以下は interatomic_distance.py の解析結果を、ドキュメント形式でまとめたものです。 1) プログラムの動作概要 - 目的 - 与えられた結晶系の格子パラメータとサイト情報に基づき、原子間距離を計算してリスト化する。距離の範囲は rmin 〜 rmax で絞り込み、近傍セル間の距離も考慮して全ての組み合わせを列挙する。 - 処理の流れ - lattice_parameters を基に格子ベクトルを計算（cal_lattice_vectors）。 - 格子のメトリックテンソル（gij など）を計算して表示。 - 容積（体積）を計算して表示。 - Reciprocal lattice（逆格子）ベクトルとそのメトリックテンソルを計算して表示。 - nxmax, nymax, nzmax を求め、近傍セルの最大オフセットを決定。 - sites に定義された原子配列を二重ループで回し、各サイト間の距離を rmin〜rmax の範囲で計算。距離 dis が範囲内であれば、リスト rlist に [dis, site0, site1, ix, iy, iz] の形で格納する。 - rlist を距離で昇順ソートし、各組み合わせの距離と原子配列位置・セルオフセットを表示。 - 出力 - 標準出力（コンソール）へ、格子パラメータ、格子ベクトル、逆格子、メトリックテンソル、体積、近傍セルの最大オフセット、そして距離リストを整形して出力する。 - 使用言語/依存 - Python 3.x - NumPy および Matplotlib（描画用のモジュールとして import あり） - 依存モジュール tkcrystalbase（ローカルの tkcrystalbase.py が必要） 2) 必要な非標準ライブラリとインストール方法 - 必須ライブラリ（非標準とみなされうるもの） - numpy - matplotlib - mpl_toolkits.mplot3d（matplotlib の一部） - tkcrystalbase（外部モジュール、同梱の tkcrystalbase.py が必要） - インストール方法 - Python 環境で以下を実行（例: pip 使用） - pip install numpy matplotlib - tkcrystalbase.py の入手 - tkcrystalbase.py が別ファイルとして入手可能な場合は、同じディレクトリに配置。公式リポジトリや提供元の指示に従って入手してください。 - 必要に応じて仮想環境を作成: - python -m venv venv - source venv/bin/activate (Windows の場合: venv\Scripts\activate) - pip install numpy matplotlib - tkcrystalbase.py を同ディレクトリへ配置 - 注意 - 本コードは tkcrystalbase.py の関数 cal_lattice_vectors、cal_metrics、cal_volume、cal_reciprocal_lattice_vectors、cal_reciprocal_lattice_parameters、distance などを前提としており、これらの関数の実装依存で動作します。tkcrystalbase.py が正しく提供されていることを確認してください。 3) 必要な入力ファイル - 本プログラムはコード内にハードコードされたデータを用いて処理します。入力ファイルを外部から読み込む設計にはなっていません。 - 依存する外部モジュール tkcrystalbase.py が正しく提供されていること、また必要な関数が期待通り動作することを前提とします。 - もし外部からデータを取り入れたい場合は、sites、lattice_parameters、rmin、rmax などの設定をファイル読み込みに変更する必要があります。 4) 実行後に生成される出力ファイル - 本コードはファイルへ出力を行わず、標準出力へテキストを表示します。 - したがって、実行によって生成される「出力ファイル」はありません（必要に応じてリダイレクトして保存可能）。 - 例: python interatomic_distance.py > results.txt などで結果をファイルに保存可能。 - なお、生成される主な情報は以下の通り： - Lattice parameters、格子ベクトル（ax, ay, az） - Reciprocal lattice parameters・格子ベクトル・メトリックテンソル - 実体積（Volume）および Reciprocal unit cell volume - nxmax, nymax, nzmax の値 - Interatomic distances の一覧（距離、サイト0/サイト1、対応する座標、オフセット ix, iy, iz、距離 dis） 5) コマンドラインでの使用例 (Usage) - 基本的な実行方法 - python interatomic_distance.py - 出力をファイルに保存する例 - python interatomic_distance.py > interatomic_distances.txt - 実行環境の例（UNIX系・Windows 共通） - 事前準備 - Ensure Python 3.x がインストールされていること - numpy, matplotlib がインストール済みで tkcrystalbase.py が同一ディレクトリにあること - 実行コマンド - cd <スクリプトのあるディレクトリ> - python interatomic_distance.py - 実行時の出力をスクリーンに表示して、必要に応じファイルへリダイレクト - python interatomic_distance.py | less - python interatomic_distance.py > distances.log 2>&1 補足 - 現在のスクリプトでは距離計算の対象として、Na と Cl のサイトが複数定義された NaCl型の仮想結晶が用意されています。格子は立方晶（abc が 5.62 Å、α=β=γ=90°）として扱われ、rmin 0.1 Å、rmax 4.5 Å の範囲で近傍距離を抽出します。 - 実際の挙動は tkcrystalbase.py の実装に依存します。距離の計算関数 distance(pos0, pos1, gij) や格子ベクトル計算関数 cal_lattice_vectors などが正しく動作する必要があります。 この文書は、コードの機能と実行方法を要点だけを分かりやすくまとめたものです。必要に応じて、使われている関数の仕様（引数・返り値・例）を tkcrystalbase.py のドキュメントと突き合わせて追加すると、より詳しい技術ドキュメントになります。