Characterizations of B and N heteroatoms as substitutional doping on structure,stability, and aromaticity of novel heterofullerenes evolved from the smallest fullerene cage C20: a density functional theory perspective

The structural stabilities and electronic properties of C₂₀ fullerene and some its incorporated boron and nitrogen derivatives are probed at B3LYP/AUG‐cc‐pVTZ level of theory. According to density functional theory results, the topology of inserted B or N heteroatoms in [20]‐fullerene perturbs strongly the stability, energy, geometry, charge, polarity, nucleus‐independent chemical shifts, aromaticity, and highest‐occupied molecular orbital and lowest‐unoccupied molecular orbital (HOMO–LUMO) gap of the resulting heterofullerenes. Vibrational frequency (υ_min) calculations show that except N₁₀C₁₀, all other B_bN_nC_{20‐(b + n)} heterofullerenes with b, and n = 0, 4, 5, 8, and 10 are true minima. The calculated band gaps (?E_HOMO–LUMO) of B₈C₁₂, and N₈C₁₂ (2.86 eV), show them the most stable heterofullerenes against electronic excitations. While 10 B substituting in equatorial position increase the conductivity of B₁₀C₁₀ through decreasing its band gaps, 10 N doping in equatorial position enhance stability of N₁₀C₁₀ against electronic excitations via increasing its band gaps. High natural bond orbital and Mulliken charge transfer on the surfaces of B atoms, especially B₅N₅C₁₀with five B–N bonds in the equatorial position, provokes further investigation on its possible application for hydrogen storage. Nucleus‐independent chemical shift values show that B₅N₅C₁₀ is the most aromatic species. The calculated heat of atomization per carbon (ΔH_at/C) of B₈C₁₂ shows it the most thermodynamic stable heterofullerenes of each. Copyright © 2016 John Wiley & Sons, Ltd.     

Theoretical study on C100 fullerenes and C96X4 (X=N,P, B,Si)

The geometrical structures and electronic properties of six fullerene isomers of C₁₀₀ were studied at the HF/6-31G^? and B3LYP/6-31G^? levels, respectively. The results of the fully optimized calculations show that three C₁₀₀ isomers 449:D₂, 425:C₁ and 442:C₂ are near isoenergetic isomers. The energies and properties of C₁₀₀ hexaanions were calculated. The C₁₀₀^6? (450:D₅) isomer is predicted to be the most stable isomer at the B3LYP/6-31G^? level, and the C₁₀₀^6? (449:D₂) isomer is 44.1 kcal/mol higher in energy. The heterofullerenes C₉₆X₄ (X=N, P, B, Si) formed from the initial C₁₀₀ (449:D₂) have also been investigated at the B3LYP/6-31G^? level. The HOMO–LUMO gaps and aromaticities show that the replacement of fullerene carbon atoms with four heteroatoms can enhance the electronic stabilization of C₁₀₀ (449:D₂).     

A comparative study on the Ge6C14 heterofullerene nanocages: a density functional survey

Ten C₁₄Ge₆ heterofullerene isomers of C₂₀ have been investigated by density functional theory (DFT) methods with Becke 3‐Parameter (Exchange), Lee, Yang and Parr (B3LYP) functional at the 6‐311 + G*, 6‐311++G** and AUG‐cc‐pVTZ levels. In contrast to identical bonds in the latter, contractions of C═C double bonds are encountered at the expense of longer C―Ge bonds in the former. Vibrational frequency analysis confirms that all of the nanocages are true minima. In contrast to the common belief, for obtaining highly silicon‐doped stable heterofullerenes, that the silicon dopants must be completely isolated from each other by means of strong C═C double bonds. Here, linking the germanium substitutions together is an applicable strategy for obtaining highly doped stable isolated heterofullerenes since it avoids weak heteroatom─heteroatom bonds. Therefore, none of the computed heterofullerenes collapses to open, to deform, or to segregate fullerenic cages. As to band gaps (ΔE_HOMO‐LUMO), and nucleus‐independent chemical shifts at cage centers (NICS (0)), C₁₄Ge₆‐2 immerges with the highest value. Hence, it is predicted to be the most stable against electronic excitation. It contains 2 Ge─Ge single bonds at the cap‐equatorial positions. On the other hand, as to zero‐point vibrational energy and heat of atomization (ΔH_at), C₁₄Ge₆‐8 appears with the lowest and highest value, respectively. It contains 6 alternating germanium atoms in the equatorial and cap positions. Thus, it is predicted to be the most thermodynamically stable. So, germanium substitution leads to a high charge distribution on the surfaces of all the isomers specially C₁₄Ge₆‐9 with +1.496 charged germanum atoms. C₁₄Ge₆ isomers seem to be a good candidate for the hydrogen storage material.     

Adsorption of NH3 and NO2 molecules on C48B6N6 heterofullerene: A DFT study on electronic properties

Adsorption of NH₃ and NO₂ molecules on the external surface of C₄₈B₆N₆ heterofullerene is investigated using DFT method. Attachment of NH₃ and NO₂ on C₄₈B₆N₆ heterofullerenes are compared with the bare C₄₈B₆N₆ model optimized at the B3LYP/6-31G^? level. The high surface binding energies indicates that ammonia undergoes chemical adsorption and could be compatible with the long recovery time but C₄₈B₆N₆ should be good NO₂ sensors with quick response as well as short recovery time. Total (TDOS) and partial (PDOS) density of state calculations is also considered to elucidate the difference in the NH₃ and NO₂ gas detection mechanism of C₄₈B₆N₆. The overlap population density of state (OPDOS) indicated that the chemical adsorption is due to the overlap of atomic orbitals below the Fermi level. The calculated results suggest that the C₄₈B₆N₆ heterofullerene is a suitable sensor material for NO₂ and is an ideal material for elimination and filtering of ammonia.     

Raman and infrared,microwave spectra,conformational stability,adjusted r0 structural parameters,and vibrational assignments of cyclopentylamine

Fourier transform microwave spectrum of cyclopentylamine, c–C₅H₉NH₂ has been recorded, and seven transitions have been assigned for the most abundant conformer, and the rotational constants have been determined: A = 4909.46(5), B = 3599.01(4), and C = 2932.94(4). From the determined microwave rotational constants and ab initio MP2(full)/6‐311 + G(d,p) predicted structural values, adjusted r₀ parameters are reported with distances (Å): rC_α–C_β = 1.529(3), rC_β–C_γ = 1.544(3), rC_γ–C_γ = 1.550(3), rC_α–N = 1.470(3), and angles (°) ∠CCN = 108.7(5), ∠C_βC_αC_β = 101.4(5), and τC_βC_αC_β′C_γ′ = 42.0(5). The infrared spectra (4000–220 cm⁻¹) of the gas have been recorded. Additionally, the variable temperature (−60 to −100 °C) Raman spectra of the sample dissolved in liquefied xenon was recorded from (3800–50 cm⁻¹). The four possible conformers have been identified, and their relative stabilities obtained with enthalpy difference relative to t‐Ax of 211 ± 21 cm⁻¹ for t‐Eq ≥ 227 ± 22 cm⁻¹ for g‐Eq ≥ 255 ± 25 cm⁻¹ for g‐Ax. The percentage of the four conformers is estimated to be 53% for the t‐Ax, 11 ± 1% for t‐Eq, 20 ± 2% for g‐Ax and 16 ± 2% for g‐Eq at ambient temperature. The conformational stabilities have been predicted from ab initio calculations by utilizing several different basis sets up to aug‐cc‐pVTZ from both MP2(full) and density functional theory calculations by the B3LYP method. Vibrational assignments have been provided for the observed bands for all four conformers, which are predicted by MP2(full)/6‐31G(d) ab initio calculations to predict harmonic force constants, wavenumbers, infrared intensities, Raman activities, and depolarization ratios for all of the conformers. The results are discussed and compared to the corresponding properties of some related molecules. Copyright © 2014 John Wiley & Sons, Ltd.     

(Ca3N2)n(n=1—4)团簇结构与性质的密度泛函理论研究

用密度泛函理论(DFT)的杂化密度泛函B3LYP方法在6-31G^*基组水平上对(Ca₃N₂)_n(n=1—4)团簇各种可能的构型进行几何结构优化,预测了各团簇的最稳定结构.并对最稳定结构的振动特性、成键特性、电荷特性和稳定性等进行了理论分析.结果表明,(Ca₃N₂)_n(n=1—4)团簇最稳定构型中N原子为3—5配位,Ca—N键长为0.231—0.251nm,Ca—Ca键长为0.295—0.358nm;N原子的自然电荷在-1.553e—-2.241e之间,Ca原子的自然电荷在1.035e—1.445e之间,Ca和N原子间相互作用呈现较强的离子性,Ca₃N₂和(Ca₃N₂)₃团簇有相对较高的动力学稳定性. 关键词： 3N₂)_n(n=1—4)团簇')" href="#">(Ca₃N₂)_n(n=1—4)团簇 密度泛函理论 结构与性质     

Electronic structure and enhanced visible‐light absorption of N,B‐codoped TiO2

We present GGA+U calculations to investigate the electronic structure and visible‐light absorption of N,B‐codoped anatase TiO₂. The N_sB_i (substitutional N, interstitial B) codoped TiO₂ produces significant Ti 3d and N 2p mid‐gap states when the distance of N and B atoms is far, whereas the N_iB_i (interstitial N and B) and N_sB_s (substitutional N and B) codoped TiO₂ prefer to form localized p states at 0.3–1.2 eV above the valence band maximum. Further, the optical band edges of the three codoped systems shift slightly to the visible region, but only the far‐distance N_sB_i codoped TiO₂ clearly shows an optical transition. These results indicate that N_sB_i codoped TiO₂ has a dominant contribution to the optical absorption of N,B‐codoped TiO₂. (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

DFT investigation on the structures and stabilities of exohedral derivatives for D<Subscript>3</Subscript> C<Subscript>32</Subscript> fullerene: C<Subscript>32</Subscript>X<Subscript>n</Subscript> (X = H and Cl)

A systematic investigation of D₃ C₃₂ fullerene and its derivatives C₃₂X_n (X = H and Cl) has been performed using B3LYP/6-31G(d) method based on the density functional theory. The geometry structures, reaction energies, relative stabilities, and electronic properties have been studied. By investigating the possible C₃₂X_n (X = H and Cl) molecules, C₃₂H₂ and C₃₂Cl₂ behave more thermodynamically accessible with respect to other derivatives. The frontier molecular orbitals and electronic density of states calculations of C₃₂X₂ system indicate that H and Cl passivation have less contribution to the electronic structures, but significantly improve the stability of D₃ C₃₂ fullerene. Finally, the ¹³C NMR chemical shifts of C₃₂H₂ and C₃₂Cl₂ have been simulated to provide helpful information for further experiment identification.     

Short‐time photodecay dynamics of meso‐tetra(p‐hydroxyphenyl)porphine in the condensed phase explored via resonance Raman spectroscopy and density functional theory calculation

Resonance Raman (RR) spectra of free‐base meso‐tetra(p‐hydroxyphenyl)porphine(THPP) were obtained with 397.9, 416 and 514 nm excitation wavelengths, and density functional calculations were carried out to help the elucidation of the photorelaxation dynamics of Soret (B_x and B_y bands) electronic transitions and the RR spectra of THPP. The RR spectrum indicates that the Franck–Condon (FC) region photorelaxation dynamics for the S₀ → S₅ excited electronic state is predominantly along the totally symmetric C_m phC stretching and the C_βC_β stretchingand simultaneously along the asymmetric (C_mC_α)as stretching, ν(phC  C)asstretching, δ(NH)s and γ(C_βH) vibrational relaxation processes, while that for S₀ → S₄ electronic state is predominantly along the C_m phC stretching and pyrrole breathing. The excited‐state structural dynamics of THPP determined from the RR spectra shows that internal conversion (IC) B_y → B_x electronic relaxation occurs in tens of femtoseconds, and the short‐time dynamics is interpreted using the time‐dependent wave packet theory and Herzberg–Teller (vibronic coupling) contributions. Copyright © 2010 John Wiley & Sons, Ltd.     

Vibrational spectroscopy and DFT calculations of N,N′‐dicyclohexylcarbodiimide

Infrared (IR) and Raman spectra were obtained for N,N′‐dicyclohexylcarbodiimide (DCC) in the solid state and in CHCl₃ solution. Structures and vibrational spectra of isolated, gas‐phase DCC molecules with C₂ and C_i symmetries, computed at the B3‐LYP/cc‐pVTZ level, show that the IR and Raman spectra provide convincing evidence for a C₂ structure in both the solid state and in CHCl₃ solution. Using a scaled quantum‐chemical force field, these density functional theory calculations have provided detailed assignments of the observed IR and Raman bands in terms of potential energy distributions. Comparison of solid‐state and solution spectra, together with a Raman study of the melting behaviour of DCC, revealed that no solid‐state effects were evident in the spectra. Copyright © 2010 John Wiley & Sons, Ltd.     

Physical image versus structure relation. Part 14—an attempt to rationalize some acidic region 13C NMR‐pH titration shifts for tetraaza macrocycles throughout the conformational GIAO DFT computational results: a pendant‐arm cyclam case

The most probable time‐averaged conformations of three polyammonium cations H_n 2 ⁿ⁺ (n = 3–5) formed from the macrocyclic pentamine ligand ( 2 , scorpiand) [derivative of 1,4,8,11‐tetraazacyclotetradecane (cyclam)] were analyzed in order to elucidate an origin of ‘wrong‐way’ amine‐protonation shifts found in some ¹³C NMR pH‐profiles determined for the acidic H₂O/D₂O solution. These NMR trends were reproduced quite well in δ_Cs computed for multicomponent shapes of related cations, which were in turn elucidated by the best fitting experimental data to those predicted by the gauge‐independent atomic orbital (GIAO) B3LYP/6‐31G* method, including the IEF‐PCM approach. A consistent DFT methodology of the treatment of such equilibrated cationic mixtures is proposed. Moreover, a few novel ONIOM2‐GIAO B3LYP/6‐31G*:STO‐3G type supermolecular calculations were performed for a simulated presence of bulk water molecules surrounding H₅ 2 ⁵⁺. Copyright © 2009 John Wiley & Sons, Ltd.     

A combined experimental and theoretical study of the alkylation of 3,5‐dithioxo‐[1,2,4]triazepines

The chemo‐ and regioselective alkylation reactions of 3,5‐dithioxo[1,2,4]triazepine 1 in a basic medium with α,ω‐dibromoalkanes 2a – c , Br(CH₂)_nBr (n = 1–3), are studied experimentally and theoretically. These alkylations, which occur at the thioxo sulfur atom in position 5, afford mainly 5‐bromomethylthio‐2,7‐dimethyl‐ 2,3‐dihydro‐ 4H[1,2,4]triazepin‐3‐one 3 for n = 1, 6,8‐dimethyl‐5‐thioxo‐2,3,4,5‐tetrahydro‐6H[1,3]thiazolo[4,5‐d][1,2,4]triazepine 4 for n = 2 and 7,9‐dimethyl‐6‐thioxo‐2,3,4,5,6,7‐hexahydro[1,3]thiazino [4,5‐d][1,2,4]triazepine 5 for n = 3. Theoretical calculations have been carried out at the B3LYP/6‐31G* and B3LYP(benzene)/6‐311+G*//B3LYP/6‐31G* levels, in order to rationalize the experimental observations. Both chemo‐ and regio‐selectivities of the alkylation reactions are analyzed. Copyright © 2008 John Wiley & Sons, Ltd.     

Structures,stabilities and electronic properties of boron-doped silicon clusters B3Sin (n=1–17) and their anions

The structures, stabilities and electronic properties of neutral and anionic B₃Si_n (n?=?1–17) clusters have been systemically investigated on the basis of density functional theory at the B3LYP/6-311?+?G(d) level and CALYPSO structure prediction method. The structural searches show that three boron atoms tend to form B₃ triangle encapsulated into Si_n cages with the increasing number of silicon atoms. Most of the lowest energy structures can be derived by using the squashed pentagonal bipyramid structure of B₃Si₄ and B₃Si₄^? as the major building unit. The relative stabilities are studied based on the calculated binding energies, second-order difference of energies and HOMO–LUMO gaps of the lowest energy structures. In addition, Hirshfeld, natural population analysis, Bader approaches and natural electronic configuration are performed to explore the charge transfer. At last, molecular orbital, magnetic properties, IR, Raman and UV–vis spectra are also, respectively, analysed for providing strong support for essential theoretical and experimental research.     

Piperidine‐appended imidazolium ionic liquids as task‐specific catalysts: computational study,synthesis, and multinuclear NMR

Imidazolium ionic liquids (IMILs) with a piperidine moiety appended via variable length methylene spacers (with n = 1–4) were studied computationally to assess their potential to act as internal base for N‐heterocyclic carbene (NHC) generation. Proton transfer energies computed by B3LYP/6‐311+G(2d,p) were least endothermic for the basic‐IL with n = 3, whose optimized structure showed the shortest C₂‐H‐‐‐‐N(piperidine) distance. Inclusion of counter anion (Cl or NTf₂) caused dramatic conformational changes to enable close contact between the acidic C₂‐H and the anions. To examine the prospect for internal C₂‐H‐‐‐‐N coordination, multinuclear NMR data (¹H, ¹⁵N, and ¹³C) were computed by gauge independent atomic orbitals–density functional theory (GIAO‐DFT) in the gas phase and in several solvents by the PCM method for comparison with the experimental NMR data for the basic ILs (with n = 2–4) synthesized in the laboratory. These studies indicate that interactions with solvent and counter ion are dominant forces that could disrupt internal C₂‐H‐‐‐‐N coordination/proton transfer, making carbene generation from these basic‐ILs unlikely without an added external base. Therefore, the piperidine‐appended IMILs appear suitable for application as dual solvent/base in organic/organometallic transformations that require the use of mild base, without the necessity to alkylate at C‐2 to prevent N‐heterocyclic carbene formation. Copyright © 2016 John Wiley & Sons, Ltd.     

Density functional investigation of metal encapsulated X@C12Si8 heterofullerene (X=Li, Na, K, Be, Mg, Ca, Al, Ga)

The stability and the possible application of our recently reported SiC heterofullerenes inspire the investigation of their further stabilization through ion encapsulation. The endohedral complexes X@C₁₂Si₈, where X=Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺, Ca²⁺, Al³⁺, and Ga³⁺, are probed at the MPWB1K/6-311G? and B3LYP/6-311G* levels of theory. The optimized geometries show the expanding or contracting capability of C₁₂Si₈ in order to accommodate metal ion guests. The inclusion energies indicate the stability of the complexes compared to the components. Meanwhile, the calculated binding energies show the stabilization of C₁₂Si₈ through the inclusion of Be²⁺, Mg²⁺, Al³⁺, and Ga³⁺. The host-guest interaction that is probed through NBO atomic charges supports the obtained results. This study refers to “metal ion encapsulation” as a strategy for stabilization of SiC heterofullerenes.     

Ammonia adsorption on the C30B15N15 heterofullerene: DFT study of nuclear magnetic shielding and electric field gradient tensors of N and B nuclei

Ammonia adsorption on the external surface of C₃₀B₁₅N₁₅ heterofullerene was studied using density functional calculations. Three models of the ammonia-attached C₃₀B₁₅N₁₅ together with the perfect model were optimized at the B3LYP/6-31G^? level. The optimization process reveals that dramatic influences occurred for the geometrical structure of C₃₀B₁₅N₁₅ after ammonia adsorption; the B atom relaxes outwardly and consequently the heterofullerene distorts from the spherical form in the adsorption sites. The chemical shielding (CS) tensors and nuclear quadrupole coupling constants of B and N nuclei were calculated at the B3LYP/6-311G^** level. Our calculations reveal that the B atom is chemically bonded to NH₃ molecule. The B atom in the NH₃-attached form has the largest chemical shielding isotropic (CSI) value among the other boron nuclei. The C_Q parameters of B nuclei at the interaction sites are significantly decreased after ammonia adsorption.     

Hydrogen incorporation into BN fullerene-like nanostructures: A first-principles study

We performed density functional theory calculations to investigate the possibility of formation of endohedrally H@(BN)_n–fullerene (n: 24, 36, 60) and H@C₆₀ complexes for potential applications in solid-state quantum-computers. Spin-polarized approach within the generalized gradient approximation with the Perdew–Burke–Ernzerhof functional was used for the total energies and structural relaxation calculations. The calculated binding energies show that H atom being incorporated into B₆₀N₆₀ nanocage can form most stable complexes while the B₂₄N₂₄ and C₆₀ nanocages might form unstable complex with positive binding energy. We have also examined the penetration of an H atom into the respective nanocages and the calculated barrier energies indicate that the H atom prefers to penetrate into the B₂₄N₂₄ and B₆₀N₆₀ nanocages with barrier energy of about 0.47 eV (10.84 kcal/mol). Furthermore the binding characteristic is rationalized by analyzing the electronic structures. Our findings reveal that the B₆₀N₆₀ nanocage has fascinating potential application in future solid-state quantum-computers.     

Density functional study of Agn-1Y (n=2-10) clusters

Employing the density functional theory, we investigate the lowest-energy geometric, the stable and the electronic properties of Ag_n-1Y (n=2-10) clusters in this paper. The structural optimization and the frequency analysis are performed at the B3LYP/LANL2DZ level. Meanwhile, the differences in geometry, stability and electronic properties between Ag_n and Ag_n-1Y (n=2-10) clusters are also studied. The results show that for the doping of the yttrium atoms, the structures and the average binding lengths of the Ag_n clusters are greatly changed. In addition, the thermodynamic stabilities of the Ag_n clusters are enhanced generally with the doping of the Y atoms. In addition, the chemical stabilities of the Ag_n-1Y clusters are still improved compared with that of the three-dimensional Ag_n clusters.     

SimCn （m+n≤7）团簇的密度泛函研究

使用密度泛函理论（DFT）的杂化密度泛函B3LYP方法在6-31G*基组水平上对Si_mC_n（m+n≤7）团簇各种可能的构型进行几何结构优化,预测了各团簇的最稳定结构. 并对最稳定结构的平均结合能（E_b）,二阶能量差分（Δ₂E）和能隙（E_g）等进行了理论研究. 结果表明,随着原子个数的增加,SiC二元团 关键词： mC_n （m+n≤7）团簇')" href="#">Si_mC_n （m+n≤7）团簇 密度泛函理论 结构与性质     

Structure and properties of α‐cyclo[N]thiophenes as potential electronic materials – a theoretical study

Density functional theory (DFT) calculations of ring‐shaped α‐cyclo[N]thiophenes with N = 2 to N = 18 have been performed for ideal structures of high symmetry (point groups C_nv and D_nh) and for optimum energy structures of lower symmetry (D_2d, C_s, C_2v, C_i or C₁). Whereas the first three members of the series behave exceptionally the higher members are typical cyclothiophenes consisting of weakly interacting thiophene rings. In contrast to neutral compounds, cations and dications of cyclothiophenes with N ≥5 exhibit pronounced electron delocalizations along the carbon backbone. However, if the functional B3LYP is replaced by BH cations of large ring‐size cations show polaron‐type charge defects. According to broken symmetry DFT calculations dications with N = 14 and N = 18 have biradical character. These structures correspond to two‐polaron‐type structures rather than to dipolarons. The calculated vertical ionization energies of cyclo[N]thiophene are comparable with those of oligo[N]thiophenes of the same number of thiophene rings but the calculated absolute energies are probably too low at large ring size. Cyclothiophenes absorb light of lower energies than the related oligothiophenes. Cyclothiophenes belong to the strongly absorbing organic chromophores. In case of high molecular symmetry some of the excited states of cyclothiophenes are degenerate. The degeneracy is lifted with lower symmetries but the general absorption feature remains. The theoretical results are discussed with respect to recent experimental findings. Copyright © 2007 John Wiley & Sons, Ltd.