From 8f8c59965c7924d64330f72157a59f01759c5a9d Mon Sep 17 00:00:00 2001
From: "pre-commit-ci[bot]"
 <66853113+pre-commit-ci[bot]@users.noreply.github.com>
Date: Tue, 8 Oct 2024 01:29:51 +0000
Subject: [PATCH 1/9] [pre-commit.ci] pre-commit autoupdate
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

updates:
- [github.com/pre-commit/pre-commit-hooks: v4.6.0 → v5.0.0](https://github.com/pre-commit/pre-commit-hooks/compare/v4.6.0...v5.0.0)
- [github.com/astral-sh/ruff-pre-commit: v0.6.7 → v0.6.9](https://github.com/astral-sh/ruff-pre-commit/compare/v0.6.7...v0.6.9)
---
 .pre-commit-config.yaml | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/.pre-commit-config.yaml b/.pre-commit-config.yaml
index 94ea09a8..58fdbcce 100644
--- a/.pre-commit-config.yaml
+++ b/.pre-commit-config.yaml
@@ -2,7 +2,7 @@
 # See https://pre-commit.com/hooks.html for more hooks
 repos:
 - repo: https://github.com/pre-commit/pre-commit-hooks
-  rev: v4.6.0
+  rev: v5.0.0
   hooks:
   - id: trailing-whitespace
   - id: end-of-file-fixer
@@ -11,7 +11,7 @@ repos:
     exclude: ^example/AlN-LDA/
 
 - repo: https://github.com/astral-sh/ruff-pre-commit
-  rev: v0.6.7
+  rev: v0.6.9
   hooks:
   - id: ruff
     args: [ "--fix", "--show-fixes" ]

From 7ed6007be45f4ffc78c716f31aef32b6b98d6214 Mon Sep 17 00:00:00 2001
From: Atsushi Togo <atz.togo@gmail.com>
Date: Wed, 9 Oct 2024 09:57:12 +0900
Subject: [PATCH 2/9] pytest with pypolymlp on github workflow

---
 .github/workflows/phono3py-pytest-conda.yml | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/.github/workflows/phono3py-pytest-conda.yml b/.github/workflows/phono3py-pytest-conda.yml
index 1df4515a..7ae69ede 100644
--- a/.github/workflows/phono3py-pytest-conda.yml
+++ b/.github/workflows/phono3py-pytest-conda.yml
@@ -24,7 +24,7 @@ jobs:
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
-        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py scipy pytest spglib alm cmake c-compiler cxx-compiler
+        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py scipy pytest spglib alm cmake c-compiler cxx-compiler pypolymlp
     - name: Install phonopy develop branch
       run: |
         conda activate test

From fa769411821df3f470d358ec89fd12c0285fa3f1 Mon Sep 17 00:00:00 2001
From: Atsushi Togo <atz.togo@gmail.com>
Date: Wed, 9 Oct 2024 10:02:26 +0900
Subject: [PATCH 3/9] pytest with pypolymlp on github workflow

---
 .github/workflows/phono3py-pytest-conda.yml | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/.github/workflows/phono3py-pytest-conda.yml b/.github/workflows/phono3py-pytest-conda.yml
index 7ae69ede..b1941f55 100644
--- a/.github/workflows/phono3py-pytest-conda.yml
+++ b/.github/workflows/phono3py-pytest-conda.yml
@@ -24,7 +24,7 @@ jobs:
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
-        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py scipy pytest spglib alm cmake c-compiler cxx-compiler pypolymlp
+        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py "numpy>=2.1" scipy pytest spglib alm cmake c-compiler cxx-compiler pypolymlp symfc
     - name: Install phonopy develop branch
       run: |
         conda activate test

From 170ebde32c839fa1e9c1953cc393cc450762b3bc Mon Sep 17 00:00:00 2001
From: Atsushi Togo <atz.togo@gmail.com>
Date: Wed, 9 Oct 2024 10:07:32 +0900
Subject: [PATCH 4/9] pytest with pypolymlp on github workflow

---
 .github/workflows/phono3py-pytest-conda.yml | 9 ++++++++-
 1 file changed, 8 insertions(+), 1 deletion(-)

diff --git a/.github/workflows/phono3py-pytest-conda.yml b/.github/workflows/phono3py-pytest-conda.yml
index b1941f55..9ff51620 100644
--- a/.github/workflows/phono3py-pytest-conda.yml
+++ b/.github/workflows/phono3py-pytest-conda.yml
@@ -20,11 +20,18 @@ jobs:
     - uses: conda-incubator/setup-miniconda@v3
       with:
         miniforge-version: latest
-    - name: Install dependent packages
+    - name: Install dependent packages for python > 3.9
+      if: ${{ matrix.python-version != 3.9 }}
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
         conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py "numpy>=2.1" scipy pytest spglib alm cmake c-compiler cxx-compiler pypolymlp symfc
+    - name: Install dependent packages for python == 3.9
+      if: ${{ matrix.python-version == 3.9 }}
+      run: |
+        conda activate test
+        conda install --yes python=${{ matrix.python-version }}
+        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py scipy pytest spglib alm cmake c-compiler cxx-compiler symfc
     - name: Install phonopy develop branch
       run: |
         conda activate test

From 9fbaaacb422f1f24e7c75af46fb54c20ac89aead Mon Sep 17 00:00:00 2001
From: Atsushi Togo <atz.togo@gmail.com>
Date: Wed, 9 Oct 2024 10:12:39 +0900
Subject: [PATCH 5/9] pytest with pypolymlp on github workflow

---
 .../phono3py-pytest-conda-mkl-phphmtblas.yml          |  2 +-
 .github/workflows/phono3py-pytest-conda-mkl-v2.yml    |  2 +-
 .github/workflows/phono3py-pytest-conda-mkl.yml       |  2 +-
 .github/workflows/phono3py-pytest-conda-numpy2.yml    |  2 +-
 .../workflows/phono3py-pytest-conda-phphmtblas.yml    |  2 +-
 .github/workflows/phono3py-pytest-conda.yml           | 11 +++++++++--
 6 files changed, 14 insertions(+), 7 deletions(-)

diff --git a/.github/workflows/phono3py-pytest-conda-mkl-phphmtblas.yml b/.github/workflows/phono3py-pytest-conda-mkl-phphmtblas.yml
index a262cdf6..49e60dd4 100644
--- a/.github/workflows/phono3py-pytest-conda-mkl-phphmtblas.yml
+++ b/.github/workflows/phono3py-pytest-conda-mkl-phphmtblas.yml
@@ -24,7 +24,7 @@ jobs:
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
-        conda install --yes matplotlib-base pyyaml "libblas=*=*mkl" mkl-include h5py scipy pytest spglib alm cmake c-compiler cxx-compiler
+        conda install --yes matplotlib-base pyyaml "libblas=*=*mkl" mkl-include h5py scipy pytest spglib alm cmake c-compiler cxx-compiler pypolymlp
     - name: Install symfc develop branch
       run: |
         conda activate test
diff --git a/.github/workflows/phono3py-pytest-conda-mkl-v2.yml b/.github/workflows/phono3py-pytest-conda-mkl-v2.yml
index a6d5c8d5..c1e95df6 100644
--- a/.github/workflows/phono3py-pytest-conda-mkl-v2.yml
+++ b/.github/workflows/phono3py-pytest-conda-mkl-v2.yml
@@ -24,7 +24,7 @@ jobs:
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
-        conda install --yes matplotlib-base pyyaml "libblas=*=*mkl" mkl-include h5py scipy pytest spglib alm cmake c-compiler cxx-compiler
+        conda install --yes matplotlib-base pyyaml "libblas=*=*mkl" mkl-include h5py scipy pytest spglib alm cmake c-compiler cxx-compiler pypolymlp
     - name: Install symfc develop branch
       run: |
         conda activate test
diff --git a/.github/workflows/phono3py-pytest-conda-mkl.yml b/.github/workflows/phono3py-pytest-conda-mkl.yml
index f42b8264..c1ae62de 100644
--- a/.github/workflows/phono3py-pytest-conda-mkl.yml
+++ b/.github/workflows/phono3py-pytest-conda-mkl.yml
@@ -24,7 +24,7 @@ jobs:
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
-        conda install --yes matplotlib-base pyyaml "libblas=*=*mkl" mkl-include h5py scipy pytest codecov pytest-cov spglib alm cmake c-compiler cxx-compiler
+        conda install --yes matplotlib-base pyyaml "libblas=*=*mkl" mkl-include h5py scipy pytest codecov pytest-cov spglib alm cmake c-compiler cxx-compiler pypolymlp
     - name: Install symfc develop branch
       run: |
         conda activate test
diff --git a/.github/workflows/phono3py-pytest-conda-numpy2.yml b/.github/workflows/phono3py-pytest-conda-numpy2.yml
index 85aee8be..c2e05745 100644
--- a/.github/workflows/phono3py-pytest-conda-numpy2.yml
+++ b/.github/workflows/phono3py-pytest-conda-numpy2.yml
@@ -24,7 +24,7 @@ jobs:
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
-        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py "numpy=2" scipy pytest cmake c-compiler cxx-compiler
+        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py "numpy=2" scipy pytest cmake c-compiler cxx-compiler pypolymlp
     - name: Install spglib develop branch
       run: |
         conda activate test
diff --git a/.github/workflows/phono3py-pytest-conda-phphmtblas.yml b/.github/workflows/phono3py-pytest-conda-phphmtblas.yml
index 5771b33b..cd4d625e 100644
--- a/.github/workflows/phono3py-pytest-conda-phphmtblas.yml
+++ b/.github/workflows/phono3py-pytest-conda-phphmtblas.yml
@@ -24,7 +24,7 @@ jobs:
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
-        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py scipy pytest spglib alm cmake c-compiler cxx-compiler
+        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py scipy pytest spglib alm cmake c-compiler cxx-compiler pypolymlp
     - name: Install symfc develop branch
       run: |
         conda activate test
diff --git a/.github/workflows/phono3py-pytest-conda.yml b/.github/workflows/phono3py-pytest-conda.yml
index 9ff51620..66bb1385 100644
--- a/.github/workflows/phono3py-pytest-conda.yml
+++ b/.github/workflows/phono3py-pytest-conda.yml
@@ -25,13 +25,20 @@ jobs:
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
-        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py "numpy>=2.1" scipy pytest spglib alm cmake c-compiler cxx-compiler pypolymlp symfc
+        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py "numpy>=2.1" scipy pytest spglib alm cmake c-compiler cxx-compiler pypolymlp
     - name: Install dependent packages for python == 3.9
       if: ${{ matrix.python-version == 3.9 }}
       run: |
         conda activate test
         conda install --yes python=${{ matrix.python-version }}
-        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py scipy pytest spglib alm cmake c-compiler cxx-compiler symfc
+        conda install --yes matplotlib-base pyyaml "libblas=*=*openblas" openblas h5py scipy pytest spglib alm cmake c-compiler cxx-compiler
+    - name: Install symfc develop branch
+      run: |
+        conda activate test
+        git clone --depth 1 https://github.com/symfc/symfc.git
+        cd symfc
+        pip install -e . -vvv
+        cd ..
     - name: Install phonopy develop branch
       run: |
         conda activate test

From 6c183fcd4b3d2ea7c0a6b3ae316eabfd5b0e2953 Mon Sep 17 00:00:00 2001
From: Atsushi Togo <atz.togo@gmail.com>
Date: Wed, 9 Oct 2024 13:21:57 +0900
Subject: [PATCH 6/9] Parse forces with displacement position check

---
 phono3py/cui/create_force_sets.py | 115 +++++++++++++++++++++++++-----
 phono3py/cui/phono3py_argparse.py |   8 +++
 phono3py/cui/settings.py          |  14 ++++
 3 files changed, 120 insertions(+), 17 deletions(-)

diff --git a/phono3py/cui/create_force_sets.py b/phono3py/cui/create_force_sets.py
index ed5b654f..bcb672ee 100644
--- a/phono3py/cui/create_force_sets.py
+++ b/phono3py/cui/create_force_sets.py
@@ -41,12 +41,18 @@
 from typing import Optional
 
 import numpy as np
-from phonopy.cui.create_force_sets import check_number_of_force_files
+from phonopy.cui.create_force_sets import (
+    check_agreement_of_supercell_positions,
+    check_agreements_of_displacements,
+    check_number_of_force_files,
+)
 from phonopy.cui.load_helper import get_nac_params
 from phonopy.cui.phonopy_script import file_exists, files_exist, print_error
 from phonopy.file_IO import is_file_phonopy_yaml, parse_FORCE_SETS, write_FORCE_SETS
 from phonopy.interface.calculator import get_calc_dataset
+from phonopy.structure.atoms import PhonopyAtoms
 
+from phono3py.cui.settings import Phono3pySettings
 from phono3py.file_IO import (
     get_length_of_first_line,
     parse_FORCES_FC2,
@@ -60,7 +66,7 @@
 
 
 def create_FORCES_FC3_and_FORCES_FC2(
-    settings,
+    settings: Phono3pySettings,
     cell_filename: Optional[str],
     log_level: int = 0,
 ):
@@ -106,11 +112,16 @@ def create_FORCES_FC3_and_FORCES_FC2(
 
     if settings.create_forces_fc3 or settings.create_forces_fc3_file:
         calc_dataset_fc3 = _get_force_sets_fc3(
-            settings, ph3py_yaml.dataset, disp_filename, interface_mode, log_level
+            settings,
+            ph3py_yaml.dataset,
+            ph3py_yaml.supercell,
+            disp_filename,
+            interface_mode,
+            log_level,
         )
         if not calc_dataset_fc3["forces"]:
             if log_level:
-                print("%s could not be created." % "FORCES_FC3")
+                print('"FORCES_FC3" could not be created.')
                 print_error()
             sys.exit(1)
 
@@ -118,13 +129,14 @@ def create_FORCES_FC3_and_FORCES_FC2(
         calc_dataset_fc2 = _get_force_sets_fc2(
             settings,
             ph3py_yaml.phonon_dataset,
+            ph3py_yaml.phonon_supercell,
             disp_filename,
             interface_mode,
             log_level,
         )
         if not calc_dataset_fc2["forces"]:
             if log_level:
-                print("%s could not be created." % "FORCES_FC2")
+                print('"FORCES_FC2" could not be created.')
                 print_error()
             sys.exit(1)
     else:
@@ -157,7 +169,7 @@ def create_FORCES_FC3_and_FORCES_FC2(
                 filename="FORCES_FC3",
             )
             if log_level:
-                print("%s has been created." % "FORCES_FC3")
+                print('"FORCES_FC3" has been created.')
 
         if settings.create_forces_fc2:
             write_FORCES_FC2(
@@ -166,7 +178,7 @@ def create_FORCES_FC3_and_FORCES_FC2(
                 filename="FORCES_FC2",
             )
             if log_level:
-                print("%s has been created." % "FORCES_FC2")
+                print('"FORCES_FC2" has been created.')
 
 
 def create_FORCES_FC2_from_FORCE_SETS(log_level):
@@ -247,7 +259,12 @@ def create_FORCE_SETS_from_FORCES_FCx(
 
 
 def _get_force_sets_fc2(
-    settings, disp_dataset, disp_filename, interface_mode, log_level
+    settings: Phono3pySettings,
+    disp_dataset: dict,
+    supercell: PhonopyAtoms,
+    disp_filename: str,
+    interface_mode: Optional[str],
+    log_level: int,
 ) -> dict:
     interface_mode = settings.calculator
     if log_level:
@@ -277,9 +294,19 @@ def _get_force_sets_fc2(
         verbose=(log_level > 0),
     )
     force_sets = calc_dataset["forces"]
+    if "points" in calc_dataset:
+        if filename := check_agreements_of_displacements(
+            supercell, disp_dataset, calc_dataset["points"], force_filenames
+        ):
+            raise RuntimeError(
+                f'Displacements don\'t match with atomic positions in "{filename}".'
+            )
 
-    if settings.subtract_forces:
-        force_filename = settings.subtract_forces
+    if settings.subtract_forces or settings.subtract_forces_fc2:
+        if settings.subtract_forces_fc2:
+            force_filename = settings.subtract_forces_fc2
+        else:
+            force_filename = settings.subtract_forces
         file_exists(force_filename, log_level=log_level)
         calc_dataset_zero = get_calc_dataset(
             interface_mode,
@@ -289,12 +316,26 @@ def _get_force_sets_fc2(
             ],
             verbose=(log_level > 0),
         )
+        if "points" in calc_dataset_zero:
+            if check_agreement_of_supercell_positions(
+                supercell, calc_dataset_zero["points"][0]
+            ):
+                raise RuntimeError(
+                    "Supercell doesn't match with atomic positions in "
+                    f'"{force_filename}".'
+                )
         force_set_zero = calc_dataset_zero["forces"][0]
         for fs in force_sets:
             fs -= force_set_zero
 
         if log_level > 0:
-            print("Forces in '{force_filename}' were subtracted from supercell forces.")
+            print(
+                f'Forces in "{force_filename}" were subtracted from supercell forces.'
+            )
+
+    if log_level > 0 and "first_atoms" not in disp_dataset:
+        if len(force_filenames) < num_disps:
+            print("** Number of supercell files is less than displacements. **")
 
     if log_level > 0:
         print("")
@@ -303,7 +344,12 @@ def _get_force_sets_fc2(
 
 
 def _get_force_sets_fc3(
-    settings, disp_dataset, disp_filename, interface_mode, log_level
+    settings: Phono3pySettings,
+    disp_dataset: dict,
+    supercell: PhonopyAtoms,
+    disp_filename: str,
+    interface_mode: Optional[str],
+    log_level: int,
 ) -> dict:
     if log_level:
         print(f'FC3 Displacement dataset was read from "{disp_filename}".')
@@ -331,12 +377,14 @@ def _get_force_sets_fc3(
         file_exists(filename, log_level=log_level)
 
     if log_level > 0:
-        print(f"  Number of displacements: {num_disps}")
         print(f"  Number of supercell files: {len(force_filenames)}")
+        print(f'  Number of displacements in "{disp_filename}": {num_disps}')
 
-    if not check_number_of_force_files(num_disps, force_filenames, disp_filename):
+    if "first_atoms" in disp_dataset and not check_number_of_force_files(
+        num_disps, force_filenames, disp_filename
+    ):  # type-1
         calc_dataset = {"forces": []}
-    else:
+    else:  # type-2
         calc_dataset = get_calc_dataset(
             interface_mode,
             num_atoms,
@@ -344,11 +392,18 @@ def _get_force_sets_fc3(
             verbose=(log_level > 0),
         )
         force_sets = calc_dataset["forces"]
+        if "points" in calc_dataset:
+            if filename := check_agreements_of_displacements(
+                supercell, disp_dataset, calc_dataset["points"], force_filenames
+            ):
+                raise RuntimeError(
+                    f'Displacements don\'t match with atomic positions in "{filename}".'
+                )
 
     if settings.subtract_forces:
         force_filename = settings.subtract_forces
         file_exists(force_filename, log_level=log_level)
-        calc_dataset = get_calc_dataset(
+        calc_dataset_zero = get_calc_dataset(
             interface_mode,
             num_atoms,
             [
@@ -356,7 +411,16 @@ def _get_force_sets_fc3(
             ],
             verbose=(log_level > 0),
         )
-        force_set_zero = calc_dataset["forces"][0]
+        force_set_zero = calc_dataset_zero["forces"][0]
+        if "points" in calc_dataset_zero:
+            if check_agreement_of_supercell_positions(
+                supercell, calc_dataset_zero["points"][0]
+            ):
+                raise RuntimeError(
+                    "Supercell doesn't match with atomic positions in "
+                    f'"{force_filename}".'
+                )
+
         for fs in force_sets:
             fs -= force_set_zero
 
@@ -365,6 +429,10 @@ def _get_force_sets_fc3(
                 f"Forces in '{force_filename}' were subtracted from supercell forces."
             )
 
+    if log_level > 0 and "first_atoms" not in disp_dataset:
+        if len(force_filenames) < num_disps:
+            print("** Number of supercell files is less than displacements. **")
+
     if log_level > 0:
         print("")
 
@@ -377,6 +445,7 @@ def _set_forces_and_nac_params(
     calc_dataset_fc3: dict,
     calc_dataset_fc2: Optional[dict],
 ):
+    """Set forces, energies and nac_params to phono3py_yaml."""
     if "first_atoms" in ph3py_yaml.dataset:
         count = len(ph3py_yaml.dataset["first_atoms"])
         for i, d1 in enumerate(ph3py_yaml.dataset["first_atoms"]):
@@ -400,6 +469,12 @@ def _set_forces_and_nac_params(
         ph3py_yaml.dataset["supercell_energies"] = np.array(
             calc_dataset_fc3["supercell_energies"], dtype="double"
         )
+        if len(ph3py_yaml.dataset["forces"]) != len(
+            ph3py_yaml.dataset["displacements"]
+        ):
+            ph3py_yaml.dataset["displacements"] = ph3py_yaml.dataset["displacements"][
+                : len(ph3py_yaml.dataset["forces"])
+            ]
 
     if settings.create_forces_fc2 and calc_dataset_fc2:
         if "first_atoms" in ph3py_yaml.phonon_dataset:
@@ -416,6 +491,12 @@ def _set_forces_and_nac_params(
             ph3py_yaml.phonon_dataset["supercell_energies"] = np.array(
                 calc_dataset_fc2["supercell_energies"], dtype="double"
             )
+            if len(ph3py_yaml.phonon_dataset["forces"]) != len(
+                ph3py_yaml.phonon_dataset["displacements"]
+            ):
+                ph3py_yaml.phonon_dataset["displacements"] = ph3py_yaml.phonon_dataset[
+                    "displacements"
+                ][: len(ph3py_yaml.phonon_dataset["forces"])]
 
     if ph3py_yaml.nac_params is None:
         nac_params = get_nac_params(primitive=ph3py_yaml.primitive)
diff --git a/phono3py/cui/phono3py_argparse.py b/phono3py/cui/phono3py_argparse.py
index 32c392f4..8ac6400f 100644
--- a/phono3py/cui/phono3py_argparse.py
+++ b/phono3py/cui/phono3py_argparse.py
@@ -159,6 +159,14 @@ def get_parser(fc_symmetry=False, is_nac=False, load_phono3py_yaml=False):
         default=None,
         help="Subtract recidual forces from supercell forces",
     )
+    parser.add_argument(
+        "--cfz-fc2",
+        "--subtract-forces-fc2",
+        metavar="FILE",
+        dest="subtract_forces_fc2",
+        default=None,
+        help="Subtract recidual forces from supercell forces for fc2",
+    )
     parser.add_argument(
         "--cfc",
         "--compact-fc",
diff --git a/phono3py/cui/settings.py b/phono3py/cui/settings.py
index 45f71814..096ee8d2 100644
--- a/phono3py/cui/settings.py
+++ b/phono3py/cui/settings.py
@@ -95,6 +95,7 @@ class Phono3pySettings(Settings):
         "solve_collective_phonon": False,
         "emulate_v2": False,
         "subtract_forces": None,
+        "subtract_forces_fc2": None,
         "temperatures": None,
         "use_ave_pp": False,
         "use_grg": False,
@@ -321,6 +322,10 @@ def set_subtract_forces(self, val):
         """Set subtract_forces."""
         self._v["subtract_forces"] = val
 
+    def set_subtract_forces_fc2(self, val):
+        """Set subtract_forces_fc2."""
+        self._v["subtract_forces_fc2"] = val
+
     def set_temperatures(self, val):
         """Set temperatures."""
         self._v["temperatures"] = val
@@ -605,6 +610,10 @@ def _read_options(self):
             if self._args.subtract_forces:
                 self._confs["subtract_forces"] = self._args.subtract_forces
 
+        if "subtract_forces_fc2" in self._args:
+            if self._args.subtract_forces_fc2:
+                self._confs["subtract_forces_fc2"] = self._args.subtract_forces_fc2
+
         if "temperatures" in self._args:
             if self._args.temperatures is not None:
                 self._confs["temperatures"] = " ".join(self._args.temperatures)
@@ -716,6 +725,7 @@ def _parse_conf(self):
                 "create_forces_fc3_file",
                 "output_yaml_filename",
                 "subtract_forces",
+                "subtract_forces_fc2",
             ):
                 self.set_parameter(conf_key, confs[conf_key])
 
@@ -1009,6 +1019,10 @@ def _set_settings(self):
         if "subtract_forces" in params:
             self._settings.set_subtract_forces(params["subtract_forces"])
 
+        # Subtract residual forces to create FORCES_FC2
+        if "subtract_forces_fc2" in params:
+            self._settings.set_subtract_forces_fc2(params["subtract_forces_fc2"])
+
         # Temperatures for scatterings
         if "temperatures" in params:
             self._settings.set_temperatures(params["temperatures"])

From fc0d833368744cb405feb502da0e4e2bfc924c49 Mon Sep 17 00:00:00 2001
From: "pre-commit-ci[bot]"
 <66853113+pre-commit-ci[bot]@users.noreply.github.com>
Date: Mon, 21 Oct 2024 23:49:17 +0000
Subject: [PATCH 7/9] [pre-commit.ci] pre-commit autoupdate
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

updates:
- [github.com/astral-sh/ruff-pre-commit: v0.6.9 → v0.7.0](https://github.com/astral-sh/ruff-pre-commit/compare/v0.6.9...v0.7.0)
---
 .pre-commit-config.yaml | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/.pre-commit-config.yaml b/.pre-commit-config.yaml
index 58fdbcce..6d3c4221 100644
--- a/.pre-commit-config.yaml
+++ b/.pre-commit-config.yaml
@@ -11,7 +11,7 @@ repos:
     exclude: ^example/AlN-LDA/
 
 - repo: https://github.com/astral-sh/ruff-pre-commit
-  rev: v0.6.9
+  rev: v0.7.0
   hooks:
   - id: ruff
     args: [ "--fix", "--show-fixes" ]

From bbaa1cdd47f11aa0f31600bd84f61b3dad89f57d Mon Sep 17 00:00:00 2001
From: Atsushi Togo <atz.togo@gmail.com>
Date: Tue, 29 Oct 2024 10:13:42 +0900
Subject: [PATCH 8/9] Update phonopy version dependenty

---
 pyproject.toml | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/pyproject.toml b/pyproject.toml
index af837355..787f70be 100644
--- a/pyproject.toml
+++ b/pyproject.toml
@@ -16,7 +16,7 @@ dependencies = [
     "matplotlib>=2.2.2",
     "h5py>=3.0",
     "spglib>=2.3",
-    "phonopy>=2.28,<2.29",
+    "phonopy>=2.29,<2.30",
 ]
 license = { file = "LICENSE" }
 

From 1a106fdee92d996bf900e93a690ab49e603fbe71 Mon Sep 17 00:00:00 2001
From: Atsushi Togo <atz.togo@gmail.com>
Date: Sun, 3 Nov 2024 15:27:50 +0900
Subject: [PATCH 9/9] Set version 3.6.0

---
 doc/changelog.md    | 4 ++++
 doc/conf.py         | 4 ++--
 phono3py/version.py | 2 +-
 3 files changed, 7 insertions(+), 3 deletions(-)

diff --git a/doc/changelog.md b/doc/changelog.md
index 71a671f2..4faa2c1f 100644
--- a/doc/changelog.md
+++ b/doc/changelog.md
@@ -2,6 +2,10 @@
 
 # Change Log
 
+## Nov-3-2024: Version 3.6.0
+
+- Maintenance release.
+
 ## Sep-24-2024: Version 3.5.2
 
 - Fix a memory leak.
diff --git a/doc/conf.py b/doc/conf.py
index f186e32a..2f64b5f2 100644
--- a/doc/conf.py
+++ b/doc/conf.py
@@ -58,9 +58,9 @@
 # built documents.
 #
 # The short X.Y version.
-version = "3.5"
+version = "3.6"
 # The full version, including alpha/beta/rc tags.
-release = "3.5.2"
+release = "3.6.0"
 
 # The language for content autogenerated by Sphinx. Refer to documentation
 # for a list of supported languages.
diff --git a/phono3py/version.py b/phono3py/version.py
index d4df8bf1..ea0e6101 100644
--- a/phono3py/version.py
+++ b/phono3py/version.py
@@ -34,4 +34,4 @@
 # ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 # POSSIBILITY OF SUCH DAMAGE.
 
-__version__ = "3.5.2"
+__version__ = "3.6.0"