Abstract:Efficiently obtaining atomic-scale thermodynamic parameters characterizing crystallization from solution is key to developing the modeling strategies needed in the quest for digital design strategies of industrial crystallization processes. Based on the thermodynamics of crystal nucleation in confined solutions, we develop a simulation framework to efficiently estimate the solubility and surface tension of organic crystals in solution from a few unbiased molecular dynamics simulations at a reference temperature. We then show that such a result can be extended with minimal computational overhead to capture the whole solubility curve. This enables an efficient and self-consistent estimate of the solubility and limit of solution stability associated with crystal nucleation in molecular systems from equilibrium molecular dynamics without requiring sophisticated free energy calculations. We apply our analysis to investigate the relative thermodynamics stability and aqueous solubility of the α and β polymorphs of L-Glutamic acid. Our analysis enables an efficient appraisal of emergent ensemble properties associated with the thermodynamics of nucleation from solutions against experimental data, demonstrating that while the absolute solubility is still far from being quantitatively captured by an off-the-shelf point charge transferable forcefield, the relative polymorphic stability and solubility obtained from finite temperature simulation are consistent with the experimentally available information on glutamic acid. We foresee the ability to efficiently obtain solubility information from a limited number of computational experiments as a key component of high-throughput computational polymorph screenings.

Polymorph-Specific Solubility Prediction of Urea Using Constant Chemical Potential Molecular Dynamics Simulations

Polymorph-Specific Solubility Prediction using Constant Chemical Potential Molecular Dynamics

Solubility Prediction of Organic Molecules with Molecular Dynamics Simulations

A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals

Studies on the Conformational Interconversion of Sulfonylurea Molecules in Solution by Molecular Dynamics Simulation Method

Comparative assessment of physics-based in silico methods to calculate relative solubilities

Free Energy Perturbation Approach for Accurate Crystalline Aqueous Solubility Predictions

Properties of compatible solutes in aqueous solution

Solution Thermodynamics of L-Glutamic Acid Polymorphs from Finite-Sized Molecular Dynamics Simulations

Organic Solubility Prediction at the Limit of Aleatoric Uncertainty

Mechanism of Urea Crystal Dissolution in Water from Molecular Dynamics Simulation.

Predicting Solubility Limits of Organic Solutes for a Wide Range of Solvents and Temperatures

Molecular Mechanism of Gas Solubility in Liquid: Constant Chemical Potential Molecular Dynamics Simulations

Study the solubility of pharmaceutical ingredients and their eutectic mixtures: An in-depth density functional theory and molecular dynamics simulations approaches

Intermolecular Interactions as a Measure of Dapsone Solubility in Neat Solvents and Binary Solvent Mixtures

Unravelling the structural and dynamical properties of concentrated aqueous ammonium nitrate solutions: MD simulation studies

Modeling the solubility of pharmaceuticals in pure solvents and solvent mixtures for drug process design

Analytical and Computational Methods for the Determination of Drug-Polymer Solubility and Miscibility

Co‐solutes Induced Changes in the Properties of Polymeric Solution and Water Dynamics

Solubility of Organic Salts in Solvent–Antisolvent Mixtures: A Combined Experimental and Molecular Dynamics Simulations Approach

Assessment of Solubility Behavior of a Copolymer in Supercritical CO 2 and Organic Solvents: Neural Network Prediction and Statistical Analysis