Magnetism and Bond Order in Diatomic Molecules Described by Semiclassical Electrons

Adventures in interdisciplinary science: a half century at the nexus between chemistry, physics and biology

A look back over five decades of research.

Emergence of Linnett's "double quartets" from a model of "Lewis dots"

Chemists routinely explicate molecular structures and chemical reactions in terms of the propensities of semiclassical valence electrons (aka "Lewis dots"). Typically, the electrons are viewed as forming spin pairs and recent efforts to translate this concise and intuitive qualitative picture into an efficient and relatable quantitative model have made good progress. But electrons are not always paired and advanced quantum calculations have shown that this is so even in small diamagnetic species such as dicarbon and benzene. Here we show that the latest semiclassical model for paired electrons can clarify the limitations on pairing simply by dissecting the elements of the interparticle potentials. Although not trained to do so, these elements produce a Linnett-like benzene, with three valence electrons in each CC bond, when the electrons are free to move singly. At the same time, sustaining higher order bonds with independently mobile electrons requires adjustments in the details of the model potentials at short distances. This is addressed with new training data and new forms for the contributions from Coulomb integrals. Although trained on hydrogen and carbon species separately, the combination applied to ethyne predicts the pairing of spins in the CH bonds and the dispersion of spins in the CC bond that is found in ab initio calculations. This adjusted force field is named LINNETT, in appreciation of Linnett's insightful double quartet interpretation of the Lewis octet.

Modeling Electronic Response Properties with an Explicit-Electron Machine Learning Potential

Explicit-electron force fields introduce electrons or electron pairs as semiclassical particles in force fields or empirical potentials, which are suitable for molecular dynamics simulations. Even though semiclassical electrons are a drastic simplification compared to a quantum-mechanical electronic wave function, they still retain a relatively detailed electronic model compared to conventional polarizable and reactive force fields. The ability of explicit-electron models to describe chemical reactions and electronic response properties has already been demonstrated, yet the description of short-range interactions for a broad range of chemical systems remains challenging. In this work, we present the electron machine learning potential (eMLP), a new explicit electron force field in which the short-range interactions are modeled with machine learning. The electron pair particles will be located at well-defined positions, derived from localized molecular orbitals or Wannier centers, naturally imposing the correct dielectric and piezoelectric behavior of the system. The eMLP is benchmarked on two newly constructed data sets: eQM7, an extension of the QM7 data set for small molecules, and a data set for the crystalline β-glycine. It is shown that the eMLP can predict dipole moments, polarizabilities, and IR-spectra of unseen molecules with high precision. Furthermore, a variety of response properties, for example, stiffness or piezoelectric constants, can be accurately reproduced.

Bending the Curve: Molecular Manifestations of Electron Antisymmetry

As very light fermions, electrons are governed by antisymmetric wave functions that lead to exchange integrals in the evaluation of the energy. Here we use the localized representation of orbitals to decompose the electronic energy in a fashion that isolates the enigmatic exchange contributions and characterizes their distinctive control over electron distributions. The key to this completely general analysis is considering the electrons in groups of three, drawing attention to the curvatures of pair potentials, rather than just their amplitudes and slopes. We show that a positive curvature at short distances is essential for the mutual distancing of electrons and a negative curvature at longer distances is essential to account for the influence of lone pairs on bond torsion. Neither curvature is available in the absence of the exchange contributions. Thus, although exchange energies are much shorter range than Coulomb energies, their influence on molecular geometry is profound and readily understood.

Exchange potentials for semi-classical electrons

Semi-classical electrons offer access to efficient and intuitive simulations of chemical reactions. As for any treatment of fermions, the greatest difficulty is in accounting for anti-symmetry effects. Semi-classical efforts to-date either reference Slater-determinants from ab initio treatments or adopt a heuristic approach inspired by density functional treatments. Here we revisit the problem with a combined approach. We conclude that semi-classical electrons need to reference a non-conventional wave function and find that (1) contrary to earlier suppositions, contributions from the electrostatic terms in the Hamiltonian are of similar magnitude to those from the kinetic terms and (2) the former point to a need to supplement pair potentials with 3-body potentials. The first result explains features of reported heuristic potentials, and the second provides a firm footing for extending the transferability of potentials across a wider range of elements and bonding scenarios.

SOMO-HOMO Level Inversion in Biologically Important Radicals

Conventionally, the singly occupied molecular orbital (SOMO) of a radical species is considered to be the highest occupied molecular orbital (HOMO), but this is not the case always. In this study, we considered a number of radicals from smallest diatomic anion radicals such as superoxide anion radical to one-electron oxidized DNA related base radicals that show the SOMO is energetically lower than one or more doubly occupied molecular orbitals (MOs) (SOMO–HOMO level inversion). The electronic configurations are calculated employing the B3LYP/6-31++G** method, with the inclusion of aqueous phase via the integral equation formalism of the polarized continuum model solvation model. From the extensive study of the electronic configurations of radicals produced by one-electron oxidation or reduction of natural-DNA bases, bromine-, sulfur-, selenium-, and aza-substituted DNA bases, as well as 20 diatomic molecules, we highlight the following important findings: (i) SOMO–HOMO level inversion is a common phenomenon in radical species. (ii) The more localized spin density in σ-orbital on a single atom (carbon, nitrogen, oxygen, sulfur, or selenium), the greater the gap between HOMO and SOMO. (iii) In species with SOMO–HOMO level inversion, one-electron oxidation takes place from HOMO not from the SOMO, which produces a molecule in its triplet ground state. Oxidation of aqueous superoxide anion producing triplet molecular oxygen is one example of many. (iv) These results are for conventional radicals and in contrast with those reported for distonic radical anions in which SOMO–HOMO gaps are smaller for more localized radicals and the orbital inversions vanish in water. Our findings yield new insights into the properties of free radical systems.

Chemistry with semi-classical electrons: reaction trajectories auto-generated by sub-atomistic force fields

For a century now, "Lewis dots" have been a mainstay of chemical thinking, teaching and communication. However, chemists have assumed that this semi-classical picture of electrons needs to be abandoned for quantitative work, and the recourse in computational simulations has been to the extremes of first principles treatments of electrons on the one hand and force fields that avoid explicit electrons on the other hand. Given both the successes and limitations of these highly divergent approaches, it seems worth considering whether the Lewis dot picture might be made quantitative after all. Here we review progress to that end, including variations that have been implemented and examples of applications, specifically the acid-base behavior of water, several organic reactions, and electron dynamics in silicon fracture. In each case, the semi-classical approach is highly efficient and generates reasonable and readily interpreted reaction trajectories in turnkey fashion (i.e., without any input about products). Avenues for further progress are also discussed.