On the Anisotropy of van der Waals Atomic Radii of O, S, Se, F, Cl, Br, and I

Analyzing Atomic Interactions in Molecules as Learned by Neural Networks

While machine learning (ML) models have been able to achieve unprecedented accuracies across various prediction tasks in quantum chemistry, it is now apparent that accuracy on a test set alone is not a guarantee for robust chemical modeling such as stable molecular dynamics (MD). To go beyond accuracy, we use explainable artificial intelligence (XAI) techniques to develop a general analysis framework for atomic interactions and apply it to the SchNet and PaiNN neural network models. We compare these interactions with a set of fundamental chemical principles to understand how well the models have learned the underlying physicochemical concepts from the data. We focus on the strength of the interactions for different atomic species, how predictions for intensive and extensive quantum molecular properties are made, and analyze the decay and many-body nature of the interactions with interatomic distance. Models that deviate too far from known physical principles produce unstable MD trajectories, even when they have very high energy and force prediction accuracy. We also suggest further improvements to the ML architectures to better account for the polynomial decay of atomic interactions.

van der Waals Radii of Free and Bonded Atoms from Hydrogen (Z = 1) to Oganesson (Z = 118)

Reliable numerical values of van der Waals (vdW) radii are required for constructing empirical force fields, vdW-inclusive density functional, and quantum-chemical methods, as well as for implicit solvent models. However, multiple definitions exist for vdW radii, involving either equilibrium or the closest contact distances between free or bonded atoms within molecules or crystals. For the paradigmatic case of the hydrogen atom, its reported vdW radius fluctuates between 2.15 and 3.70 Bohr depending on the definition, leading to a high uncertainty in calculations and different conceptual interpretations of noncovalent interactions. In this work, we systematically review different definitions and methodologies to establish the free and bonded vdW radii for hydrogen, based on equilibrium vdW distances in noncovalently bonded molecules, enveloping electron density cutoffs, noncovalent positron bonds in hydrogen anion dimer, vacuum virtual photon cloud caused by the hydrogen atom, and atomic dipole polarizability. By doing so, we show that the vdW radius of the free hydrogen atom is 3.16 ± 0.06 Bohr. By employing the most general and elegant definition of atomic vdW radius as a function of the atomic polarizability, we tabulate consistent values of vdW radii for all atoms in the periodic table up to Z = 118.

Competition of interactions and a new high-temperature phase of selenourea

The aggregation of molecules is usually associated with a specific type of interaction, which can be altered by thermodynamic conditions. Under normal conditions, the crystal structure of selenourea, SeC(NH₂)₂, phase α is trigonal, space group P3₁, Z = 27. Its large number of independent molecules (Z_α' = 9) can be associated with the formation of an NH...N hydrogen bond substituting one of 36 independent NH...Se hydrogen bonds, which prevail among intermolecular interactions. Phase α approximates the trigonal structure with a threefold smaller unit cell (Z = 9), which in turn approximates another still threefold smaller unit cell (Z = 3). The temperature-induced transformations of selenourea have been characterized by calorimetry and by performing 21 single-crystal X-ray diffraction structural determinations as a function of temperature. At 381.0 K, phase α undergoes a first-order displacive transition to phase γ, with space group P3₁21 and Z reduced to 9, when the NH...N bond is broken and an NH...Se bond is formed in its place. Previously, an analogous competition was observed between NH...N and NH...O hydrogen bonds in high-pressure phase III of urea. The lattice vectors along the (001) plane in low- and high-temperature phases of selenourea are related by a similarity rule, while the lattice dimensions along direction c are not affected. This similarity rule also applies to the structures of phase γ and hypothetical phase δ (Z = 3). The thermally controlled transition between enantiomorphic phases of selenourea contrasts with its high-pressure transition at 0.21 GPa to a centrosymmetric phase β, where both the NH...Se and NH...N bonds are present. The compression and heating reduce the number of independent molecules from Z' = 9 in phase α, to Z' = 2 in phase β and to Z' = 1.5 in phase γ.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

The π-hole revisited

It follows from the Schrödinger equation that the forces operating within molecules and molecular complexes are Coulombic, which necessarily entails both electrostatics and polarization. A common and important class of molecular complexes is due to π-holes. These are molecular regions of low electronic density that are perpendicular to planar portions of the molecular frameworks. π-Holes often have positive electrostatic potentials associated with them, which result in mutually polarizing attractive forces with negative sites such as lone pairs, π electrons or anions. In many molecules, π-holes correspond to a flattening of the electronic density surface but in benzene derivatives and in polyazines the π-holes are craters above and below the rings. The interaction energies of π-hole complexes can be expressed quite well in terms of regression relationships that account for both the electrostatics and the polarization. There is a marked gradation in the interaction energies, from quite weak (about -2 kcal mol-1) to relatively strong (about -40 kcal mol-1). Gradations are also evident in the ratios of the intermolecular separations to the sums of the respective van der Waals radii and in the gradual transition of the π-hole atoms from trigonal to quasi-tetrahedral configurations. These trends are consistent with the concept that chemical interactions form a continuum, from very weak to very strong.

Can Counter-Intuitive Halogen Bonding Be Coulombic?

We use the term "counter-intuitive" to describe an intermolecular interaction in which the electrostatic potentials of the interacting regions of the ground-state molecules have the same sign, both positive or both negative. In the present work, we consider counter-intuitive halogen bonding with nitrogen bases, in which both the halogen σ-hole and the nitrogen lone pair have negative potentials on their molecular surfaces. We show that these interactions can be treated as Coulombic despite the apparent repulsion between the ground-state molecules, provided that both electrostatics and polarization are explicitly taken into account. We demonstrate first that the energies of 20 counter-intuitive interactions with four nitrogen bases can be expressed very well in terms of just two molecular properties: the electrostatic potential of the halogen σ-hole and the average polarizability of the nitrogen base. Then we show that the same two properties can also represent the energies of an expanded data base that includes the 20 counter-intuitive plus an additional 20 weak and moderately-strong intuitive halogen bonding interactions (in which the σ-hole potentials are now positive).

Structural biology of DNA abasic site protection by SRAP proteins

Abasic (AP) sites are one of the most frequently occurring types of DNA damage. They lead to DNA strand breaks, interstrand DNA crosslinks, and block transcription and replication. Mutagenicity of AP sites arises from translesion synthesis (TLS) by error-prone bypass polymerases. Recently, a new cellular response to AP sites was discovered, in which the protein HMCES (5-hydroxymethlycytosine (5hmC) binding, embryonic stem cell-specific) forms a stable, covalent DNA-protein crosslink (DPC) to AP sites at stalled replication forks. The stability of the HMCES-DPC prevents strand cleavage by endonucleases and mutagenic bypass by TLS polymerases. Crosslinking is carried out by a unique SRAP (SOS Response Associated Peptidase) domain conserved across all domains of life. Here, we review the collection of recently reported SRAP crystal structures from human HMCES and E. coli YedK, which provide a unified basis for SRAP specificity and a putative chemical mechanism of AP site crosslinking. We discuss the structural and chemical basis for the stability of the SRAP DPC and how it differs from covalent protein-DNA intermediates in DNA lyase catalysis of strand scission.

Interaction and Polarization Energy Relationships in σ-Hole and π-Hole Bonding

We demonstrate that a wide range of σ- and π-hole interaction energies can be related to (a) the electrostatic potentials and electric fields of the σ- and π-hole molecules at the approximate positions of the negative sites and (b) the electrostatic potentials and polarizabilities of the latter. This is consistent with the Coulombic nature of these interactions, which should be understood to include both electrostatics and polarization. The energies associated with polarization were estimated and were shown to overall be greater for the stronger interactions; no new factors need be introduced to account for these. All of the interactions can be treated in the same manner.

A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

The founding in 1965 of what is now called the Cambridge Structural Database (CSD) has reaped dividends in numerous and diverse areas of chemical research. Each of the million or so crystal structures in the database was solved for its own particular reason, but collected together, the structures can be reused to address a multitude of new problems. In this Review, which is focused mainly on the last 10 years, we chronicle the contribution of the CSD to research into molecular geometries, molecular interactions, and molecular assemblies and demonstrate its value in the design of biologically active molecules and the solid forms in which they are delivered. Its potential in other commercially relevant areas is described, including gas storage and delivery, thin films, and (opto)electronics. The CSD also aids the solution of new crystal structures. Because no scientific instrument is without shortcomings, the limitations of CSD research are assessed. We emphasize the importance of maintaining database quality: notwithstanding the arrival of big data and machine learning, it remains perilous to ignore the principle of garbage in, garbage out. Finally, we explain why the CSD must evolve with the world around it to ensure it remains fit for purpose in the years ahead.

σ-Holes vs. Buildups of Electronic Density on the Extensions of Bonds to Halogen Atoms

Our discussion focuses upon three possible features that a bonded halogen atom may exhibit on its outer side, on the extension of the bond. These are (1) a region of lower electronic density (a σ-hole) accompanied by a positive electrostatic potential with a local maximum, (2) a region of lower electronic density (a σ-hole) accompanied by a negative electrostatic potential that also has a local maximum, and (3) a buildup of electronic density accompanied by a negative electrostatic potential that has a local minimum. In the last case, there is no σ-hole. We show that for diatomic halides and halogen-substituted hydrides, the signs and magnitudes of these maxima and minima can be expressed quite well in terms of the differences in the electronegativities of the halogen atoms and their bonding partners, and the polarizabilities of both. We suggest that the buildup of electronic density and absence of a σ-hole on the extension of the bond to the halogen may be an operational indication of ionicity.

Comparison between Tetrel Bonded Complexes Stabilized by σ and π Hole Interactions

The σ-hole tetrel bonds formed by a tetravalent molecule are compared with those involving a π-hole above the tetrel atom in a trivalent bonding situation. The former are modeled by TH₄, TH₃F, and TH₂F₂ (T = Si, Ge, Sn) and the latter by TH₂=CH₂, THF=CH₂, and TF₂=CH₂, all paired with NH₃ as Lewis base. The latter π-bonded complexes are considerably more strongly bound, despite the near equivalence of the σ and π-hole intensities. The larger binding energies of the π-dimers are attributed to greater electrostatic attraction and orbital interaction. Each progressive replacement of H by F increases the strength of the tetrel bond, whether σ or π. The magnitudes of the maxima of the molecular electrostatic potential in the two types of systems are not good indicators of either the interaction energy or even the full Coulombic energy. The geometry of the Lewis acid is significantly distorted by the formation of the dimer, more so in the case of the σ-bonded complexes, and this deformation intensifies the σ and π holes.

The σ-hole revisited

A covalently-bonded atom typically has a region of lower electronic density, a "σ-hole," on the side of the atom opposite to the bond, along its extension. There is frequently a positive electrostatic potential associated with this region, through which the atom can interact attractively but noncovalently with negative sites. This positive potential reflects not only the lower electronic density of the σ-hole but also contributions from other portions of the molecule. These can significantly influence both the value and also the angular position of the positive potential, causing it to deviate from the extension of the covalent bond. We have surveyed these effects, and their consequences for the directionalities of subsequent noncovalent intermolecular interactions, for atoms of Groups IV-VII. The overall trends are that larger deviations of the positive potential result in less linear intermolecular interactions, while smaller deviations lead to more linear interactions. We find that the deviations of the positive potentials and the nonlinearities of the noncovalent interactions tend to be greatest for atoms of Groups V and VI. We also present arguments supporting the use of the 0.001 a.u. contour of the electronic density as the molecular surface on which to compute the electrostatic potential.

Cooperative intermolecular S–Cl⋯O and F⋯F associations in the crystal packing of α,ω-di(sulfonyl chloride) perfluoroalkanes, ClSO2(CF2)nSO2Cl, where n = 4, 6

Halogen bonding between neighboring sulfonyl chloride groups and short fluorine–fluorine contacts supports crystal formation in the title compounds.

σ-Hole Interactions: Perspectives and Misconceptions

After a brief discussion of the σ-hole concept and the significance of molecular electrostatic potentials in noncovalent interactions, we draw attention to some common misconceptions that are encountered in that context: (1) Since the electrostatic potential reflects the contributions of both the nuclei and the electrons, it cannot be assumed that negative potentials correspond to “electron-rich” regions and positive potentials to “electron-poor” ones; (2) The electrostatic potential in a given region is determined not only by the electrons and nuclei in that region, but also by those in other portions of the molecule, especially neighboring ones; (3) A σ-hole is a region of lower electronic density on the extension of a covalent bond, not an electrostatic potential; (4) Noncovalent interactions are between positive and negative regions, which are not necessarily associated with specific atoms, so that “close contacts” between atoms do not always indicate the actual interactions.

On the role of charge transfer in halogen bonding

We have quantified the energetic contribution of charge transfer to halogen bonding to be about 10% of the interaction energy.

Alkyl groups as electron density donors in π-hole bonding

A combined structural and computational analysis has demonstrated that alkyl groups can act as Lewis bases in π-hole bonding.

Beyond Born-Mayer: Improved Models for Short-Range Repulsion in ab Initio Force Fields

Short-range repulsion within intermolecular force fields is conventionally described by either Lennard-Jones (A/r(12)) or Born-Mayer (A exp(-Br)) forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of intermolecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, and robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Finally, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach.