Lattice energies of molecular crystals can also be obtained with CrystalExplorer17, using the same CE-B3LYP interaction energies that underpin the calculation of energy frameworks.
Using four different benchmark sets of molecular crystals, we have established the level of confidence for lattice energies estimated using CE-B3LYP model energies and experimental crystal structures. They compare very well with available benchmark estimates derived from sublimation enthalpies, and in many cases they are comparable with, and sometimes better than, more computationally demanding approaches, such as those based on periodic DFT+dispersion methodologies. The performance over the complete set of 110 crystals indicates a mean absolute deviation from benchmark energies of only 6.6 kJ/mol. That reference also discussed applications to polymorphic crystals and some much larger molecules.
At present the calculation of lattice energies has not been automated, but it is a relatively straightforward procedure. CE-B3LYP lattice energies are computed by direct summation of interaction energies in CrystalExplorer over molecules interacting with a central molecule until is converged to better than ~1 kJ/mol, using a cutoff based on the separation of molecular centroids, . Sums over Edis and Erep converge rapidly, but the sum over terms is much longer range, with highly dipolar compounds (e.g., HCN, formamide, zwitterionic amino acids) requiring much greater cutoff distances than nondipolar molecules (e.g., anthracene, coronene). For molecules separated by more than 12 Å, is computed using a distributed multipole model of the electron density.
The first step is to construct a cluster of sufficient size around a molecule, typically with the Generate Atoms Within Radius button, making sure that all molecule fragments are completed . For non-dipolar molecules a radius of 20 Å should be adequate, but more dipolar molecules (e.g., amino acids, zwitterions) may require a cluster radius of 30 Å or more.
With the central molecule still selected, clicking on Calculate Energies will open the dialog box for energies. Select Accurate [B3LYP/6-31G(d,p)]. This is the only energy model that has been calibrated against experimental sublimation enthalpies. Clicking on OK will initiate wavefunction calculation, followed by computation of all unique pair energies of the central molecule with other molecules in the cluster. This (naturally) takes a little time, and progress of the computation can be followed by the dialog and progress bar at the bottom of the graphics window.
The Information window will open automatically when all unique pair energies have been computed. An example for formamide is shown on the right. The text in this window can be copied and pasted into a spreadsheet application, and from this the lattice energy obtained as one-half the product of the two columns labelled N and Etot, where N is the number of molecule pairs in the cluster with that particular interaction energy. In the same manner the lattice energy can be broken down into separate contributions, electrostatic, dispersion, etc. (Note that the number of decimal places for energy quantities in the Information window can be modified using the Preferences dialog).
Pay careful attention to whether the lattice energy estimate you have obtained is adequately converged. You can do this by sorting the energies by separation of molecular pairs, and computing partial sums, or by extending the cluster by several Å, computing the additional energies and the revised estimate of the lattice energy, then repeating until convergence. Suitably converged lattice energies should show consistent changes of less than 1 kJ mol–1.
The procedure here is similar to that above, but more complicated because separate clusters need to be created for each unique molecule.
For each unique molecule you need to do the following:
For unit cells with nonzero dipole moments (i.e., polar space groups such as , , and ) the simple lattice summation described above is not representative of the entire crystal, and an additional consideration is required, the cell dipole energy, dipole. Numerous examples of the additional term can be found in the Supporting Information to the lattice energy paper. Calculation of this term is straightforward, but it requires an estimate of the unit cell dipole moment. This can be obtained by summing the components along the polar direction of all molecular dipole moments in the unit cell, and that information is easily obtained from the Gaussian output files for each molecule.
The lattice energy can be related to the sublimation enthalpy by taking into account the electronic relaxation energy associated with the crystal to gas transition, as well the difference between vibrational energies for the two phases,
and this can be (reasonably well) approximated by
The important point here is that for crystals where little molecular geometry change occurs between crystal and gas phase, the sublimation enthalpy is a positive quantity, whose magnitude at room temperature is ~5 kJ/mol less than that of the lattice energy.
Conformational polymorphism is common in pharmaceutical compounds, a consequence of their considerable internal flexibility, and any attempt to estimate relative polymorph stability must also incorporate the energy differences between the different conformations found in the polymorphs. We have discussed one way of estimating , with application to the ROY polymorphs.