MultiOptPy 1.20.4__tar.gz → 1.20.6__tar.gz

{multioptpy-1.20.4 → multioptpy-1.20.6/MultiOptPy.egg-info}/PKG-INFO

@@ -1,21 +1,21 @@
 Metadata-Version: 2.4
 Name: MultiOptPy
-Version: 1.20.4
+Version: 1.20.6
 Summary: Multifunctional geometry optimization tools for quantum chemical calculations.
 Author-email: ss0832 <highlighty876@gmail.com>
-License: GPLv3
+License-Expression: GPL-3.0-or-later
 Requires-Python: >=3.12
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: numpy~=2.2.0
+Requires-Dist: numpy>=2.2.0
 Requires-Dist: scipy>=1.13.0
-Requires-Dist: matplotlib~=3.10.0
+Requires-Dist: matplotlib>=3.10.0
 Requires-Dist: torch~=2.6.0
-Requires-Dist: pyscf~=2.9.0
-Requires-Dist: tblite~=0.4.0
+Requires-Dist: pyscf>=2.9.0
+Requires-Dist: tblite>=0.4.0
 Requires-Dist: ase~=3.26.0
 Requires-Dist: fairchem-core~=2.7.0
-Requires-Dist: sympy~=1.13.0
+Requires-Dist: sympy>=1.13.0
 Dynamic: license-file
 
 # MultiOptPy
@@ -29,7 +29,7 @@ Dynamic: license-file
 [![Buy Me A Coffee](https://img.shields.io/badge/Buy%20Me%20a%20Coffee-ffdd00?style=for-the-badge&logo=buy-me-a-coffee&logoColor=black)](https://buymeacoffee.com/ss0832)
 
 [![PyPI Downloads](https://static.pepy.tech/personalized-badge/multioptpy?period=total&units=INTERNATIONAL_SYSTEM&left_color=BLACK&right_color=GREEN&left_text=downloads)](https://pepy.tech/projects/multioptpy)
-[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.17970774.svg)](https://doi.org/10.5281/zenodo.17970774)
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.17973395.svg)](https://doi.org/10.5281/zenodo.17973395)
 
 If this tool helped your studies, education, or saved your time, I'd appreciate a coffee!
 Your support serves as a great encouragement for this personal project and fuels my next journey.
@@ -80,9 +80,9 @@ conda create -n test_mop python=3.12.7
 conda activate test_mop
 
 ## 3. Download and install MultiOptPy:
-wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.3.zip
-unzip v1.20.3.zip
-cd MultiOptPy-1.20.3
+wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.5.zip
+unzip v1.20.5.zip
+cd MultiOptPy-1.20.5
 pip install -r requirements.txt
 
 ## 4. Copy the test configuration file and run the AutoTS workflow:
@@ -107,10 +107,10 @@ python run_autots.py aldol_rxn.xyz -cfg config_autots_run_xtb_test.json
 # Installation via pip (Linux)
 conda create -n <env-name> python=3.12 pip
 conda activate <env-name>
-pip install git+https://github.com/ss0832/MultiOptPy.git@v1.20.3
-wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.3.zip
-unzip v1.20.3.zip
-cd MultiOptPy-1.20.3
+pip install git+https://github.com/ss0832/MultiOptPy.git@v1.20.5
+wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.5.zip
+unzip v1.20.5.zip
+cd MultiOptPy-1.20.5
 
 ## 💻 Command Line Interface (CLI) Functionality (v1.20.2)
 # The following eight core functionalities are available as direct executable commands in your terminal after installation:
@@ -157,7 +157,7 @@ References are given in the source code.
 After downloading the repository using git clone or similar commands, move to the generated directory and run the following:
 python command
 ```
-python optmain.py SN2.xyz -ma 150 1 6 -pyscf -elec 0 -spin 0 -opt rsirfo_block_fsb -modelhess
+python optmain.py SN2.xyz -ma 150 1 6 -pyscf -elec -1 -spin 0 -opt rsirfo_block_fsb -modelhess
 ```
 CLI command (arbitrary directory)
 ```
@@ -412,13 +412,13 @@ If you use MultiOptPy in your research, please cite it as follows:
   month        = dec,
   year         = 2025,
   publisher    = {Zenodo},
-  version      = {v1.20.3},
-  doi          = {10.5281/zenodo.17970774},
-  url          = {https://doi.org/10.5281/zenodo.17970774}
+  version      = {v1.20.4},
+  doi          = {10.5281/zenodo.17973395},
+  url          = {https://doi.org/10.5281/zenodo.17973395}
 }
 ```
 ```
-ss0832. (2025). MultiOptPy: Multifunctional geometry optimization tools for quantum chemical calculations (v1.20.3). Zenodo. https://doi.org/10.5281/zenodo.17970774
+ss0832. (2025). MultiOptPy: Multifunctional geometry optimization tools for quantum chemical calculations (v1.20.4). Zenodo. https://doi.org/10.5281/zenodo.17973395
 ```
 
 ## Setting Up an Environment for Using NNP(UMA) on Windows 11

{multioptpy-1.20.4 → multioptpy-1.20.6}/MultiOptPy.egg-info/SOURCES.txt

@@ -34,6 +34,7 @@ multioptpy/Calculator/ase_tools/__init__.py
 multioptpy/Calculator/ase_tools/fairchem.py
 multioptpy/Calculator/ase_tools/gamess.py
 multioptpy/Calculator/ase_tools/gaussian.py
+multioptpy/Calculator/ase_tools/gxtb_dev.py
 multioptpy/Calculator/ase_tools/mace.py
 multioptpy/Calculator/ase_tools/mopac.py
 multioptpy/Calculator/ase_tools/nwchem.py
@@ -175,6 +176,7 @@ multioptpy/OtherMethod/dimer.py
 multioptpy/OtherMethod/elastic_image_pair.py
 multioptpy/OtherMethod/modelfunction.py
 multioptpy/OtherMethod/newton_traj.py
+multioptpy/OtherMethod/spring_pair_method.py
 multioptpy/OtherMethod/twopshs.py
 multioptpy/PESAnalyzer/__init__.py
 multioptpy/PESAnalyzer/calc_irc_curvature.py

multioptpy-1.20.6/MultiOptPy.egg-info/requires.txt

@@ -0,0 +1,9 @@
+numpy>=2.2.0
+scipy>=1.13.0
+matplotlib>=3.10.0
+torch~=2.6.0
+pyscf>=2.9.0
+tblite>=0.4.0
+ase~=3.26.0
+fairchem-core~=2.7.0
+sympy>=1.13.0

{multioptpy-1.20.4/MultiOptPy.egg-info → multioptpy-1.20.6}/PKG-INFO

@@ -1,21 +1,21 @@
 Metadata-Version: 2.4
 Name: MultiOptPy
-Version: 1.20.4
+Version: 1.20.6
 Summary: Multifunctional geometry optimization tools for quantum chemical calculations.
 Author-email: ss0832 <highlighty876@gmail.com>
-License: GPLv3
+License-Expression: GPL-3.0-or-later
 Requires-Python: >=3.12
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: numpy~=2.2.0
+Requires-Dist: numpy>=2.2.0
 Requires-Dist: scipy>=1.13.0
-Requires-Dist: matplotlib~=3.10.0
+Requires-Dist: matplotlib>=3.10.0
 Requires-Dist: torch~=2.6.0
-Requires-Dist: pyscf~=2.9.0
-Requires-Dist: tblite~=0.4.0
+Requires-Dist: pyscf>=2.9.0
+Requires-Dist: tblite>=0.4.0
 Requires-Dist: ase~=3.26.0
 Requires-Dist: fairchem-core~=2.7.0
-Requires-Dist: sympy~=1.13.0
+Requires-Dist: sympy>=1.13.0
 Dynamic: license-file
 
 # MultiOptPy
@@ -29,7 +29,7 @@ Dynamic: license-file
 [![Buy Me A Coffee](https://img.shields.io/badge/Buy%20Me%20a%20Coffee-ffdd00?style=for-the-badge&logo=buy-me-a-coffee&logoColor=black)](https://buymeacoffee.com/ss0832)
 
 [![PyPI Downloads](https://static.pepy.tech/personalized-badge/multioptpy?period=total&units=INTERNATIONAL_SYSTEM&left_color=BLACK&right_color=GREEN&left_text=downloads)](https://pepy.tech/projects/multioptpy)
-[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.17970774.svg)](https://doi.org/10.5281/zenodo.17970774)
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.17973395.svg)](https://doi.org/10.5281/zenodo.17973395)
 
 If this tool helped your studies, education, or saved your time, I'd appreciate a coffee!
 Your support serves as a great encouragement for this personal project and fuels my next journey.
@@ -80,9 +80,9 @@ conda create -n test_mop python=3.12.7
 conda activate test_mop
 
 ## 3. Download and install MultiOptPy:
-wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.3.zip
-unzip v1.20.3.zip
-cd MultiOptPy-1.20.3
+wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.5.zip
+unzip v1.20.5.zip
+cd MultiOptPy-1.20.5
 pip install -r requirements.txt
 
 ## 4. Copy the test configuration file and run the AutoTS workflow:
@@ -107,10 +107,10 @@ python run_autots.py aldol_rxn.xyz -cfg config_autots_run_xtb_test.json
 # Installation via pip (Linux)
 conda create -n <env-name> python=3.12 pip
 conda activate <env-name>
-pip install git+https://github.com/ss0832/MultiOptPy.git@v1.20.3
-wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.3.zip
-unzip v1.20.3.zip
-cd MultiOptPy-1.20.3
+pip install git+https://github.com/ss0832/MultiOptPy.git@v1.20.5
+wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.5.zip
+unzip v1.20.5.zip
+cd MultiOptPy-1.20.5
 
 ## 💻 Command Line Interface (CLI) Functionality (v1.20.2)
 # The following eight core functionalities are available as direct executable commands in your terminal after installation:
@@ -157,7 +157,7 @@ References are given in the source code.
 After downloading the repository using git clone or similar commands, move to the generated directory and run the following:
 python command
 ```
-python optmain.py SN2.xyz -ma 150 1 6 -pyscf -elec 0 -spin 0 -opt rsirfo_block_fsb -modelhess
+python optmain.py SN2.xyz -ma 150 1 6 -pyscf -elec -1 -spin 0 -opt rsirfo_block_fsb -modelhess
 ```
 CLI command (arbitrary directory)
 ```
@@ -412,13 +412,13 @@ If you use MultiOptPy in your research, please cite it as follows:
   month        = dec,
   year         = 2025,
   publisher    = {Zenodo},
-  version      = {v1.20.3},
-  doi          = {10.5281/zenodo.17970774},
-  url          = {https://doi.org/10.5281/zenodo.17970774}
+  version      = {v1.20.4},
+  doi          = {10.5281/zenodo.17973395},
+  url          = {https://doi.org/10.5281/zenodo.17973395}
 }
 ```
 ```
-ss0832. (2025). MultiOptPy: Multifunctional geometry optimization tools for quantum chemical calculations (v1.20.3). Zenodo. https://doi.org/10.5281/zenodo.17970774
+ss0832. (2025). MultiOptPy: Multifunctional geometry optimization tools for quantum chemical calculations (v1.20.4). Zenodo. https://doi.org/10.5281/zenodo.17973395
 ```
 
 ## Setting Up an Environment for Using NNP(UMA) on Windows 11

{multioptpy-1.20.4 → multioptpy-1.20.6}/README.md

@@ -9,7 +9,7 @@
 [![Buy Me A Coffee](https://img.shields.io/badge/Buy%20Me%20a%20Coffee-ffdd00?style=for-the-badge&logo=buy-me-a-coffee&logoColor=black)](https://buymeacoffee.com/ss0832)
 
 [![PyPI Downloads](https://static.pepy.tech/personalized-badge/multioptpy?period=total&units=INTERNATIONAL_SYSTEM&left_color=BLACK&right_color=GREEN&left_text=downloads)](https://pepy.tech/projects/multioptpy)
-[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.17970774.svg)](https://doi.org/10.5281/zenodo.17970774)
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.17973395.svg)](https://doi.org/10.5281/zenodo.17973395)
 
 If this tool helped your studies, education, or saved your time, I'd appreciate a coffee!
 Your support serves as a great encouragement for this personal project and fuels my next journey.
@@ -60,9 +60,9 @@ conda create -n test_mop python=3.12.7
 conda activate test_mop
 
 ## 3. Download and install MultiOptPy:
-wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.3.zip
-unzip v1.20.3.zip
-cd MultiOptPy-1.20.3
+wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.5.zip
+unzip v1.20.5.zip
+cd MultiOptPy-1.20.5
 pip install -r requirements.txt
 
 ## 4. Copy the test configuration file and run the AutoTS workflow:
@@ -87,10 +87,10 @@ python run_autots.py aldol_rxn.xyz -cfg config_autots_run_xtb_test.json
 # Installation via pip (Linux)
 conda create -n <env-name> python=3.12 pip
 conda activate <env-name>
-pip install git+https://github.com/ss0832/MultiOptPy.git@v1.20.3
-wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.3.zip
-unzip v1.20.3.zip
-cd MultiOptPy-1.20.3
+pip install git+https://github.com/ss0832/MultiOptPy.git@v1.20.5
+wget https://github.com/ss0832/MultiOptPy/archive/refs/tags/v1.20.5.zip
+unzip v1.20.5.zip
+cd MultiOptPy-1.20.5
 
 ## 💻 Command Line Interface (CLI) Functionality (v1.20.2)
 # The following eight core functionalities are available as direct executable commands in your terminal after installation:
@@ -137,7 +137,7 @@ References are given in the source code.
 After downloading the repository using git clone or similar commands, move to the generated directory and run the following:
 python command
 ```
-python optmain.py SN2.xyz -ma 150 1 6 -pyscf -elec 0 -spin 0 -opt rsirfo_block_fsb -modelhess
+python optmain.py SN2.xyz -ma 150 1 6 -pyscf -elec -1 -spin 0 -opt rsirfo_block_fsb -modelhess
 ```
 CLI command (arbitrary directory)
 ```
@@ -392,13 +392,13 @@ If you use MultiOptPy in your research, please cite it as follows:
   month        = dec,
   year         = 2025,
   publisher    = {Zenodo},
-  version      = {v1.20.3},
-  doi          = {10.5281/zenodo.17970774},
-  url          = {https://doi.org/10.5281/zenodo.17970774}
+  version      = {v1.20.4},
+  doi          = {10.5281/zenodo.17973395},
+  url          = {https://doi.org/10.5281/zenodo.17973395}
 }
 ```
 ```
-ss0832. (2025). MultiOptPy: Multifunctional geometry optimization tools for quantum chemical calculations (v1.20.3). Zenodo. https://doi.org/10.5281/zenodo.17970774
+ss0832. (2025). MultiOptPy: Multifunctional geometry optimization tools for quantum chemical calculations (v1.20.4). Zenodo. https://doi.org/10.5281/zenodo.17973395
 ```
 
 ## Setting Up an Environment for Using NNP(UMA) on Windows 11

{multioptpy-1.20.4 → multioptpy-1.20.6}/multioptpy/Calculator/ase_calculation_tools.py

@@ -22,6 +22,8 @@ from multioptpy.Calculator.ase_tools.fairchem import ASE_FAIRCHEM
 from multioptpy.Calculator.ase_tools.mace import ASE_MACE
 from multioptpy.Calculator.ase_tools.mopac import ASE_MOPAC
 from multioptpy.Calculator.ase_tools.pygfn0 import ASE_GFN0
+from multioptpy.Calculator.ase_tools.gxtb_dev import ASE_gxTB_Dev
+
 
 """
 referrence:
@@ -89,7 +91,9 @@ class Calculation:
         if self.software_type == "gaussian":
             print("Calculating exact Hessian using Gaussian...")
             exact_hess = calc_obj.calc_analytic_hessian()  # in hartree/Bohr^2
-        
+        elif self.software_type == "orca":
+            hess_path = calc_obj.run_frequency_analysis()
+            exact_hess = calc_obj.get_hessian_matrix(hess_path)
         else:
             vib = Vibrations(atoms=calc_obj.atom_obj, delta=0.001, name="z_hess_"+timestamp)
             vib.run()
@@ -124,7 +128,7 @@ class Calculation:
             file_list = glob.glob(file_directory+"/*_[0-9].xyz")
     
         for num, input_file in enumerate(file_list):
-            try:
+            if True:#try:
                 if geom_num_list is None:
                     positions, _, electric_charge_and_multiplicity = xyz2list(input_file, electric_charge_and_multiplicity)
                 else:
@@ -155,11 +159,11 @@ class Calculation:
                 elif iter % self.FC_COUNT == 0 or self.hessian_flag:
                     # exact numerical hessian
                     _ = self.calc_exact_hess(calc_obj, positions, element_list)             
-            except Exception as error:
-                print(error)
-                print("This molecule could not be optimized.")
-                finish_frag = True
-                return np.array([0]), np.array([0]), np.array([0]), finish_frag 
+            #except Exception as error:
+            #    print(error)
+            #    print("This molecule could not be optimized.")
+            #    finish_frag = True
+            #    return np.array([0]), np.array([0]), np.array([0]), finish_frag 
 
         positions /= self.bohr2angstroms
         self.energy = e
@@ -286,7 +290,10 @@ class ASEEngine(CalculationEngine):
                     timestamp = datetime.datetime.now().strftime("%Y_%m_%d_%H_%M_%S_%f")[:-2]
                     if software_type == "gaussian":
                         exact_hess = calc_obj.calc_analytic_hessian()  # in hartree/Bohr^2
-                        exact_hess = exact_hess 
+                   
+                    elif software_type == "orca":
+                        hess_path = calc_obj.run_frequency_analysis()
+                        exact_hess = calc_obj.get_hessian_matrix(hess_path)
                     else:
                         vib = Vibrations(atoms=calc_obj.atom_obj, delta=0.001, name="z_hess_"+timestamp)
                         vib.run()
@@ -431,6 +438,12 @@ def setup_calculator(atom_obj, software_type, electric_charge_and_multiplicity,
                              input_file=input_file)
         return calc_obj
 
+    if software_type == "gxtb_dev":
+        calc_obj = ASE_gxTB_Dev(atom_obj=atom_obj,
+                             electric_charge_and_multiplicity=electric_charge_and_multiplicity,
+                             software_type=software_type)
+        return calc_obj
+    
     # Unknown software type
     raise ValueError(f"Unsupported software type: {software_type}")
 

multioptpy-1.20.6/multioptpy/Calculator/ase_tools/gxtb_dev.py

@@ -0,0 +1,41 @@
+class ASE_gxTB_Dev:
+    """
+    Wrapper class to set up and run g-xTB (preliminary version) calculations via ASE.
+    $ pip install pygxtb==0.7.0
+    
+    """
+    def __init__(self, **kwargs):
+        
+        self.atom_obj = kwargs.get('atom_obj', None)
+        self.electric_charge_and_multiplicity = kwargs.get('electric_charge_and_multiplicity', None)
+        self.software_path = kwargs.get('software_path', None)
+        self.software_type = kwargs.get('software_type', None)
+        from pygxtb import PygxTB
+        self.pygxtb = PygxTB
+        
+    def set_calculator(self):
+        """
+        Sets the ASE calculator object based on software_type.
+        """
+        
+        charge = 0  # Default charge
+        if self.electric_charge_and_multiplicity is not None:
+            try:
+                # Get charge from [charge, multiplicity] list
+                charge = int(self.electric_charge_and_multiplicity[0])
+            except (IndexError, TypeError, ValueError):
+                print(f"Warning: Could not parse charge from {self.electric_charge_and_multiplicity}. Defaulting to 0.")
+                pass
+        
+        # Instantiate GFN0 class and pass the charge
+        gxtb_calc = self.pygxtb(charge=charge)
+        return gxtb_calc
+        
+
+    def run(self):
+        """
+        Attaches the calculator to the atoms object and returns it.
+        """
+        calc_obj = self.set_calculator()
+        self.atom_obj.calc = calc_obj
+        return self.atom_obj

multioptpy-1.20.6/multioptpy/Calculator/ase_tools/orca.py

@@ -0,0 +1,236 @@
+import os
+import shutil
+import subprocess
+import numpy as np
+from ase.calculators.orca import ORCA, OrcaProfile
+from ase.data import atomic_numbers
+
+"""
+Please specify the absolute path to the ORCA executable in software_path.conf using the format orca::<path>. For Linux, provide the path to the binary (e.g., /absolute/path/to/orca), and for Windows, provide the path to the executable file (e.g., C:\absolute\path\to\orca.exe).
+"""
+
+
+class ASE_ORCA:
+    TARGET_ORCA_VERSION = '6.1.0'
+
+    def __init__(self, **kwargs):
+        self.atom_obj = kwargs.get('atom_obj', None)
+        self.electric_charge_and_multiplicity = kwargs.get('electric_charge_and_multiplicity', None)
+        
+        # NOTE: self.input_file is updated in _setup_calculator to enforce 'orca.inp' in CWD.
+        raw_input_file = kwargs.get('input_file', 'orca.inp')
+        self.input_file = os.path.abspath(raw_input_file).replace('\\', '/')
+        
+        self.orca_path = kwargs.get('orca_path', None)
+        self.functional = kwargs.get('functional', 'B3LYP')
+        self.basis_set = kwargs.get('basis_set', 'def2-SVP')
+
+        # Optional ORCA input blocks (raw string)
+        self.orca_blocks = kwargs.get('orca_blocks', '')
+
+        # Auto-fix for unsupported Pople basis sets on heavy elements
+        self.auto_fix_basis = kwargs.get('auto_fix_basis', True)
+        self.heavy_atom_basis = kwargs.get('heavy_atom_basis', 'def2-SVP')
+
+    def _resolve_orca_exe(self, provided_path):
+        """
+        Resolve the absolute path to the ORCA executable.
+        Handles directories, stripping whitespace from config files, and WSL/Windows paths.
+        """
+        if not provided_path:
+            return None
+        
+        # CRITICAL FIX: Strip whitespace/newlines that might come from config parsing
+        clean_path = provided_path.strip()
+        # Expand ~ to home directory if present
+        clean_path = os.path.expanduser(clean_path)
+        clean_path = os.path.normpath(clean_path)
+
+        candidates = []
+        # If the path is a directory, look for the executable inside it
+        if os.path.isdir(clean_path):
+            candidates.append(os.path.join(clean_path, 'orca'))
+            candidates.append(os.path.join(clean_path, 'orca.exe'))
+        else:
+            # If it's a file path (or doesn't exist yet), use it as is
+            candidates.append(clean_path)
+
+        for candidate in candidates:
+            # Check if file exists
+            if os.path.exists(candidate) and os.path.isfile(candidate):
+                return os.path.abspath(candidate)
+            
+            # Check system PATH
+            resolved = shutil.which(candidate)
+            if resolved and os.path.exists(resolved):
+                return os.path.abspath(resolved)
+
+        # Use repr() in error message to reveal hidden characters like \n
+        candidate_reprs = [repr(c) for c in candidates]
+        raise FileNotFoundError(f"Cannot locate ORCA executable. Checked: {', '.join(candidate_reprs)}")
+
+    def _is_pople_basis(self, basis_name):
+        if not basis_name: return False
+        b = basis_name.strip().lower()
+        return b.startswith("6-31") or b.startswith("6-311")
+
+    def _get_heavy_elements_for_pople(self):
+        if self.atom_obj is None: return []
+        symbols = self.atom_obj.get_chemical_symbols()
+        return sorted({s for s in symbols if atomic_numbers.get(s, 0) > 30})
+
+    def _build_orca_blocks(self):
+        blocks = (self.orca_blocks or "").strip()
+        heavy_elements = []
+        if self._is_pople_basis(self.basis_set):
+            heavy_elements = self._get_heavy_elements_for_pople()
+
+        if heavy_elements:
+            if not self.auto_fix_basis:
+                raise ValueError("Unsupported Pople basis for heavy elements.")
+
+            basis_lines = ["%basis"]
+            for elem in heavy_elements:
+                # Per user instruction: No ECP, just GTO
+                basis_lines.append(f'  NewGTO {elem} "{self.heavy_atom_basis}"')
+            basis_lines.append("end")
+
+            if blocks:
+                blocks = blocks + "\n" + "\n".join(basis_lines)
+            else:
+                blocks = "\n".join(basis_lines)
+
+        return blocks if blocks else ""
+
+    def _print_orca_output_on_error(self):
+        """Helper to print ORCA output file content if it exists."""
+        if not hasattr(self, 'input_file'): return
+        
+        out_file = os.path.splitext(self.input_file)[0] + ".out"
+        if os.path.exists(out_file):
+            print(f"\n--- ORCA OUTPUT ({out_file}) ---")
+            try:
+                with open(out_file, 'r', encoding='utf-8', errors='replace') as f:
+                    content = f.read()
+                    print(content[-3000:] if len(content) > 3000 else content)
+            except Exception as e:
+                print(f"Failed to read output file: {e}")
+            print("--- END ORCA OUTPUT ---\n")
+        else:
+            print(f"\n--- ORCA OUTPUT NOT FOUND ({out_file}) ---\n")
+
+    def _setup_calculator(self, task_keyword):
+        # Force usage of Current Working Directory + "orca.inp"
+        cwd = os.getcwd()
+        label_path = os.path.join(cwd, 'orca').replace('\\', '/')
+        self.input_file = label_path + '.inp'
+        
+        print(f"DEBUG: ASE Label Path : {label_path}")
+
+        simple_input = f"{self.functional} {self.basis_set} {task_keyword}"
+        
+        profile_obj = None
+        if self.orca_path:
+            real_exe_path = self._resolve_orca_exe(self.orca_path)
+            orca_dir = os.path.dirname(real_exe_path)
+            path_env = os.environ.get('PATH', '')
+            if orca_dir not in path_env:
+                os.environ['PATH'] = orca_dir + os.pathsep + path_env
+
+            ase_safe_path = real_exe_path.replace('\\', '/')
+            profile_obj = OrcaProfile(ase_safe_path)
+            print(f"DEBUG: ORCA Executable: {ase_safe_path}")
+
+        orca_blocks = self._build_orca_blocks()
+        
+        calc = ORCA(
+            label=label_path,
+            profile=profile_obj,
+            charge=int(self.electric_charge_and_multiplicity[0]),
+            mult=int(self.electric_charge_and_multiplicity[1]),
+            orcasimpleinput=simple_input,
+            orcablocks=orca_blocks
+        )
+        
+        self.atom_obj.calc = calc
+        return self.atom_obj
+
+    def run(self):
+        self._setup_calculator("EnGrad")
+        print(f"--- Starting Gradient Calculation (ORCA {self.TARGET_ORCA_VERSION}) ---")
+        try:
+            forces = self.atom_obj.get_forces()
+            potential_energy = self.atom_obj.get_potential_energy()
+            print("Gradient calculation completed.")
+            return forces, potential_energy
+        except subprocess.CalledProcessError as e:
+            print(f"CRITICAL: ORCA execution failed with exit code {e.returncode}")
+            self._print_orca_output_on_error()
+            raise e
+        except Exception as e:
+            print(f"CRITICAL: An unexpected error occurred: {e}")
+            self._print_orca_output_on_error()
+            raise e
+
+    def run_frequency_analysis(self):
+        print(f"--- Starting Frequency Calculation (ORCA {self.TARGET_ORCA_VERSION}) ---")
+        self._setup_calculator("Freq")
+        try:
+            self.atom_obj.get_potential_energy()
+            print("Frequency calculation completed.")
+            # Use self.input_file to construct hess_path instead of relying on calc.label
+            hess_path = os.path.splitext(self.input_file)[0] + ".hess"
+            return hess_path
+        except subprocess.CalledProcessError as e:
+            print(f"CRITICAL: ORCA Frequency execution failed with exit code {e.returncode}")
+            self._print_orca_output_on_error()
+            raise e
+        except Exception as e:
+            print(f"CRITICAL: An unexpected error occurred during frequency analysis: {e}")
+            self._print_orca_output_on_error()
+            raise e
+
+    def get_hessian_matrix(self, hess_file_path=None):
+        if hess_file_path is None:
+            # Default to orca.hess in the same dir as input_file
+            input_dir = os.path.dirname(self.input_file)
+            hess_file_path = os.path.join(input_dir, 'orca.hess')
+
+        if not os.path.exists(hess_file_path):
+            raise FileNotFoundError(f"Hessian file not found: {hess_file_path}")
+
+        print(f"Reading Hessian from: {hess_file_path}")
+        
+        with open(hess_file_path, 'r', encoding='utf-8') as f:
+            lines = f.readlines()
+
+        hessian_matrix = None
+        iterator = iter(lines)
+        for line in iterator:
+            if "$hessian" in line:
+                dim_line = next(iterator).strip().split()
+                # Parse dimensions: square matrix usually provides just one dimension
+                if len(dim_line) == 1:
+                    n_rows = int(dim_line[0])
+                    n_cols = n_rows
+                else:
+                    n_rows, n_cols = map(int, dim_line[:2])
+                
+                hessian_matrix = np.zeros((n_rows, n_cols))
+                col_pointer = 0
+                while col_pointer < n_cols:
+                    header = next(iterator).strip()
+                    if not header: break
+                    col_indices = [int(c) for c in header.split()]
+                    for r in range(n_rows):
+                        row_data = next(iterator).strip().split()
+                        row_idx = int(row_data[0])
+                        values = [float(x) for x in row_data[1:]]
+                        for i, val in enumerate(values):
+                            hessian_matrix[row_idx, col_indices[i]] = val
+                    col_pointer += len(col_indices)
+                break
+
+        if hessian_matrix is None:
+            raise ValueError("Could not find $hessian block in the file.")
+        return hessian_matrix