2,4-二硫代尿嘧啶的紫外吸收光譜和共振拉曼光譜

金穎淳 鄭旭明

(浙江理工大學化學系，杭州 310018)

2,4-二硫代尿嘧啶的紫外吸收光譜和共振拉曼光譜

金穎淳 鄭旭明＊

(浙江理工大學化學系，杭州 310018)

硫代嘧啶堿基是光動力療法潛在的重要光敏劑，其最低單重激發(fā)態(tài)的光物理研究已有廣泛報道。然而，其較高激發(fā)態(tài)的躍遷性質和反應動力學研究較為稀少。因此，本文采用共振拉曼光譜和密度泛函理論計算方法研究2,4-二硫代尿嘧啶的紫外光譜和幾個較高單重激發(fā)態(tài)的短時結構動力學。首先，基于共振拉曼光譜強度與電子吸收帶振子強度f的關系，將紫外光譜去卷積成四個吸收帶，分別為358 nm (f =0.0336)中等強度吸收帶(A帶)，338 nm (f = 0.1491)、301 nm (f = 0.1795)和278 nm (f = 0.3532)強而寬的吸收帶(B、C和D帶)。這一結果既吻合密度泛函理論計算結果，又符合共振拉曼光譜強度模式對紫外光譜帶的預期。據此，去卷積得到的四個吸收帶被分別指認為S0→ S2躍遷、S0→ S6躍遷、S0→ S7躍遷和S0→ S8躍遷。同時，分別對B，C和D帶共振拉曼光譜進行了詳細的指認，獲得了短時動力學信息。結果表明，S8態(tài)短時動力學的顯著特征是在Franck-Condon區(qū)域或附近發(fā)生了S8(ππ＊)/S(nπ＊)勢能面交叉引發(fā)的、伴隨超快結構扭轉的非絕熱過程。S7和S6態(tài)短時動力學的主要特征是反應坐標的多維性，它們分別沿C5C6/C2S8/C4S10/N2C3+ C4N3H9/N1C2N3/C2N1C6/C6N1H7/C5C6H12和C5C6/N3C2/C4S10/C2S8+C6N1H7/C5C6H12/C5C6N1/C5C6H12/C2N1C6/N1C2N3/C4N3H9/N1C2N3等內坐標演化。

2,4-二硫代尿嘧啶；激發(fā)態(tài)結構動力學；紫外吸收光譜；共振拉曼光譜；密度泛函理論

1 引 言

天然核酸堿基受光誘導的激發(fā)態(tài)衰變動力學是一個由 S2(ππ*)/S0勢能面交叉主導的超快內轉換過程1–5。這一過程能使激發(fā)態(tài)分子來不及發(fā)生化學反應就回到了基態(tài)，從而保障了生命體中的DNA分子免受紫外輻射影響，使遺傳信息復錄得以正常進行。與天然核酸堿基不同，硫代核酸堿基的低能態(tài)激發(fā)態(tài)衰變過程主要受系間竄越過程的控制6–16。例如，2-硫代胸腺嘧啶初始布居單重態(tài)的主要弛豫通道為系間竄越到 T1(ππ*)6–9，2-硫代胸腺嘧啶6,8,9、4-硫代胸腺嘧啶8,10–13和硫代鳥嘌呤14等分子經系間竄越到到T1(ππ*)態(tài)的量子產率很高，有些接近1。硫代堿基具有較高的光敏化單線態(tài)氧產率，具體與結構有關8。這使得它們在光動力療法等領域倍受關注。

在動力學方面，納秒、皮秒和飛秒時間分辨瞬態(tài)吸收光譜技術揭示了2-硫代胸腺嘧啶T1(ππ*)形成的時間尺度9,11–13。理論計算研究洞察了硫代尿嘧啶和硫代鳥嘌呤系間竄越過程的微觀機制16,17。借助于全活化空間自洽場(CASSCF)理論及其二階微擾(CASPT2)方法所獲得的2-硫代尿嘧啶激發(fā)態(tài)勢能面和旋-軌耦合信息，崔剛龍和方維海提出了從初始布居態(tài) S2(ππ*)產生最低三重態(tài) T1(ππ*)的三條互相競爭、高效的非絕熱通道17。Pollum and Crespo-Hernández18采用飛秒時間分辨寬帶瞬態(tài)吸收實驗評估了崔剛龍和方維海提出的三條ISC 通道，提出了 S2(ππ*) → S1(nπ*) → T1(ππ*)過程是主要通道。這一結論與前期Gobbo的研究結果一致19。目前認為，S1(nπ*)態(tài)是完成 S2(ππ*) →T1(ππ*)系間竄越的門態(tài)，決定著 T1(ππ*)量子產率大小20,21。姜杰等22采用共振拉曼光譜結合CASSCF方法研究了 2-硫代尿嘧啶激發(fā)態(tài)衰變動力學，揭示經 S2(ππ*)/S1(nπ*)勢能面交叉從 S2(ππ*)態(tài)弛豫到S1(nπ*)態(tài)是產生T1(ππ*)態(tài)的主要通道。最近，崔剛龍等采用 CASPT2//CASSCF和QM(CASPT2// CASSCF)/MM方法提出了2,4-二硫代胸腺嘧啶(2,4-diiothymine)在真空、微溶劑化和水環(huán)境下的激發(fā)態(tài)衰變機制和溶劑效應23。

2,4-二硫代尿嘧啶(簡稱 24DTU)是尿嘧啶分子中兩個羰基氧原子均被硫原子取代的產物。相對于單硫代堿基，二硫代堿基的電子態(tài)結構較為復雜，紫外光區(qū)多個光吸收態(tài)的存在使得紫外光譜指認變得不確定。因此，本文擬借助共振拉曼光譜的強度模式對不同電子態(tài)的對應關系，通過紫外吸收光譜和共振拉曼光譜實驗，對比密度泛函理論計算結果，研究24DTU分子激發(fā)態(tài)電子結構，指認紫外吸收光譜和共振拉曼光譜，以獲得電子態(tài)耦合信息和短時結構動力學結果。

2 實驗部分

2.1 試劑和儀器

24DTU (2,4-dithiouracil) (99.9%，北京百靈威科技有限公司)，乙腈(99.9%，斯百全化學有限公司，甲醇(99.9%，斯百全化學有限公司，純水(杭州娃哈哈集團有限公司)。

紫外可見分光光度計(Cary 50，瓦里安公司，美國)，傅里葉紅外光譜儀(Nicolet Avatar 370，尼高力公司，美國)；傅里葉拉曼光譜儀(Nicolet Raman 960，尼高力公司，美國)；共振拉曼光譜儀(自制)。

2.2 實驗方法

24DTU在乙腈、甲醇和水中的濃度約為0.90–1.0 × 10?3mol·L?1。共振拉曼光譜實驗裝置和方法參見文獻24。24DTU的共振拉曼光譜的激發(fā)波長(266.0、273.9、282.4、299.1、309.1、319.9、341.5和354.7 nm)由10 Hz Nd:YAG激光器的二、三、四倍頻激光(266.0、355和532 nm)和它們的氫受激拉曼線提供。樣品溶液由循環(huán)泵輸送至噴嘴口，噴嘴口可以擠壓溶液使得通過的溶液變成薄膜狀。此時，通過棱鏡折射后的激光打在溶液薄膜上，橢球鏡則收集拉曼散射的信號，最后Princeton Instruments PyLoN:2KBUV型CCD檢測器(2048 ×512前感光液氮冷卻)收集信號并將它們顯示在計算機上。噴嘴口流出的溶液則通過漏斗收集，由導管通入進錐形瓶，產生溶液流動的循環(huán)。共振拉曼光譜的一次采集時間由強度或整體信噪比決定，一般在100–150 s。50–100次采集數據的累加結果即為該激發(fā)波長下的共振拉曼光譜原圖。原圖經像素-波數轉換(cm?1)、靈敏度和自吸收等校正、溶劑峰扣除等處理后得到樣品的共振拉曼光譜圖。

圖1 (a) 24DTU在乙腈、甲醇和水中的紫外吸收光譜；(b) 在PCM溶劑模型下由B3LYP/6-31+G(d)計算獲得的24DTU的幾何結構示意圖Fig.1 (a) UV absorption spectra of 24DTU in acetonitrile, methanol and water; (b) Schematic diagram of the geometry structure of 24DTU.

3 理論計算

電子基態(tài)幾何結構和振動頻率采用密度泛函理論在極化連續(xù)介質模型(PCM)溶劑模型和B3LYP/6-31+G(d)計算水平下計算獲得。電子躍遷能和振子強度在 B3LYP-TD/6-311++G(3df,3pd)水平下計算獲得。本文所有量子化學計算均采用Gaussian 09W程序包完成25。簡正模式分析采用VEDA4程序完成26。

4 結果與討論

4.1 紫外光譜研究

圖1示出24DTU在乙腈、甲醇和水溶液中的紫外吸收光譜，圖中箭頭上方數字為共振拉曼光譜實驗所用的激光波長。24DTU紫外光譜展示兩個寬吸收帶，分別位于282和346 nm，其中346 nm較強吸收帶的不對稱性預示該帶由某些強吸收帶組合。表 1列出 B3LYP-TD/6-311++G(3df,3pd)計算水平和PCM溶劑模型下獲得的24DTU在乙腈中的電子躍遷能、躍遷軌道和振子強度。表 1可見，在 > 240 nm光譜區(qū)域，計算結果給出四個允許躍遷，最大吸收波長和振子強度f分別為336 nm(f = 0.0451)，297 nm (f = 0.1472)，272 nm (f =0.3851)和253 nm (f = 0.4066)。圖2為24DTU在乙腈溶液中的紫外吸收光譜的去卷積圖。圖 2和表 1可見，去卷積獲得的四個主要吸收峰的最大吸收波長和振子強度 f與理論計算結果又較好的相關性。因此，根據圖1、圖2和表1的結果，圖2中的358 nm (f = 0.0336)中等強度吸收帶(A帶)被指認為S0→ S2躍遷, 338 nm (f = 0.1491)、301 nm (f = 0.1795)和278 nm (f = 0.3532)強而寬的吸收帶(B、C 和 D 帶)被分別指認為 S0→ S6、S0→ S7和S0→ S8躍遷。

值得指出的是，為了排除去卷積過程存在的隨意性，綜合考慮了不同激發(fā)波長下的共振拉曼光譜與振子強度的關系，即拉曼光譜強度與振子強度f的平方成正比。因此，把338 nm去卷積帶(f =0.1491)指認為S0→ S6吸收帶(f = 0.1789)而非S0→S2吸收帶(f = 0.0541)。這得到了341.5和354.7 nm共振拉曼光譜強度模式的有力支持。此外，我們注意到，構成S0→ S6、S0→ S7和S0→ S7躍遷的主要軌道基本相同，區(qū)別在于相應的躍遷軌道系數大小或符號不同，這表明S6、S7、和S8三個態(tài)之間存在態(tài)態(tài)耦合。溶劑對紫外吸收光譜的結構和帶寬影響不大。根據朗伯比爾定律(Lambert-Beer)定律，24DTU在乙腈、甲醇和水溶液中的最大吸收波長處的摩爾消光系數εmax/λmax分別為：2.29 × 104L·mol?1·cm?1(280.0 nm)、2.21 × 104L·mol?1·cm-1(283.0 nm)和2.20 × 104L·mol?1·cm?1(281.0 nm)。

4.2 振動光譜指認

圖3示出24DTU的傅里葉變換拉曼和紅外光譜，以及由B3LYP/6-31+G(d)計算獲得并由GaussView軟件輸出的拉曼光譜(參數為半峰寬2 cm?1，步長1 cm?1)27。圖3中，計算拉曼光譜頻率已按實驗光譜頻率做了校正。表2列出實驗和計算振動頻率。表2可見，24DTU有30個振動模，其中屬于A′和A′不可約表示的振動模分別有21個和9個。圖3和表2可見，拉曼光譜中主要譜帶的實驗拉曼位移(頻率)與計算值基本吻合，表明密度泛函理論計算結果總體合理。例如，實驗拉曼光譜中的1549、1491、1253、1228、1188、1118、1077、683、460、443和229 cm?1等譜帶與計算光譜中的ν6、ν7、ν10、ν11、ν12、ν13、ν14、ν17、ν18、ν19、ν21等有良好的一一對應關系。基于密度泛函理論計算結果開展了振動光譜的簡正模式分析，獲得的勢能分布(potential energy distribution)列于表2。表2可見，除NH/CH伸縮振動外，絕大多數振動模的振動形式都由多個內坐標的復合運動構成。

圖2 24DTU在乙腈溶液中的紫外吸收光譜的去卷積曲線(點線)Fig.2 Deconvoluted curves (dotted line) of the UV absorption spectrum of 24DTU in acetonitrile.

表1 在PCM溶劑模型下由B3LYP-TD/6-311++G(3df,3pd)計算獲得的24DTU的電子躍遷能(△E)、躍遷軌道和振子強度(f)Table 1 Electronic transition energies (△E), transition orbitals and Oscillator strengths (f) for 24DTU computed by B3LYP-TD/6-311++G(3df,3pd) and polarization continum model (PCM).

4.3 共振拉曼光譜和短時結構動力學

圖 4示出 266.0、273.9、282.4、299.1、309.1、319.9、341.5和 354.7 nm 激發(fā)下 500–1800 cm?1光譜區(qū)域內的共振拉曼光譜，以展現(xiàn)基頻信息。

圖3 24DTU的傅里葉變換紅外和拉曼光譜以及B3LYP/6-31+G(d)計算拉曼光譜Fig.3 Comparison of FT-IR, FT-Raman and B3LYP/6-31+G(d) computed Raman spectra of 24DTU.

表2 24DTU在B3LYP/6-31+G(d)計算水平下的計算頻率和實驗觀察到的傅里葉變換紅外和拉曼光譜的振動頻率和指認Table 2 Vibrational frequencies and assignments of 24DTU calculated by using B3LYP/6-31+G(d) and PCM model and observed by experimental FT-IR and FT-Raman spectra.

圖5展示 500–3500 cm?1光譜區(qū)域內的266.0、319.9和341.5 nm共振拉曼光譜，以全面展現(xiàn)倍頻和組合頻信息。借助圖3和表2，圖4和圖5中的共振拉曼光譜可被初步指認為由 ν17、ν16、ν15、ν14、ν13、ν11、ν9、ν8、ν7、ν6和 ν5等 11 個振動基頻以及它們的倍頻和組合頻構成。為了排除因溶劑扣減引起的瑕疵峰對共振拉曼光譜的干擾和考察溶劑性質對24DTU共振拉曼光譜的強度模式影響，獲取了在乙腈、甲醇和水中的共振拉曼光譜。圖4可見，溶劑扣減的瑕疵峰和溶劑性質對共振拉曼光譜影響不大。

continued Table 2

圖4 2.4-二硫代尿嘧啶在乙腈(a)、甲醇(b)和水(c)溶劑中的不同激發(fā)波長下的共振拉曼光譜Fig.4 Resonance Raman spectra of 24DTU in acetonitrile (a), methanol (b) and water (c) solvents.

圖 4可見，在乙腈和甲醇中不同激發(fā)波長下的共振拉曼光譜存在良好的一一對應關系，表明24DTU的 B、C和 D帶吸收躍遷的能量和Franck-Condon區(qū)域勢能面構造在兩種溶劑環(huán)境中基本相同。但水對24DTU的C吸收帶產生明顯的藍移作用，并且相應的共振拉曼光譜強度模式與在乙腈和甲醇中的有一定的差異，表明C吸收帶的躍遷能和對應勢能面在 Franck-Condon區(qū)域的構造在水環(huán)境中發(fā)生了明顯變化。

眾所周知，共振拉曼光譜的強度模式對激發(fā)態(tài)性質很敏感。激發(fā)態(tài)性質不同，所對應的共振拉曼光譜的強度模式往往也不同。因此，通過改變激發(fā)波長使之與不同激發(fā)態(tài)共振，能獲得不同激發(fā)態(tài)上的短時結構動力學信息。以水中共振拉曼光譜為例，圖4(c)中341.5和354.7 nm共振拉曼光譜的強度模式基本相同，266.0、273.9和282.4 nm共振拉曼光譜的強度模式基本相同，而 309.1和319.9 nm共振拉曼光譜的強度模式基本相同，但與前兩者都不同。顯然，341.5和354.7 nm激發(fā)波長與S6態(tài)共振，因此，其共振拉曼光譜的強度模式體現(xiàn)了波包沿著 S6態(tài)的Franck-Condon區(qū)域勢能面演化的超快動力學過程。同理，266.0、273.9和282.4 nm激發(fā)波長與S8態(tài)共振，309.1和319.9 nm激發(fā)波長與S7態(tài)共振，它們對應的共振拉曼光譜強度模式體現(xiàn)了在相應激發(fā)態(tài)的Franck-Condon區(qū)域的超快動力學過程。

在共振拉曼光譜中，由于高電子激發(fā)態(tài)的非共振(off-resonance)或預共振(pre-resoncne)效應28、多電子態(tài)間躍遷矩的非一致性相干28或振動-電子耦合等效應29，高電子激發(fā)態(tài)會對低激發(fā)態(tài)共振拉曼光譜的基頻產生額外的增強或減弱效應。圖 4可見，266.0、273.9和282.4 nm共振拉曼光譜中最強的拉曼峰為ν13和ν23，其中ν23峰強度延伸至309.1 nm共振拉曼光譜中，而ν13峰強度一直延伸至354.7 nm共振拉曼光譜中。這表明S8態(tài)對309.1和319.9 nm共振拉曼光譜的ν23峰產生明顯的預共振增強效應，對309.1、319.9、341.5和354.7 nm共振拉曼光譜的ν13峰產生明顯的預共振增強效應。同理，S7對341.5和354.7 nm共振拉曼光譜的ν5和ν11基頻也會有一定的預共振增強效應。然而，高電子激發(fā)態(tài)會對低激發(fā)態(tài)共振拉曼光譜的倍頻和組合頻一般不產生預共振效應或影響不大。30因此，扣除預共振效應、排除溶劑扣減的瑕疵峰和可能的激光線影響后，從各共振拉曼光譜中能得到如下信息：

圖5 24DTU在乙腈溶劑中的266.0 nm (上)、319.9 nm(中)和 341.5 nm (下)共振拉曼光譜Fig.5 266.0 nm (top), 319.9 nm (middle) and 341.5 nm(bottom) resonance Raman spectra of 24DTU in acetonitrile.

S8態(tài)短時結構動力學：266.0、273.9和282.4 nm 共振拉曼光譜可指認為 ν23、ν17、ν16、ν15、ν14、ν13、ν12、ν11、ν9、ν8、ν7、ν6和 ν5等 13 個振動基頻以及它們的倍頻和合頻，其中，N1C6C5H11/C4C5C6H12扭轉振動模 ν23、C2S8/C4S10伸縮 +C2N1C6面內彎曲振動模 ν16、C2S8/C4S10伸縮 +C4N3H9/N1C2N3面內彎曲振動模ν13、N1C2伸縮 +C2N1C6/N1C2N3/C4N3C2面內彎曲振動模 ν17、N1C6/N3C4伸縮+C6C3H11/C5C6H12面內彎曲振動模ν12和C5C6N1/C5C6H12/ C2N1C6/N1C2N3面內彎曲 +N1C6伸縮振動模ν15振動模的基頻、倍頻和合頻，及其它們與其他振動模的合頻強度占據共振拉曼光譜總強度的絕大部分。這表明，S8態(tài)短時結構動力學主要沿 N1C6C5H11/C4C5C6H12、C2S8/C4S10/N1C6/N3C4和C2N1C6/C4N3H9/N1C2N3/C5C6N1/C5C6H12/C2N1C6等多維內坐標展開。ν23歸屬 A′不可約表示，在S8(ππ*)態(tài)絕熱動力學過程中為非拉曼活性的。因此，根據拉曼強度的振動電子耦合理論31,32和非絕熱體系共振拉曼光譜研究結果33–36，強ν23(A′)振動模的出現(xiàn)預示 24DTU 在 S8(ππ*)態(tài)的Franck-Condon區(qū)域或附近發(fā)生了與其它勢能面交叉的非絕熱過程，其結構動力學涉及極快速的結構扭轉運動。

S7態(tài)短時結構動力學：309.1和319.9 nm共振拉曼光譜可指認為 ν17、ν16、ν15、ν14、ν13、ν12、ν11、ν9、ν8、ν7、ν6和 ν5等 12 個振動基頻以及某些較強基頻的倍頻和組合頻。與S8態(tài)短時結構動力學不同，S7(ππ*)態(tài)短時動力學不含 A′不可約表示的振動坐標，表明在Franck-Condon區(qū)域及附近沒有S(nπ*)態(tài)的強力耦合。重要的是，309.1和319.9 nm共振拉曼光譜的強度模式與266.0、273.9和282.4 nm共振拉曼光譜的明顯不同，其強度模式主要由C5C6伸縮振動模 ν5、C2S8/C4S10伸縮 + C4N3H9/N1C2N3面內彎曲振動模 ν13、C2S8/C4S10伸縮 +C2N1C6面內彎曲振動模 ν16和 N2C3/C2S8/C4S10伸縮 + C6N1H7/C5C6H12面內彎曲振動模ν11等構成，說明 S7(ππ*)態(tài)短時動力學主要沿 C5C6/C2S8/C4S10/N2C3和 C4N3H9/N1C2N3/C2N1C6/C6N1H7/C5C6H12等內坐標演化。

S6態(tài)短時結構動力學：341.5和354.7 nm共振拉曼光譜可指認為 ν17、ν16、ν15、ν14、ν13、ν10、ν9、ν8、ν7、ν6和 ν5等 11 個振動基頻以及它們的倍頻和合頻。它們的強度模式與與S7或S8態(tài)共振拉曼光譜強度模式都不同，其短時結構動力學主要沿 C5C6/N3C2/C4S10/C2S8和 C6N1H7/C5C6H12/C5C6N1/C5C6H12/C2N1C6/N1C2N3/C4N3H9/N1C2N3等內坐標演化。

5 結 論

采用共振拉曼光譜技術結合密度泛函理論計算方法研究了 24DTU 的電子激發(fā)和Franck-Condon區(qū)域結構動力學，獲得如下結論：

(1) 共振拉曼光譜強度模式與去卷積獲得的四個紫外吸收帶的振子強度 f間存在一一對應關系。中等強度的350 nm吸收帶被指認為S0→ S2躍遷(A帶吸收), 338 nm (f = 0.1491)、301 nm (f =0.1795)和278 nm (f = 0.3532)強而寬的吸收帶(B、C 和 D帶)被分別指認為 S0→ S6、S0→ S7和 S0→S8躍遷。

(2) B帶、C帶和D帶共振拉曼光譜可被分別指認為由11、12和13個振動基頻以及它們的倍頻和合頻。在乙腈和甲醇中，341.5和354.7 nm共振拉曼光譜主要歸屬于 B帶吸收或 S0→ S6躍遷；309.1和319.9 nm共振拉曼光譜主要歸屬于C帶吸收或S0→ S7躍遷；282.4、273.9和266.0 nm共振拉曼光譜主要歸屬于D帶吸收中或S0→ S8躍遷。

(3) S6、S7和S8態(tài)短時動力學主要沿多維內坐標展開。S8態(tài)短時結構動力學最重要的特征是在Franck-Condon區(qū)域或附件發(fā)生了S8(ππ*)/S(nπ*)勢能面交叉引發(fā)的、伴隨超快結構扭轉運動的非絕熱過程。S7和S6激發(fā)態(tài)短時動力學分別沿C5C6/C2S8/C4S10/N2C3+ C4N3H9/N1C2N3/C2N1C6/C6N1H7/C5C6H12和C5C6/N3C2/C4S10/C2S8+ C6N1H7/C5C6H12/C5C6N1/C5C6H12/C2N1C6/N1C2N3/C4N3H9/N1C2N3等內坐標演化。

(1) Crespo-Hernández, C. E.; Cohen, B.; Hare, P. M.; Kohler, B. Chem.Rev. 2004, 104 (4), 1977. doi: 10.1021/cr0206770

(2) Crespo-Hernández, C. E.; Cohen, B.; Kohler, B. Nature. 2005, 436(7054), 1141. doi: 10.1038/nature03933

(3) Middleton, C. T.; De La Harpe, K.; Su, C.; Law, Y. K.;Crespo-Hernández, C. E.; B. Kohler. Ann. Rev. Phys. Chem. 2009,60, 217. doi: 10.1146/annurev.physchem.59.032607.093719

(4) Buchvarov, I.; Wang, Q.; Raytchev, M.; Trifonov, A.; Fiebig, T.Proc. Nat. Acad. Sci. 2007, 104 (12), 4794.doi: 10.1073/pnas.0606757104

(5) Markovitsi, D.; Onidas, D.; Gustavsson, T.; Talbot, F.; Lazzarotto, E.J. Am. Chem. Soc. 2005, 127 (49), 17130. doi: 10.1021/ja054955z

(6) Kuramochi, H.; Kobayashi, T.; Suzuki, T.; Ichimura, T. J. Phys.Chem. A 2010, 114 (26), 8782-8789. doi: 10.1021/jp102067t

(7) Pollum, M.; Crespo-Hernández, C. E.; J. Phys. Chem. 2014, 140,07110. doi: 10.1063/1.4866447

(8) Pollum, M.; Jockusch, S.; Crespo-Hernández, C. E. J. Am. Chem.Soc. 2014, 136 (52), 17930. doi: 10.1021/ja510611j

(9) Taras-Go?lińska, K.; Burdziński, G.; Wenska, G. J. Photochem.Photobiol. A 2014, 275, 89. doi: 10.1016/j.jphotochem.2013.11.003

(10) Harada, Y.; Suzuki, T.; Ichimura, T.; Xu, Y. J. Phys. Chem. B 2007,111 (19), 5518. doi: 10.1021/jp0678094

(11) Harada, Y.; Okabe, C.; Kobayashi, T.; Suzuki, T.; Ichimura, T.;Nishi, N.; Xu, Y. Z. J. Phys. Chem. Lett. 2009, 1 (2), 480.doi: 10.1021/jz900276x

(12) Reichardt, C.; Crespo-Hernández, C. E. J. Phys. Chem. Lett. 2010, 1(15), 2239. doi: 10.1021/jz100729w

(13) Reichardt, C.; Crespo-Hernández, C. E. Chem. Commun. 2010, 46(32), 5963. doi: 10.1039/C0CC01181A

(14) Zhang, Y.; Zhu, X.; Smith, J.; Haygood, M.; Gao, R. J. Phys. Chem.B 2011, 115 (8), 1889. doi: 10.1021/jp109590t

(15) Reichardt, C.; Guo, C.; Crespo-Hernández, C. E; J. Phys. Chem. B 2011, 115 (12), 3263-3270. doi: 10.1021/jp112018u

(16) Martínez-Fernández, L.; González, L.; Corral, I. Chem. Commun.2012, 48 (15), 2134. doi: 10.1039/C2CC15775F

(17) Cui, G.; Fang, W. J. Chem. Phys. 2013, 138 (4), 044315.doi: 10.1063/1.4776261

(18) Pollum, M,;Martínez-Fernández, L.; Crespo-Hernández, C. E.Photoinduced Phenomena in Nucleic Acid I; Springer: NY, US,1970; p 33. doi: 10.1007/128_2014_554

(19) Gobbo, J. P.; Borin, A. C.; Serrano-Andrés, L. J. Phys. Chem. B.2011, 115 (19), 6243. doi: 10.1021/jp200297z

(20) Kobayashi, T.; Kuramochi, H.; Harada, Y.; Suzuki, T.; Ichimura T.J. Phys. Chem. A 2009, 113 (44), 12088. doi: 10.1021/jp905433s

(21) Kobayashi, T.; Harada, Y.; Suzuki, T.; Ichimura, T. J. Phys. Chem. A.2008, 112 (51), 13308. doi: 10.1021/jp803096j

(22) Jiang, J.; Zhang, T. S.; Xue, J. D.; Zheng, X. M. J. Chem. Phys.2015, 143 (17), 11B605_1. doi: 10.1063/1.4935047

(23) Xie, B. B.; Wang, Q.; Guo, W. W.; Cui, G. L. Phys. Chem. Chem.Phys. 2017, 19 (11), 7689. doi: 10.1039/C7CP00478H

(24) Li, M. J.; Liu, M. X.; Zheng, X. J. Acta Phys.?Chim. Sin. 2013, 29(5), 903. [李明娟, 劉明霞, 鄭旭明. 物理化學學報, 2013, 29 (5),903.] doi: 10.3866/PKU.WHXB201302272

(25) Frisch, M. J.; Treucks, G. W.; Schlegel, H. B.; et al. Gauss 09[M];Gaussian Inc.: Wallingford, CT, 2009.

(26) Jamróz, M. H. Spectrochim. Acta A 2013, 114, 220.doi: 10.1016/j.saa.2013.05.096

(27) Dennington, R.; Keith, T.; Millam, J. GaussView, version 5.Semichem Inc.: Shawnee Mission, KS, 2009.

(28) Galica, G. E.; Johnson, B. R.; Kinsey, J. L.; Hale, M. O. J. Phys.Chem. 1991, 95, 7994. doi: 10.1021/j100174a003

(29) Phillips, D. L.; Myers, A. B. J. Chem. Phys. 1991, 95, 226.doi: 10.1063/1.461479

(30) Biswas, N.; Umapathy, S. J. Phys. 1997, 48 (4), 937.doi: 10.1007/BF02845597

(31) Tang, J.; Albrecht, A. C. J. Chem. Phys. 1968, 49 (3), 1144.doi: 10.1063/1.1670202

(32) Tang, J, Albrecht, A. C. Raman Spectrosc. 1970, 33.doi: 10.1007/978-1-4684-3027-1_2

(33) Li, M. J.; Liu, M. X.; Zhao, Y. Y.; Pei, K. M.; Wang, H. G. J. Phys.Chem. B 2013, 117 (39), 11660. doi: 10.1021/jp403798d

(34) Fang, W. X.; Zheng, X. M.; Wang, H. G.; Zhao, Y. Y.; Guang, X.G.; Phillips, D. L.; Chen, X. B.; Fang, W. H. J. Phys. Chem. 2008,133, 134507. doi: 10.1021/jp510396y

(35) Yang, Y.; Pan, S.; Xue, J. D.; Zheng, X. D.; Phillips, D. L.; Fang, W.H. J. Raman Spectrosc. 2014, 45 (1), 105. doi: 10.1002/jrs.4420

(36) Liu, M X.; Xie, B. B.; Li, M. J.; Zhao, Y. Y.; Pei, K. M.;Wang, H. G.; Zheng, X. M. J. Raman Spectrosc. 2013, 44 (3),440. doi: 10.1002/jrs.4213

UV Absorption and Resonance Raman Spectra of 2,4-Dithiouracil

JIN Ying-Chun ZHENG Xu-Ming＊
(Department of Chemistry, Zhejiang Sci-Tech University, Hangzhou 310018, P. R. China)

2,4-Dithiouracil is potentially an important photosensitizer for use in photodynamic therapy. Its photophysics when populated in the lowest excited state has been studied extensively. However, its higher light absorbing excited states and the corresponding reaction dynamics have not been investigated sufficiently. Herein, the resonance Raman spectroscopy and density functional theory were adopted to clarify the electronic transitions associated with the UV absorptions in the far-UV region and the short-time structural dynamics corresponding to the higher light absorbing excited states. The UV absorption spectrum in acetonitrile was deconvoluted into four bands∶ the moderate intense absorption band at 358 nm (f = 0.0336) (A band), the intense broad absorption bands at 338 nm (f = 0.1491), 301 nm (f = 0.1795),and 278 nm (f = 0.3532) (B, C, and D bands) respectively, on the basis of the relationship between the resonance Raman intensities and the oscillator strength f. The result was consistent with the predictions made using the time-dependent density functional theory calculations and the resonance Raman intensitypatterns. Thus, the four bands resulted from the deconvolution are assigned as the S0→ S2, S0→ S6,S0→ S7and S0→ S8transitions, respectively. The resonance Raman spectra of the corresponding B, C,and D bands are assigned and the qualitative short-time structural dynamics are obtained. The major character in the short-time structural dynamics of 2,4-dithiouracil in the S8excited state is that a non-adiabatic process via S8(ππ＊)/S(nπ＊) curve-crossing, accompanied with ultrafast structural distortion,takes place in or near the Franck-Condon region, while the major character in the short-time structural dynamics in the S7and S6excited state appears in the multidimensional reaction coordinates, which are mostly along the C5C6/C2S8/C4S10/N2C3bond lengths + C4N3H9/N1C2N3/C2N1C6/C6N1H7/C5C6H12bond angles for the S7excited state and the C5C6/N3C2/C4S10/C2S8bond lengths + C6N1H7/C5C6H12/C5C6N1/C5C6H12/C2N1C6/ N1C2N3/C4N3H9/N1C2N3bond angles for the S6excited state.

2,4-Dithiouracil; Excited state structural dynamics; UV absorption spectrum;Resonance Raman spectrum; Density functional calculation

April 11, 2017; Revised: May 8, 2017; Published online: May 17, 2017.

O641；O643

10.3866/PKU.WHXB201705175 www.whxb.pku.edu.cn

＊Corresponding author. Email: zxm@zstu.edu.cn; Tel: +86-571-86843627.

The project was supported by the National Natural Science Foundation of China (21473163) and National Key Basic Research Program of China (973)(2013CB834607).

國家自然科學基金(21473163)和國家重點基礎研究發(fā)展規(guī)劃項目(973) (2013CB834607)資助

? Editorial office of Acta Physico-Chimica Sinica