An ab initio theoretical study of electronic structure and properties of 2'-deoxyguanosine in gas phase and aqueous media

Molecular geometries of two structural forms of 2'-deoxyguanosine (keto-N9R and keto-N7R, R = the sugar moiety) considering both the C2'-endo and C3'-endo conformations of the sugar ring and those of the complexes of these species with two water molecules each were optimized employing the ab initio RHF procedure. A mixed basis set consisting of the 6-311+G* basis set for the nitrogen atom of the amino group and the 4-31G basis set for all the other atoms was used. The RHF calculations were followed by correlation correction of the total energy at the MP2 level. Both the structural forms of 2'-deoxyguanosine were solvated using the polarized continuum model (PCM) of the self-consistent reaction field (SCRF) theory and the corresponding RHF optimized geometries at the RHF and MP2 levels. Geometry optimization was also performed in aqueous media using the Onsager model at the RHF level using the above-mentioned mixed basis set, and subsequently, using the reoptimized geometries, single-point MP2 calculations were performed. It is found that both the keto-N9R and keto-N7R forms of 2'-deoxyguanosine as well as their complexes with two water molecules each would occur, particularly at the water-air interface. Though the normal Watson-Crick-type base pairing would not be possible with the keto-N7R form of 2'-deoxyguanosine(G*), two other (G*-C and G*-T) base pairing schemes may occur with this form of the nucleoside, which may cause mutation. The present calculated geometry of the keto-N9R form of the anti-conformation of 2'-deoxyguanosine including the dihedral angle chi(CN) agree satisfactorily with the available crystallographic results. The present results also agree satisfactorily with those obtained by other authors earlier for the keto-N9R form of 2'-deoxyguanosine using B3LYP and MP2 methods employing the 6-31G* basis set. Using transition state calculations, it is shown that tautomerism of guanine and other similar molecules where the tautomers would coexist would be facilitated by the occurrence of the H(+) and OH(-) fragments of water molecules. Further, this coexistence of the two tautomers appears to make the C8 carbon atom located between the N7 and N9 nitrogen atoms susceptible to attack by the OH(-) group. Thus, an explanation is obtained for the efficient formation of the reaction product 8-hydroxy-2'-deoxyguanosine, which serves as a biomarker for oxidative damage to DNA in biological systems.     

Electronic spectra of adenine and 2-aminopurine: an ab initio study of energy level diagrams of different tautomers in gas phase and aqueous solution

Ground and lowest two singlet excited state geometries of four tautomeric forms (N9H, N7H, N3H and N1H) of each of adenine and 2-aminopurine (2AP) were optimized using an ab initio approach employing a mixed basis set (6-311 + G* on the nitrogen atom of the amino group and 4-31G basis set on the other atoms). Excited states were generated employing configuration interaction involving single electron excitations (CIS). Subsequently, the different species were solvated in water employing the self-consistent reaction field (SCRF) approach along with the corresponding gas phase optimized geometries. Thus the observed absorption and fluorescence spectra of adenine and 2AP have been explained successfully. It is concluded that both the N9H and N7H forms of 2AP would contribute to absorption and fluorescence spectra. Further, the fluorescence of 2AP would be absorbed by its cation in which both the N9 and N7 atoms are protonated, the fluorescence of which can have an anti-Stokes component. Among the different tautomers of adenine, the N9H form would be present dominantly in the ground state in aqueous solutions but the N7H form would be produced by energy transfer and subsequent fluorescence. The N3H form of adenine appears to be responsible for the observed absorption near 300 nm by its solutions intermittently exposed to ultraviolet radiation. The rings of the different species related to 2AP and adenine remain almost planar in the pi-pi* and n-pi* singlet excited states as in the ground state. The pyramidal character of the amino group is usually less in the pi-pi* excited states than that in the corresponding ground or n-pi* excited states. Molecular electrostatic potential (MEP) maps of the molecules provide useful clues regarding phototautomerism.     

An ab initio study of excited states of guanine in the gas phase and aqueous media: Electronic transitions and mechanism of spectral oscillations

Ground state geometries of the four tautomeric forms keto‐N9H, keto‐N7H, enol‐N9H, and enol‐N7H of guanine were optimized in the gas phase at the RHF level using a mixed basis set consisting of the 4‐31G basis set for all the atoms except the nitrogen atom of the amino group for which the 6‐311+G* basis set was used. These calculations were also extended to hydrogen‐bonded complexes of three water molecules with each of the keto‐N9H (G9‐3W) and keto‐N7H (G7‐3W) forms of guanine. Relative stabilities of the four above‐mentioned tautomers of guanine as well as those of G9‐3W and G7‐3W complexes in the ground state in the gas phase were studied employing the MP2 correlation correction. In aqueous solution, relative stabilities of these systems were studied using the MP2 correlation correction and polarized continuum model (PCM) or the isodensity surface polarized continuum model (IPCM) of the self‐consistent reaction field (SCRF) theory. Geometry optimization in the gas phase at the RHF level using the 6‐31+G* basis set for all atoms and the solvation calculations in water at the MP2 level using the same basis set were also carried out for the nonplanar keto‐N9H and keto‐N7H forms of guanine. Thus, it is shown that among the different tautomers of guanine, the keto‐N7H form is most stable in the gas phase, while the keto‐N9H form is most stable in aqueous solution. It appears that both the keto‐N9H and keto‐N7H forms of guanine would be present in the ground state, particularly near the aqueous solution–air interface. Vertical excitation and excited state geometry optimization calculations were performed using configuration interaction involving single electron excitation (CIS). It is found that the absorption spectrum of guanine would arise mainly due to its keto‐N9H form but the keto‐N7H form of the same would also make some contribution to it. The enol‐N9H and enol‐N7H forms of the molecule are not expected to occur in appreciable abundance in the gas phase or aqueous media. The normal fluorescence spectrum of guanine in aqueous solution with a peak near 332 nm seems to originate from the lowest singlet excited state of the keto‐N7H form of the molecule while the fluorescence of oxygen‐rich aqueous solutions of guanine with a peak near 450 nm appears to originate from the lowest singlet excited state of the keto‐N9H form of the molecule. The origin of the slow damped spectral oscillation observed in the absorption spectrum of guanine has been explained. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 826–846, 2000     

The half projected Hartree-Fock function for determining singlet excited states. Application to cyclobutanone and 3-cyclopenten-1-one

The Half-Projected-Hartree-Fock procedure (HPHF) for determining singlet ground states is briefly described and extended to the direct determination of singlet excited states. The procedure is applied, using a [7s,3p/2s,1p] basis set, to determine the optimal geometry of two relatively large molecules, to which large CI calculations cannot be easily applied. These two molecules are cyclobutanone and 3-cyclopenten-1-one in their lowest singlet (n → π) excited state. Both molecules are found to exhibit in their excited state a pyramidal structure with the carbonyl oxygen atom pointing outward from the molecular plane. RHF calculations for the singlet ground state were also performed for comparison. The theoretical geometrical parameters compare well with the experimental data.     

Theoretical study of some 1,3-substituted (1,3-diketonato)borondifluorides. Potential energy curve relevant to the barrier crossing reaction

Ab initio Hartree-Fock calculations utilising STO-3G, 3-21G^* and 6-31G^* basis sets have been performed on three neutral and highly polar molecules, (diformylmethine)borondifluoride, (acetylacetonato)borondifluoride and (dibenzylmethine) borondifluoride. The calculated and experimental structures are well correlated when using the HF/3-21G^* basis set, except for the structure parameters involving the boron atom. The HF/6-31G^* basis set does not improve the accuracy in structure calculations. The conformational analysis is in agreement with the experimentally observed C_2v symmetrical structures, where the boron atom is tetrahedrally coordinated. The calculations support a one-dimensional ground state barrier crossing reaction for (dibenzylmethine)borondifluoride, where the phenyl torsion is the most likely reaction coordinate. Both HF/6-31G^* calculations and the second-order Møller-Plesset correction with the 3-21G^* basis set suggest an activation energy of the ground state reaction of about 30 kJ mol⁻¹. The ground state barrier crossing reaction kinetics is evaluated by the Kramers theory. The calculated ground state parameters relevant to the barrier crossing reaction are compared with the experimentally observed excited state values.     

Hydration-dependent structural deformation of guanine in the electronic singlet excited state

Theoretical study was performed to investigate how the degree of hydration affects the structures and properties of the canonical form (keto-N9H) of guanine in the ground and lowest singlet pipi* excited state. This work is the continuation of our earlier work where we have studied the hydration of guanine in the first solvation shell with one, three, five, and six water molecules. In the present investigation, we have considered 7-13 water molecules in hydrating guanine. Ground-state geometries were optimized at the Hartree-Fock level, whereas the configuration interaction-singles (CIS) method was used for the excited-state geometry optimization. The 6-311G(d,p) basis set was used in all calculations. The harmonic vibrational frequency analysis was used to determine the nature of the optimized ground- and excited-state potential energy surfaces; all geometries were found to be minima at the respective potential surfaces. It was found that the degree of hydration has a significant influence on the excited-state structural nonplanarity of guanine. It is expected that excited-state dynamics of guanine will depend on the degree of hydration. Ground- and excited-state geometries of selected hydrated species were also optimized in the bulk water solution using the polarizable continuum model (PCM). It was found that bulk water solution generally does not have significant influence on the structure of the hydrated species. Effects of hydration on different stretching vibrations in the ground and excited states are also discussed.     

1. 5-环辛二烯-3, 7-二炔(C~8H~4)结构和光谱性质的理论研究

采用量子化学从头算方法研究了1,5-环辛二烯-3,7-二炔(C~8H~4)的结构和光谱性质,根据等键反应分析和自然键轨道方法研究了它的稳定性、成键情况和共轭性。结果表明1,5-环辛二烯-3,7-二炔(C~8H~4)分子为平面刚性结构,可能稳定存在。分子中C≡C键与C=C键存在一定程度的共轭,可能具有芳香性。     

去氢抗坏血酸分子振动光谱的理论研究

采用RHF, MP2, DFT(B3LYP)方法, 以6-311++G**为基组研究了去氢抗坏血酸分子(DHA)的平衡几何构型和振动光谱. 计算结果表明, 采用RHF, B3LYP以及MP2 方法优化得到的几何结构以及频率值是一致的. 采用B3LYP/6-311++G**计算了DHA分子平衡构型下的谐振动力场﹑振动频率和振动强度. 使用Wilson的GF矩阵方法对DHA分子进行了简正坐标分析, 依据所得的势能分布对DHA分子的振动基频进行了合理的理论归属.     

Comparison of ground and excited state polarizabilities of thiophene,cyclopentadiene, and fulvene oligomers and their cyano substituted derivatives—Ab initio study

Ab Initio study of the ground and excited state polarizabilities of thiophene, fulvene, and cyclopentadiene based conducting oligomers and their cyano derivatives have been performed using the restricted Hartree–Fock (RHF) and the configuration interaction singles (CIS) approaches, respectively, with 3‐21G* basis set. For comparison purposes, for some small oligomers (monomers and dimers), higher basis sets (6‐31G*, 6‐31+G*, aug‐cc‐pVTZ) were also employed in the computations of polarizabilities. The trends in polarizability as a function of oligomer length were investigated. For all systems, the RHF polarizability increases as n^1.2–1.3 as n gets larger and the CIS polarizability increases as n^1.4–1.6 for n less than seven or eight rings and then increases approximately linearly with n for larger n. For the thiophene based systems the dependence of the polarizability on bond length alternation (BLA) along the backbone of the oligomers was also investigated using the RHF, density functional (DFT), and CIS theories (with 3‐21G* basis set). For thiophene dimer, we also performed RHF/aug‐cc‐pVTZ calculations of polarizabilities versus BLA. We found that the polarizability is largest when BLA is near zero (for both ground and excited states), which correlates with the lowest excitation energy. Comparison with experimental results has been made where possible. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1983–1995, 2007     

Excited state proton transfer in guanine in the gas phase and in water solution: a theoretical study

Theoretical investigations were performed to study the phenomena of ground and electronic excited state proton transfer in the isolated and monohydrated forms of guanine. Ground and transition state geometries were optimized at both the B3LYP/6-311++G(d,p) and HF/6-311G(d,p) levels. The geometries of tautomers including those of transition states corresponding to the proton transfer from the keto to the enol form of guanine were also optimized in the lowest singlet pipi* excited state using the configuration interaction singles (CIS) method and the 6-311G(d,p) basis set. The time-dependent density function theory method augmented with the B3LYP functional (TD-B3LYP) and the 6-311++G(d,p) basis set was used to compute vertical transition energies using the B3LYP/6-311++G(d,p) geometries. The TD-B3LYP/6-311++G(d,p) calculations were also performed using the CIS/6-311G(d,p) geometries to predict the adiabatic transition energies of different tautomers and the excited state proton transfer barrier heights of guanine tautomerization. The effect of the bulk aqueous environment was considered using the polarizable continuum model (PCM). The harmonic vibrational frequency calculations were performed to ascertain the nature of potential energy surfaces. The excited state geometries including that of transition states were found to be largely nonplanar. The nonplanar fragment was mostly localized in the six-membered ring. Geometries of the hydrated transition states in the ground and lowest singlet pipi* excited states were found to be zwitterionic in which the water molecule is in the form of hydronium cation (H3O(+)) and guanine is in the anionic form, except for the N9H form in the excited state where water molecule is in the hydroxyl anionic form (OH(-)) and the guanine is in the cationic form. It was found that proton transfer is characterized by a high barrier height both in the gas phase and in the bulk water solution. The explicit inclusion of a water molecule in the proton transfer reaction path reduces the barrier height drastically. The excited state barrier height was generally found to be increased as compared to that in the ground state. On the basis of the current theoretical calculation it appears that the singlet electronic excitation of guanine may not facilitate the excited state proton transfer corresponding to the tautomerization of the keto to the enol form.     

尿嘧啶和5-氯尿嘧啶1S0→1S2跃迁动态结构的共振拉曼光谱

采用含时量子波包理论的简单模型对5-氯尿嘧啶和尿嘧啶的共振拉曼光谱开展了强度分析拟合, 获得了1(π, π*)激发态的几何结构变化动态特征. 结果表明, 尿嘧啶1S0→1S2跃迁的动态结构特征因5-位氯原子取代而改变. 5-氯尿嘧啶的动态结构特征主要沿C5=C6伸缩振动+C6H12 弯曲振动和N3H9/N1H7弯曲振动+N1C6伸缩振动反应坐标展开, 而尿嘧啶的动态结构特征主要沿嘧啶环的伸缩振动+C5H11/C6H12/N1H7弯曲振动和C4=O10伸缩振动反应坐标展开. π和π*轨道中氯原子的pz电子参与嘧啶环的p-π共轭作用导致了在1(π, π*)激发态上5-氯尿嘧啶的振动重组能更多地配分给嘧啶环的弯曲振动模式和C5=C6伸缩振动模式. 尿嘧啶在甲醇中的激发态动态结构特征与在水中的基本一致, 但波包沿C5H11/C6H12/N1H7弯曲振动+N1C6伸缩振动(υ12)和环呼吸振动(υ17)反应坐标的运动明显增强.     

A gas phase ab initio excited state geometry optimization study of thymine, cytosine and uracil

Molecular geometries of the nucleic acid bases thymine, cytosine and uracil in the ground and the lowest two singlet excited states were optimized using the ab initio approach employing the 4-31G basis set for all the atoms except the amino group of cytosine for which the 6-311+G^* basis set was used. The excited state calculations were performed employing configuration interaction involving singly excited configurations (CIS). Vibrational frequencies were computed in order to examine the nature of the stationary points on the potential energy surfaces obtained by geometry optimization. While the ground state geometries of uracil and thymine (except the methyl group hydrogens) are planar, the corresponding excited state geometries were found to be significantly nonplanar. In the case of cytosine, the amino group is pyramidal and the rest of the molecule is only slightly nonplanar in the ground state, but the excited state geometries are appreciably nonplanar. In particular, consequent to the S₂(n–π^*) excitation of cytosine, the amino group plane is strongly rotated. While thymine is stable in the S₂(π–π^*) excited state, uracil appears to be dissociative in the corresponding excited state.     

Dibenzo-p-dioxin. An ab initio CASSCF/CASPT2 study of the pi-pi* and n-pi* valence excited states

The pi-pi* and n-pi* valence excited states of dibenzo-p-dioxin (DD) were studied via the complete active space SCF and multiconfigurational second-order perturbation theory employing the cc-pVDZ basis set and the full pi-electron active spaces of 16 electrons in 14 active orbitals. The geometry and harmonic vibrational wavenumbers of the ground state correlate well with the experimental and other theoretical data. In particular, significant improvements over previously reported theoretical results are observed for the excitation energies. All of the pi-pi* excited states exhibit planar D(2h)minima. Thus no evidence was found for a C(2v) butterfly-like relaxation, although the wavenumbers of the b(3u) butterfly flapping mode proved exceedingly low in both the ground S(0)((1)A(g)) and the lowest dipole allowed excited S(1)((1)B(2u)) state. The calculations of oscillator strengths established the 2(1)B(2u) <-- 1(1)A(g) and 2(1)B(1u) <-- 1(1)A(g) transitions as by far the most intense, whereas the only allowed of the n-pi* transitions ((1)B(3u)) should possess only a modest intensity. Studies into dependence of the oscillator strengths on the extent of the butterfly-like folding showed that the electronic spectrum is more consistent with a folded equilibrium geometry assumed by DD in solution.     

Excited state properties of 7-hydroxy-4-methylcoumarin in the gas phase and in solution. A theoretical study

TDDFT/B3LYP and RI-CC2 calculations with different basis sets have been performed for vertical and adiabatic excitations and emission properties of the lowest singlet states for the neutral (enol and keto), protonated and deprotonated forms of 7-hydroxy-4-methylcoumarin (7H4MC) in the gas phase and in solution. The effect of 7H4MC-solvent (water) interactions on the lowest excited and fluorescence states were computed using the Polarizable Continuum Method (PCM), 7H4MC-water clusters and a combination of both approaches. The calculations revealed that in aqueous solution the pi pi* energy is the lowest one for excitation and fluorescence transitions of all forms of 7H4MC studied. The calculated excitation and fluorescence energies in aqueous solution are in good agreement with experiment. It was found that, depending on the polarity of the medium, the solvent shifts vary, leading to a change in the character of the lowest excitation and fluorescence transition. The dipole-moment and electron-density changes of the excited states relative to the ground state correlate with the solvation effect on the singlet excited states and on transition energies, respectively. The calculations show that, in contrast to the ground state, the keto form has a lower energy in the pi pi* state as compared to enol, demonstrating from this point of view the energetic possibility of proton transfer from the enol to the keto form in the excited state.     

Effect of hydration on the lowest singlet PiPi* excited-state geometry of guanine: a theoretical study

An ab-initio computational study was performed to investigate the effect of explicit hydration on the ground and lowest singlet PiPi* excited-state geometry and on the selected stretching vibrational frequencies corresponding to the different NH sites of the guanine acting as hydrogen-bond donors. The studied systems consisted of guanine interacting with one, three, five, six, and seven water molecules. Ground-state geometries were optimized at the HF level, while excited-state geometries were optimized at the CIS level. The 6-311G(d,p) basis set was used in all calculations. The nature of potential energy surfaces was ascertained via the harmonic vibrational frequency analysis; all structures were found minima at the respective potential energy surfaces. The changes in the geometry and the stretching vibrational frequencies of hydrogen-bond-donating sites of the guanine in the ground and excited state consequent to the hydration are discussed. It was found that the first solvation shell of the guanine can accommodate up to six water molecules. The addition of the another water molecule distorts the hydrogen-bonding network by displacing other neighboring water molecules away from the guanine plane.     

Excited-state proton transfer behavior of 7-hydroxy-4-methylcoumarin in AOT reverse micelles

The excited-state proton transfer and phototautomerization of 7-hydroxy-4-methylcoumarin (7H4MC) dye has been studied in the confined water pools of AOT reverse micelles using steady-state and time-resolved fluorescence measurements. In the "dry" reverse micelles ([water]/[AOT], w(0) = 0), only the neutral form of the dye is present both in the ground and the excited states. At higher w(0) values, three prototropic forms, namely, neutral, anionic, and tautomeric, can be identified in the excited state, although only the neutral form of the dye is present in the ground state. From steady-state fluorescence results and time-resolved area-normalized emission spectra (TRANES), it is indicated that the anionic and tautomeric forms of the dye are the excited-state reaction products and that they arise apparently independently from the excited neutral form of the dye. In bulk water, however, there is no evidence of the tautomeric species and only the anionic form is observed in the excited state. The fluorescence quenching results of the three forms of 7H4MC by the different quenchers, potassium iodide, aniline, and N, N-dimethylaniline, suggest that the distribution of 7H4MC molecules in the reverse micelles is not diverse but that the different prototropic forms arise from the same population of the excited dye in the interfacial region.     

Electronic transitions of guanine tautomers, their stacked dimers, trimers and sodium complexes

Planar and nonplanar geometries of the keto-N9H and keto-N7H tautomers of the guanine base of DNA as well as the hydrogen bonded complexes of these species with three water molecules each were optimized using the density functional theory at the B3LYP/6-31G** level. Geometries of the isolated bases were also optimized using the ab initio approach at the MP2/6-31G** level. The isolated keto-N9H and keto-N7H tautomers as well as their hydrogen bonded complexes with three water molecules each were solvated in bulk water employing the polarized continuum model (PCM) of the self-consistent reaction field theory (SCRF). Stacked dimers and trimers of both the tautomers of guanine were generated by placing the planar forms of the species at a fixed distance of 3.5 A from the neighboring one and rotating one molecule with respect to the other by 110 degrees for the keto-N9H form and 90 degrees for the keto-N7H form which corresponded to total energy minima at the B3LYP/6-31G** level. Geometry optimization for the cation of the monomer of guanine was performed at the same level of theory, and its solvation in bulk water was treated using the PCM model of the SCRF theory. The geometries of complexes of the two tautomers of guanine with a Na+ ion each were optimized at the B3LYP/6-31G** level, and the Na+ ion is predicted to bind with the keto-N9H tautomer preferentially. While the complex of the keto-N7H form of guanine with three water molecules in gas phase is slightly more stable than the corresponding complex of the keto-N9H form of guanine, the reverse is true in bulk water. Stacking interactions enhance the relative stability of the keto-N9H tautomer over that of the keto-N7H tautomer, suggesting that in bulk solutions, the former would be dominant. Electronic spectra of the isolated tautomers of guanine, those of their complexes with three water molecules each, the (keto-n9h and keto-n7h) cation of guanine, the complexes of the tautomers with a Na+ ion each, the stacked dimers and trimers of the two tautomers were calculated using configuration interaction involving single electron excitations (CIS). The relative absorption intensities of the two tautomers of guanine near 275 and 248 nm in the monomer, dimer, and trimer are predicated to be in the opposite order. Thus the absorption intensity oscillation observed using a guanine aqueous solution can be explained in terms of oscillation of relative populations of the two tautomers of the molecule. The 248 nm absorption peak would be appreciably red-shifted on formation of the cation of guanine. Binding of the Na+ ion with the two tautomers of guanine reduces intensities of their transitions appreciably and also it causes large red-shifts in the same.     

Solid-state 25Mg NMR spectroscopic and computational studies of organic compounds. square-pyramidal magnesium(II) ions in aqua(magnesium) phthalocyanine and chlorophyll a

We report a solid-state (25)Mg NMR spectroscopic study of two magnesium-containing organic compounds: monopyridinated aqua(magnesium) phthalocyanine (MgPc.H(2)O.Py) and chlorophyll a (Chla). Each of these compounds contains a Mg(II) ion coordinating to four nitrogen atoms and a water molecule in a square-pyramidal geometry. Solid-state (25)Mg NMR spectra for MgPc.H(2)O.Py were obtained at 11.7 T (500 MHz for (1)H) for a (25)Mg-enriched sample (99.1% (25)Mg atom) using both Hahn-echo and quadrupole Carr-Purcell Meiboom-Gill (QCPMG) pulse sequences. Solid-state (25)Mg NMR spectra for Chla were recorded at (25)Mg natural abundance (10.1%) at 19.6 T (830 MHz for (1)H). The (25)Mg quadrupole parameters were determined from spectral analyses: MgPc.H(2)O.Py, C(Q) = 13.0 +/- 0.1 MHz and eta(Q) = 0.00 +/- 0.05; Chla, C(Q) = 12.9 +/- 0.1 MHz and eta(Q) = 1.00 +/- 0.05. This work represents the first time that Mg(II) ions in a square-pyramidal geometry have been characterized by solid-state (25)Mg NMR spectroscopy. Extensive quantum mechanical calculations for electric-field-gradient (EFG) and chemical shielding tensors were performed at restricted Hartee-Fock (RHF), density functional theory (DFT), and second-order M?ller-Plesset perturbation theory (MP2) levels for both compounds. Computed (25)Mg nuclear quadrupole coupling constants at the RHF and MP2 levels show a reasonable basis-set convergence at the cc-pV5Z basis set (within 7% of the experimental value); however, B3LYP results display a drastic divergence beyond the cc-pVTZ basis set. A new crystal structure for MgPc.H(2)O.Py is also reported.     

磷的d轨道在H~3PO分子中的作用

用分子片轨道在分子环境中发生极化的概念研究d轨道在H~3PO分子中的作用。H~3PO分子被分为两个分子片---H~3P和O.在RHF/6-31G^*水平上计算出分子环境中的极化了的分子片轨道(FOM)。再剔除d函数为主的FOM,用剩余的FOM为基进行构型优化,得到与RHF/6-31G^*相近的结果。这一结果说明磷原子的d函数在H~3PO分子中仅仅起一个极化函数的作用,而不是起价轨道作用。     

二氧化钚分子的多体展式势能函数

从导出基态PuO2分子的电子状态X5Σ g正确地判断其离解极限出发,采用MP2方法,应用相对论有效原子实模型(RECP)优化出PuO2(X5Σ g)分子稳定构型为线性OPuO(D∞h),其平衡核间距Re=0.18004nm.同时也计算出振动频率,并优化出存在亚稳态的Pu-O-O(C∞v)构型.使用多体项展式理论方法,导出了基态PuO2分子的分析势能函数.该势能表面准确地再现了O-Pu-O(D∞h)平衡结构和亚稳态的Pu-O-O(C∞v)构型.然后根据势能函数等值图讨论了O(3Pg) PuO反应的势能面静态特征.