DFT-ZigとQuantum ESPRESSO (QE) のPAW計算比較。silicon-pawの応力・力・バンドを比較
DFT-Zig の PAW 計算結果を Quantum ESPRESSO と比較する。
external/qe/conda_env/bin/pw.xDYLD_LIBRARY_PATH=external/qe/conda_env/lib が必要な場合ありexternal/qe/qe-7.4/build/bin/pw.x は LAPACK (ZHEGV) エラーあり — 使用不可.venv/ (numpy, matplotlib)just buildリグレッションテスト（DFT-Zig実行 + QEベースラインとの比較）:
just test-paw
QEを直接実行して比較する場合:
cd benchmarks/silicon_paw && source ../../.venv/bin/activate
QE=../../external/qe/conda_env/bin/pw.x
$QE < qe_scf.in > qe_scf.out 2>&1
DYLD_LIBRARY_PATH=../../external/qe/conda_env/lib:$DYLD_LIBRARY_PATH $QE < qe_band.in > qe_band.out 2>&1
QE入力ファイル: benchmarks/silicon_paw/
qe_scf.in — SCF (ecut=44Ry, 4x4x4)qe_band.in — バンドqe_force_test.in — 力テスト応力比較: benchmarks/silicon_paw/qe_stress/
| 項目 | QE出力 | DFT-Zig出力 |
|---|---|---|
| バンドギャップ | grep "highest occupied" | band_energies.csv |
| 応力 | grep "total stress" | scf.log |
| 力 | grep "Forces acting" | scf.log |
| エネルギー | 絶対値はG=0規約差あり。差分で比較 |
$ARGUMENTS