Natural Science
Vol.6 No.4(2014), Article ID:43346,7 pages DOI:10.4236/ns.2014.64026

The helix surface of sugar ring for tetrahydrofuran*

Jin Li1, Zhenghe Zhu2#

1College of Materials and Chemical Engineering, Hainan University, Haikou, China

2Institutes of Atomic and Molecular Physics, Sichuan University, Chengdu, China; #Corresponding Author: zhuxm@scu.edu.cn

Copyright © 2014 Jin Li, Zhenghe Zhu. 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. In accordance of the Creative Commons Attribution License all Copyrights © 2014 are reserved for SCIRP and the owner of the intellectual property Jin Li, Zhenghe Zhu. All Copyright © 2014 are guarded by law and by SCIRP as a guardian.

Received 11 January 2014; revised 11 February 2014; accepted 18 February 2014

Keywords:Coupled Cluster; Electronic Structure; Excited States

ABSTRACT

The present calculations of three geometry conformation for tetrahydrofuran (THF, C4H8O) using density functional theory lead to the energy level sequence E(C1) wang#Bracket## E(Cs) < E(C2v), where, both the geometry symmetry C2v and Cs are probably a transition structure. The C1 symmetry is the equilibrium conformation, for its energy is the lowest one and it is only possible without the negative harmonic frequency. In C1 conformation, the 4 of carbon atoms C and oxygen atom O constitute a curved surface or the helix surface. The tetrahydrofuran molecule can be as a simple prototype of deoxyribose and has been reviving a quite deal of interest in radiation damage research. It is well known that the sugar ring of deoxyribose is considered as a plane in most of the references. If the equilibrium conformation of THF is C1 symmetry with the helix surface, and for the same reason, the sugar ring of deoxyribose would be also the helix surface; this probably causes the double helix of DNA. It is hardly believed that the sugar ring is a plane in deoxyribose.

1. INTRODUCTION

Tetrahydrofuran (THF) also called diethylene oxide or oxacyclopentane is a good solvent, and the study of its structure is a current focus for it can be as a moiety of nucleosides. It is well-known that a large amount of energy deposited in living cells by ionizing radiation is channeled into the production of low-energy secondary electrons. Therefore, investigation of radiation damage in living tissue upon exposure to high-energy radiation must include and apprehend these subsequent processes describing reaction between secondary electrons and cell constitutes. The dissociation of small DNA’s basic components as a function of incident energy of bombarding electrons has been shown to correlate to the measured DNA damage [1]. Being a simple prototype to investigate electron-induced decomposition of deoxyribose ring, tetrahydrofuran molecule (C4H8O) has revived a quite deal of interest in radiation damage research in recent years.

Where (a) is the tetrahydrofuran molecule (C4H8O), (b) is the deoxyribose ring.

Several papers on electron interaction with either gaseous or surface deposited THF have been recently published, which consider electron spectroscopy of resonance-enhanced vibration excitation [2], electron-stimulated desorption yields of H from thin films [3,4] and electron induced damage of solid THF film [5]. The earlier work on THF has also been summarized [5].

Since the discovery of the DNA double helix in 1953 [6], the structure of tetrahydrofuran is a kind of scientific controversy, for there are various descriptions. One of equilibrium geometries of THF is a twisted C2 group, which may be a transition structure [7,8]. However, ab initio calculation predicts that the equilibrium structure of THF is an envelope conformation of Cs [9], in which 4 of carbon atoms in a plane and oxygen atom out of this plane. Many of the works, specially, pay attention the structure of liquid tetrahydrofuran. Hydrogen/deuterium isotopic substitution neutron diffraction techniques have been used to measure the structural correlation functions of liquid tetrahydrofuran at room temperature [10]. The molecular interactions of pyrazine with THF are the weak Hydrogen bond N-C-H…O [11]. The tetrahydrofuran clathrate hydrate rapidly densifies at 77 K and ∼1.3 GPa. The densified sample is then heated to 150 K at 1.5 Gpa, and recovered at 1 bar and 77 K. Its X-ray and Raman measurements show that the densified sample is amorphous. The molecular arrangement of water in this sample appears to be similar to that in high-density amorphous ice [12].

2. DENSITY FUNCTIONAL THEORY

In DFT, electronic energy E is divided into several components,

, (1)

where ET is the electronic kinetic energy, EV is the potential energy including the electron-nucleus attraction and nucleus-nucleus repulsion, EJ is the electron-electron repulsion, and EXC is the exchange-correlation. Except nucleus-nucleus repulsion, each term can be expressed as the function of electron density

, (2)

,(3)

, (4)

where all the three terms are the functions of electron density. EX(ρ) and EC(ρ) are the exchange and correlation functional, respectively, the former is a local functional only related with the electron density, and the later is a gradient-corrected functional related with the electron density and its gradient Ñρ. In 1988, Becke [13] derived the local exchange functional as follows:

, (5)

, (6)

where ρ is a function of r, , γ is a parameter fitted with the exchange energy of noble atoms (γ = 0.0042 a.u.). Similarly, in 1992, Perdew and Wang [14]

suggested a corrected functional as

, (7)

where

Here rs is the density parameter and ζ is the correlated spin polarization. DFT is a joint calculation with the exchange and correlation functional. The B3LYP method is to combine Becke’s exchange functional and the correlated functional of Lee, Yang and Parr, with both are gradient-corrected. The local correlation functional is from Vosko, Wilk and Nusair(VWN). Then the so-called Becke’s three parameter functional [15] is

(8)

Using G1 method, Becke has obtained: c0 = 0.20, cX = 0.72 and cC = 0.81 by optimizing of exchange and correlation energy. The Eq.9 with the lower energy calculated will be used for the presented DFT-SCF calculations.

(9)

3. THE THEORETICAL DERIVATION OF MOLECULAR ELECTRONIC STATE

The molecular electronic state is the irreducible representation of molecular group, therefore, before the calculation, it is necessary to derive the possible electronic states based on the resolution, reduction and product of group representation [16-19]. Now, we derive the possible electronic states for some of possible symmetry groups, for example, tetrahydrofuran (THF) will be possible group C2v, Cs, or C1.

For simplicity, we firstly derive the electronic state of furan (C4H4O), its ground state is known as of C2v. It is supposed that we have the following reaction

All the irreducible representation for given the reactive molecules with given state are resolved into the irreducible representation of molecular group C2v of C4H4O, as follows. Oxygen O atom must be the singlet excited state, in order to result the singlet ground state of furan based on the quantization vector addition of spin. The resolution of the irreducible representation, i.e. the electronic state for 2 of C2 and 2 of H2, all are the same as

And the resolution of of O atom is

The direct product will be

There will be component for furan (C4H4O) i.e. the allowed electronic state.

For tetrahydrofuran (THF), we have the following reaction

Here, we only list these results by a table without the details of derivation.

The molecular property without indicating the symmetry and electronic state is meaningless. The electronic states derived here are proved by the following calculations.

4. CALCULATION AND COMPARISON

The present calculation is based on the density functional theory (DFT) described above and used the B3P86 method of Eq.9 with basis set 6-311G**.The calculated energy with B3P86 is usually the lowest one. The GAUSSIAN03 [20] package of program were employed to carry out all the calculations reported in this work.

The calculated geometry of C2v is shown in Figure 1, where, 4 of carbon atoms C and oxygen atom O constitute a plane, i.e. of C2v and C2 is z-axis passing through oxygen atom O and bisecting the opposite bond C-C. The electronic state is and triplet state is in a high level by 6.569 eV than. Each carbon atoms C is not to be sp2 hybridization, for each carbon atoms C has four bonds to connect four of atoms, then, it will be impossible to have plane structure. In addition, the lowest two of 33 harmonic frequencies (cm1) are negative, i.e. B1 = –222.7; A2 = –207.4 (To see the end of Table 1), this is an implicit assent of that C2v geometry would be a transition state, which is agreement with reference [7-9]. This is why to have energy level sequence C1 < Cs < C2v.

The calculated geometry of Cs is shown in Figure 2, where, 4 of carbon atoms C and oxygen atom O are not on a plane, atom O out of this plane. This is similar to the envelope Cs conformation [3]. The energy of Cs is a little lower than that of C2v. The lowest two of 33 harmonic frequencies (cm1) are also negative, i.e. A" = –222.8; A" = –207.7 (To see the end of Table 2). Therefore, Cs would be a transition state, and is not a ground state which is contrary to that of [3].

The calculated geometry of C1 is shown in Figure 3 and Table 3. The energy of this state is lower by 0.16705 ev or 1347.3 cm−1 than that of C2v with. The harmonic frequencies (cm−1) for the C1 geometry are A = 53.91, A = 266.1, A = 576.1,···, A = 3080, A = 3113.8 and A = 3122.3. Moreover, the important thing is no negative value only for this conformation, Therefore, the ground electronic state is for tetrahydrofuran (THF) in which 4 of carbon atoms C and oxygen atom O do not constitute a plane, but, like a “wave plane”, or the “helix surface” as shown in Figure 3(c), where the hydrogen atoms have been neglected.

5. CONCLUSION

The calculations with method B3P86/6-311G** show that the energy level sequence is for tetrahydrofuran (THF). Although the symmetry of C2v is relatively high, the energy of is not the lowest, or the highest one. And there are two of negative harmonic frequencies. Therefore, the C2v conformation with state is probably a transition structure.

For Cs conformation with state, there are also two of negative harmonic frequencies and it is not possible to be the equilibrium conformation of THF. This would be contrary to Reference [9].

The C1 symmetry with ground state () is the equilibrium conformation of THF, we believe, for its energy is the lowest one and it is only possible without the negative harmonic frequency. In this conformation the 4 of carbon atoms C and oxygen atom O constitute a curved surface or the helix surface shown in Figure 3(c).

In the C2v conformation, 4 of carbon atoms C and oxygen atom O all are on a plane, and for the Cs conformation, 4 of carbon atoms C are on a plane, both of them would be under the much molecular internal stress-strain, because of no sp2 hybridization for each carbon atoms C as described above.

Table 1. The results of theoretical derivation for the possible electronic state of THF.

Table 2. The geometry of C2v and Cs for THF with method (B3P86/6-311G**).

As mentioned before, tetrahydrofuran (THF) molecule (C4H8O) can be as a simple prototype of deoxyribose and has been reviving a quite deal of interest in radiation damage research. It is well known that the sugar ring of deoxyribose is considered as a plane in most of the references such as [21]. If the equilibrium conformation of THF is C1 symmetry with the helix surface, and the sugar ring of deoxyribose would be the helix surface, this

Table 3. The geometry of C1 for THF with method (B3P86/6-311G**).

(a)(b)

Figure 1. The schematic drawing of geometry for THF: (a) Top view; (b) Side view.

(a)(b)

Figure 2. The schematic drawing of geometry for THF: (a) Top view; (b) Side view.

(a)(b)(c)

Figure 3. The schematic drawing of geometry for THF: (a) Side view; (b) Top view; (c) the helix surface where all the hydrogen atoms neglected.

probably causes the double helix of DNA. It is hardly believed that the sugar ring is a plane in deoxyribose ring. Many of properties for molecular biology, for example, the sugar moiety itself in DNA being active in the initial molecular processes leading to single strand breaks, have been studied [22], however, there are a few researches on the molecular level.

ACKNOWLEDGEMENTS

This work is supported by the Natural Science Foundation of Hainan (Grant No.110001).

REFERENCES

  1. Boudaiffa, B., Cloutier, P., Hunting, D., et al. (2000) Resonant formation of DNA strand breaks by low-energy (3 to 20 eV) electrons. Science, 287, 1658-1660. http://dx.doi.org/10.1126/science.287.5458.1658
  2. Lepage, M., Letarte, S., Michaud, M., et al. (1998) Electron spectroscopy of resonance-enhanced vibrational excitations of gaseous and solid tetrahydrofuran. Journal of Chemical Physics, 109, 5980-5986. http://dx.doi.org/10.1063/1.477223
  3. Antic, D., Parenteau, L., Lepage, M., et al. (1999) Lowenergy electron damage to condensed-phase deoxyribose analogues investigated by electron stimulated desorption of H and electron energy loss spectroscopy. Journal of Chemical Physics B, 103, 6611-6619. http://dx.doi.org/10.1021/jp990686l
  4. Antic, D., Parenteau, L. and Sanche, L. (2000) Electronstimulated desorption of H from condensed phase deoxyrib ose analogues: Dissociative electron attachment versus resonance decay into dipolar dissociation. Journal of Chemical Physics B, 104, 4711-4716. http://dx.doi.org/10.1021/jp000206m
  5. Breton, S.-P., Michaud, M., Jaggle, C., et al. (2004) Damage induced by low-energy electrons in solid films of tetrahydrofuran. Journal of Chemical Physics, 121, 11240- 11249. http://dx.doi.org/10.1063/1.1814632
  6. Watson, J.D. and Crick, F.H.C. (1953) Molecular structure of nucleic acids, a structure for deoxyribose nucleic acid. Nature, 171, 737-738. http://dx.doi.org/10.1038/171737a0
  7. Mamleev, A.H., Gunderova, L.N. and Gallev, R.V. (2001) Microwave spectrum and hindered pseudorotation of tetrahydrofuran. Journal of Structural Chemistry, 42, 365- 370. http://dx.doi.org/10.1023/A:1012448618502
  8. Melnik, D.G., Gopalakrishnan, S., Miller, T.A., et al. (2003) The absorption spectroscopy of the lowest pseudorotational states of tetrahydrofuran. Journal of Chemical Physics, 118, 3589-3599. http://dx.doi.org/10.1063/1.1538241
  9. Victor, M.R. and Jose, A.S. (2005) Pseudorotation motion in tetrahydrofuran: An ab initio study. Journal of Chemical Physics, 122, 204303-1-8.
  10. Bowron, D.T., Finney, J.L. and Soper, A.K. (2006) The structure of liquid tetrahydrofuran. Journal of the American Chemical Society, 128, 5119-5126. http://dx.doi.org/10.1021/ja0583057
  11. Wang, R., Zhong, Q.L., Wu, R., et al. (2008) Molecular interactions between pyrazine and n-propanol, chloroform, or tetrahydrofuran. Spectrochimica Acta, Part A, 70, 793- 798. http://dx.doi.org/10.1016/j.saa.2007.09.014
  12. Yoshiharu, S. (2004) Evidence of pressure-induced amorphization of tetrahydrofuran clathrate hydrate. Physical Review B, 70, 172108-1-4.
  13. Becke, A.D. (1988) Density-functional exchange-energy approximation with correct asymptotic behavior. Physical Reviews, 38, 3098-3100. http://dx.doi.org/10.1103/PhysRevA.38.3098
  14. Perdew, J.P. and Wang, Y. (1992) Accurate and simple analytic representation of the electron-gas correlation energy. Physical Reviews B, 45, 13244-13249. http://dx.doi.org/10.1103/PhysRevB.45.13244
  15. Becke, A.D. (1993) Density functional thermo chemistry III, the role of exact exchange. Journal of Chemical Physics B, 98, 5648-5652.
  16. Zhu, Z.H. (1985) The general principle for determining the “microscopic” processes involving diatomic moleculeions. In the Page 14-15 of Takayanagi K. Atomic and Molecular Physics (Japanese). The Institute of Space and Astronautically Science, Tokyo.
  17. Zhu, Z.H. (2007) Atomic and molecular reaction statics. Science in China Series G, 50, 581-590. http://dx.doi.org/10.1007/s11433-007-0054-6
  18. Zhu, Z.H. (2007) Atomic and molecular reaction statics. Science Press, Beijing.
  19. Klumpp, G.W. (1982) Reactivity in organic chemistry. A Wiley-Interscience Publication.
  20. Frisch, M.J., Trucks, G.W., Schlegel, H.B., et al. (2003) Gaussian03. Gaussian, Inc., Pittsburgh.
  21. Cheng, C.X. (1998) Biophysical chemistry. Science and Technology Press, Beijing, 387.
  22. Philipp, S., Sylwia, P., Fabio, Z., et al. (2006) Dissociative electron attachment to furan, tetrahydrofuran, and fructose. Journal of Chemical Physics, 125, 044304-1-6.

NOTES

*PACC numbers: 3120C; 3130; 3150.