Related Experiment Video
Updated: May 22, 2026
Construction and Systematical Symmetric Studies of a Series of Supramolecular Clusters with Binary or Ternary Ammonium Triphenylacetates
Published on: February 15, 2016
The Structure, Thermodynamics and Solubility of Organic Crystals from Simulation with a Polarizable Force Field
Michael J Schnieders1, Jonas Baltrusaitis, Yue Shi
1Department of Biomedical Engineering, University of Texas at Austin, Austin, TX 78712.
Abstract:
An important unsolved problem in materials science is prediction of the thermodynamic stability of organic crystals and their solubility from first principles. Solubility can be defined as the saturating concentration of a molecule within a liquid solvent, where the physical picture is of solvated molecules in equilibrium with their solid phase. Despite the importance of solubility in determining the oral bioavailability of pharmaceuticals, prediction tools are currently limited to quantitative structure-property relationships that are fit to experimental solubility measurements. For the first time, we describe a consistent procedure for the prediction of the structure, thermodynamic stability and solubility of organic crystals from molecular dynamics simulations using the polarizable multipole AMOEBA force field. Our approach is based on a thermodynamic cycle that decomposes standard state solubility into the sum of solid-vapor sublimation and vapor-liquid solvation free energies [Formula: see text], which are computed via the orthogonal space random walk (OSRW) sampling strategy. Application to the n-alkylamides series from aeetamide through octanamide was selected due to the dependence of their solubility on both amide hydrogen bonding and the hydrophobic effect, which are each fundamental to protein structure and solubility. On average, the calculated absolute standard state solubility free energies are accurate to within 1.1 kcal/mol. The experimental trend of decreasing solubility as a function of n-alkylamide chain length is recapitulated by the increasing stability of the crystalline state and to a lesser degree by decreasing favorability of solvation (i.e. the hydrophobic effect). Our results suggest that coupling the polarizable AMOEBA force field with an orthogonal space based free energy algorithm, as implemented in the program Force Field X, is a consistent procedure for predicting the structure, thermodynamic stability and solubility of organic crystals.
Related Concept Videos
03:00Solubility
A solution is a homogeneous mixture composed of a solvent, the major component, and a solute, the minor component. The physical state of a solution—solid, liquid, or gas—is typically the same as that of the solvent. Solute concentrations are often described with qualitative terms such as dilute (of relatively low concentration) and concentrated (of relatively high concentration).
In a solution, the solute particles (molecules, atoms, and/or ions)...
02:05Entropy and Solvation
01:21Polymer Classification: Crystallinity
Crystalline domains are the regions where polymer chains are aligned in an orderly manner and held together in proximity by intermolecular forces. For example, chains in the crystalline domains of polyethylene and nylon are bound together by van der Waals...
03:37Molecular Shape and Polarity
02:54Molecular and Ionic Solids
Molecular Solids
Molecular crystalline solids, such as ice, sucrose (table sugar), and iodine, are solids that are composed of neutral molecules as their constituent units. These molecules are held together by weak intermolecular forces such as London dispersion forces, dipole-dipole interactions, or hydrogen bonds, which...
02:28Intermolecular Forces in Solutions
When the strengths of the intermolecular forces of attraction between solute and solvent species in a solution are no different than those present in the separated components, the solution is formed with no accompanying energy change. Such a solution is called an ideal solution. A mixture of ideal gases (or gases such as helium and argon,...