Note
Go to the end to download the full example code.
Finding Reaction Paths with EON and a Metatomic Potential¶
- Authors:
Rohit Goswami @HaoZeke, Hanna Tuerk @HannaTuerk, Arslan Mazitov @abmazitov, Michele Ceriotti @ceriottim
This example describes how to find the reaction pathway for oxadiazole formation from N₂O and ethylene. We will use the PET-MAD metatomic model to calculate the potential energy and forces.
The primary goal is to contrast a standard Nudged Elastic Band (NEB) calculation using the atomic simulation environment (ASE) with more sophisticated methods available in the EON package. For even a relatively simple reaction like this, a basic NEB implementation can struggle to converge or may time out. We will show how EON’s advanced features, such as energy-weighted springs and mixing in single-ended dimer search steps, can efficiently locate and refine the transition state along the path.
Our approach will be:
Set up the PET-MAD metatomic calculator.
Use ASE to generate an initial IDPP reaction path.
Illustrate the limitations of a standard NEB calculation in ASE.
Refine the path and locate the transition state saddle point using EON’s optimizers, including energy-weighted springs and the dimer method.
Visualize the final converged pathway.
Demonstrate endpoind relaxation with EON
Importing Required Packages¶
First, we import all the necessary python packages for this task.
By convention, all import
statements are at the top of the file.
import contextlib
import os
import shutil
import subprocess
import sys
from pathlib import Path
from typing import List, Optional, Union
import ase
import ase.io as aseio
import ira_mod
import matplotlib.image as mpimg
import matplotlib.pyplot as plt
import requests
from ase.mep import NEB
from ase.optimize import LBFGS
from ase.visualize import view
from ase.visualize.plot import plot_atoms
from metatomic.torch.ase_calculator import MetatomicCalculator
# sphinx_gallery_thumbnail_number = 4
Obtaining the Foundation Model - PET-MAD¶
PET-MAD is an instance of a point edge transformer model trained on
the diverse MAD dataset
which learns equivariance through data driven measures
instead of having equivariance baked in [1]. In turn, this enables
the PET model to have greater design space to learn over. Integration in
Python and the C++ EON client occurs through the metatomic software [2],
which in turn relies on the atomistic machine learning toolkit build
over metatensor. Essentially using any of the metatomic models involves
grabbing weights off of HuggingFace and loading them with
metatomic before running the
engine of choice.
repo_id = "lab-cosmo/pet-mad"
tag = "v1.0.2"
url_path = f"models/pet-mad-{tag}.ckpt"
fname = Path(url_path).name
url = f"https://huggingface.co/{repo_id}/resolve/{tag}/{url_path}"
if not Path(fname).exists():
response = requests.get(url)
if response.ok:
with open(fname, "wb") as f:
f.write(response.content)
print(f"Downloaded {fname} from tag {tag}.")
else:
print("Failed to download:", response.status_code)
else:
print(f"{fname} from tag {tag} already present.")
if not Path(fname.replace("ckpt", "pt")).exists():
subprocess.run(
[
"mtt",
"export",
"pet-mad-v1.0.2.ckpt", # noqa: E501
]
)
Downloaded pet-mad-v1.0.2.ckpt from tag v1.0.2.
Nudged Elastic Band (NEB)¶
Given two known configurations on a potential energy surface (PES), often one wishes to determine the path of highest probability between the two. Under the harmonic approximation to transition state theory, connecting the configurations (each point representing a full molecular structure) by a discrete set of images allows one to evolve the path under an optimization algorithm, and allows approximating the reaction to three states: the reactant, product, and transition state.
The location of this transition state (≈ the point with the highest energy along this path) determines the barrier height of the reaction. This saddle point can be found by transforming the second derivatives (Hessian) to step along the softest mode. However, an approximation which is free from explicitly finding this mode involves moving the highest image of a NEB path: the “climbing” image.
Mathematically, the saddle point has zero first derivatives and a single negative eigenvalue. The climbing image technique moves the highest energy image along the reversed NEB tangent force, avoiding the cost of full Hessian diagonalization used in single-ended methods [3].
Concretely, in this example, we will consider a reactant and product state, for oxadiazole formation, namely N₂O and ethylene.
reactant = aseio.read("data/min_reactant.con")
product = aseio.read("data/min_product.con")
We can visualize these structures using ASE.
fig, (ax1, ax2) = plt.subplots(1, 2)
plot_atoms(reactant, ax1, rotation=("-90x,0y,0z"))
plot_atoms(product, ax2, rotation=("-90x,0y,0z"))
ax1.text(0.3, -1, "reactant")
ax2.text(0.3, -1, "product")
ax1.set_axis_off()
ax2.set_axis_off()
Endpoint minimization¶
For finding reaction pathways, the endpoints should be minimized. We provide initial configurations which are already minimized, but in order to see how to relax endpoints with EON, please have a look at the end of this tutorial.
Initial path generation¶
To begin an NEB method, an initial path is required, the optimal construction of which still forms an active area of research. The simplest approximation to an initial path for NEB methods linearly interpolate between the two known configurations building on intuition developed from “drag coordinate” methods. This may break bonds or otherwise also unphysically pass atoms through each other, similar to the effect of incorrect permutations. To ameliorate this effect, the NEB algorithm is often started from the linear interpolation but then the path is optimized on a surrogate potential energy surface, commonly something cheap and analytic, like the IDPP (Image dependent pair potential, [5]) which provides a surface based on bond distances, and thus preventing atom-in-atom collisions.
Here we use the IDPP from ASE to setup the initial path. You can find more information about this method in the corresponding ASE tutorial or in the original IDPP publication by S. Smidstru et al. . A brief pedagogical discussion of the transition state methods may be found on the Rowan blog, though the software is proprietary there.
N_INTERMEDIATE_IMGS = 10
# total includes the endpoints
TOTAL_IMGS = N_INTERMEDIATE_IMGS + 2
images = [reactant]
images += [reactant.copy() for img in range(N_INTERMEDIATE_IMGS)]
images += [product]
neb = NEB(images)
neb.interpolate("idpp")
We don’t cover subtleties in setting the number of images, typically too many intermediate images may cause kinks but too few will be unable to resolve the tangent to any reasonable quality.
For EON, we write the initial path to a file called idppPath.dat.
output_dir = "path"
os.makedirs(output_dir, exist_ok=True)
output_files = [f"{output_dir}/{num:02d}.con" for num in range(TOTAL_IMGS)]
for outfile, img in zip(output_files, images):
ase.io.write(outfile, img)
print(f"Wrote {len(output_files)} IDPP images to '{output_dir}/'.")
summary_file_path = "idppPath.dat"
with open(summary_file_path, "w") as f:
for filepath in output_files:
abs_path = os.path.abspath(filepath)
f.write(f"{abs_path}\n")
print(f"Wrote absolute paths to '{summary_file_path}'.")
Wrote 12 IDPP images to 'path/'.
Wrote absolute paths to 'idppPath.dat'.
Running NEBs¶
We will now consider actually running the Nudged Elastic Band with different codes.
ASE and Metatomic¶
We first consider using metatomic with the ASE calculator.
# define the calculator
def mk_mta_calc():
return MetatomicCalculator(
"pet-mad-v1.0.2.pt",
device="cpu",
non_conservative=False,
uncertainty_threshold=0.001,
)
# set calculators for images
ipath = [reactant] + [reactant.copy() for img in range(10)] + [product]
for img in ipath:
img.calc = mk_mta_calc()
print(img.calc._model.capabilities().outputs)
neb = NEB(ipath, climb=True, k=5, method="improvedtangent")
neb.interpolate("idpp")
# store initial path guess for plotting
initial_energies = [img.get_potential_energy() for img in ipath]
# setup the NEB clalculation
optimizer = LBFGS(neb, trajectory="A2B.traj", logfile="opt.log")
conv = optimizer.run(fmax=0.01, steps=100)
print("Check if calculation has converged:", conv)
if conv:
print(neb)
final_energies = [i.get_potential_energy() for i in ipath]
# Plot initial and final path
plt.figure(figsize=(8, 6))
# Initial Path (Blue)
plt.plot(
initial_energies - initial_energies[0],
"o-",
label="Initial Path (IDPP)",
color="xkcd:blue",
)
# Final Path (Orange)
plt.plot(
final_energies - initial_energies[0],
"o-",
label="Optimized Path (LBFGS)",
color="xkcd:orange",
)
# Metadata
plt.xlabel("Image number")
plt.ylabel("Potential Energy (eV)")
plt.legend()
plt.grid(True, alpha=0.3)
plt.title("NEB Path Evolution")
plt.show()
{'energy': <torch.ScriptObject object at 0x5595462d0ed0>, 'energy_ensemble': <torch.ScriptObject object at 0x559549eb19b0>, 'energy_uncertainty': <torch.ScriptObject object at 0x559549eb1780>, 'features': <torch.ScriptObject object at 0x559547d32930>, 'mtt::aux::energy_last_layer_features': <torch.ScriptObject object at 0x559548f00020>}
/home/runner/work/atomistic-cookbook/atomistic-cookbook/.nox/eon-pet-neb/lib/python3.13/site-packages/ase/calculators/calculator.py:519: UserWarning: Some of the atomic energy uncertainties are larger than the threshold of 0.001 eV. The prediction is above the threshold for atoms [0 1 2 3 4 5 6 7 8].
self.calculate(atoms, [name], system_changes)
Check if calculation has converged: False
In the 100 NEB steps we took, the structure did unfortunately not converge. The metatomic calculator for PET-MAD v1.0.2 provides LLPR based energy uncertainities. As we obtain a warning that the uncertainty of the path structure sampled is very high, we stop after 100 steps. The ASE algorithm with LBFGS optimizer does not find good intermediate structures and does not converge at all. Our test showed that the FIRE optimizer works better in this context, but still takes over 500 steps to converge, and since second order methods are faster, we consider the LBFGS routine throughout this notebook.
We thus want to look at a different code, which manages to compute a NEB for this simple system more efficiently.
EON and Metatomic¶
EON has two improvements to accurately locate the saddle point.
Energy weighting for improving tangent resolution near the climbing image
The Hybrid MMF-NEB-CI which involves iteratively switching to the dimer method for faster convergence by the climbing image.
To use EON, we setup a function that writes the desired EON input for us and
runs the eonclient binary. Since we are in a notebook environment, we will
use several abstractions over raw subprocess calls. In practice, writing
and using EON is much simpler.
def write_neb_eon_config(
run_dir: Path, model_path: Path, ninterm: int = N_INTERMEDIATE_IMGS
):
"""
Writes the config.ini file for an EON NEB job,
using the user's provided template.
"""
config_content = f"""[Main]
job=nudged_elastic_band
random_seed = 706253457
[Potential]
potential = Metatomic
[Metatomic]
model_path = {Path(str(model_path).replace("ckpt", "pt")).absolute()}
[Nudged Elastic Band]
images={ninterm}
energy_weighted=true
ew_ksp_min = 0.972
ew_ksp_max = 9.72
initial_path_in = idppPath.dat
minimize_endpoints = false
climbing_image_method = true
climbing_image_converged_only = true
ci_after = 1.5
ci_mmf = true
ci_mmf_after = 0.5
ci_mmf_nsteps = 20
# for PET-MAD 1.1
# ci_after = 0.5
# ci_mmf_nsteps = 10
[Optimizer]
max_iterations = 100
opt_method = lbfgs
max_move = 0.1
converged_force = 0.01
[Debug]
write_movies=true
"""
config_path = run_dir / "config.ini"
with open(config_path, "w") as f:
f.write(config_content)
print(f"Wrote EON NEB config to '{config_path}'")
Which now let’s us write out the final triplet of reactant, product, and configuration of the EON-NEB.
write_neb_eon_config(Path("."), Path(fname).absolute())
aseio.write("reactant.con", reactant)
aseio.write("product.con", product)
Wrote EON NEB config to 'config.ini'
Now we can finally define helpers for easier visualization.
def _run_command_live(
cmd: Union[str, List[str]],
*,
check: bool = True,
timeout: Optional[float] = None,
capture: bool = False,
encoding: str = "utf-8",
) -> subprocess.CompletedProcess:
"""
Internal: run command and stream stdout/stderr live to current stdout.
If capture=True, also collect combined output and return it
in CompletedProcess.stdout.
"""
shell = isinstance(cmd, str)
cmd_str = cmd if shell else cmd[0]
# If list form, ensure program exists before trying to run
if not shell and shutil.which(cmd_str) is None:
raise FileNotFoundError(f"{cmd_str!r} is not on PATH")
# Start the process
# We combine stderr into stdout so we only have one stream to read
proc = subprocess.Popen(
cmd,
stdout=subprocess.PIPE,
stderr=subprocess.STDOUT,
text=True,
encoding=encoding,
shell=shell,
bufsize=1, # Line buffered
)
collected = [] if capture else None
try:
assert proc.stdout is not None
for line in proc.stdout:
# Stream into notebook or terminal live
print(line, end="")
sys.stdout.flush()
if capture:
collected.append(line)
# Wait for the process to actually exit after stream closes
returncode = proc.wait(timeout=timeout)
except subprocess.TimeoutExpired:
proc.kill()
proc.wait()
raise
except KeyboardInterrupt:
proc.kill()
proc.wait()
raise
finally:
if proc.stdout:
proc.stdout.close()
if check and returncode != 0:
output_str = "".join(collected) if capture else ""
raise subprocess.CalledProcessError(returncode, cmd, output=output_str)
return subprocess.CompletedProcess(
cmd, returncode, stdout="".join(collected) if capture else None
)
def run_command_or_exit(
cmd: Union[str, List[str]], capture: bool = False, timeout: Optional[float] = 300
) -> subprocess.CompletedProcess:
"""
Helper wrapper to run commands, stream output, and exit script/notebook
cleanly on failure so sphinx-gallery sees the errors appropriately.
"""
try:
return _run_command_live(cmd, check=True, capture=capture, timeout=timeout)
except FileNotFoundError as e:
print(f"Executable not found: {e}", file=sys.stderr)
sys.exit(2)
except subprocess.CalledProcessError as e:
print(f"Command failed with exit code {e.returncode}", file=sys.stderr)
sys.exit(e.returncode)
except subprocess.TimeoutExpired:
print("Command timed out", file=sys.stderr)
sys.exit(124)
Run the main C++ client¶
This runs ‘eonclient’ and streams output live. If it fails, the notebook execution stops here.
run_command_or_exit(["eonclient"], capture=True, timeout=300)
EON Client
VERSION: f01406f
BUILD DATE: Tue Dec 02 01:53:51 AM GMT 2025
OS: linux
Arch: x86_64
Hostname: runnervmoqczp
PID: 2507
DIR: /home/runner/work/atomistic-cookbook/atomistic-cookbook/examples/eon-pet-neb
Loading parameter file config.ini
* [Main] job: nudged_elastic_band
* [Main] random_seed: 706253457
* [Potential] potential: Metatomic
* [Debug] write_movies: true
* [Optimizer] opt_method: lbfgs
* [Optimizer] converged_force: 0.01
* [Optimizer] max_iterations: 100
* [Optimizer] max_move: 0.1
* [Metatomic] model_path: /home/runner/work/atomistic-cookbook/atomistic-cookbook/examples/eon-pet-neb/pet-mad-v1.0.2.pt
* [Nudged Elastic Band] images: 10
* [Nudged Elastic Band] climbing_image_method: true
* [Nudged Elastic Band] climbing_image_converged_only: true
* [Nudged Elastic Band] energy_weighted: true
* [Nudged Elastic Band] ew_ksp_min: 0.972
* [Nudged Elastic Band] ew_ksp_max: 9.72
* [Nudged Elastic Band] initial_path_in: idppPath.dat
* [Nudged Elastic Band] minimize_endpoints: false
* [Nudged Elastic Band] ci_after: 1.5
* [Nudged Elastic Band] ci_mmf: true
* [Nudged Elastic Band] ci_mmf_after: 0.5
* [Nudged Elastic Band] ci_mmf_nsteps: 20
nebMinimEP == false: not minimizing endpoints.
Nudged elastic band calculation started.
iteration step size ||Force|| max image max energy
---------------------------------------------------------------
1 0.0000e+00 6.4530e+00 7 1.476
2 2.7135e-02 3.4314e+00 7 1.313
3 2.3174e-02 2.9734e+00 7 1.246
4 7.9659e-02 4.1153e+00 7 1.144
5 4.2119e-02 3.4505e+00 7 1.055
6 1.3126e-02 2.5218e+00 8 1.027
7 4.3104e-02 2.1791e+00 8 1.007
8 6.5064e-02 2.2961e+00 8 0.9603
9 7.1163e-02 2.1858e+00 8 0.9309
10 1.3802e-02 9.2779e-01 8 0.9
11 1.0319e-02 9.8000e-01 8 0.8936
12 1.8940e-02 1.3644e+00 8 0.8813
13 5.1184e-02 1.7403e+00 8 0.8499
14 3.3102e-02 1.0731e+00 8 0.8264
15 5.3943e-03 8.5662e-01 8 0.8209
16 1.4975e-02 9.7261e-01 8 0.8166
17 1.7579e-02 1.0942e+00 8 0.8065
18 4.1575e-02 1.0257e+00 8 0.7855
19 6.1904e-02 8.6704e-01 8 0.7693
20 2.4118e-02 6.7097e-01 8 0.7669
21 1.3001e-02 6.1148e-01 8 0.7647
22 3.0034e-02 5.7133e-01 8 0.7647
23 4.7631e-02 8.4101e-01 8 0.7684
24 4.1914e-02 1.7999e+00 8 0.7753
25 2.7705e-02 7.0061e-01 8 0.7605
26 8.5950e-03 4.6628e-01 8 0.7593
Saddle point search started from reactant with energy -53.27388000488281 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.7511 3.792 2.39 1
[Dimer] 1 0.0186833 -0.0015 8.88839e-01 -17.7511 3.792 2.393 1
[IDimerRot] ----- --------- ---------- ------------------ -17.8759 6.108 ------ ----
[Dimer] 2 0.0165101 -0.0059 4.08662e-01 -17.8759 3.054 3.674 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7504 3.630 2.10 1
[Dimer] 3 0.0139743 -0.0087 2.10910e-01 -17.7504 3.630 2.101 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6027 4.685 ------ ----
[Dimer] 4 0.0088965 -0.0097 2.04202e-01 -17.6027 2.342 3.472 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6127 4.699 ------ ----
[Dimer] 5 0.0186002 -0.0098 3.14657e-01 -17.6127 2.350 3.799 0
[IDimerRot] ----- --------- ---------- ------------------ -18.1543 4.708 ------ ----
[Dimer] 6 0.0065046 -0.0100 1.71341e-01 -18.1543 2.354 3.694 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8074 4.594 ------ ----
[Dimer] 7 0.0040005 -0.0102 4.36882e-02 -17.8074 2.297 3.675 0
[IDimerRot] ----- --------- ---------- ------------------ -17.9121 4.537 ------ ----
[Dimer] 8 0.0007495 -0.0102 3.45538e-02 -17.9121 2.269 3.609 0
[IDimerRot] ----- --------- ---------- ------------------ -17.9171 4.530 ------ ----
[Dimer] 9 0.0025785 -0.0103 4.86666e-02 -17.9171 2.265 3.603 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8873 4.550 ------ ----
[Dimer] 10 0.0033890 -0.0104 5.63030e-02 -17.8873 2.275 3.624 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8446 4.515 ------ ----
[Dimer] 11 0.0074655 -0.0104 7.51035e-02 -17.8446 2.257 3.605 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7884 4.327 ------ ----
[Dimer] 12 0.0045799 -0.0104 9.79015e-02 -17.7884 2.163 3.467 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7174 4.251 ------ ----
[Dimer] 13 0.0026098 -0.0105 1.99273e-02 -17.7174 2.125 3.420 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7402 4.292 ------ ----
[Dimer] 14 0.0003734 -0.0105 1.47679e-02 -17.7402 2.146 3.449 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7279 4.278 ------ ----
[Dimer] 15 0.0015581 -0.0105 2.75431e-02 -17.7279 2.139 3.440 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6926 4.253 ------ ----
[Dimer] 16 0.0023492 -0.0105 4.26412e-02 -17.6926 2.127 3.427 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6535 4.207 ------ ----
[Dimer] 17 0.0066182 -0.0105 6.18012e-02 -17.6535 2.104 3.398 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5841 4.114 ------ ----
[Dimer] 18 0.0106054 -0.0106 6.60917e-02 -17.5841 2.057 3.336 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5231 4.004 ------ ----
[Dimer] 19 0.0285095 -0.0107 1.78217e-01 -17.5231 2.002 3.259 0
[IDimerRot] ----- --------- ---------- ------------------ -17.3765 3.918 ------ ----
[Dimer] 20 0.0082968 -0.0108 7.05261e-02 -17.3765 1.959 3.216 0
27 5.0021e-02 6.6456e-02 8 0.7485
Saddle point search started from reactant with energy -53.28470230102539 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.4712 3.807 2.30 1
[Dimer] 1 0.0005027 -0.0000 3.91019e-02 -17.4712 3.807 2.300 1
[IDimerRot] ----- --------- ---------- ------------------ -17.4542 3.656 ------ ----
[Dimer] 2 0.0002679 -0.0000 3.66216e-02 -17.4542 1.828 2.688 0
[IDimerRot] ----- --------- ---------- ------------------ -17.4563 3.655 ------ ----
[Dimer] 3 0.0049187 -0.0001 4.55772e-02 -17.4563 1.827 2.988 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5819 3.524 ------ ----
[Dimer] 4 0.0045667 -0.0002 4.99926e-02 -17.5819 1.762 2.862 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6530 3.657 ------ ----
[Dimer] 5 0.0033963 -0.0002 4.24197e-02 -17.6530 1.828 2.956 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7338 3.784 ------ ----
[Dimer] 6 0.0006650 -0.0002 2.52607e-02 -17.7338 1.892 3.045 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7051 3.848 ------ ----
[Dimer] 7 0.0018353 -0.0002 2.15064e-02 -17.7051 1.924 3.101 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6821 4.038 ------ ----
[Dimer] 8 0.0016934 -0.0003 2.36644e-02 -17.6821 2.019 3.257 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7007 4.146 ------ ----
[Dimer] 9 0.0065684 -0.0003 5.61226e-02 -17.7007 2.073 3.340 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7781 4.266 ------ ----
[Dimer] 10 0.0040765 -0.0003 2.82622e-02 -17.7781 2.133 3.421 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8319 4.239 ------ ----
[Dimer] 11 0.0012096 -0.0003 2.50769e-02 -17.8319 2.120 3.389 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8399 4.108 ------ ----
[Dimer] 12 0.0056081 -0.0004 3.33256e-02 -17.8399 2.054 3.284 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8419 3.860 ------ ----
[Dimer] 13 0.0119066 -0.0005 5.63237e-02 -17.8419 1.930 3.087 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8616 3.739 ------ ----
[Dimer] 14 0.0322391 -0.0007 9.39375e-02 -17.8616 1.869 2.988 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8062 3.640 ------ ----
[Dimer] 15 0.0189463 -0.0008 5.57000e-02 -17.8062 1.820 2.918 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8234 3.894 ------ ----
[Dimer] 16 0.0238557 -0.0009 4.64804e-02 -17.8234 1.947 3.117 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7842 4.410 ------ ----
[Dimer] 17 0.0242690 -0.0010 5.36686e-02 -17.7842 2.205 3.534 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7350 4.851 ------ ----
[Dimer] 18 0.0462534 -0.0012 7.98999e-02 -17.7350 2.426 3.894 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6298 5.583 ------ ----
[Dimer] 19 0.0630176 -0.0014 1.14882e-01 -17.6298 2.791 4.499 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5742 3.186 1.73 1
[Dimer] 20 0.0572431 -0.0014 1.87180e-01 -17.5742 3.186 1.728 1
28 1.5729e-01 6.0795e-01 8 0.7487
29 9.4161e-02 8.9647e-01 8 0.7518
30 3.4622e-02 4.2468e-01 8 0.7472
Saddle point search started from reactant with energy -53.28594970703125 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.0447 4.583 5.30 1
[IDimerRot] ----- --------- ---------- ------------------ -17.3874 5.129 5.75 2
[IDimerRot] ----- --------- ---------- ------------------ -17.3873 6.163 ------ ----
[Dimer] 1 0.0030546 -0.0006 1.72727e-01 -17.3873 3.082 3.914 2
[IDimerRot] ----- --------- ---------- ------------------ -17.5742 3.301 1.29 1
[Dimer] 2 0.0011993 -0.0007 1.00714e-01 -17.5742 3.301 1.292 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6405 2.702 ------ ----
[Dimer] 3 0.0020557 -0.0007 5.80249e-02 -17.6405 1.351 2.015 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6746 2.792 ------ ----
[Dimer] 4 0.0015384 -0.0008 5.80295e-02 -17.6746 1.396 2.258 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6769 2.866 ------ ----
[Dimer] 5 0.0052372 -0.0009 6.58701e-02 -17.6769 1.433 2.317 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6520 3.231 ------ ----
[Dimer] 6 0.0086470 -0.0012 1.44669e-01 -17.6520 1.616 2.615 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6546 3.442 ------ ----
[Dimer] 7 0.0039242 -0.0014 8.83451e-02 -17.6546 1.721 2.784 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6483 3.471 ------ ----
[Dimer] 8 0.0139955 -0.0016 9.49049e-02 -17.6483 1.736 2.808 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6909 3.258 ------ ----
[Dimer] 9 0.0153976 -0.0018 9.56474e-02 -17.6909 1.629 2.631 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7818 2.999 ------ ----
[Dimer] 10 0.0191170 -0.0020 9.31785e-02 -17.7818 1.500 2.410 0
[IDimerRot] ----- --------- ---------- ------------------ -17.9805 3.383 ------ ----
[Dimer] 11 0.0115168 -0.0020 9.20679e-02 -17.9805 1.691 2.687 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8820 3.168 ------ ----
[Dimer] 12 0.0055241 -0.0021 3.88238e-02 -17.8820 1.584 2.531 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8995 3.175 ------ ----
[Dimer] 13 0.0017167 -0.0021 3.32521e-02 -17.8995 1.587 2.534 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8821 3.183 ------ ----
[Dimer] 14 0.0046314 -0.0022 4.48691e-02 -17.8821 1.592 2.543 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8063 3.192 ------ ----
[Dimer] 15 0.0062909 -0.0022 5.80797e-02 -17.8063 1.596 2.561 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7693 3.206 ------ ----
[Dimer] 16 0.0056639 -0.0022 2.46150e-02 -17.7693 1.603 2.577 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7177 3.252 ------ ----
[Dimer] 17 0.0010379 -0.0023 1.72468e-02 -17.7177 1.626 2.622 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7322 3.253 ------ ----
[Dimer] 18 0.0014562 -0.0023 2.34185e-02 -17.7322 1.626 2.620 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7488 3.267 ------ ----
[Dimer] 19 0.0027184 -0.0023 3.51576e-02 -17.7488 1.633 2.629 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7450 3.273 ------ ----
[Dimer] 20 0.0062004 -0.0023 4.40626e-02 -17.7450 1.637 2.635 0
31 1.3745e-01 1.0119e+00 8 0.7531
32 2.1750e-02 4.5098e-01 8 0.7482
Saddle point search started from reactant with energy -53.28496170043945 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -16.7645 7.254 5.35 1
[IDimerRot] ----- --------- ---------- ------------------ -17.5518 7.131 8.78 2
[IDimerRot] ----- --------- ---------- ------------------ -17.7578 4.072 3.54 3
[Dimer] 1 0.0058540 -0.0013 2.25586e-01 -17.7578 4.072 3.541 3
[IDimerRot] ----- --------- ---------- ------------------ -17.9623 3.833 ------ ----
[Dimer] 2 0.0026683 -0.0014 1.48721e-01 -17.9623 1.917 2.529 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8766 3.661 ------ ----
[Dimer] 3 0.0007327 -0.0015 8.05034e-02 -17.8766 1.831 2.923 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8397 3.816 ------ ----
[Dimer] 4 0.0015095 -0.0016 7.74910e-02 -17.8397 1.908 3.053 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7567 4.121 ------ ----
[Dimer] 5 0.0026335 -0.0017 1.10506e-01 -17.7567 2.061 3.310 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6588 4.561 ------ ----
[Dimer] 6 0.0061106 -0.0019 1.41959e-01 -17.6588 2.281 3.680 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5304 5.230 ------ ----
[Dimer] 7 0.0113424 -0.0024 1.52377e-01 -17.5304 2.615 4.242 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5142 5.383 ------ ----
[Dimer] 8 0.0269133 -0.0023 3.82310e-01 -17.5142 2.692 4.368 0
[IDimerRot] ----- --------- ---------- ------------------ -17.2801 4.656 ------ ----
[Dimer] 9 0.0037183 -0.0029 1.01693e-01 -17.2801 2.328 3.836 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6056 4.172 ------ ----
[Dimer] 10 0.0023577 -0.0030 6.89033e-02 -17.6056 2.086 3.378 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6890 3.576 ------ ----
[Dimer] 11 0.0065721 -0.0032 9.22760e-02 -17.6890 1.788 2.886 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7212 2.918 ------ ----
[Dimer] 12 0.0112717 -0.0035 1.06833e-01 -17.7212 1.459 2.353 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6557 2.726 ------ ----
[Dimer] 13 0.0286049 -0.0039 1.31719e-01 -17.6557 1.363 2.207 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7469 3.194 ------ ----
[Dimer] 14 0.0172980 -0.0038 1.68076e-01 -17.7469 1.597 2.571 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5554 4.121 ------ ----
[Dimer] 15 0.0089408 -0.0040 3.02141e-02 -17.5554 2.061 3.348 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6044 4.080 ------ ----
[Dimer] 16 0.0006806 -0.0040 1.96092e-02 -17.6044 2.040 3.305 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6147 4.240 ------ ----
[Dimer] 17 0.0012700 -0.0040 2.32826e-02 -17.6147 2.120 3.432 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6142 4.424 ------ ----
[Dimer] 18 0.0022913 -0.0041 3.38714e-02 -17.6142 2.212 3.579 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6325 4.479 ------ ----
[Dimer] 19 0.0050936 -0.0041 3.19423e-02 -17.6325 2.240 3.620 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6296 4.494 ------ ----
[Dimer] 20 0.0043822 -0.0041 4.93339e-02 -17.6296 2.247 3.632 0
33 6.1520e-02 2.2993e-01 8 0.7445
Saddle point search started from reactant with energy -53.28866958618164 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.0445 6.393 5.36 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6203 6.506 7.18 2
[IDimerRot] ----- --------- ---------- ------------------ -17.6199 7.806 ------ ----
[Dimer] 1 0.0022962 -0.0003 1.05137e-01 -17.6199 3.903 4.938 2
[IDimerRot] ----- --------- ---------- ------------------ -17.6285 3.905 2.13 1
[Dimer] 2 0.0009249 -0.0003 3.91022e-02 -17.6285 3.905 2.127 1
[IDimerRot] ----- --------- ---------- ------------------ -17.5999 3.697 ------ ----
[Dimer] 3 0.0003285 -0.0003 2.75477e-02 -17.5999 1.849 2.725 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6220 3.682 ------ ----
[Dimer] 4 0.0009295 -0.0003 2.69741e-02 -17.6220 1.841 2.982 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6620 3.667 ------ ----
[Dimer] 5 0.0013890 -0.0003 3.58414e-02 -17.6620 1.833 2.963 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6946 3.665 ------ ----
[Dimer] 6 0.0036469 -0.0004 4.18327e-02 -17.6946 1.833 2.956 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7231 3.583 ------ ----
[Dimer] 7 0.0071206 -0.0004 1.06688e-01 -17.7231 1.791 2.886 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7671 3.708 ------ ----
[Dimer] 8 0.0013238 -0.0005 4.35987e-02 -17.7671 1.854 2.979 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7250 3.647 ------ ----
[Dimer] 9 0.0022879 -0.0005 3.16250e-02 -17.7250 1.824 2.937 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6740 3.689 ------ ----
[Dimer] 10 0.0033851 -0.0005 4.36517e-02 -17.6740 1.845 2.979 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6409 3.779 ------ ----
[Dimer] 11 0.0065942 -0.0006 4.57707e-02 -17.6409 1.889 3.057 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6453 3.872 ------ ----
[Dimer] 12 0.0192952 -0.0008 6.23592e-02 -17.6453 1.936 3.130 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5417 3.633 ------ ----
[Dimer] 13 0.0445432 -0.0006 1.86832e-01 -17.5417 1.817 2.956 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5289 3.690 ------ ----
[Dimer] 14 0.0209166 -0.0008 3.02960e-02 -17.5289 1.845 3.004 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6110 3.648 ------ ----
[Dimer] 15 0.0029158 -0.0008 2.29262e-02 -17.6110 1.824 2.956 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6474 3.719 ------ ----
[Dimer] 16 0.0048682 -0.0009 2.55737e-02 -17.6474 1.859 3.007 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7013 3.798 ------ ----
[Dimer] 17 0.0040658 -0.0009 5.17058e-02 -17.7013 1.899 3.062 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7138 3.857 ------ ----
[Dimer] 18 0.0047134 -0.0009 2.57320e-02 -17.7138 1.929 3.107 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6958 3.828 ------ ----
[Dimer] 19 0.0182250 -0.0010 3.55781e-02 -17.6958 1.914 3.086 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6217 3.842 ------ ----
[Dimer] 20 0.0096075 -0.0010 4.28253e-02 -17.6217 1.921 3.111 0
34 6.3226e-02 3.6384e-02 8 0.7435
Saddle point search started from reactant with energy -53.2896728515625 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -16.5917 6.852 8.10 1
[IDimerRot] ----- --------- ---------- ------------------ -17.4231 8.292 8.89 2
[IDimerRot] ----- --------- ---------- ------------------ -17.6258 4.408 3.23 3
[Dimer] 1 0.0010842 -0.0000 4.23503e-02 -17.6258 4.408 3.227 3
[IDimerRot] ----- --------- ---------- ------------------ -17.6381 3.583 ------ ----
[Dimer] 2 0.0008253 -0.0000 3.89795e-02 -17.6381 1.791 2.610 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6489 3.637 ------ ----
[Dimer] 3 0.0108179 -0.0001 9.55080e-02 -17.6489 1.819 2.942 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7245 4.014 ------ ----
[Dimer] 4 0.0014020 -0.0002 8.80723e-02 -17.7245 2.007 3.230 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7277 3.978 ------ ----
[Dimer] 5 0.0016016 -0.0002 2.77570e-02 -17.7277 1.989 3.200 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7158 3.889 ------ ----
[Dimer] 6 0.0002632 -0.0002 2.10972e-02 -17.7158 1.944 3.132 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7143 3.836 ------ ----
[Dimer] 7 0.0027060 -0.0002 3.54648e-02 -17.7143 1.918 3.090 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7035 3.568 ------ ----
[Dimer] 8 0.0031177 -0.0003 4.59413e-02 -17.7035 1.784 2.877 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6939 3.422 ------ ----
[Dimer] 9 0.0083902 -0.0003 5.89878e-02 -17.6939 1.711 2.762 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6905 3.296 ------ ----
[Dimer] 10 0.0150948 -0.0004 9.83018e-02 -17.6905 1.648 2.661 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6460 3.198 ------ ----
[Dimer] 11 0.0032537 -0.0004 4.16833e-02 -17.6460 1.599 2.589 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6882 3.427 ------ ----
[Dimer] 12 0.0032497 -0.0005 4.00265e-02 -17.6882 1.713 2.767 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7132 3.809 ------ ----
[Dimer] 13 0.0021409 -0.0005 3.81935e-02 -17.7132 1.904 3.068 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7169 3.886 ------ ----
[Dimer] 14 0.0155190 -0.0006 1.11311e-01 -17.7169 1.943 3.130 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6627 3.759 ------ ----
[Dimer] 15 0.0092842 -0.0006 2.96548e-02 -17.6627 1.880 3.037 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7083 3.778 ------ ----
[Dimer] 16 0.0013749 -0.0006 1.84338e-02 -17.7083 1.889 3.044 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6829 3.680 ------ ----
[Dimer] 17 0.0039840 -0.0006 2.03833e-02 -17.6829 1.840 2.970 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6650 3.573 ------ ----
[Dimer] 18 0.0084938 -0.0007 4.01791e-02 -17.6650 1.787 2.887 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6044 3.540 ------ ----
[Dimer] 19 0.0114645 -0.0007 4.83708e-02 -17.6044 1.770 2.871 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5852 3.534 ------ ----
[Dimer] 20 0.0361547 -0.0009 8.19751e-02 -17.5852 1.767 2.869 0
35 6.5778e-02 7.3923e-02 8 0.7426
Saddle point search started from reactant with energy -53.290592193603516 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -16.0971 7.563 10.07 1
[IDimerRot] ----- --------- ---------- ------------------ -17.2874 9.976 11.07 2
[IDimerRot] ----- --------- ---------- ------------------ -17.5631 5.185 3.67 3
[Dimer] 1 0.0022988 -0.0000 8.68632e-02 -17.5631 5.185 3.665 3
[IDimerRot] ----- --------- ---------- ------------------ -17.7232 3.990 ------ ----
[Dimer] 2 0.0017628 -0.0001 5.43971e-02 -17.7232 1.995 2.854 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7420 4.187 ------ ----
[Dimer] 3 0.0029208 -0.0001 7.37342e-02 -17.7420 2.094 3.365 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7368 4.343 ------ ----
[Dimer] 4 0.0008361 -0.0002 4.01829e-02 -17.7368 2.171 3.490 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7522 4.415 ------ ----
[Dimer] 5 0.0004649 -0.0002 3.99705e-02 -17.7522 2.208 3.544 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7439 4.363 ------ ----
[Dimer] 6 0.0035331 -0.0002 4.84407e-02 -17.7439 2.182 3.505 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6752 4.006 ------ ----
[Dimer] 7 0.0039528 -0.0003 5.03209e-02 -17.6752 2.003 3.233 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6472 3.765 ------ ----
[Dimer] 8 0.0074953 -0.0003 6.64698e-02 -17.6472 1.883 3.045 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5914 3.413 ------ ----
[Dimer] 9 0.0019434 -0.0003 4.00183e-02 -17.5914 1.707 2.771 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6274 3.486 ------ ----
[Dimer] 10 0.0043269 -0.0004 5.44886e-02 -17.6274 1.743 2.824 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6846 3.715 ------ ----
[Dimer] 11 0.0064612 -0.0005 6.11097e-02 -17.6846 1.857 2.998 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7117 3.945 ------ ----
[Dimer] 12 0.0160476 -0.0006 6.73859e-02 -17.7117 1.972 3.177 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7116 4.178 ------ ----
[Dimer] 13 0.0216730 -0.0001 3.21701e-01 -17.7116 2.089 3.363 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5527 4.661 ------ ----
[Dimer] 14 0.0113152 -0.0007 2.41866e-02 -17.5527 2.331 3.782 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6822 4.267 ------ ----
[Dimer] 15 0.0018666 -0.0007 1.53560e-02 -17.6822 2.133 3.440 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6768 4.224 ------ ----
[Dimer] 16 0.0026477 -0.0007 1.80750e-02 -17.6768 2.112 3.407 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6718 4.174 ------ ----
[Dimer] 17 0.0031901 -0.0007 3.01389e-02 -17.6718 2.087 3.368 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6657 4.122 ------ ----
[Dimer] 18 0.0075265 -0.0007 5.44092e-02 -17.6657 2.061 3.327 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6585 4.030 ------ ----
[Dimer] 19 0.0128904 -0.0008 7.72842e-02 -17.6585 2.015 3.255 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6478 3.980 ------ ----
[Dimer] 20 0.0221173 -0.0009 1.02437e-01 -17.6478 1.990 3.217 0
36 1.0000e-01 5.6394e-01 8 0.7429
37 6.1011e-02 4.0354e+00 8 0.7926
38 3.7416e-02 1.4381e+00 8 0.7476
39 1.9157e-02 6.8016e-01 8 0.7437
40 8.6507e-02 1.2682e+00 8 0.7467
41 5.8430e-02 1.0812e+00 8 0.744
42 6.5042e-02 9.6822e-01 8 0.7431
43 1.3735e-02 4.7441e-01 8 0.7412
Saddle point search started from reactant with energy -53.291969299316406 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -16.7674 5.917 6.40 1
[IDimerRot] ----- --------- ---------- ------------------ -17.0437 4.890 4.22 2
[Dimer] 1 0.0030486 -0.0004 1.92566e-01 -17.0437 4.890 4.221 2
[IDimerRot] ----- --------- ---------- ------------------ -17.1845 6.099 ------ ----
[Dimer] 2 0.0011055 -0.0002 1.59725e-01 -17.1845 3.049 4.391 0
[IDimerRot] ----- --------- ---------- ------------------ -17.3272 3.103 1.78 1
[Dimer] 3 0.0071266 0.0003 6.09842e-02 -17.3272 3.103 1.777 1
[IDimerRot] ----- --------- ---------- ------------------ -17.5644 5.246 ------ ----
[Dimer] 4 0.0026741 0.0003 1.21334e-01 -17.5644 2.623 3.538 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6184 5.749 ------ ----
[Dimer] 5 0.0013762 0.0003 1.54177e-02 -17.6184 2.874 4.633 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5934 5.562 ------ ----
[Dimer] 6 0.0001752 0.0003 1.46934e-02 -17.5934 2.781 4.491 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5901 5.545 ------ ----
[Dimer] 7 0.0061044 0.0002 4.77200e-02 -17.5901 2.772 4.478 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5149 5.451 ------ ----
[Dimer] 8 0.0085683 0.0001 7.06336e-02 -17.5149 2.725 4.422 0
[IDimerRot] ----- --------- ---------- ------------------ -17.4694 5.420 ------ ----
[Dimer] 9 0.0195888 0.0000 8.72174e-02 -17.4694 2.710 4.409 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5022 5.156 ------ ----
[Dimer] 10 0.0235357 0.0001 1.30143e-01 -17.5022 2.578 4.190 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6635 5.302 ------ ----
[Dimer] 11 0.0093507 -0.0001 3.44241e-02 -17.6635 2.651 4.268 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6298 5.232 ------ ----
[Dimer] 12 0.0032069 -0.0001 2.09172e-02 -17.6298 2.616 4.220 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6419 5.295 ------ ----
[Dimer] 13 0.0015523 -0.0001 2.52486e-02 -17.6419 2.648 4.267 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6490 5.386 ------ ----
[Dimer] 14 0.0034716 -0.0001 3.55007e-02 -17.6490 2.693 4.337 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6491 5.477 ------ ----
[Dimer] 15 0.0100648 -0.0002 5.79383e-02 -17.6491 2.738 4.410 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6139 5.542 ------ ----
[Dimer] 16 0.0128877 -0.0002 3.31817e-02 -17.6139 2.771 4.470 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5705 5.619 ------ ----
[Dimer] 17 0.0036306 -0.0002 2.22766e-02 -17.5705 2.809 4.542 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5702 5.587 ------ ----
[Dimer] 18 0.0116286 -0.0003 4.75761e-02 -17.5702 2.793 4.517 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5400 5.483 ------ ----
[Dimer] 19 0.0139421 -0.0003 7.92040e-02 -17.5400 2.741 4.442 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5147 5.380 ------ ----
[Dimer] 20 0.0335641 -0.0004 1.34018e-01 -17.5147 2.690 4.366 0
44 7.1147e-02 1.5175e-01 8 0.7408
Saddle point search started from reactant with energy -53.29235076904297 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.2395 5.613 4.88 1
[Dimer] 1 0.0012462 -0.0001 7.60688e-02 -17.2395 5.613 4.883 1
[IDimerRot] ----- --------- ---------- ------------------ -17.4879 4.621 4.05 1
[Dimer] 2 0.0005151 -0.0001 4.42320e-02 -17.4879 4.621 4.046 1
[IDimerRot] ----- --------- ---------- ------------------ -17.4786 6.413 ------ ----
[Dimer] 3 0.0009196 -0.0001 4.30234e-02 -17.4786 3.206 4.054 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5805 3.216 1.46 1
[Dimer] 4 0.0015382 -0.0002 5.70119e-02 -17.5805 3.216 1.458 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6308 3.557 ------ ----
[Dimer] 5 0.0031501 -0.0002 5.83977e-02 -17.6308 1.778 2.504 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7143 3.594 ------ ----
[Dimer] 6 0.0059079 -0.0003 6.68993e-02 -17.7143 1.797 2.897 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6899 3.564 ------ ----
[Dimer] 7 0.0151078 -0.0002 2.16318e-01 -17.6899 1.782 2.876 0
[IDimerRot] ----- --------- ---------- ------------------ -17.9165 3.734 ------ ----
[Dimer] 8 0.0047909 -0.0004 3.89769e-02 -17.9165 1.867 2.975 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6978 3.573 ------ ----
[Dimer] 9 0.0009653 -0.0005 3.28469e-02 -17.6978 1.787 2.882 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6749 3.557 ------ ----
[Dimer] 10 0.0059571 -0.0005 5.12552e-02 -17.6749 1.778 2.873 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6555 3.509 ------ ----
[Dimer] 11 0.0067618 -0.0006 6.08301e-02 -17.6555 1.755 2.838 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6899 3.523 ------ ----
[Dimer] 12 0.0171549 -0.0006 1.49710e-01 -17.6899 1.762 2.843 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6630 3.552 ------ ----
[Dimer] 13 0.0023283 -0.0007 4.23069e-02 -17.6630 1.776 2.870 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7490 3.629 ------ ----
[Dimer] 14 0.0023298 -0.0008 2.88841e-02 -17.7490 1.814 2.918 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7600 3.660 ------ ----
[Dimer] 15 0.0020412 -0.0008 2.60318e-02 -17.7600 1.830 2.942 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7259 3.683 ------ ----
[Dimer] 16 0.0055852 -0.0008 5.83544e-02 -17.7259 1.841 2.965 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6929 3.703 ------ ----
[Dimer] 17 0.0049592 -0.0008 2.25393e-02 -17.6929 1.851 2.987 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6534 3.709 ------ ----
[Dimer] 18 0.0038781 -0.0009 1.75767e-02 -17.6534 1.854 2.998 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6421 3.703 ------ ----
[Dimer] 19 0.0043744 -0.0009 2.11259e-02 -17.6421 1.852 2.996 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6420 3.731 ------ ----
[Dimer] 20 0.0041962 -0.0009 2.51169e-02 -17.6420 1.866 3.018 0
45 5.2000e-02 2.0492e-02 8 0.7399
Saddle point search started from reactant with energy -53.293243408203125 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.3519 5.622 4.45 1
[Dimer] 1 0.0006987 -0.0000 2.71440e-02 -17.3519 5.622 4.447 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6070 4.568 4.05 1
[Dimer] 2 0.0005624 -0.0000 2.43417e-02 -17.6070 4.568 4.054 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6068 5.973 ------ ----
[Dimer] 3 0.0069810 -0.0000 1.06123e-01 -17.6068 2.986 3.860 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6459 6.136 ------ ----
[Dimer] 4 0.0007810 -0.0001 4.00513e-02 -17.6459 3.068 4.932 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6549 5.929 ------ ----
[Dimer] 5 0.0005969 -0.0001 2.51947e-02 -17.6549 2.964 4.766 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6435 5.977 ------ ----
[Dimer] 6 0.0007323 -0.0001 2.06435e-02 -17.6435 2.988 4.806 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6246 6.015 ------ ----
[Dimer] 7 0.0034097 -0.0001 4.44849e-02 -17.6246 3.008 4.842 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5941 6.098 ------ ----
[Dimer] 8 0.0064479 -0.0002 5.92689e-02 -17.5941 3.049 4.916 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5803 6.161 ------ ----
[Dimer] 9 0.0112059 -0.0002 6.30787e-02 -17.5803 3.080 4.969 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5999 6.142 ------ ----
[Dimer] 10 0.0290662 -0.0001 2.30890e-01 -17.5999 3.071 4.949 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6247 3.205 1.36 1
[Dimer] 11 0.0056721 -0.0003 4.42215e-02 -17.6247 3.205 1.363 1
[IDimerRot] ----- --------- ---------- ------------------ -17.7252 3.485 ------ ----
[Dimer] 12 0.0027816 -0.0003 2.59118e-02 -17.7252 1.742 2.358 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7401 3.615 ------ ----
[Dimer] 13 0.0030502 -0.0004 3.13948e-02 -17.7401 1.808 2.909 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7513 3.749 ------ ----
[Dimer] 14 0.0087372 -0.0004 9.63984e-02 -17.7513 1.875 3.014 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6448 3.771 ------ ----
[Dimer] 15 0.0033204 -0.0005 3.56212e-02 -17.6448 1.886 3.050 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6887 3.738 ------ ----
[Dimer] 16 0.0071304 -0.0005 2.63685e-02 -17.6887 1.869 3.015 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6902 3.657 ------ ----
[Dimer] 17 0.0133123 -0.0005 4.45056e-02 -17.6902 1.829 2.951 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6645 3.559 ------ ----
[Dimer] 18 0.0134047 -0.0005 4.90727e-02 -17.6645 1.780 2.877 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6465 3.511 ------ ----
[Dimer] 19 0.0211755 -0.0006 6.33420e-02 -17.6465 1.756 2.841 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6356 3.576 ------ ----
[Dimer] 20 0.0358747 -0.0007 1.39034e-01 -17.6356 1.788 2.895 0
46 7.9178e-02 1.3608e-01 8 0.7392
Saddle point search started from reactant with energy -53.29399490356445 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.5113 5.640 4.14 1
[Dimer] 1 0.0018230 -0.0001 5.23020e-02 -17.5113 5.640 4.136 1
[IDimerRot] ----- --------- ---------- ------------------ -17.5830 4.182 3.45 1
[Dimer] 2 0.0006216 -0.0001 4.55288e-02 -17.5830 4.182 3.447 1
[IDimerRot] ----- --------- ---------- ------------------ -17.5930 5.292 ------ ----
[Dimer] 3 0.0053472 -0.0003 5.63808e-02 -17.5930 2.646 3.492 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6254 5.330 ------ ----
[Dimer] 4 0.0051930 -0.0002 2.41292e-01 -17.6254 2.665 4.299 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6462 5.064 ------ ----
[Dimer] 5 0.0018575 -0.0003 4.12787e-02 -17.6462 2.532 4.083 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6200 5.314 ------ ----
[Dimer] 6 0.0003325 -0.0003 3.20612e-02 -17.6200 2.657 4.288 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6198 5.339 ------ ----
[Dimer] 7 0.0023856 -0.0004 5.30913e-02 -17.6198 2.669 4.308 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6292 5.401 ------ ----
[Dimer] 8 0.0024710 -0.0004 6.02418e-02 -17.6292 2.701 4.355 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6584 5.435 ------ ----
[Dimer] 9 0.0056706 -0.0005 5.61435e-02 -17.6584 2.717 4.374 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6678 5.461 ------ ----
[Dimer] 10 0.0102061 -0.0005 7.82805e-02 -17.6678 2.730 4.392 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7111 5.494 ------ ----
[Dimer] 11 0.0033309 -0.0006 3.97699e-02 -17.7111 2.747 4.408 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6964 5.488 ------ ----
[Dimer] 12 0.0024755 -0.0006 3.87575e-02 -17.6964 2.744 4.407 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6713 5.505 ------ ----
[Dimer] 13 0.0054800 -0.0007 4.10462e-02 -17.6713 2.752 4.427 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6640 5.537 ------ ----
[Dimer] 14 0.0113737 -0.0007 8.96727e-02 -17.6640 2.769 4.454 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6078 5.631 ------ ----
[Dimer] 15 0.0036000 -0.0008 3.65381e-02 -17.6078 2.816 4.543 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6405 5.617 ------ ----
[Dimer] 16 0.0048106 -0.0008 2.26379e-02 -17.6405 2.809 4.523 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6580 5.629 ------ ----
[Dimer] 17 0.0051004 -0.0008 2.21227e-02 -17.6580 2.815 4.528 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6436 5.676 ------ ----
[Dimer] 18 0.0039644 -0.0008 6.64481e-02 -17.6436 2.838 4.569 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6163 5.680 ------ ----
[Dimer] 19 0.0016364 -0.0009 1.85967e-02 -17.6163 2.840 4.579 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6188 5.666 ------ ----
[Dimer] 20 0.0006485 -0.0009 1.36509e-02 -17.6188 2.833 4.567 0
47 4.3263e-02 1.3157e-02 8 0.7383
Saddle point search started from reactant with energy -53.294864654541016 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.4148 5.636 3.95 1
[Dimer] 1 0.0019554 -0.0000 3.47009e-02 -17.4148 5.636 3.949 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6644 4.326 3.75 1
[Dimer] 2 0.0017770 -0.0000 2.92930e-02 -17.6644 4.326 3.752 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6949 5.372 ------ ----
[Dimer] 3 0.0026394 -0.0000 3.51614e-02 -17.6949 2.686 3.602 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7084 5.416 ------ ----
[Dimer] 4 0.0026682 -0.0001 3.68958e-02 -17.7084 2.708 4.347 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7031 5.417 ------ ----
[Dimer] 5 0.0086244 -0.0001 4.27433e-02 -17.7031 2.708 4.349 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6804 5.439 ------ ----
[Dimer] 6 0.0110322 -0.0001 1.09448e-01 -17.6804 2.719 4.372 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5864 5.567 ------ ----
[Dimer] 7 0.0008071 -0.0002 4.18014e-02 -17.5864 2.783 4.497 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6370 5.514 ------ ----
[Dimer] 8 0.0020352 -0.0002 2.55514e-02 -17.6370 2.757 4.442 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6720 5.491 ------ ----
[Dimer] 9 0.0053278 -0.0002 3.67761e-02 -17.6720 2.746 4.415 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7110 5.439 ------ ----
[Dimer] 10 0.0097348 -0.0003 4.02431e-02 -17.7110 2.719 4.365 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7361 5.415 ------ ----
[Dimer] 11 0.0258366 -0.0004 5.79182e-02 -17.7361 2.707 4.339 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6987 5.414 ------ ----
[Dimer] 12 0.0353756 -0.0003 1.47229e-01 -17.6987 2.707 4.348 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6992 5.413 ------ ----
[Dimer] 13 0.0134993 -0.0005 3.20373e-02 -17.6992 2.707 4.347 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6496 5.450 ------ ----
[Dimer] 14 0.0014537 -0.0005 2.42081e-02 -17.6496 2.725 4.388 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6346 5.440 ------ ----
[Dimer] 15 0.0108303 -0.0006 4.66644e-02 -17.6346 2.720 4.384 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5856 5.477 ------ ----
[Dimer] 16 0.0115998 -0.0006 5.96404e-02 -17.5856 2.738 4.425 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5753 5.494 ------ ----
[Dimer] 17 0.0170947 -0.0007 5.26898e-02 -17.5753 2.747 4.442 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5857 5.483 ------ ----
[Dimer] 18 0.0517138 -0.0008 1.03790e-01 -17.5857 2.741 4.430 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5276 5.524 ------ ----
[Dimer] 19 0.0496689 -0.0011 9.91694e-02 -17.5276 2.762 4.477 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5962 5.597 ------ ----
[Dimer] 20 0.1148644 -0.0017 1.64317e-01 -17.5962 2.798 4.518 0
48 1.6933e-01 3.6113e-01 8 0.7367
Saddle point search started from reactant with energy -53.29645538330078 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -15.7596 8.091 9.18 1
[IDimerRot] ----- --------- ---------- ------------------ -16.9228 9.534 10.60 2
[IDimerRot] ----- --------- ---------- ------------------ -17.3157 5.939 4.88 3
[Dimer] 1 0.0033396 -0.0006 1.28283e-01 -17.3157 5.939 4.877 3
[IDimerRot] ----- --------- ---------- ------------------ -17.5357 5.443 ------ ----
[Dimer] 2 0.0011015 -0.0007 6.48642e-02 -17.5357 2.721 3.554 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5415 5.334 ------ ----
[Dimer] 3 0.0005531 -0.0007 5.51835e-02 -17.5415 2.667 4.322 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5328 5.335 ------ ----
[Dimer] 4 0.0033139 -0.0008 5.65086e-02 -17.5328 2.668 4.326 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5509 5.323 ------ ----
[Dimer] 5 0.0037999 -0.0009 5.34274e-02 -17.5509 2.661 4.311 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5802 5.316 ------ ----
[Dimer] 6 0.0059193 -0.0009 1.07179e-01 -17.5802 2.658 4.299 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6652 5.420 ------ ----
[Dimer] 7 0.0011464 -0.0010 3.64614e-02 -17.6652 2.710 4.361 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6356 5.413 ------ ----
[Dimer] 8 0.0010649 -0.0010 2.98751e-02 -17.6356 2.706 4.362 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6246 5.408 ------ ----
[Dimer] 9 0.0035534 -0.0010 3.65925e-02 -17.6246 2.704 4.361 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6315 5.411 ------ ----
[Dimer] 10 0.0044505 -0.0011 3.52658e-02 -17.6315 2.705 4.362 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6772 5.379 ------ ----
[Dimer] 11 0.0106824 -0.0011 7.73657e-02 -17.6772 2.689 4.325 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7067 5.392 ------ ----
[Dimer] 12 0.0044664 -0.0012 3.21949e-02 -17.7067 2.696 4.328 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7202 5.362 ------ ----
[Dimer] 13 0.0021347 -0.0012 3.23380e-02 -17.7202 2.681 4.302 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7119 5.326 ------ ----
[Dimer] 14 0.0060329 -0.0013 2.92472e-02 -17.7119 2.663 4.275 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6968 5.334 ------ ----
[Dimer] 15 0.0096878 -0.0013 6.09647e-02 -17.6968 2.667 4.285 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6878 5.311 ------ ----
[Dimer] 16 0.0082394 -0.0014 2.98130e-02 -17.6878 2.656 4.269 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6096 5.477 ------ ----
[Dimer] 17 0.0033329 -0.0014 1.08183e-02 -17.6096 2.739 4.420 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6385 5.435 ------ ----
[Dimer] 18 0.0006404 -0.0014 9.79839e-03 -17.6385 2.718 4.380 0
49 8.5049e-02 1.4403e-02 8 0.7353
Saddle point search started from reactant with energy -53.29783248901367 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -16.1554 8.328 17.75 1
[IDimerRot] ----- --------- ---------- ------------------ -17.3889 11.443 10.06 2
[IDimerRot] ----- --------- ---------- ------------------ -17.6775 5.387 3.54 3
[Dimer] 1 0.0002814 0.0000 1.29894e-02 -17.6775 5.387 3.538 3
[IDimerRot] ----- --------- ---------- ------------------ -17.6892 4.706 ------ ----
[Dimer] 2 0.0001902 0.0000 1.23134e-02 -17.6892 2.353 3.262 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6929 4.721 ------ ----
[Dimer] 3 0.0017193 -0.0000 2.46431e-02 -17.6929 2.361 3.800 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7094 4.844 ------ ----
[Dimer] 4 0.0021163 -0.0000 2.17872e-02 -17.7094 2.422 3.894 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7076 4.966 ------ ----
[Dimer] 5 0.0047722 -0.0001 2.90908e-02 -17.7076 2.483 3.991 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7221 5.045 ------ ----
[Dimer] 6 0.0194302 -0.0001 9.76366e-02 -17.7221 2.522 4.050 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6095 5.130 ------ ----
[Dimer] 7 0.0034557 -0.0002 2.93643e-02 -17.6095 2.565 4.144 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6664 5.007 ------ ----
[Dimer] 8 0.0045577 -0.0002 2.72607e-02 -17.6664 2.503 4.033 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6579 4.814 ------ ----
[Dimer] 9 0.0096527 -0.0002 3.07033e-02 -17.6579 2.407 3.881 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6155 4.657 ------ ----
[Dimer] 10 0.0121764 -0.0003 4.44298e-02 -17.6155 2.329 3.765 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5790 4.573 ------ ----
[Dimer] 11 0.0323106 -0.0003 4.43615e-02 -17.5790 2.287 3.706 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5127 4.577 ------ ----
[Dimer] 12 0.0171764 -0.0004 8.14639e-02 -17.5127 2.288 3.722 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5649 4.662 ------ ----
[Dimer] 13 0.0039748 -0.0004 4.24255e-02 -17.5649 2.331 3.780 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6247 4.847 ------ ----
[Dimer] 14 0.0066524 -0.0005 2.79617e-02 -17.6247 2.423 3.914 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7022 5.007 ------ ----
[Dimer] 15 0.0024747 -0.0005 3.17884e-02 -17.7022 2.503 4.025 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7241 5.003 ------ ----
[Dimer] 16 0.0252535 -0.0006 6.13975e-02 -17.7241 2.501 4.016 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7738 4.963 ------ ----
[Dimer] 17 0.0427684 -0.0007 8.79025e-02 -17.7738 2.481 3.974 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8035 4.907 ------ ----
[Dimer] 18 0.0263639 -0.0008 6.85145e-02 -17.8035 2.454 3.923 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8044 4.934 ------ ----
[Dimer] 19 0.1223202 -0.0011 9.38192e-02 -17.8044 2.467 3.945 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8388 5.339 ------ ----
[Dimer] 20 0.0986618 -0.0015 1.82899e-01 -17.8388 2.670 4.256 0
50 1.9483e-01 3.0641e+01 7 1.536
51 5.4677e-02 4.3515e+00 8 0.7344
52 3.6745e-02 2.3970e+00 8 0.7339
53 2.9026e-02 1.9361e+00 8 0.7309
54 6.5330e-02 2.4608e+00 8 0.7159
55 2.0116e-02 2.0944e+00 8 0.6984
56 6.4388e-02 4.9255e+00 8 0.6134
57 6.7568e-02 1.2881e+01 6 0.6063
58 4.4988e-02 2.2388e+00 8 0.6044
59 8.8201e-03 1.0773e+00 8 0.6013
60 2.0476e-02 1.6285e+00 8 0.5888
61 7.9713e-02 1.3163e+00 8 0.5709
62 5.0149e-02 1.8050e+00 8 0.5275
63 6.4346e-02 1.6669e+00 9 0.5479
64 5.5440e-02 1.4293e+00 9 0.6265
65 9.5589e-03 1.9552e+00 9 0.6199
66 8.4602e-02 1.6324e+00 9 0.7103
67 5.7304e-02 8.3646e-01 9 0.6438
68 1.3138e-02 1.8234e+00 9 0.6259
69 6.3265e-03 5.5313e-01 9 0.6339
70 1.2540e-02 1.7881e+00 9 0.6354
71 7.2947e-02 1.9618e+00 9 0.6807
72 1.3641e-02 1.2585e+00 9 0.6614
73 1.7967e-02 1.6433e+00 9 0.6485
74 2.6259e-02 7.0828e-01 9 0.6558
75 1.9809e-02 1.5554e+00 9 0.6575
76 2.6490e-02 9.0793e-01 9 0.6675
77 6.4910e-03 1.5080e+00 9 0.6656
78 2.7223e-02 1.0678e+00 9 0.6709
79 1.8411e-02 1.4784e+00 9 0.6709
80 3.6177e-02 1.4638e+00 9 0.6764
81 2.9488e-02 1.5086e+00 9 0.6796
82 4.1936e-02 9.2091e-01 9 0.684
83 2.0034e-02 1.3410e+00 9 0.6847
84 4.2011e-02 1.2793e+00 9 0.6913
85 2.1815e-02 1.2376e+00 9 0.6934
86 7.2402e-03 1.0445e+00 9 0.7087
87 1.1378e-02 7.6578e-01 9 0.728
88 2.9580e-02 1.1498e+00 9 0.7393
89 8.1096e-03 5.0772e-01 9 0.7351
90 6.1899e-03 2.7515e-01 9 0.7375
Saddle point search started from reactant with energy -53.29566192626953 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -18.4078 6.064 5.74 1
[IDimerRot] ----- --------- ---------- ------------------ -18.5945 3.947 3.28 2
[Dimer] 1 0.0067456 -0.0003 3.88448e-01 -18.5945 3.947 3.282 2
[IDimerRot] ----- --------- ---------- ------------------ -18.8184 5.006 ------ ----
[Dimer] 2 0.0055525 -0.0010 2.52693e-01 -18.8184 2.503 2.950 0
[IDimerRot] ----- --------- ---------- ------------------ -18.8600 4.932 ------ ----
[Dimer] 3 0.0125527 -0.0024 1.84014e-01 -18.8600 2.466 3.724 0
[IDimerRot] ----- --------- ---------- ------------------ -18.7337 5.426 ------ ----
[Dimer] 4 0.0092984 -0.0035 1.88039e-01 -18.7337 2.713 4.120 0
[IDimerRot] ----- --------- ---------- ------------------ -18.6534 5.951 ------ ----
[Dimer] 5 0.0312919 -0.0051 2.58389e-01 -18.6534 2.976 4.532 0
[IDimerRot] ----- --------- ---------- ------------------ -18.0694 3.816 1.70 1
[Dimer] 6 0.0079605 -0.0053 1.51212e-01 -18.0694 3.816 1.703 1
[IDimerRot] ----- --------- ---------- ------------------ -18.0253 4.295 ------ ----
[Dimer] 7 0.0020084 -0.0054 7.27090e-02 -18.0253 2.147 2.824 0
[IDimerRot] ----- --------- ---------- ------------------ -17.9528 4.157 ------ ----
[Dimer] 8 0.0033482 -0.0055 6.34534e-02 -17.9528 2.078 3.302 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7957 4.025 ------ ----
[Dimer] 9 0.0033912 -0.0056 6.86270e-02 -17.7957 2.012 3.226 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7005 3.975 ------ ----
[Dimer] 10 0.0066711 -0.0057 5.30532e-02 -17.7005 1.987 3.203 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6700 3.909 ------ ----
[Dimer] 11 0.0054548 -0.0055 2.27131e-01 -17.6700 1.955 3.156 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6305 4.272 ------ ----
[Dimer] 12 0.0032319 -0.0057 1.52195e-02 -17.6305 2.136 3.454 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7167 4.023 ------ ----
[Dimer] 13 0.0004327 -0.0057 1.01858e-02 -17.7167 2.012 3.239 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7351 4.039 ------ ----
[Dimer] 14 0.0007016 -0.0057 1.00529e-02 -17.7351 2.020 3.248 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7468 4.101 ------ ----
[Dimer] 15 0.0008676 -0.0057 1.14342e-02 -17.7468 2.050 3.295 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7510 4.102 ------ ----
[Dimer] 16 0.0021336 -0.0058 2.17523e-02 -17.7510 2.051 3.296 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7502 4.094 ------ ----
[Dimer] 17 0.0020158 -0.0058 1.40137e-02 -17.7502 2.047 3.289 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7535 4.054 ------ ----
[Dimer] 18 0.0004753 -0.0058 9.24782e-03 -17.7535 2.027 3.257 0
91 6.7926e-02 1.6959e-01 9 0.7322
Saddle point search started from reactant with energy -53.30095672607422 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -16.6053 7.030 6.78 1
[IDimerRot] ----- --------- ---------- ------------------ -17.2659 7.596 7.27 2
[IDimerRot] ----- --------- ---------- ------------------ -17.4071 3.505 2.64 3
[Dimer] 1 0.0010673 -0.0001 6.27775e-02 -17.4071 3.505 2.644 3
[IDimerRot] ----- --------- ---------- ------------------ -17.4177 3.392 ------ ----
[Dimer] 2 0.0003432 -0.0001 5.72816e-02 -17.4177 1.696 2.439 0
[IDimerRot] ----- --------- ---------- ------------------ -17.4111 3.388 ------ ----
[Dimer] 3 0.0069198 -0.0003 6.14517e-02 -17.4111 1.694 2.778 0
[IDimerRot] ----- --------- ---------- ------------------ -17.4194 3.466 ------ ----
[Dimer] 4 0.0044598 -0.0003 1.25898e-01 -17.4194 1.733 2.841 0
[IDimerRot] ----- --------- ---------- ------------------ -17.4174 3.446 ------ ----
[Dimer] 5 0.0009221 -0.0004 5.53167e-02 -17.4174 1.723 2.825 0
[IDimerRot] ----- --------- ---------- ------------------ -17.5114 3.543 ------ ----
[Dimer] 6 0.0023111 -0.0005 2.24385e-02 -17.5114 1.772 2.888 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6444 3.648 ------ ----
[Dimer] 7 0.0017957 -0.0005 2.19531e-02 -17.6444 1.824 2.951 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7161 3.718 ------ ----
[Dimer] 8 0.0023844 -0.0005 1.88001e-02 -17.7161 1.859 2.995 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7850 3.804 ------ ----
[Dimer] 9 0.0038464 -0.0005 6.80745e-02 -17.7850 1.902 3.052 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7428 3.860 ------ ----
[Dimer] 10 0.0011593 -0.0005 7.17333e-03 -17.7428 1.930 3.104 0
92 1.9614e-02 3.8861e-01 9 0.7315
Saddle point search started from reactant with energy -53.30168914794922 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.0598 6.287 7.32 1
[IDimerRot] ----- --------- ---------- ------------------ -17.3476 5.196 4.10 2
[Dimer] 1 0.0023570 -0.0002 1.40543e-01 -17.3476 5.196 4.103 2
[IDimerRot] ----- --------- ---------- ------------------ -17.3764 4.957 ------ ----
[Dimer] 2 0.0008629 -0.0002 1.26085e-01 -17.3764 2.479 3.814 0
[IDimerRot] ----- --------- ---------- ------------------ -17.4270 5.068 ------ ----
[Dimer] 3 0.0083748 0.0002 1.30180e-01 -17.4270 2.534 4.137 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6455 3.248 2.21 1
[Dimer] 4 0.0024208 0.0003 1.35435e-01 -17.6455 3.248 2.213 1
[IDimerRot] ----- --------- ---------- ------------------ -17.8527 4.989 ------ ----
[Dimer] 5 0.0013833 0.0002 1.27117e-02 -17.8527 2.494 2.919 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7320 5.114 ------ ----
[Dimer] 6 0.0001255 0.0002 1.09239e-02 -17.7320 2.557 4.103 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7360 5.128 ------ ----
[Dimer] 7 0.0009371 0.0002 1.71921e-02 -17.7360 2.564 4.113 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7500 5.085 ------ ----
[Dimer] 8 0.0012917 0.0002 1.82971e-02 -17.7500 2.542 4.076 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7583 5.069 ------ ----
[Dimer] 9 0.0019536 0.0002 2.80647e-02 -17.7583 2.534 4.061 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7426 5.108 ------ ----
[Dimer] 10 0.0008442 0.0002 9.21547e-03 -17.7426 2.554 4.095 0
93 1.5317e-02 3.4760e-01 9 0.7316
Saddle point search started from reactant with energy -53.30158615112305 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.3730 6.192 6.29 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6548 4.425 4.48 2
[Dimer] 1 0.0020058 -0.0002 9.74716e-02 -17.6548 4.425 4.478 2
[IDimerRot] ----- --------- ---------- ------------------ -17.6130 3.848 ------ ----
[Dimer] 2 0.0005768 -0.0002 8.85026e-02 -17.6130 1.924 2.912 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6225 3.946 ------ ----
[Dimer] 3 0.0054785 0.0000 8.81213e-02 -17.6225 1.973 3.194 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7525 4.895 ------ ----
[Dimer] 4 0.0022256 0.0000 9.61552e-02 -17.7525 2.448 3.925 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6771 5.210 ------ ----
[Dimer] 5 0.0014405 0.0000 6.59491e-03 -17.6771 2.605 4.192 0
94 8.3856e-03 6.2215e-02 9 0.7315
Saddle point search started from reactant with energy -53.30162048339844 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.4824 6.078 5.46 1
[IDimerRot] ----- --------- ---------- ------------------ -17.7486 4.314 4.33 2
[Dimer] 1 0.0004226 0.0000 2.99555e-02 -17.7486 4.314 4.329 2
[IDimerRot] ----- --------- ---------- ------------------ -17.7399 4.797 ------ ----
[Dimer] 2 0.0001944 0.0000 2.66313e-02 -17.7399 2.398 3.384 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7434 4.805 ------ ----
[Dimer] 3 0.0018740 0.0000 3.60346e-02 -17.7434 2.402 3.855 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7610 4.885 ------ ----
[Dimer] 4 0.0006778 0.0000 2.88690e-02 -17.7610 2.443 3.915 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7902 4.838 ------ ----
[Dimer] 5 0.0003530 0.0000 8.00247e-03 -17.7902 2.419 3.872 0
95 9.7215e-03 3.9727e-02 9 0.7315
Saddle point search started from reactant with energy -53.301666259765625 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.6301 5.992 4.68 1
[Dimer] 1 0.0008308 0.0000 4.31939e-02 -17.6301 5.992 4.683 1
[IDimerRot] ----- --------- ---------- ------------------ -17.7949 4.222 3.94 1
[Dimer] 2 0.0006057 0.0000 2.56175e-02 -17.7949 4.222 3.940 1
[IDimerRot] ----- --------- ---------- ------------------ -17.7753 5.421 ------ ----
[Dimer] 3 0.0007222 0.0000 1.04089e-02 -17.7753 2.710 3.657 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7751 5.432 ------ ----
[Dimer] 4 0.0002352 0.0000 9.67876e-03 -17.7751 2.716 4.344 0
96 1.4725e-02 1.3620e-01 9 0.7314
Saddle point search started from reactant with energy -53.301727294921875 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.6869 5.632 4.13 1
[Dimer] 1 0.0008875 -0.0000 5.56824e-02 -17.6869 5.632 4.131 1
[IDimerRot] ----- --------- ---------- ------------------ -17.8646 3.984 3.37 1
[Dimer] 2 0.0002931 -0.0000 4.68744e-02 -17.8646 3.984 3.372 1
[IDimerRot] ----- --------- ---------- ------------------ -17.8458 5.347 ------ ----
[Dimer] 3 0.0026252 0.0000 1.68900e-02 -17.8458 2.674 3.515 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7748 5.574 ------ ----
[Dimer] 4 0.0004773 0.0000 3.04919e-02 -17.7748 2.787 4.456 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7509 5.626 ------ ----
[Dimer] 5 0.0002662 0.0000 6.84089e-03 -17.7509 2.813 4.502 0
97 1.5321e-02 1.9075e-01 9 0.7315
Saddle point search started from reactant with energy -53.30166244506836 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.6996 5.464 3.64 1
[Dimer] 1 0.0010913 -0.0001 3.68170e-02 -17.6996 5.464 3.638 1
[IDimerRot] ----- --------- ---------- ------------------ -17.7153 7.561 ------ ----
[Dimer] 2 0.0001918 -0.0001 3.21716e-02 -17.7153 3.781 4.915 0
[IDimerRot] ----- --------- ---------- ------------------ -17.8275 3.769 1.76 1
[Dimer] 3 0.0019253 -0.0001 2.34463e-02 -17.8275 3.769 1.761 1
[IDimerRot] ----- --------- ---------- ------------------ -17.7860 5.043 ------ ----
[Dimer] 4 0.0006522 -0.0001 2.65289e-02 -17.7860 2.521 3.462 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7663 5.051 ------ ----
[Dimer] 5 0.0001817 -0.0001 1.31186e-02 -17.7663 2.526 4.045 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7735 5.006 ------ ----
[Dimer] 6 0.0003298 -0.0001 1.05159e-02 -17.7735 2.503 4.008 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7722 4.999 ------ ----
[Dimer] 7 0.0009807 -0.0001 1.58063e-02 -17.7722 2.500 4.003 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7552 4.966 ------ ----
[Dimer] 8 0.0012014 -0.0001 1.43548e-02 -17.7552 2.483 3.981 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7357 4.953 ------ ----
[Dimer] 9 0.0018016 -0.0001 1.41579e-02 -17.7357 2.477 3.975 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7279 4.929 ------ ----
[Dimer] 10 0.0030930 -0.0001 2.49768e-02 -17.7279 2.465 3.957 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7040 4.845 ------ ----
[Dimer] 11 0.0005617 -0.0001 9.34264e-03 -17.7040 2.423 3.896 0
98 3.5307e-02 5.4306e-01 9 0.7322
99 9.0356e-03 1.1793e-01 9 0.7314
Saddle point search started from reactant with energy -53.30179214477539 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.4503 5.120 4.17 1
[Dimer] 1 0.0007901 -0.0000 4.86714e-02 -17.4503 5.120 4.165 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6160 4.069 3.14 1
[Dimer] 2 0.0002834 -0.0000 4.05183e-02 -17.6160 4.069 3.139 1
[IDimerRot] ----- --------- ---------- ------------------ -17.6043 4.642 ------ ----
[Dimer] 3 0.0023568 -0.0001 2.70710e-02 -17.6043 2.321 3.206 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6054 4.717 ------ ----
[Dimer] 4 0.0014455 -0.0001 3.32475e-02 -17.6054 2.359 3.815 0
[IDimerRot] ----- --------- ---------- ------------------ -17.6293 4.736 ------ ----
[Dimer] 5 0.0026278 -0.0001 2.30855e-02 -17.6293 2.368 3.825 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7184 4.656 ------ ----
[Dimer] 6 0.0029157 -0.0001 2.65869e-02 -17.7184 2.328 3.742 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7489 4.645 ------ ----
[Dimer] 7 0.0017354 -0.0001 3.02649e-02 -17.7489 2.322 3.727 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7806 4.607 ------ ----
[Dimer] 8 0.0006492 -0.0001 9.49622e-03 -17.7806 2.304 3.691 0
100 8.4827e-03 4.8794e-02 9 0.7312
Saddle point search started from reactant with energy -53.30195617675781 eV.
[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots
[IDimerRot] ----- --------- ---------- ------------------ -17.4715 5.177 4.19 1
[Dimer] 1 0.0013349 0.0000 5.08536e-02 -17.4715 5.177 4.192 1
[IDimerRot] ----- --------- ---------- ------------------ -17.7672 3.859 3.10 1
[Dimer] 2 0.0009469 0.0000 2.10970e-02 -17.7672 3.859 3.103 1
[IDimerRot] ----- --------- ---------- ------------------ -17.7763 4.593 ------ ----
[Dimer] 3 0.0005979 0.0000 1.59646e-02 -17.7763 2.297 3.040 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7702 4.552 ------ ----
[Dimer] 4 0.0002152 0.0000 1.52322e-02 -17.7702 2.276 3.650 0
[IDimerRot] ----- --------- ---------- ------------------ -17.7617 4.644 ------ ----
[Dimer] 5 0.0000849 0.0000 8.17779e-03 -17.7617 2.322 3.724 0
Final state:
Nudged elastic band, too many iterations.
timing information:
real 28.322 seconds
user 54.499 seconds
sys 0.865 seconds
CompletedProcess(args=['eonclient'], returncode=0, stdout='EON Client\nVERSION: f01406f\nBUILD DATE: Tue Dec 02 01:53:51 AM GMT 2025\n\nOS: linux\nArch: x86_64\nHostname: runnervmoqczp\nPID: 2507\nDIR: /home/runner/work/atomistic-cookbook/atomistic-cookbook/examples/eon-pet-neb\nLoading parameter file config.ini\n* [Main] job: nudged_elastic_band\n* [Main] random_seed: 706253457\n* [Potential] potential: Metatomic\n* [Debug] write_movies: true\n* [Optimizer] opt_method: lbfgs\n* [Optimizer] converged_force: 0.01\n* [Optimizer] max_iterations: 100\n* [Optimizer] max_move: 0.1\n* [Metatomic] model_path: /home/runner/work/atomistic-cookbook/atomistic-cookbook/examples/eon-pet-neb/pet-mad-v1.0.2.pt\n* [Nudged Elastic Band] images: 10\n* [Nudged Elastic Band] climbing_image_method: true\n* [Nudged Elastic Band] climbing_image_converged_only: true\n* [Nudged Elastic Band] energy_weighted: true\n* [Nudged Elastic Band] ew_ksp_min: 0.972\n* [Nudged Elastic Band] ew_ksp_max: 9.72\n* [Nudged Elastic Band] initial_path_in: idppPath.dat\n* [Nudged Elastic Band] minimize_endpoints: false\n* [Nudged Elastic Band] ci_after: 1.5\n* [Nudged Elastic Band] ci_mmf: true\n* [Nudged Elastic Band] ci_mmf_after: 0.5\n* [Nudged Elastic Band] ci_mmf_nsteps: 20\nnebMinimEP == false: not minimizing endpoints.\nNudged elastic band calculation started.\n iteration step size ||Force|| max image max energy\n---------------------------------------------------------------\n\n 1 0.0000e+00 6.4530e+00 7 1.476\n 2 2.7135e-02 3.4314e+00 7 1.313\n 3 2.3174e-02 2.9734e+00 7 1.246\n 4 7.9659e-02 4.1153e+00 7 1.144\n 5 4.2119e-02 3.4505e+00 7 1.055\n 6 1.3126e-02 2.5218e+00 8 1.027\n 7 4.3104e-02 2.1791e+00 8 1.007\n 8 6.5064e-02 2.2961e+00 8 0.9603\n 9 7.1163e-02 2.1858e+00 8 0.9309\n 10 1.3802e-02 9.2779e-01 8 0.9\n 11 1.0319e-02 9.8000e-01 8 0.8936\n 12 1.8940e-02 1.3644e+00 8 0.8813\n 13 5.1184e-02 1.7403e+00 8 0.8499\n 14 3.3102e-02 1.0731e+00 8 0.8264\n 15 5.3943e-03 8.5662e-01 8 0.8209\n 16 1.4975e-02 9.7261e-01 8 0.8166\n 17 1.7579e-02 1.0942e+00 8 0.8065\n 18 4.1575e-02 1.0257e+00 8 0.7855\n 19 6.1904e-02 8.6704e-01 8 0.7693\n 20 2.4118e-02 6.7097e-01 8 0.7669\n 21 1.3001e-02 6.1148e-01 8 0.7647\n 22 3.0034e-02 5.7133e-01 8 0.7647\n 23 4.7631e-02 8.4101e-01 8 0.7684\n 24 4.1914e-02 1.7999e+00 8 0.7753\n 25 2.7705e-02 7.0061e-01 8 0.7605\n 26 8.5950e-03 4.6628e-01 8 0.7593\nSaddle point search started from reactant with energy -53.27388000488281 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7511 3.792 2.39 1\n[Dimer] 1 0.0186833 -0.0015 8.88839e-01 -17.7511 3.792 2.393 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8759 6.108 ------ ----\n[Dimer] 2 0.0165101 -0.0059 4.08662e-01 -17.8759 3.054 3.674 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7504 3.630 2.10 1\n[Dimer] 3 0.0139743 -0.0087 2.10910e-01 -17.7504 3.630 2.101 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6027 4.685 ------ ----\n[Dimer] 4 0.0088965 -0.0097 2.04202e-01 -17.6027 2.342 3.472 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6127 4.699 ------ ----\n[Dimer] 5 0.0186002 -0.0098 3.14657e-01 -17.6127 2.350 3.799 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -18.1543 4.708 ------ ----\n[Dimer] 6 0.0065046 -0.0100 1.71341e-01 -18.1543 2.354 3.694 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8074 4.594 ------ ----\n[Dimer] 7 0.0040005 -0.0102 4.36882e-02 -17.8074 2.297 3.675 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.9121 4.537 ------ ----\n[Dimer] 8 0.0007495 -0.0102 3.45538e-02 -17.9121 2.269 3.609 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.9171 4.530 ------ ----\n[Dimer] 9 0.0025785 -0.0103 4.86666e-02 -17.9171 2.265 3.603 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8873 4.550 ------ ----\n[Dimer] 10 0.0033890 -0.0104 5.63030e-02 -17.8873 2.275 3.624 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8446 4.515 ------ ----\n[Dimer] 11 0.0074655 -0.0104 7.51035e-02 -17.8446 2.257 3.605 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7884 4.327 ------ ----\n[Dimer] 12 0.0045799 -0.0104 9.79015e-02 -17.7884 2.163 3.467 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7174 4.251 ------ ----\n[Dimer] 13 0.0026098 -0.0105 1.99273e-02 -17.7174 2.125 3.420 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7402 4.292 ------ ----\n[Dimer] 14 0.0003734 -0.0105 1.47679e-02 -17.7402 2.146 3.449 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7279 4.278 ------ ----\n[Dimer] 15 0.0015581 -0.0105 2.75431e-02 -17.7279 2.139 3.440 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6926 4.253 ------ ----\n[Dimer] 16 0.0023492 -0.0105 4.26412e-02 -17.6926 2.127 3.427 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6535 4.207 ------ ----\n[Dimer] 17 0.0066182 -0.0105 6.18012e-02 -17.6535 2.104 3.398 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5841 4.114 ------ ----\n[Dimer] 18 0.0106054 -0.0106 6.60917e-02 -17.5841 2.057 3.336 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5231 4.004 ------ ----\n[Dimer] 19 0.0285095 -0.0107 1.78217e-01 -17.5231 2.002 3.259 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.3765 3.918 ------ ----\n[Dimer] 20 0.0082968 -0.0108 7.05261e-02 -17.3765 1.959 3.216 0\n\n 27 5.0021e-02 6.6456e-02 8 0.7485\nSaddle point search started from reactant with energy -53.28470230102539 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4712 3.807 2.30 1\n[Dimer] 1 0.0005027 -0.0000 3.91019e-02 -17.4712 3.807 2.300 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4542 3.656 ------ ----\n[Dimer] 2 0.0002679 -0.0000 3.66216e-02 -17.4542 1.828 2.688 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4563 3.655 ------ ----\n[Dimer] 3 0.0049187 -0.0001 4.55772e-02 -17.4563 1.827 2.988 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5819 3.524 ------ ----\n[Dimer] 4 0.0045667 -0.0002 4.99926e-02 -17.5819 1.762 2.862 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6530 3.657 ------ ----\n[Dimer] 5 0.0033963 -0.0002 4.24197e-02 -17.6530 1.828 2.956 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7338 3.784 ------ ----\n[Dimer] 6 0.0006650 -0.0002 2.52607e-02 -17.7338 1.892 3.045 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7051 3.848 ------ ----\n[Dimer] 7 0.0018353 -0.0002 2.15064e-02 -17.7051 1.924 3.101 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6821 4.038 ------ ----\n[Dimer] 8 0.0016934 -0.0003 2.36644e-02 -17.6821 2.019 3.257 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7007 4.146 ------ ----\n[Dimer] 9 0.0065684 -0.0003 5.61226e-02 -17.7007 2.073 3.340 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7781 4.266 ------ ----\n[Dimer] 10 0.0040765 -0.0003 2.82622e-02 -17.7781 2.133 3.421 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8319 4.239 ------ ----\n[Dimer] 11 0.0012096 -0.0003 2.50769e-02 -17.8319 2.120 3.389 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8399 4.108 ------ ----\n[Dimer] 12 0.0056081 -0.0004 3.33256e-02 -17.8399 2.054 3.284 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8419 3.860 ------ ----\n[Dimer] 13 0.0119066 -0.0005 5.63237e-02 -17.8419 1.930 3.087 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8616 3.739 ------ ----\n[Dimer] 14 0.0322391 -0.0007 9.39375e-02 -17.8616 1.869 2.988 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8062 3.640 ------ ----\n[Dimer] 15 0.0189463 -0.0008 5.57000e-02 -17.8062 1.820 2.918 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8234 3.894 ------ ----\n[Dimer] 16 0.0238557 -0.0009 4.64804e-02 -17.8234 1.947 3.117 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7842 4.410 ------ ----\n[Dimer] 17 0.0242690 -0.0010 5.36686e-02 -17.7842 2.205 3.534 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7350 4.851 ------ ----\n[Dimer] 18 0.0462534 -0.0012 7.98999e-02 -17.7350 2.426 3.894 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6298 5.583 ------ ----\n[Dimer] 19 0.0630176 -0.0014 1.14882e-01 -17.6298 2.791 4.499 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5742 3.186 1.73 1\n[Dimer] 20 0.0572431 -0.0014 1.87180e-01 -17.5742 3.186 1.728 1\n\n 28 1.5729e-01 6.0795e-01 8 0.7487\n 29 9.4161e-02 8.9647e-01 8 0.7518\n 30 3.4622e-02 4.2468e-01 8 0.7472\nSaddle point search started from reactant with energy -53.28594970703125 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.0447 4.583 5.30 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.3874 5.129 5.75 2\n[IDimerRot] ----- --------- ---------- ------------------ -17.3873 6.163 ------ ----\n[Dimer] 1 0.0030546 -0.0006 1.72727e-01 -17.3873 3.082 3.914 2\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5742 3.301 1.29 1\n[Dimer] 2 0.0011993 -0.0007 1.00714e-01 -17.5742 3.301 1.292 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6405 2.702 ------ ----\n[Dimer] 3 0.0020557 -0.0007 5.80249e-02 -17.6405 1.351 2.015 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6746 2.792 ------ ----\n[Dimer] 4 0.0015384 -0.0008 5.80295e-02 -17.6746 1.396 2.258 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6769 2.866 ------ ----\n[Dimer] 5 0.0052372 -0.0009 6.58701e-02 -17.6769 1.433 2.317 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6520 3.231 ------ ----\n[Dimer] 6 0.0086470 -0.0012 1.44669e-01 -17.6520 1.616 2.615 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6546 3.442 ------ ----\n[Dimer] 7 0.0039242 -0.0014 8.83451e-02 -17.6546 1.721 2.784 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6483 3.471 ------ ----\n[Dimer] 8 0.0139955 -0.0016 9.49049e-02 -17.6483 1.736 2.808 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6909 3.258 ------ ----\n[Dimer] 9 0.0153976 -0.0018 9.56474e-02 -17.6909 1.629 2.631 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7818 2.999 ------ ----\n[Dimer] 10 0.0191170 -0.0020 9.31785e-02 -17.7818 1.500 2.410 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.9805 3.383 ------ ----\n[Dimer] 11 0.0115168 -0.0020 9.20679e-02 -17.9805 1.691 2.687 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8820 3.168 ------ ----\n[Dimer] 12 0.0055241 -0.0021 3.88238e-02 -17.8820 1.584 2.531 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8995 3.175 ------ ----\n[Dimer] 13 0.0017167 -0.0021 3.32521e-02 -17.8995 1.587 2.534 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8821 3.183 ------ ----\n[Dimer] 14 0.0046314 -0.0022 4.48691e-02 -17.8821 1.592 2.543 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8063 3.192 ------ ----\n[Dimer] 15 0.0062909 -0.0022 5.80797e-02 -17.8063 1.596 2.561 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7693 3.206 ------ ----\n[Dimer] 16 0.0056639 -0.0022 2.46150e-02 -17.7693 1.603 2.577 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7177 3.252 ------ ----\n[Dimer] 17 0.0010379 -0.0023 1.72468e-02 -17.7177 1.626 2.622 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7322 3.253 ------ ----\n[Dimer] 18 0.0014562 -0.0023 2.34185e-02 -17.7322 1.626 2.620 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7488 3.267 ------ ----\n[Dimer] 19 0.0027184 -0.0023 3.51576e-02 -17.7488 1.633 2.629 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7450 3.273 ------ ----\n[Dimer] 20 0.0062004 -0.0023 4.40626e-02 -17.7450 1.637 2.635 0\n\n 31 1.3745e-01 1.0119e+00 8 0.7531\n 32 2.1750e-02 4.5098e-01 8 0.7482\nSaddle point search started from reactant with energy -53.28496170043945 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -16.7645 7.254 5.35 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.5518 7.131 8.78 2\n[IDimerRot] ----- --------- ---------- ------------------ -17.7578 4.072 3.54 3\n[Dimer] 1 0.0058540 -0.0013 2.25586e-01 -17.7578 4.072 3.541 3\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.9623 3.833 ------ ----\n[Dimer] 2 0.0026683 -0.0014 1.48721e-01 -17.9623 1.917 2.529 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8766 3.661 ------ ----\n[Dimer] 3 0.0007327 -0.0015 8.05034e-02 -17.8766 1.831 2.923 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8397 3.816 ------ ----\n[Dimer] 4 0.0015095 -0.0016 7.74910e-02 -17.8397 1.908 3.053 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7567 4.121 ------ ----\n[Dimer] 5 0.0026335 -0.0017 1.10506e-01 -17.7567 2.061 3.310 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6588 4.561 ------ ----\n[Dimer] 6 0.0061106 -0.0019 1.41959e-01 -17.6588 2.281 3.680 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5304 5.230 ------ ----\n[Dimer] 7 0.0113424 -0.0024 1.52377e-01 -17.5304 2.615 4.242 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5142 5.383 ------ ----\n[Dimer] 8 0.0269133 -0.0023 3.82310e-01 -17.5142 2.692 4.368 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.2801 4.656 ------ ----\n[Dimer] 9 0.0037183 -0.0029 1.01693e-01 -17.2801 2.328 3.836 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6056 4.172 ------ ----\n[Dimer] 10 0.0023577 -0.0030 6.89033e-02 -17.6056 2.086 3.378 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6890 3.576 ------ ----\n[Dimer] 11 0.0065721 -0.0032 9.22760e-02 -17.6890 1.788 2.886 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7212 2.918 ------ ----\n[Dimer] 12 0.0112717 -0.0035 1.06833e-01 -17.7212 1.459 2.353 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6557 2.726 ------ ----\n[Dimer] 13 0.0286049 -0.0039 1.31719e-01 -17.6557 1.363 2.207 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7469 3.194 ------ ----\n[Dimer] 14 0.0172980 -0.0038 1.68076e-01 -17.7469 1.597 2.571 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5554 4.121 ------ ----\n[Dimer] 15 0.0089408 -0.0040 3.02141e-02 -17.5554 2.061 3.348 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6044 4.080 ------ ----\n[Dimer] 16 0.0006806 -0.0040 1.96092e-02 -17.6044 2.040 3.305 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6147 4.240 ------ ----\n[Dimer] 17 0.0012700 -0.0040 2.32826e-02 -17.6147 2.120 3.432 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6142 4.424 ------ ----\n[Dimer] 18 0.0022913 -0.0041 3.38714e-02 -17.6142 2.212 3.579 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6325 4.479 ------ ----\n[Dimer] 19 0.0050936 -0.0041 3.19423e-02 -17.6325 2.240 3.620 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6296 4.494 ------ ----\n[Dimer] 20 0.0043822 -0.0041 4.93339e-02 -17.6296 2.247 3.632 0\n\n 33 6.1520e-02 2.2993e-01 8 0.7445\nSaddle point search started from reactant with energy -53.28866958618164 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.0445 6.393 5.36 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.6203 6.506 7.18 2\n[IDimerRot] ----- --------- ---------- ------------------ -17.6199 7.806 ------ ----\n[Dimer] 1 0.0022962 -0.0003 1.05137e-01 -17.6199 3.903 4.938 2\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6285 3.905 2.13 1\n[Dimer] 2 0.0009249 -0.0003 3.91022e-02 -17.6285 3.905 2.127 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5999 3.697 ------ ----\n[Dimer] 3 0.0003285 -0.0003 2.75477e-02 -17.5999 1.849 2.725 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6220 3.682 ------ ----\n[Dimer] 4 0.0009295 -0.0003 2.69741e-02 -17.6220 1.841 2.982 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6620 3.667 ------ ----\n[Dimer] 5 0.0013890 -0.0003 3.58414e-02 -17.6620 1.833 2.963 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6946 3.665 ------ ----\n[Dimer] 6 0.0036469 -0.0004 4.18327e-02 -17.6946 1.833 2.956 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7231 3.583 ------ ----\n[Dimer] 7 0.0071206 -0.0004 1.06688e-01 -17.7231 1.791 2.886 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7671 3.708 ------ ----\n[Dimer] 8 0.0013238 -0.0005 4.35987e-02 -17.7671 1.854 2.979 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7250 3.647 ------ ----\n[Dimer] 9 0.0022879 -0.0005 3.16250e-02 -17.7250 1.824 2.937 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6740 3.689 ------ ----\n[Dimer] 10 0.0033851 -0.0005 4.36517e-02 -17.6740 1.845 2.979 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6409 3.779 ------ ----\n[Dimer] 11 0.0065942 -0.0006 4.57707e-02 -17.6409 1.889 3.057 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6453 3.872 ------ ----\n[Dimer] 12 0.0192952 -0.0008 6.23592e-02 -17.6453 1.936 3.130 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5417 3.633 ------ ----\n[Dimer] 13 0.0445432 -0.0006 1.86832e-01 -17.5417 1.817 2.956 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5289 3.690 ------ ----\n[Dimer] 14 0.0209166 -0.0008 3.02960e-02 -17.5289 1.845 3.004 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6110 3.648 ------ ----\n[Dimer] 15 0.0029158 -0.0008 2.29262e-02 -17.6110 1.824 2.956 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6474 3.719 ------ ----\n[Dimer] 16 0.0048682 -0.0009 2.55737e-02 -17.6474 1.859 3.007 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7013 3.798 ------ ----\n[Dimer] 17 0.0040658 -0.0009 5.17058e-02 -17.7013 1.899 3.062 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7138 3.857 ------ ----\n[Dimer] 18 0.0047134 -0.0009 2.57320e-02 -17.7138 1.929 3.107 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6958 3.828 ------ ----\n[Dimer] 19 0.0182250 -0.0010 3.55781e-02 -17.6958 1.914 3.086 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6217 3.842 ------ ----\n[Dimer] 20 0.0096075 -0.0010 4.28253e-02 -17.6217 1.921 3.111 0\n\n 34 6.3226e-02 3.6384e-02 8 0.7435\nSaddle point search started from reactant with energy -53.2896728515625 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -16.5917 6.852 8.10 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.4231 8.292 8.89 2\n[IDimerRot] ----- --------- ---------- ------------------ -17.6258 4.408 3.23 3\n[Dimer] 1 0.0010842 -0.0000 4.23503e-02 -17.6258 4.408 3.227 3\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6381 3.583 ------ ----\n[Dimer] 2 0.0008253 -0.0000 3.89795e-02 -17.6381 1.791 2.610 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6489 3.637 ------ ----\n[Dimer] 3 0.0108179 -0.0001 9.55080e-02 -17.6489 1.819 2.942 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7245 4.014 ------ ----\n[Dimer] 4 0.0014020 -0.0002 8.80723e-02 -17.7245 2.007 3.230 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7277 3.978 ------ ----\n[Dimer] 5 0.0016016 -0.0002 2.77570e-02 -17.7277 1.989 3.200 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7158 3.889 ------ ----\n[Dimer] 6 0.0002632 -0.0002 2.10972e-02 -17.7158 1.944 3.132 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7143 3.836 ------ ----\n[Dimer] 7 0.0027060 -0.0002 3.54648e-02 -17.7143 1.918 3.090 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7035 3.568 ------ ----\n[Dimer] 8 0.0031177 -0.0003 4.59413e-02 -17.7035 1.784 2.877 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6939 3.422 ------ ----\n[Dimer] 9 0.0083902 -0.0003 5.89878e-02 -17.6939 1.711 2.762 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6905 3.296 ------ ----\n[Dimer] 10 0.0150948 -0.0004 9.83018e-02 -17.6905 1.648 2.661 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6460 3.198 ------ ----\n[Dimer] 11 0.0032537 -0.0004 4.16833e-02 -17.6460 1.599 2.589 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6882 3.427 ------ ----\n[Dimer] 12 0.0032497 -0.0005 4.00265e-02 -17.6882 1.713 2.767 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7132 3.809 ------ ----\n[Dimer] 13 0.0021409 -0.0005 3.81935e-02 -17.7132 1.904 3.068 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7169 3.886 ------ ----\n[Dimer] 14 0.0155190 -0.0006 1.11311e-01 -17.7169 1.943 3.130 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6627 3.759 ------ ----\n[Dimer] 15 0.0092842 -0.0006 2.96548e-02 -17.6627 1.880 3.037 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7083 3.778 ------ ----\n[Dimer] 16 0.0013749 -0.0006 1.84338e-02 -17.7083 1.889 3.044 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6829 3.680 ------ ----\n[Dimer] 17 0.0039840 -0.0006 2.03833e-02 -17.6829 1.840 2.970 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6650 3.573 ------ ----\n[Dimer] 18 0.0084938 -0.0007 4.01791e-02 -17.6650 1.787 2.887 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6044 3.540 ------ ----\n[Dimer] 19 0.0114645 -0.0007 4.83708e-02 -17.6044 1.770 2.871 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5852 3.534 ------ ----\n[Dimer] 20 0.0361547 -0.0009 8.19751e-02 -17.5852 1.767 2.869 0\n\n 35 6.5778e-02 7.3923e-02 8 0.7426\nSaddle point search started from reactant with energy -53.290592193603516 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -16.0971 7.563 10.07 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.2874 9.976 11.07 2\n[IDimerRot] ----- --------- ---------- ------------------ -17.5631 5.185 3.67 3\n[Dimer] 1 0.0022988 -0.0000 8.68632e-02 -17.5631 5.185 3.665 3\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7232 3.990 ------ ----\n[Dimer] 2 0.0017628 -0.0001 5.43971e-02 -17.7232 1.995 2.854 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7420 4.187 ------ ----\n[Dimer] 3 0.0029208 -0.0001 7.37342e-02 -17.7420 2.094 3.365 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7368 4.343 ------ ----\n[Dimer] 4 0.0008361 -0.0002 4.01829e-02 -17.7368 2.171 3.490 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7522 4.415 ------ ----\n[Dimer] 5 0.0004649 -0.0002 3.99705e-02 -17.7522 2.208 3.544 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7439 4.363 ------ ----\n[Dimer] 6 0.0035331 -0.0002 4.84407e-02 -17.7439 2.182 3.505 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6752 4.006 ------ ----\n[Dimer] 7 0.0039528 -0.0003 5.03209e-02 -17.6752 2.003 3.233 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6472 3.765 ------ ----\n[Dimer] 8 0.0074953 -0.0003 6.64698e-02 -17.6472 1.883 3.045 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5914 3.413 ------ ----\n[Dimer] 9 0.0019434 -0.0003 4.00183e-02 -17.5914 1.707 2.771 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6274 3.486 ------ ----\n[Dimer] 10 0.0043269 -0.0004 5.44886e-02 -17.6274 1.743 2.824 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6846 3.715 ------ ----\n[Dimer] 11 0.0064612 -0.0005 6.11097e-02 -17.6846 1.857 2.998 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7117 3.945 ------ ----\n[Dimer] 12 0.0160476 -0.0006 6.73859e-02 -17.7117 1.972 3.177 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7116 4.178 ------ ----\n[Dimer] 13 0.0216730 -0.0001 3.21701e-01 -17.7116 2.089 3.363 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5527 4.661 ------ ----\n[Dimer] 14 0.0113152 -0.0007 2.41866e-02 -17.5527 2.331 3.782 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6822 4.267 ------ ----\n[Dimer] 15 0.0018666 -0.0007 1.53560e-02 -17.6822 2.133 3.440 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6768 4.224 ------ ----\n[Dimer] 16 0.0026477 -0.0007 1.80750e-02 -17.6768 2.112 3.407 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6718 4.174 ------ ----\n[Dimer] 17 0.0031901 -0.0007 3.01389e-02 -17.6718 2.087 3.368 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6657 4.122 ------ ----\n[Dimer] 18 0.0075265 -0.0007 5.44092e-02 -17.6657 2.061 3.327 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6585 4.030 ------ ----\n[Dimer] 19 0.0128904 -0.0008 7.72842e-02 -17.6585 2.015 3.255 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6478 3.980 ------ ----\n[Dimer] 20 0.0221173 -0.0009 1.02437e-01 -17.6478 1.990 3.217 0\n\n 36 1.0000e-01 5.6394e-01 8 0.7429\n 37 6.1011e-02 4.0354e+00 8 0.7926\n 38 3.7416e-02 1.4381e+00 8 0.7476\n 39 1.9157e-02 6.8016e-01 8 0.7437\n 40 8.6507e-02 1.2682e+00 8 0.7467\n 41 5.8430e-02 1.0812e+00 8 0.744\n 42 6.5042e-02 9.6822e-01 8 0.7431\n 43 1.3735e-02 4.7441e-01 8 0.7412\nSaddle point search started from reactant with energy -53.291969299316406 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -16.7674 5.917 6.40 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.0437 4.890 4.22 2\n[Dimer] 1 0.0030486 -0.0004 1.92566e-01 -17.0437 4.890 4.221 2\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.1845 6.099 ------ ----\n[Dimer] 2 0.0011055 -0.0002 1.59725e-01 -17.1845 3.049 4.391 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.3272 3.103 1.78 1\n[Dimer] 3 0.0071266 0.0003 6.09842e-02 -17.3272 3.103 1.777 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5644 5.246 ------ ----\n[Dimer] 4 0.0026741 0.0003 1.21334e-01 -17.5644 2.623 3.538 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6184 5.749 ------ ----\n[Dimer] 5 0.0013762 0.0003 1.54177e-02 -17.6184 2.874 4.633 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5934 5.562 ------ ----\n[Dimer] 6 0.0001752 0.0003 1.46934e-02 -17.5934 2.781 4.491 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5901 5.545 ------ ----\n[Dimer] 7 0.0061044 0.0002 4.77200e-02 -17.5901 2.772 4.478 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5149 5.451 ------ ----\n[Dimer] 8 0.0085683 0.0001 7.06336e-02 -17.5149 2.725 4.422 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4694 5.420 ------ ----\n[Dimer] 9 0.0195888 0.0000 8.72174e-02 -17.4694 2.710 4.409 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5022 5.156 ------ ----\n[Dimer] 10 0.0235357 0.0001 1.30143e-01 -17.5022 2.578 4.190 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6635 5.302 ------ ----\n[Dimer] 11 0.0093507 -0.0001 3.44241e-02 -17.6635 2.651 4.268 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6298 5.232 ------ ----\n[Dimer] 12 0.0032069 -0.0001 2.09172e-02 -17.6298 2.616 4.220 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6419 5.295 ------ ----\n[Dimer] 13 0.0015523 -0.0001 2.52486e-02 -17.6419 2.648 4.267 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6490 5.386 ------ ----\n[Dimer] 14 0.0034716 -0.0001 3.55007e-02 -17.6490 2.693 4.337 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6491 5.477 ------ ----\n[Dimer] 15 0.0100648 -0.0002 5.79383e-02 -17.6491 2.738 4.410 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6139 5.542 ------ ----\n[Dimer] 16 0.0128877 -0.0002 3.31817e-02 -17.6139 2.771 4.470 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5705 5.619 ------ ----\n[Dimer] 17 0.0036306 -0.0002 2.22766e-02 -17.5705 2.809 4.542 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5702 5.587 ------ ----\n[Dimer] 18 0.0116286 -0.0003 4.75761e-02 -17.5702 2.793 4.517 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5400 5.483 ------ ----\n[Dimer] 19 0.0139421 -0.0003 7.92040e-02 -17.5400 2.741 4.442 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5147 5.380 ------ ----\n[Dimer] 20 0.0335641 -0.0004 1.34018e-01 -17.5147 2.690 4.366 0\n\n 44 7.1147e-02 1.5175e-01 8 0.7408\nSaddle point search started from reactant with energy -53.29235076904297 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.2395 5.613 4.88 1\n[Dimer] 1 0.0012462 -0.0001 7.60688e-02 -17.2395 5.613 4.883 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4879 4.621 4.05 1\n[Dimer] 2 0.0005151 -0.0001 4.42320e-02 -17.4879 4.621 4.046 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4786 6.413 ------ ----\n[Dimer] 3 0.0009196 -0.0001 4.30234e-02 -17.4786 3.206 4.054 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5805 3.216 1.46 1\n[Dimer] 4 0.0015382 -0.0002 5.70119e-02 -17.5805 3.216 1.458 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6308 3.557 ------ ----\n[Dimer] 5 0.0031501 -0.0002 5.83977e-02 -17.6308 1.778 2.504 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7143 3.594 ------ ----\n[Dimer] 6 0.0059079 -0.0003 6.68993e-02 -17.7143 1.797 2.897 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6899 3.564 ------ ----\n[Dimer] 7 0.0151078 -0.0002 2.16318e-01 -17.6899 1.782 2.876 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.9165 3.734 ------ ----\n[Dimer] 8 0.0047909 -0.0004 3.89769e-02 -17.9165 1.867 2.975 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6978 3.573 ------ ----\n[Dimer] 9 0.0009653 -0.0005 3.28469e-02 -17.6978 1.787 2.882 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6749 3.557 ------ ----\n[Dimer] 10 0.0059571 -0.0005 5.12552e-02 -17.6749 1.778 2.873 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6555 3.509 ------ ----\n[Dimer] 11 0.0067618 -0.0006 6.08301e-02 -17.6555 1.755 2.838 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6899 3.523 ------ ----\n[Dimer] 12 0.0171549 -0.0006 1.49710e-01 -17.6899 1.762 2.843 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6630 3.552 ------ ----\n[Dimer] 13 0.0023283 -0.0007 4.23069e-02 -17.6630 1.776 2.870 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7490 3.629 ------ ----\n[Dimer] 14 0.0023298 -0.0008 2.88841e-02 -17.7490 1.814 2.918 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7600 3.660 ------ ----\n[Dimer] 15 0.0020412 -0.0008 2.60318e-02 -17.7600 1.830 2.942 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7259 3.683 ------ ----\n[Dimer] 16 0.0055852 -0.0008 5.83544e-02 -17.7259 1.841 2.965 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6929 3.703 ------ ----\n[Dimer] 17 0.0049592 -0.0008 2.25393e-02 -17.6929 1.851 2.987 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6534 3.709 ------ ----\n[Dimer] 18 0.0038781 -0.0009 1.75767e-02 -17.6534 1.854 2.998 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6421 3.703 ------ ----\n[Dimer] 19 0.0043744 -0.0009 2.11259e-02 -17.6421 1.852 2.996 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6420 3.731 ------ ----\n[Dimer] 20 0.0041962 -0.0009 2.51169e-02 -17.6420 1.866 3.018 0\n\n 45 5.2000e-02 2.0492e-02 8 0.7399\nSaddle point search started from reactant with energy -53.293243408203125 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.3519 5.622 4.45 1\n[Dimer] 1 0.0006987 -0.0000 2.71440e-02 -17.3519 5.622 4.447 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6070 4.568 4.05 1\n[Dimer] 2 0.0005624 -0.0000 2.43417e-02 -17.6070 4.568 4.054 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6068 5.973 ------ ----\n[Dimer] 3 0.0069810 -0.0000 1.06123e-01 -17.6068 2.986 3.860 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6459 6.136 ------ ----\n[Dimer] 4 0.0007810 -0.0001 4.00513e-02 -17.6459 3.068 4.932 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6549 5.929 ------ ----\n[Dimer] 5 0.0005969 -0.0001 2.51947e-02 -17.6549 2.964 4.766 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6435 5.977 ------ ----\n[Dimer] 6 0.0007323 -0.0001 2.06435e-02 -17.6435 2.988 4.806 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6246 6.015 ------ ----\n[Dimer] 7 0.0034097 -0.0001 4.44849e-02 -17.6246 3.008 4.842 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5941 6.098 ------ ----\n[Dimer] 8 0.0064479 -0.0002 5.92689e-02 -17.5941 3.049 4.916 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5803 6.161 ------ ----\n[Dimer] 9 0.0112059 -0.0002 6.30787e-02 -17.5803 3.080 4.969 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5999 6.142 ------ ----\n[Dimer] 10 0.0290662 -0.0001 2.30890e-01 -17.5999 3.071 4.949 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6247 3.205 1.36 1\n[Dimer] 11 0.0056721 -0.0003 4.42215e-02 -17.6247 3.205 1.363 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7252 3.485 ------ ----\n[Dimer] 12 0.0027816 -0.0003 2.59118e-02 -17.7252 1.742 2.358 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7401 3.615 ------ ----\n[Dimer] 13 0.0030502 -0.0004 3.13948e-02 -17.7401 1.808 2.909 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7513 3.749 ------ ----\n[Dimer] 14 0.0087372 -0.0004 9.63984e-02 -17.7513 1.875 3.014 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6448 3.771 ------ ----\n[Dimer] 15 0.0033204 -0.0005 3.56212e-02 -17.6448 1.886 3.050 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6887 3.738 ------ ----\n[Dimer] 16 0.0071304 -0.0005 2.63685e-02 -17.6887 1.869 3.015 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6902 3.657 ------ ----\n[Dimer] 17 0.0133123 -0.0005 4.45056e-02 -17.6902 1.829 2.951 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6645 3.559 ------ ----\n[Dimer] 18 0.0134047 -0.0005 4.90727e-02 -17.6645 1.780 2.877 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6465 3.511 ------ ----\n[Dimer] 19 0.0211755 -0.0006 6.33420e-02 -17.6465 1.756 2.841 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6356 3.576 ------ ----\n[Dimer] 20 0.0358747 -0.0007 1.39034e-01 -17.6356 1.788 2.895 0\n\n 46 7.9178e-02 1.3608e-01 8 0.7392\nSaddle point search started from reactant with energy -53.29399490356445 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5113 5.640 4.14 1\n[Dimer] 1 0.0018230 -0.0001 5.23020e-02 -17.5113 5.640 4.136 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5830 4.182 3.45 1\n[Dimer] 2 0.0006216 -0.0001 4.55288e-02 -17.5830 4.182 3.447 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5930 5.292 ------ ----\n[Dimer] 3 0.0053472 -0.0003 5.63808e-02 -17.5930 2.646 3.492 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6254 5.330 ------ ----\n[Dimer] 4 0.0051930 -0.0002 2.41292e-01 -17.6254 2.665 4.299 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6462 5.064 ------ ----\n[Dimer] 5 0.0018575 -0.0003 4.12787e-02 -17.6462 2.532 4.083 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6200 5.314 ------ ----\n[Dimer] 6 0.0003325 -0.0003 3.20612e-02 -17.6200 2.657 4.288 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6198 5.339 ------ ----\n[Dimer] 7 0.0023856 -0.0004 5.30913e-02 -17.6198 2.669 4.308 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6292 5.401 ------ ----\n[Dimer] 8 0.0024710 -0.0004 6.02418e-02 -17.6292 2.701 4.355 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6584 5.435 ------ ----\n[Dimer] 9 0.0056706 -0.0005 5.61435e-02 -17.6584 2.717 4.374 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6678 5.461 ------ ----\n[Dimer] 10 0.0102061 -0.0005 7.82805e-02 -17.6678 2.730 4.392 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7111 5.494 ------ ----\n[Dimer] 11 0.0033309 -0.0006 3.97699e-02 -17.7111 2.747 4.408 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6964 5.488 ------ ----\n[Dimer] 12 0.0024755 -0.0006 3.87575e-02 -17.6964 2.744 4.407 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6713 5.505 ------ ----\n[Dimer] 13 0.0054800 -0.0007 4.10462e-02 -17.6713 2.752 4.427 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6640 5.537 ------ ----\n[Dimer] 14 0.0113737 -0.0007 8.96727e-02 -17.6640 2.769 4.454 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6078 5.631 ------ ----\n[Dimer] 15 0.0036000 -0.0008 3.65381e-02 -17.6078 2.816 4.543 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6405 5.617 ------ ----\n[Dimer] 16 0.0048106 -0.0008 2.26379e-02 -17.6405 2.809 4.523 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6580 5.629 ------ ----\n[Dimer] 17 0.0051004 -0.0008 2.21227e-02 -17.6580 2.815 4.528 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6436 5.676 ------ ----\n[Dimer] 18 0.0039644 -0.0008 6.64481e-02 -17.6436 2.838 4.569 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6163 5.680 ------ ----\n[Dimer] 19 0.0016364 -0.0009 1.85967e-02 -17.6163 2.840 4.579 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6188 5.666 ------ ----\n[Dimer] 20 0.0006485 -0.0009 1.36509e-02 -17.6188 2.833 4.567 0\n\n 47 4.3263e-02 1.3157e-02 8 0.7383\nSaddle point search started from reactant with energy -53.294864654541016 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4148 5.636 3.95 1\n[Dimer] 1 0.0019554 -0.0000 3.47009e-02 -17.4148 5.636 3.949 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6644 4.326 3.75 1\n[Dimer] 2 0.0017770 -0.0000 2.92930e-02 -17.6644 4.326 3.752 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6949 5.372 ------ ----\n[Dimer] 3 0.0026394 -0.0000 3.51614e-02 -17.6949 2.686 3.602 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7084 5.416 ------ ----\n[Dimer] 4 0.0026682 -0.0001 3.68958e-02 -17.7084 2.708 4.347 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7031 5.417 ------ ----\n[Dimer] 5 0.0086244 -0.0001 4.27433e-02 -17.7031 2.708 4.349 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6804 5.439 ------ ----\n[Dimer] 6 0.0110322 -0.0001 1.09448e-01 -17.6804 2.719 4.372 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5864 5.567 ------ ----\n[Dimer] 7 0.0008071 -0.0002 4.18014e-02 -17.5864 2.783 4.497 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6370 5.514 ------ ----\n[Dimer] 8 0.0020352 -0.0002 2.55514e-02 -17.6370 2.757 4.442 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6720 5.491 ------ ----\n[Dimer] 9 0.0053278 -0.0002 3.67761e-02 -17.6720 2.746 4.415 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7110 5.439 ------ ----\n[Dimer] 10 0.0097348 -0.0003 4.02431e-02 -17.7110 2.719 4.365 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7361 5.415 ------ ----\n[Dimer] 11 0.0258366 -0.0004 5.79182e-02 -17.7361 2.707 4.339 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6987 5.414 ------ ----\n[Dimer] 12 0.0353756 -0.0003 1.47229e-01 -17.6987 2.707 4.348 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6992 5.413 ------ ----\n[Dimer] 13 0.0134993 -0.0005 3.20373e-02 -17.6992 2.707 4.347 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6496 5.450 ------ ----\n[Dimer] 14 0.0014537 -0.0005 2.42081e-02 -17.6496 2.725 4.388 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6346 5.440 ------ ----\n[Dimer] 15 0.0108303 -0.0006 4.66644e-02 -17.6346 2.720 4.384 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5856 5.477 ------ ----\n[Dimer] 16 0.0115998 -0.0006 5.96404e-02 -17.5856 2.738 4.425 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5753 5.494 ------ ----\n[Dimer] 17 0.0170947 -0.0007 5.26898e-02 -17.5753 2.747 4.442 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5857 5.483 ------ ----\n[Dimer] 18 0.0517138 -0.0008 1.03790e-01 -17.5857 2.741 4.430 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5276 5.524 ------ ----\n[Dimer] 19 0.0496689 -0.0011 9.91694e-02 -17.5276 2.762 4.477 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5962 5.597 ------ ----\n[Dimer] 20 0.1148644 -0.0017 1.64317e-01 -17.5962 2.798 4.518 0\n\n 48 1.6933e-01 3.6113e-01 8 0.7367\nSaddle point search started from reactant with energy -53.29645538330078 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -15.7596 8.091 9.18 1\n[IDimerRot] ----- --------- ---------- ------------------ -16.9228 9.534 10.60 2\n[IDimerRot] ----- --------- ---------- ------------------ -17.3157 5.939 4.88 3\n[Dimer] 1 0.0033396 -0.0006 1.28283e-01 -17.3157 5.939 4.877 3\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5357 5.443 ------ ----\n[Dimer] 2 0.0011015 -0.0007 6.48642e-02 -17.5357 2.721 3.554 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5415 5.334 ------ ----\n[Dimer] 3 0.0005531 -0.0007 5.51835e-02 -17.5415 2.667 4.322 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5328 5.335 ------ ----\n[Dimer] 4 0.0033139 -0.0008 5.65086e-02 -17.5328 2.668 4.326 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5509 5.323 ------ ----\n[Dimer] 5 0.0037999 -0.0009 5.34274e-02 -17.5509 2.661 4.311 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5802 5.316 ------ ----\n[Dimer] 6 0.0059193 -0.0009 1.07179e-01 -17.5802 2.658 4.299 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6652 5.420 ------ ----\n[Dimer] 7 0.0011464 -0.0010 3.64614e-02 -17.6652 2.710 4.361 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6356 5.413 ------ ----\n[Dimer] 8 0.0010649 -0.0010 2.98751e-02 -17.6356 2.706 4.362 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6246 5.408 ------ ----\n[Dimer] 9 0.0035534 -0.0010 3.65925e-02 -17.6246 2.704 4.361 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6315 5.411 ------ ----\n[Dimer] 10 0.0044505 -0.0011 3.52658e-02 -17.6315 2.705 4.362 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6772 5.379 ------ ----\n[Dimer] 11 0.0106824 -0.0011 7.73657e-02 -17.6772 2.689 4.325 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7067 5.392 ------ ----\n[Dimer] 12 0.0044664 -0.0012 3.21949e-02 -17.7067 2.696 4.328 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7202 5.362 ------ ----\n[Dimer] 13 0.0021347 -0.0012 3.23380e-02 -17.7202 2.681 4.302 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7119 5.326 ------ ----\n[Dimer] 14 0.0060329 -0.0013 2.92472e-02 -17.7119 2.663 4.275 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6968 5.334 ------ ----\n[Dimer] 15 0.0096878 -0.0013 6.09647e-02 -17.6968 2.667 4.285 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6878 5.311 ------ ----\n[Dimer] 16 0.0082394 -0.0014 2.98130e-02 -17.6878 2.656 4.269 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6096 5.477 ------ ----\n[Dimer] 17 0.0033329 -0.0014 1.08183e-02 -17.6096 2.739 4.420 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6385 5.435 ------ ----\n[Dimer] 18 0.0006404 -0.0014 9.79839e-03 -17.6385 2.718 4.380 0\n\n 49 8.5049e-02 1.4403e-02 8 0.7353\nSaddle point search started from reactant with energy -53.29783248901367 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -16.1554 8.328 17.75 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.3889 11.443 10.06 2\n[IDimerRot] ----- --------- ---------- ------------------ -17.6775 5.387 3.54 3\n[Dimer] 1 0.0002814 0.0000 1.29894e-02 -17.6775 5.387 3.538 3\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6892 4.706 ------ ----\n[Dimer] 2 0.0001902 0.0000 1.23134e-02 -17.6892 2.353 3.262 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6929 4.721 ------ ----\n[Dimer] 3 0.0017193 -0.0000 2.46431e-02 -17.6929 2.361 3.800 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7094 4.844 ------ ----\n[Dimer] 4 0.0021163 -0.0000 2.17872e-02 -17.7094 2.422 3.894 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7076 4.966 ------ ----\n[Dimer] 5 0.0047722 -0.0001 2.90908e-02 -17.7076 2.483 3.991 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7221 5.045 ------ ----\n[Dimer] 6 0.0194302 -0.0001 9.76366e-02 -17.7221 2.522 4.050 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6095 5.130 ------ ----\n[Dimer] 7 0.0034557 -0.0002 2.93643e-02 -17.6095 2.565 4.144 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6664 5.007 ------ ----\n[Dimer] 8 0.0045577 -0.0002 2.72607e-02 -17.6664 2.503 4.033 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6579 4.814 ------ ----\n[Dimer] 9 0.0096527 -0.0002 3.07033e-02 -17.6579 2.407 3.881 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6155 4.657 ------ ----\n[Dimer] 10 0.0121764 -0.0003 4.44298e-02 -17.6155 2.329 3.765 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5790 4.573 ------ ----\n[Dimer] 11 0.0323106 -0.0003 4.43615e-02 -17.5790 2.287 3.706 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5127 4.577 ------ ----\n[Dimer] 12 0.0171764 -0.0004 8.14639e-02 -17.5127 2.288 3.722 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5649 4.662 ------ ----\n[Dimer] 13 0.0039748 -0.0004 4.24255e-02 -17.5649 2.331 3.780 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6247 4.847 ------ ----\n[Dimer] 14 0.0066524 -0.0005 2.79617e-02 -17.6247 2.423 3.914 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7022 5.007 ------ ----\n[Dimer] 15 0.0024747 -0.0005 3.17884e-02 -17.7022 2.503 4.025 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7241 5.003 ------ ----\n[Dimer] 16 0.0252535 -0.0006 6.13975e-02 -17.7241 2.501 4.016 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7738 4.963 ------ ----\n[Dimer] 17 0.0427684 -0.0007 8.79025e-02 -17.7738 2.481 3.974 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8035 4.907 ------ ----\n[Dimer] 18 0.0263639 -0.0008 6.85145e-02 -17.8035 2.454 3.923 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8044 4.934 ------ ----\n[Dimer] 19 0.1223202 -0.0011 9.38192e-02 -17.8044 2.467 3.945 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8388 5.339 ------ ----\n[Dimer] 20 0.0986618 -0.0015 1.82899e-01 -17.8388 2.670 4.256 0\n\n 50 1.9483e-01 3.0641e+01 7 1.536\n 51 5.4677e-02 4.3515e+00 8 0.7344\n 52 3.6745e-02 2.3970e+00 8 0.7339\n 53 2.9026e-02 1.9361e+00 8 0.7309\n 54 6.5330e-02 2.4608e+00 8 0.7159\n 55 2.0116e-02 2.0944e+00 8 0.6984\n 56 6.4388e-02 4.9255e+00 8 0.6134\n 57 6.7568e-02 1.2881e+01 6 0.6063\n 58 4.4988e-02 2.2388e+00 8 0.6044\n 59 8.8201e-03 1.0773e+00 8 0.6013\n 60 2.0476e-02 1.6285e+00 8 0.5888\n 61 7.9713e-02 1.3163e+00 8 0.5709\n 62 5.0149e-02 1.8050e+00 8 0.5275\n 63 6.4346e-02 1.6669e+00 9 0.5479\n 64 5.5440e-02 1.4293e+00 9 0.6265\n 65 9.5589e-03 1.9552e+00 9 0.6199\n 66 8.4602e-02 1.6324e+00 9 0.7103\n 67 5.7304e-02 8.3646e-01 9 0.6438\n 68 1.3138e-02 1.8234e+00 9 0.6259\n 69 6.3265e-03 5.5313e-01 9 0.6339\n 70 1.2540e-02 1.7881e+00 9 0.6354\n 71 7.2947e-02 1.9618e+00 9 0.6807\n 72 1.3641e-02 1.2585e+00 9 0.6614\n 73 1.7967e-02 1.6433e+00 9 0.6485\n 74 2.6259e-02 7.0828e-01 9 0.6558\n 75 1.9809e-02 1.5554e+00 9 0.6575\n 76 2.6490e-02 9.0793e-01 9 0.6675\n 77 6.4910e-03 1.5080e+00 9 0.6656\n 78 2.7223e-02 1.0678e+00 9 0.6709\n 79 1.8411e-02 1.4784e+00 9 0.6709\n 80 3.6177e-02 1.4638e+00 9 0.6764\n 81 2.9488e-02 1.5086e+00 9 0.6796\n 82 4.1936e-02 9.2091e-01 9 0.684\n 83 2.0034e-02 1.3410e+00 9 0.6847\n 84 4.2011e-02 1.2793e+00 9 0.6913\n 85 2.1815e-02 1.2376e+00 9 0.6934\n 86 7.2402e-03 1.0445e+00 9 0.7087\n 87 1.1378e-02 7.6578e-01 9 0.728\n 88 2.9580e-02 1.1498e+00 9 0.7393\n 89 8.1096e-03 5.0772e-01 9 0.7351\n 90 6.1899e-03 2.7515e-01 9 0.7375\nSaddle point search started from reactant with energy -53.29566192626953 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -18.4078 6.064 5.74 1\n[IDimerRot] ----- --------- ---------- ------------------ -18.5945 3.947 3.28 2\n[Dimer] 1 0.0067456 -0.0003 3.88448e-01 -18.5945 3.947 3.282 2\n\n[IDimerRot] ----- --------- ---------- ------------------ -18.8184 5.006 ------ ----\n[Dimer] 2 0.0055525 -0.0010 2.52693e-01 -18.8184 2.503 2.950 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -18.8600 4.932 ------ ----\n[Dimer] 3 0.0125527 -0.0024 1.84014e-01 -18.8600 2.466 3.724 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -18.7337 5.426 ------ ----\n[Dimer] 4 0.0092984 -0.0035 1.88039e-01 -18.7337 2.713 4.120 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -18.6534 5.951 ------ ----\n[Dimer] 5 0.0312919 -0.0051 2.58389e-01 -18.6534 2.976 4.532 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -18.0694 3.816 1.70 1\n[Dimer] 6 0.0079605 -0.0053 1.51212e-01 -18.0694 3.816 1.703 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -18.0253 4.295 ------ ----\n[Dimer] 7 0.0020084 -0.0054 7.27090e-02 -18.0253 2.147 2.824 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.9528 4.157 ------ ----\n[Dimer] 8 0.0033482 -0.0055 6.34534e-02 -17.9528 2.078 3.302 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7957 4.025 ------ ----\n[Dimer] 9 0.0033912 -0.0056 6.86270e-02 -17.7957 2.012 3.226 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7005 3.975 ------ ----\n[Dimer] 10 0.0066711 -0.0057 5.30532e-02 -17.7005 1.987 3.203 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6700 3.909 ------ ----\n[Dimer] 11 0.0054548 -0.0055 2.27131e-01 -17.6700 1.955 3.156 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6305 4.272 ------ ----\n[Dimer] 12 0.0032319 -0.0057 1.52195e-02 -17.6305 2.136 3.454 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7167 4.023 ------ ----\n[Dimer] 13 0.0004327 -0.0057 1.01858e-02 -17.7167 2.012 3.239 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7351 4.039 ------ ----\n[Dimer] 14 0.0007016 -0.0057 1.00529e-02 -17.7351 2.020 3.248 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7468 4.101 ------ ----\n[Dimer] 15 0.0008676 -0.0057 1.14342e-02 -17.7468 2.050 3.295 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7510 4.102 ------ ----\n[Dimer] 16 0.0021336 -0.0058 2.17523e-02 -17.7510 2.051 3.296 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7502 4.094 ------ ----\n[Dimer] 17 0.0020158 -0.0058 1.40137e-02 -17.7502 2.047 3.289 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7535 4.054 ------ ----\n[Dimer] 18 0.0004753 -0.0058 9.24782e-03 -17.7535 2.027 3.257 0\n\n 91 6.7926e-02 1.6959e-01 9 0.7322\nSaddle point search started from reactant with energy -53.30095672607422 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -16.6053 7.030 6.78 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.2659 7.596 7.27 2\n[IDimerRot] ----- --------- ---------- ------------------ -17.4071 3.505 2.64 3\n[Dimer] 1 0.0010673 -0.0001 6.27775e-02 -17.4071 3.505 2.644 3\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4177 3.392 ------ ----\n[Dimer] 2 0.0003432 -0.0001 5.72816e-02 -17.4177 1.696 2.439 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4111 3.388 ------ ----\n[Dimer] 3 0.0069198 -0.0003 6.14517e-02 -17.4111 1.694 2.778 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4194 3.466 ------ ----\n[Dimer] 4 0.0044598 -0.0003 1.25898e-01 -17.4194 1.733 2.841 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4174 3.446 ------ ----\n[Dimer] 5 0.0009221 -0.0004 5.53167e-02 -17.4174 1.723 2.825 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.5114 3.543 ------ ----\n[Dimer] 6 0.0023111 -0.0005 2.24385e-02 -17.5114 1.772 2.888 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6444 3.648 ------ ----\n[Dimer] 7 0.0017957 -0.0005 2.19531e-02 -17.6444 1.824 2.951 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7161 3.718 ------ ----\n[Dimer] 8 0.0023844 -0.0005 1.88001e-02 -17.7161 1.859 2.995 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7850 3.804 ------ ----\n[Dimer] 9 0.0038464 -0.0005 6.80745e-02 -17.7850 1.902 3.052 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7428 3.860 ------ ----\n[Dimer] 10 0.0011593 -0.0005 7.17333e-03 -17.7428 1.930 3.104 0\n\n 92 1.9614e-02 3.8861e-01 9 0.7315\nSaddle point search started from reactant with energy -53.30168914794922 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.0598 6.287 7.32 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.3476 5.196 4.10 2\n[Dimer] 1 0.0023570 -0.0002 1.40543e-01 -17.3476 5.196 4.103 2\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.3764 4.957 ------ ----\n[Dimer] 2 0.0008629 -0.0002 1.26085e-01 -17.3764 2.479 3.814 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4270 5.068 ------ ----\n[Dimer] 3 0.0083748 0.0002 1.30180e-01 -17.4270 2.534 4.137 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6455 3.248 2.21 1\n[Dimer] 4 0.0024208 0.0003 1.35435e-01 -17.6455 3.248 2.213 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8527 4.989 ------ ----\n[Dimer] 5 0.0013833 0.0002 1.27117e-02 -17.8527 2.494 2.919 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7320 5.114 ------ ----\n[Dimer] 6 0.0001255 0.0002 1.09239e-02 -17.7320 2.557 4.103 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7360 5.128 ------ ----\n[Dimer] 7 0.0009371 0.0002 1.71921e-02 -17.7360 2.564 4.113 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7500 5.085 ------ ----\n[Dimer] 8 0.0012917 0.0002 1.82971e-02 -17.7500 2.542 4.076 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7583 5.069 ------ ----\n[Dimer] 9 0.0019536 0.0002 2.80647e-02 -17.7583 2.534 4.061 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7426 5.108 ------ ----\n[Dimer] 10 0.0008442 0.0002 9.21547e-03 -17.7426 2.554 4.095 0\n\n 93 1.5317e-02 3.4760e-01 9 0.7316\nSaddle point search started from reactant with energy -53.30158615112305 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.3730 6.192 6.29 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.6548 4.425 4.48 2\n[Dimer] 1 0.0020058 -0.0002 9.74716e-02 -17.6548 4.425 4.478 2\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6130 3.848 ------ ----\n[Dimer] 2 0.0005768 -0.0002 8.85026e-02 -17.6130 1.924 2.912 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6225 3.946 ------ ----\n[Dimer] 3 0.0054785 0.0000 8.81213e-02 -17.6225 1.973 3.194 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7525 4.895 ------ ----\n[Dimer] 4 0.0022256 0.0000 9.61552e-02 -17.7525 2.448 3.925 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6771 5.210 ------ ----\n[Dimer] 5 0.0014405 0.0000 6.59491e-03 -17.6771 2.605 4.192 0\n\n 94 8.3856e-03 6.2215e-02 9 0.7315\nSaddle point search started from reactant with energy -53.30162048339844 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4824 6.078 5.46 1\n[IDimerRot] ----- --------- ---------- ------------------ -17.7486 4.314 4.33 2\n[Dimer] 1 0.0004226 0.0000 2.99555e-02 -17.7486 4.314 4.329 2\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7399 4.797 ------ ----\n[Dimer] 2 0.0001944 0.0000 2.66313e-02 -17.7399 2.398 3.384 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7434 4.805 ------ ----\n[Dimer] 3 0.0018740 0.0000 3.60346e-02 -17.7434 2.402 3.855 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7610 4.885 ------ ----\n[Dimer] 4 0.0006778 0.0000 2.88690e-02 -17.7610 2.443 3.915 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7902 4.838 ------ ----\n[Dimer] 5 0.0003530 0.0000 8.00247e-03 -17.7902 2.419 3.872 0\n\n 95 9.7215e-03 3.9727e-02 9 0.7315\nSaddle point search started from reactant with energy -53.301666259765625 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6301 5.992 4.68 1\n[Dimer] 1 0.0008308 0.0000 4.31939e-02 -17.6301 5.992 4.683 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7949 4.222 3.94 1\n[Dimer] 2 0.0006057 0.0000 2.56175e-02 -17.7949 4.222 3.940 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7753 5.421 ------ ----\n[Dimer] 3 0.0007222 0.0000 1.04089e-02 -17.7753 2.710 3.657 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7751 5.432 ------ ----\n[Dimer] 4 0.0002352 0.0000 9.67876e-03 -17.7751 2.716 4.344 0\n\n 96 1.4725e-02 1.3620e-01 9 0.7314\nSaddle point search started from reactant with energy -53.301727294921875 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6869 5.632 4.13 1\n[Dimer] 1 0.0008875 -0.0000 5.56824e-02 -17.6869 5.632 4.131 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8646 3.984 3.37 1\n[Dimer] 2 0.0002931 -0.0000 4.68744e-02 -17.8646 3.984 3.372 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8458 5.347 ------ ----\n[Dimer] 3 0.0026252 0.0000 1.68900e-02 -17.8458 2.674 3.515 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7748 5.574 ------ ----\n[Dimer] 4 0.0004773 0.0000 3.04919e-02 -17.7748 2.787 4.456 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7509 5.626 ------ ----\n[Dimer] 5 0.0002662 0.0000 6.84089e-03 -17.7509 2.813 4.502 0\n\n 97 1.5321e-02 1.9075e-01 9 0.7315\nSaddle point search started from reactant with energy -53.30166244506836 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6996 5.464 3.64 1\n[Dimer] 1 0.0010913 -0.0001 3.68170e-02 -17.6996 5.464 3.638 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7153 7.561 ------ ----\n[Dimer] 2 0.0001918 -0.0001 3.21716e-02 -17.7153 3.781 4.915 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.8275 3.769 1.76 1\n[Dimer] 3 0.0019253 -0.0001 2.34463e-02 -17.8275 3.769 1.761 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7860 5.043 ------ ----\n[Dimer] 4 0.0006522 -0.0001 2.65289e-02 -17.7860 2.521 3.462 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7663 5.051 ------ ----\n[Dimer] 5 0.0001817 -0.0001 1.31186e-02 -17.7663 2.526 4.045 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7735 5.006 ------ ----\n[Dimer] 6 0.0003298 -0.0001 1.05159e-02 -17.7735 2.503 4.008 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7722 4.999 ------ ----\n[Dimer] 7 0.0009807 -0.0001 1.58063e-02 -17.7722 2.500 4.003 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7552 4.966 ------ ----\n[Dimer] 8 0.0012014 -0.0001 1.43548e-02 -17.7552 2.483 3.981 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7357 4.953 ------ ----\n[Dimer] 9 0.0018016 -0.0001 1.41579e-02 -17.7357 2.477 3.975 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7279 4.929 ------ ----\n[Dimer] 10 0.0030930 -0.0001 2.49768e-02 -17.7279 2.465 3.957 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7040 4.845 ------ ----\n[Dimer] 11 0.0005617 -0.0001 9.34264e-03 -17.7040 2.423 3.896 0\n\n 98 3.5307e-02 5.4306e-01 9 0.7322\n 99 9.0356e-03 1.1793e-01 9 0.7314\nSaddle point search started from reactant with energy -53.30179214477539 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4503 5.120 4.17 1\n[Dimer] 1 0.0007901 -0.0000 4.86714e-02 -17.4503 5.120 4.165 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6160 4.069 3.14 1\n[Dimer] 2 0.0002834 -0.0000 4.05183e-02 -17.6160 4.069 3.139 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6043 4.642 ------ ----\n[Dimer] 3 0.0023568 -0.0001 2.70710e-02 -17.6043 2.321 3.206 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6054 4.717 ------ ----\n[Dimer] 4 0.0014455 -0.0001 3.32475e-02 -17.6054 2.359 3.815 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.6293 4.736 ------ ----\n[Dimer] 5 0.0026278 -0.0001 2.30855e-02 -17.6293 2.368 3.825 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7184 4.656 ------ ----\n[Dimer] 6 0.0029157 -0.0001 2.65869e-02 -17.7184 2.328 3.742 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7489 4.645 ------ ----\n[Dimer] 7 0.0017354 -0.0001 3.02649e-02 -17.7489 2.322 3.727 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7806 4.607 ------ ----\n[Dimer] 8 0.0006492 -0.0001 9.49622e-03 -17.7806 2.304 3.691 0\n\n 100 8.4827e-03 4.8794e-02 9 0.7312\nSaddle point search started from reactant with energy -53.30195617675781 eV.\n[Dimer] Step Step Size Delta E ||Force|| Curvature Torque Angle Rots\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.4715 5.177 4.19 1\n[Dimer] 1 0.0013349 0.0000 5.08536e-02 -17.4715 5.177 4.192 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7672 3.859 3.10 1\n[Dimer] 2 0.0009469 0.0000 2.10970e-02 -17.7672 3.859 3.103 1\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7763 4.593 ------ ----\n[Dimer] 3 0.0005979 0.0000 1.59646e-02 -17.7763 2.297 3.040 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7702 4.552 ------ ----\n[Dimer] 4 0.0002152 0.0000 1.52322e-02 -17.7702 2.276 3.650 0\n\n[IDimerRot] ----- --------- ---------- ------------------ -17.7617 4.644 ------ ----\n[Dimer] 5 0.0000849 0.0000 8.17779e-03 -17.7617 2.322 3.724 0\n\nFinal state: \nNudged elastic band, too many iterations.\ntiming information:\nreal 28.322 seconds\nuser 54.499 seconds\nsys 0.865 seconds\n')
Visual interpretation¶
rgpycrumbs is a visualization toolkit designed to bridge the gap between raw computational output and physical intuition, mapping high-dimensional NEB trajectories onto interpretable 1D energy profiles and 2D RMSD landscapes.
We check both the standard 1D profile against the path reaction coordinate, or the distance between intermediate images:
# Clean env to prevent backend conflicts in notebooks
os.environ.pop("MPLBACKEND", None)
oneDprof_oxad = [
sys.executable,
"-m",
"rgpycrumbs.cli",
"eon",
"plt_neb",
"--theme",
"ruhi",
"--con-file",
"neb.con",
"--plot-structures",
"crit_points",
"--facecolor",
"white",
"--plot-type",
"profile",
"--ase-rotation=90x,0y,0z",
"--title=NEB Path Optimization",
"--output-file",
"1D_oxad.png",
]
# Run the 1D plotting command using the helper
run_command_or_exit(oneDprof_oxad, capture=False, timeout=60)
# Display the result
img = mpimg.imread("1D_oxad.png")
plt.figure(figsize=(8, 6))
plt.imshow(img)
plt.axis("off")
plt.show()
--> Dispatching to: uv run /home/runner/work/atomistic-cookbook/atomistic-cookbook/.nox/eon-pet-neb/lib/python3.13/site-packages/rgpycrumbs/eon/plt_neb.py --theme ruhi --con-file neb.con --plot-structures crit_points --facecolor white --plot-type profile --ase-rotation=90x,0y,0z --title=NEB Path Optimization --output-file 1D_oxad.png
Downloading matplotlib (8.3MiB)
Downloading kiwisolver (1.4MiB)
Downloading numpy (15.8MiB)
Downloading pygments (1.2MiB)
Downloading pillow (6.7MiB)
Downloading scipy (34.0MiB)
Downloading fonttools (4.7MiB)
Downloading polars-runtime-32 (42.4MiB)
Downloading ase (2.8MiB)
Downloaded kiwisolver
Downloaded pygments
Downloaded fonttools
Downloaded pillow
Downloaded matplotlib
Downloaded numpy
Downloaded ase
Downloaded polars-runtime-32
Downloaded scipy
Installed 21 packages in 29ms
[12/10/25 16:29:50] INFO INFO - Overriding theme facecolor with white
INFO INFO - Setting global rcParams for ruhi theme
WARNING WARNING - Font 'Atkinson Hyperlegible' not found.
Falling back to 'sans-serif'.
WARNING WARNING - For custom fonts to work, they must be
installed on your system and recognized by
matplotlib.
WARNING WARNING - You may need to clear the matplotlib
cache: /home/runner/.cache/matplotlib
INFO INFO - Using figsize: width=10.00", height=7.00"
INFO INFO - Applying axis-specific theme for ruhi
INFO INFO - Reading structures from neb.con
INFO INFO - Loaded 12 structures.
INFO INFO - Searching for files with pattern:
'neb_*.dat'
INFO INFO - Found 101 file(s).
[12/10/25 16:29:51] INFO INFO - Saving plot to 1D_oxad.png
Also, the PES 2D landscape profile as a function of the RMSD [3] which shows the relative distance between the endpoints as the optimization takes place:
twoDprof_oxad = [
sys.executable,
"-m",
"rgpycrumbs.cli",
"eon",
"plt_neb",
"--theme",
"ruhi",
"--con-file",
"neb.con",
"--plot-structures",
"crit_points",
"--facecolor",
"white",
"--rc-mode",
"path",
"--ase-rotation=90x,0y,0z",
"--title",
"NEB-RMSD Surface",
"--fontsize-base",
"16",
"--landscape-mode",
"surface",
"--landscape-path",
"all",
"--plot-type",
"landscape",
"--show-pts",
"--surface-type",
"rbf",
"--output-file",
"2D_oxad.png",
]
# Run the 2D plotting command using the helper
run_command_or_exit(twoDprof_oxad, capture=False, timeout=60)
# Display the result
img = mpimg.imread("2D_oxad.png")
plt.figure(figsize=(8, 6))
plt.imshow(img)
plt.axis("off")
plt.show()
--> Dispatching to: uv run /home/runner/work/atomistic-cookbook/atomistic-cookbook/.nox/eon-pet-neb/lib/python3.13/site-packages/rgpycrumbs/eon/plt_neb.py --theme ruhi --con-file neb.con --plot-structures crit_points --facecolor white --rc-mode path --ase-rotation=90x,0y,0z --title NEB-RMSD Surface --fontsize-base 16 --landscape-mode surface --landscape-path all --plot-type landscape --show-pts --surface-type rbf --output-file 2D_oxad.png
[12/10/25 16:29:53] INFO INFO - Overriding theme font size with 16
INFO INFO - Overriding theme facecolor with white
INFO INFO - Setting global rcParams for ruhi theme
WARNING WARNING - Font 'Atkinson Hyperlegible' not found.
Falling back to 'sans-serif'.
WARNING WARNING - For custom fonts to work, they must be
installed on your system and recognized by
matplotlib.
WARNING WARNING - You may need to clear the matplotlib
cache: /home/runner/.cache/matplotlib
INFO INFO - Using figsize: width=10.00", height=7.00"
INFO INFO - Applying axis-specific theme for ruhi
INFO INFO - Reading structures from neb.con
INFO INFO - Loaded 12 structures.
INFO INFO - Found 101 numerically-matched dat/con step
indices.
INFO INFO - Computing Landscape data...
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
[12/10/25 16:29:54] INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Calculating landscape coordinates using
ira.match...
INFO INFO - Landscape coordinate calculation complete.
INFO INFO - Saving Landscape cache to
.neb_landscape.parquet...
INFO INFO - Generating interpolated RBF 2D surface...
INFO INFO - Saving plot to 2D_oxad.png
Where each black dot represents a configuration at which the energy and forces are calculated during the NEB optimization run, and the RMSD is the “permutation corrected” distance from the reactant and product. The interpolated relative energy surface is generated by a Gaussian approximation with a radial basis function kernel interpolation over the energy and forces of all the black dots. See [3, Chapter 4] for more details of this visualization.
Relaxing the endpoints with EON¶
In this final part we come back to an essential point of performing NEB calculations, and that is the relaxation of the initial states. In the tutorials above we used directly relaxed structures, and here we are demonstrating how these can be relaxed. We first load structures which are not relaxed.
reactant = aseio.read("data/reactant.con")
product = aseio.read("data/product.con")
# For compatibility with EON, we also need to provide
# a unit cell
def center_cell(atoms):
atoms.set_cell([20, 20, 20])
atoms.pbc = True
atoms.center()
return atoms
reactant = center_cell(reactant)
product = center_cell(product)
The resulting reactant has a larger box:
fig, (ax1, ax2) = plt.subplots(1, 2)
plot_atoms(reactant, ax1, rotation=("-90x,0y,0z"))
plot_atoms(product, ax2, rotation=("-90x,0y,0z"))
ax1.text(0.3, -1, "reactant")
ax2.text(0.3, -1, "product")
ax1.set_axis_off()
ax2.set_axis_off()
def write_eon_min_config(run_dir: Path, model_path: Path):
"""
Writes the config.ini file for an EON NEB job,
using the user's provided template.
"""
config_content = f"""[Main]
job=minimization
random_seed = 706253457
[Potential]
potential = Metatomic
[Metatomic]
model_path = {Path(str(model_path).replace("ckpt", "pt")).absolute()}
[Optimizer]
max_iterations = 2000
opt_method = lbfgs
max_move = 0.1
converged_force = 0.01
"""
config_path = run_dir / "config.ini"
with open(config_path, "w") as f:
f.write(config_content)
print(f"Wrote EON NEB config to '{config_path}'")
# Reactant setup
dir_reactant = Path("min_reactant")
dir_reactant.mkdir(exist_ok=True)
aseio.write(dir_reactant / "pos.con", reactant)
write_eon_min_config(dir_reactant, Path(fname).absolute())
# Product setup
dir_product = Path("min_product")
dir_product.mkdir(exist_ok=True)
aseio.write(dir_product / "pos.con", product)
write_eon_min_config(dir_product, Path(fname).absolute())
Wrote EON NEB config to 'min_reactant/config.ini'
Wrote EON NEB config to 'min_product/config.ini'
Run the minimization¶
The ‘eonclient’ will use the correct configuration within the folder.
@contextlib.contextmanager
def work_in_dir(path: Path):
"""
Context manager to safely change directory and return to previous
one afterwards. Crucial for notebooks to avoid path drift.
"""
prev_cwd = Path.cwd()
os.chdir(path)
try:
yield
finally:
os.chdir(prev_cwd)
with work_in_dir(dir_reactant):
run_command_or_exit(["eonclient"], capture=True, timeout=300)
with work_in_dir(dir_product):
run_command_or_exit(["eonclient"], capture=True, timeout=300)
EON Client
VERSION: f01406f
BUILD DATE: Tue Dec 02 01:53:51 AM GMT 2025
OS: linux
Arch: x86_64
Hostname: runnervmoqczp
PID: 2570
DIR: /home/runner/work/atomistic-cookbook/atomistic-cookbook/examples/eon-pet-neb/min_reactant
Loading parameter file config.ini
* [Main] job: minimization
* [Main] random_seed: 706253457
* [Potential] potential: Metatomic
* [Optimizer] opt_method: lbfgs
* [Optimizer] converged_force: 0.01
* [Optimizer] max_iterations: 2000
* [Optimizer] max_move: 0.1
* [Metatomic] model_path: /home/runner/work/atomistic-cookbook/atomistic-cookbook/examples/eon-pet-neb/pet-mad-v1.0.2.pt
Beginning minimization of pos.con
[Matter] Iter Step size ||Force|| Energy
[Matter] 0 0.00000e+00 3.35994e+00 -54.11823
[Matter] 1 2.13300e-02 2.13510e+00 -54.16295
[Matter] 2 8.64217e-03 5.69126e-01 -54.17617
[Matter] 3 2.05505e-03 4.85098e-01 -54.17807
[Matter] 4 5.61483e-03 3.34296e-01 -54.18124
[Matter] 5 3.06752e-03 1.15048e-01 -54.18205
[Matter] 6 1.70797e-03 1.14860e-01 -54.18222
[Matter] 7 2.58184e-03 1.65259e-01 -54.18243
[Matter] 8 6.45005e-03 1.83993e-01 -54.18282
[Matter] 9 1.08150e-02 1.73501e-01 -54.18347
[Matter] 10 5.57214e-02 3.04151e-01 -54.18568
[Matter] 11 8.22242e-02 4.28157e-01 -54.18802
[Matter] 12 3.75054e-02 2.71154e-01 -54.19065
[Matter] 13 7.89329e-02 2.31768e-01 -54.19278
[Matter] 14 9.24889e-03 2.62950e-01 -54.19288
[Matter] 15 5.55577e-03 7.74731e-02 -54.19317
[Matter] 16 1.04434e-03 5.50555e-02 -54.19323
[Matter] 17 7.39721e-03 6.94049e-02 -54.19339
[Matter] 18 3.39791e-03 6.10213e-02 -54.19351
[Matter] 19 4.05978e-03 6.29989e-02 -54.19357
[Matter] 20 1.18729e-03 2.24259e-02 -54.19360
[Matter] 21 9.99609e-04 2.10560e-02 -54.19361
[Matter] 22 1.18117e-03 2.88133e-02 -54.19363
[Matter] 23 3.44700e-03 5.17056e-02 -54.19365
[Matter] 24 2.21968e-03 3.24301e-02 -54.19368
[Matter] 25 1.96420e-03 2.98051e-02 -54.19369
[Matter] 26 5.39770e-03 4.47915e-02 -54.19376
[Matter] 27 2.08058e-02 6.42344e-02 -54.19399
[Matter] 28 4.42303e-02 2.10412e-01 -54.19418
[Matter] 29 1.36149e-02 1.13144e-01 -54.19456
[Matter] 30 3.26525e-02 5.34401e-02 -54.19494
[Matter] 31 3.43590e-02 1.59670e-01 -54.19478
[Matter] 32 2.39536e-02 6.66176e-02 -54.19505
[Matter] 33 4.71602e-03 5.01543e-02 -54.19508
[Matter] 34 1.31916e-03 2.87957e-02 -54.19510
[Matter] 35 3.46048e-03 3.00308e-02 -54.19514
[Matter] 36 1.31145e-03 3.16942e-02 -54.19517
[Matter] 37 1.00937e-02 1.68636e-01 -54.19505
[Matter] 38 6.99502e-03 1.89155e-02 -54.19521
[Matter] 39 6.15134e-04 1.22553e-02 -54.19521
[Matter] 40 6.94035e-04 1.44800e-02 -54.19522
[Matter] 41 1.17171e-03 3.44055e-02 -54.19522
[Matter] 42 8.91710e-04 2.54493e-02 -54.19523
[Matter] 43 8.40842e-03 5.23743e-02 -54.19530
[Matter] 44 1.51423e-02 1.12806e-01 -54.19540
[Matter] 45 3.00541e-02 2.01424e-01 -54.19556
[Matter] 46 1.99334e-02 2.06066e-01 -54.19582
[Matter] 47 1.47601e-02 2.86510e+00 -54.16262
[Matter] 48 1.35789e-02 6.35039e-02 -54.19606
[Matter] 49 2.97498e-04 6.39530e-02 -54.19608
[Matter] 50 2.70639e-04 5.64622e-02 -54.19609
[Matter] 51 2.04382e-03 2.41977e-02 -54.19614
[Matter] 52 4.57460e-04 3.60070e-02 -54.19614
[Matter] 53 2.02473e-04 1.44476e-02 -54.19615
[Matter] 54 1.45797e-04 1.09404e-02 -54.19615
[Matter] 55 8.63914e-04 2.26838e-02 -54.19617
[Matter] 56 1.30279e-03 2.99261e-02 -54.19618
[Matter] 57 3.12584e-03 3.62484e-02 -54.19621
[Matter] 58 1.03842e-02 1.01008e-01 -54.19624
[Matter] 59 5.60303e-03 4.54559e-02 -54.19632
[Matter] 60 6.42112e-03 3.25697e-02 -54.19638
[Matter] 61 1.29644e-02 5.63956e-02 -54.19647
[Matter] 62 2.13468e-02 1.35651e-01 -54.19655
[Matter] 63 1.94048e-02 8.93571e-02 -54.19668
[Matter] 64 1.91308e-02 9.22966e-02 -54.19687
[Matter] 65 1.68858e-02 1.03039e-01 -54.19702
[Matter] 66 1.10646e-02 7.96517e-02 -54.19716
[Matter] 67 4.78943e-02 2.97596e-01 -54.19691
[Matter] 68 3.86362e-02 3.15286e-02 -54.19726
[Matter] 69 6.67822e-03 3.09921e-02 -54.19729
[Matter] 70 3.61675e-02 1.63606e-01 -54.19728
[Matter] 71 8.01655e-03 8.70214e-02 -54.19736
[Matter] 72 1.40798e-03 2.75092e-02 -54.19739
[Matter] 73 7.64475e-03 3.19429e-02 -54.19741
[Matter] 74 2.53417e-02 4.29464e-02 -54.19745
[Matter] 75 2.07332e-02 4.59619e-02 -54.19753
[Matter] 76 3.00587e-02 8.62458e-02 -54.19761
[Matter] 77 5.22281e-03 1.68464e-01 -54.19762
[Matter] 78 3.23256e-02 8.30238e-02 -54.19778
[Matter] 79 3.34794e-03 5.57963e-02 -54.19785
[Matter] 80 3.15576e-02 7.91883e-02 -54.19806
[Matter] 81 3.78017e-02 1.38264e-01 -54.19843
[Matter] 82 4.23751e-02 1.76429e-01 -54.19909
[Matter] 83 2.43435e-03 3.83794e-01 -54.19906
[Matter] 84 1.08371e-03 1.13490e-01 -54.19940
[Matter] 85 2.28149e-04 8.69504e-02 -54.19946
[Matter] 86 1.93757e-03 6.26631e-02 -54.19968
[Matter] 87 2.20366e-03 5.41573e-02 -54.19980
[Matter] 88 3.05952e-03 1.04606e-01 -54.19992
[Matter] 89 2.60609e-03 7.02966e-02 -54.20004
[Matter] 90 1.67229e-03 4.72268e-02 -54.20010
[Matter] 91 1.14044e-03 3.94763e-02 -54.20014
[Matter] 92 1.95462e-03 3.44234e-02 -54.20018
[Matter] 93 5.02763e-03 4.69987e-02 -54.20027
[Matter] 94 9.70108e-03 6.45405e-02 -54.20042
[Matter] 95 2.15335e-02 1.33208e-01 -54.20064
[Matter] 96 2.19392e-02 1.10219e-01 -54.20095
[Matter] 97 1.35240e-02 7.84058e-02 -54.20131
[Matter] 98 2.17364e-02 6.55860e-02 -54.20161
[Matter] 99 4.37914e-03 8.01132e-02 -54.20166
[Matter] 100 2.98891e-03 3.73772e-02 -54.20171
[Matter] 101 6.61838e-03 2.84772e-02 -54.20176
[Matter] 102 5.65942e-03 3.29427e-02 -54.20179
[Matter] 103 1.05077e-02 3.81033e-02 -54.20186
[Matter] 104 1.62336e-02 7.25631e-02 -54.20200
[Matter] 105 3.01170e-02 1.14619e-01 -54.20216
[Matter] 106 1.66715e-02 9.21741e-02 -54.20236
[Matter] 107 1.06290e-02 8.82887e-02 -54.20256
[Matter] 108 2.87583e-02 1.24912e-01 -54.20274
[Matter] 109 3.35153e-03 1.32399e-01 -54.20298
[Matter] 110 3.27277e-02 1.78752e-01 -54.20332
[Matter] 111 3.32625e-03 1.04020e-01 -54.20353
[Matter] 112 2.53603e-02 6.67521e-02 -54.20376
[Matter] 113 1.55836e-02 1.21973e-01 -54.20386
[Matter] 114 1.15801e-02 6.37750e-02 -54.20393
[Matter] 115 3.68519e-03 6.09228e-02 -54.20399
[Matter] 116 7.57257e-03 5.59562e-02 -54.20407
[Matter] 117 7.52381e-03 5.68448e-02 -54.20414
[Matter] 118 1.09631e-02 6.80684e-02 -54.20423
[Matter] 119 4.58184e-03 6.22465e-02 -54.20433
[Matter] 120 9.90518e-03 7.90822e-02 -54.20450
[Matter] 121 2.40655e-02 1.42548e-01 -54.20491
[Matter] 122 4.93923e-02 2.19583e-01 -54.20581
[Matter] 123 5.40878e-02 3.27944e-01 -54.20691
[Matter] 124 6.39374e-03 7.15219e-01 -54.20540
[Matter] 125 4.19793e-03 1.16278e-01 -54.20789
[Matter] 126 7.03100e-04 9.06617e-02 -54.20803
[Matter] 127 3.76818e-03 2.89006e-01 -54.20815
[Matter] 128 1.14111e-03 9.95391e-02 -54.20832
[Matter] 129 4.59101e-04 4.93038e-02 -54.20834
[Matter] 130 2.19537e-04 4.05293e-02 -54.20836
[Matter] 131 1.39114e-03 6.68066e-02 -54.20843
[Matter] 132 1.67480e-03 7.07079e-02 -54.20848
[Matter] 133 3.62236e-03 1.36642e-01 -54.20852
[Matter] 134 1.04774e-03 4.47788e-02 -54.20859
[Matter] 135 5.49334e-04 2.68529e-02 -54.20860
[Matter] 136 5.04334e-04 2.99504e-02 -54.20862
[Matter] 137 1.53520e-03 5.65289e-02 -54.20865
[Matter] 138 3.05052e-03 7.80733e-02 -54.20869
[Matter] 139 5.43370e-03 8.06951e-02 -54.20876
[Matter] 140 1.51532e-02 1.16888e-01 -54.20895
[Matter] 141 3.70510e-02 1.85510e-01 -54.20929
[Matter] 142 2.74276e-02 2.08571e-01 -54.20972
[Matter] 143 3.59814e-02 1.31082e-01 -54.21060
[Matter] 144 3.65769e-02 1.93250e-01 -54.21094
[Matter] 145 2.16360e-02 1.53911e-01 -54.21121
[Matter] 146 2.98955e-02 1.91438e-01 -54.21148
[Matter] 147 5.90424e-03 8.54403e-02 -54.21165
[Matter] 148 1.53044e-02 7.26897e-02 -54.21180
[Matter] 149 3.30745e-02 1.23511e-01 -54.21203
[Matter] 150 5.20897e-02 4.71689e-01 -54.21162
[Matter] 151 3.13438e-02 2.30814e-01 -54.21222
[Matter] 152 1.61378e-02 2.41111e-01 -54.21227
[Matter] 153 3.70307e-03 5.59423e-02 -54.21248
[Matter] 154 1.23226e-02 7.64374e-02 -54.21257
[Matter] 155 2.71652e-02 9.98032e-02 -54.21275
[Matter] 156 3.32433e-02 2.04274e-01 -54.21303
[Matter] 157 2.17996e-02 1.27652e-01 -54.21339
[Matter] 158 2.07666e-02 1.48198e-01 -54.21346
[Matter] 159 1.00259e-02 1.73915e-01 -54.21350
[Matter] 160 3.94405e-03 6.95972e-02 -54.21366
[Matter] 161 7.19456e-03 5.86682e-02 -54.21374
[Matter] 162 1.66473e-02 1.00232e-01 -54.21392
[Matter] 163 9.74421e-03 9.52959e-02 -54.21410
[Matter] 164 3.45415e-02 1.43938e-01 -54.21429
[Matter] 165 7.34192e-03 6.17666e-02 -54.21444
[Matter] 166 1.14150e-02 4.52182e-02 -54.21453
[Matter] 167 1.85403e-02 6.25672e-02 -54.21461
[Matter] 168 6.49744e-03 5.99286e-02 -54.21467
[Matter] 169 1.34825e-02 5.74583e-02 -54.21478
[Matter] 170 3.53771e-02 1.36717e-01 -54.21478
[Matter] 171 1.93409e-02 3.93742e-02 -54.21488
[Matter] 172 4.50900e-03 2.82102e-02 -54.21490
[Matter] 173 7.82294e-03 1.69509e-02 -54.21492
[Matter] 174 8.04884e-04 1.57452e-02 -54.21493
[Matter] 175 2.57630e-03 3.76793e-02 -54.21495
[Matter] 176 2.12580e-03 3.75752e-02 -54.21497
[Matter] 177 3.40618e-03 6.29038e-02 -54.21502
[Matter] 178 1.29041e-02 1.08045e-01 -54.21517
[Matter] 179 1.79418e-02 1.05324e-01 -54.21534
[Matter] 180 4.22479e-02 4.22966e-01 -54.21499
[Matter] 181 8.59048e-03 8.82492e-02 -54.21562
[Matter] 182 2.66002e-03 4.50519e-02 -54.21569
[Matter] 183 7.75794e-03 4.28520e-02 -54.21577
[Matter] 184 5.78282e-03 4.96211e-02 -54.21581
[Matter] 185 1.34927e-02 8.15384e-02 -54.21581
[Matter] 186 5.93939e-03 2.26017e-02 -54.21585
[Matter] 187 3.36634e-03 3.19088e-02 -54.21585
[Matter] 188 6.36630e-03 2.18895e-02 -54.21587
[Matter] 189 7.25423e-03 3.11477e-02 -54.21588
[Matter] 190 5.19965e-03 1.16910e-02 -54.21589
[Matter] 191 8.49927e-04 9.22649e-03 -54.21589
Minimization converged within tolerence
Saving result to min.con
Final Energy: -54.21588897705078
timing information:
real 7.545 seconds
user 13.577 seconds
sys 0.205 seconds
EON Client
VERSION: f01406f
BUILD DATE: Tue Dec 02 01:53:51 AM GMT 2025
OS: linux
Arch: x86_64
Hostname: runnervmoqczp
PID: 2574
DIR: /home/runner/work/atomistic-cookbook/atomistic-cookbook/examples/eon-pet-neb/min_product
Loading parameter file config.ini
* [Main] job: minimization
* [Main] random_seed: 706253457
* [Potential] potential: Metatomic
* [Optimizer] opt_method: lbfgs
* [Optimizer] converged_force: 0.01
* [Optimizer] max_iterations: 2000
* [Optimizer] max_move: 0.1
* [Metatomic] model_path: /home/runner/work/atomistic-cookbook/atomistic-cookbook/examples/eon-pet-neb/pet-mad-v1.0.2.pt
Beginning minimization of pos.con
[Matter] Iter Step size ||Force|| Energy
[Matter] 0 0.00000e+00 2.03521e+00 -54.72154
[Matter] 1 2.30092e-02 1.54286e+00 -54.75799
[Matter] 2 1.14317e-02 4.98111e-01 -54.77272
[Matter] 3 4.37960e-03 2.89456e-01 -54.77559
[Matter] 4 4.27196e-03 1.55479e-01 -54.77702
[Matter] 5 3.15189e-03 1.01303e-01 -54.77753
[Matter] 6 3.79740e-03 9.20883e-02 -54.77778
[Matter] 7 1.42435e-03 4.40399e-02 -54.77789
[Matter] 8 4.17118e-04 2.66551e-02 -54.77792
[Matter] 9 7.49504e-04 2.74143e-02 -54.77794
[Matter] 10 1.29938e-03 3.68717e-02 -54.77798
[Matter] 11 2.77309e-03 7.09546e-02 -54.77800
[Matter] 12 1.30622e-03 3.65371e-02 -54.77805
[Matter] 13 1.62314e-03 3.98186e-02 -54.77809
[Matter] 14 3.87701e-03 7.53303e-02 -54.77817
[Matter] 15 1.09979e-02 1.52479e-01 -54.77837
[Matter] 16 4.37323e-02 4.02710e-01 -54.77907
[Matter] 17 6.13801e-02 7.21846e-01 -54.77945
[Matter] 18 6.09176e-03 5.70980e-01 -54.78160
[Matter] 19 8.00970e-02 4.25474e-01 -54.79069
[Matter] 20 1.27584e-02 1.53876e+00 -54.78309
[Matter] 21 9.32641e-03 2.14231e-01 -54.79369
[Matter] 22 1.17776e-03 1.30921e-01 -54.79409
[Matter] 23 3.34097e-03 9.50682e-02 -54.79457
[Matter] 24 2.48103e-03 9.23926e-02 -54.79483
[Matter] 25 8.04328e-03 1.83141e-01 -54.79518
[Matter] 26 3.74425e-03 1.01847e-01 -54.79548
[Matter] 27 1.22234e-03 8.89367e-02 -54.79564
[Matter] 28 9.94251e-03 1.30721e-01 -54.79638
[Matter] 29 2.18906e-02 3.25293e-01 -54.79702
[Matter] 30 9.10332e-03 1.85220e-01 -54.79797
[Matter] 31 1.64896e-02 1.53094e-01 -54.79901
[Matter] 32 1.10536e-02 1.32567e-01 -54.79959
[Matter] 33 1.53104e-02 1.07156e-01 -54.80007
[Matter] 34 4.80920e-03 1.88434e-01 -54.79992
[Matter] 35 3.30005e-03 2.84513e-02 -54.80020
[Matter] 36 6.26490e-04 2.36519e-02 -54.80021
[Matter] 37 2.71879e-03 3.31060e-02 -54.80026
[Matter] 38 1.47254e-03 3.82295e-02 -54.80027
[Matter] 39 4.91832e-04 1.32723e-02 -54.80028
[Matter] 40 2.80982e-04 9.26102e-03 -54.80028
Minimization converged within tolerence
Saving result to min.con
Final Energy: -54.8002815246582
timing information:
real 3.095 seconds
user 4.837 seconds
sys 0.132 seconds
Additionally, the relative ordering must be preserved, for which we use IRA [4].
reactant = aseio.read(dir_reactant / "min.con")
product = aseio.read(dir_product / "min.con")
ira = ira_mod.IRA()
# Default value
kmax_factor = 1.8
nat1 = len(reactant)
typ1 = reactant.get_atomic_numbers()
coords1 = reactant.get_positions()
nat2 = len(product)
typ2 = product.get_atomic_numbers()
coords2 = product.get_positions()
print("Running ira.match to find rotation, translation, AND permutation...")
# r = rotation, t = translation, p = permutation, hd = Hausdorff distance
r, t, p, hd = ira.match(nat1, typ1, coords1, nat2, typ2, coords2, kmax_factor)
print(f"Matching complete. Hausdorff Distance (hd) = {hd:.6f} Angstrom")
# Apply rotation (r) and translation (t) to the original product coordinates
# This aligns the product's orientation to the reactant's
coords2_aligned = (coords2 @ r.T) + t
# Apply the permutation (p)
# This re-orders the aligned product atoms to match the reactant's atom order
# p[i] = j means reactant atom 'i' matches product atom 'j'
# So, the new coordinate array's i-th element should be coords2_aligned[j]
coords2_aligned_permuted = coords2_aligned[p]
# Save the new aligned-and-permuted structure
# CRUCIAL: Use chemical symbols from the reactant,
# because we have now permuted the product coordinates to match the reactant order.
product = reactant.copy()
product.positions = coords2_aligned_permuted
Running ira.match to find rotation, translation, AND permutation...
Matching complete. Hausdorff Distance (hd) = 1.390119 Angstrom
Finally we can visualize these with ASE.
view(reactant, viewer="x3d")
view(product, viewer="x3d")
fig, (ax1, ax2) = plt.subplots(1, 2)
plot_atoms(reactant, ax1, rotation=("-90x,0y,0z"))
plot_atoms(product, ax2, rotation=("-90x,0y,0z"))
ax1.text(0.3, -1, "reactant")
ax2.text(0.3, -1, "product")
ax1.set_axis_off()
ax2.set_axis_off()
References¶
Mazitov, A.; Bigi, F.; Kellner, M.; Pegolo, P.; Tisi, D.; Fraux, G.; Pozdnyakov, S.; Loche, P.; Ceriotti, M. PET-MAD, a Universal Interatomic Potential for Advanced Materials Modeling. arXiv March 18, 2025. https://doi.org/10.48550/arXiv.2503.14118.
Bigi, F.; Abbott, J. W.; Loche, P.; Mazitov, A.; Tisi, D.; Langer, M. F.; Goscinski, A.; Pegolo, P.; Chong, S.; Goswami, R.; Chorna, S.; Kellner, M.; Ceriotti, M.; Fraux, G. Metatensor and Metatomic: Foundational Libraries for Interoperable Atomistic Machine Learning. arXiv August 21, 2025. https://doi.org/10.48550/arXiv.2508.15704.
Goswami, R. Efficient Exploration of Chemical Kinetics. arXiv October 24, 2025. https://doi.org/10.48550/arXiv.2510.21368.
Gunde, M.; Salles, N.; Hémeryck, A.; Martin-Samos, L. IRA: A Shape Matching Approach for Recognition and Comparison of Generic Atomic Patterns. J. Chem. Inf. Model. 2021, 61 (11), 5446–5457. https://doi.org/10.1021/acs.jcim.1c00567.
Smidstrup, S.; Pedersen, A.; Stokbro, K.; Jónsson, H. Improved Initial Guess for Minimum Energy Path Calculations. J. Chem. Phys. 2014, 140 (21), 214106. https://doi.org/10.1063/1.4878664.
Total running time of the script: (1 minutes 29.870 seconds)