A Relativistic Density Functional Study of the U2F6 Molecule
Total Page:16
File Type:pdf, Size:1020Kb
Journal of Modern Physics, 2012, 3, 865-869 http://dx.doi.org/10.4236/jmp.2012.38113 Published Online August 2012 (http://www.SciRP.org/journal/jmp) A Relativistic Density Functional Study of * the U2F6 Molecule Yunguang Zhang School of Science, Xi’an University of Posts and Telecommunications, Xi’an, China Email: [email protected] Received May 26, 2012; revised June 24, 2012; accepted July 19, 2012 ABSTRACT All-electronic relativistic density functional theory (DFT) method has been used to study the U2F6 molecule. Results from calculations predict the existence of U2F6 molecule, which has been found to be stable with a multiply bonded U2 unit. The calculations also predict that D3d symmetry of U2F6 is more stable than D3h. The optimized geometries, vibra- tional frequencies and infrared intensities are reported for D3d symmetry of U2F6 from Becke Three-parameter Lee-Yang-Parr (B3LYP) function with triple zeta valence plus polarization functions basis set (TZP). The bond disso- ciation energy (BDE) for U-U bond in the U2F6 was obtained using the same method. In addition, the entropies of U2F6 have been investigated at temperature rang from 0 to 3000 K in 10 steps using the B3LYP method. Keywords: Molecular Structure; Vibrational Frequency; Bond Dissociation Energy 1. Introduction try could occur in solution, analogous to the widely studied transition metal dimmer chemistry. In theory, Cl U Cl The uranium compounds bring a prodigious challenge for 3 2 3 molecule has been studied by Complete-Active-Space investigation from theory and experiment. Still, uranium compounds are not easy to manage in the laboratory. It is Self-Consistent-Field (CASSCF)/Complete-Active-Space very important that we could gain a good understanding Second-Order-Perturbation Theory (CASPT2) method, of their chemical properties in many areas. Quantum the results from calculation predict to be stable with a chemistry is now mature enough to compute uranium multiply bonded U2 unit [4]. In the paper, All-electronic compounds properties with good accuracy. Gagliardi and DFT has been used to show that diuranium compound Roos have explored how two uranium atoms bind to F3U2F3 is stable with a multiple U-U bond. The calcula- tions also predict that D symmetry of U F is more sta- form the U2 molecule by multiconfigurational quantum 3d 2 6 ble than D symmetry. The optimized geometries, BDE chemistry methods [1,2]. They find that the U2 chemical 3h bond is more complex than any other known diatomic and vibrational frequencies are reported for D3d symme- bond, with summation of the bonding electrons suggesting try of U2F6 from B3LYP function with TZP basis set [5]. a quintuple bond. In other words, the uranium atom has In addition, the entropies of U2F6 have been investigated the ground-state electronic configuration 567f 312ds at temperature rang from 0 to 3000 K in 10 steps using (six valence electrons), and the U-U bond of U2 molecule the same method. is composed of three normal two-electron bonds, four Methods of calculation are given in the next section, electrons in different bonding orbitals, and two non- and the Result and Discussion are then presented in Sec- bonded electrons leading to a quintuple bond between the tion 3. The article ends with the Summary. two uranium atoms. In addition, its multi-radical nature may support chemical bonds to the framework of diura- 2. Computational Details nium chemistry, such as in F3U2F3. Our calculations were performed using the Amsterdam Not much diuranium compounds have been found ex- density functional (ADF) program. The underlying the- perimentally. So far, only the molecules OU2O and H2U2H2 ory of the program is the Kohn-Sham approach to the [3] have been detected. However, no real inorganic chem- DFT. Kohn-Sham DFT is an important computational istry based on the U2 unit has been observed. This chemis- method to predict chemical properties accurately and to analyze and interpret these in convenient and simple *Project supported by the young teacher scientific research foundation of Xi’an University of Posts and Telecommunications (Grant No. chemical terms [6]. DFT gives superior accuracy to Har- 0001294). tree-Fock theory and semi-empirical approaches, and it is Copyright © 2012 SciRes. JMP 866 Y. G. ZHANG well suited for molecules containing heavy metal atoms. B3LYP level of theory using the TZP basis set. The most In this method, scalar relativistic effects are introduced stable spin-free electronic state was found to be a singlet via the zeroth-order regular approximation (ZORA) ap- state. The U-U bond distance is 2.425 Å, the U-F dis- proach to the Dirac equation [7-9], and that spin-orbit tance is 1.998 Å, and the U-U-F angle is 120.0˚C In ad- coupling is not included in the present study. It will have dition, we know that scalar relativistic effects are intro- only a minor effect on the computed structure for singlet duced via the ZORA to the Dirac equation, and that spin- electronic state. Certainly, the spin-orbit coupling effect orbit coupling is not included in the present study. Cer- will further lower the energy but has a negligible effect tainly, the spin-orbit coupling effect will further lower on the geometry. For the close-shell systems, the spin the energy but has a negligible effect on the geometry. restricted DFT scheme is used. In this paper, a triple-ζ In the U2F6, the formal charge of uranium atom is +3. Slater type basis set augmented with one set of polariza- Thus, 6 of 12 valence electrons remain, and a triple bond tion functions, i.e. the triple-zeta polarized(TZP) basis set can be formed. The electronic configuration for this is employed for all atoms(F and U). Moreover, the cal- molecule including all the electrons on U and F atoms is: culations were performed using GGA level to [10,11] 2424 2 DFT: B3LYP function in all-electronic method. 2AA12ug 2 18 E 1 u 18 E 1 g 20 A 2 u (1) ADF has a fragment oriented approach: the poly- 4 2 19EA 21 atomic system to be computed is conceptually built up 11ug from fragments, the molecular orbitals are calculated as where the highest occupied molecular orbitals of each linear combinations of fragment orbitals, and final analy- symmetry are indicated. sis of the bonding energy are in terms of fragment prop- erties. The fragments may be single atoms or larger 3.2. Bonding Energy Analysis of F3U-UF3 moieties. When we compute a system in terms of its con- stituent fragments, these fragments must have been com- The energy decomposition between the two radical mo- puted before and their properties must be passed on to lecular fragments UF3 in the compound F3U-UF3 have the current calculation. This is done by attaching frag- been analyzed with the program package ADF which is ment files, which contain the necessary information. In found on the methods of Morokuma [12,13] and Ziegler and Rauk [14,15]. The overall bond energy E be- the paper, the BDE for U-U bond in the U2F6 between the e tween the two fragments UF3 is given by the summation F3U and UF3 radical fragments were obtained using the fragment approach. of Eprep and Eint . 3. Results and Discussion 3.1. Structure and Stability In a previous study of U2 molecule [1,2], it was shown 2 4 that the ground state configuration is 76sd guπ 111111 6655ddfgggu fπ 55ffug and that the chemical bond is thus more complex than any other known diatomic bond, with summation of the bonding electrons suggesting a quintuple bond. Its multi- radical nature (six electrons are available for binding) Figure 1. Staggered conformation of U2F6. may support chemical bonds to F3U2F3 molecule. Here we have considered the diuranium compound of U2F6. It can have either an eclipsed or a staggered conformation (see Figures 1 and 2). Preliminary calculations have in- dicated that the staggered conformation is lower in en- ergy than the eclipsed form, and the calculations were done under limiting symmetry. Moreover, when the mo- lecular structure of U2F6 was optimized under D3d and D3h symmetry respectively (unlimiting symmetry), the final optimum structure was found to be staggered. So the stability of D3d symmetry is better than D3h. We will thus focus our analysis on the staggered structure. Optimization of the geometry was performed at the Figure 2. Eclipsed conformation of U2F6. Copyright © 2012 SciRes. JMP Y. G. ZHANG 867 EEe prep E int (2) repulsive term EPauli , the attractive electrostatic term Eelstat , orbital interaction term Eorb , and preparation The preparation energy Eprep is the sum total of en- energy term Eprep are 250.20 kcal/mol, –63.37 kcal/ ergy which is necessary associated to the structure distor- mol, –205.33 kcal/mol and 1.24 kcal/mol, respectively. tion of both fragments from their equilibrium geometry From the Equations (2) and (4), the calculational bond to the geometry they have when the complex is formed. dissociation energy of U-U bond is 17.26 kcal/mol. In this EE(distorted fragments) molecule, there is a net flow of electrons from the metal prep total E (fragments in the equilibrium structure) to substituent F3, resulting in a positive atomic charge on total U according to both the VDD (0.270) and the Hirshfeld (3) (0.722) [19,20]. In addition, we also give the breakdown Eint in Equation (2) is the instantaneous interaction of the E term into contributions of E (–340.48 orb A1 energy between the two distorted fragments in the mole- kcal/mol), E (–2.65 kcal/mol) and E (137.83 A2 E1 cule.