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

Matteo Salvalaglio,Fabienne Bachtiger,Lunna Li,Aliff Rahimee

DOI: https://doi.org/10.26434/chemrxiv-2024-v335b

2024-07-09

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.

Chemistry

What problem does this paper attempt to address?

This paper mainly discusses the thermodynamic properties of two polymorphs (α and β) of L-Glutamic acid in solution, especially their solubility and surface tension. Based on the thermodynamics of crystal nucleation in restricted solution, a simulation framework is proposed to estimate the solubility and surface tension of organic crystals in solution using a finite number of unbiased molecular dynamics simulations at a reference temperature, and further extended to the whole solubility curve. This approach avoids complex free energy calculations and improves efficiency.

The paper applies this analytical method to study the relative thermodynamic stability and water solubility of the α and β polymorphs of L-Glutamic acid, and compares them with experimental data. The results show that while there is still a discrepancy in absolute solubility predictions, the relative stability and solubility of the polymorphs obtained through finite temperature simulations are consistent with experimental information. The authors believe that being able to efficiently obtain solubility information from limited computational experiments is a crucial step in high-throughput polymorph screening processes.

In addition, the paper introduces the theoretical background, including classical nucleation theory and the thermodynamics of crystal-solution equilibrium in finite size systems. Through a series of molecular dynamics simulations, the solubility curve, nucleation barrier, and solution stability limits at different temperatures are investigated. Finally, the paper provides a stability map of L-Glutamic acid under restricted conditions, revealing how limiting conditions affect the stability of supersaturated solutions.

Overall, this paper aims to develop a new computational method to more effectively predict and understand the nucleation and dissolution behavior of organic crystals in solution, which is of significant importance in fields such as drug discovery, materials science, and chemical engineering.