[molpro-user] To understand the reason for difference betweeen two outputs for the same situation developed by MOLPRO

[Sham K]

Hello Users,

One simple possibility is to check nuclear repulsion energy to satisfy if
two geometries are the same.

There is another nice way to check if two structures (with same atoms and
the order in which they are listed is same) with different xyz data are
same. This comes from an elegant algorithm known as Kabsch algorithm - for
details of how it works see - https://en.wikipedia.org/wiki/Kabsch_algorithm

It calculates an RMSD between two structures. For the case here, using the
calculate_rmsd program from github, download it from
https://raw.githubusercontent.com/charnley/rmsd/master/rmsd/calculate_rmsd.py

sham at shree:~/junk$ ./calculate_rmsd.py a.xyz b.xyz
Normal RMSD: 0.305543115084
Kabsch RMSD: 7.14474667163e-06
Quater RMSD: 7.14474667163e-06

Showing that two structures are indeed the same.

-Sham
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20171116/07e4af8e/attachment.html>