## scientific commentaries

## Looking at local classical and quantum forces in stable crystals using multipole-model refined electron densities and orbital-free DFT approximations

^{a}CNR-SCITEC, via Golgi section, via Golgi 19, Milano, Italy 20134, Italy^{*}Correspondence e-mail: c.gatti@scitec.cnr.it

Two distinct approaches, that of *energy* and that of *force*, are adopted in quantum mechanics to get insights on chemical processes. In the second one, the net forces acting on the electrons and nuclei in a system (Ehrnefest and Hellmann–Feynman forces, respectively) are determined and a *local* version of the approach, in terms of force density fields rather than forces, has also been proposed for electrons. This is the path followed by Tsirelson & Stash (2020) in this issue of *Acta Crystallographica Section B*, to study for the first time the spatial distribution of the electronic forces of different nature acting in stable crystals. Interestingly, by relying on approximations taken from orbital-free DFT, all components of the inner-crystal force can be easily retrieved from multipole-model refined experimental electron densities and their derivatives. No less important is that these calculations are becoming easily doable for any X-ray density crystallographer thanks to a new version of the computer program *WinXPRO*, purposely developed in the study which is discussed in this commentary.

Keywords: forces in crystals; orbital-free DFT; quantum crystallography; atoms in molecules topology.

In their article, Tsirelson and Stash (T&S) describe a method to obtain a local mapping of electrostatic and quantum forces acting on electrons in stable crystals. To achieve their goal, T&S have carefully selected and ingeniously combined several ingredients that might not be immediately familiar to crystallographers. The ingredients include the quantum mechanical *force approach* in chemistry (Deb, 1973), the theory of *gradient dynamical systems* (Abraham & Shaw, 1992), *Orbital-Free Density Functional Theory* (Wesolowski & Wang, 2013) and the old, but recently broadened and revitalized field of *Quantum Crystallography* (Genoni *et al.*, 2018). Each of these bits and pieces of science deserves a brief introduction on its own, before analysing the role they play when combined together (Fig. 1).

In quantum chemistry, rather than using the more customary energy approach, chemical bonding and molecular processes may be studied using the force approach, that is in terms of the forces of various nature acting on the electrons and nuclei in a system. In the *local force* form adopted by T&S, forces acting on electrons are expressed in terms of force density fields (Deb, 1979; Bader *et al.*, 2007; Pendás & Hernández-Trujillo, 2012; Pendás *et al.*, 2016). The latter may be obtained either from the expectation values of the gradients of the corresponding potential operators or, as performed by T&S, as gradients of the associated scalar potential fields.

These scalar potentials and their gradient fields (the forces) are analysed by T&S in terms of their topology via the theory of *gradient dynamical systems*, using its typical entities (attractors, repellers, basins, gradient paths, critical points) and concepts (homeomorphism). It is a tool common to most of the Quantum Chemical Topology methods, such as Bader's Quantum Theory of Atoms in Molecules (QTAIM) (Bader, 1990). *Attractors* and *repellers* are those local potential maxima and local potential minima where gradient paths (the force lines) converge to and are expelled from, leading to a partitioning of chemical systems into space-exhaustive pieces, named *basins*. Each basin is traversed by all ascending gradient paths (GPs) terminating at its own 3D attractor (an alternative partitioning in terms of descending GPs and 3D repellers is equally possible). In most cases the basins bear a chemical meaning. Those with 3D attractors are often atom-like in shape and in QTAIM they are identified as the (quantum) atoms in molecules. The *critical points* (CPs) are those points where the gradient of the potential, the force, vanishes and are characterized by their signature. The signature is given by the dimensionality of the space traversed by the ascending GPs expelled from the CP minus that of the space traversed by the ascending GPs with terminus at the CP. Signatures −3 and +3 correspond to attractors (local maxima) and repellers (local minima), while the −1 and +1 signatures denote first- and second-order saddle points, respectively. If two potential fields exhibit the same number of CPs for each signature, the fields are said to be *homeomorphic* within the system where this occurs. Basins sharing a portion of their bounding surface imply a CP of −1 signature lying on their common surface and being the potential maximum for it. Two unique ascending GPs originate at this CP. One GP ends up at the attractor of one of the two basins sharing the surface and the other at the attractor of the remaining basin. The union of such paths is a line joining these two attractors and with maximum potential value relative to any lateral displacement from the line. In QTAIM it is known as a *bond path*, a line of maximum electron density ρ linking the nuclei of a pair of atoms. Accumulation of ρ along the bond path has been associated with chemical bonding (Bader, 1990). T&S discuss the analogy and dissimilarity of the topology of the scalar potential and force fields they analyse with that typical of QTAIM (where ρ is taken as the potential and its gradient ∇ρ as the corresponding force).

Orbital-free density functional theory (OF-DFT) and the methods derived thereof represent a fascinating challenge to go beyond the established routes in computational DFT. In principle, Hohenberg–Kohn theorems make it possible to describe the ground state of a multi-electron system exactly without constructing the wavefunction explicitly. Yet, the extreme difficulty of building sufficiently accurate approximations to the electron density (ED) functional for the ). Besides, the and the electrostatic, exchange and densities, all entering with their corresponding functionals in OF-DFT approaches, were demonstrated to be of great importance for interpreting key chemical concepts (Bader, 1990; Gatti, 2005; Wesolowski & Wang, 2013; Pendás *et al.*, 2016). Using information from experimental EDs and their derivatives, Tsirelson (2007) has in the past pioneered the study of bonding in crystals in terms of these various local OF-DFT energies. T&S now take a significant step forward. Rather than the energies, they study the OF-DFT potentials and forces derived thereof, using again experimentally derived densities and examining which forces behave as *heterotropic* and which as *homotropic* (Revici, 1961). The latter tend to preserve homogenous disorganization of the matter, while the former favour its local organization, such as the electronic shell structure around nuclei.

Next to `orbital-free', the little known discipline of Quantum Crystallography (QCr) appears into the title of T&S's paper before the aim of `viewing forces in crystals' is declared. QCr was originally defined as the operation of using `crystallographic information to enhance quantum mechanical calculations and the information derived from them' (Massa *et al.*, 1995) or, conversely, of that of performing `*ad hoc* quantum mechanical calculations to enhance the accuracy of the crystallographic (Huang *et al.*, 1999). None of these definitions strictly applies to the work performed by T&S, yet it surely does to the more recent, highly broadened QCr view (Genoni *et al.*, 2018), embracing many other active and emerging research areas involving quantum mechanics and scattering experiments. Significantly, in 2017 the IUCr Commission on Charge and Spin Momentum Densities was renamed as the Commission on Quantum Crystallography. Orbital-free QCr can then be seen as a successful attempt of using precious information obtained from crystallographic methods for evaluating potential and forces in crystals within a fundamental quantum mechanical approach.

After illustrating their multi-faceted scientific framework, the results obtained by T&S deserve a brief overview and concluding comments. Within DFT, and considering a number of simplifications (Tsirelson & Stash, 2020), the local one-electron potential μ(**r**) is given by the sum of a kinetic potential, ν_{kin}(**r**), and of the so-called potential acting on an electron in a molecule (PAEM; Zhao & Yang, 2014), ν_{PAEM}(**r**). For a stable system, μ(**r**) vanishes everywhere. Thereby, the two electron potentials balance themselves everywhere, ν_{kin}(**r**) = −ν_{PAEM}(**r**), and so do their corresponding forces −*f*_{kin}(**r**) = *f*_{PAEM}(**r**) [with ∇ν_{kin}(**r**) = −*f*_{kin}(**r**) and −∇ν_{PAEM}(**r**) = *f*_{PAEM}(**r**)]. Since ν_{kin} and ν_{PAEM} differ in sign, they have the same total number of CPs and placed at the same locations, but each with reversed signature. For instance, the local maxima of one potential become the local minima of the other and *vice versa*. The kinetic potential includes two terms, one due to the electron fluctuations arising from the uncertainty principle and the other to the corrections arising from the wavefunction antisymmetry requirements. Both terms are clearly quantum mechanical in nature. Instead, PAEM is composed of the electrostatic, the exchange and the correlation potential contributions, the first one being purely classical.

T&S have analysed all ν_{kin} and ν_{PAEM} components for two simple stable crystals, diamond and rock salt, using ED and its derivatives obtained from high-resolution X-ray diffraction experiments, followed by an accurate multipole-model of the electronic structure factors. Many are the intriguing results of their study. First, it is found that the forces of kinetic origin draw electrons from the nuclei and are heterotropic in nature, while all *f*_{PAEM} components behave as homotropic. The latter push electrons to the nuclei that therefore act as sinks for the PAEM field, rather than as local maxima for ν_{kin}. Hence, only the kinetic potential and its associated forces are able to create the electronic shells yielding atomic-like structures and basins resembling those of ρ in the QTAIM, yet with noticeable differences (Tsirelson & Stash, 2020). Both ν_{kin} and ν_{PAEM} fields concur to the accumulation of electronic charge along curved lines (GPs) joining nearest atoms. However, not all atomic pairs linked by bond paths in the ρ field (QTAIM) are found to be joined by the corresponding GPs in the ν_{kin} field. This fact has an impact on the topologies of the two fields that are therefore not always homeomorphic. Even more importantly, it opens the question on the physical meaning of both the lack of a shared portion of bounding surface and of a common CP of signature −1 in some circumstances (in the present case, the Cl⋯Cl interaction in the NaCl crystal, which is missing in the ν_{kin} field and, instead, present in the ρ field). Interestingly, the missing CPs in the ν_{kin} field topologically correspond to the QTAIM bond critical points (bcps), whose meaning is and has been widely debated in the literature [see for instance, Pendás *et al.* (2007)]. Understanding the reasons behind the possible dissimilarities of the ρ and the ν_{kin} fields could shed further light on the physical interpretation of bcps. It will certainly be a challenge for forthcoming contributions in the area pioneered by T&S in their paper.

### Supporting information

Originally published version of this article. DOI: https://doi.org//10.1107/S2052520620012895/bm5134sup1.pdf

### References

Abraham, R. H. & Shaw, C. D. (1992). *Dynamics: The Geometry of Behavior.* Redwood City, CA: Addison Wesley. Google Scholar

Bader, R. F. W. (1990). *Atoms in Molecules: A Quantum Theory.* New York: Oxford University Press. Google Scholar

Bader, R. F. W., Hernández-Trujillo, J. & Cortés-Guzmán, F. (2007). *J. Comput. Chem.* **28**, 4–14. Web of Science CrossRef PubMed CAS Google Scholar

Deb, B. M. (1973). *Rev. Mod. Phys.* **45**, 22–43. CrossRef CAS Web of Science Google Scholar

Deb, B. M. (1979). *J. Phys. B At. Mol. Phys.* **12**, 3857–3871. CrossRef CAS Web of Science Google Scholar

Gatti, C. (2005). *Z. Kristallogr.* **220**, 399–457. Web of Science CrossRef CAS Google Scholar

Genoni, A., Bučinský, L., Claiser, N., Contreras-García, J., Dittrich, B., Dominiak, P. M., Espinosa, E., Gatti, C., Giannozzi, P., Gillet, J.-M., Jayatilaka, D., Macchi, P., Madsen, A. Ø., Massa, L., Matta, C. F., Merz, K. M. Jr, Nakashima, P. N. H., Ott, H., Ryde, U., Schwarz, K., Sierka, M. & Grabowsky, S. (2018). *Chem. Eur. J.* **24**, 10881–10905. Web of Science CrossRef CAS PubMed Google Scholar

Huang, L., Massa, L. & Karle, J. (1999). *Int. J. Quant. Chem.* **73**, 439–450. CrossRef CAS Google Scholar

Massa, L., Huang, L. & Karle, J. (1995). *Int. J. Quantum Chem.* **56**, 371–384. CrossRef Google Scholar

Pendás, A. M., Francisco, E., Blanco, M. A. & Gatti, C. (2007). *Chem. Eur. J.* **13**, 9362–9371. Web of Science PubMed Google Scholar

Pendás, A. M., Francisco, E., Gallo Bueno, A., Guevara Vela, J. M. & Costales, A. (2016). *Applications of Topological Methods in Molecular Chemistry*, edited by R. Chauvin, C. Lepetit, B. Silvi and E. Alikhani, pp. 131–150. Springer International Publishing, Switzerland. Google Scholar

Pendás, A. M. & Hernández-Trujillo, J. (2012). *J. Chem. Phys.* **137**, 134101. PubMed Google Scholar

Revici, E. (1961). *Research in Physiopathology as Basis of Guided Chemotherapy with Special Application to Cancer.* New York, London and Toronto: Van Nostrand Co Inc. (See https://www.biri.org/resources/revici-book.). Google Scholar

Tsirelson, V. & Stash, A. (2020). *Acta Cryst.* B**76**, 769–778. CrossRef IUCr Journals Google Scholar

Tsirelson, V. G. (2007). *The Quantum Theory of Atoms in Molecules: From Solid State to DNA and Drug Design*, edited by C. Matta and R. Boyd, pp. 259–283. Weinheim: Wiley-VCH. Google Scholar

Wesolowski, T. A. & Wang, Y. A. (2013). Editors. *Recent Progress in Orbital-Free Density Functional Theory.* World Scientific. Google Scholar

Zhao, D. X. & Yang, Z. Z. (2014). *J. Phys. Chem. A*, **118**, 9045–9057. CrossRef CAS PubMed Google Scholar

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.