【DFT+實驗】張林娟團隊Nature子刊:雙配體空穴作為LiNiO2的催化活性位點 2024年3月22日 上午10:36 ? 計算 ? 閱讀 143 研究背景 氧釋放反應(yīng)(OER)的電催化作用是許多能量轉(zhuǎn)換裝置的核心,如水和二氧化碳電解器,以及金屬-空氣電池,因此設(shè)計高效的OER催化劑具有重要意義。陰離子氧化還原允許形成O-O鍵,并提供比傳統(tǒng)金屬位點更高的OER活性。 研究背景 近日,中國科學(xué)院上海應(yīng)用物理研究所的張林娟教授團隊在高氧壓力下成功地制備了3d8 L構(gòu)型的LiNiO2,由于在NiIII氧化物的O 2p軌道上發(fā)生了一個電子去除,因此在OER下實現(xiàn)了一個雙配體空穴3d8 L2。LiNiO2在LiMO2、RMO3中表現(xiàn)出超高效的OER活性。理論計算表明,NiIV導(dǎo)致晶格氧之間的直接O-O耦合,*O中間體加速OER活性。這些發(fā)現(xiàn)突出了一種設(shè)計晶格氧氧化還原的新方法,在OER過程中產(chǎn)生足夠的配體空穴。 圖文導(dǎo)讀 圖1. 具有正式的3d7填充的NiIII氧化物的單電子去除光譜示意圖。 NiIII氧化物本質(zhì)上具有陰離子-氧化還原活性,正式的三價Ni離子處于高度共價狀態(tài),其中氧2p配體中的空穴發(fā)揮了重要作用。NiIII型氧化物,如RNiO3、RSrNiO4(R=稀土元素)和LiNiO2,都處于負電荷轉(zhuǎn)移能量狀態(tài)(圖1a),導(dǎo)致3d8 L(L表示位于由O 2p軌道形成的分子軌道上的配體空穴)作為主要的基態(tài)構(gòu)型(圖1b)。眾所周知,當(dāng)NiIII被氧化為NiIIV時,這種電荷轉(zhuǎn)移系統(tǒng)中的電子去除應(yīng)該發(fā)生在O 2p帶內(nèi)而不是Ni 3d帶,,得到一個3d8 L2構(gòu)型,即雙O 2p空穴態(tài)(圖1a)。 基于晶格氧氧化機理,作者提出了雙配體空穴態(tài)催化劑具有較高的OER活性。為了驗證預(yù)測,考慮到LiNiO2具有一個邊緣的Ni-O-Ni網(wǎng)絡(luò),作者采用LiNiO2作為雙配體空穴態(tài)的預(yù)催化劑。此外,NiIII很容易被解離和氧化形成NiIV離子。高活性NiIV的穩(wěn)定歸因于LiNiO2的高熵電荷樣玻璃狀態(tài)。LiNiO2的電子結(jié)構(gòu)具有許多近簡并態(tài),由[3d8]八面體組成,具有3d8、3d8 L和3d8 L2三種構(gòu)型。LiNiO2系統(tǒng)可以在這些狀態(tài)之間波動,這樣氧配體上的空穴態(tài)就被自由能的熵項所穩(wěn)定了。 LiNiO2的結(jié)構(gòu)表征 圖2. 所制備的LiNiO2的結(jié)構(gòu)表征。 作者用同步x射線衍射(XRD)以及透射電子顯微鏡(TEM)(圖2)對制備的LiNiO2的晶相和形態(tài)進行了表征。LiNiO2的XRD圖與R-3m LiNiO2的標(biāo)準(zhǔn)圖匹配良好,沒有任何明顯的晶體雜質(zhì)(圖2a)。圖2a中作者通過RietveLd精修得到了制備的LiNiO2的晶體結(jié)構(gòu)參數(shù)證明了樣品的高純度。作者又通過TEM技術(shù)表明了制備的LiNiO2的單晶性質(zhì)(圖2b-d)。LiNiO2的一個典型的選擇區(qū)電子衍射(SAED)模式顯示出清晰的衍射點,可以沿[010]區(qū)軸指向LiNiO2的R-3m結(jié)構(gòu)(圖2b)。LiNiO2薄部分的高分辨率TEM圖像分別顯示了R-3m結(jié)構(gòu)的(003)和(104)的均勻晶格條紋(圖2c)。作者記錄了校正的高角度環(huán)形暗場掃描透射電子顯微鏡圖像以直接可視化Ni原子沿邊緣的排列,與R-3m LiNiO2的(003)晶格平面很匹配(圖2d)。透射電鏡、掃描電子顯微鏡(SEM)和原子力顯微鏡(AFM)表明,所制備的LiNiO2的形貌為納米薄片,橫向尺寸為200?400nm,層厚度~10nm(圖2e)。 作者通過擴展的x射線吸收精細結(jié)構(gòu)(EXAFS)光譜得到過渡金屬離子的局部配位環(huán)境。LiNiO2的Ni局部配位與一氧化鎳和Ni(OH)2相似,在~1.8和~2.8 ?處有兩個優(yōu)勢峰。這分別是歸因于Ni-O和Ni-Ni的散射路徑(圖3a)。對于這些峰,LiNiO2的徑向距離比Ni(OH)2和一氧化鎳的徑向距離更短,由于Ni具有較高的形式氧化態(tài),支持了較短的Ni-O鍵長。 圖3. 所制備的LiNiO2的XAS表征。 從離子圖來看,低自旋NiIII離子在八面體晶體場中采用了(t2g)6(eg)1的結(jié)構(gòu),并具有Jahn-TeLLer活性。因此,可以期望一個[NiO6]八面體的強局部畸變來消除例如軌道的簡并度。這里存在[NiIIIO6]單元的局部畸變,通過擬合Rspace EXAFS譜,LiNiO2的第一配位殼和較大的無序因子證明了這一點(圖3a)。作者提出了兩種畸變機制。一個是傳統(tǒng)的局部Jahn-TeLLer畸變,在[NiIIIO6]八面體中有四個短的Ni-O(~1.920 ?)和2個長的Ni-O(~2.069 ?)鍵。然而,純Jahn-TeLLer畸變并不能解釋LiNiO2的中子對分布函數(shù)和電子結(jié)構(gòu),因此,作者提出了一種高尺寸的不成比例化(SD)機制。其中Ni?Oshort和Ni?OLong的數(shù)相等。基于這兩種可能的結(jié)構(gòu)模型,作者利用兩個Ni?O配位殼層進行了EXAFS擬合。與僅使用一個Ni-O散射路徑相比,兩種模型都給出了類似但更好的EXAFS數(shù)據(jù)擬合,表明在原始LiNiO2中存在扭曲的[NiIIIO6]。 圖3b中作者用x射線吸收近邊緣結(jié)構(gòu)光譜法測定了所制備的LiNiO2的氧化態(tài)。在Ni K邊緣,LiNiO2的能量位置位于Ni2+和NiIV之間,表明LiNiO2的Ni處于Ni2+和NiIV之間的中間氧化態(tài)。 通過在L3邊緣的電子總產(chǎn)生軟x射線吸收光譜(TEY-sXAS)可以得到表面區(qū)域的詳細電子結(jié)構(gòu)。三維元素的多重光譜特征和L3邊緣的能量位置是高度敏感的價態(tài)和局部環(huán)境。圖3c顯示了LiNiO2的Ni L3邊緣sXAS,以及一氧化鎳和γ-NiOOH的邊緣分別作為Ni2+和NiIII參考。Ni-L3邊緣的多重光譜特征和能量位置與γ-NiOOH非常相似,說明LiNiO2的Ni處于正式的+3氧化態(tài)。γ-NiOOH和LiNiO2的主峰均位于~855 eV,遠高于一氧化鎳的~853 eV,進一步觀察發(fā)現(xiàn),LiNiO2的低能峰定位的比圖3c中一氧化鎳的優(yōu)勢峰高~0.3eV。低能峰是NiIII的固有光譜特征。 作者用O-K sXAS進一步研究了Ni的價態(tài)和Ni-O的共價性。在533 eV以下的前邊緣峰可以被分配到未被占用的3d,如與O 2p軌道雜化。三維元素的價態(tài)的增加與前邊向較低能量的轉(zhuǎn)移以及由于強共價而引起的光譜強度的增強有關(guān)。LiNiO2的O-K sXAS在528.2 eV處有一個強峰,而一氧化鎳在531.9 eV處沒有(圖3d)。因此,在Ni L3邊緣和OK邊緣的TEY-sXAS表明,LiNiO2的表面Ni氧化態(tài)為NiIII,不含Ni2+雜質(zhì)。 電化學(xué)測試 作者利用旋轉(zhuǎn)環(huán)圓盤電極研究了LiNiO2的OER活性。在整個測量過程中,鉑環(huán)的電位保持在0.42 V,這樣一旦在盤上產(chǎn)生的O2可以被收集并在環(huán)上電化學(xué)還原(圖4a)。這樣就可以確定OER的起始電位了。作者采用線性掃描伏安法(LSV)在飽和和無鐵1M氫氧化鉀中估計LiNiO2的OER活性,并與基準(zhǔn)IrO2催化劑和NiOOH進行了比較。從圖4b的結(jié)果上來看,LiNiO2的OER幾何活性最好。此外,LiNiO2具有類似于IrO2的Tafel斜率,并與傳統(tǒng)的NiOOH(圖4c)有所不同,表明LiNiO2具有較快的OER動力學(xué)。與LiMO2、RMO3和其他一元三維過渡金屬催化劑的過電位值相比,想要達到相同的電流密度,LiNiO2需要更少的過電位(圖4d)。 圖4. LiNiO2的電化學(xué)測量和表面化學(xué)。 從RRDE的環(huán)電流來看LiNiO2、IrO2和NiOOH的OER起始電位分別為~1.437V、~1.438V和~1.495V(圖4b),表明LiNiO2在比NiOOH低得多的過電位下就可以激活水的氧化。因此,從過電位、TafeL斜率和起始電位的角度來看,LiNiO2對OER的電催化性能與最先進的IrO2一樣好。在10 mA cm?2條件下,作者用時電位法測定LiNiO2的穩(wěn)定性(圖4e)。在前18h過電位約為10mV,LiNiO2在很大程度上保持了其初始電位,然后逐漸增加,在隨后的30小時內(nèi)保持不變。LiNiO2的高穩(wěn)定性也表明,在48 h穩(wěn)定性試驗前后,LSV的活性惡化程度很?。▓D4e的插圖)。 Ni2+-基氧化物向NiOOH的相互轉(zhuǎn)化通常通過1.2-1.4V電位窗口中的Ni2+/NiIII氧化還原峰來表征。因此,作者采用循環(huán)伏安法研究了LiNiO2生成的NiOOH。如圖4f所示,電位下限從1.41 V逐漸降低到1.11 V,接近于以往研究中通常使用的值。當(dāng)電壓下限設(shè)置為1.41 V時,OER電流是最優(yōu)的,并在多個CV循環(huán)中保持穩(wěn)定。然而,降低到~1.21和~1.11 V的極限會顯著惡化OER電流,伴隨著在-1.36V和~1.26V處的氧化還原峰的出現(xiàn)并逐漸增加,表明表面形成氧化還原活性,但OER相是非活性的。作者使用拉曼光譜對發(fā)生了Ni2+/NiIII氧化還原的LiNiO2樣品在1.56 V下進行了表征(圖4g)。其拉曼光譜在拉曼波段的位置和相對強度上有明顯的差異,而整體的拉曼散射強度也有顯著的增加。與原始的LiNiO2相比,具有相同的電位。所有這些光譜特征都可以歸于β-NiOOH,表明LiNiO2的NiOOH物種經(jīng)過Ni2+/NiIII氧化還原循環(huán)后的結(jié)構(gòu)轉(zhuǎn)變。在未純化的1M氫氧化鉀中也進行了相同的伏安法測試,盡管在相近電位范圍內(nèi)出現(xiàn)了Ni2+/NiIII氧化還原峰,但OER電流這表明Fe雜質(zhì)顯著提高了NiOOH的性能。 綜上循環(huán)伏安分析表明,LiNiO2的OER活性和電化學(xué)性質(zhì)受到電解液中電勢下限和鐵雜質(zhì)存在的影響。在進行Ni2+/NiIII氧化還原時,LiNiO2的表面結(jié)構(gòu)轉(zhuǎn)化并降解為NiOOH,在沒有Fe的情況下,其OER活性較差。相比之下,原始和結(jié)構(gòu)未改變的LiNiO2,表現(xiàn)出更好的內(nèi)在OER活性。 原位-XRD,XAS和拉曼研究 從上述結(jié)構(gòu)和電化學(xué)測量中,LiNiO2的原始性質(zhì)已從結(jié)晶相、局部配位環(huán)境、氧化狀態(tài)等方面得到了表征,證實了LiNiO2的OER活性優(yōu)于類似的NiOOH。然而,LiNiO2的活性物種卻不同于NiOOH。此外,鋰離子電池典型的陰極材料LiNiO2可以脫氫,增加Ni氧化態(tài),誘導(dǎo)結(jié)構(gòu)轉(zhuǎn)變,因此,在氫氧化鉀水溶液中的OER期間,剝離將變得類似,甚至更容易。為了建立LiNiO2的結(jié)構(gòu)與其OER活性之間的關(guān)系,作者進行了基于同步加速器的原位XRD和XAS,以及拉曼光譜實驗。 圖5. OER條件下的LiNiO2脫鋰。 圖5a顯示了LiNiO2在1M氫氧化鉀中在1.41-1.81V的電位范圍內(nèi)的原位同步加速器XRD。這些模式由三個部分組成:位于~12°和18°的電池和材料的非晶態(tài)背景,碳紙石墨碳衍射峰(箭頭所示)以及LiNiO2及其分解衍生物的衍射模式。圖5b顯示了1.41-1.81V電位范圍內(nèi)的XANES譜。隨著外加電位的增加,上升邊緣向更高的能量移動,表明NiIII的氧化與OER的發(fā)生相一致。圖6c顯示LiNiO2的平均Ni氧化態(tài)是應(yīng)用勢的函數(shù)。氧化態(tài)在1.41 V時的氧化態(tài)為~3.0,在1.56 V時急劇增加到~3.3,在電位>1.56 V時逐漸增加,達到~3.4的平臺。此外,前邊緣峰值強度隨著應(yīng)用電位的增加而增加(圖5b),對應(yīng)于八面體對稱性的降低。結(jié)合Ni的電位依賴性氧化態(tài)與LiNiO2的電化學(xué)反應(yīng),可以發(fā)現(xiàn)NiIV離子的數(shù)量與OER電流呈正相關(guān)(圖6c),,表明游離LixNiO2中的NiIV可以被認為是OER活性中心的一部分。 圖6. 在OER條件下,LiNiO2的局部相互作用。 作者利用原位EXAFS探測了Ni局部配位環(huán)境與電位的關(guān)系(圖6a,b)。與非原位數(shù)據(jù)相似,Ni配位殼層以Ni-O和Ni-Ni散射路徑為主。在~1.45A時,Ni-O的振幅隨著應(yīng)用電位的增加而增加,從1.41V增加到1.81V(圖6b)。由于原始LiNiO2的Ni第一殼層配位被6個氧原子飽和,Ni-O的配位數(shù)不能進一步增加,這也意味著Ni-O鍵隨著應(yīng)用勢的變化而變得越來越有序。隨著振幅的增加,Ni-O峰邊界的降低也表明了Ni-O畸變的降低(圖6b)。除振幅外,Ni-O的等高線圖在相應(yīng)的較低的R區(qū)域的貢獻逐漸增加(圖6b)。作者通過EXAFS擬合確定了Ni-O鍵長從1.41 V時的1.955 ?收縮到1.81 V時的1.909? (圖6c),證實了在OER條件下NiIV的存在。綜上所述,在OER過程中發(fā)現(xiàn)了LixNiO2中的NiIV和氧空位,這歸因于高OER活性。 作者進一步利用原位拉曼光譜追蹤了在大于1.4 V的高電位下LiNiO2的中間種。與上述圖4g中的拉曼實驗不同,LiNiO2現(xiàn)在沒有被迫進行Ni2+/NiIII氧化還原,但是直接被施加1.46V的電壓。圖6d中LiNiO2拉曼的光譜輪廓僅在OER過程中,隨著應(yīng)用電位的增加,存在輕微的紅移帶,與脫氫LiNiO2的形成一致。這支持了上述原位數(shù)XRD和XAS的結(jié)果。通過對圖6d和圖4g的比較,作者推斷LiNiO2的晶體結(jié)構(gòu)、中間體和OER活性很大程度上依賴于所應(yīng)用的最低電位。 γ-NiOOH的存在和形成被認為是許多鎳基氧化物中OER活性的起源。然而,根據(jù)原位拉曼光譜(圖6d),LiNiO2的活性物種不是γ-NiOOH。此外,在原位數(shù)據(jù)中沒有觀察到LiNiO2(4.72 ?)向γ-NiOOH(≥7 ?)的轉(zhuǎn)換導(dǎo)致[NiO2]層間間距的大幅擴展。 DFT計算和OER機制的DEMS測試 圖7. LiNiO2的OER機理。 作者采用第一性原理DFT計算,研究了脫鋰LiNiO2增強OER活性的潛在機理。對涉及四個質(zhì)子-電子轉(zhuǎn)移步驟的OER機制,作者考慮了三種機制:金屬位點吸附質(zhì)演化(MAE)機制,晶格氧空位點(LOV)機制以及金屬和晶格氧空位位(MLOV)機制。在這些機制中,反應(yīng)步驟的吸附位點是不同的:MAE中只有Ni位點,LOV中是氧位點, MLOV中是兩個位點(圖7a)。作者建立了完美LiNiO2和脫鋰LixNiO2模型,研究OER過程中產(chǎn)生的不飽和NiIV狀態(tài)的影響。 對Li0.5NiO2和LiNiO2,通過三種OER機制得到的每個OER步的吉布斯自由能差(ΔG)分別如圖7b所示。首先,無論OER途徑如何,Li0.5NiO2的電位限制步驟都比LiNiO2具有更低的過電位,表明LixNiO2為OER的實際活性相。其次,在Li0.5NiO2的三種OER場景中,MLOV路徑顯示最低的過電位(0.69 V),表明通過MLOV的OER占脫Li0.5NiO2,且電位限制步驟是對*OH的吸附(步驟3)。第三,對于LiNiO2上的OER,MLOV路線也占優(yōu)勢,但潛在限制步驟轉(zhuǎn)變?yōu)?sup>*OO中間體和氧空位的形成(步驟1)。 LiNiO2在脫鋰前后的限制步驟的轉(zhuǎn)變可以歸因于Li0.5NiO2中NiIV離子的存在引起的雙配體空穴態(tài)。在八面體對稱[NiO6]體系的電子結(jié)構(gòu)中,從O 2p到Ni 3d軌道的電荷轉(zhuǎn)移能量隨著Ni離子價態(tài)的增加而減小。LiNiO2體系已經(jīng)處于負電荷轉(zhuǎn)移狀態(tài),因此,Ni-O相互作用是高度共價的,其中NiIII采用接近3d8 L的平均構(gòu)型,在O 2p軌道上留下一個配體空穴。當(dāng)部分NiIII通過LiNiO2脫鋰氧化為NiIV時,共價進一步增加,得到的NiIV被賦予3d8 L2的平均構(gòu)型,此時,O 2p中處于雙配體空穴態(tài)。LixNiO2中具有O 2p配體特征的高未占據(jù)密度狀態(tài)激活晶格氧,有利于與含氧中間體直接耦合并形成氧空位,這是LiNiO2的電勢限制步驟。對OER過程中氧空位形成的預(yù)測與XAS原位測試過程中的實驗觀察結(jié)果很一致(圖6c)。因此,密度泛函理論表明,在OER條件下產(chǎn)生的雙配體空穴很容易激活晶格氧,加速了O-O結(jié)合和氧空位的形成,雙配體空穴態(tài)的存在與脫鋰的LixNiO2的高OER活性密切相關(guān)。 為了證實晶格氧在OER過程中的參與,作者使用原位差分電化學(xué)質(zhì)譜(DEMS)進行了18O同位素標(biāo)記實驗。實驗首先在富含18O的氫氧化鉀電解質(zhì)中的未標(biāo)記的LiNiO2上進行,在1.4 V到1.8 V的電壓范圍內(nèi)進行6圈CV測試。圖7c顯示了第一圈CV和對應(yīng)的質(zhì)譜循環(huán)伏安圖(MSCV)。在MSCV曲線中,16O18O的質(zhì)量信號m/z=34,主要是由LiNiO2晶格的16O與電解質(zhì)的18O耦合引起的,顯示出顯著的離子電流,以36O2為主要分子氧同位素。測量到的16O18O質(zhì)量信號與理論上預(yù)期的18O富含電解質(zhì)的氧化曲線比較(圖7c的插圖),顯示出16O18O離子電荷。這證明了晶格氧在OER過程中O2的形成中起著直接的作用。18O同位素的量隨著循環(huán)次數(shù)的增加而逐漸減少,表明標(biāo)記的LiNiO2消耗了18O。因此,雙向同位素標(biāo)記實驗證實了LiNiO2晶格氧參與了OER,證實了預(yù)測的MLOV機制和雙配體孔穴的活化晶格氧的存在。 Fe3+雜質(zhì)效應(yīng) 人們普遍認為可以將Fe3+摻入Ni(氧)氫氧化物中,通過形成氧橋接的Ni和Fe位點的協(xié)同活性中心來促進OER活性,且其活性可以通過結(jié)構(gòu)排序進一步增強。結(jié)構(gòu)完整的LixNiO2具有穩(wěn)定的雙配體空穴,可以作為一個比通常無定形的NiOOH更好的形成這種協(xié)同中心的平臺。 圖8a比較了在純化的1M氫氧化鉀中,用Fe3+滴定的LiNiO2的線性掃描伏安圖。與報道的Ni(oxy)氫氧化物相似,LiNiO2的OER活性與Fe3+濃度的E和TafeL斜率密切相關(guān)(圖8a-c):(1)E隨著Fe3+濃度的降低而降低,在100μM Fe3+存在下在1.500 V時飽和8a;(2)在只有1μM Fe3+存在的情況下,TafeL斜率急劇下降,在Fe3+濃度為200μM時基本保持不變(圖8b)。在如此低的Fe3+濃度下,TafeL斜率下降,表明OER切換到活性較高的Fe電催化位點,而這種新的活性位點的數(shù)量隨著Fe3+濃度的增加而增加,導(dǎo)致了依賴Fe濃度的E。為了進行比較,作者在含有100μM Fe3+的1M氫氧化鉀中也測量了NiOOH。盡管整體活動如預(yù)期的那樣得到了高度的提升,但NiOOH-Fe3+的E比LiNiO2-Fe3+高36 mV,OER起始電位高25 mV(圖8d),證明了LiNiO2是形成鎳鐵活性位點的良好平臺。此外,在鎳鐵活性位點中高價鎳的存在也促進了電催化的耐久性。 總結(jié)展望 本文中作者通過對電化學(xué)性質(zhì)和原位光譜的綜合對比,揭示了OER的高熵LiNiO2的原始活性、實際活性物種和反應(yīng)中心。當(dāng)防止Ni2+/NiIII的電化學(xué)氧化還原和鐵雜質(zhì)的引入時,LiNiO2獲得了最佳和真實的OER活性,其活性與最先進的IrO2的活性相當(dāng),并且在過電位方面遠遠高于傳統(tǒng)的NiOOH。本研究闡明了鋰化金屬氧化物的電催化作用,并強調(diào)了精細化的電化學(xué)測量和原位技術(shù)的重要性。 文獻信息 HaoLiang Huang, Yu-Chung Chang, Yu-Cheng Huang, LiLi Li, ALexander C. Komarek, Liu Hao Tjeng, Yuki Orikasa, Chih-Wen Pao,Ting-Shan Chan, Jin-Ming Chen, Shu-Chih Haw, Jing Zhou, Yifeng Wang,Hong-Ji Lin, Chien-Te Chen, Chung-Li Dong, Chang-Yang Kuo,Jian-Qiang Wang, Zhiwei Hu & Linjuan Zhang, UnusuaL doubLe Ligand hoLes as cataLytic active sites in LiNiO2.Nature Communications,14, 2112 (2023). 原創(chuàng)文章,作者:計算搬磚工程師,如若轉(zhuǎn)載,請注明來源華算科技,注明出處:http://www.xiubac.cn/index.php/2024/03/22/2fd189d321/ 贊 (0) 0 生成海報 相關(guān)推薦 國家級“四青”人才、Science一作!張雷教授線上報告:氫燃料電池超低Pt催化劑! 2023年10月6日 計算成果,8篇頂刊:Nature子刊、Angew、JACS、EES、AM、AFM等! 2023年10月17日 喬世璋/焦研JACS:理論計算+機器學(xué)習(xí),揭示Mxene上C-N電偶聯(lián)反應(yīng)的活性和選擇性路線 2023年10月30日 6個案例,3個比方,講清倒空間、倒格子、第一布里淵區(qū)、高對稱k點、k path等概念!11月7日發(fā)車! 2023年11月6日 PRX:機器學(xué)習(xí)預(yù)測相干聲子熱傳導(dǎo)最小化 2024年3月17日 【純計算】JACS:用多級量子力學(xué)/分子動力學(xué)模擬探討水溶液中Mg和Ca離子的pH依賴脫水動力學(xué) 2024年1月16日