The electronic structure of benzene from a tiling of the correlated 126-dimensional wavefunction

A Unique Hybridization Route to Access Hydrazylnaphthalimidols as Novel Structural Scaffolds of Multitargeting Broad-Spectrum Antifungal Candidates

This study developed a class of novel structural antifungal hydrazylnaphthalimidols (HNs) with multitargeting broad-spectrum potential via multicomponent hybridization to confront increasingly severe fungal invasion. Some prepared HNs exhibited considerable antifungal potency; especially nitrofuryl HN 4a (MIC = 0.001 mM) exhibited a potent antifungal activity against Candida albicans, which is 13-fold higher than that of fluconazole. Furthermore, nitrofuryl HN 4a displayed low cytotoxicity, hemolysis and resistance, as well as a rapid fungicidal efficacy. Preliminary mechanistic investigations revealed that nitrofuryl HN 4a could inhibit lactate dehydrogenase to decrease metabolic activity and promote the accumulation of reactive oxygen species, leading to oxidative stress. Moreover, nitrofuryl HN 4a did not exhibit membrane-targeting ability; it could embed into DNA to block DNA replication but could not cleave DNA. These findings implied that HNs are promising as novel structural scaffolds of potential multitargeting broad-spectrum antifungal candidates for treating fungal infection.

Persistence of atoms in molecules: there is room beyond electron densities

Evidence that the electronic structure of atoms persists in molecules to a much greater extent than has been usually admitted is presented. This is achieved by resorting to N-electron real-space descriptors instead of one- or at most two-particle projections like the electron or exchange-correlation densities. Here, the 3N-dimensional maxima of the square of the wavefunction, the so-called Born maxima, are used. Since this technique is relatively unknown to the crystallographic community, a case-based approach is taken, revisiting first the Born maxima of atoms in their ground state and then some of their excited states. It is shown how they survive in molecules and that, beyond any doubt, the distribution of electrons around an atom in a molecule can be recognized as that of its isolated, in many cases excited, counterpart, relating this fact with the concept of energetic promotion. Several other cases that exemplify the applicability of the technique to solve chemical bonding conflicts and to introduce predictability in real-space analyses are also examined.

Separation of benzene and toluene associated with vapochromic behaviors by hybrid[4]arene-based co-crystals

The combination of macrocyclic chemistry with co-crystal engineering has promoted the development of materials with vapochromic behaviors in supramolecular science. Herein, we develop a macrocycle co-crystal based on hybrid[4]arene and 1,2,4,5-tetracyanobenzene that is able to construct vapochromic materials. After the capture of benzene and toluene vapors, activated hybrid[4]arene-based co-crystal forms new structures, accompanied by color changes from brown to yellow. However, when hybrid[4]arene-based co-crystal captures cyclohexane and pyridine, neither structures nor colors change. Interestingly, hybrid[4]arene-based co-crystal can separate benzene from a benzene/cyclohexane equal-volume mixture and allow toluene to be removed from a toluene/ pyridine equal-volume mixture with purities reaching 100%. In addition, the process of adsorptive separation can be visually monitored. The selectivity of benzene from a benzene/cyclohexane equal-volume mixture and toluene from a toluene/ pyridine equal-volume mixture is attributed to the different changes in the charge-transfer interaction between hybrid[4]arene and 1,2,4,5-tetracyanobenzene when hybrid[4]arene-based co-crystal captures different vapors. Moreover, hybrid[4]arene-based co-crystal can be reused without losing selectivity and performance. This work constructs a vapochromic material for hydrocarbon separation.

Towards the ground state of molecules via diffusion Monte Carlo on neural networks

Diffusion Monte Carlo (DMC) based on fixed-node approximation has enjoyed significant developments in the past decades and become one of the go-to methods when accurate ground state energy of molecules and materials is needed. However, the inaccurate nodal structure hinders the application of DMC for more challenging electronic correlation problems. In this work, we apply the neural-network based trial wavefunction in fixed-node DMC, which allows accurate calculations of a broad range of atomic and molecular systems of different electronic characteristics. Our method is superior in both accuracy and efficiency compared to state-of-the-art neural network methods using variational Monte Carlo (VMC). We also introduce an extrapolation scheme based on the empirical linearity between VMC and DMC energies, and significantly improve our binding energy calculation. Overall, this computational framework provides a benchmark for accurate solutions of correlated electronic wavefunction and also sheds light on the chemical understanding of molecules.

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.

Atomic shell structure from Born probabilities: Comparison to other shell descriptors and persistence in molecules

Real space chemical bonding descriptors, such as the electron localization function or the Laplacian of the electron density, have been widely used in electronic structure theory thanks to their power to provide chemically intuitive spatial images of bonded and non-bonded interactions. This capacity stems from their ability to display the shell structure of atoms and its distortion upon molecular formation. Here, we examine the spatial position of the N electrons of an atom at the maximum of the square of the wavefunction, the so-called Born maximum, as a shell structure descriptor for ground state atoms with Z = 1-36, comparing it to other available indices. The maximization is performed with the help of variational quantum Monte Carlo calculations. We show that many electron effects (mainly Pauli driven) are non-negligible, that Born shells are closer to the nucleus than any other of the examined descriptors, and that these shells are very well preserved in simple molecules.

The Hitchhiker's Guide to the Wave Function∥

The electronic wave function of molecules is 3N-dimensional and inseparable in the coordinates of the N electrons. Whereas molecular orbitals are often invoked to visualize the electronic structure, they are nonunique, with the same 3N-dimensional wave function being represented by an infinite number of 3-D, one-electron functions (orbitals). Furthermore, multireference wave functions cannot be described by an antisymmetrized product of a single set of occupied orbitals. What is required is a way to visualize the full dimensionality of the wave function, including the effects of correlation, as a 3N-dimensional being would be able to do. In the past 5 years, we have been developing a way to analyze and visualize highly dimensional wave functions by focusing on the structure of the repeating unit demanded by fermionic behavior. This 3N-dimensional repeating unit, the wave function "tile", can be projected onto the three dimensions of each electron, in turn, to reveal the complete electronic structure. It is found that the tile reproduces canonical chemical motifs such as core-electrons, single bonds and lone pairs. Multiple bonds emerge as the "banana" bonds favored by Pauling. As a function of the reaction coordinate, electron motions are visualized that correspond to the curly arrow notation of organic chemists. Excited states can also be inspected. Analyzing a wave function in terms of fermionic tiling allows for insight not facilitated by the inspection of orbitals or configuration interaction vectors: The wave function tiles of resonance structures reveal that electron correlation in benzene pushes opposing spin electrons to occupy alternate Kekulé structures, and in C₂, the emerging structure supports the notion of a triply bonded structure with a weak, fourth bonding contribution.

Chemical bonding theories as guides for self-interaction corrected solutions: Multiple local minima and symmetry breaking

Fermi-Löwdin orbitals (FLOs) are a special set of localized orbitals, which have become commonly used in combination with the Perdew-Zunger self-interaction correction (SIC) in the FLO-SIC method. The FLOs are obtained for a set of occupied orbitals by specifying a classical position for each electron. These positions are known as Fermi-orbital descriptors (FODs), and they have a clear relation to chemical bonding. In this study, we show how FLOs and FODs can be used to initialize, interpret, and justify SIC solutions in a common chemical picture, both within FLO-SIC and in traditional variational SIC, and to locate distinct local minima in either of these approaches. We demonstrate that FLOs based on Lewis theory lead to symmetry breaking for benzene-the electron density is found to break symmetry already at the symmetric molecular structure-while ones from Linnett's double-quartet theory reproduce symmetric electron densities and molecular geometries. Introducing a benchmark set of 16 planar cyclic molecules, we show that using Lewis theory as the starting point can lead to artifactual dipole moments of up to 1 D, while Linnett SIC dipole moments are in better agreement with experimental values. We suggest using the dipole moment as a diagnostic of symmetry breaking in SIC and monitoring it in all SIC calculations. We show that Linnett structures can often be seen as superpositions of Lewis structures and propose Linnett structures as a simple way to describe aromatic systems in SIC with reduced symmetry breaking. The role of hovering FODs is also briefly discussed.

Fluoromaticity: The Molecular Orbital Contributions of Fluorine Substituents to the π-Systems of Aromatic Rings

The addition of fluorine atoms to an aromatic ring brings about an additional set of π-bonding and antibonding orbitals culminating after the addition of the sixth fluorine with a new set of π-aromatic-like orbitals that affect the molecule in a way that we will refer to hereafter as "fluoromaticity". Depending on the number and position of the fluorine atoms, the contributed π-orbitals can even further stabilize the ring leading to smaller bond lengths within the ring and higher resistance to addition reactions. This added ring stability partially explains the high thermostability and chemical resistance found in polymers containing fluorinated aromatics in their architecture. A similar molecular orbital effect is seen with the addition of other halogen atoms to aromatic rings, though to a much smaller degree and not resulting in the additional ring stability.

Real space electron delocalization, resonance, and aromaticity in chemistry

Chemists explaining a molecule’s stability and reactivity often refer to the concepts of delocalization, resonance, and aromaticity. Resonance is commonly discussed within valence bond theory as the stabilizing effect of mixing different Lewis structures. Yet, most computational chemists work with delocalized molecular orbitals, which are also usually employed to explain the concept of aromaticity, a ring delocalization in cyclic planar systems which abide certain number rules. However, all three concepts lack a real space definition, that is not reliant on orbitals or specific wave function expansions. Here, we outline a redefinition from first principles: delocalization means that likely electron arrangements are connected via paths of high probability density in the many-electron real space. In this picture, resonance is the consideration of additional electron arrangements, which offer alternative paths. Most notably, the famous 4n + 2 Hückel rule is generalized and derived from nothing but the antisymmetry of fermionic wave functions.

The concept of delocalization, resonance and aromaticity are commonly discussed within electronic structure frameworks relying on specific wave function expansions. Here the authors propose a redefinition of these concepts from first-principles by investigating saddle points of the all-electron probability density.

The Shape of Full Configuration Interaction to Come

We present a Perspective on what the future holds for full configuration interaction (FCI) theory, with an emphasis on conceptual rather than technical details. Upon revisiting the early history of FCI, a number of its key contemporary approximations are compared on as equal a footing as possible, using a recent blind challenge on the benzene molecule as a testbed [Eriksen et al., J. Phys. Chem. Lett., 2020 11, 8922]. In the process, we review the scope of applications for which FCI continues to prove indispensable, and the required traits in terms of robustness, efficacy, and reliability its modern approximations must satisfy are discussed. We close by conveying a number of general observations on the merits offered by the state-of-the-art alongside some of the challenges still faced to this day. While the field has altogether seen immense progress over the years-the past decade, in particular-it remains clear that our community as a whole has a substantial way to go in enhancing the overall applicability of near-exact electronic structure theory for systems of general composition and increasing size.

Synchronizing chemistry, quantum mechanics and radioactivity in a revolutionary renewed atom model. Part 1: the elements where Z is 1–10 (H, He, Li, Be, B, C, N, O, F, Ne)

The alliance between the reigning quantum mechanical atom model and chemistry still is a difficult one when it comes to an adequate explanation for e.g. the covalent bond, inversion, chirality, or hydrogen bonds. Overruling Rutherford's extrapolation from gold to hydrogen, an atom model is described that provides improved answers to these phenomena while the hybridization principle and the covalent bond are re-defined by giving neutrons a much more prominent role than they have in the reigning quantum mechanical model. It is postulated that a neutron is not just there to assist the strong force in surpassing the repulsive coulombic forces between the protons in the nucleus, but the neutron is the modus operandi of molecular geometry, and as such plays a part in chemical reactivity, bond length and bond strength.

A renewed atom model is described that provides improved answers to unsolved phenomena like inversion, chirality, hydrogen bonds, the hybridization principle and covalent bonds by giving neutrons a much more prominent role.

Hybrid Boron-Carbon Chemistry

The recently proved one-to-one structural equivalence between a conjugated hydrocarbon CnHm and the corresponding borane BnHm+n is applied here to hybrid systems, where each C=C double bond in the hydrocarbon is consecutively substituted by planar B(H2)B moieties from diborane(6). Quantum chemical computations with the B3LYP/cc-pVTZ method show that the structural equivalences are maintained along the substitutions, even for non-planar systems. We use as benchmark aromatic and antiaromatic (poly)cyclic conjugated hydrocarbons: cyclobutadiene, benzene, cyclooctatetraene, pentalene, benzocyclobutadiene, naphthalene and azulene. The transformation of these conjugated hydrocarbons to the corresponding boranes is analyzed from the viewpoint of geometry and electronic structure.

The Ground State Electronic Energy of Benzene

We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground-state energy of the benzene molecule in a standard correlation-consistent basis set of double-ζ quality. As a broad international endeavor, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around -863 mE_H. However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 mE_H), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward.

On the connection between probability density analysis, QTAIM, and VB theory

Classification of bonds is essential for understanding and predicting the reactivity of chemical compounds. This classification mainly manifests in the bond order and the contribution of different Lewis resonance structures. Here, we outline a first principles approach to obtain these orders and contributions for arbitrary wave functions in a manner that is both, related to the quantum theory of atoms in molecules and consistent with valence bond theory insight: the Lewis structures arise naturally as attractors of the all-electron probability density |Ψ|2. Doing so, we introduce a valence bond weight definition that does not collapse in the basis set limit.