In this section
@article{Zhou:2026:10.1021/acs.jctc.6c00750,
author = {Zhou, H and Mallia, G and Harrison, NM},
doi = {10.1021/acs.jctc.6c00750},
journal = {Journal of Chemical Theory and Computation},
title = {Predictive Modeling of Paracetamol Polymorphs: LinkingBasis Set and Hamiltonian Effects to an Optimized Quasi-Harmonic LatticeDynamics Workflow},
url = {http://dx.doi.org/10.1021/acs.jctc.6c00750},
year = {2026}
}
TY - JOUR
AB - <jats:title>Abstract</jats:title> <jats:p>Thermodynamic properties of molecular crystals are the cornerstone of a broad range of applications, spanning pharmaceutics, organic semiconductors, and many more. The reliable modeling of these properties is usually hindered by collective influences of delicate dispersion interactions and insufficient prior knowledge. Providing reasonable accuracy and minimal reliance on experimental data, lattice dynamics (LD) based on the quasi-harmonic approximation (QHA) and density functional theory (DFT) has been widely and successfully adopted for crystal structure prediction (CSP). However, for polymorphic molecular crystals, it remains challenging to accurately establish the relative stability among multiple polymorphs that are distinguished by subtle energy differences. Despite recent intense efforts to achieve calculations that can reproduce observations in well-characterized systems, a comprehensive study elucidating how predictions are governed by errors within the theoretical framework and details of implementation is lacking. This undermines confidence in the transferability and efficiency of QHA LD workflows proposed in the literature. The need for reliable and predictive calculations necessitates systematic investigation and disentangling of the influences from key parameters such as the choice of basis set and Hamiltonian, so the cost and precision of DFT-based QHA LD calculations can be rationally balanced in practical workflows. In this study, paracetamol polymorphs are adopted as the prototype system, as they are well-characterized and display representative chemical interactions, where covalent bonds, hydrogen bonds, and London dispersions are present with sufficient complexity. The intermixed sources of error underlying free-energy predictions are resolved by comparing the structural, vibrational, and thermochemical properties of paracetamol polymorphs with high-quality references. A
AU - Zhou,H
AU - Mallia,G
AU - Harrison,NM
DO - 10.1021/acs.jctc.6c00750
PY - 2026///
SN - 1549-9618
TI - Predictive Modeling of Paracetamol Polymorphs: LinkingBasis Set and Hamiltonian Effects to an Optimized Quasi-Harmonic LatticeDynamics Workflow
T2 - Journal of Chemical Theory and Computation
UR - http://dx.doi.org/10.1021/acs.jctc.6c00750
UR - https://doi.org/10.1021/acs.jctc.6c00750
ER -
