The helix surface of sugar ring for tetrahydrofuran *

The present calculations of three geometry conformation for tetrahydrofuran (THF, C4H8O) using density functional theory lead to the energy level sequence E(C1) < 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.


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 chan-neled 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 (C 4 H 8 O) has revived a quite deal of interest in radiation damage research in recent years.
Where (a) is the tetrahydrofuran molecule (C 4 H 8 O), (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 C 2 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 C s [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].

DENSITY FUNCTIONAL THEORY
In DFT, electronic energy E is divided into several components, where E T is the electronic kinetic energy, E V is the potential energy including the electron-nucleus attraction and nucleus-nucleus repulsion, E J is the electron-electron repulsion, and E XC is the exchange-correlation. Except nucleus-nucleus repulsion, each term can be expressed as the function of electron density ρ ( )( ) ( ) where all the three terms are the functions of electron density.E X (ρ) and E C (ρ) 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: ( ) where ρ is a function of r, 4 3 x ρ ρ − = ▽ , γ is a parame- ter fitted with the exchange energy of noble atoms (γ = 0.0042 a.u.).Similarly, in 1992, Perdew and Wang [14] suggested a corrected functional as ( ) Here r s 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 Using G1 method, Becke has obtained: c 0 = 0.20, c X = 0.72 and c C = 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.

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][17][18][19].Now, we derive the possible electronic states for some of possible symmetry groups, for example, tetrahydrofuran (THF) will be possible group C 2v , C s , or C 1 .
For simplicity, we firstly derive the electronic state of furan (C 4 H 4 O), its ground state is known as Χ Α  of furan based on the quantization vector addition of spin.The resolution of the irreducible representation, i.e. the electronic state for 2 of C 2 and 2 of H 2 , all are the same as The direct product will be 2 There will be component 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.

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 C 2v is shown in Figure 1, where, 4 of carbon atoms C and oxygen atom O constitute a plane, i.e.

( )
of C 2v and C 2 is z-axis passing through oxygen atom O and bisecting the opposite bond C-C.The electronic state is 1   1 Α and triplet state 3   1 Β is in a high level by 6.569 eV than 1 1 Α .Each carbon atoms C is not to be sp 2 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 (cm -1 ) are negative, i.e.B 1 = -222.7;A 2 = -207.4(To see the end of Table 1), this is an implicit assent of that C 2v geometry would be a transition state, which is agreement with reference [7][8][9].This is why to have energy level sequence C 1 < C s < C 2v .
The calculated geometry of C s 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 C s conformation [3].The energy of C s is a little lower than that of C 2v .The lowest two of 33 harmonic frequencies (cm -1 ) are also negative, i.e.A" = -222.8;A" = -207.7 (To see the end of Table 2).Therefore, C s would be a transition state, and is not a ground state which is contrary to that of [3].
The calculated geometry of C 1 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 C 2v with 1   1   Α .The harmonic frequencies (cm −1 ) for the C 1 geometry are A = 53.91,A = 266.1,A = 576.1,•••,A = 3080, A = 3113.8and 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.

CONCLUSION
The calculations with method B3P86/6-311G** show that the energy level sequence is ( ) for tetrahydrofuran (THF).Although the symmetry of C 2v is relatively high, the energy of ( ) C , Α is not the lowest, or the highest one.And there are two of negative harmonic frequencies.Therefore, the C 2v conformation with state 1 1 Α is probably a transition structure.
For C s conformation with state 1 ′ Α , 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 C 1 symmetry with ground state ( 1 Χ Α  ) 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 C 2v conformation, 4 of carbon atoms C and oxygen atom O all are on a plane, and for the C s 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 sp 2 hybridization for each carbon atoms C as described above.
Resolution and direct product Resolution and direct product Resolution and direct product Α ⊗ Α ⊗ Α = Α for ( ) As mentioned before, tetrahydrofuran (THF) molecule (C 4 H 8 O) 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 C 1 symmetry with the helix surface, and the sugar ring of deoxyribose would be the helix surface, this Table 3.The geometry of C 1 for THF with method (B3P86/6-311G**).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.

Figure 3 . 1 C ,Χ Α 1 
Figure 3.The schematic drawing of geometry for THF ( ) 1 C ,Χ Α 1  : (a) Side view; (b) Top view; (c) the helix surface All the irreducible representation for given the reactive molecules with given state are resolved into the irreducible representation of molecular group C 2v of C 4 H 4 O, as follows.Oxygen O atom must be the singlet excited state

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