[CP2K-user] [CP2K:16318] Error reading coordinates from restart file in QM/MM force mixing calculations

[Pedro Maximiano]

Dear CP2K users and developers

I'm trying to do QM/MM simulations of a condensation reaction between two 
orthosilicic acid molecules, using the force mixing method to model the 
first solvation layer at the QM level. 

I have noticed that, when starting the QM/MM simulations from a restart 
file of a previous fully-MM run or when restarting the QM/MM simulation 
itself, the coordinates of the restart file are not read (even though 
RESTART_POS .TRUE. is activated in the &EXT_RESTART section), yielding the 
following error message:
 *******************************************************************************
 *   
___                                                                        
                                                     *
 *  /   
\                                                                           
                                                  *
 * 
[ABORT]                                                                    
                                                  *
 *  \___/                 No atomic coordinates have been 
found!                                 *
 *    
|                                                                                                                                
*
 *  
O/|                                                                         
                                                     *
 * /| 
|                                                                                                                                
*
 * / \                                                          
topology.F:595                                            *
 *******************************************************************************

This issue happens only if force mixing is enabled. However, if I copy the 
coordinates from the &COORD section of the restart file into an external 
file and then load that file in the section &FORCE_EVAL / &SUBSYS / 
&TOPOLOGY the simulation runs normally.

On the other hand, if I load the same coordinates as a pdb file and not a 
coordinate file in CP2K format, the simulation runs but the atoms appear 
almost frozen in the trajectory. It is not just a visualization artifact, 
since the fluctuations in temperature and energies change completely and 
the RDF of water oxygen atoms around the Si atoms is completely wrong. This 
is true even when the QM/MM simulation is restarted in this way.

The calculations were done with the xTB model for the QM region, although 
the same issues occurred with PM6 and PM6-FM. The coordinates issue also 
occurred when starting a DFT GPW simulation from a restart file. I have 
tried both the 7.1 and 8.2 versions of the code.
Are these issues related to a bug or am I missing something?

Best regards
Pedro Maximiano

-- 
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 on the web visit https://groups.google.com/d/msgid/cp2k/6aef7662-7c08-40e0-9bd8-cfb746484c74n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20211129/bcac6c56/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MM-MD-1-last.pdb
Type: application/x-pdb
Size: 474582 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20211129/bcac6c56/attachment.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MM-MD-1.restart
Type: application/octet-stream
Size: 789591 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20211129/bcac6c56/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MD-1.coord
Type: application/octet-stream
Size: 380564 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20211129/bcac6c56/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: interface.prm
Type: application/octet-stream
Size: 201402 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20211129/bcac6c56/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: box.psf
Type: application/x-font-linux-psf
Size: 665372 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20211129/bcac6c56/attachment-0001.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: md.inp
Type: chemical/x-gamess-input
Size: 3807 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20211129/bcac6c56/attachment.inp>