[CP2K-user] [CP2K:19080] Re: Ehrenfest dynamics, inconsistent results from versions 9.1 and 2023.1

Guillaume Le Breton guillaume_le_breton at live.fr
Fri Jul 7 09:16:36 UTC 2023


Dear Natalia,

I have modified recently the calculation of the electron-field energy for 
real-time field. Previously, this calculation was made on the grid, now it 
is done as an expectation value of the dipole operator. This was done in 
order to compare the total electronic energy with the velocity gauge 
approach (now available for Real-Time Propagation, not yet for Ehrenfest 
MD). 

This modification of the electronic energy calculation should have an 
effect only when the electric field is on. If you continue this calculation 
after the pulse end you should recover the same energy as for the previous 
version of CP2K. 
Moreover, the other quantities such as the dipole moment or the nuclei 
displacement should be the same for this version and the 9.1 throughout the 
whole simulation. 

Please let us know if this is not the case.

Best,
Guillaume 


Le jeudi 6 juillet 2023 à 12:45:27 UTC+2, Natalia K a écrit :

> Dear community,
>
> I am performing test calculations for a system of 2 silver atoms under an 
> external field. Using the same input files (attached), I get completely 
> different energies from the two versions of the code (see images attached). 
> The results from cp2k-9.1 look correct to me. Is this the code problem or 
> should I modify the input file? It seems consistent when compared to the 
> files from the /tests/regtests.. files. 
>
> P.S. the xyz file is obtained from the geometry optimization with the 
> version 9.1, but I wouldn't think this is the problem.
>

-- 
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/5add1d9b-e6b9-45a1-87f7-5ceddb254869n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230707/dc00be28/attachment.htm>


More information about the CP2K-user mailing list