Complexes with a Single Metal–Metal Bond as a Sensitive Probe of Quality of Exchange-Correlation Functionals

A comprehensive electron wavefunction analysis toolbox for chemists, Multiwfn

Analysis of electron wavefunction is a key component of quantum chemistry investigations and is indispensable for the practical research of many chemical problems. After more than ten years of active development, the wavefunction analysis program Multiwfn has accumulated very rich functions, and its application scope has covered numerous aspects of theoretical chemical research, including charge distribution, chemical bond, electron localization and delocalization, aromaticity, intramolecular and intermolecular interactions, electronic excitation, and response property. This article systematically introduces the features and functions of the latest version of Multiwfn and provides many representative examples. Through this article, readers will be able to fully understand the characteristics and recognize the unique value of Multiwfn. The source code and precompiled executable files of Multiwfn, as well as the manual containing a detailed introduction to theoretical backgrounds and very rich tutorials, can all be downloaded for free from the Multiwfn website (http://sobereva.com/multiwfn).

B[double bond, length as m-dash]P double bonds relieved from steric encumbrance: matrix-isolation infrared spectroscopy of the phosphaborene F2B-P[double bond, length as m-dash]BF and the triradical B[double bond, length as m-dash]PF3

Free phosphaborenes have a labile boron-phosphorus double bond and therefore require extensive steric shielding by bulky substituents to prevent isomerisation and oligomerisation. In the present work, the small free phosphaborene F2B-P[double bond, length as m-dash]BF was isolated by matrix-isolation techniques and was characterised by infrared spectroscopy in conjunction with quantum-chemical methods. In contrast to its sterically hindered relatives, this small phosphaborene exhibits an acute BPB angle of 83° at the CCSD(T) level. An alternative orbital structure for the B[double bond, length as m-dash]P double bond is found in the triradical B[double bond, length as m-dash]PF3, the direct adduct of laser-ablated atomic B and PF3. The single-bonded isomer F2B-PF and the dimer F3P-B[triple bond, length as m-dash]B-PF3 are also tentatively assigned.

Oxidative addition of an alkyl halide to form a stable Cu(III) product

The step that cleaves the carbon-halogen bond in copper-catalyzed cross-coupling reactions remains ill defined because of the multiple redox manifolds available to copper and the instability of the high-valent copper product formed. We report the oxidative addition of α-haloacetonitrile to ionic and neutral copper(I) complexes to form previously elusive but here fully characterized copper(III) complexes. The stability of these complexes stems from the strong Cu-CF3 bond and the high barrier for C(CF3)-C(CH2CN) bond-forming reductive elimination. The mechanistic studies we performed suggest that oxidative addition to ionic and neutral copper(I) complexes proceeds by means of two different pathways: an SN2-type substitution to the ionic complex and a halogen-atom transfer to the neutral complex. We observed a pronounced ligand acceleration of the oxidative addition, which correlates with that observed in the copper-catalyzed couplings of azoles, amines, or alkynes with alkyl electrophiles.

Reductive Al-B σ-Bond Formation in Alumaboranes: Facile Scission of Polar Multiple Bonds

We present facile access to an alumaborane species with electron precise Al-B σ-bond. The reductive rearrangement of 1-(AlI2 ), 8-(BMes2 ) naphthalene (Mes=2,4,6-Me3 C6 H2 ) affords the alumaborane species cyclo-(1,8-C10 H6 )-[1-Al(Mes)(OEt2 )-8-B(Mes)] with a covalent Al-B σ-bond. The Al-B σ-bond performs the reductive scission of multiple bonds: S=C(NiPrCMe)2 affords the naphthalene bridged motif B-S-Al(NHC), NHC=N-heterocyclic carbene, while O=CPh2 is deoxygenated to afford an B-O-Al bridged species with incorporation of the remaining ≡CPh2 fragment into the naphthalene scaffold. The reaction with isonitrile Xyl-N≡C (Xyl=2,6-Me2 C6 H4 ) proceeds via a proposed (amino boryl) carbene species; which adds a second equivalent of isonitrile to ultimately form the Al-N-B bridged species cyclo-(1,8-C10 H6 )-[1-Al(Mes)-N(Xyl)-8-B{C(Mes)=C-N-Xyl}] with complete scission of the C≡N triple bond. The latter reaction is supported with isolated intermediates and by DFT calculations.

Structural control in the nanoassembly of the tungsten and molybdenum dithiolene complex analog

A strategy for precisely tuning the self-assembly of tungsten and molybdenum dithiolene complexes to nanoflowers and nanopolyhedra is put forward.

Insights into Molecular Magnetism in Metal-Metal Bonded Systems as Revealed by a Spectroscopic and Computational Analysis of Diiron Complexes

A pair of bimetallic compounds featuring Fe-Fe bonds, [Fe(ⁱPrNPPh₂)₃FeR] (R = PMe₃, ≡N^tBu), have been investigated using High-Frequency Electron Paramagnetic Resonance (HFEPR) as well as field- and temperature-dependent ⁵⁷Fe nuclear γ resonance (Mössbauer) spectroscopy. To gain insight into the local site electronic structure, we have concurrently studied a compound containing a single Fe(II) in a geometry analogous to that of one of the dimer sites. Our spectroscopic studies have allowed for the assessment of the electronic structure via the determination of the zero-field splitting and ⁵⁷Fe hyperfine parameters for the entire series. We also report on our efforts to correlate structure with physical properties in metal-metal bonded systems using ligand field theory guided by quantum chemical calculations. Through the insight gained in this study, we discuss strategies for the design of single-molecule magnets based on polymetallic compounds linked via direct metal-metal bonds.

Hydrogen release from a single water molecule on Vn+ (3 ≤ n ≤ 30)

Water and its interactions with metals are closely bound up with human life, and the reactivity of metal clusters with water is of fundamental importance for the understanding of hydrogen generation. Here a prominent hydrogen evolution reaction (HER) of single water molecule on vanadium clusters V_n⁺ (3 ≤ n ≤ 30) is observed in the reaction of cationic vanadium clusters with water at room temperature. The combined experimental and theoretical studies reveal that the wagging vibrations of a V-OH group give rise to readily formed V-O-V intermediate states on V_n⁺ (n ≥ 3) clusters and allow the terminal hydrogen to interact with an adsorbed hydrogen atom, enabling hydrogen release. The presence of three metal atoms reduces the energy barrier of the rate-determining step, giving rise to an effective production of hydrogen from single water molecules. This mechanism differs from dissociative chemisorption of multiple water molecules on aluminium cluster anions, which usually proceeds by dissociative chemisorption of at least two water molecules at multiple surface sites followed by a recombination of the adsorbed hydrogen atoms.

A density functional theory study on the electronic and adsorption characteristics of cyclo M9N9 (M = B and Al)

On the basis of experimental and theoretical calculations conducted for cyclo C₁₈, we predicted two novel inorganic cyclo M₉N₉ (M = B and Al) molecules. Because of the significant difference in electronegativity between M and N atoms, M-N bonds were ionic. Furthermore, the interaction of cyclo M₉N₉ with cyclopropylpiperazine (CPPP) was investigated. In cyclo M₉N₉, each M atom could adsorb one CPPP molecule. The CPPP molecules exhibited a preference to remain outside cyclo M₉N₉ molecules. Depending on the structural characteristics of CPPP molecules, the exciting part is that up to four CPPP molecules could be adsorbed on the exterior surface of cyclo M₉N₉. We calculated adsorption energies and analyzed the main structural parameters in the process. The research results indicated that adsorption on cyclo Al₉N₉ was energetically more favorable than that on cyclo B₉N₉. The cyclo M₉N₉ have considerable potential in the future. Graphical abstract.

Unveiling the Relationship between Energy Transfer and the Triplet Energy Level by Tuning Diarylethene within Europium(III) Complexes

Luminescence performance and photoisomerization control of sensitized energy transfer in a series of Eu(acac)₃De complexes that contain photochromic diarylethene (De) as the ligand are studied by theoretical methods. Both the open-ring and closed-ring isomers exhibit a consistent coordination mode between the Eu^III ion and De. An unneglected weak interaction originating from electrostatic attraction is found in the region of the coordinate bond Eu-N. The open-ring isomer has higher triplet energy levels than ⁵D₁ and ⁵D₀ of the Eu^III ion, which facilitates forward energy transfer from De to the Eu^III ion. The closed-ring isomer, for the extended conjugated system formed in cyclization, has a much lower triplet energy level than ⁵D₀ of the Eu^III ion. The energy-gap deficit makes energy transfer unavailable. By utilization of this phenomenon, regulation of energy transfer and reversible on/off luminescence switching of the europium(III) complex can be achieved. The forward and backward energy-transfer rates in different channels are also calculated for the series of complexes. A statistics diagram is obtained to exhibit the change trend of energy-transfer rates in the forward and backward directions as a function of the triplet energy level, which indicates the contribution of different channels to energy transfer in each level region and figures out that the optimal triplet energy level should be in the range of 21740-19532 cm^-1.

Importance of the iron-sulfur component and of the siroheme modification in the resting state of sulfite reductase

The active site of sulfite reductase (SiR) consists of an unusual siroheme-Fe₄S₄ assembly coupled via a cysteinate sulfur, and serves for multi-electron reduction reactions. Clear explanations have not been demonstrated for the reasons behind the choice of siroheme (vs. other types of heme) or for the single-atom coupling to an Fe₄S₄ center (as opposed to simple adjacency or to coupling via chains consisting of more than one atom). Possible explanations for these choices have previously been invoked, relating to the control of the spin state of the substrate-binding (siro)heme iron, modulation of the trans effect of the (Fe₄S₄-bound) cysteinate, or modulation of the redox potential. Reported here is a density functional theory (DFT) investigation of the structural interplay (in terms of geometry, molecular orbitals and magnetic interactions) between the siroheme and the Fe₄S₄ center as well as the importance of the covalent modifications within siroheme compared to the more common heme b, aiming to verify the role of the siroheme modification and of the Fe₄S₄ cluster at the SiR active site, with focus on previously-formulated hypotheses (geometrical/sterics, spin state, redox and electron-transfer control). A calibration of various DFT methods/variants for the correct description of ground state spin multiplicity is performed using a set of problematic cases of bioinorganic Fe centers; out of 11 functionals tested, M06-L and B3LYP offer the best results - though none of them correctly predict the spin state for all test cases. Upon examination of the relative energies of spin states, reduction potentials, energy decomposition (electrostatic, exchange-repulsion, orbital relaxation, correlation and dispersion interactions) and Mayer bond indices in SiR models, the following main roles of the siroheme and cubane are identified: (1) the cubane cofactor decreases the reduction potential of the siroheme and stabilizes the siroheme-cysteine bond interaction, and (2) the siroheme removes the quasi-degeneracy between the intermediate and high-spin states found in ferrous systems by preserving the latter as ground state; the higher-spin preference and the increased accessibility of multiple spin states are likely to be important in selective binding of the substrate and of the subsequent reaction intermediates, and in efficient changes in redox states throughout the catalytic cycle.

Dysprosium complexes bearing unsupported DyIII-GeII/SnII metal-metal bonds as single-ion magnets

Two dysprosium complexes bearing unsupported Dy-Ge/Sn metal-metal bonds are reported here, wherein the Dy-Ge and Dy-Sn bonds both contain relatively large covalency. The complexes exhibit slow relaxation of magnetization at zero field with energy barriers of 485 and 620 K, respectively, and the blocking temperature of 6 K.

Electronic Structure of Cubane-Like Vanadium–Nitrogen Cationic Clusters [V4N4]+ and [V6N6]+

Density Functional Theory and Complete Active Space Self-Consistent Field (CASSCF) methodologies are used to explore the electronic structure of the cationic V–N clusters, [V4N4]+ and [V6N6]+, that have been identified in recent mass spectrometric experiments. Our calculations indicate that both clusters are based on cubane-like fragments of the rock-salt lattice. In the smaller [V4N4]+ cluster, the V–V bonding is delocalized over the tetrahedron, with net bond orders of 1/3 per V–V bond. In [V6N6]+, in contrast, the V–V bonding is strongly localized in the central V2N2 unit, which has a short V=V double bond. CASSCF calculations reveal that both localized and delocalized V–V bonds are highly multi-configurational.

Probing Fe-V Bonding in a C3-Symmetric Heterobimetallic Complex

Direct metal-metal bonding of two distinct first-row transition metals remains relatively unexplored compared to their second- and third-row heterobimetallic counterparts. Herein, a recently reported Fe-V triply bonded species, [V( ⁱPrNPPh₂)₃FeI] (1; Kuppuswamy, S.; Powers, T. M.; Krogman, J. P.; Bezpalko, M. W.; Foxman, B. M.; Thomas, C. M. Vanadium-iron complexes featuring metal-metal multiple bonds. Chem. Sci. 2013, 4, 3557-3565), is investigated using high-frequency electron paramagnetic resonance, field- and temperature-dependent ⁵⁷Fe nuclear gamma resonance (Mössbauer) spectroscopy, and high-field electron-electron double resonance detected nuclear magnetic resonance. From the use of this suite of physical methods, we have assessed the electronic structure of 1. These studies allow us to establish the effective g̃ tensors as well as the Fe/V electro-nuclear hyperfine interaction tensors of the spin S = ¹/₂ ground state. We have rationalized these tensors in the context of ligand field theory supported by quantum chemical calculations. This theoretical analysis suggests that the S = ¹/₂ ground state originates from a single unpaired electron predominately localized on the Fe site.

Unsaturated trinuclear iron fluoroborylene complexes

The unsaturated trinuclear iron fluoroborylene complexes Fe₃(BF)₃(CO) _n (n = 7, 6) have been studied using density functional theory (DFT). Relatively complicated potential energy surfaces are found with nine and eight structures within 15 kcal mol^-1 of the lowest energy structures for the Fe₃(BF)₃(CO)₇ and Fe₃(BF)₃(CO)₆ systems, respectively. In each of these low-energy structures all three BF groups are either edge-bridging or face-bridging but never terminal groups. Some, but not all, of the low-energy structures also have edge-bridging and/or face-bridging CO groups leading to some structures with as many as five bridging groups. The relatively narrow range of Fe-Fe distances in the central Fe₃ triangles of the Fe₃(BF)₃(CO) _n (n = 7, 6) structures, mainly between 2.37 and 2.55 Å, suggests considerable delocalization in these unsaturated systems. Graphical Abstract The lowest energy Fe₃(BF)₃(CO)₇ and Fe₃(BF)₃(CO)₆ structures have a face-bridging μ₃-BF group with the two remaining BF groups bridging edges. The lowest energy Fe₃(BF)₃(CO)₆ structure also has one four-electron donor bridging η²-μ-CO group.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Solid-state (185/187)Re NMR and GIPAW DFT study of perrhenates and Re2(CO)10: chemical shift anisotropy, NMR crystallography, and a metal-metal bond

Advances in solid-state nuclear magnetic resonance (SSNMR) methods, such as dynamic nuclear polarization (DNP), intricate pulse sequences, and increased applied magnetic fields, allow for the study of systems which even very recently would be impractical. However, SSNMR methods using certain quadrupolar probe nuclei (i.e., I > 1/2), such as (185/187)Re remain far from fully developed due to the exceedingly strong interaction between the quadrupole moment of these nuclei and local electric field gradients (EFGs). We present a detailed high-field (B0 = 21.1 T) experimental SSNMR study on several perrhenates (KReO4, AgReO4, Ca(ReO4)2·2H2O), as well as ReO3 and Re2(CO)10. We propose solid ReO3 as a new rhenium SSNMR chemical shift standard due to its reproducible and sharp (185/187)Re NMR resonances. We show that for KReO4, previously poorly understood high-order quadrupole-induced effects (HOQIE) on the satellite transitions can be used to measure the EFG tensor asymmetry (i.e., ηQ) to nearly an order-of-magnitude greater precision than competing SSNMR and nuclear quadrupole resonance (NQR) approaches. Samples of AgReO4 and Ca(ReO4)2·2H2O enable us to comment on the effects of counter-ions and hydration upon Re(vii) chemical shifts. Calcium-43 and (185/187)Re NMR tensor parameters allow us to conclude that two proposed crystal structures for Ca(ReO4)2·2H2O, which would be considered as distinct, are in fact the same structure. Study of Re2(CO)10 provides insights into the effects of Re-Re bonding on the rhenium NMR tensor parameters and rhenium oxidation state on the Re chemical shift value. As overtone NQR experiments allowed us to precisely measure the (185/187)Re EFG tensor of Re2(CO)10, we were able to measure rhenium chemical shift anisotropy (CSA) for the first time in a powdered sample. Experimental observations are supported by gauge-including projector augmented-wave (GIPAW) density functional theory (DFT) calculations, with NMR tensor calculations also provided for NH4ReO4, NaReO4 and RbReO4. These calculations are able to reproduce many of the experimental trends in rhenium δiso values and EFG tensor magnitudes. Using KReO4 as a prototypical perrhenate-containing system, we establish a correlation between the tetrahedral shear strain parameter (|ψ|) and the nuclear electric quadrupolar coupling constant (CQ), which enables the refinement of the structure of ND4ReO4. Shortcomings in traditional DFT approaches, even when including relativistic effects via the zeroth-order regular approximation (ZORA), for calculating rhenium NMR tensor parameters are identified for Re2(CO)10.

Fluoride-free Hiyama coupling by palladium abnormal N-heterocyclic carbene complexes

A series of Pd complexes of the 1,2,3-triazole based abnormal NHC ligands of the type (a-NHC)PdI₂(L) [L = NC₅H₅ and PPh₃] successfully catalyzed the fluoride-free Hiyama coupling in air.