01引言

大位阻烯烴分子馬達(dá)可以在光/熱交替控制下實(shí)現(xiàn)單向旋轉(zhuǎn)，從而能有效地將光能轉(zhuǎn)化為機(jī)械能。其中熱翻轉(zhuǎn)的能壘較高，影響了馬達(dá)的旋轉(zhuǎn)效率。全光控的分子馬達(dá)可以規(guī)避熱反應(yīng)路徑，其旋轉(zhuǎn)效率僅依賴于光吸收，相比于傳統(tǒng)馬達(dá)具有明顯優(yōu)勢(shì)，然而探究全光控條件需要借助高精度量子化學(xué)方法計(jì)算其復(fù)雜的激發(fā)態(tài)勢(shì)能面，另外需要通過非絕熱動(dòng)力學(xué)模擬揭示其動(dòng)力學(xué)過程。近年來，基于機(jī)器學(xué)習(xí)的勢(shì)能面應(yīng)用于加速分子動(dòng)力學(xué)模擬，但擬合激發(fā)態(tài)勢(shì)能面仍然充滿挑戰(zhàn)。

02成果簡(jiǎn)介

針對(duì)近期實(shí)驗(yàn)報(bào)道的一種全光控的分子馬達(dá)，我們基于自旋翻轉(zhuǎn)的含時(shí)密度泛函理論（Spin-Flip TDDFT）計(jì)算了基態(tài)/激發(fā)態(tài)勢(shì)能面，從而揭示了其光異構(gòu)化反應(yīng)過程，探究了其實(shí)現(xiàn)全光控旋轉(zhuǎn)的機(jī)理。同時(shí)，我們結(jié)合經(jīng)典分子力場(chǎng)，利用機(jī)器學(xué)習(xí)構(gòu)建了具有高精度的解析勢(shì)能面，用于加速模擬分子馬達(dá)的光致異構(gòu)反應(yīng)，得到了與實(shí)驗(yàn)一致的異構(gòu)化產(chǎn)率。這套方法還可以用于探究其它復(fù)雜體系的光異構(gòu)化反應(yīng)。

03圖文導(dǎo)讀

1、靜態(tài)計(jì)算研究

利用SF-TDDFT方法對(duì)BHHLYP、B3LYP和CAM-B3LYP等三種泛函進(jìn)行標(biāo)定，比較了不同泛函計(jì)算激發(fā)態(tài)勢(shì)能面的差異。如圖1所示，我們發(fā)現(xiàn)SF-BHHLYP方法能夠準(zhǔn)確描述CI區(qū)域的勢(shì)能面形貌。在扭轉(zhuǎn)角達(dá)到-90°時(shí)出現(xiàn)能量兼并區(qū)域，即S0/S1電子態(tài)的強(qiáng)耦合區(qū)域。隨后使用SF-BHHLYP方法計(jì)算二維勢(shì)能剖面，并分析了光異構(gòu)化路徑。如圖2所示，通過對(duì)二維勢(shì)能面的分析，我們發(fā)現(xiàn)了馬達(dá)轉(zhuǎn)動(dòng)的光異構(gòu)化反應(yīng)路徑：[a]→[CI-1]→[b]→[CI-2]→[d]。最后比較了TDDFT和SF-TDDFT兩種方法計(jì)算光譜的差異，并分析了S0→S1的軌道躍遷類型，表明其躍遷類型是π→π*。另外光譜結(jié)果表明SF-BHHLYP方法得到的垂直激發(fā)能與實(shí)驗(yàn)光譜的吸收峰位置符合。

圖1 SF-BHHLYP/cc-pVDZ計(jì)算的一維勢(shì)能面

圖2 SF-BHHLYP計(jì)算的（a）二維絕熱勢(shì)能面和（b）電子態(tài)能量差的等高繪圖

圖3 TDDFT和SF-TDDFT計(jì)算的吸收光譜及其對(duì)應(yīng)的π→π*軌道躍遷

2、非絕熱動(dòng)力學(xué)模擬

利用機(jī)器學(xué)習(xí)模型構(gòu)建基態(tài)/激發(fā)態(tài)勢(shì)能面，如圖4所示，ML模型在測(cè)試集上的MAE均小于2.0 kcal/mol（內(nèi)插圖），預(yù)測(cè)效果接近化學(xué)精度?；赪igner取樣生成初始構(gòu)型和速度，用于Surface Hopping非絕熱動(dòng)力學(xué)模擬，并統(tǒng)計(jì)電子態(tài)布局?jǐn)?shù)和反應(yīng)坐標(biāo)的演化，計(jì)算反應(yīng)速率和產(chǎn)率。如圖5所示，結(jié)果表明對(duì)稱折疊的異構(gòu)體[b’]是異構(gòu)化反應(yīng)的主要產(chǎn)物（48%），而反向折疊異構(gòu)體[b]的轉(zhuǎn)化效率為8%。計(jì)算得到反應(yīng)時(shí)間常數(shù)為235 fs，與以往的二代分子馬達(dá)理論計(jì)算結(jié)果一致。

圖4基態(tài)（a）和激發(fā)態(tài)（b）勢(shì)能面擬合的交叉驗(yàn)證結(jié)果

圖5非絕熱動(dòng)力學(xué)模擬中的幾何參數(shù)演化：（a）烯烴平面的金字塔化（b）扭轉(zhuǎn)角和（c）環(huán)翻轉(zhuǎn)

04小結(jié)

我們使用鴻之微量子化學(xué)計(jì)算軟件BDF，基于SF-TDDFT方法研究了一種新型全光控分子馬達(dá)的異構(gòu)化反應(yīng)，表明其在連續(xù)光照條件下可以實(shí)現(xiàn)單向旋轉(zhuǎn)。另外我們利用機(jī)器學(xué)習(xí)模型構(gòu)建了基態(tài)/激發(fā)態(tài)勢(shì)能面，并對(duì)光異構(gòu)化過程進(jìn)行非絕熱動(dòng)力學(xué)模擬。計(jì)算結(jié)果表明，實(shí)現(xiàn)全光控轉(zhuǎn)動(dòng)的關(guān)鍵因素是四個(gè)異構(gòu)體具有接近的垂直激發(fā)能，這有利于馬達(dá)在相同光照條件下連續(xù)激發(fā)。并且在不穩(wěn)定異構(gòu)體之間存在著光反應(yīng)通道，這有利于馬達(dá)的單向轉(zhuǎn)動(dòng)。另外，定子的環(huán)翻轉(zhuǎn)對(duì)光異構(gòu)化有顯著影響，導(dǎo)致其主要產(chǎn)物為對(duì)稱折疊異構(gòu)體。我們利用機(jī)器學(xué)習(xí)模型開發(fā)了一套高效的非絕熱動(dòng)力學(xué)模擬方法，用于光異構(gòu)化反應(yīng)體系的研究，并給出和實(shí)驗(yàn)結(jié)果符合的反應(yīng)速率和產(chǎn)率。該模型還可以推廣到其他的復(fù)雜光化學(xué)反應(yīng)體系。








審核編輯：劉清