MD Equilibration Issue – Temperature Fluctuations in CP2K

132 views
Skip to first unread message

Hasan Tunçer

unread,
Feb 3, 2025, 11:24:24 AM2/3/25
to cp2k

Hello,

I am running a molecular dynamics  simulation using CP2K, and I am encountering difficulties in achieving equilibration. The temperature fluctuations remain significant, indicating that my system has not yet reached a stable state.

I have attached my CP2K input file for reference.

I would appreciate any insights or suggestions on how to stabilize the temperature and ensure proper equilibration. Thank you in advance for your time and assistance!

Best regards,

Hasantemperature.png

input.inp

Jürg Hutter

unread,
Feb 4, 2025, 3:46:32 AM2/4/25
to cp...@googlegroups.com
Hi

fluctuations in Temperature will be proportional to 1/SQRT(Natoms) and also
depending on you thermostat settings. I guess you are showing us a temperature
trajectory of a rather small molecule?

regards
JH

________________________________________
From: cp...@googlegroups.com <cp...@googlegroups.com> on behalf of Hasan Tunçer <tuncer...@gmail.com>
Sent: Monday, February 3, 2025 5:06 PM
To: cp2k
Subject: [CP2K:21103] MD Equilibration Issue – Temperature Fluctuations in CP2K
Hasan[temperature.png]

--
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+uns...@googlegroups.com<mailto:cp2k+uns...@googlegroups.com>.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/4938f957-d3a7-4cc7-a3d2-79a6910f0c0en%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/4938f957-d3a7-4cc7-a3d2-79a6910f0c0en%40googlegroups.com?utm_medium=email&utm_source=footer>.
temperature.png

Thomas Kühne

unread,
Feb 4, 2025, 4:55:23 AM2/4/25
to 'Dorothea Golze' via cp2k
Dear Hasan Tuncer,

let me add that by default REGION GLOBAL, i.e. only 1 or 3 degrees of freedom (DOFs) of the total system
are coupled to the thermostat. Please use REGION MASSIVE in which case 3 DOFs per atom are coupled
leading to much fast equipartition.

Best,
Thomas Kühne
> To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+uns...@googlegroups.com.
> To view this discussion visit https://groups.google.com/d/msgid/cp2k/ZR0P278MB07599B5E6D8ACF4650A07F439FF42%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM.
> <temperature.png>

Hasan Tunçer

unread,
Feb 6, 2025, 12:06:44 AM2/6/25
to cp2k
Hello,

Thank you for your response! I appreciate the suggestion and will give it a try.

Regards,
Hasan


4 Şubat 2025 Salı tarihinde saat 12:55:23 UTC+3 itibarıyla Thomas Kühne şunları yazdı:

Hasan Tunçer

unread,
Feb 6, 2025, 12:06:53 AM2/6/25
to cp2k

Thank you for your response. Yes, the data corresponds to a temperature trajectory of a 76-atom system. The thermostat settings are included in the input file attached to my previous email, but I am unsure if they are appropriately configured.

I would appreciate your feedback on whether any adjustments are needed.

Best regards,
Hasan

4 Şubat 2025 Salı tarihinde saat 11:46:32 UTC+3 itibarıyla Jürg Hutter şunları yazdı:

Hasan Tunçer

unread,
Feb 16, 2025, 1:48:05 AM2/16/25
to cp2k
Hello,

I actually did what you suggested me, however still it the temperature is flactuating so much. Do you have any other suggestion?

Thanks,
Hasan

4 Şubat 2025 Salı tarihinde saat 12:55:23 UTC+3 itibarıyla Thomas Kühne şunları yazdı:
Dear Hasan Tuncer,
Reply all
Reply to author
Forward
0 new messages