腺嘌呤和质子化腺嘌呤的结构和振动光谱(完整版)实用资料.doc

《腺嘌呤和质子化腺嘌呤的结构和振动光谱(完整版)实用资料.doc》由会员分享，可在线阅读，更多相关《腺嘌呤和质子化腺嘌呤的结构和振动光谱(完整版)实用资料.doc（229页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、腺嘌呤和质子化腺嘌呤的结构和振动光谱（完整版）实用资料(可以直接使用，可编辑 完整版实用资料，欢迎下载）Article物理化学学报(Wuli Huaxue Xuebao Acta Phys.-Chim.Sin .,2021,25(8:1605-1610August Received:March 5,2021;Revised:April 29,2021;Published on Web:June 11,2021.*Corresponding author.Email:dywuxmu.edu ;Tel:+86-592-2189023.The project was supported by the

2、 National Natural Science Foundation of China (20573087,National Key Basic Research Program of China (973(2007CB815303,2021CB930703and Program for New Century Excellent Talent (NCETin Universities of Fujian Province,China.国家自然科学基金(20573087、国家重点基础研究发展计划(973项目(2007CB815303、2021CB930703和福建省高等学校新世纪人才支持计

3、划项目(NCET资助鬁Editorial office of Acta Physico -Chimica Sinica腺嘌呤和质子化腺嘌呤的结构和振动光谱梁晓静崔丽吴德印*田中群(固体表面物理化学国家重点实验室,厦门大学化学化工学院化学系,福建厦门361005摘要:用密度泛函方法B3LYP/aug -cc -pVTZ 分析了腺嘌呤和质子化腺嘌呤的低能稳定异构体的结构和振动光谱.结果发现,对于中性腺嘌呤分子,腺嘌呤的异构体N9H 比N7H 的能量低32.76kJmol -1(在极化连续模型下为6.28kJ mol -1.基于标度量子力场方法所得到的势能分布,对异构体N9H 的部分振动基频重新进行

4、了归属.在极化连续模型下,质子化腺嘌呤分子有5种低能稳定构型,其中N1位质子化的9-位氢腺嘌呤最为稳定.基于振动模式分析,对这种最稳定构型的振动基频进行了归属,并对腺嘌呤在pH=1的高氯酸溶液中的实验拉曼光谱进行了指认.关键词:腺嘌呤;密度泛函理论;振动光谱;质子化腺嘌呤中图分类号:O641Structures and Vibrational Spectra of Adenine and Protonated AdenineLIANG Xiao -JingCUI LiWU De -Yin *TIAN Zhong -Qun(State Key Laboratory of Physical Che

5、mistry of Solid Surfaces,Dopartment of Chemistry,College of Chemistry and ChemicalEngineering,Xiamen University,Xiamen 361005,Fujian Province,P.R.China Abstract :The structures and vibrational spectra of neutral and protonated adenine molecules were calculated at the B3LYP/aug -cc -pVTZ level.For ne

6、utral adenine,the N9H adenine configuration is more stable in energy of about 32.76kJ mol -1(6.28kJ mol -1by using the polarization continuous model (PCMthan N7H adenine.Based on the potential energy distribution (PEDcalculated using the scaled quantum mechanical field (SQMFprocedure,we corrected th

7、e assignments of some N9H adenine fundamental vibrations.There are five stable configurations for protonated adenine and the isomer of adenine that is protonated at the N1position is the most stable.Based on vibration analysis,we assigned the fundamental vibrations of this configuration and analyzed

8、 the Raman spectra of adenine in the HClO 4(pH=1solution.Key Words :Adenine;Density functional theory;Vibrational spectrum;Protonated adenine由于腺嘌呤及其衍生物是核酸和储能分子ATP 的重要组成部分,同时很多重要的生物过程诸如碱基错配、核苷酸的自聚作用以及由DNA 中互补碱基对的氢键质子隧道所引起的自发点突变等都与核酸碱基质子化有关,尤其涉及腺嘌呤质子化的参与1.虽然已经有较多有关该分子的振动光谱的实验和理论研究2-8,但至今对质子化腺嘌呤分子的振动光谱

9、的归属仅限于理论上HF/STO -3G 低水平指认1以及实验上的经验指认9,10.因此,进一步研究质子化腺嘌呤的分子结构以及质子化对其振动光谱的影响在振动光谱应用于生命科学上有着重要意义1.1605Acta Phys.-Chim.Sin.,2021Vol.25至今有关腺嘌呤分子的振动光谱实验数据有:腺嘌呤在低温氩基质下的红外光谱3,腺嘌呤晶体的红外1和拉曼光谱11,质子化腺嘌呤在固相9和气相10的红外光谱.由于该分子在低温氩基质下红外光谱的谱峰信噪比相对较好,而且分子间相互作用较小,因此本文同时参照了该实验的研究工作3,10.本文采用B3LYP 12-16/aug -cc -pVTZ 17方法

10、计算腺嘌呤分子和质子化腺嘌呤的可能低能构型,并基于标度量子力场方法以及结合理论计算和实验测量谱频率和强度对最稳定结构的振动基频进行了指认.通过比较前人对该分子基频归属的分析2-4,我们对腺嘌呤分子的个别谱峰的归属进行了校正,对质子化谱峰提出了新的指认.1理论计算方法和实验条件采用Gaussian 03程序包18中的B3LYP 方法得到了9-位氢腺嘌呤和7-位氢腺嘌呤的优化构型(图1;在相同的计算水平上,采用极化连续模型以水作溶剂优化了质子化腺嘌呤的5种低能构型(图1.所有计算中的C 、N 和H 原子均采用aug -cc -pVTZ 基组,其中cc -pVTZ 是Dunning 基组,aug 表

11、示包含弥散函数.为了进行振动分析,采用Pulay 标度量子力场方法19.首先在相同理论水平上计算了腺嘌呤和质子化腺嘌呤的直角坐标下的谐性力常数,然后将直角坐标下的力常数转化为内坐标下的力常数矩阵F DFT ij ,再对力常数进行校正:F scaled ij =(C i C j 1/2F DFT ij其中,C i 为第i 个坐标所对应力常数的标度因子,F DFT ij 是B3LYP 方法计算得到的内坐标力常数矩阵,F scaled ij为标度力常数矩阵.基于现有的计算方法和基组,本文中除对具有高频的C H 、N H 和NH 2的伸缩振动采用标度因子0.935外,对腺嘌呤分子的内坐标采用的标度因子

12、均为0.96320;由于我们是采用PCM 模型水作溶剂的条件下计算质子化腺嘌呤的优化结构和振动光谱的,比较其与腺嘌呤分子在酸性溶液中得到的拉曼光谱的数值,我们对其内坐标所采用的标度因子为0.990.然后按照GF (Wilson GF matrix矩阵方法得到标度后的振动频率和势能分布(PED,对振动模式进行了归属,振动分析中腺嘌呤分子的内对称坐标定义与Xue 等2定义的相同,而质子化腺嘌呤分子的定义则为s 1-s 16,s 18-s 28,s 30-s 33,s 35-s 42,分别与Xue 等2定义的s 1-s 16,s 17-s 27,s 28-s 31,s 32-s 39相同,与质子相关

13、的三个内坐标则分别定义为s 17=L (1,16,s 29=(2-1/2(2,1,16-(6,1,16,s 34=D (1,6,2,16.这里si 是内对称坐标,i 代表内对称坐标的编号,详细参照文献2.实验观测到的拉曼散射强度正比于微分拉曼散射截面.在谐振子模型近似下,微分拉曼散射截面可以表示为20-22:d i=(2445h 82c 軌i(軌0-軌i 41-exp(-hc 軌i/k BT S i上式中,c 、h 、k B 和T 分别是真空中的光速、Planck 常数、Boltzmann 常数和体系所处的绝对温度(298K,軌i 是第i 个振动模式的振动频率(cm -1,軌0是激发光的频率(

14、cm -1,而拉曼散射因子S i 来自B3LYP 计算.本论文中对于固体腺嘌呤分子(C 5H 5N 5,纯度99%,Alfa Aesar 公司拉曼光谱的测量采用法国Dilor 公司生产的LabRam I 型共焦显微拉曼光谱仪,使用He -Ne 激光器,激发线波长为632.8nm ,选择透镜为50倍,数值孔径(NA为0.55,工作距离(WD为8mm;对于腺嘌呤分子在pH=1的高氯酸(优级纯溶液中拉曼光谱的测量采用英国Renishaw 公司生产的UV -Vis R 1000共聚焦显微拉曼光谱仪,使用He -Cd 激光器,激发线波长为325.0nm,选择15倍紫外透镜,数值孔径为0.32,工作距离为

15、8.5mm.仪器的详细细节可参照文献23,24.2结果与讨论2.1相对稳定性和分子结构在B3LYP/aug -cc -pVTZ 水平下,腺嘌呤分子的两种构型(N9H 和N7H中,N9H 异构体比N7H 的能量低32.76kJ mol -1(表1.在同样水平下在极化连 图1腺嘌呤以及质子化腺嘌呤的构型以及原子标号Fig.1Structures and atomic numberings of N9H,N7H and five low -energy configurations of protonatedadenine1606No.8梁晓静等:腺嘌呤和质子化腺嘌呤的结构和振动光谱续模型以水作溶剂

16、的条件下得到了这两种构型在水溶液中的稳定结构,同样N9H 的能量更低,但与N7H 的能差减小为6.28kJ mol -1(表1.这主要是由于N7H 和N9H 之间偶极矩相差较大,其分别为2.3010-29和7.9110-30C m,导致在考虑溶剂化的条件下N7H 和N9H 之间的能差减少了26.48kJ mol -1(表1.因此,依据Boltzmann 分布,这两种异构体在水溶液中的比例约为12125,26,而在低温气相氩基质中,腺嘌呤主要以N9H 构型存在2-4.采用水作溶剂的极化连续模型得到腺嘌呤5种质子化异构体的构型如图1.表1列出了5种构型在气相和考虑PCM 模型时计算的相对能量.结果

17、表明在酸性溶液中N1位质子化的N9H 异构体(图1,P1N9是最稳定的,这与文献27的结果是一致的.表2列出的是计算得到的孤立的N9H 、水溶剂化下的N9H 和水溶剂化下的P1N9以及实验测到的腺嘌呤晶体中的结构参数,由于得到的实验数据的局限性,我们得到自由分子的C6N10键的键长比实验测得的长了0.0013nm.我们发现溶剂化效应或N1位质子化均使该键的键长缩短.理论计算该分子转动常数(A =2361.9,B =1576.1,C =947.3与在氩基质中测到的实验值(A =2371.9,B =1573.4,C =946.329接近.因此,理论结构与实验的偏差主要是由腺嘌呤晶体的分子间氢键作用

18、所引起的.质子化使腺嘌呤分子的结构发生了较大变化.从表2中的实验结果来看,质子化使N9H 结构中的共轭大键的电子云分布变形而使N1C2和N1C6键的成键作用削弱,从而键长增加,同时其力常数也发生了相应的减小.但与之相邻的C2N3和C6N10键长缩短,力常数增大.(C2N1C6以及其相邻的两个角(N1C2N3和(N1C6C5在质子化前后的变化最大(分别增加了4.4以及减少了3.5和3.3.质子化后P1N9氨基与嘌呤环几乎处于同一平面(D (N1C6N10H15=0.1,其相应的剪式振动模式的力常数(从106N m -1减小到48N m -1以及摇摆振动模式的力常数(从92N m -1减小到58N

19、 m -1也发生了明显减小.质子化在N1位使N9H 面外弯曲模式的力常数也增大,因此N1位质子化更有利于N9H 构型的稳定.2.2腺嘌呤分子的振动光谱对于腺嘌呤分子,我们对大部分基频的指认都与以前的研究工作相一致2-8.图2a 和2b 分别是我们理论模拟和实验观测的腺嘌呤分子的拉曼光谱.基于标度量子力场方法,我们在基频的归属上与Giese 和McNaughton 在2002年的归属取得了一致8.计算得到的基频与在低温氩基质下得到的实验红外光谱相比3,存在四个偏差超过10cm -1的谱峰,它们分别是与NH 2的剪式振动或者摇摆振动相关的谱峰,实验频率为1633、1599、1229和1005cm

20、-1,它们相应的理论计算标度频率分别是1623、1579、1219及990cm -1.我们认为其偏差主要是我们得到的C6N10键偏长引起理论上低估氮原子上孤对电表1在B3LYP/aug -cc -pVTZ 级别下计算得到的7种构型的相对能量值Table 1The relative energies of neutral and protonated adenine molecules calculated atB3LYP/aug -cc -pVTZ levelaThe unit of the energy is in kJ mol -1;busing PCM model;cThe total

21、energies of N9H and P3N7are -467.49634and -467.87011a.u.,respectively;dThe total energies of N9H and P1N9are -467.52400and -467.98089a.u.using the PCM model,respectively.Energy aNeutral adenine Protonated adenineN9H P1N9P3N7P3N9P7N9P1N7E 0c1.510c6.8435.5845.78E s b0d 0d1.647.113.749.59表2中性腺嘌呤以及质子化腺嘌

22、呤的优化几何结构aTable 2Optimized geometries of neutral adenine and protonated adenine aabond length in nm,bond angle and dihedral angle in degree;bB3LYP/aug -cc -pVTZ;cB3LYP/aug -cc -pVTZ combining the PCM model;dfrom Ref.28,and they are average values from 21crystal structures of adenine;efrom Ref.28,and

23、they are average values from13crystal structures of protonated adenine.L (N1C2L (C2N3L (N1C6L (C6N10(N1C2N3(N1C6C5(C2N1C6D (N1C6N10H14D (N1C6N10H15calc.1b0.13380.13320.13490.1337129.0117.6118.8protonated adenine calc.c0.1expt.e0.13620.13060.13600.1322125.5114.3123.2neutral adenine1607Acta Phys.-Chim

24、.Sin.,2021Vol.25 子与苯环的共轭作用,从而引起的NH 2剪式振动以及NH 2摇摆振动的力常数偏低导致的.同时峰1612cm -1(主要是由键C2N3伸缩振动(PED 32引起的和802cm -1(主要是由六元环扭转1(PED 45以及五元环扭转1(PED 32引起的也与理论值(1600和813cm -1有一定的偏差,其余的谱峰均与实验吻合.我们对于腺嘌呤在低温气相氩基质中红外光谱的指认与Nowak 等3在1996年、Xue 等2在2000年的指认进行比较,主要存在差别有如下五点:(1我们认为实验得到的红外峰927cm -1主要是由五元环畸变I(PED 76引起的,这与Xue 等

25、一致;而Nowak 等认为其是由五元环畸变II 、五元环畸变I 以及C4C5伸缩振动引起的.考虑到该分子的9-位氮以及7-位氮的同位素位移分别为7.7和12.3cm -1,该频率的同位素位移与实验得到的同位素位移7.5和12.5cm -13相符,因此,将该峰的指认为五元环畸变I 更为合理;(2我们认为实验得到的红外峰655cm -1是由六元环扭转II(PED 62贡献的,这与Xue 等2的一致,而Nowak 等3认为其是由六元环扭转I 和六元环扭转II 引起的;(3实验得到的两个较强红外峰583和591cm -1,Xue 等2认为由基质效应30引起分裂的红外峰是由NH 2扭转振动的基频,而我们

26、并没有相应基频与之相对应,并且我们得到的由NH 2扭转引起的振动峰红外强度非常弱,因此该谱峰产生的原因需要更进一步的研究来确认;(4Xue 等2认为实验得到的红外峰566cm -1是由N9H 面外弯曲以及六元环扭转引起的,而我们认为除以上之外,五元环扭转(PED 30应更为重要;(5基频峰522和521cm -1的振动频率非常接近,我们将其指认为N9H 面外弯曲和六元环畸变II 贡献,但二者相应PED 值分别是39和62及40和34.我们将它们分别指认为红外谱峰513cm -1和拉曼谱峰536cm -1,这与Xue 2和Nowak 3等的指认均不一致.2.3质子化腺嘌呤的振动光谱前人对腺嘌呤分

27、子在酸性条件下可能存在的质子化结构进行了计算27,也有人从理论1和实验9方面对酸性介质中腺嘌呤分子振动光谱进行指认,但由于受当时理论计算水平限制,质子化腺嘌呤的振动分析仍有待改进.我们分析了5种质子化构型,并对能量最低的构型(P1N9进行振动分析,结果如表3所示.表3总结了利用振动分析所得到在500-1600cm -1的31个基频及其势能分布.实验红外频率11和拉曼频率(图2(c 、d也列于表3中.比较频率变化可以发现质子化对1000cm -1以上的谱峰影响较大,N1H 面内弯曲(s 29振动模式分散到该区域中的其中4个基频中,而且基频都发生了重新组合.但是质子化对1000cm -1以下的谱峰

28、影响较小,N1H 面外弯曲(s 34振动模式基频是721cm -1.实验谱中与该振动模相对应的谱峰常被环呼吸振动引起的很强的红外(722cm -1和拉曼谱谱(720cm -1所掩盖.以下是我们对得到的振动光谱进行分区归属.在区域1700-1550cm -1区域中理论上有4个基频振动,实验得到的红外光谱中则发现了5个很强的红外峰,预测到有4个基频峰与其中的4个峰对应.而对于观测峰1650cm -1,文献11认为其是由构型P3N7的N7H 伸缩振动引起的,而我们通过振动分析P3N7构型预测到基频峰1657cm -1具有很强的红外活性,且其是由C5C6、C6N10伸缩振动引起的,因此我们确认酸性介质

29、中有少量的其他构型(P3N7和P3N9的存在31.而对于观测到的拉曼光谱而言,由于在稀的水溶液中对于拉曼光谱较低的分辨,在这个区间主要有1616和1584cm -1两个基频峰,其中峰1616cm -1主要来自CN 键伸缩振动和NH 2剪式振动,而峰1584cm -1是由NH 2剪式振动以及N1H 面外弯曲引起的,与中性分子晶体的峰1596cm -1相比有一定的红移.在该区域1550-1300cm -1有7个基频振动,其中在实验得到的红外谱中观测到6个基频峰,而在实验拉曼光谱中仅观测到3个基频峰.拉曼观测峰1502cm -1是由N3C4伸缩、N7C8伸缩以及C8H 面内弯曲图2腺嘌呤和质子化腺嘌

30、呤分子的理论以及实验拉曼光谱Fig.2Calculated and experimental Raman spectra of adenine and protonated adenine molecules(acalculated using B3LYP/aug -cc -pVTZ for adenine molecule;(bsolid powder with excitation at 632.8nm for adenine molecule;(ccalculated using B3LYP/aug -cc -pVTZ and PCM model for protonated adeni

31、ne molecule P1N9;(dRaman spectrumof protonated adenine molecules dissolved in HClO 4(0.1molL -1with excitation at 325.0nm1608No.8梁晓静等:腺嘌呤和质子化腺嘌呤的结构和振动光谱振动引起的,与中性腺嘌呤晶体的拉曼观测峰1481cm -1相比有一定的蓝移.观测峰1352cm -1主要是由C8H 面内弯曲和C8N7伸缩振动引起的,与中性腺嘌呤晶体的拉曼观测峰1470cm -1相比发生显著的红移.峰1313cm -1主要是由C5N7和C6N1伸缩振动引起的,与中性腺嘌呤晶体的

32、拉曼观测峰1331cm -1相表3计算以及观察到的质子化腺嘌呤的振动频率及其指认Table 3Theoretical and experimental vibrational frequencies and assignments of protonated adenineaIR of protonated adenine in gas from Ref.11;bRaman Shift of adenine in HClO 4solution (pH=1from Fig.2c;cthe calculated frequencies ofvibrational spectrum;dthe cal

33、culated IR intensity;ethe calculated Raman scattering factor;fthe calculated frequencies of vibrational spectrum using the scaled quantum mechanical field (SQMFprocedure;gabbreviations used:str,stretching;bend,bending;def,deformation;sciss,scissoring;rock,rocking;tor,torsion;o -o -p,out -of -plane;s

34、,symmetric;as,antisymmetric.IR a /cm -1Raman b /cm -1Theo c /cm -1I d /(km mol -1104S i e /(nm 4amu -1scaled /cm -1PED Mode description g -3480283.798.63329s 16(98NH 2as str -3370632.3495.83224s 15(99NH 2s str -3229616.8410.13098s 14(99NH str -315629.7285.03053s 13(99CH str -3138624.5284.93013s 12(9

35、4CH str -311417.0212.63000s 17(93NH str 1690-16911366.7 4.81682s 11(34,s 5(32CN str+CH bend 1416-144035.5119.71433s 8(25,s 10(24,s 3(16,s 2+s 4+s 22(27CN str1226-125162.08.51244s 25(31,s 8(19,s 6(10CH bend+CN str 1172-117945.615.11173s 1(25,s 27(18,s 29(12CN str+NH 2rock 1128-114312.028.31137s 1(21,

36、s 3(10,s 10(10CN str1066-108310.126.31077s 9(50,s 26(33,s 25(10CN str+NH bend 1004-1002 1.912.2997s 27(52,s 6(23,s 1(8NH 2rock -954 1.9 1.7949s 31(92CH o -o -p 947-94533.528.6940s 21(70,s 4(11ring -str -8998.7 3.4894s 32(109,s 39(27NH o -o -p 902-89737.817.7892s 18(53,s 20(13ring -str 812-81710.0 3.

37、4813s 38(43,s 35(40,s 30(23,s 36(14ring -tor -72492.60.9721s 34(63,s 30(24,s 35(14N +H o -o -p 72272071313.789.7709s 3(20,s 4+s 5+s 6+s 7+s 10+s 20+s 21(48CN str 697-68056.40.4676s 39(48,s 30+s 34+s 37(24ring -tor 646-65046.30.2647s 39(50,s 30(28,s 33(21ring -tor 622616621148.90.2618s 33(99,s 38(46,

38、s 39(21N +H o -o -p 580-5967.135.0593s 22(39,s 20(28,s 5(22ring -str 566-568 2.4 2.1565s 37(13,s 33(12,s 35(12ring -tor 53052652827.07.7525s 20(22,s 23(20ring -str -51313.30.2510s 41(92NH 2tor 509503511 1.27.7508s 19(66ring -str -384317.6 2.8382s 42(107NH 2o -o -p -2860.20.1284s 36(35,s 40(30,s 37(1

39、3,s 34(10ring -tor -27113.28.0270s 23(62,s 20(19ring -str -19950.40.0198s 40(53,s 36(23,s 37(11ring -tor -15219.00.8151s 37(64,s 35(39,s 36(18ring -tor16091610 Acta Phys. -Chim. Sin., 2021 Vol.25 比有一定的红移. 同时我们发现该区域谱峰的相对 强度与理论相差较大, 这是由于实验 激 发 光 能 量 (325 nm很靠近质子化腺嘌呤的基态与第一激发态 能差所产生的预共振效应造成的31,32. 在区域 1

40、300-1000 cm -1 有 4 个基频振动, 其中 在红外谱中观测到 4 个红外基频峰, 但没有较强的 拉曼谱峰与之对应. 值得注意的是在红外谱中 1290 cm -1 峰属于很强的谱峰, 但在 P1N9 中没有相应的 基频与之相对应. 但我们在对 P3N7 结构的振动分 析得到了 1297 cm-1 且具有较强的红外活性. 最后, 在区域 1000-500 cm -1 有 16 个基频, 目 前在红外光谱中观测到 12 个红外基频峰, 而在拉曼 光谱中仅观测到 4 个拉曼基频峰. 在该区域的红外 峰强度较弱, 且主要以面外振动模式居多. 图 2d 中 拉曼谱峰 935 cm-1 是来自

41、溶液中的高氯酸根的全对 称振动 , 它掩盖了由五元环畸变所引起的一个较 弱的基频峰 940 cm -1. 拉曼谱峰 720 cm -1 和红外谱 24 2 Xue, Y.; Xie, D. Q.; Yan, G. S. Int. J. Quantum Chem., 2000, 76: 686 3 Nowak, M. J.; Lapinski, L.; Kwiatkowski, J. S.; Leszczynski. J. J. Phys. Chem., 1996, 100: 3527 4 Nowak, M. J.; Rostkowska, H.; Lapinski, L.; Kwiatkows

42、ki, J. S.; Leszczynski. J. Spectrochim. Acta A, 1994, 50: 1081 5 Stepanian, S. G.; Sheina, G. G.; Radchenko, E. D.; Blagoi, Y. P. J. Mol. Struct., 1985, 131: 333 6 Nowak, M. J.; Lapinski, L.; Kwiatkowski, J. S.; Leszczynski. J. Spectrochim. Acta A, 1991, 47: 87 7 Wirkiewicz-Kuczera, J.; Karplus, M.

43、J. Am. Chem. Soc., 1990, 112: 5324 8 Giese, B.; McNaughton, D. J. J. Phys. Chem. B, 2002, 106: 101 9 Ivanova, B. B. Spectrosc. Lett., 2005, 38: 635 10 11 12 13 14 15 16 17 18 Dkhissi, A.; Houben, L.; Adamowicz, L.; Maes, G. Chem. Phys. Lett., 2004, 387: 362 Majoube, M. J. Raman Specrrosc., 1985, 16:

44、 98 Pople, J. A.; Gordon, M. H.; Fox, D. J.; Raghavachari, K.; Curtiss, L. A. J. Chem. Phys., 1989, 90: 5622 Curtiss, L. A.; Jones, C.; Trucks, G. W.; Raghavachari, K.; Pople, J. A. J. Chem. Phys., 1990, 93: 2537 Becke, A. D. J. Chem. Phys., 1993, 98: 5648 Lee, C.; Yang, W.; Parr, R. G. Phys. Rev. B

45、, 1988, 37: 785 Miehlich, B.; Savin, A.; Stoll, H.; Preuss. H. Chem. Phys. Lett., 1989, 157: 200 Woon, D. E.; Dunning Jr., T. H. J. Chem. Phys., 1993, 98: 1358 Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; et al. Gaussian 03, Revision C.02. Wallingford, CT: Gaussian Inc., 2004 19 20 21 22 23 24 25

46、26 27 28 29 30 31 Pulay, P.; Fogarasi, G.; Pang, F.; Boggs, J. E. J. Am. Chem. Soc., 1979, 101: 2550 Wu, D. Y.; Ren, B.; Jiang, Y. X.; Xu, X.; Tian, Z. Q. J. Phys. Chem. A, 2002, 106: 9042 Wu, D. Y.; Liu, X. M.; Xu, Y. C. J. Lightscatting, 2006, 4: 323 Wu, D. Y.; Hayashi, M.; Tian, Z. Q.; Lin, S. H.

47、 Spectrochim. Acta A, 2004, 60: 137 Cui, L.; Wang, A.; Wu, D. Y.; Ren, B.; Tian, Z. Q. J. Phys. Chem. C, 2021, 112: 17618 Cui, L. Ph. D. Dissertation. Xiamen: Xiamen University, 2021 崔 丽. 博士学位论文. 厦门: 厦门大学, 2021 Chenon, M. T.; Pugmire, R. J.; Grant, D. M.; Panzica, R. P.; Townsend, L. B. J. Am. Chem.

48、 Soc., 1975, 97: 4636 Eastman, J. W. Ber. Bunsen-Ges. Phys. Chem., 1969, 73: 407 Taylor, R.; Kennard, O. J. Mol. Struct., 1982, 78: 1 Gu, J. D.; Jerzy, L. J. Phys. Chem. A, 1999, 103: 2744 Brown, R. D.; Godfrey, P. D.; McNaughton, D.; Pierlot, A. P. Chem. Phys. Lett., 1989, 156: 61 Sponer, J.; Hobza, D. J. Phys. Chem., 1994, 98: 3161 Marian, C.; Nolting, D.; Weinkauf, R. Phys. Chem. Chem. Phys., 2005,