Prediction of stable radon fluoride molecules and geometry optimization using first

Sep 29, 2023

Scientific Reports volume 13, Article number: 2898 (2023) Cite this article

804 Accesses

3 Citations

4 Altmetric

Metrics details

Noble gases possess extremely low reactivity because their valence shells are closed. However, previous studies have suggested that these gases can form molecules when they combine with other elements with high electron affinity, such as fluorine. Radon is a naturally occurring radioactive noble gas, and the formation of radon-fluorine molecules is of significant interest owing to its potential application in future technologies that address environmental radioactivity. Nevertheless, because all isotopes of radon are radioactive and the longest radon half-life is only 3.82 days, experiments on radon chemistry have been limited. Here, we study the formation of radon molecules using first-principles calculations; additionally, possible compositions of radon fluorides are predicted using a crystal structure prediction approach. Similar to xenon fluorides, di-, tetra-, and hexafluorides are found to be stabilized. Coupled-cluster calculations reveal that RnF6 stabilizes with Oh point symmetry, unlike XeF6 with C3v symmetry. Moreover, we provide the vibrational spectra of our predicted radon fluorides as a reference. The molecular stability of radon di-, tetra-, and hexafluoride obtained through calculations may lead to advances in radon chemistry.

A noble gas with a fully filled outer shell is not reactive. When argon was first discovered, chemists found that it was not reactive with the other elements in the periodic table, and noble gases were generally considered to be unreactive. In the 1930s, Pauling1 predicted that xenon (Xe) would be able to form compounds with fluorine. In related experiments, they succeeded only in corroding the walls of a quartz flask and were unaware of the presence of new compounds in it2. After several trials, xenon and fluorine readily reacted and formed the solid XeF4, which is stable even at room temperature3. The structures of xenon difluoride (XeF2) and xenon tetrafluoride (XeF4) were identified using their vibrational spectra4. However, the structure of xenon hexafluoride (XeF6) is controversial regarding the stereo activity of valence electron lone pairs5. Experimental evidence obtained using electron diffraction6 and vibrational spectroscopy7 suggests that XeF6 forms a distorted octahedral symmetry8.

Radon (Rn) is a noble gas like Xe; it is a naturally occurring radioactive material (NORM) found in the underground environment. It has a high density (9.73 g/L at standard temperature and pressure) and solubility (230 cm3/L at 20 °C)9 among the noble gases. When gaseous Rn is inhaled, it either emits alpha rays directly or decays into daughter radionuclides, which may cause lung cancer10. In addition, owing to its short half-life (only 3.82 d), experiments and studies involving radon have been limited thus far. Heavy noble gases, such as Xe and Rn, are both rare elements and highly radioactive gases; therefore, they present various challenges to experimental researchers. Radon difluoride (RnF2) was synthesized by Fields11; it formed nonvolatile RnF2 when exposed to fluorine and heated to 400 °C. RnF2 is currently the only known radon molecule. Owing to high temperature and pressure requirements, obtaining a complete solid RnF2 sample in a natural environment is challenging.

Computational studies based on first-principles density functional theory (DFT) are regarded as a practical alternative for the study of these gases as DFT has been successfully applied for the detailed analysis of a vast number of materials. The DFT approach can be used to provide qualitative predictions of geometric features and determine various chemical and physical properties. However, computational studies based on first-principles about interactions of Rn atoms with environments (e. g., Rn·H2O complex)12 or formation of Rn molecules are very rare.

A previous study13 used redundancy analysis to determine the Rn content in groundwater, and they found a positive correlation with fluorine concentration. This raised the questions of whether Rn can interact with fluorine, and how the radon fluoride reaction should be elucidated. Only a few theoretical studies14 have demonstrated the relationship between Rn and fluorine, but most of this literature is focused on Xe or Kr, and Rn was not treated as the main subject. Malli15 and Filatov16 suggested the atomization energy and RnF bond length of radon hexafluoride (RnF6) at the Hartree–Fock (HF) and MP2 levels of theory. However, they assumed only octahedral RnF6 geometry as the initial structure and did not consider other possibilities. Our study aims to understand the possible chemical bonding between Rn and fluorine and its geometry through first-principles calculations. Additionally, the vibrational frequencies of radon di-, tetra-, and hexafluorides were calculated at the coupled-cluster level, suggesting as a reference. Herein, first-principles calculations were performed for Xe and Rn fluorides with various basis sets at different methods to compare geometrical parameters. Subsequently, the most stable combination describing the Xe and Rn fluorides was provided for geometry optimizations.

Gillespie17 and Hedberg18 predicted that the geometry of ground state XeF6 is a distorted octahedron with C3v point symmetry. According to prior studies, two characteristic bond lengths of Xe-F are experimentally determined to be 1.85 and 1.94 ± 0.036 Å19, and the∠F-Xe-F of 114.9° and 81.0° were obtained at the self-consistent field (SCF) level20. To calculate the geometries of XeF6, initial symmetries were set to C2v, C3v, and octahedral (Oh); subsequently, the bond lengths and angles were fully relaxed. These three symmetries of XeF6 were calculated using the DFT, Møller–Plesset second-order perturbation theory (MP2)21, and coupled cluster singles and doubles (CCSD)22,23,24,25 methods. The calculated structures by DFT and MP2 are described in Supplementary Tables S1 and S2 of the Supplementary Information (SI).

Certain calculated structures converged equivalently to the Oh molecular geometry, even though they started at different initial point groups such as C2v or C3v, according to the DFT and MP2 results. Other calculated structures stay at their initial symmetry and they tended to slightly overestimate the bond length compared with the experimentally obtained values by ~ 0.06 Å. Note that bond lengths determined using the MP2 method were rather closer to the experimental values compared with those provided by the DFT method. The angles obtained from these two methods tend to be underestimated compared with other calculated values, and the C3v structure was stable at the lowest energies for all basis sets, except in the case of the Oh-converged structure. The energy differences for other symmetries compared with that of C3v are denoted as ΔE. The bond angles calculated using the DFT and MP2 methods are deviated noticeably from other calculated values up to 9.55%.

In case DFT and MP2 gave the same answer, we additionally applied CC method as a cross-check to validate our calculation results. The XeF6 structure was determined at the CCSD level, as presented in Supplementary Table S3. At the CCSD level, the relativistic effect of Xe was not considered, and three initial symmetries were the same as those in the DFT and MP2 results. The initial symmetries were maintained during the relaxation and C3v symmetry is predicted to be the ground state for all basis sets. For the LANL2DZ and CEP-31G basis sets, the calculated lengths have a good agreement with experimental data, with a small difference of the order of two decimal places in angstrom. The angles for the C3v symmetry were different from other calculated values only by approximately 1° for these basis sets.

In summary, the bond lengths and bond angles determined using the DFT and MP2 methods were overestimated and underestimated respectively, implying stabilizing the metastable stereo-isomers was challenging. When it comes to XeF6, CCSD level calculations with non-relativistic effect gave quantitative agreement with experimental data, which is consistent with several previous studies26,27.

In the case of many-electron atoms, the relativistic contraction of inner-shell orbitals by screening affects the outer-shell orbitals28. This may result in a significant impact on the chemical and physical properties of heavy inert gases in the lower half of the periodic table29. Therefore, the relativistic effect of Xe was further considered at the CCSD level using the DKH Hamiltonian. Supplementary Table S4 summarizes the bond length and total energies of XeF6 for the three geometries with DK3 basis sets. The relativistic effect did not significantly affect geometric parameters. The C3v structure possesses lower total energy than the Oh structure and by C2v structure by 18.88 kcal/mol and 8.28 kcal/mol, respectively using the cc-pVTZ-DK3 basis set. For aug-cc-VTZ-DK3, the overall energy differences, \(\Delta \mathrm{E}\) were lower than those of cc-VTZ-DK3, although there was a tendency to overestimate the bond length by approximately 2.06%. In the DKH calculation, the angle was reduced by approximately 5º compared with the case employing the LANL2DZ basis set in the NR calculation. Although the NR calculation produced values that were in better agreement with the experimental results for XeF6, the DKH calculation also exhibited small differences that were within 0.01 Å. Even when relativistic effects were not considered, the CCSD level was found to be consistent with the experiment for XeF6 (C3v). Our calculation of XeF6 parameters can be used for the following computational study to determine the effective computational level depending on initial symmetry in Rn-F chemistry.

The atomization energies of argon, krypton, and xenon difluoride indicate that these molecules tend to become stable with increasing atomic number30. Following this trend, radon difluoride is expected to be more stable than the corresponding xenon molecules31. Therefore, the simpler structures of radon difluoride (RnF2) and radon tetrafluoride (RnF4) were first determined. However, the only reported experimental values for radon molecules are those for RnF2, which were proposed by Fields11. Therefore, first, the stability was investigated by calculating the formation energies of radon fluorides. The convex hull diagram in Supplementary Fig. S1 indicated that the lowest energy can occur between the Rn and binary F2 phases. The energy convex hull graph helps to understand that RnF6 is the most stable molecule along the thermodynamic route while RnF2 and RnF4 may exist as metastable forms. Speculatively, as there is an experimental observation of RnF2, we can easily speculate the possible formation of RnF4 and RnF6 as well. For structural optimization, the geometry of RnF4 and RnF6 were determined by CCSD, with relativistic effects (DKH). As presented in Table 1, the bond lengths of RnF2 and RnF4 were estimated to be 2.04–2.05 Å and 2.00–2.01 Å, respectively, depending on the basis sets employed. The optimized structures of RnF2 and RnF4 were obtained as linear and square planar, as in the case of xenon fluoride, respectively. The bond lengths of radon fluorides were longer than those of xenon fluorides, which is a predictable result considering the larger atomic size of radon. In addition, radon di- and tetrafluoride retained a similar structure to xenon fluorides.

The symmetrically optimized structures and energies for C2v, C3v, and Oh obtained using DFT are summarized in Supplementary Table S5. Several previous studies15,16,20 predicted the bond lengths of RnF6, but they exhibited large differences in bond length depending on the type of calculation employed. For the DFT simulations, all basis sets converged equivalently to Oh, and the bond lengths of RnF6 differed by up to 0.09 Å compared with the experimental XeF6 structures19. In general, the stronger the attraction between bonding atoms, the shorter the bond, and the greater the atom size, the longer the bond. As expected, the bond lengths of RnF6 were greater than those of xenon, although the difference between these was only 0.09 Å. As Rn has an atomic radius that is approximately 0.12 Å larger than that of Xe, the calculated bond lengths for radon fluoride were predicted to be shorter, as expected. Using the enthalpy of RnF2 formation14,32, Rn was theoretically demonstrated to form relatively strong covalent bonds compared with other noble gases33. In that study, Rn was expected to form shorter and more stable bonds with fluorine than XeF6, and these results are consistent with our calculations.

The calculated structures and relative energies for RnF6 for the MP2 method are described in Supplementary Table S6. For all converged Oh structures, the RnF6 bond length varied between 1.98 and 2.03 Å, and these values tended to be slightly underestimated compared to the DFT results.

Unlike XeF6, all converged geometries exhibited Oh symmetry in the case of RnF6, irrespective of initial symmetry. Although stabilization of the metastable C2v and C3v state was attempted, this was only possible using a fixed bond angle, as shown in Fig. 1. These results reveal that the energy barrier between the Oh ground state and the other metastable states (C3v and C2v) was particularly small, and therefore, the initial states easily relaxed into the Oh symmetry state. The bond lengths and total energies that were fully relaxed at the CCSD level are presented in Table 2. The radon fluoride bond lengths obtained using CCSD were slightly shorter for all basis sets compared with the other two methods (DFT and MP2). Unlike the DFT and MP2 results, the CEP-31G basis set was used to predict that C3v symmetry was the most stable structure at the CCSD level, and this differed from the second most stable structure of C2v by approximately 5 kcal/mol.

Comparison of the total energy (keV) for RnF6 calculated by the CCSD method. aanglefix is the fixed given angles and f symmfollow is the freely relaxed bond length.

Additionally, DKH calculations were conducted considering the relativistic effect on the structural parameters of Rn. Table 3 summarizes the bond length and total energies of RnF6 for the three initial geometries (C2v, C3v, and Oh). All initial geometries with different basis sets were converged to Oh, and the bond lengths were estimated to be approximately 1.97 Å. Note that the relativistic effect on heavy atoms destabilizes the d and f orbitals because the inner (core) orbitals of the atom are strongly attracted to the nucleus. However, the outer orbitals are subject to minimal relativistic effects because the effective nuclear charge felt by the electrons decreases due to the screening effect. This usually causes the s and p orbitals to contract further and the d and f orbitals to expand, thereby stabilizing the 6p orbitals slightly but the 6s orbital strongly for Rn16. Therefore, this effect reduces the covalent radius of Rn and consequently induces a shorter radon fluoride bond length16, following the typical trend of bond shrinkage in molecules15. The relativistic effect leads to the presence of shorter bond lengths because of the stabilization of Rn atom inner orbitals.

Also, it is worth to note that the ground-state geometry of XeF6 is rather different from that of RnF6 when computed using the same methods and basis sets. To compare relativistic effects between Xe and Rn molecules, the molecular electrostatic potential (MEP) of RnF6 and XeF6 was mapped, as shown in Fig. 2. The molecular electrostatic potential surfaces describe the charge distribution of molecules in three dimensions34. Figure 2a shows a sigma hole (σ-hole), a region with a positive surface electrostatic potential35 in the MEP of XeF6(C3v). The lone pair of Xe experiences an electron shielding effect that shields a partial positive charge in the Xe nucleus, inducing a σ-hole due to the repulsive force between the Xe lone pair and the three F atoms36. In contrast, RnF6 has evenly distributed σ-holes in the middle of each face of the octahedral37, as shown in Fig. 2b. The bond lengths of RnF6 calculated using NR and DKH calculations and the CCSD method are shown in Fig. 3. Evidently, RnF6 converged equivalently to Oh symmetry for all basis sets, except for CEP-31G, even in the case of the DKH calculations.

Molecular electrostatic potential (MEP) surface mapped (a) XeF6 (C3v) and (b) RnF6 (Oh).

The bond lengths of RnF6 including all calculations by the CCSD method. *The calculated values29 with relativistic effect (Rel) and non-relativistic effect (NR).

Infrared (IR) vibrational spectra were obtained at the CCSD level, considering the relativistic effect, for the optimized structures of radon fluorides with different stoichiometries. The calculated vibrational spectra obtained with the aug-cc-pVTZ-DK3 basis set are tabulated in Table 4, wherein the calculated IR spectra are plotted against the wavenumbers. As no studies conducted thus far have reported the IR spectra of Rn molecules, these could not be compared with experimentally obtained frequencies. For RnF6, the all given initial structures converged to Oh symmetry, and therefore, a frequency comparison for each structure was meaningless. In the case of XeF6, the C2v structure had one imaginary frequency, the Oh structure had a three-fold degenerate imaginary frequency, whereas Rn had no imaginary frequency. The presence of an imaginary frequency implies that the optimized structure under investigation is not stable26. This suggests that the optimized structures of radon fluoride were harmonically stable on the potential energy surface of the molecules. Note that the calculated frequencies at 648.28 cm−1, consistent with the eu (double degeneracy) mode, accounted for most of the IR active frequency of RnF4. The calculated vibrational spectra for RnF6 indicated that the t1u bending mode (at 663.93 cm−1) contributed to the greatest extent to the intensity. As RnF4 and RnF6 have not been characterized experimentally, our predicted vibrational spectra can play as a critical reference for future studies.

The dissociation energy of radon fluorides was predicted at CCSD levels with DKH calculations. Previous studies described the dissociation energies of xenon fluoride, which were obtained using classical thermodynamic equilibrium measurements with predicted entropies38. Likewise, this study described the thermodynamic reaction and dissociation energy of MF6 (M = Xe or Rn) using the following expressions.

Previous studies have described the heat of formation calculations for XeF6, and obtained – 62.1 ± 1.4 kcal/mol at 0 K and − 64.0 ± 1.4 kcal/mol at 298 K26. In this study, the dissociation energy of XeF6 (C3v) → Xe + 3F2 for the most stable structure, which was obtained consistently with the measured values, was 5.04 kcal/mol smaller than the experimentally obtained equilibrium value38, and 26.94 kcal/mol lower than that from the photoionization experiment39, as shown in Supplementary Fig. S2. Several previous studies reported errors in the experimental equilibrium values of 1–2 kcal/mol, which depended on structural differences8. However, the difference in the total dissociation energy in this study was relatively large, depending on symmetry. For example, the total dissociation energy of the CF4+ ion exhibited a variation of approximately 50 kcal/mol according to structural differences at the same level owing to the Jahn–teller effect40. Therefore, providing experimental evidence that addresses the question of energy differences according to XeF6 structures is challenging. The calculation results for XeF6 (C3v) using aug-cc-pVTZ-DK3A may be reasonable as these are rather similar to the equilibrium experimental values.

Dissociation energy was calculated according to various reactions of radon fluoride, and the results are presented in Supplementary Table S7; additionally, the dissociation energy of RnF6 → Rn + 3F2 is shown in Fig. 4. In general, the loss of F2 from RnF2 is slightly harder than the F2 loss from RnF4, and the loss of F2 in RnF4 is much harder than that in RnF6. This is consistent with the presence of a growing steric community as more fluoride is added to the central Rn8. Because convergence to Oh was achieved for all basis sets, no difference in the dissociation energy between Rn molecule structures was observed. There are no significant differences between previous studies16 with the MP2 method and our calculated values (aug-cc-pVTZ-DK3A: 139.62 kcal/mol and cc-pVTZ-DK3A: 122.85 kcal/mol, see Fig. 4).

The dissociation energy of RnF6 → Rn + 3F2 with relativistic calculation by CCSD method. *The calculated values16 with relativistic effect (Rel) and non-relativistic effect (NR).

Despite the necessity of studying the behavior of radioactive materials, such as Rn, which naturally occur in the environment, only radon difluoride (RnF2) and its complexes are known. Both a short half-life and radiological risks challenge the experimental assessment of Rn, and the existence of Rn molecules remains controversial. This study described the molecular geometry, vibrational spectra, and dissociation energies of the RnF2, RnF4, and RnF6 family of molecules using ab initio calculations.

Detailed parameters employed to describe xenon fluoride were used to obtain an optimized computational approach. For the XeF6 structure, the geometric parameters obtained at the CCSD level indicated a difference between calculated bond lengths of approximately 0.98% and angles that differed from other calculated values by only 1°. In addition, the dissociation energies for Xe + 3F2 were in good agreement with the experimental38,39 and other calculation references26. Using these calculations as a base, the geometry, dissociation energies, and vibrational spectra of RnF2, RnF4, and RnF6 were predicted.

CCSD-level calculations with relativistic effects provided radon fluoride bond lengths of 2.04 and 2.00 Å for radon di- and tetrafluoride, respectively. Unlike XeF6, which maintained C3v symmetry, all possible initial group symmetries of RnF6 converged to an octahedral structure, with estimated bond lengths of approximately 1.97 Å. The dissociation energies for Rn and 3F2 were approximately 139.62 kcal/mol, which is higher than that of Xe. This may be assumed to be an accurate value, considering the high reactivity of Rn. This study can be extended by adding other anionic elements that can be possibly encountered in nature to the studied system, such as oxygen and chlorine, particularly in an aqueous environment. In addition, this study was the first to address the vibrational frequencies of radon di-, tetra-, and hexafluoride, and this information can be used as a basis for future radon-based experiments.

The initial molecular structures and formation energies of radon di-, tetra-, and hexafluoride were obtained via ab initio random structure searching (AIRSS)41 using the projected augmented plane-wave method (PAW) embedded in the Vienna ab initio simulation package (VASP)42,43,44. For structure searching, we placed a supercell with a volume of 10 × 10 × 10 Å3 and randomly generated Rn and F atoms for various Rn/F ratios. We searched ~ 5000 trial structures and 17 compositions. For formation energy calculations, the Perdew−Burke−Ernzerhof (GGA-PBE) exchange-correlation functional45 was used with a plane-wave basis set energy cutoff of 550 eV and the self-consistent-field energy convergence criteria of 1E−6 eV. Only a gamma-point k-mesh and a cubic unit cell with a lattice constant of 15 Å were used in this study.

All ab initio calculations for structural optimization and vibrational spectra were performed using the Gaussian09 simulation package46. The optimized molecular geometries of Xe and Rn fluorides were obtained using DFT and two different levels of the post Hartree–Fock (HF) theory method: MP221, and the coupled-cluster (CC) method. In this study, the B3LYP hybrid function, which is one of the hybrid functionals provided by the combination of Becke's three-parameter functional47 with the Lee, Yang, and Parr correlation functional48, was used for DFT calculations. In addition, MP2 was developed by Moller and Plesset21 to address HF theory in the case of many-electron systems. This approach was used in this study because it is known to correct minor deficiencies in the HF method and has been widely applied to molecular energy correlation. The CC method is one of the most prevalent post-HF methods and applies to sufficiently large molecules. A feature of the CC method is that, as opposed to other calculation methods, it can be systematically improved by including a higher excitation operator49.

The CCSD22,23,24,25 was used in our study, with higher excitations approximated by a product of lower excitations. Further, Gaussian09 was used at all levels in these calculations, with an SCF energy cutoff of 1E–6 Hartree.

Three basis sets were used for Xe and Rn. For nonrelativistic calculations, all methods (DFT, MP2, and CCSD) were used with a variety of basis sets, namely LANL2DZ50, def2-SVP51, and CEP-31G52,53,54 for Xe. Furthermore, the aug-cc-pVDZ-PP55, def2-SVP51, and CEP-31G52,53,54 basis sets were defined for Rn. For comparison, using the same basis sets for both elements produced better results. However, LANL2DZ could not be used for radon because the largest atom was defined as bismuth (Bi).

Relativistic effects were considered within the 2nd-order Douglas–Kroll–Hess (DKH) Hamiltonian implemented in the Gaussian09 software. Although this method does not involve spin–orbit coupling (SOC), the inclusion of SOC usually does not alter the molecular geometry qualitatively significantly56. The aug-cc-pVTZ-DK3 and cc-pVTZ-DK3 basis sets were selected for Xe and Rn, which is the only possible basis set both can be used for Xe and Rn. In addition, the effective core potential (ECP) was employed with both basis sets used to describe Rn. These extended ECP valence basis sets are denoted as "DK3A". In each calculation that involved NR and DKH, the 6-31G++57 basis set with polarization and diffuse functions was used for fluorine in the molecules.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Pauling, L. The formulas of antimonic acid and the antimonates. J. Am. Chem. Soc. 55, 1895–1900. https://doi.org/10.1021/ja01332a016 (1933).

Article CAS Google Scholar

Labinger, J. Why isn't noble gas chemistry 30 years older? The failed (?) 1933 experiment of Yost and Kaye. Bull. Hist. Chem. 40, 29–36 (2015).

CAS Google Scholar

Claassen, H. H., Selig, H. & Malm, J. G. Xenon tetrafluoride. J. Am. Chem. Soc. 84, 3593–3593 (1962).

Article CAS Google Scholar

Hyman, H. H. The chemistry noble gas compounds. Science 145, 773–783. https://doi.org/10.1126/science.145.3634.773 (1964).

Article ADS CAS PubMed Google Scholar

Seppelt, K. Nonoctahedral structures. Acc. Chem. Res. 36, 147–153. https://doi.org/10.1021/ar020052o (2003).

Article CAS PubMed Google Scholar

Bartell, L. S. & Gavin, R. M. Molecular structure of XeF6. II. Internal motion and mean geometry deduced by electron diffraction. J. Chem. Phys. 48, 2466–2483. https://doi.org/10.1063/1.1669471 (1968).

Article ADS CAS Google Scholar

Claassen, H. H., Goodman, G. L. & Kim, H. Spectral observations on molecular XeF6: Raman scattering and infrared, visible and ultraviolet absorption in the vapor and in matrix isolation. J. Chem. Phys. 56, 5042–5053. https://doi.org/10.1063/1.1676986 (1972).

Article ADS CAS Google Scholar

Dixon, D. A., de Jong, W. A., Peterson, K. A., Christe, K. O. & Schrobilgen, G. J. Heats of formation of xenon fluorides and the fluxionality of XeF6 from high level electronic structure calculations. J. Am. Chem. Soc. 127, 8627–8634. https://doi.org/10.1021/ja0423116 (2005).

Article CAS PubMed Google Scholar

Greenwood, N. N. & Earnshaw, A. Chemistry of the Elements (Elsevier, 2012).

Google Scholar

Samet, J. M. Radon and lung cancer. JNCI 81, 745–758 (1989).

Article CAS PubMed Google Scholar

Fields, P. R., Stein, L. & Zirin, M. H. Radon fluoride. J. Am. Chem. Soc. 84, 4164–4165. https://doi.org/10.1021/ja00880a048 (1962).

Article Google Scholar

Lee, E. P. F. & Wright, T. G. Interaction energy of the radon-water (Rn, H2O) complex. J. Phys. Chem. A 103, 7843–7847 (1999).

Article CAS Google Scholar

Kang, J., Singh, B. K. & Um, W. Efficient radon removal using fluorine-functionalized natural zeolite. J. Environ. Radioact. 233, 106607. https://doi.org/10.1016/j.jenvrad.2021.106607 (2021).

Article CAS PubMed Google Scholar

Liao, M.-S. & Zhang, Q.-E. Chemical bonding in XeF2, XeF4, KrF2, KrF4, RnF2, XeCl2, and XeBr 2: From the gas phase to the solid state. J. Phys. Chem. A 102, 10647–10654 (1998).

Article CAS Google Scholar

Malli, G. L. Relativistic all-electron Dirac-Fock calculations on RnF6 and its ions. J. Mol. Struct. (Thoechem.) 537, 71–77. https://doi.org/10.1016/S0166-1280(00)00663-1 (2001).

Article CAS Google Scholar

Filatov, M. & Cremer, D. Bonding in radon hexafluoride: An unusual relativistic problem?. Phys. Chem. Chem. Phys. 5, 1103–1105. https://doi.org/10.1039/B212460M (2003).

Article CAS Google Scholar

Gillespie, R. J. & Hargittai, I. The VSEPR Model of Molecular Geometry (Courier Corporation, 1991).

Google Scholar

Hedberg, K., Peterson, S. H., Ryan, R. R. & Weinstock, B. On the structure of gaseous XeF6. J. Chem. Phys. 44, 1726–1726. https://doi.org/10.1063/1.1726928 (1966).

Article ADS CAS Google Scholar

Pitzer, K. S. & Bernstein, L. S. Molecular structure of XeF6. J. Chem. Phys. 63, 3849. https://doi.org/10.1063/1.431880 (1975).

Article ADS CAS Google Scholar

Kaupp, M., van Wüllen, C., Franke, R., Schmitz, F. & Kutzelnigg, W. The structure of XeF6 and of compounds isoelectronic with it. A challenge to computational chemistry and to the qualitative theory of the chemical bond. J. Am. Chem. Soc. 118, 11939–11950. https://doi.org/10.1021/ja9621556 (1996).

Article CAS Google Scholar

Møller, C. & Plesset, M. S. Note on an approximation treatment for many-electron systems. Phys. Rev. 46, 618–622. https://doi.org/10.1103/PhysRev.46.618 (1934).

Article ADS MATH Google Scholar

Purvis, G. D. & Bartlett, R. J. A full coupled-cluster singles and doubles model: The inclusion of disconnected triples. J. Chem. Phys. 76, 1910–1918. https://doi.org/10.1063/1.443164 (1982).

Article ADS CAS Google Scholar

Scuseria, G. E., Janssen, C. L. & Schaefer, H. F. An efficient reformulation of the closed-shell coupled cluster single and double excitation (CCSD) equations. J. Chem. Phys. 89, 7382–7387. https://doi.org/10.1063/1.455269 (1988).

Article ADS CAS Google Scholar

Scuseria, G. E. & Schaefer, H. F. Is coupled cluster singles and doubles (CCSD) more computationally intensive than quadratic configuration interaction (QCISD)?. J. Chem. Phys. 90, 3700–3703. https://doi.org/10.1063/1.455827 (1989).

Article ADS CAS Google Scholar

Cížek, J. Advances in Chemical Physics Vol. 14 (Wiley Interscience, 1969).

Google Scholar

Peterson, K. A., Dixon, D. A. & Stoll, H. The use of explicitly correlated methods on XeF6 predicts a C3v minimum with a sterically active, free valence electron pair on Xe. J. Phys. Chem. A 116, 9777–9782. https://doi.org/10.1021/jp3084259 (2012).

Article CAS PubMed Google Scholar

Cheng, L., Gauss, J. & Stanton, J. F. Relativistic coupled-cluster calculations on XeF6: Delicate interplay between electron-correlation and basis-set effects. J. Chem. Phys. 142, 224309 (2015).

Article ADS PubMed Google Scholar

Desclaux, J. P. & Kim, Y.-K. Relativistic effects in outer shells of heavy atoms. J. Phys. B 8, 1177–1182. https://doi.org/10.1088/0022-3700/8/8/007 (1975).

Article ADS CAS Google Scholar

Wilson, A. K. Heaviest element has unusual shell structure. Physics 11, 10 (2018).

Article Google Scholar

Braïda, B. & Hiberty, P. C. The essential role of charge-shift bonding in hypervalent prototype XeF2. Nat. Chem. 5, 417–422. https://doi.org/10.1038/nchem.1619 (2013).

Article CAS PubMed Google Scholar

Avrorin, V. V., Krasikova, R. N., Nefedov, V. D. & Toropova, M. A. The chemistry of radon. Russ. Chem. Rev. 51, 12–20. https://doi.org/10.1070/rc1982v051n01abeh002787 (1982).

Article ADS Google Scholar

Han, Y.-K. & Lee, Y. S. Structures of RgF n (Rg= Xe, Rn, and element 118. n= 2, 4.) calculated by two-component spin-orbit methods. A spin-orbit induced isomer of (118) F4. J. Phys. Chem. A 103, 1104–1108 (1999).

Article CAS Google Scholar

Grandinetti, F. Noble Gas Chemistry: Structure, Bonding, and Gas-Phase Chemistry (Wiley, 2018).

Book Google Scholar

Lakshminarayanan, S., Jeyasingh, V., Murugesan, K., Selvapalam, N. & Dass, G. Molecular electrostatic potential (MEP) surface analysis of chemo sensors: An extra supporting hand for strength, selectivity & non-traditional interactions. J. Photochem. Photobiol. 6, 100022. https://doi.org/10.1016/j.jpap.2021.100022 (2021).

Article Google Scholar

Wang, H., Wang, W. & Jin, W. J. σ-hole bond vs π-hole bond: A comparison based on halogen bond. Chem. Rev. 116, 5072–5104. https://doi.org/10.1021/acs.chemrev.5b00527 (2016).

Article CAS PubMed Google Scholar

Haner, J., Matsumoto, K., Mercier, H. P. A. & Schrobilgen, G. J. Nature of the XeVI−N bonds in F6XeNCCH3 and F6Xe(NCCH3)2 and the stereochemical activity of their xenon valence electron lone pairs. Chemistry 22, 4833–4842. https://doi.org/10.1002/chem.201504904 (2016).

Article CAS PubMed Google Scholar

Frontera, A. Noble gas bonding interactions involving xenon oxides and fluorides. Molecules 25, 3419 (2020).

Article CAS PubMed PubMed Central Google Scholar

Weinstock, B., Weaver, E. E. & Knop, C. P. The xenon-fluorine system. Inorg. Chem. 5, 2189–2203. https://doi.org/10.1021/ic50046a026 (1966).

Article CAS Google Scholar

Woolf, A. A. Advances in Inorganic Chemistry and Radiochemistry Vol. 24, 1–55 (Academic Press, 1981).

Google Scholar

de la Vega, J. M. G. & Fabián, E. S. Jahn-Teller effect and dissociation from the ground state of CF4+. Chem. Phys. 151, 335–342. https://doi.org/10.1016/0301-0104(91)80019-E (1991).

Article Google Scholar

Pickard, C. J. & Needs, R. High-pressure phases of silane. Phys. Rev. Lett. 97, 045504 (2006).

Article ADS PubMed Google Scholar

Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994).

Article ADS CAS Google Scholar

Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).

Article CAS Google Scholar

Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169 (1996).

Article ADS CAS Google Scholar

Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).

Article ADS CAS PubMed Google Scholar

Gaussian09, R. et al. Gaussian 09, Revision E. 01, Gaussian, Inc, Wallingford, CT, 2004. Inc., Wallingford CT 121, 150–166 (2009).

Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98, 5648–5652. https://doi.org/10.1063/1.464913 (1993).

Article ADS CAS Google Scholar

Calhorda, M. J., Pregosin, P. S. & Veiros, L. F. Geometry optimization of a Ru(IV) allyl dicationic complex: A DFT failure?. J. Chem. Theory Comput. https://doi.org/10.1021/ct600341r (2007).

Article PubMed Google Scholar

Bartlett, R. J. & Musiał, M. Coupled-cluster theory in quantum chemistry. Rev. Mod. Phys. 79, 291–352. https://doi.org/10.1103/RevModPhys.79.291 (2007).

Article ADS CAS Google Scholar

Wadt, W. R. & Hay, P. J. Ab initio effective core potentials for molecular calculations. Potentials for main group elements Na to Bi. J. Chem. Phys. 82, 284–298. https://doi.org/10.1063/1.448800 (1985).

Article ADS CAS Google Scholar

Weigend, F. & Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 7, 3297–3305. https://doi.org/10.1039/B508541A (2005).

Article CAS PubMed Google Scholar

Stevens, W. J., Basch, H. & Krauss, M. Compact effective potentials and efficient shared-exponent basis sets for the first-and second-row atoms. J. Chem. Phys. 81, 6026–6033 (1984).

Article ADS Google Scholar

Stevens, W. J., Krauss, M., Basch, H. & Jasien, P. G. Relativistic compact effective potentials and efficient, shared-exponent basis sets for the third-, fourth-, and fifth-row atoms. Can. J. Chem. 70, 612–630 (1992).

Article CAS Google Scholar

Cundari, T. R. & Stevens, W. J. Effective core potential methods for the lanthanides. J. Chem. Phys. 98, 5555–5565 (1993).

Article ADS CAS Google Scholar

Peterson, K. A. Systematically convergent basis sets with relativistic pseudopotentials. I. Correlation consistent basis sets for the post-d group 13–15 elements. J. Chem. Phys. 119, 11099–11112. https://doi.org/10.1063/1.1622923 (2003).

Article ADS CAS Google Scholar

Xiao, L. & Wang, L. From planar to three-dimensional structural transition in gold clusters and the spin–orbit coupling effect. Chem. Phys. Lett. 392, 452–455. https://doi.org/10.1016/j.cplett.2004.05.095 (2004).

Article ADS CAS Google Scholar

Clark, T., Chandrasekhar, J., Spitznagel, G. W. & Schleyer, P. V. R. Efficient diffuse function-augmented basis sets for anion calculations. III. The 3–21+G basis set for first-row elements, Li–F. J. Comput. Chem. 4, 294–301. https://doi.org/10.1002/jcc.540040303 (1983).

Article CAS Google Scholar

Download references

This work was supported by the Institute for Korea Spent Nuclear Fuel (iKSNF) and National Research Foundation of Korea (NRF) grant funded by the Korea government (Ministry of Science and ICT, MSIT) (No.2021M2E1A1085202, 2021R1A2C2010972, 2020R1A5A1019141, and 2020R1A2C1005236). D. Y. K. also acknowledges the support of the National Natural Science Foundation of China (U1930401).

Division of Advanced Nuclear Engineering (DANE), Pohang University of Science and Technology (POSTECH), 77 Cheongam–ro, Nam–Gu, Pohang, Gyeongbuk, 790–784, Republic of Korea

Jaeeun Kang & Wooyong Um

Department of Chemistry, Pohang University of Science and Technology (POSTECH), 77 Cheongam–ro, Nam–Gu, Pohang, Gyeongbuk, 790–784, Republic of Korea

Ina Park & Ji Hoon Shim

Center for High Pressure Science & Technology Advanced Research, Shanghai, People's Republic of China

Duck Young Kim

Division of Environmental Sciences and Engineering (DESE), Pohang University of Science and Technology (POSTECH), 77 Chongam–ro, Nam–Gu, Pohang, 790–784, Republic of Korea

Wooyong Um

Nuclear Environmental Technology Institute (NETI), Pohang University of Science and Technology (POSTECH), Pohang, Gyeongbuk, 790–784, Republic of Korea

Wooyong Um

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

J.K. and I.P. conducted the computational work and wrote the manuscript. J.S. assisted with the data analyses and provided helpful comments on the manuscript. D.Y.K. and W.U. conceived and designed this study. All authors have contributed to the interpretation of results and the preparation of the manuscript, which was drafted by J.K.

Correspondence to Duck Young Kim or Wooyong Um.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Kang, J., Park, I., Shim, J.H. et al. Prediction of stable radon fluoride molecules and geometry optimization using first-principles calculations. Sci Rep 13, 2898 (2023). https://doi.org/10.1038/s41598-023-29313-5

Download citation

Received: 17 October 2022

Accepted: 02 February 2023

Published: 18 February 2023

DOI: https://doi.org/10.1038/s41598-023-29313-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.