From Open‐Shell Singlet Diradicaloid to Closed‐Shell Global Antiaromatic Macrocycles

A dithieno[a,h]‐s‐indacene‐ (DTI‐) based diradicaloid DTI‐2Br was synthesized and its open‐shell singlet diradical character was validated by magnetic measurements. On the other hand, its macrocyclic trimer DTI‐MC3 and tetramer DTI‐MC4 turned out to be closed‐shell compounds with global antiaromaticity, which was supported by X‐ray crystallographic analysis and NMR spectroscopy, assisted by ACID and 2D‐ICSS calculations. Such change can be explained by a subtle balance between two types of antiferromagnetic spin–spin coupling along the π‐conjugated macrocycles. The dications of DTI‐MC3 and DTI‐MC4 turned out to be open‐shell singlet diradical dications, with a singlet–triplet energy gap of ?2.90 and ?2.60 kcal mol^?1, respectively. At the same time, they are both global aromatic. Our studies show that intramolecular spin–spin interactions play important roles on electronic properties of π‐conjugated macrocycles.     

Revitalizing Spin Natural Orbital Analysis: Electronic Structures of Mixed-Valence Compounds,Singlet Biradicals,and Antiferromagnetically Coupled Systems

Chemical systems with open-shell electronic structure have been gaining attention these days. Their potential applications in first-row transition metal catalysis, molecular wires, photovoltaics and other potential applications have urged the adoption of a simple analysis tool to better understand their open-shell electronic structures, especially the role played by the unpaired electrons. Despite its lack of popularity, spin natural orbital (SNO) analysis is a tool we found to well-suit this purpose. We have therefore re-examined how the SNO could help us analyze some interesting open-shell systems, including mixed-valence compounds, singlet biradicals, and antiferromagnetically coupled systems. We found that some interesting patterns emerge from SNO analysis, especially those associated with exchange interaction. © 2019 Wiley Periodicals, Inc.     

Electronic Stabilization Effect of a Spin‐Delocalized Neutral Radical: Synthesis of an 8‐Cyano‐6‐oxophenalenoxyl Derivative and Quantitative Evaluation of the Electronic Spin Structure in terms of Resonance Structures

A new 2,5‐di‐tert‐butyl‐6‐oxophenalenoxyl (6OPO) derivative with a cyano group at the 8‐position, where a large spin density resides, has been synthesized. This neutral radical exhibits high stability in the solid state in air despite the low steric protection on the 8‐position; the stability is comparable to that of a corresponding 8‐tert‐butylated 6OPO derivative. EPR/¹H‐ENDOR/TRIPLE (electron paramagnetic resonance/¹H‐electron‐nuclear double resonance/TRIPLE) spectroscopy and cyclic voltammetry showed an extended spin delocalization on the cyano group and a significant increase in electron‐accepting ability relative to that of the 8‐tert‐butylated 6OPO derivative. DFT calculations indicated the extension of a singly occupied molecular orbital (SOMO) onto the cyano group and the lower‐lying SOMO and LUMO in comparison with those of the 8‐tert‐butylated 6OPO derivative, which was consistent with experimental results. Furthermore, the extended nature of π conjugation onto the cyano group was quantitatively evaluated by calculating the contributing weights of resonance structures in terms of a molecular orbital (MO)‐based valence‐bond (VB) method. Herein, the synthesis and physical properties of the 8‐cyano‐6OPO derivative are described, emphasizing that the high stability arises from the electronic effect of the cyano group. Also, the usefulness of the quantitative resonance structure analysis is shown.     

Charge‐Shift Bonding: A New and Unique Form of Bonding

Charge‐shift bonds (CSBs) constitute a new class of bonds different than covalent/polar‐covalent and ionic bonds. Bonding in CSBs does not arise from either the covalent or the ionic structures of the bond, but rather from the resonance interaction between the structures. This Essay describes the reasons why the CSB family was overlooked by valence‐bond pioneers and then demonstrates that the unique status of CSBs is not theory‐dependent. Thus, valence bond (VB), molecular orbital (MO), and energy decomposition analysis (EDA), as well as a variety of electron density theories all show the distinction of CSBs vis‐à‐vis covalent and ionic bonds. Furthermore, the covalent–ionic resonance energy can be quantified from experiment, and hence has the same essential status as resonance energies of organic molecules, e.g., benzene. The Essay ends by arguing that CSBs are a distinct family of bonding, with a potential to bring about a Renaissance in the mental map of the chemical bond, and to contribute to productive chemical diversity.     

NMR spin–spin coupling constants in hydrogen‐bonded glycine clusters

The influence of the hydrogen bond formation on the NMR spin–spin coupling constants (SSCC), including the Fermi contact (FC), the diamagnetic spin‐orbit, the paramagnetic spin‐orbit, and the spin dipole term, has been investigated systematically for the homogeneous glycine cluster, in gas phase, containing up to three monomers. The one‐bond and two‐bond SSCCs for several intramolecular (through covalent bond) and intermolecular (across the hydrogen‐bond) atomic pairs are calculated employing the density functional theory with B3LYP and KT3 functionals and different types of extended basis sets. The ab initio SOPPA(CCSD) is used as benchmark for the SSCCs of the glycine monomer. The hydrogen bonding is found to cause significant variations in the one‐bond SSCCs, mostly due to contribution from electronic interactions. However, the nature of variation depends on the type of oxygen atom (proton‐acceptor or proton‐donor) present in the interaction. Two‐bond intermolecular coupling constants vary more than the corresponding one‐bond constants when the size of the cluster increases. Among the four Ramsey terms that constitute the total SSCC, the FC term is the most dominant contributor followed by the paramagnetic spin‐orbit term in all one‐bond interaction.     

Spin uncoupling in ethylene activation by palladium and platinum atoms

Spin‐dependent effects in complex formation reactions of the ethylene molecule with palladium and platinum atoms were studied by electron correlation calculations with account of spin–orbit coupling. Simple correlation diagrams illustrating spin‐uncoupling mechanisms were obtained, showing that the low spin state of the transition‐metal atom or the transition‐metal atom complex is always more reactive than are the high spin states because of the involvement of the triplet excited molecule in the chemical activation. Spin–orbit coupling calculations of the reaction between a platinum atom and ethylene explain the high‐spin Pt(³D) reactivity as due to an effective spin flip at the stage of the weak triplet complex formation. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 581–596, 1999     

Spin–orbit coupling of spin‐frustrated systems

We have investigated the effects of spin–orbit (SO) interactions on noncollinear molecular magnetism by combining the classical Dzyaloshinsky–Moriya (DM) model and ab initio generalized spin orbital (GSO) method. We have derived an estimation scheme of the magnetic anisotropy energy (MAE) and the Dzyaloshinsky vector based on the SO first‐order perturbation theory (SOPT1) for GSO Hartree–Fock (GHF) solutions. We found that the fundamental results of GHF‐SOPT1 method can be reproduced by diagonalizing the core Hamiltonian plus SO terms, and that the spin topologies of odd‐ring systems can be determined by the topological indices of the singly occupied molecular orbitals. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Magnitude of Finite‐Nucleus‐Size Effects in Relativistic Density Functional Computations of Indirect NMR Nuclear Spin–Spin Coupling Constants

A spherical Gaussian nuclear charge distribution model has been implemented for spin‐free (scalar) and two‐component (spin–orbit) relativistic density functional calculations of indirect NMR nuclear spin–spin coupling (J‐coupling) constants. The finite nuclear volume effects on the hyperfine integrals are quite pronounced and as a consequence they noticeably alter coupling constants involving heavy NMR nuclei such as W, Pt, Hg, Tl, and Pb. Typically, the isotropic J‐couplings are reduced in magnitude by about 10 to 15 % for couplings between one of the heaviest NMR nuclei and a light atomic ligand, and even more so for couplings between two heavy atoms. For a subset of the systems studied, viz. the Hg atom, Hg₂²⁺, and Tl? X where X=Br, I, the basis set convergence of the hyperfine integrals and the coupling constants was monitored. For the Hg atom, numerical and basis set calculations of the electron density and the 1s and 6s orbital hyperfine integrals are directly compared. The coupling anisotropies of TlBr and TlI increase by about 2 % due to finite‐nucleus effects.     

Generalized population analysis of three‐center two‐electron bonding

Modern valence bond theory, in its spin‐coupled (SC) form, is used to assess the utility of generalized population analysis, using correlated density matrices, for identifying patterns of three‐center two‐electron bonding in B₂H₆, B₃H, and CH₂Li₂. The effect on the relevant three‐center indices of including electron correlation, via the SC method, is assessed by means of comparison with analogous Hartree–Fock values. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

The Local Kinetic Energy Profile of an Inverted CarbonCarbon Bond Reveals and Refines its Charge‐Shift Character

Analysis of the kinetic energy density within a molecule identifies patterns in its electronic structure that are linked to the concept of charge‐shift bonding. This is illustrated in a detailed study of twelve molecules, possessing carbon‐carbon covalent as well as carbon‐carbon charge‐shift bonds in various degrees of orders, including propellanes and heteropropellanes. Regions of slow electrons are fundamental for such a correlation, and a RoSE (region of slow electrons) indicator ν_±, based on the positive definite kinetic energy density τ, is employed to characterize classes of charge‐shift bonds in terms of its full topology of all critical points of rank three.     

Ab initio calculation of the Dzyaloshinskii–Moriya parameters: Spin–orbit GSO‐HF,DFT, and CI approaches

We present ab initio methods to determine the Dzyaloshinskii–Moriya (DM) parameter, which provides the anisotropic effects of noncollinear spin systems. For this purpose, we explore various general spin orbital (GSO) approaches, such as Hartree–Fock (HF), density functional theory (DFT), and configuration interaction (CI), with one‐electron spin–orbit coupling (SOC1). As examples, two simple D_3h‐symmetric models, H₃ and B(CH₂)₃, are examined. Implications of the computational results are discussed in relation to as isotropic and anisotropic interactions of molecular‐based magnets. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A theoretical structural analysis of the factors that affect 1JNH, 1hJNH and 2hJNN in N–H···N hydrogen‐bonded complexes

Calculations of ¹ J_NH, ^1h J_NH and ^2h J_NN spin–spin coupling constants of 27 complexes presenting N–H·N hydrogen bonds have allowed to analyze these through hydrogen‐bond coupling as a function of the hybridization of both nitrogen atoms and the charge (+1, 0, ? 1) of the complex. The main conclusions are that the hybridization of N atom of the hydrogen bond donor is much more important than that of the hydrogen bond acceptor. Positive and negative charges (cationic and anionic complexes) exert opposite effects while the effect of the transition states ‘proton‐in‐the‐middle’ is considerable. Copyright © 2008 John Wiley & Sons, Ltd.     

Organic–inorganic hybrid materials. I. Characterization and degradation of poly(imide–silica) hybrids

Poly(imide–silica) hybrid materials with covalent bonds were prepared by (3-aminopropyl)methyldiethoxysilane (APrMDEOS) terminated amic acid, water, and tetramethoxysilane (TMOS) via a sol–gel technique. Infrared (IR), ²⁹Si and ¹³C CP/MAS nuclear magnetic resonance (NMR) spectroscopy, and thermogravimetric analysis (TGA) were used to study hybrids containing various proportions of TMOS and hydrolysis ratios. The microstructure and chain mobility of hybrids were investigated by proton spin–spin relaxation T₂ measurements. The apparent activation energy E_a for degradation of hybrids in air was studied by the van Krevelen method. © 1999 John Wiley & Sons, Inc. J Polym Sci A: Polym Chem 37: 2275–2284, 1999     

Theoretical study of bifurcated hydrogen bonding effects on the 1J(N,H), 1hJ(N,H), 2hJ(N,N) couplings and 1H, 15N shieldings in model pyrroles

According to the density functional theory calculations, the X···H···N (X?N, O) intramolecular bifurcated (three‐centered) hydrogen bond with one hydrogen donor and two hydrogen acceptors causes a significant decrease of the ^1hJ(N,H) and ^2hJ(N,N) coupling constants across the N? H···N hydrogen bond and an increase of the ¹J(N,H) coupling constant across the N? H covalent bond in the 2,5‐disubsituted pyrroles. This occurs due to a weakening of the N? H···N hydrogen bridge resulting in a lengthening of the N···H distance and a decrease of the hydrogen bond angle at the bifurcated hydrogen bond formation. The gauge‐independent atomic orbital calculations of the shielding constants suggest that a weakening of the N? H···N hydrogen bridge in case of the three‐centered hydrogen bond yields a shielding of the bridge proton and deshielding of the acceptor nitrogen atom. The atoms‐in‐molecules analysis shows that an attenuation of the ^1hJ(N,H) and ^2hJ(N,N) couplings in the compounds with bifurcated hydrogen bond is connected with a decrease of the electron density ρ_H···N at the hydrogen bond critical point and Laplacian of this electron density ?²ρ_H···N. The natural bond orbital analysis suggests that the additional N? H···X interaction partly inhibits the charge transfer from the nitrogen lone pair to the σ*_{N? H} antibonding orbital across hydrogen bond weakening of the ^1hJ(N,H) and ^2hJ(N,N) trans‐hydrogen bond couplings through Fermi‐contact mechanism. An increase of the nitrogen s‐character percentage of the N? H bond in consequence of the bifurcated hydrogen bonding leads to an increase of the ¹J(N,H) coupling constant across the N? H covalent bond and deshielding of the hydrogen donor nitrogen atom. Copyright © 2010 John Wiley & Sons, Ltd.     

Singlet and Triplet Excited States and Intersystem Crossing in Free‐Base Porphyrin: TDDFT and DFT/MRCI Study

Extensive time-dependent DFT (TDDFT) and DFT/multireference configuration interaction (MRCI) calculations are performed on the singlet and triplet excited states of free-base porphyrin, with emphasis on intersystem crossing processes. The equilibrium geometries, as well as the vertical and adiabatic excitation energies of the lowest singlet and triplet excited states are determined. Single and double proton-transfer reactions in the first excited singlet state are explored. Harmonic vibrational frequencies are calculated at the equilibrium geometries of the ground state and of the lowest singlet and triplet excited states. Furthermore, spin–orbit coupling matrix elements of the lowest singlet and triplet states and their numerical derivatives with respect to nuclear displacements are computed. It is shown that opening of an unprotonated pyrrole ring as well as excited-state single and double proton transfer inside the porphyrin cavity lead to crossings of the potential energy curves of the lowest singlet and triplet excited states. It is also found that displacements along out-of-plane normal modes of the first excited singlet state cause a significant increase of the 2|H_so|S₁>, 1|H_so|S₁>, and 1|H_so|S₀> spin–orbit coupling matrix elements. These phenomena lead to efficient radiationless deactivation of the lowest excited states of free-base porphyrin via intercombination conversion. In particular, the S₁→T₁ population transfer is found to proceed at a rate of ≈10⁷ s⁻¹ in the isolated molecule.     

The Nature of the Fourth Bond in the Ground State of C2: The Quadruple Bond Conundrum

Does, or doesn’t C₂ break the glass ceiling of triple bonding? This work provides an overview on the bonding conundrum in C₂ and on the recent discussions regarding our proposal that it possesses a quadruple bond. As such, we focus herein on the main point of contention, the 4th bond of C₂, and discuss the main views. We present new data and an overview of the nature of the 4th bond—its proposed antiferromagnetically coupled nature, its strength, and a derivation of its bond energy from experimentally based thermochemical data. We address the bond‐order conundrum of C₂ arising from generalized VB (GVB) calculations by comparing it to HC?CH, and showing that the two molecules behave very similarly, and C₂ is in no way an exception. We analyse the root cause of the deviation of C₂ from the Badger Rule, and demonstrate that the reason for the smaller force constant (FC) of C₂ relative to HC?CH has nothing to do with the bond energies, or with the number of bonds in the two molecules. The FC is determined primarily by the bond length, which is set by the balance between the bond length preferences of the σ‐ versus π‐bonds in the two molecules. This interplay in the case of C₂ clearly shows the fingerprints of the 4th bond. Our discussion resolves the points of contention and shows that the arguments used to dismiss the quadruple bond nature of C₂ are not well founded.     

A Serendipitous Rendezvous with a Four‐Center Two‐Electron Bonded Intermediate in the Aerial Oxidation of Hydrazine

Oxidation by dioxygen has a rich repertoire of mechanistic intricacies. Herein, we report a hitherto unknown paradigm of dioxygen activation reaction which propagates through a four center two electron (4c–2e) bound species. Using static DFT and ab initio quantum chemical techniques we have unraveled the oxidation pathway for hydrazine and its methylated analogues by dioxygen which involves formation of this unconventional 4c–2e bonded species en route to the oxidation products. Inconvertible evidence in favor of such an unprecedented dioxygen activation route is provided by capturing the events of formation of the 4c–2e species in aqueous phase for hydrazine and its congeners and the experimentally observed products from the respective 4c–2e species, like H₂O₂ and N₂H₂, diazene in the case of hydrazine using Car–Parrinello molecular dynamics simulations.