第 3 週工作週報
Overview · 本週總覽
這禮拜做了什麼
4個二元合金完成 slab 收斂
66個 CuPt 吸附位已收斂
2組方法學對照定案
收斂配方推廣四合金
Research
CuPt 配方複製到 CoPt／FePt／NiPt，全數收斂
CuPt 吸附位大批送算
Research
96 位點分批送出，top 16／bridge 48 已全數收斂
方法學對照定案
Research
偶極修正差 1 meV、鎖 3 vs 4 層差 2 meV
研究是主軸（3 頁詳細講）；這週 side project 都按兵不動。
Research 1 · 收斂配方推廣
兩階段 LDIPOL 配方，從 CuPt 複製到整個 Pt-skin 家族
- 階段一：先關
LDIPOL讓結構鬆弛、SCF 穩定 - 階段二：拿
CONTCAR當起點、開LDIPOL、補質心座標 - CoPt 第一段撞牆鐘被中斷，auto-chain 自動接續後收斂
- FePt／NiPt 磁矩差很多，原擔心難收斂，結果一樣順利
學到：沒收斂不等於設定錯 —— 要先分清是「物理沒收斂」還是「被系統掐斷」。驗證過的配方可直接複用到相近的二元合金。
Research 2 · 吸附位大批送算
CuPt 96 個吸附位：分批送出，已收斂 66 個
- 位點共 96 組（top 16／bridge 48／hollow 32），依教授的「逐類型」順序分批送
- 已收斂 66 組：top 16 全數、bridge 48 全數，加 hollow 2 組先導
踩到的坑
- 一次投 76 個 job 撞 **QOS 上限**，十幾個被拒
- 誤把排隊中約 60 個 job 一起 `scancel`
- `bridge_17` 跑兩段仍未收斂
怎麼解
- 改用 `.submitted` 標記的 **refill 腳本**，撞上限就停、隊列空了自動接續
- OSZICAR 的 CG 試探步暴衝假值，要往前找真值
剩 `hollow_3~32` 共 30 組待送；ΔG_H* 目前用 PPT 的暫定 ½E(H₂) 值算，**要自算 H₂ 參考態後才是定值**。其他三個合金（NiPt／CoPt／FePt）的吸附位尚未開始。
Research 3 · 方法學對照
兩組對照實驗：把「要不要多花算力」問清楚
偶極修正開 vs 關
- 同一個 top 位點，`LDIPOL` 開與關各算一次
- 能量淨差只有 **1 meV**
- 結論：production 用**關**的，省算力不失精度
底層鎖 3 層 vs 4 層
- 鎖 3 層（基準）與鎖 4 層對照
- 能量差 **2 meV**
- 結論：**鎖 3 層就夠**，不必多鎖一層
另一組對照（鎖 2 層）因為單一 job 跑超過一天太慢，與老師討論後直接取消，不硬等。
這兩組對照的意義：先用少量計算把參數選擇「釘死」，後面 96 組（乃至四個合金共約 384 組）才不會用錯設定重算。
Research 4 · 初步結果與下一步
66 組能量出來了，但 ΔG_H* 還不是定值
為什麼還不能下結論
- ΔG_H* ＝ E(slab+H) − E(slab) − **½E(H₂)** ＋ 0.24 eV
- 目前 **½E(H₂) 用的是文獻／簡報上的暫用值**，還沒自己算
- 基準一改，**所有位點的 ΔG_H\* 會整批平移**
目前看到的趨勢（暫定）
- 已有 66 組收斂能量可以代入
- 相當一部分位點落在 **|ΔG_H*| < 0.1 eV** 的理想區間
- 趨勢（哪類位點較好）可信；**絕對值待 H₂ 自算後定案**
下一步照順序：①自算 ½E(H₂) 參考態 → ②送 `hollow_3~32` 補滿 96 組 → ③重投未收斂的 `bridge_17` → ④CuPt 收齊後，才推到 NiPt／CoPt／FePt。
Build · 週報自動化管線
「收工 → 彙整 → 發佈」這條線繼續補防呆
掃描來源擴充
修正 EXCLUDE 規則
發佈時一併清草稿
專題日誌自動彙整
project-roots.mjs掃描來源加入 codex 工作區，之後新專案會自動被收工／週報涵蓋- 一度誤把桌面「網站更新建議」資料夾整個排除，隔天發現後改回每日自動更新
publish-weekly.mjs：發佈週報時一併把草稿檔刪除，不留殘留- 每天自動把「專題日誌」寫成部落格文章（07/21–07/25 共 4 篇），這份週報也是拿它們拼出研究進度
這些都是修 bug／補防呆，讓自動化這條線更穩，不是新功能。