Dear Leili,<div><br /></div><div>Sorry for the late reply. Two files for you version 2022.1 are attached. You have to make substitution into src/motion and recompile the code: </div><div>make ARCH=local VERSION="psmp".</div><div><br /></div><div>To specify the range of atoms to shift from equilibrium positions, just add the keywords START_ATOM_ID and END_ATOM_ID:</div><div><br /></div><div>&VIBRATIONAL_ANALYSIS<br />DX 0.005<br /><span style="background-color: white;"><font color="#008000">START_ATOM_ID 1<br />END_ATOM_ID 100</font></span><br />.....................</div><div>.....................</div><div><br /></div><div>For instance, you can split the task of 600 atoms into 6 parts 1-100,101-200, etc and then gather the force derivatives into the dynamical matrix.</div><div>If you need further assistance, please let me know.</div><div><br /></div><div>Best regards,</div><div>Mikhail</div><div><br /><br /></div><div><div dir="auto">On Wednesday, January 3, 2024 at 4:47:41 AM UTC+1 Leili Rassouli wrote:<br /></div><blockquote style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Thank you, Eric. I made a new inquiry.<br /><br /><br /><div><div dir="auto">On Tuesday, January 2, 2024 at 1:50:00 PM UTC-5 Eric Patterson wrote:<br /></div><blockquote style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Dear Lili,<div><br /></div><div>I suggest you make your TDDFT inquiry a new thread since it is a separate issue. I have no experience with TDDFT in CP2K, so I’m afraid I’m of no use in this case..</div><div><br /></div><div>Best,</div><div>Eric</div><div><br /><div><br /><blockquote type="cite"></blockquote></div></div></div><div><div><div><blockquote type="cite"><div>On Jan 2, 2024, at 1:34 PM, Leili Rassouli <<a rel="nofollow">rassoul...@gmail.com</a>> wrote:</div><br /></blockquote></div></div></div><div><div><div><blockquote type="cite"><div>Dear Mikhail and Eric,<br /><br />Thank you for your responses. My job is completed after 10 days, but to enhance the computational efficiency, I am considering the use of a mode-selective method or the dividing the full range into several parts to potentially expedite the calculation process. Mikhail, I would be appreciative if you could provide me access to the source files, given that I am currently utilizing CP2K v2022.1.<br /><br />I'm currently facing an issue with negative frequencies. While using the ground state wavefunction in the ground state geometry, frequencies are positive and within the expected range. However, switching to the first excited wavefunction yields a significant number of negative frequencies.<br /><br />Enclosed is my input file for 221 supercells of Fe2O3 with 120 atoms. Your insights and recommendations on this matter would be greatly valued.<br /><br />I've explored various parameters, including:<br /><br />Adjusting the CUTOFF value of TDDFPT from 300 to 600 and 700<br />Refining the optimized geometry<br />Increasing the TDDFT convergence from 1.0e-5 to 1.0e-6 and 1.0e-7<br />Despite these efforts, a notable proportion of frequencies persist as negative values. <br /><br />I express my gratitude in advance for your time and assistance.<br /><br />Best regards,<br />Lili<br /><br /><br /><br /><br /><br /><br /><br /><br /><br /><div><div dir="auto">On Tuesday, January 2, 2024 at 11:21:34 AM UTC-5 Eric Patterson wrote:<br /></div><blockquote style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Hello Lili,<div><br /></div><div>I did do a little bit of testing of the mode-selective method vs the normal method. The results are comparable if I choose all atoms and restart until all modes are converged (or at least nearly so — I didn’t care vey much about low energy modes).</div><div><br /></div><div>However, I’ve also managed to compile a working MPI version on my HPC, so I am able to run on enough nodes to get through a default vibrational analysis before the queue kills the job. </div><div><br /></div><div>Best,</div><div>Eric</div><div><br /><div><br /><blockquote type="cite"></blockquote></div></div></div><div><div><div><blockquote type="cite"><div>On Dec 31, 2023, at 3:24 PM, Leili Rassouli <<a rel="nofollow">rassoul...@gmail.com</a>> wrote:</div><br /></blockquote></div></div></div><div><div><div><blockquote type="cite"><div>Hi Eric,<br /><br />I've faced a similar issue and was considering the approach you suggested. Have you had a chance to try it, and does it work well in your case? Your insights would be greatly appreciated.<br /><br />Thank you for your assistance.<br />Lili Rassouli<br /><br /><div><div dir="auto">On Friday, February 10, 2023 at 11:32:16 AM UTC-5 Eric Patterson wrote:<br /></div><blockquote style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>On a closely-related note…<div><br /></div><div>Since normal vibrational analysis is not restartable, but mode-selective vibrational analysis is… Is there any drawback to using mode-selective and choosing all atoms, and then restarting as necessary until all modes have converged? Given the queue time limit imposed by my HPC resources, this seems to me the only way to get the job done. Are there any pitfalls one should bear in mind?</div><div><br /></div><div>Thanks,</div><div>Eric</div><div><br /><div><br /><blockquote type="cite"></blockquote></div></div></div><div><div><div><blockquote type="cite"><div>On Feb 9, 2023, at 6:03 AM, Krack Matthias <<a rel="nofollow">matthia...@psi.ch</a>> wrote:</div><br /></blockquote></div></div></div><div><div><div><blockquote type="cite"><div><div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><div style="margin: 0cm; font-size: 10pt; font-family: Calibri, sans-serif;"><span lang="EN-US" style="font-size: 11pt;">It is not possible to restart a vibrational analysis run. Such a functionality is not implemented.<u></u><u></u></span></div><div style="margin: 0cm; font-size: 10pt; font-family: Calibri, sans-serif;"><span lang="EN-US" style="font-size: 11pt;"><u></u> <u></u></span></div><div style="border-style: solid none none; border-top-width: 1pt; border-top-color: rgb(181, 196, 223); padding: 3pt 0cm 0cm;"><p style="margin: 0cm 0cm 12pt 36pt; font-size: 10pt; font-family: Calibri, sans-serif;"><b><span style="font-size: 12pt;">From:<span> </span></span></b><span style="font-size: 12pt;"><a style="color: blue; text-decoration: underline;" rel="nofollow">cp...@googlegroups.com</a><span> </span><<a style="color: blue; text-decoration: underline;" rel="nofollow">cp...@googlegroups.com</a>> on behalf of Zac Smith <<a style="color: blue; text-decoration: underline;" rel="nofollow">qingxi...@gmail.com</a>><br /><b>Date:<span> </span></b>Thursday, 9 February 2023 at 02:54<br /><b>To:<span> </span></b>cp2k <<a style="color: blue; text-decoration: underline;" rel="nofollow">cp...@googlegroups.com</a>><br /><b>Subject:<span> </span></b>[CP2K:18425] Help: cp2k vibration analysis stopped in the middle, how to renew the calculation?<u></u><u></u></span></p></div><div><div style="margin: 0cm 0cm 0cm 36pt; font-size: 10pt; font-family: Calibri, sans-serif;"><span style="font-size: 11pt;">Hello cp2k users, <u></u><u></u></span></div></div><div style="margin: 0cm 0cm 0cm 36pt; font-size: 10pt; font-family: Calibri, sans-serif;"><span style="font-size: 11pt;">I used a virtual machine to do vibration analysis with cp2k, but the task was killed when the hard disk was full, and the output file only contains out, molden, cube, and the eig file of size 0. Can the calculation be continued in this case?<u></u><u></u></span></div><div style="margin: 0cm 0cm 0cm 36pt; font-size: 10pt; font-family: Calibri, sans-serif;"><span style="font-size: 11pt;">--<span> </span><br />You received this message because you are subscribed to the Google Groups "cp2k" group.<br />To unsubscribe from this group and stop receiving emails from it, send an email to<a style="color: blue; text-decoration: underline;" rel="nofollow">cp2k+uns...@googlegroups.com</a>.<br />To view this discussion on the web visit<span> </span><a href="https://groups.google.com/d/msgid/cp2k/9d914e4f-bcb1-42ae-9b02-6c14f6fec8a5n%40googlegroups.com?utm_medium=email&utm_source=footer" style="color: blue; text-decoration: underline;" rel="nofollow" target="_blank">https://groups.google.com/d/msgid/cp2k/9d914e4f-bcb1-42ae-9b02-6c14f6fec8a5n%40googlegroups.com</a>.<u></u><u></u></span></div></div><div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;"><br /></div><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; float: none; display: inline;">--<span> </span></span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;" /><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; float: none; display: inline;">You received this message because you are subscribed to the Google Groups "cp2k" group.</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;" /><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; float: none; display: inline;">To unsubscribe from this group and stop receiving emails from it, send an email to<span> </span></span><a style="color: blue; text-decoration: underline; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;" rel="nofollow">cp2k+uns...@googlegroups.com</a><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; float: none; display: inline;">.</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;" /></div></blockquote></div></div></div><div><div><div><blockquote type="cite"><div><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; float: none; display: inline;">To view this discussion on the web visit<span> </span></span><a href="https://groups.google.com/d/msgid/cp2k/ZRAP278MB082784950ADD367A217A8C43F4D99%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer" style="color: blue; text-decoration: underline; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;" rel="nofollow" target="_blank">https://groups.google.com/d/msgid/cp2k/ZRAP278MB082784950ADD367A217A8C43F4D99%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; float: none; display: inline;">.</span></div></blockquote></div><br /></div></div></blockquote></div><div><br /></div>
-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp2k+uns...@googlegroups.com</a>.<br /></div></blockquote></div></div></div><div><div><div><blockquote type="cite"><div>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0d3984ac-a0a6-4040-b500-6a5d67f5c556n%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank">https://groups.google.com/d/msgid/cp2k/0d3984ac-a0a6-4040-b500-6a5d67f5c556n%40googlegroups.com</a>.<br />
</div></blockquote></div><br /></div></div></blockquote></div><div><br /></div>
-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp2k+uns...@googlegroups.com</a>.<br /></div></blockquote></div></div></div><div><div><div><blockquote type="cite"><div>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/3e2997be-cac2-48d0-aeca-78080fa420f9n%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank">https://groups.google.com/d/msgid/cp2k/3e2997be-cac2-48d0-aeca-78080fa420f9n%40googlegroups.com</a>.<br />
<span><Fe2O3_vib.inp></span></div></blockquote></div><br /></div></div></blockquote></div></blockquote></div>
<p></p>
-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ced8a300-d1f5-4d76-8c66-7cd3cfed6bd2n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ced8a300-d1f5-4d76-8c66-7cd3cfed6bd2n%40googlegroups.com</a>.<br />