Agent-almanac formulate-quantum-problem

install
source · Clone the upstream repo
git clone https://github.com/pjt222/agent-almanac
Claude Code · Install into ~/.claude/skills/
T=$(mktemp -d) && git clone --depth=1 https://github.com/pjt222/agent-almanac "$T" && mkdir -p ~/.claude/skills && cp -r "$T/i18n/ja/skills/formulate-quantum-problem" ~/.claude/skills/pjt222-agent-almanac-formulate-quantum-problem-253f1e && rm -rf "$T"
manifest: i18n/ja/skills/formulate-quantum-problem/SKILL.md
source content

量子問題の定式化

関連する自由度を特定し、ハミルトニアンと状態空間を構築し、境界条件を指定し、適切な近似法を選択し、既知の極限に対して定式化を検証することで、物理系を適切に定式化された量子力学の問題に変換する。

使用タイミング

  • 量子力学の問題を解析的または数値的な解法のために設定する時
  • 量子化学計算(分子軌道、電子構造)を定式化する時
  • 物理的シナリオをディラックまたはシュレーディンガー形式に変換する時
  • 摂動論、変分法、DFT、または厳密対角化の中から選択する時
  • 実験的な分光データまたは散乱データとの比較のために理論モデルを準備する時

入力

  • 必須: 物理系の説明(原子、分子、固体、場など)
  • 必須: 関心のある観測量(エネルギースペクトル、遷移率、基底状態の性質)
  • 任意: 一致させるべき実験的制約またはデータ(スペクトル線、束縛エネルギー)
  • 任意: 望ましい精度レベルまたは計算予算
  • 任意: 好みの形式(波動力学、行列力学、第二量子化、経路積分)

手順

ステップ1: 物理系と関連する自由度の特定

方程式を書く前に系を完全に特性化する:

  1. 粒子内容: すべての粒子(電子、原子核、光子、フォノン)とその量子数(スピン、電荷、質量)をリストアップする。
  2. 対称性: 空間対称性(球、円筒、並進、結晶群)、内部対称性(スピン回転、ゲージ)、離散対称性(パリティ、時間反転)を特定する。
  3. エネルギースケール: 関連するエネルギースケールを決定し、どの自由度が活性でどれを凍結または断熱的に扱えるかを判断する。
  4. 自由度の縮約: 原子核と電子のタイムスケールが分離する場合、ボルン-オッペンハイマー近似を適用する。多体簡約化が適用可能な場合、集団座標を特定する。
## System Characterization
- **Particles**: [list with quantum numbers]
- **Active degrees of freedom**: [coordinates, spins, fields]
- **Frozen degrees of freedom**: [and justification for freezing]
- **Symmetry group**: [continuous and discrete]
- **Energy scale hierarchy**: [e.g., electronic >> vibrational >> rotational]

期待結果: 粒子、量子数、対称性の完全な目録と、活性自由度と凍結自由度の根拠のある選択。

失敗時: エネルギースケールの階層が不明な場合、最初にすべての自由度を保持し、スケール解析の必要性をフラグする。早期の切り捨ては定性的に誤った物理につながる。

ステップ2: ハミルトニアンと状態空間の構築

ステップ1で特定された自由度から数学的枠組みを構築する:

  1. ヒルベルト空間: 状態空間を定義する。有限次元系では基底を指定する(例えばスピン1/2基底|up>、|down>)。無限次元系では関数空間を指定する(例えば3Dの単一粒子にL2(R^3))。
  2. 運動項: 各粒子の運動エネルギー演算子を書く。位置表示ではT = -hbar^2/(2m) nabla^2。
  3. ポテンシャル項: すべての相互作用ポテンシャル(クーロン、調和、スピン-軌道、外部場)を書く。関数形式と結合定数を明示する。
  4. 合成ハミルトニアン: H = T + Vを組み立て、相互作用の種類ごとにグループ化する。多粒子系では交換項と相関項を含めるか、近似を通じてどこで導入されるかを注記する。
  5. 演算子代数: ハミルトニアンがエルミートであることを検証する。問題をブロック対角化するために使用できる運動の定数([H, O] = 0)を特定する。
## Hamiltonian Structure
- **Hilbert space**: [definition and basis]
- **H = T + V decomposition**:
  - T = [kinetic terms]
  - V = [potential terms, grouped by type]
- **Constants of motion**: [operators commuting with H]
- **Symmetry-adapted basis**: [if block diagonalization is possible]

期待結果: すべての項が明示的に記述された完全なエルミートハミルトニアン、定義されたヒルベルト空間、特定された運動の定数。

失敗時: ハミルトニアンが明らかにエルミートでない場合、欠落した共役項やゲージ依存の位相を確認する。ヒルベルト空間が曖昧な場合(例えば相対論的粒子の場合)、形式を明示的に指定し問題点を注記する。

ステップ3: 境界条件と初期条件の指定

一意の解を持つように問題を制約する:

  1. 境界条件: 束縛状態の問題では規格化可能性を要求する(無限遠でpsi -> 0)。散乱問題では入射波の境界条件を指定する。周期系ではブロッホ条件またはボルン-フォン・カルマン条件を適用する。
  2. 領域の制限: 空間領域を指定する。箱の中の粒子では壁を定義する。水素原子では動径と角度の領域を定義する。格子モデルでは格子とそのトポロジーを定義する。
  3. 初期状態(時間依存問題): t=0での状態をエネルギー固有基底での展開として、または指定された中心と幅を持つ波束として定義する。
  4. 制約方程式: 不可分粒子では対称化(ボソン)または反対称化(フェルミオン)を強制する。ゲージ理論ではゲージ固定条件を課す。
## Boundary and Initial Conditions
- **Spatial domain**: [definition]
- **Boundary type**: [Dirichlet / Neumann / periodic / scattering]
- **Normalization**: [condition]
- **Particle statistics**: [bosonic / fermionic / distinguishable]
- **Initial state** (if time-dependent): [specification]

期待結果: 物理的に動機付けされ、ハミルトニアンの定義域と数学的に整合し、一意の解(または明確に定義された散乱行列)を決定するのに十分な境界条件。

失敗時: 境界条件が過剰決定または不足決定の場合、選択された定義域上でのハミルトニアンの自己共役性を確認する。非自己共役なハミルトニアンでは不足指数の慎重な取り扱いが必要である。

ステップ4: 近似法の選択

問題の構造に適した解法戦略を選択する:

  1. 厳密可解性の評価: 問題が既知の厳密に解ける模型(調和振動子、水素原子、イジング模型など)に簡約されるか確認する。該当する場合、厳密解を主要な結果として使用し、補正に摂動論を用いる。

  2. 摂動論(弱結合):

    • H = H0 + lambda Vに分割する(H0は厳密に解ける)
    • lambda VがH0のレベル間隔に比べて小さいことを検証する
    • 縮退を確認する; 必要に応じて縮退摂動論を使用する
    • 適する場合: 相互作用が弱い、少数体系、解析的結果が必要
  3. 変分法(基底状態に焦点):

    • 調整可能なパラメータを持つ試行波動関数を選択する
    • 試行関数が境界条件と対称性を満たすことを確認する
    • 適する場合: 基底状態エネルギーが主要な目標、多体系
  4. 密度汎関数理論(多電子系):

    • 交換相関汎関数を選択する(LDA、GGA、ハイブリッド)
    • 基底セットを定義する(平面波、ガウス型、数値原子軌道)
    • 適する場合: 多電子系、基底状態の密度とエネルギーが必要
  5. 数値的厳密法(小さな系、ベンチマーク):

    • 小さなヒルベルト空間に対する厳密対角化
    • 基底状態サンプリングのための量子モンテカルロ
    • 1次元または擬1次元系に対するDMRG
    • 適する場合: 高精度が必要で系が十分小さい
## Approximation Method Selection
- **Method chosen**: [name]
- **Justification**: [why this method fits the problem structure]
- **Expected accuracy**: [order of perturbation, variational bound quality, DFT functional accuracy]
- **Computational cost**: [scaling with system size]
- **Alternatives considered**: [and why they were rejected]

期待結果: 期待される精度と計算コストの明確な記述、および検討した代替案の文書化を伴う、根拠のある近似法の選択。

失敗時: 単一の方法が明確に適切でない場合、2つの方法で問題を定式化し結果を比較する。方法間の不一致は問題の困難さを明らかにし、さらなる改良を導く。

ステップ5: 既知の極限に対する定式化の検証

解く前に、定式化が既知の物理を再現することを検証する:

  1. 古典極限: hbar -> 0(または大きな量子数)を取り、ハミルトニアンが正しい古典力学に簡約されることを検証する。
  2. 非相互作用極限: 結合定数をゼロに設定し、解が単一粒子状態の直積であることを検証する。
  3. 対称性の極限: 定式化が特定されたすべての対称性を尊重していることを検証する。ハミルトニアンが対称群の下で正しく変換されることを確認する。
  4. 次元解析: ハミルトニアンのすべての項がエネルギーの単位を持つことを検証する。特性長さ、エネルギー、時間スケールが物理的に妥当であることを確認する。
  5. 既知の厳密結果: 系が特殊な場合に既知の厳密解を持つ場合(例えばZ=1の水素原子、二次ポテンシャルの調和振動子)、定式化がそれらを再現することを検証する。
## Validation Checks
| Check | Expected Result | Status |
|-------|----------------|--------|
| Classical limit (hbar -> 0) | [classical Hamiltonian] | [Pass/Fail] |
| Non-interacting limit | [product states] | [Pass/Fail] |
| Symmetry transformation | [correct representation] | [Pass/Fail] |
| Dimensional analysis | [all terms in energy units] | [Pass/Fail] |
| Known exact case | [reproduced result] | [Pass/Fail] |

期待結果: すべての検証チェックが合格すること。定式化が自己整合的で解法の準備ができていること。

失敗時: 検証チェックの失敗は、ハミルトニアンの構築または境界条件の誤りを示す。特定の項または条件まで失敗を遡り、解法に進む前に修正する。

バリデーション

  • すべての粒子と量子数が明示的にリストされた
  • 明確な基底を持つヒルベルト空間が定義された
  • ハミルトニアンがエルミートですべての項が正しい単位を持つ
  • 運動の定数が特定され簡約化に使用された
  • 境界条件が物理的に動機付けされ数学的に十分である
  • 粒子統計(ボソン的/フェルミオン的)が正しく強制された
  • 近似法の選択が期待精度の記述とともに根拠が示された
  • 古典極限、非相互作用極限、対称性の極限が確認された
  • 特殊な場合に既知の厳密結果が再現された
  • 定式化が他の研究者が実装できるほど完全である

よくある落とし穴

  • 自由度の早期省略: エネルギースケールの階層を確認せずに自由度を凍結すると、定性的に重要な物理を見逃す可能性がある。常にエネルギースケールの議論ですべての簡約を正当化する。
  • 非エルミートなハミルトニアン: スピン-軌道結合や複素ポテンシャルでの共役項の忘却。常にH = H-daggerを明示的に検証する。
  • 散乱の境界条件の誤り: 散乱問題に束縛状態の境界条件(規格化可能性)を使用すると、連続スペクトルが完全に失われる。物理的な問題に境界条件を一致させる。
  • 摂動論での縮退の無視: 縮退レベルに非縮退摂動論を適用すると発散する補正が生じる。展開の前に常に縮退を確認する。
  • 単一の近似法への過度の依存: 異なる方法は相補的な失敗モードを持つ。変分法は上界を与えるが励起状態を見逃す可能性がある。摂動論は強結合で発散する。可能な場合はクロスバリデーションを行う。
  • 次元の不整合: 自然単位(hbar = 1)とSI単位を同じ表現で混在させること。最初に一貫した単位系を採用し明示的に述べる。

関連スキル

  • derive-theoretical-result
    -- 定式化された問題から解析的結果を導出する
  • survey-theoretical-literature
    -- 類似の量子系に関する先行研究を見つける