肖宿點點頭,又問:“自洽場迭代的收斂判據是多少?”
“預設的,10的負8次方。”
肖宿翻到某一頁,手指輕點公式:“這個地方的交換相關泛函,用的是雜化型的?”
萬匯楊眼睛一亮:“對,B3LYP。你怎麼看出來的?”
肖宿並未回答,繼續往下翻閱。
方清在一旁看著,心中不禁有些驚訝。
這些問題問得十分專業,絕非隨意翻閱就能問出來的。
但轉念一想,搞計算化學的人也都會問這些,關鍵還得看能否找到問題的根源。
他忍不住開口:“肖宿同學,其實我們遇到的主要問題,在於處理分子對稱性時。一些高對稱性的構型,計算結果總是偏離實驗值。萬老師嘗試將你的度量方法加進去,卻似乎不太相容。”
這話雖委婉,但潛臺詞很明顯:你的方法我們試過了,沒什麼效果。
肖宿抬起頭,目光平靜地看向方清。
那眼神波瀾不驚,卻讓方清莫名感到一絲壓力。
“你們的實現方式存在問題。”肖宿直言。
方清一愣:“什麼問題?”
肖宿並未直接回應,而是從書包裡掏出筆,在列印稿的空白處畫了幾筆。
“你剛才提到,分子屬於C2v點群?”
方清點頭:“沒錯。空間群是C2v,有四個對稱操作:恆等操作、繞z軸的二重旋轉、以及兩個鏡面反射。我們在構造哈密頓量時,已將這些對稱性考慮在內。”
肖宿在幾筆旁邊寫下一個公式:“你們的模型中,這個係數應當滿足對稱性約束。但在處理奇點附近時,你們破壞了這種約束。”
方清盯著那個公式,眉頭緊鎖。
那個係數他自然知曉,是哈密頓量裡某個關鍵項的權重因子。
按照對稱性要求,在C2v群的作用下,它應保持不變。
可肖宿所畫奇點附近的區域,他仔細思索後發現,似乎在數值離散化時,為處理奇異性,他們做的某種近似處理……
難道真的破壞了對稱性?
萬匯楊湊過來,盯著肖宿畫的圖看了許久,突然一拍大腿。
“我明白了!”
他猛地站起身,在辦公室裡來回踱步,嘴裡唸唸有詞:
“對對對,我們在奇點附近採用了非結構網格,網格劃分時未考慮對稱性約束,致使那個係數在對稱操作下不再相等。所以計算出的能量本徵值才會出現偏差!”
方清愣住了。
。轉運速飛腦大,式公個那看看又,筆幾的畫宿肖向看他
……此如實確乎似,說一麼這師老萬經
。除解誤錯被度併簡的級能些某者或,級能的在存該不本原出裂分會便譜徵本的量頓哈,壞破被稱對若
。關相稱對子分與律規差偏且,符不值驗實與量能的出算計致導將這
——是可
”。束約稱對援支不並生格網其,庫元限有的現是的用們我且況。高提幅大會度雜複的分劃格網,稱對持保格嚴要,格網構結非是域區個那“,道說住不忍清方”,師老萬“
”。格網分劃新重需無“:淡平氣語,他向看宿肖








