Physics Research International

Physics Research International / 2008 / Article

Research Letter | Open Access

Volume 2008 |Article ID 879017 |

D. Rodr铆guez, D. Soto, E. Ram铆rez, A. Cruz, A. Santana, G. E. L贸pez, "Study of Magnesium Diboride Clusters Using Hybrid Density Functional Theory", Physics Research International, vol. 2008, Article ID 879017, 4 pages, 2008.

Study of Magnesium Diboride Clusters Using Hybrid Density Functional Theory

Academic Editor: Angel Rubio Secades
Received25 Sep 2007
Accepted31 Oct 2007
Published09 Dec 2007


Using hybrid density functional theory and a relatively large basis set, the lowest energy equilibrium structure, vibrational spectrum, and natural orbital analysis were obtained for magnesium diboride clusters [(MgB2)x for 饾懃 = 1 , 2 , and 3]. For comparison, boron clusters [ B 饾懃 f o r 饾懃 = 2 , 4 , a n d 6 ] were also considered. The MgB2 and (MgB2)2 showed equilibrium structures with the boron atoms in arrangements similar to what was obtained for pure boron atoms, whereas, for (MgB2)3 a different arrangement of boron was obtained. From the population analysis, large electron density in the boron atoms forming the clusters was observed.

1. Introduction

Materials composed of diborides have attracted the attention of researchers in recent years because of its superconductivity properties [1, 2]. In particular, it has been found that MgB2 behaves as a superconductor at relatively high critical temperature, 饾憞 c = 4 0 K . Although 饾憞 c for this material is not as high as other oxides with superconducting properties, its preparation is less expensive and it has a lower anisotropy, making smoother grain boundaries and hence increasing current flow. The origin of the superconducting properties of MgB2 is still a topic of debate [3]. However, the structure of MgB2 might have a significant importance in this property. This material shows a simple hexagonal AlB2-type structure typical of borides. It is formed of layers of boron similar to the basal plane of graphite, separated by closed-packed layers of magnesium. The magnesium atoms are located at the center of the hexagons formed by the boron atom and donate the electron density to the boron layers. It has been postulated [4] that this electron donation might be related to the superconducting nature of this material. Also, electronic structure calculations of bulk MgB2 have shown [5] a metallic character of the material which causes a strong electron-phonon interaction and hence superconducting properties. Experimental evidence has suggested [4] that an increment in the 饾憞 c temperature can be obtained if nanosized particles of MgB2 are embedded in inert materials. This increment in 饾憞 c seems to be caused by the change in the electronic density of states. The main objective of this study is to use cluster models to understand the electronic problem of MgB2. These models have been extensively used to understand the behavior of bulk materials because of the great theoretical detail that can be obtained [6, 7]. The lowest energy equilibrium structure was obtained and corroborated with a harmonic vibrational analysis. Natural orbital analysis was performed in order to characterize the electron density of the system. Boron clusters were also considered for comparison purposes with the results obtained being in agreement with recent studies [8].

2. Methods

The theoretical methods used in this study were based on ab-initio methods within the density functional theory approximation. Specific details of the procedure followed can be found elsewhere [9]. Specifically, the three parameter hybrid functional B3PW91 was used with the Dunning basis set cc-pvtz augmented with polarization functions. Hence, the method was termed B3PW91/aug-cc-pvtz. For all systems, a full geometry optimization was performed and various initial geometries were used to guarantee the determination of the lowest energy equilibrium structure. Harmonic vibrational analysis was performed for each system not only to obtain the vibrational frequencies, but also to characterize the nature of the structure obtained in the potential energy surface (PES). All calculations were performed using the Gaussian 03 software package [10] installed on an SGI 200 workstation. In order to analyze the charge distribution in the system, natural bond orbital (NBO) analysis was performed on the optimized wave function using the NBO program of the Gaussian 03 package. This analysis provides an improvement to the Mulliken population analysis usually used in the description of the charge distribution.

3. Results and Discussion

The equilibrium bond lengths of all the optimized structures are presented in Table 1 and the structures are shown in Figures 1 and 2. Figure 1(a) shows the lowest energy structure of MgB2. As in previous studies [11, 12] a nonlinear structure was obtained as the lowest energy equilibrium structure with the boron atoms at an equilibrium distance 0.09鈥壝 shorter than in the B2 molecule (see Figure 1(b)). From the NBO analysis shown in Figure 1(b), it can be observed that the electronic density is concentrated at the boron atoms, that is, the magnesium is donating electron density to the diborine. Table 2 shows the harmonic vibrational spectrum of MgB2 and B2 and it can be observed that all frequencies obtained have a positive sign which means that the global minimum was obtained. In the case of MgB2, the frequency that corresponds to the diboride vibration, frequency 1, is 16% higher than in the B2 molecule because the bond shortening. Figures 1(c) and 1(d) show the lowest energy for equilibrium structure of (MgB2)2 and B4, respectively. Both structures show a planar arrangement of atoms, with the boron atoms in a two-dimensional unsymmetrical rhombic arrangement for (MgB2)2 and a perfect square for B4. Specifically, two of the boron-boron bonds in the molecular species do not change and the other two were stretched by 9%. The magnesium atoms were accommodated at a relatively long distance when compared to the separation of the boron atoms. The vibrational frequency analysis shown in Table 2 demonstrates that the structure is a minimum in the PES. A direct comparison between the modes in the molecular species and the pure boron clusters was not possible. Figure 1(c) shows the NBO density analysis for the dimer and it can be observed that all the boron atoms forming the 2D rhombic arrangement have a negative charge. The boron atoms, which are closer to the Mg atoms, showed larger value of electron density than the other atoms which has the magnesium at a larger distance. Hence, the electron density is being concentrated on the boron atoms in a nonsymmetrical fashion and is being donated by the Mg atoms. Figure 2 shows the lowest energy equilibrium structure of (MgB2)3 and B6. Boron atoms form a nonsymmetrical nonplanar arrangement. Although the B6 do not show a planar structure, the arrangement of atoms is different to the one obtained in the molecular species. The vibrational analysis shows all positive frequencies and similar to the previous system, it is not possible to correlate the frequencies between the B6 and the molecular system. The NBO analysis for the molecular systems shows a similar electronic distribution to the one observed for the dimer, namely, the electronic distribution is located on the boron atoms. Two of the boron atoms show a larger electronic distribution (鈭0.853) because of their proximity to the magnesium atoms.


饾憻 1 = 1.551 饾憻 1 = 1.726 饾憻 1 = 1.545 饾憻 1 = 1.641 饾憻 1 = 1.522 饾憻 1 = 1.797
饾憻 2 = 2.197 饾憻 2 = 1.581 饾憻 2 = 1.585 饾憻 2 = 1.586
饾憻 3 = 2.157 饾憻 3 = 1.730 饾憻 3 = 1.573
饾憻 4 = 2.254 饾憻 4 = 1.700 饾憻 4 = 1.541
饾憻 5 = 2.247
饾憻 6 = 2.323
饾憻 7 = 3.352
饾憻 8 = 2.800


(1) 452(1) 64(1) 70(1) 962(1) 313(1) 184
(2) 496(2) 191(2) 138(2) 352(2) 305
(3) 1132(3) 248(3) 146(3) 1063(3) 346
(4) 253(4) 175(4) 1182(4) 381
(5) 324(5) 254(5) 1276(5) 560
(6) 428(6) 293(6) 1277(6) 563
(7) 448(7) 306(7) 713
(8) 578(8) 336(8) 838
(9) 606(9) 380(9) 1074
(10) 889(10) 431(10) 1120
(11) 1080(11) 460(11) 1261
(12) 1104(12) 489(12) 1320
(13) 573
(14) 581
(15) 696
(16) 725
(17) 819
(18) 973
(19) 1078
(20) 1193
(21) 1221

4. Conclusions

In all cases considered, the lowest energy equilibrium structure for MgB2 clusters showed an arrangement of boron atoms that are two dimensional or quasi-two dimensional. The electron density is concentrated on the boron atoms and hence the magnesium atoms are acting as electronic donating species. These results are in accordance to what is observed for MgB2 in the bulk [13]. If the superconductivity behavior of this material is increased by the larger electronic density in a specific plane, then it can be expected that this clusters should exhibit a similar superconducting behavior.


This work is supported by grants from the NIH-SCORE (Grant No. GM008103-03) and NIH-COBRE (Grant No. P2RR16439-01) programs.


  1. S. L. Bud'ko, G. Lapertot, C. Petrovic, C. E. Cunningham, N. Anderson, and P. C. Canfield, 鈥淏oron isotope effect in superconducting MgB2,鈥 Physical Review Letters, vol. 86, no. 9, pp. 1877鈥1880, 2001. View at: Publisher Site | Google Scholar
  2. J. Kortus, 鈥淐urrent progress in the theoretical understanding of MgB2,鈥 Physica C, vol. 456, no. 1-2, pp. 54鈥62, 2007. View at: Publisher Site | Google Scholar
  3. H. J. Choi, D. Roundy, H. Sun, M. L. Cohen, and S. G. Louie, 鈥淭he origin of the anomalous superconducting properties of MgB2,鈥 Nature, vol. 418, no. 6899, pp. 758鈥760, 2002. View at: Publisher Site | Google Scholar
  4. H. Y. Zhai, H. M. Christen, L. Zhang et al., 鈥淪uperconducting magnesium diboride films on Si with Tc0~24 K grown via vacuum annealing from stoichiometric precursors,鈥 Applied Physics Letters, vol. 79, no. 16, pp. 2603鈥2605, 2001. View at: Publisher Site | Google Scholar
  5. G. L. Zhao, A. Pullen, and D. Bagayoko, 鈥淭he metallic nature of boron layers in magnesium diboride,鈥 International Journal of Modern Physics B, vol. 17, no. 31-32, pp. 5905鈥5910, 2003. View at: Publisher Site | Google Scholar
  6. J. Zhao and R.-H. Xie, 鈥淐luster-assembled materials based on Na6Pb,鈥 Physical Review B, vol. 68, no. 3, Article ID 035401, 5 pages, 2003. View at: Publisher Site | Google Scholar
  7. X. Gu, S. Bulusu, X. Li et al., 鈥Au34-: a fluxional core-shell cluster,鈥 Journal of Physical Chemistry C, vol. 111, no. 23, pp. 8228鈥8232, 2007. View at: Publisher Site | Google Scholar
  8. A. N. Alexandrova, A. I. Boldyrev, H.-J. Zhai, and L.-S. Wang, 鈥淎ll-boron aromatic clusters as potential new inorganic ligands and building blocks in chemistry,鈥 Coordination Chemistry Reviews, vol. 250, no. 21-22, pp. 2811鈥2866, 2006. View at: Publisher Site | Google Scholar
  9. M. C. Vicens and G. E. L贸pez, 鈥淒ensity functional studies of cation-water complexes,鈥 Journal of Computational Chemistry, vol. 21, no. 1, pp. 63鈥68, 2000. View at: Publisher Site | Google Scholar
  10. M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., 2004, Gaussian Inc Wallingford CT.
  11. S. Zhang, P. Yang, and S. Li, 鈥淕eometries and properties of adenine-BH3 complex: an investigation with density functional theory (DFT) method,鈥 Journal of Molecular Structure: THEOCHEM, vol. 677, no. 1鈥3, pp. 161鈥166, 2004. View at: Publisher Site | Google Scholar
  12. E. P. F. Lee and T. G. Wright, 鈥淭he ground state and heat of formation of MgB2,鈥 Journal of Physical Chemistry A, vol. 108, no. 36, pp. 7424鈥7428, 2004. View at: Publisher Site | Google Scholar
  13. 艦. Erko莽, 鈥淪tructural and electronic properties of magnesium diboride,鈥 International Journal of Modern Physics C, vol. 14, no. 4, pp. 483鈥489, 2003. View at: Publisher Site | Google Scholar

Copyright © 2008 D. Rodríguez et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.