Iron bispentazole Fe(eta5-N5)2, a theoretically predicted high-energy compound: structure, bonding analysis, metal-ligand bond strength and a comparison with the isoelectronic ferrocene

Chemical Bonding in [Fe(η4-P4)2]2- and Related Complexes

Quantum chemical calculations of the six valence isoelectronic complexes [FeL2]2-, [CoL2]-, and NiL2 with L = η4-P4, η4-C4H4 using density functional theory have been carried out. The molecular structures were investigated with a variety of methods. The analysis of the electronic structure in [Fe(η4-P4)2]2- shows that the bonding situation is very similar to valence isoelectronic Ni(η4-C4H4)2. The orbital interactions in the 18 electron complexes [TML2]q (TMq = Fe2-, Co-, Ni) come mainly from TM(dπ)→L2 backdonation, enhanced by smaller contributions from TM(dδ)→L2 backdonation and TM(s)←L2 donation. Calculations of the six TML2 species indicate that all of them are viable candidates for synthetic work. The bonding situation is very similar and can straightforwardly be explained with the Dewar-Chatt-Duncanson bonding model in terms of dative bonding between d10 metal atoms and the ligands in the electronic singlet state. EDA-NOCV calculations using the ligands and the metal atoms with different charges and electronic states indicate that the metal-ligand bonds in the charged complexes [FeL2]2- and [CoL2]- are best described with fragments in the electronic triplet state between the metal atoms with d8 configuration and triplet ligands. The singlet fragments give the degenerate TM(dπ)→L2 π backdonation as the strongest component, whereas the triplet fragments have the related electron-sharing TMq (dπ)-(L2)2- π bonding as the major component, differing only by the assignment of the bonded two electrons to one or both fragments. The calculations of the charge distribution using the Hirshfeld and Voronoi partitioning methods suggest that the metal atoms are nearly neutral or carry small negative charges in all complexes. The NBO method gives erratic charges, because of the neglect of the 4p AOs of the transition metals as genuine valence orbitals.

Revisiting the origin of the bending in group 2 metallocenes AeCp2 (Ae = Be-Ba)

Metallocenes are well-established compounds in organometallic chemistry, and can exhibit either a coplanar structure or a bent structure according to the nature of the metal center (E) and the cyclopentadienyl ligands (Cp). Herein, we re-examine the chemical bonding to underline the origins of the geometry and stability observed experimentally. To this end, we have analysed a series of group 2 metallocenes [Ae(C₅R₅)₂] (Ae = Be-Ba and R = H, Me, F, Cl, Br, and I) with a combination of computational methods, namely energy decomposition analysis (EDA), polarizability model (PM), and dispersion interaction densities (DIDs). Although the metal-ligand bonding nature is mainly an electrostatic interaction (65-78%), the covalent character is not negligible (33-22%). Notably, the heavier the metal center, the stronger the d-orbital interaction with a 50% contribution to the total covalent interaction. The dispersion interaction between the Cp ligands counts only for 1% of the interaction. Despite that orbital contributions become stronger for heavier metals, they never represent the energy main term. Instead, given the electrostatic nature of the metallocene bonds, we propose a model based on polarizability, which faithfully predicts the bending angle. Although dispersion interactions have a fair contribution to strengthen the bending angle, the polarizability plays a major role.

Stress-responsive properties of metallocenes in metallopolymers

This review article combines the field of metallopolymers and stress-responsiveness on a molecular level, namely, metallocenes, as emerging stress-responsive building blocks for materials.

Theoretical Insights on the High Pressure Behavior of Pentazolate Anion Complex [Co(H2O)4(N5)2]·4H2O

Periodic dispersion corrected density functional theory (DFT) calculations were carried out to examine the Hirshfeld surface, two dimensional (2D) fingerprint plots, crystal structure, molecular structure and density of state of all-nitrogen pentazolate anion complex [Co(H2O)4(N5)2]·4H2O under hydrostatic pressure from 0 to 20 GPa. The GGA/PW91-OBS method was applied in the present study. The intercontacts in [Co(H2O)4(N5)2]·4H2O were analyzed by Hirshfeld surfaces and 2D fingerprint plots. With ascending pressure, the lattice constants, compression rates, bond lengths, bond angles, and density of states change irregularly. Under 11.5, 13.0 and 15.8 GPa, covalent interaction competition is obvious between Co-N and Co-O bonds. It is possible to achieve orderly modification and regulation of the internal structure of [Co(H2O)4(N5)2]·4H2O by applied pressure. This is in accordance with the results from density of states analysis. The external compression causes the nonuniformity of electron density and the differential covalent interaction between pentazolate anion, coordinated water and atom Co. It is of great significance to interpret inter/intramolecular interaction and structural stability of [Co(H2O)4(N5)2]·4H2O and provide theoretical guidance for the design of metal complexes of all-nitrogen pentazolate anion.

Generalizing metallocene mechanochemistry to ruthenocene mechanophores

Recent reports have shown that ferrocene displays an unexpected combination of force-free stability and mechanochemical activity, as it acts as the preferred site of chain scission along the backbone of highly extended polymer chains. This observation raises the tantalizing question as to whether similar mechanochemical activity might be present in other metallocenes, and, if so, what features of metallocenes dictate their relative ability to act as mechanophores. In this work, we elucidate polymerization methodologies towards main-chain ruthenocene-based polymers and explore the mechanochemistry of ruthenocene. We find that ruthenocene, in analogy to ferrocene, acts as a highly selective site of main chain scission despite the fact that it is even more inert. A comparison of ruthenocene and ferrocene reactivity provides insights as to the possible origins of metallocene mechanochemistry, including the relative importance of structural and thermodynamic parameters such as bond length and bond dissociation energy. These results suggest that metallocenes might be privileged mechanophores through which highly inert coordination complexes can be made dynamic in a stimuli-responsive fashion, offering potential opportunities in dynamic metallo-supramolecular materials and in mechanochemical routes to reactive intermediates that are otherwise difficult to obtain.

Structural evolution of LiNn+ (n = 2, 4, 6, 8, and 10) clusters: mass spectrometry and theoretical calculations

Mixed nitrogen-lithium cluster cations LiN_n⁺ were generated by laser vaporization and analyzed by time-of-flight mass spectrometry. It is found that LiN₈⁺ has the highest ion abundance among the LiN_n⁺ ions in the mass spectrum. Density functional calculations were conducted to search for the stable structures of the Li–N clusters. The theoretical results show that the most stable isomers of LiN_n⁺ clusters are in the form of Li⁺(N₂)_n/2, and the order of their calculated binding energies is consistent with that of Li–N₂ bond lengths. The most stable structures of LiN_n⁺ evolve from one-dimensional linear type (C_∞v, n = 2; D_∞h, n = 4), to two-dimensional branch type (D_3h, n = 6), then to three-dimensional tetrahedral (T_d, n = 8) and square pyramid (C_4v, n = 10) types. Further natural bond orbital analyses show that electrons are transferred from the lone pair on N_α of every N₂ unit to the empty orbitals of lithium atom in LiN_2–8⁺, while in LiN₁₀⁺, electrons are transferred from the bonding orbital of the Li–N_α bonds to the antibonding orbital of the other Li–N_α bonds. In both cases, the N₂ units become dipoles and strongly interact with Li⁺. The average second-order perturbation stabilization energy for LiN₈⁺ is the highest among the observed LiN_n⁺ clusters. For neutral LiN_2–8 clusters, the most stable isomers were also formed by a Li atom and n/2 number of N₂ units, while that of LiN₁₀ is in the form of Li⁺(N₂)₃(η¹-N₄).

LiN_n⁺ clusters were generated by laser ablation and the LiN₈⁺ with tetrahedral Li⁺(N₂)₄ structure has the highest ion abundance.

A carbon-free inorganic-metal complex consisting of an all-nitrogen pentazole anion, a Zn(ii) cation and H2O

A carbon-free inorganic-metal complex [Zn(H₂O)₄(N₅)₂]·4H₂O was synthesized by the ion metathesis of [Na(H₂O)(N₅)]·2H₂O solution with Zn(NO₃)₂·6H₂O. The complex was well characterized by IR and Raman spectroscopy, elemental analysis (EA), powder X-ray diffraction (PXRD), and differential scanning calorimetry (DSC). The structure of the complex was confirmed by single-crystal X-ray crystallography and a Zn(ii) ion is coordinated in a quadrilateral bipyramid environment in which the axial position is formed by two nitrogen atoms (N1) from two pentazolate rings (cyclo-N₅^-) and the equatorial plane is formed by four oxygen atoms (O1) from four coordinated water molecules. The thermal analysis of [Zn(H₂O)₄(N₅)₂]·4H₂O reveals that although water plays an important role in stabilizing cyclo-N₅^-, dehydration does not cause immediate decomposition of the anion. However, cyclo-N₅^- decomposed into N₃^- and N₂ gas at 107.9 °C (onset). Based on its chemical compatibility and stability, the complex exhibits promising potential as a modern environmentally-friendly energetic material.

Maximum bonding fragment orbitals for deciphering complex chemical interactions

An optimal set of fragment orbitals is proposed as a simple and powerful tool for analyzing complex bonding interactions.

Theoretically predicted ferrocene analogues with triplet aromatic CB5H5 ligands

Three ferrocene analogues, D _5h (η⁵-CB₅H₅)₂M (M = Fe^2-, Co^-, and Ni), with triplet aromatic CB₅H₅ ligands have been predicted at TPSSh/6-311+G(d,p) level. We find that the M atom interacts drastically with the two CB₅H₅ ligands through a nearly fully-filled 3d subshell, which is different from (η⁵-C₅H₅)₂Fe. The natural population analyses suggest that (η⁵-CB₅H₅)₂M have an unconventional charge distribution, i.e., the M atom is negatively charged, while the two boron rings are positively charged. The analyses of the electronic and dynamic stabilities indicate that (η⁵-CB₅H₅)₂Co^- is the most stable among (η⁵-CB₅H₅)₂M. Thus, we theoretically confirm that the triplet aromatic CB₅H₅ cluster can be regarded as a potential new ligand. Our theoretical predictions are awaiting future experimental verification.

Computational insights into novel dicobalt polynitrogen: structure, stability, intermolecular interaction, and application

Previous studies have suggested that polynitrogen species are significant as potential candidates for superior energetic material. In this paper, the polynitrogen species of Co2(N5)4 were reasonably designed and studied by the density functional theory (DFT), and five isomers of Co2(N5)4 were selected. These species were explored in detail, including structure, stability, intermolecular interaction, and application. The five isomers, each with its own special structure feature, were stable enough based on the analysis of bond energy, chemical hardness, and aromaticity. Furthermore, the intermolecular interactions suggested the presence of a covalent interaction in the Co–Co and N–N bonds, the electronic delocalization in cyclo-N5, and the ionic feature in the Co–N bond. In addition, all of the title species held high-energy content. Compared with the known high energy density materials of HB(N5)3Be2(N5)3BH, energetic material of nitromethane, and famous nitramine explosive HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane), Co2(N5)4 holds a stronger advantage. The five Co2(N5)4 species were located at 27.8–35.8 kcal/mol per N2 unit, their energy densities were about 2.73 × 104 MJ/kg, and their mass densities were in the range of 2.60–2.74 g/cm3. Significantly, the 4-1 was the most stable, and its density was also the greatest among the five species. Thus, it has the most potential as a high energy density material.

Choosing the right precursor for thermal decomposition solution-phase synthesis of iron nanoparticles: tunable dissociation energies of ferrocene derivatives

Pathways to low-temperature thermal dissociation of ferrocene derivatives as iron nanoparticle precursors.

Theoretical investigations on stability of pyridylpentazoles, pyridazylpentazoles, triazinylpentazoles, tetrazinylpentazoles, and pentazinylpentazole searching for a replacement of phenylpentazole as N5 (-) source

Stabilities of pyridylpentazoles, pyridazylpentazoles, triazinylpentazoles, tetrazinylpentazoles, and pentazinylpentazole were studied using density functional theory to assess their potentials as the source of pentazole anion (N5 (-)) for replacement of phenylpentazole (PhN 5 ). Replacing the aryl group of PhN 5 by six-member heterocycle weakens pentazole ring. Compared to PhN 5 , title molecules have longer N-N bonds and lower activation energy (E a,1) needed for the N5 ring breaking. E a,1 decreases with the increasing number of nitrogen atoms of heterocycle. The ortho nitrogen of heterocycle most obviously lowers the stability of pentazole. The central C-N bond dissociation energies (BDEs) of title molecules are lower than that of PhN 5 . For the molecule with 0~1 ortho-nitrogen, H rearrangement happens during the central C-N bond breaking. The energy (E a,2) required for H rearrangement is considerably smaller than the corresponding BDE. ΔE a,2 (E a,2(PhN5) - E a,2 = 7.5~35.7 kJ mol(-1)) is larger than ΔE a,1 (E a,2(PhN5) - E a,2 = 4.6~15.5 kJ mol(-1)), while ΔE a,2/E a,2(PhN5) (2~9.5 %) is smaller than ΔE a,1/E a,1(PhN5) ( 4.4~15.0 %). The larger ΔE a,1/E a,1(PhN5) suggests that title molecules can not be the better N5 (-) than PhN 5 .

Ferrocene Orientation Determined Intramolecular Interactions Using Energy Decomposition Analysis

Two very different quantum mechanically based energy decomposition analyses (EDA) schemes are employed to study the dominant energy differences between the eclipsed and staggered ferrocene conformers. One is the extended transition state (ETS) based on the Amsterdam Density Functional (ADF) package and the other is natural EDA (NEDA) based in the General Atomic and Molecular Electronic Structure System (GAMESS) package. It reveals that in addition to the model (theory and basis set), the fragmentation channels more significantly affect the interaction energy terms (ΔE) between the conformers. It is discovered that such an interaction energy can be absorbed into the pre-partitioned fragment channels so that to affect the interaction energies in a particular conformer of Fc. To avoid this, the present study employs a complete fragment channel—the fragments of ferrocene are individual neutral atoms. It therefore discovers that the major difference between the ferrocene conformers is due to the quantum mechanical Pauli repulsive energy and orbital attractive energy, leading to the eclipsed ferrocene the energy preferred structure. The NEDA scheme further indicates that the sum of attractive (negative) polarization (POL) and charge transfer (CL) energies prefers the eclipsed ferrocene. The repulsive (positive) deformation (DEF) energy, which is dominated by the cyclopentadienyle (Cp) rings, prefers the staggered ferrocene. Again, the cancellation results in a small energy residue in favour of the eclipsed ferrocene, in agreement with the ETS scheme. Further Natural Bond Orbital (NBO) analysis indicates that all NBO energies, total Lewis (no Fe) and lone pair (LP) deletion all prefer the eclipsed Fc conformer. The most significant energy preferring the eclipsed ferrocene without cancellation is the interactions between the donor lone pairs (LP) of the Fe atom and the acceptor antibond (BD*) NBOs of all C–C and C–H bonds in the ligand, LP(Fe)-BD*(C–C & C–H), which strongly stabilizes the eclipsed (D_5h) conformation by −457.6 kcal·mol⁻¹.

The d-electrons of Fe in ferrocene: the excess orbital energy spectrum (EOES)

The EOES (Δε_i=εE-Fci −εS-Fci) shows that the orbitals with significantly excess energies are Fe d-electron dominant.

Structure, stability and intramolecular interaction of M(N5)2(M = Mg, Ca, Sr and Ba) : a theoretical study

The potential energetic materials, alkaline earth metal complexes of the pentazole anion (M(N₅)₂, M = Mg²⁺, Ca²⁺, Sr²⁺and Ba²⁺), were studied using the density functional theory.

Pyridylpentazole and its derivatives: a new source of N5−?

Pyridylpentazole (PyN₅) and its derivatives with 1–2 electron withdrawing groups were studied using the density functional theory to assess their potentials as the source of pentazole anion N₅⁻ for replacement of phenylpentazole (PhN₅).

Ferrocene analogues of sandwich M(CrB₆H₆)₂: a theoretical investigation

The stability and electronic structures of the new sandwich compounds M(CrB6H6)2 (M = Cr, Mn(+), Fe(2+)) are investigated by density functional theory. All the investigated sandwich complexes are in D(6d) symmetry and all of them are thermodynamically stable according to the large HOMO-LUMO gap, binding energy, vertical ionization potential and vertical electron affinity analyses, as well as Fe(C5H5)2 and Cr(C6H6)2, following the 18-electron principle. The natural bond orbital, detailed molecular orbitals and adaptive natural density partitioning analyses suggest that the spd-π interaction plays an important role in the sandwich compounds. This work challenges the traditional chemical bonding of the inorganic metal compound, investigates first the bonding style of the d atom orbital interacting with the π MO which was formed by p-d atomic orbitals, and indicates that the metal-doped borane ring can also be an ideal type π-electron donor ligand to stabilize the transition metal.

Chain or ring: which one is favorable in nitrogen-rich molecules N6XHm, N8XHm, and N10XHm (X = B, Al, Ga, m = 1 and X = C, Si, Ge, m = 2)?

A series of nitrogen-rich molecules N6XHm, N8XHm, and N10XHm (X = B, Al, Ga, m = 1 and X = C, Si, Ge, m = 2) consisting of N3 and N5 radicals, are systematically investigated by using B3LYP and B3PW91 DFT methods. It is found that for the nitrogen-rich molecules, the structures with N3-chains (N5-ring) are more stable than those containing a N3-ring (N5-chain). This result could be well-explained by the intrinsic stability of the N3 and N5 radicals and their charge distribution in nitrogen-rich molecules. The dissociation energies further indicate that the B-doped and C-doped structures are the most stable among the molecules with three elements of group 13 and 14, respectively. Energy decomposition analysis shows the bond of boron-nitrogen is stronger than that of carbon-nitrogen. Detailed bonding analysis demonstrates that the B-N bond is determined by σ and π interactions between the B and N atoms, whereas C-N bonds by only σ interactions. These results imply that the boron atom is more suitable than the carbon atom for building the nitrogen-rich molecules studied in this article.

CHARACTERIZATIONS OF NOVEL BINUCLEAR ALKALINE-EARTH METAL COMPOUNDS: M2(ηn-N5)2 (M=Be AND Mg, n = 1, 2; M=Ca, n = 2, 5)

The first structural characterization for the twelve binuclear alkaline-earth metal compounds M 2(ηn- N 5)2 ( M=Be , Mg ; n = 1, 2) and Ca 2(ηn- N 5)2(n = 2, 5) have been optimized with local energy minimum by density functional theory (DFT). The most energetically favored structures in M 2(ηn- N 5)2( M=Be , Mg , Ca ) are of D2d symmetry Be 2(η1- N 5)2, Mg 2(η2- N 5)2 and Ca 2(η2- N 5)2 and the metal–metal distances are 2.03 Å for Be-Be , 2.77 Å for Mg-Mg and 3.72 Å for Ca-Ca , which are significantly shorter than the experiment values of weakly bound bare diatomic Be 2, Mg 2 and Ca 2.1,2 Ca 2(η5- N 5)2 (D5d or D5h) is the only stable specie with sandwiched structure, bearing an even shorter Ca-Ca distance of 3.66 Å, and lying 24 kcal/mol higher in energy than the D2d structure. The dissociation enthalpies of the twelve M 2(ηn- N 5)2( M=Be , Mg , Ca ) to two M (ηn- N 5) fragments are predicted to be 72.6–73.1, 41.2–43.8, and 27.4–29.7 kcal/mol, respectively, implying a substantial metal–metal bonding. Natural bond orbital (NBO) analysis suggests that metal–metal bonds are of σ-bond. The natural charge of the alkali earth metal atom in the twelve M 2(ηn- N 5)2 species is larger than +0.88, which is consistent with the +1 oxidation state of the metal atoms. Nucleus-independent chemical shift (NICS) values confirm that the planar [Formula: see text] exhibits characteristics of aromaticity for these M 2(ηn- N 5)2 species.