<div dir="ltr">Dear Juerg Hutter,<div> </div><div>Thank you for taking the time to answer. </div><div>From your replies I gather that I may be using the wrong trajectory file as I am using the ions+centres file and that of course means there are different number of particles in the topology file and the trajectory file.</div><div>This probably should be addressed in the Travis tutorial in the cp2k web page as it appears is not clear for newbies. Maybe seeing this thread some other people will avoid hitting the same stone.</div><div> </div><div>Ion El martes, 5 de diciembre de 2017, 7:32:44 (UTC+1), jgh escribió:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi I am not a TRAVIS user, but looking at your files it seems to me that you have used the wrong 'trajectory' file. CP2K writes a trajectory file during MD that is usually call ....<a href="http://pos-x.xyz" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fpos-x.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEGro4YeFgDBl2j_UA_CeXGm3GTtQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fpos-x.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEGro4YeFgDBl2j_UA_CeXGm3GTtQ';return true;">pos-x.xyz</a>. This file contains the coordinates of all particles. If you calculate Wannier functions during a MD you can also request a file ions+wannier-centers. This last file should only be used for analysis, not as a REFTRAJ! as it contains Wannier center information. If this distinction is not clear from the TRAVIS documentation, tell the TRAVIS developers to upgrade the manual. regards Juerg -------------------------------------------------------------- Juerg Hutter Phone : ++41 44 635 4491 Institut für Chemie C FAX : ++41 44 635 6838 Universität Zürich E-mail: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="eBcjGOr5BgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">hut...@chem.uzh.ch</a> Winterthurerstrasse 190 CH-8057 Zürich, Switzerland --------------------------------------------------------------- -----<a href="javascript:" target="_blank" gdf-obfuscated-mailto="eBcjGOr5BgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a> wrote: ----- To: cp2k <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="eBcjGOr5BgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>> From: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="eBcjGOr5BgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">nju...@gmail.com</a> Sent by: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="eBcjGOr5BgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a> Date: 12/04/2017 04:06PM Subject: [CP2K:9772] Re: RAMAN spectrum from AIMD Dear developers and users, I’m plaguing with the same problem and it’s outputting the same error. Could anyone give some advice? Thank you very much. 在 2017年11月17日星期五 UTC+8下午7:33:13，ION SUCH写道： Dear users and developers, My name is Ion Such and I am a recent CP2K user. I am interested in getting raman spectra of adsorbed molecules. As an exercise I am following one of the HOWTOs on the cp2k web: "How to calculate infrared and Raman spectra from MD with CP2K and TRAVIS". In my case I am simulating a molecule of para-nitroaniline in a box. At the time of writing I have made the first Molecular Dynamics run and have obtained the infrared spectrum and a cp2k input file for the calculation of the polarizabilities by using TRAVIS but I run into an error when trying to execute the polarizabilities calculation that I don't understand and need some help to go ahead. The error in the output file reads MD| Molecular Dynamics Protocol MD| Ensemble Type REFTRAJ MD| Number of Time Steps 2000 MD| Time Step [fs] 0.50 MD| Temperature [K] 300.00 MD| Temperature tolerance [K] 0.00 MD| Print MD information every 1 step(s) MD| File type Print frequency[steps] File names MD| Coordinates 1 <a href="http://polarizability-pos-1.xyz" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fpolarizability-pos-1.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNF6GU56AGZ6B7slUQ1tcNMKn8F9ZQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fpolarizability-pos-1.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNF6GU56AGZ6B7slUQ1tcNMKn8F9ZQ';return true;">polarizability-pos-1.xyz</a> MD| Velocities 1 <a href="http://polarizability-vel-1.xyz" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fpolarizability-vel-1.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNE_BLJGBGOvG4G779F1DyV5FIBlEA';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fpolarizability-vel-1.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNE_BLJGBGOvG4G779F1DyV5FIBlEA';return true;">polarizability-vel-1.xyz</a> MD| Energies 1 polarizability-1.ener MD| Dump 1 polarizability-1.restart ROT| Rotational Analysis Info ROT| Principal axes and moments of inertia in atomic units: ROT| 1 2 3 ROT| EIGENVALUES 0.853999142E+06 0.406672822E+07 0.491400783E+07 ROT| X -0.022761946 -0.821027156 -0.570435187 ROT| Y 0.067435865 0.568022568 -0.820245552 ROT| Z 0.997463933 -0.057138175 0.042437378 ROT| Numer of Rotovibrational vectors: 6 Calculation of degrees of freedom Number of atoms: 16 Number of Intramolecular constraints: 0 Number of Intermolecular constraints: 0 Invariants(translation + rotations): 3 Degrees of freedom: 45 Restraints Information Number of Intramolecular restraints: 0 Number of Intermolecular restraints: 0 ************************** Velocities initialization ************************** Initial Temperature 300.00 K COM velocity: 0.000000000000 -0.000000000000 0.000000000000 ******************************************************************************* ******************************************************************************* * ___ * * / \ * * [ABORT] * * \___/ CPASSERT failed * * | * * O/| * * /| | * * / \ motion/integrator.F:1541 * ******************************************************************************* ===== Routine Calling Stack ===== 3 velocity_verlet 2 qs_mol_dyn_low 1 CP2K The relevant part of the TRAVIS generated input is the following: &MOTION &MD ENSEMBLE REFTRAJ STEPS 2000 TIMESTEP 0.500000 &REFTRAJ EVAL_ENERGY_FORCES TRAJ_FILE_NAME <a href="http://polarizability_reftraj.xyz" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fpolarizability_reftraj.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFcfKtN3lfL_F4cVAioM0hGkKsOlA';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fpolarizability_reftraj.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFcfKtN3lfL_F4cVAioM0hGkKsOlA';return true;">polarizability_reftraj.xyz</a> LAST_SNAPSHOT 2001 STRIDE 1 &END REFTRAJ &END MD &PRINT &RESTART &EACH MD 1 &END EACH &END RESTART &END PRINT &END MOTION THANKS Ion -- 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 <a href="javascript:" target="_blank" gdf-obfuscated-mailto="eBcjGOr5BgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.com</a>. To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="eBcjGOr5BgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>. Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/group/cp2k</a>. For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/optout</a>. </blockquote></div></div>