A Systemic Insight into Exohedral Actinides and Endohedral Borospherenes: An&Bm and An@Bn (An=U,Np, Pu; m = 28, 32, 34, 36, 38, 40; n = 36, 38, 40)

A series of exohedral actinide borospherenes, An&B_m, and endohedral borospherenes, An@B_n (An=U, Np, Pu; m = 28, 32, 34, 36, 38, 40; n = 36, 38, 40), have been characterized by density functional theory calculations. The electronic structures, chemical bond topological properties and spectra have been systematically investigated. It was found that An@B_n is more stable than An&B_n in terms of structure and energy, and UB₃₆ in an aqueous solution is the most stable molecular in this research. The IR and UV-vis spectra of An&B_m and An@B_n are computationally predicted to facilitate further experimental investigations. Charge-transfer spectroscopy decomposes the total UV-Vis absorption curve into the contributions of different excitation features, allowing insight into what form of electronic excitation the UV–Vis absorption peak is from the perspective of charge transfer between the An atoms and borospherenes.     

Endohedral and Exohedral Metalloborospherenes: M@B40 (M=Ca,Sr) and M&B40 (M=Be,Mg)

The recent discovery of the all‐boron fullerenes or borospherenes, D_2d B₄₀^?/0, paves the way for borospherene chemistry. Here we report a density functional theory study on the viability of metalloborospherenes: endohedral M@B₄₀ (M=Ca, Sr) and exohedral M&B₄₀ (M=Be, Mg). Extensive global structural searches indicate that Ca@B₄₀ ( 1 , C_2v, ¹A₁) and Sr@B₄₀ ( 3 , D_2d, ¹A₁) possess almost perfect endohedral borospherene structures with a metal atom at the center, while Be&B₄₀ ( 5 , C_s, ¹A′) and Mg&B₄₀ ( 7 , C_s, ¹A′) favor exohedral borospherene geometries with a η⁷‐M atom face‐capping a heptagon on the waist. Metalloborospherenes provide indirect evidence for the robustness of the borospherene structural motif. The metalloborospherenes are characterized as charge‐transfer complexes (M²⁺B₄₀^2?), where an alkaline earth metal atom donates two electrons to the B₄₀ cage. The high stability of endohedral Ca@B₄₀ ( 1 ) and Sr@B₄₀ ( 3 ) is due to the match in size between the host cage and the dopant. Bonding analyses indicate that all 122 valence electrons in the systems are delocalized as σ or π bonds, being distributed evenly on the cage surface, akin to the D_2d B₄₀ borospherene.     

Linear uranium complexes X2UL5 with L=cyanide, isocyanate: DFT evidence for similarities between uranyl (X=O) and uranocene (X=Cp) derivatives

A DFT study of the isostructural compounds [UO2L5](n-) with n=3-5 and linear [Cp2UL5](m-) with m=1-3 has been carried out for two different anionic ligands. Structurally stable structures are obtained for all systems. The coordination competition between cyanide (CN(-)) and isocyanide (NC(-)) as well as between cyanate (OCN(-)) and isocyanate (NCO(-)) has been studied in the uranyl case. A clear preference for cyanide and isocyanate complexes is reported. The coordination of five ligands in the equatorial plane is rationalised by the analysis of the MO diagram of both systems. Moreover, the qualitative comparison of the two MO diagrams shows a high similarity in agreement with the isolobality concept. The existence of linear [Cp2UL5](-) organometallic U(VI) complexes is thus proposed, as well as the possibility of obtaining complexes of both types for U(VI) and U(V) with OCN(-) ligands. In addition, the U(IV) linear metallocene is calculated to be stable for the latter ligand.     

σ Aromaticity Dominates in the Unsaturated Three‐Membered Ring of Cyclopropametallapentalenes from Groups 7–9: A DFT Study

Aromaticity, an old but still fantastic topic, has long attracted considerable interest of chemists. Generally, π aromaticity is described by π‐electron delocalization in closed circuits of unsaturated compounds whereas σ‐electron delocalization in saturated rings leads to σ aromaticity. Interestingly, our recent study shows that σ aromaticity can be dominating in an unsaturated three‐membered ring (3MR) of cyclopropaosmapentalene. An interesting question is raised: Can the σ aromaticity, which is dominant in the unsaturated 3MR, be extended to other cyclopropametallapentalenes? If so, how could the metal centers, ligands, and substituents affect the σ aromaticity? Here, we report a thorough theoretical study on these issues. The nucleus‐independent chemical shift calculations and the anisotropy of the current‐induced density plots reveal the dominant σ aromaticity in these unsaturated 3MRs. In addition, our calculations show that substituents on the 3MRs have significant effects on the σ aromaticity, whereas the ligand effect is particularly small.     

Packed to the Rafters: Filling up C60 with Rare Gas Atoms

How many rare gas atoms can be placed into a fullerene cage until the pressure becomes large enough to break the C₆₀ framework? The answer given by density functional and ab initio computations is surprising and underlines the high stability of this unique carbon structure.     

A DFT Study of the Modulation of the Antiaromatic and Open-Shell Character of Dibenzo[a,f]pentalene by Employing Three Strategies: Additional Benzoannulation,BN/CC Isosterism,and Substitution

Dibenzo[a,f]pentalene ( [a , f ]DBP ) is a highly antiaromatic molecule having appreciable open-shell singlet character in its ground state. In this work, DFT calculations at the B3LYP/6-311+G(d,p) level of theory were performed to explore the efficiency of three strategies, that is, BN/CC isosterism, substitution, and (di)benzoannulation of [a , f ]DBP , in controlling its electronic state and (anti)aromaticity. To evaluate the type and extent of the latter, the harmonic oscillator model of aromaticity (HOMA) and aromatic fluctuation (FLU) indices were used, along with the nucleus-independent chemical shift NICS-XY-scan procedure. The results suggest that all three strategies could be employed to produce either the closed-shell system or open-shell species, which may be in the singlet or triplet ground state. Triplet states have been characterized as aromatic, which is in accordance with Baird's rule. All the singlet states were found to have weaker global paratropicity than [a , f ]DBP . Additional (di)benzo fusion adds local aromatic subunit(s) and mainly retains the topology of the paratropic ring currents of the basic molecule. The substitution of two carbon atoms by the isoelectronic BN pair, or the introduction of substituents, results either in the same type and very similar topology of ring currents as in the parent compound, or leads to (anti)aromatic and nonaromatic subunits. The triplet states of all the examined compounds are also discussed.     

Aggregation behavior and non-covalent functionalization of borofullerenes B28, B38, and B40: A density functional theory investigation

Fullerene-like boron clusters (borofullerenes) are rising stars in the field of cluster chemistry. In this work, density functional theory calculations revealed that the recently reported small borofullerenes B_n (n = 28, 38, and 40) are all highly reactive and tend to form dimers and even trimers spontaneously. In addition, the non-covalent modification of these borofullerenes by various cycloparaphenylene nanorings can form stable host-guest systems with substantial intermolecular charge transfer at both ground and excited states. Our results demonstrate that the borofullerenes are versatile platform for exohedral functionalization, and are very promising candidates for the design of novel nanomaterials with desirable properties.     

A DFT study on 1,4‐dihydro‐1,4‐azaborinine annulated linear polyacenes: Absorption spectra,singlet‐triplet energy gap,aromaticity, and HOMO–LUMO energy modulation

Linear polyacene (LPA) mimics containing multiple heterocycles have been computationally designed by annulating 1,4‐dihydro‐1,4‐azaborinine moieties to benzene (aB₁–aB₅), naphthalene (aN₁–aN₅), anthracene (aA₁–aA₅), and tetracene (aT₁–aT₅) cores. DFT studies conducted on them using M06L/6‐311++G(d,p) method reveal a perfect planar structure for all and suggest the utilization of nitrogen lone pairs for aromatic π‐electron delocalization. The computed values of aromaticity indices such as HOMA, NICS, and dehydrogenation energy (E _dh) of heterocycles support strong aromatic character for each six‐membered ring in the LPA mimics. On the basis of the minimum value of the molecular electrostatic potential (V _min) observed on each LPA unit in the LPA mimics, the extended delocalization of π‐electrons is verified. The energetic parameter E _dh showed strong linear correlation with HOMA, NICS and V _min parameters, which strongly supports the multidimensional character of aromaticity in LPA mimics. The electronic property modification is shown by the theoretical absorption spectra data and singlet‐triplet energy gap (ΔE _ST). The bandgap and ΔE _ST tunings are achieved for LPA mimics by selecting appropriate number of azaborinine type units and the size of LPA core used for annulation. © 2017 Wiley Periodicals, Inc.     

In search of fullerene-based superacids: synthesis, X-ray structure, and DFT study of C60(C2F5)(5)H

Simply super! The perfluoroalkylfullerene C(60)(C(2)F(5))(5) H is the first structurally characterized perfluoroalkylated fullerene-based acid and is also predicted to be the first gas-phase fullerene-based superacid.     

Understanding the Reactivity of Supported Late Transition Metals on a Bare Anatase (101) Surface: A Periodic Conceptual DFT Investigation

The rapidly growing interest for new heterogeneous catalytic systems providing high atomic efficiency along with high stability and reactivity triggered an impressive progress in the field of single-atom catalysis. Nevertheless, unravelling the factors governing the interaction strength between the support and the adsorbed metal atoms remains a major challenge. Based on periodic density functional theory (DFT) calculations, this paper provides insight into the adsorption of single late transition metals on a defect-free anatase surface. The obtained adsorption energies fluctuate, with the exception of Pd, between −3.11 and −3.80 eV and are indicative of a strong interaction. Depending on the considered transition metal, we could attribute the strength of this interaction with the support to i) an electron transfer towards anatase (Ru, Rh, Ni), ii) s-d orbital hybridisation effects (Pt), or iii) a synergistic effect between both factors (Fe, Co, Os, Ir). The driving forces behind the adsorption were also found to be strongly related to Klechkowsky's rule for orbital filling. In contrast, the deviating behaviour of Pd is most likely associated with the lower dissociation enthalpy of the Pd−O bond. Additionally, the reactivity of these systems was evaluated using the Fermi weighted density of states approach. The resulting softness values can be clearly related to the electron configuration of the catalytic systems as well as with the net charge on the transition metal. Finally, these indices were used to construct a model that predicts the adsorption strength of CO on these anatase-supported d-metal atoms. The values obtained from this regression model show, within a 95 % probability interval, a correlation of 84 % with the explicitly calculated CO adsorption energies.     

C68 Fullerene Isomers,Anions, and their Metallofullerenes: Charge‐Stabilizing Different Isomers

The complete set of 6332 classical isomers of the fullerene C₆₈ as well as several non‐classical isomers is investigated by PM3, and the data for some of the more stable isomers are refined by the DFT‐based methods HCTH and B3LYP. C₂:0112 possesses the lowest energy of all the neutral isomers and it prevails in a wide range of temperatures. Among the fullerene ions modeled, C₆₈^2?, C₆₈^4? and C₆₈^6?, the isomers C₆₈^2?(C_s:0064), C₆₈^4?(C_2v:0008), and C₆₈^6?(D₃:0009) respectively, are predicted to be the most stable. This reveals that the pentagon adjacency penalty rule (PAPR) does not necessarily apply to the charged fullerene cages. The vertical electron affinities of the neutral C_s:0064, C_2v:0008, and D₃:0009 isomers are 3.41, 3.29, and 3.10 eV, respectively, suggesting that they are good electron acceptors. The predicted complexation energy, that is, the adiabatic binding energy between the cage and encapsulated cluster, of Sc₂C₂@C₆₈(C_2v:0008) is ?6.95 eV, thus greatly releasing the strain of its parent fullerene (C_2v:0008). Essentially, C₆₈ fullerene isomers are charge‐stabilized. Thus, inducing charge facilitates the isolation of the different isomers. Further investigations show that the steric effect of the encaged cluster should also be an important factor to stabilize the C₆₈ fullerenes effectively.     

Pentapnictogen heterocyclic monoanions: Structure,stability, and aromaticity

Inorganic planar ring-shape molecules with 4n + 2 π electrons are always the focus of experimental synthesis and theoretical research due to their potential aromaticity and stability. In this work, the whole series of five-membered heterocycle monoanions X _nY_5-n⁻ (X, Y = group 15 elements; n = 1-4) were thoroughly investigated by means of density functional theory calculations. They all have large formation energies and HOMO-LUMO gap energies, suggesting the potential thermodynamic and kinetic stability. Their aromaticities are comparable to that of typical aromatic hydrocarbons. Their thermal stabilities were firmly established by the ab initio molecular dynamics simulations. As most of them are predicted for the first time, their various spectra were simulated for experimental characterization. Furthermore, we demonstrate that these five-membered cyclic anions can be employed as η⁵-ligand to construct novel all-inorganic metallocenes, which may serve as the building blocks of low-dimensional nanomaterials.     

Nature of Bonding in Bowl‐Like B36 Cluster Revisited: Concentric (6π+18π) Double Aromaticity and Reason for the Preference of a Hexagonal Hole in a Central Location

The bowl‐shaped C_6v B₃₆ cluster with a central hexagon hole is considered an ideal molecular model for low‐dimensional boron‐based nanosystems. Owing to the electron deficiency of boron, chemical bonding in the B₃₆ cluster is intriguing, complicated, and has remained elusive despite a couple of papers in the literature. Herein, a bonding analysis is given through canonical molecular orbitals (CMOs) and adaptive natural density partitioning (AdNDP), further aided by natural bond orbital (NBO) analysis and orbital composition calculations. The concerted computational data establish the idea of concentric double π aromaticity for the B₃₆ cluster, with inner 6π and outer 18π electron counting, which both conform to the (4n+2) Hückel rule. The updated bonding picture differs from existing knowledge of the system. A refined bonding model is also proposed for coronene, of which the B₃₆ cluster is an inorganic analogue. It is further shown that concentric double π aromaticity in the B₃₆ cluster is retained and spatially fixed, irrespective of the migration of the hexagonal hole; the latter process changes the system energetically. The hexagonal hole is a destabilizing factor for σ/π CMOs. The central hexagon hole affects substantially fewer CMOs, thus making the bowl‐shaped C_6v B₃₆ cluster the global minimum.     

Sandwich Compounds of Transition Metals with Cyclopolyenes and Isolobal Boron Analogues

A series of sandwich compounds of transition metals (M=Ni, Fe, Cr) with cyclic hydrocarbon (M(CH)_n) and borane (M(BH₂)_n), ligands (including mixed hydrocarbon/borane sandwiches) has been studied using density functional theory (B3LYP/6‐311+G(df,p)). Multicenter bonding between the central metal atom and basal cycloborane rings provides stabilization to planar cycloborane species. Large negative NICS values allude to aromatic character in the cycloboranes similar to the analogous cyclic hydrocarbons. The ability of cycloborane sandwiches to stabilize attached carbocations, radicals and carbanions is also assessed.     

Adsorption of Naphthalene and Quinoline on Pt,Pd and Rh: A DFT Study

The adsorption of naphthalene and quinoline on Pt(111), Pd(111) and Rh(111) surfaces is studied using density functional theory. The metal surfaces are simulated by means of large confined clusters and for Pt by means of a slab with periodic boundary conditions (PBC). Calculation parameters such as basis set convergence, basis set superposition error and effects of cluster relaxation and size are analyzed in order to assess the aptness of the cluster model. For all the metals, the preferred sites of adsorption are analyzed, thus revealing their different behaviors concerning structure and stability of adsorption modes. On Pt, the molecules have the richest theoretical configurational variety. Naphthalene and quinoline are found to adsorb preferentially on di‐bridge[7] sites on the three metals, and Rh exhibits higher adsorption energies than Pt and Pd. Structural features of the adsorbed molecules are correlated to the calculated adsorption energies. The di‐bridge[7] adsorption modes are studied in deeper detail decomposing the adsorption energies in two terms arising from molecular distortion and binding interaction to the metal. Molecular distortion is correlated to the HOMO–LUMO energy gap. The larger adsorption energies found for interactions with Rh result from the lower contribution of the distortion term. Binding interactions are described by analyzing the wave functions of naphthalene and quinoline adsorbed on a subunit of the large clusters in order to reduce the complexity of the analysis. Molecular orbitals are studied using concepts of Frontier Molecular Orbitals theory. This approach reveals that in the adsorption of naphthalene and quinoline on Pt and Pd, an antibonding state lies below the Fermi energy, while on Rh all antibonding states are empty, in agreement with the larger interaction energies. In addition, further insight is gained by projecting the density of states on the d band of the clean surfaces and of the adsorbed systems. This results in the rationalization of the structural features in terms of the concepts of electronic structure theory. The distributions of electronic density are described by means of Hirshfeld charges and isosurfaces of differential electron density. The net electron transfer from the metals to the molecules for most of the sites correlates with the trends of the adsorption energies.