[CP2K-user] [CP2K:20968] Wrong distance between geometrical centres and RMSD values reported in metadynLog

[Edwin Helbert Aponte Angarita]

Hi Albert,

You are right, the ratio between the distances measured on PyMOL and the 
COLVAR values reported in metadynLog is roughly 0.529, which is the 
conversion factor from Bohr and Angstrom. The PyMOL RMSD to COLVAR RMSD 
ratio tends towards 0.007 (see attached plot) but I can't figure out the 
units of RMSD in metadynLog.

Regards,
Edwin.

On 12/15/24 10:39 AM, Albert Makhmudov wrote:
> Hi Edwin,
> 
> To the best of my knowledge, by default CP2K outputs everything in a.u. but coordinates. Keeping that in mind, the distance between two geometrical centers is c. 9-10 bohr which is roughly 4-5 angstrom. I’m not sure about the RMSD values, though.
> 
> 
> Best,
> 
> Albert
> 
>> On Dec 14, 2024, at 21:27, Edwin Helbert Aponte Angarita <helbert2a at gmail.com> wrote:
>>
>> Dear All,
>>
>> I have been trying to calculate a free energy surface in the same way as presented in this tutorial:
>> https://www.cp2k.org/exercises:2020_uzh_acpc2:ex03
>> However, as collective variables I am specifying the distance between two geometrical centres and fitted RMSD to a reference structure. The calculation apparently runs fine but the values reported in metadynLog are wrong.
>>
>> I repeated the calculation using Clofenotane (DDT) as a test system so that I can post here my input files, which are attached. One of the CVs is the distance between the geometrical centre of the rings and the other is fitted RMSD to the initial structure. The distance reported in metadynLog is around 9angs while the one measured on PyMOL never goes above 5.5angs. RMSD values reported in metadynLog go up to around 160angs^2 while the value calculated on PyMOL never surpasses 1.7angs^2.
>>
>> Thank you in advance,
>>
>> Edwin.
>>
>> -- 
>> You received this message because you are subscribed to the Google Groups "cp2k" group.
>> To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
>> To view this discussion visit https://groups.google.com/d/msgid/cp2k/0b7d866a-d3f3-4fda-904a-d8633db2ba2a%40gmail.com.
>> <6WT_ideal.xyz><Clofenotane-COLVAR.metadynLog><FES-Clofenotane.inp>
> 

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/5fe47457-cbc9-40ae-a409-0e04bad8e5e1%40gmail.com.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Ratio_PyMOL-COLVAR.png
Type: image/png
Size: 62100 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20241216/bd3a7a1b/attachment-0001.png>