Theoretical investigation of LaC3n+ (n = 0, 1, 2) clusters by density functional theory

LaC₃ⁿ⁺ (n=0, 1, 2) clusters have been studied using B3LYP (Becke 3-parameter–Lee-Yang-Parr) density functional method. The basis set is Dunning/Huzinaga valence double zeta for carbon and [2s2p2d] for lanthanum, denoted LANL1DZ. Four isomers are presented for each cluster; two of them are edge binding isomers with C_2v symmetry, the other two are linear chains with C_∞v symmetry. Meanwhile, two spin states for each isomer, that is, singlet and triplet for LaC₃⁺, doublet and quartet for LaC₃ and LaC₃²⁺, respectively, are also considered. Geometries, vibrational frequencies, infrared intensities, and other quantities are reported and discussed. The results indicate that at some spin states; the C_2v symmetry isomers are the dominant structures, while for the other spin states, linear isomers are energetically favored. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 301–307, 1998     

Structures and electronic properties of the small rubidium‐doped silicon RbSin (n = 1–12) clusters

The geometries, relative stabilities, and electronic properties of small rubidium‐doped silicon clusters RbSi_n (n = 1–12) have been systematically investigated using the density functional theory at the B3LYP/GENECP level. The optimized structures show that lowest‐energy isomers of RbSi_n are similar with the ground state isomers of pure Si_n clusters and prefer the three‐dimensional for n = 3–12. The relative stabilities of RbSi_n clusters have been analyzed on the averaged binding energy, fragmentation energy, second‐order energy difference, and highest occupied molecular orbital‐lowest unoccupied molecular orbital energy gap. The calculated results indicate that the doping of Rb atom enhances the chemical activity of Si_n frame and the magic number is RbSi₂. The Mulliken population analysis reveals that the charges in the corresponding RbSi_n clusters transfer from the Rb atom to Si atoms. The partial density of states and chemical hardness are also discussed. © 2014 Wiley Periodicals, Inc.     

A theoretical study on the stability difference of the borane and carborane C2Bn−2Hn (5 ≤ n ≤ 7) clusters

In order to study the electronic structure and structural stability of borane and carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters, especially the stability difference between the borane and carborane C₂B₃H₅. The frontier orbital energy levels of the borane and carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters are calculated at CCSD(T)/aug‐cc‐pVXZ//B3LYP/def2‐TZVPP level. The results are further analyzed by qualitative frontier orbital method based on the cap–ring interaction. The results reveal that: (1) the larger E_gap(HOMO‐LUMO energy gap) of carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters than borane (5 ≤ n ≤ 7) clusters originates from the more effective cap–ring orbital overlap of carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters than that of borane (5 ≤ n ≤ 7) clusters; (2) the smallest E_gap of the borane results from the highest energy level of the ring symmetry‐adapted linear combination orbital of cluster; and (3) the largest E_gap of the carborane C₂B₃H₅ is induced by the most effective cap–ring orbital interaction of C₂B₃H₅ cluster. © 2014 Wiley Periodicals, Inc.     

Global Minimum‐Energy Structure and Spectroscopic Properties of I2.−⋅n H2O Clusters: A Monte Carlo Simulated Annealing Study

The vibrational (IR and Raman) and photoelectron spectral properties of hydrated iodine‐dimer radical‐anion clusters, I₂^.? ? n H₂O (n=1–10), are presented. Several initial guess structures are considered for each size of cluster to locate the global minimum‐energy structure by applying a Monte Carlo simulated annealing procedure including spin–orbit interaction. In the Raman spectrum, hydration reduces the intensity of the I? I stretching band but enhances the intensity of the O? H stretching band of water. Raman spectra of more highly hydrated clusters appear to be simpler than the corresponding IR spectra. Vibrational bands due to simultaneous stretching vibrations of O? H bonds in a cyclic water network are observed for I₂^.? ? n H₂O clusters with n≥3. The vertical detachment energy (VDE) profile shows stepwise saturation that indicates closing of the geometrical shell in the hydrated clusters on addition of every four water molecules. The calculated VDE of finite‐size small hydrated clusters is extrapolated to evaluate the bulk VDE value of I₂^.? in aqueous solution as 7.6 eV at the CCSD(T) level of theory. Structure and spectroscopic properties of these hydrated clusters are compared with those of hydrated clusters of Cl₂^.? and Br₂^.?.     

Theoretical study of silicon–sulfur clusters (SiS2)n− (n=1–6)

The possible geometrical structures and relative stability of silicon–sulfur clusters (SiS₂) (n=1–6) are explored by means of density functional theory (DFT) quantum chemical calculations. We also compare DFT with second‐order Møller–Plesset (MP2) and Hartree–Fock (HF) methods. The effects of polarization functions, diffuse functions, and electron correlation are included in MP2 and B3LYP quantum chemical calculations, and B3LYP is effective in larger cluster structure optimization, so we can conclude that the DFT approach is useful in establishing trends. The electronic structures and vibrational spectra of the most stable geometrical structures of (SiS₂)_n⁻ are analyzed by B3LYP. As a result, the regularity of the (SiS₂)_n⁻ cluster growing is obtained, and the calculation may predict the formation mechanism of the (SiS₂)_n⁻ cluster. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 280–290, 2001     

Analysis of Why Boron Avoids sp2 Hybridization and Classical Structures in the BnHn+2 Series

We performed global minimum searches for the B_nH_n+2 (n=2‐5) series and found that classical structures composed of 2c–2e B? H and B? B bonds become progressively less stable along the series. Relative energies increase from 2.9 kcal mol^?1 in B₂H₄ to 62.3 kcal mol^?1 in B₅H₇. We believe this occurs because boron atoms in the studied molecules are trying to avoid sp² hybridization and trigonal structure at the boron atoms, as in that case one 2p‐AO is empty, which is highly unfavorable. This affinity of boron to have some electron density on all 2p‐AOs and avoiding having one 2p‐AO empty is a main reason why classical structures are not the most stable configurations and why multicenter bonding is so important for the studied boron–hydride clusters as well as for pure boron clusters and boron compounds in general.     

Stability,Infrared Spectra and Electronic Structures of (ZrO2)n(n=3-6)Clusters:DFT Study

The stability, infrared spectra and electronic structures of (ZrO₂)_n (n=3–6) clusters have been investigated by using density‐functional theory (DFT) at B3LYP/6‐31G* level. The lowest‐energy structures have been recognized by considering a number of structural isomers for each cluster size. It is found that the lowest‐energy (ZrO₂)₅ cluster is the most stable among the (ZrO₂)_n (n=3–6) clusters. The vibration spectra of Zr? O stretching motion from terminal oxygen atom locate between 900 and 1000 cm^?1, and the vibrational band of Zr? O? Zr? O four member ring is obtained at 600–700 cm^?1, which are in good agreement with the experimental results. Mulliken populations and NBO charges of (ZrO₂)_n clusters indicate that the charge transfers occur between 4d orbital of Zr atoms and 2p orbital of O atoms. HOMO‐LUMO gaps illustrate that chemical stabilities of the lowest‐energy (ZrO₂)_n (n=3–6) clusters display an even‐odd alternating pattern with increasing cluster size.     

Structure and Bonding in Stannadiphospholes and their Dianions SnC2P2R2m (R=H,tBu m=0, −2): A Comparative Study with C5H5+ and C5H5− Analogues

The potential energy surfaces of both neutral and dianionic SnC₂P₂R₂ (R=H, tBu) ring systems have been explored at the B3PW91/LANL2DZ (Sn) and 6‐311+G* (other atoms) level. In the neutral isomers the global minimum is a nido structure in which a 1,2‐diphosphocyclobutadiene ring (1,2‐DPCB) is capped by the Sn. Interestingly, the structure established by X‐ray diffraction analysis, for R=tBu, is a 1,3‐DPCB ring capped by Sn and it is 2.4 kcal mol^?1 higher in energy than the 1,2‐DPCB ring isomer. This is possibly related to the kinetic stability of the 1,3‐DPCB ring, which might originate from the synthetic precursor ZrCp₂tBu₂C₂P₂. In the case of the dianionic isomers we observe only a 6π‐electron aromatic structure as the global minimum, similarly to the cases of our previously reported results with other types of heterodiphospholes. 1 , 4 , 19 The existence of large numbers of cluster‐type isomers in neutral and 6π‐planar structures in the dianions SnC₂P₂R₂^2? (R=H, tBu) is due to 3D aromaticity in neutral clusters and to 2D π aromaticity of the dianionic rings. Relative energies of positional isomers mainly depend on: 1) the valency and coordination number of the Sn centre, 2) individual bond strengths, and 3) the steric effect of tBu groups. A comparison of neutral stannadiphospholes with other structurally related C₅H₅⁺ analogues indicates that Sn might be a better isolobal analogue to P⁺ than to BH or CH⁺. The variation in global minima in these C₅H₅⁺ analogues is due to characteristic features such as 1) the different valencies of C, B, P and Sn, 2) the electron deficiency of B, 3) weaker pπ–pπ bonding by P and Sn atoms, and 4) the tendency of electropositive elements to donate electrons to nido clusters. Unlike the C₅H₅⁺ systems, all C₅H₅^? analogues have 6π‐planar aromatic structures as global minima. The differences in the relative ordering of the positional isomers and ligating properties are significant and depend on 1) the nature of the π orbitals involved, and 2) effective overlap of orbitals.     

Structures and properties of (H2GaN3)n (n = 1–4) clusters: A DFT study

Computations on the systems of (H₂GaN₃)_n (n = 1–4) are performed using the density functional theory (DFT)/B3LYP method with different basis sets. (H₂GaN₃)₂ possessing D_2h symmetry is found to exhibit the planar Ga₂N₂ ring structure. (H₂GaN₃)₃ involving a six‐membered Ga₃N₃ ring is found to exhibit two minima with very similar binding energies (ca. −235 ∼ −231 kJ · mol⁻¹). One minimum is the newly found boat‐like conformation possessing C_s symmetry. Another minimum possessing C_3v symmetry is the chair‐like conformation. (H₂GaN₃)₄ occurs in several structures with Ga₄N₄ eight‐membered ring structures that correspond to minima with slight energy differences among them. The structural changes of the clusters are large compared with the monomer. Frequency calculations are carried out on each optimized structure, and their infrared (IR) spectra are discussed. Thermodynamic properties demonstrate that the systems of H₂GaN₃ occur at dimer–trimer–tetramer equilibrium, and the trimer is the main component. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Cationic Closo Carboranes—Promising Weakly Coordinating Ions

The unexplored carbon rich cationic closo carboranes, C₃B_n?3H_n⁺¹ (n=5, 6, 7, 10, 12) are investigated theoretically. The position isomers were calculated at the B3LYP/6‐31G* level, and the charge distribution in the cluster is estimated by NBO analysis. The criterion of ring‐cap orbital overlap compatibility along with the number of B? C, C? C, and B? B bonds help in explaining the stability order in each category. The most stable isomer is the one with maximum ring‐cap orbital overlap and largest number of B? C bonds. The order of relative stability among the trigonal bipyramid is 1c > 1b > 1a ′, where the stability is proportional to the number of CH caps over the small three‐membered ring. The C₃B₃H₆⁺ isomer with the one allyl C3 group ( 2b ) is more favorable than the one with a cyclopropenyl group ( 2a ). Among the C₃B₄H₇⁺ isomers the stability order is 3e > 3d > 3c > 3b > 3a , which mostly depends on the ring‐cap orbital overlap. In the bicapped square antiprism (4) where there is large number of isomers, the order follows the rule of ring cap compatibility and the number of B? C bonds. The order of 5e > 5d > 5c > 5b > 5a obtained from the calculations is in perfect agreement with the above sited rules. Equations (1) – (5) devised for estimating the stability of isomers of C₃B_n?3H_n⁺ indicate an increase in stability with cage size. The mono‐positive charge of the isomers is distributed throughout the cage, making them suitable candidates as weakly electrophillic cations. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1542–1551, 2001     

Geometry,stability, and isomerization of BnN2 (n = 1−6) isomers

We perform a systematic study on the geometry, stability, nature of bonding, and potential energy surface of low‐lying isomers of planar and cyclic B_nN₂ (n = 1?6) at the CCSD(T)/6‐311+G(d)//B3LYP/6‐311+G(d) level. B_nN₂ (n = 2?4) clusters are structurally similar to pure boron clusters. The evolution of the binding energy per atom, incremental binding energy, and second‐order difference of total energy with the size of B_nN₂ reveals that the lowest energy isomer of B₃N₂ has high stability. B₅N₂ and B₆N₂ possess π‐aromaticity according to Hückel (4n + 2) rule. The aromaticity of some isomers of B₄N₂ and B₆N₂ is examined based on their valence molecular orbitals. At the CCSD(T)/6‐311+G(d)//B3LYP/6‐311+G(d) level, several B₂N₂, B₃N₂, B₄N₂, and B₅N₂ isomers are predicted to be stable both thermodynamically and kinetically, and detectable in future experiments. © 2013 Wiley Periodicals, Inc.     

Structures and Reactivity of Oxygen‐Rich Scandium Cluster Anions ScO3–5−

Oxygen‐rich scandium cluster anions ScO_3–5^? are prepared by laser ablation and allowed to react with n‐butane in a fast‐flow reactor. A time‐of‐flight mass spectrometer is used to detect the cluster distribution before and after the reactions. The ScO₃^? and ScO₄^? clusters can react with n‐butane to produce ScO₃H^?, ScO₃H₂^?, and ScO₄H^?, while the more oxygen‐rich cluster ScO₅^? is inert. The experiment suggests that unreactive cluster isomers of ScO₃^? and ScO₄^? are also present in the cluster source. Density functional theory and ab initio methods are used to calculate the structures and reaction mechanisms of the clusters. The theoretical results indicate that the unreactive and reactive cluster isomers of ScO_3,4^? contain peroxides (O₂^2?) and oxygen‐centered radicals (O^.?), respectively. The mechanisms and energetics for conversion of unreactive O₂^2? to reactive O^.? species are also theoretically studied.     

The Haloboration of Tri‐tert‐butylazadiboriridine NB2R3: Ring Opening versus NB3 Cluster Formation

The azadiboriridine [–BR–NR–BR–] ( 1 ; R = tBu) is bromoborated at the B–B bond by alkyldibromoboranes R′BBr₂ to give the products Br–BR–NR=BR–BR′–Br ( 8 a – g : R′ = Me, Bu, iBu, Bzl, CH₂CHEt₂, CH₂Cy, CH₂(4‐C₆H₄tBu)). Two isomers of each of the products 8 a – g are formed and attributed to a cis/trans isomerism at the BN double bond; the isomerization is followed thermodynamically and kinetically by NMR methods with 8 a – d . The analogous chloroboration of 1 with BCl₃ yields Cl–BR–NR=BR–BCl₂ ( 8 h ), which at ambient temperature undergoes a degenerate exchange of the ligands Cl and BCl₂ along the B–N–B skeleton. At room temperature, the isomer Cl–BR–NR=BCl–BR–Cl ( 8 h ′) is slowly formed by an irreversible exchange of R and Cl along the B–B bond of 8 h . Different from BCl₃, the chloroborane BH₂Cl is simply added to the B–B bond of 1 under formation of the aza‐nido‐tetraborane NB₃R₃H₂Cl ( 2 b ). The chloroborane BHCl₂ gives a mixture of 8 h ′ and 2 b upon addition to 1 , apparently according to a preceding dismutation into BCl₃ and BH₂Cl. The configuration at the B3 atom of the nido‐clusters NB₃R₃H₂X (X = H, Cl) is discussed on the basis of the corresponding model molecules NB₃Me₃H₂X, whose structure and NMR signals are computed by the B3LYP method. The boranes 8 b – g can be debrominated with Li in the presence of tmen on applying ultrasound. The products are found to be the B‐borylated azadiboriridines [–BR–NR–B(BRR′)–] ( 9 b – g ). The 2‐borylazadiboriridines NB₃H₄ ( 9 h ) and NB₃Me₄ ( 9 i ) were found as local minima on the energy hyperface by the B3LYP method, but minima for structural isomers with lower energy were also found; the tetrahedral clusters NB₃R₄ give high‐energy minima with triplet ground states. Computations of the ¹¹B NMR shifts of 9 h and 9 i support the proposed structures of 9 b – g .     

Regiodirected Synthesis and Stereochemistry of 2,4,8‐Trialkyl‐3‐thia‐1,5‐diazabicyclo[3.2.1]octanes and α,ω‐Bis(2,4,6‐trialkyl‐1,3,5‐dithiazinane‐5‐yl)alkanes

2,4,8‐Trialkyl‐3‐thia‐1,5‐diazabicyclo[3.2.1]octanes have been obtained by the regioselective and stereoselective cyclocondensation of 1,2‐ethanediamine with aldehydes RCHO (R═Me, Et, Prⁿ, Buⁿ, Pentⁿ) and H₂S at molar ratio 1:3:2 at 0°C. The increase in molar ratio of thiomethylation mixture RCHO–H₂S (6:4) at 40°C resulted in selective formation of bis‐(2,4,6‐trialkyl‐1,3,5‐dithiazinane‐5‐yl)ethanes. Cyclothiomethylation of aliphatic α,ω‐diamines with aldehydes RCHO (R═Me, Et) and H₂S at molar ratio 1:6:4 and at 40°С led to α,ω‐bis(2,4,6‐trialkyl‐1,3,5‐dithiazinane‐5‐yl)alkanes. Stereochemistry of 2,4,8‐trialkyl‐3‐thia‐1,5‐diazabicyclo[3.2.1]octanes have been determined by means of ¹H and ¹³С NMR spectroscopy and further supported by DFT calculations at the B3LYP/6‐31G(d,p) level. The structure of α,ω‐bis(2,4,6‐trialkyl‐1,3,5‐dithiazinane‐5‐yl)alkanes was confirmed by single‐crystal X‐ray diffraction study.     

Theoretical prediction for the structures of gas phase lithium oxide clusters: (Li2O)n (n = 1–8)

We apply genetic algorithm combining directly with density functional method to search the potential energy surface of lithium‐oxide clusters (Li₂O)_n up to n = 8. In (Li₂O)_n (n = 1–8) clusters, the planar structures are found to be global minimum up to n = 2, and the global minimum structures are all three‐dimensional at n ≥ 3. At n ≥ 4, the tetrahedral unit (TU) is found in most of the stable structures. In the TU, the central Li is bonded with four O atoms in sp³ interactions, which leads to unusual charge transformation, and the probability of the central Li participating in the bonding is higher by adaptive natural density partitioning analysis, so the central Li is in particularly low positive charge. At large cluster size, distortion of structures is viewed, which breaks the symmetry and may make energy higher. The global minimum structures of (Li₂O)₂, (Li₂O)₆, and (Li₂O)₇ clusters are the most stable magic numbers, where the first one is planar and the later both have stable structural units of tetrahedral and C_4v. © 2012 Wiley Periodicals, Inc.     

Opening of the Diamondoid Cage upon Ionization Probed by Infrared Spectra of the Amantadine Cation Solvated by Ar,N2, and H2O

Radical cations of diamondoids, a fundamental class of very stable cyclic hydrocarbon molecules, play an important role in their functionalization reactions and the chemistry of the interstellar medium. Herein, we characterize the structure, energy, and intermolecular interaction of clusters of the amantadine radical cation (Ama⁺, 1-aminoadamantane) with solvent molecules of different interaction strength by infrared photodissociation (IRPD) spectroscopy of mass-selected Ama⁺L_n clusters, with L=Ar (n≤3) and L=N₂ and H₂O (n=1), and dispersion-corrected density functional theory calculations (B3LYP−D3/cc-pVTZ). Three isomers of Ama⁺ generated by electron ionization are identified by the vibrational properties of their rather different NH₂ groups. The ligands bind preferentially to the acidic NH₂ protons, and the strength of the NH…L ionic H-bonds are probed by the solvation-induced red-shifts in the NH stretch modes. The three Ama⁺ isomers include the most abundant canonical cage isomer ( I ) produced by vertical ionization, which is separated by appreciable barriers from two bicyclic distonic iminium ions obtained from cage-opening (primary radical II ) and subsequent 1,2 H-shift (tertiary radical III ), the latter of which is the global minimum on the Ama⁺ potential energy surface. The effect of solvation on the energetics of the potential energy profile revealed by the calculations is consistent with the observed relative abundance of the three isomers. Comparison to the adamantane cation indicates that substitution of H by the electron-donating NH₂ group substantially lowers the barriers for the isomerization reaction.     

Formation of positive cluster ions LinBr (n = 2–7) and ionization energies studied by thermal ionization mass spectrometry

Clusters of the type Li_nX (X = halides) can be considered as potential building blocks of cluster‐assembly materials. In this work, Li_nBr (n = 2–7) clusters were obtained by a thermal ionization source of modified design and selected by a magnetic sector mass spectrometer. Positive ions of the Li_nBr (n = 4–7) cluster were detected for the first time. The order of ion intensities was Li₂Br⁺ > Li₄Br⁺ > Li₅Br⁺ > Li₆Br⁺ > Li₃Br⁺. The ionization energies (IEs) were measured and found to be 3.95 ± 0.20 eV for Li₂Br, 3.92 ± 0.20 eV for Li₃Br, 3.93 ± 0.20 eV for Li₄Br, 4.08 ± 0.20 eV for Li₅Br, 4.14 ± 0.20 eV for Li₆Br and 4.19 ± 0.20 eV for Li₇Br. All of these clusters have a much lower ionization potential than that of the lithium atom, so they belong to the superalkali class. The IEs of Li_nBr (n = 2–4) are slightly lower than those in the corresponding small Li_n or Li_nH clusters, whereas the IEs of Li_nBr are very similar to those of Li_n or Li_nH for n = 5 and 6. The thermal ionization source of modified design is an important means for simultaneously obtaining and measuring the IEs of Li_nBr (n = 2–7) clusters (because their ions are thermodynamically stable with respect to the loss of lithium atoms in the gas phase) and increasingly contributes toward the development of clusters for practical applications. Copyright © 2012 John Wiley & Sons, Ltd.     

Stable isomers for trifluoroacetic acid (TFA) pentahydrates obtained from density functional calculations

Stable isomers of trifluoroacetic acid (TFA) pentahydrate clusters, TFA-(H₂O)₅, have been explored by using density functional theory calculations. As done for TFA-(H₂O)₄ (Ito, 2013), structure optimization and vibrational calculations were performed for 70 isomeric structures (68 for neutral and 2 for ion-pair species, respectively) at the B971/6–311++G(3df,3pd) level. We found that the edge-sharing bicyclic isomer is at the global minimum and that three other isomers lie energetically within 100 cm⁻¹. Two types of ion-pair species were found to be unstable by 1100 cm⁻¹ in comparison with the global minimum. The results were compared with infrared spectra observed in nitrogen matrix.     

Theoretical study of structure and vibrational properties of MgnOn (n = 3–10) clusters

The structure and harmonic vibrations of Mg_nO_n (n = 3–10) clusters have been investigated using density functional theory. All structures are found to be cumulenic D_nh rings (equal bonds, alternating angles), with one intense out‐of‐plane mode and three infrared (IR)‐active degenerate modes, of which the highest one is extremely intense and increases asymptotically to 1000 cm^?1 for n = 10 at the B3LYP/6‐311++G(2d,2p) level. Comparisons with C_2n clusters show that B_nN_n and Be_nO_n clusters, the structure and bonding type for the Mg_nO_n clusters are consistent with those of the C_2n (n = 3, 5, 7,…) clusters B_nN_n(n = 3–10) and Be_nO_n(n = 3–10) clusters. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

An experimental and theoretical approach to the molecular structure of 3‐{[4‐(3‐Methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐hydrazono}‐1,3‐dihydro‐ indol‐2‐one

The title molecule, 3‐{[4‐(3‐methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐hydrazono}‐1,3‐dihydro‐indol‐2‐one (C₂₂H₂₀N₄O₁S₁), was prepared and characterized by ¹H NMR, ¹³C NMR, IR, UV–visible, and single‐crystal X‐ray diffraction. The compound crystallizes in the monoclinic space group P21 with a = 8.3401(5), b = 5.6976(3), c = 20.8155(14) Å, and β = 95.144(5)°. Molecular geometry from X‐ray experiment and vibrational frequencies of the title compound in the ground state has been calculated using the Hartree–Fock with 6‐31G(d, p) and density functional method (B3LYP) with 6‐31G(d, p) and 6‐311G(d, p) basis sets, and compared with the experimental data. The calculated results show that optimized geometries can well reproduce the crystal structural parameters, and the theoretical vibrational frequencies values show good agreement with experimental data. Density functional theory calculations of the title compound and thermodynamic properties were performed at B3LYP/6‐31G(d, p) level of theory. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012