Hi Prof. Hutter,<div><br /></div><div>Thanks for solving this issue, it is very much appreciated!</div><div><br /></div><div>Best,</div><div>Léon<br /><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, 30 June 2023 at 13:33:21 UTC+2 Jürg Hutter wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi
<br>
<br>I generalized the CP2K CIF parser, so that _atom_site_label is no longer required to be the first of
<br>the _atom_site descriptors.
<br>
<br>Thanks for reporting.
<br>
<br>JH
<br>
<br>________________________________________
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a href data-email-masked rel="nofollow">Leon.luntad...@uantwerpen.be</a>>
<br>Sent: Thursday, June 29, 2023 11:41 AM
<br>To: cp2k
<br>Subject: Re: [CP2K:19030] Re: TiO2 slab not converging
<br>
<br>So I dug a bit deeper into why the CIF from ASE could not be properly read by CP2K even though it looks correct and apparently something goes wrong when CP2K reads in the CIF file. I noticed this because the same problem happened when I took a molecule in a 25x25x25 Å box from a CIF (created by Critic2, see attachment 3mppa.cif) which could be read by CP2K and used ASE to alter the vacuum size and center the molecule and suddenly the resulting CIF could no longer be properly read by CP2K (i.e., I got errors stating atoms are too close). I noticed that even when just reading in and writing out the CIF with ASE (see attachment 3mppa-ase.cif) already showed atoms too close to one another when read into CP2K:
<br>
<br>   Atom Kind Element         X             Y             Z       Z(eff)     Mass
<br>      1    1 H     1      7.070300      9.682957      5.879925   1.0000   1.0079
<br>      2    2 O     8      7.070300      9.423014      6.500614   6.0000  15.9994
<br>      3    1 H     1      7.070300      5.945227      5.744394   1.0000   1.0079
<br>      4    2 O     8      7.070300      6.657800      6.398468   6.0000  15.9994
<br>      5    3 P    15      7.070300      8.447808      5.896924   5.0000  30.9738
<br>      6    2 O     8      7.070300      8.566346      4.554013   6.0000  15.9994
<br>      7    1 H     1      7.070300      9.631506      8.443233   1.0000   1.0079
<br>      8    1 H     1      7.070300      7.722120      8.375391   1.0000   1.0079
<br>      9    4 C     6      7.070300      8.978791      8.136181   4.0000  12.0107
<br>     10    4 C     6      7.070300      9.143231      6.759885   4.0000  12.0107
<br>     11    1 H     1      7.070300     10.392672      6.492497   1.0000   1.0079
<br>     12    5 S    16      7.070300      9.397121     10.419543   6.0000  32.0650
<br>     13    1 H     1      7.070300      8.473533      6.425727   1.0000   1.0079
<br>     14    4 C     6      7.070300      9.608639      8.784895   4.0000  12.0107
<br>     15    1 H     1      7.070300     10.868673      8.560234   1.0000   1.0079
<br>     16    1 H     1      7.070300      8.934066      8.505103   1.0000   1.0079
<br>     17    1 H     1      7.070300      9.995695     10.760287   1.0000   1.0079
<br>
<br>
<br> *** WARNING in particle_methods.F:677 :: The distance between the atoms 2 ***
<br> *** and 10 is only 0.381 angstrom and thus smaller than the threshold of  ***
<br> *** 0.500 angstrom                                                        ***
<br>
<br>
<br> *** WARNING in particle_methods.F:677 :: The distance between the atoms 7 ***
<br> *** and 14 is only 0.342 angstrom and thus smaller than the threshold of  ***
<br> *** 0.500 angstrom                                                        ***
<br>
<br>
<br> *** WARNING in particle_methods.F:677 :: The distance between the atoms 9 ***
<br> *** and 16 is only 0.372 angstrom and thus smaller than the threshold of  ***
<br> *** 0.500 angstrom                                                        ***
<br>
<br>Notice that all the x-coordinates are the same, even though the fractional x-coordinates in the CIF are different for all atoms... I noticed that by switching the order of the lines "_atom_site_type_symbol" and "_atom_site_label" in the CIF (see 3mppa-ase-edited.cif) CP2K could correctly read in the CIF file without a hitch.
<br>
<br>   Atom Kind Element         X             Y             Z       Z(eff)     Mass
<br>      1    1 H     1      8.143430      6.455697      4.999966   1.0000   1.0079
<br>      2    2 O     8      7.924816      7.137164      5.661084   6.0000  15.9994
<br>      3    3 H     1      4.999975      6.306894      6.438132   1.0000   1.0079
<br>      4    4 O     8      5.599253      7.025016      6.709042   6.0000  15.9994
<br>      5    5 P    15      7.104662      6.474360      6.880102   5.0000  30.9738
<br>      6    6 O     8      7.204353      4.999949      6.935846   6.0000  15.9994
<br>      7    7 H     1      8.100160      9.270007      7.237998   1.0000   1.0079
<br>      8    8 H     1      6.494353      9.195522      7.971859   1.0000   1.0079
<br>      9    9 C     6      7.551222      8.932889      8.129596   4.0000  12.0107
<br>     10   10 C     6      7.689517      7.421824      8.281361   4.0000  12.0107
<br>     11   11 H     1      8.740305      7.128253      8.426234   1.0000   1.0079
<br>     12   12 S    16      7.903040     11.439841      9.116297   6.0000  32.0650
<br>     13   13 H     1      7.126297      7.054944      9.153051   1.0000   1.0079
<br>     14   14 C     6      8.080929      9.645126      9.361478   4.0000  12.0107
<br>     15   15 H     1      9.140625      9.398466      9.518603   1.0000   1.0079
<br>     16   16 H     1      7.513608      9.337936     10.250780   1.0000   1.0079
<br>     17   17 H     1      8.406445     11.813951     10.314334   1.0000   1.0079
<br>
<br>Perhaps someone with more knowledge about the process of reading in a CIF file in CP2K can take a look at the code and where/why it does not work if the two labels are in the original order and maybe fix this in an upcoming version?
<br>
<br>Best,
<br>Léon
<br>On Wednesday, 28 June 2023 at 16:52:05 UTC+2 Léon Luntadila Lufungula wrote:
<br>Hi Matthias,
<br>
<br>Thanks a lot, now it does not throw the error anymore. I'll keep you updated if this fixes the convergence issues but I suspect it will. Again, thanks for the help!
<br>
<br>All the best,
<br>
<br>Léon
<br>
<br>On Wednesday, 28 June 2023 at 16:37:28 UTC+2 Krack Matthias wrote:
<br>Hi Léon
<br>
<br>Try
<br>
<br>&CELL
<br> A   10.2988004684         0.0000000000         0.0000000000
<br> B    0.0000000000        <a href="tel:(354)%20999-5422" value="+13549995422" target="_blank" rel="nofollow">11.3549995422</a><tel:(354)%20999-5422>         0.0000000000
<br> C    0.0000000000         0.0000000000        49.3899002075
<br>&END CELL
<br>
<br>&COORD
<br>O    9.300228691         0.000000000        19.999946640
<br> O    9.300228691         3.784962059        19.999946640
<br> O    9.300228691         7.570037822        19.999946640
<br>. . . further 69 O atoms
<br>36 Ti atoms
<br>&END COORD
<br>
<br>HTH
<br>
<br>Matthias
<br>
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a href data-email-masked rel="nofollow">Leon.luntad...@uantwerpen.be</a>>
<br>Date: Wednesday, 28 June 2023 at 16:29
<br>To: cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>>
<br>Subject: Re: [CP2K:19022] Re: TiO2 slab not converging
<br>Hi Matthias,
<br>
<br>So I opened the CIF file in VESTA and did Export Data > VASP file > Save > Cartesian coordinates > Output current unit cell and this resulted in the file shown below, but I don't see how this will "easily" cut&paste into a CP2K input file as the atom labels are missing...
<br>
<br>"O72 Ti36"
<br>1.0
<br>       10.2988004684         0.0000000000         0.0000000000
<br>        0.0000000000        <a href="tel:(354)%20999-5422" value="+13549995422" target="_blank" rel="nofollow">11.3549995422</a><tel:(354)%20999-5422>         0.0000000000
<br>        0.0000000000         0.0000000000        49.3899002075
<br>    O   Ti
<br>   72   36
<br>Cartesian
<br>     9.300228691         0.000000000        19.999946640
<br>     9.300228691         3.784962059        19.999946640
<br>     9.300228691         7.570037822        19.999946640
<br>     4.151034406         1.892992021        20.000439738
<br>     4.151034406         5.677953911        20.000439738
<br>     4.151034406         9.463029336        20.000439738
<br>     7.074658142         1.892992021        20.879580124
<br>     7.074658142         5.677953911        20.879580124
<br>     7.074658142         9.463029336        20.879580124
<br>     1.923610010         0.000000000        20.880567792
<br>     1.923610010         3.784962059        20.880567792
<br>     1.923610010         7.570037822        20.880567792
<br>     5.583803656         0.000000000        21.468802556
<br>     5.583803656         3.784962059        21.468802556
<br>     5.583803656         7.570037822        21.468802556
<br>     0.433682486         1.892992021        21.469295653
<br>     0.433682486         5.677953911        21.469295653
<br>     0.433682486         9.463029336        21.469295653
<br>     3.357409004         1.892992021        22.348436040
<br>     3.357409004         5.677953911        22.348436040
<br>     3.357409004         9.463029336        22.348436040
<br>     8.505058213         0.000000000        22.349423707
<br>     8.505058213         3.784962059        22.349423707
<br>     8.505058213         7.570037822        22.349423707
<br>     0.392487277         0.000000000        23.519964286
<br>     0.392487277         3.784962059        23.519964286
<br>     0.392487277         7.570037822        23.519964286
<br>     5.542093345         1.892992021        23.520458855
<br>     5.542093345         5.677953911        23.520458855
<br>     5.542093345         9.463029336        23.520458855
<br>     8.465716774         1.892992021        24.399597769
<br>     8.465716774         5.677953911        24.399597769
<br>     8.465716774         9.463029336        24.399597769
<br>     3.314668949         0.000000000        24.400586909
<br>     3.314668949         3.784962059        24.400586909
<br>     3.314668949         7.570037822        24.400586909
<br>     6.974965416         0.000000000        24.988818729
<br>     6.974965416         3.784962059        24.988818729
<br>     6.974965416         7.570037822        24.988818729
<br>     1.824741540         1.892992021        24.989807868
<br>     1.824741540         5.677953911        24.989807868
<br>     1.824741540         9.463029336        24.989807868
<br>     4.748467944         1.892992021        25.868455157
<br>     4.748467944         5.677953911        25.868455157
<br>     4.748467944         9.463029336        25.868455157
<br>     9.896220587         0.000000000        25.869935922
<br>     9.896220587         3.784962059        25.869935922
<br>     9.896220587         7.570037822        25.869935922
<br>     1.783546255         0.000000000        27.040476500
<br>     1.783546255         3.784962059        27.040476500
<br>     1.783546255         7.570037822        27.040476500
<br>     6.933152591         1.892992021        27.040476500
<br>     6.933152591         5.677953911        27.040476500
<br>     6.933152591         9.463029336        27.040476500
<br>     9.856776020         1.892992021        27.919615415
<br>     9.856776020         5.677953911        27.919615415
<br>     9.856776020         9.463029336        27.919615415
<br>     4.705727889         0.000000000        27.921099124
<br>     4.705727889         3.784962059        27.921099124
<br>     4.705727889         7.570037822        27.921099124
<br>     8.366024662         0.000000000        28.508839318
<br>     8.366024662         3.784962059        28.508839318
<br>     8.366024662         7.570037822        28.508839318
<br>     3.215800326         1.892992021        28.509825513
<br>     3.215800326         5.677953911        28.509825513
<br>     3.215800326         9.463029336        28.509825513
<br>     6.139527190         1.892992021        29.388472802
<br>     6.139527190         5.677953911        29.388472802
<br>     6.139527190         9.463029336        29.388472802
<br>     0.988478905         0.000000000        29.389953567
<br>     0.988478905         3.784962059        29.389953567
<br>     0.988478905         7.570037822        29.389953567
<br>     7.442016174         0.000000000        20.734373862
<br>     7.442016174         3.784962059        20.734373862
<br>     7.442016174         7.570037822        20.734373862
<br>     2.292821888         1.892992021        20.734868432
<br>     2.292821888         5.677953911        20.734868432
<br>     2.292821888         9.463029336        20.734868432
<br>     5.215518701         1.892992021        21.614008818
<br>     5.215518701         5.677953911        21.614008818
<br>     5.215518701         9.463029336        21.614008818
<br>     0.065397382         0.000000000        21.614996485
<br>     0.065397382         3.784962059        21.614996485
<br>     0.065397382         7.570037822        21.614996485
<br>     8.833075420         0.000000000        24.254391507
<br>     8.833075420         3.784962059        24.254391507
<br>     8.833075420         7.570037822        24.254391507
<br>     3.683880827         1.892992021        24.254886077
<br>     3.683880827         5.677953911        24.254886077
<br>     3.683880827         9.463029336        24.254886077
<br>     6.606680461         1.892992021        25.134024991
<br>     6.606680461         5.677953911        25.134024991
<br>     6.606680461         9.463029336        25.134024991
<br>     1.456456432         0.000000000        25.135014131
<br>     1.456456432         3.784962059        25.135014131
<br>     1.456456432         7.570037822        25.135014131
<br>    10.224134052         0.000000000        27.774409152
<br>    10.224134052         3.784962059        27.774409152
<br>    10.224134052         7.570037822        27.774409152
<br>     5.074939767         1.892992021        27.774903722
<br>     5.074939767         5.677953911        27.774903722
<br>     5.074939767         9.463029336        27.774903722
<br>     7.997739707         1.892992021        28.654537206
<br>     7.997739707         5.677953911        28.654537206
<br>     7.997739707         9.463029336        28.654537206
<br>     2.847515371         0.000000000        28.655031776
<br>     2.847515371         3.784962059        28.655031776
<br>     2.847515371         7.570037822        28.655031776
<br>
<br>On Wednesday, 28 June 2023 at 16:15:00 UTC+2 Krack Matthias wrote:
<br>Hi Léon
<br>
<br>You can export cell and coordinates as a POSCAR file in VESTA and from that file you can easily cut&paste CP2K’s &CELL and &COORD section.
<br>
<br>HTH
<br>
<br>Matthias
<br>
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a href data-email-masked rel="nofollow">Leon.luntad...@uantwerpen.be</a>>
<br>Date: Wednesday, 28 June 2023 at 16:10
<br>To: cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>>
<br>Subject: Re: [CP2K:19020] Re: TiO2 slab not converging
<br>Hi Matthias,
<br>
<br>Thanks a lot for the tip! The output when using the pre-compiled binary shows that some atoms are indeed too close together, which is probably the cause of the bad convergence. Is there an easy way to fix this?
<br>
<br>I made the slabs with ASE so I'm quite surprised that it did not make the slabs correctly. Here is the script I used to make them with ASE, perhaps you can see where I made a mistake, resulting in the faulty CIF files...
<br>
<br>from <a href="http://ase.io" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=http://ase.io&source=gmail&ust=1688458718677000&usg=AOvVaw356CxpjQ3V1LUv27TnERYV">ase.io</a><<a href="http://ase.io" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=http://ase.io&source=gmail&ust=1688458718677000&usg=AOvVaw356CxpjQ3V1LUv27TnERYV">http://ase.io</a>> import read, write
<br>from ase.visualize import view
<br>from ase.build.surfaces_with_termination import surfaces_with_termination
<br>from ase.build.tools import sort
<br>
<br>ana = read('anatase_opt_tzvp.pdb')
<br>
<br>for N in range(2,13):
<br>    s = surfaces_with_termination(ana, (1,0,1), N, vacuum=20, termination='O')
<br>    s = s[1]
<br>    s = s.repeat((1,3,1))
<br>    s = sort(s, tags=s.positions[:, 2])
<br>    name = f'a101-{N}-layer'
<br>    write(f'{name}.cif', s)
<br>
<br>Léon
<br>On Wednesday, 28 June 2023 at 15:58:25 UTC+2 Krack Matthias wrote:
<br>Hi Léon
<br>
<br>You are using CP2K version 8.2. Newer CP2K versions perform a distance check for CP2K/Quickstep by default (see also this input keyword<<a href="https://manual.cp2k.org/cp2k-2023_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/PRINT/INTERATOMIC_DISTANCES.html#CHECK_INTERATOMIC_DISTANCES" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://manual.cp2k.org/cp2k-2023_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/PRINT/INTERATOMIC_DISTANCES.html%23CHECK_INTERATOMIC_DISTANCES&source=gmail&ust=1688458718677000&usg=AOvVaw3qdPjeVc_J5vhl0k_bP617">https://manual.cp2k.org/cp2k-2023_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/PRINT/INTERATOMIC_DISTANCES.html#CHECK_INTERATOMIC_DISTANCES</a>>) when the system is smaller than 2000 atoms and prints a warning if a small distance is found.
<br>You can download a pre-compiled binary of CP2K 2023.1 here<<a href="https://github.com/cp2k/cp2k/releases/download/v2023.1/cp2k-2023.1-Linux-gnu-x86_64.ssmp" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://github.com/cp2k/cp2k/releases/download/v2023.1/cp2k-2023.1-Linux-gnu-x86_64.ssmp&source=gmail&ust=1688458718677000&usg=AOvVaw1p9JC9NpVO2ckNN7EbSg4o">https://github.com/cp2k/cp2k/releases/download/v2023.1/cp2k-2023.1-Linux-gnu-x86_64.ssmp</a>>, which is sufficient for that initial check.
<br>
<br>HTH
<br>
<br>Matthias
<br>
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a href data-email-masked rel="nofollow">Leon.luntad...@uantwerpen.be</a>>
<br>Date: Wednesday, 28 June 2023 at 15:46
<br>To: cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>>
<br>Subject: Re: [CP2K:19017] Re: TiO2 slab not converging
<br>
<br>Hi Matthias,
<br>
<br>
<br>
<br>Thanks for you reply!
<br>
<br>
<br>
<br>I looked at my CIF file in various visualization software and it looked okay there. Is there a way to check this via CP2K to be sure?
<br>
<br>I also have some preliminary results of some ongoing calculations with DIIS/FULL_ALL and DIIS/FULL_SINGLE_INVERSE and these seem to be a bit faster than CG but reaching convergence is still a big issue (see attachments).
<br>
<br>Best,
<br>Léon
<br>
<br>On Wednesday, 28 June 2023 at 15:35:32 UTC+2 Krack Matthias wrote:
<br>Hi Léon
<br>
<br>Did you check the set of atomic coordinates read from the cif file for close contacts possibly caused by the periodic boundary conditions?
<br>
<br>Best
<br>
<br>Matthias
<br>
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a href data-email-masked rel="nofollow">Leon.luntad...@uantwerpen.be</a>>
<br>Date: Wednesday, 28 June 2023 at 14:16
<br>To: cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>>
<br>Subject: [CP2K:19014] Re: TiO2 slab not converging
<br>I forgot to add the cif file which is read with the input, so here it is if someone needs it.
<br>On Wednesday, 28 June 2023 at 13:54:37 UTC+2 Léon Luntadila Lufungula wrote:
<br>Update:
<br>
<br>I tried to optimize a 3-layer slab with altered SCF cycle settings (inner loop: 15 steps, outer loop: 45 steps) and still it does not converge. The convergence goes towards 2E-7 after 22 outer SCF cycles, but then doesn't seem to get lower (see attachment a101-3-layer-jobID_188.out). I tried to restart the calculation with altered EPS_SCF (= 1E-6) and then the restarted SCF cycle converges as the convergence was already below 1E-6. Unfortunately, after that the SCF again seems to only very slowly converge, after 27 outer SCF steps it's still only at 8E-5 convergence (see attachment a101-3-layer-jobID_189.out). Can anybody please help me out? I feel like the calculation is not converging properly, as after nearly 2 hours it still hasn't converged one SCF cycle but I don't know what might be wrong orhow to fix it... Requiring 300+ steps to reach SCF convergence seems quite high to me and makes me start to think there is something fundamentally wrong with my input. If anyone could look at my inputs and tell me what's wrong with them, I would be forever thankful!
<br>
<br>Kind regards,
<br>Léon
<br>
<br>P.S. Could the TZVP basis set pose a problem? I often see people use DZVP-SR for similar calculations, but my promotor and co-promotor both told me that TZVP is required for proper energetics and should be used in production runs.
<br>On Wednesday, 28 June 2023 at 10:04:07 UTC+2 Léon Luntadila Lufungula wrote:
<br>Dear all,
<br>
<br>I am trying to determine the ideal slab thickness for an anatase (101) slab but my SCF cycles are not converging. I optimized the cutoffs for the anatase bulk and did a CELL_OPT calculation on the unit cell to get the converged structure. From this optimized unit cell, I then used ASE to build the (101) slabs  with a 20 Å vacuum width and thicknesses varying from 2 to 12 layers. During the GEO_OPT calculation I keep the bottom layer frozen to act as the bulk. I started the 2-layer calculation yesterday but when I looked at the output today, I see that the SCF cycle fails at every step of the geometry optimization (see output in attachment). Is there something wrong with my setup (too tight convergence criteria, wrong SCF setup, ...), or is it because a 2-layer structure is perhaps too small? The changes in energy go to 1E-5 - 1E-6 but never reach the convergence criteria of 1E-7 even after 300 steps (inner cycle: 10 steps, outer cycle 30 steps). I would guess that it could converge if I would increase the maximum number of steps, but 300 steps already seems quite a lot to me, although I could be wrong...
<br>
<br>If anyone could share their thoughts on the situation so I can continue my work, I would greatly appreciate the help! Thanks in advance!
<br>
<br>Best regards,
<br>Léon
<br>--
<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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/12b3aec6-38d8-404f-b68f-aa681ccd6266n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/12b3aec6-38d8-404f-b68f-aa681ccd6266n%2540googlegroups.com&source=gmail&ust=1688458718678000&usg=AOvVaw3ZYrYyM8IW6MMV4q77cXuv">https://groups.google.com/d/msgid/cp2k/12b3aec6-38d8-404f-b68f-aa681ccd6266n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/12b3aec6-38d8-404f-b68f-aa681ccd6266n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/12b3aec6-38d8-404f-b68f-aa681ccd6266n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1688458718678000&usg=AOvVaw2ETDWkV_DIWUS8oI9BfBrE">https://groups.google.com/d/msgid/cp2k/12b3aec6-38d8-404f-b68f-aa681ccd6266n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br>--
<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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/94c9ad9e-58b3-4af1-844e-34b82c2fbed6n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/94c9ad9e-58b3-4af1-844e-34b82c2fbed6n%2540googlegroups.com&source=gmail&ust=1688458718678000&usg=AOvVaw16ExUVlDKVXvz5PQ-Qr9SF">https://groups.google.com/d/msgid/cp2k/94c9ad9e-58b3-4af1-844e-34b82c2fbed6n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/94c9ad9e-58b3-4af1-844e-34b82c2fbed6n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/94c9ad9e-58b3-4af1-844e-34b82c2fbed6n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1688458718678000&usg=AOvVaw3FYwAO30e3Qgyhv4d8igR5">https://groups.google.com/d/msgid/cp2k/94c9ad9e-58b3-4af1-844e-34b82c2fbed6n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br>--
<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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/2fb6de53-e6de-44ce-84d4-34ebbe1fb14bn%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/2fb6de53-e6de-44ce-84d4-34ebbe1fb14bn%2540googlegroups.com&source=gmail&ust=1688458718678000&usg=AOvVaw0_u5trdkPO7hC_4AYFTVXa">https://groups.google.com/d/msgid/cp2k/2fb6de53-e6de-44ce-84d4-34ebbe1fb14bn%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/2fb6de53-e6de-44ce-84d4-34ebbe1fb14bn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/2fb6de53-e6de-44ce-84d4-34ebbe1fb14bn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1688458718678000&usg=AOvVaw3enImdJc5vHC8B7ZjsF3jY">https://groups.google.com/d/msgid/cp2k/2fb6de53-e6de-44ce-84d4-34ebbe1fb14bn%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br>--
<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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/85703e25-7ad2-4fc4-9a82-43261d13ab81n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/85703e25-7ad2-4fc4-9a82-43261d13ab81n%2540googlegroups.com&source=gmail&ust=1688458718678000&usg=AOvVaw0TDPSFLbtKwkxQvC3BRK-A">https://groups.google.com/d/msgid/cp2k/85703e25-7ad2-4fc4-9a82-43261d13ab81n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/85703e25-7ad2-4fc4-9a82-43261d13ab81n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/85703e25-7ad2-4fc4-9a82-43261d13ab81n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1688458718678000&usg=AOvVaw1I73ka-_EsFm5PJxaDlRRm">https://groups.google.com/d/msgid/cp2k/85703e25-7ad2-4fc4-9a82-43261d13ab81n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br>
<br>--
<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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/f8c23295-cf28-4f92-a65a-aff594ea3c03n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/f8c23295-cf28-4f92-a65a-aff594ea3c03n%2540googlegroups.com&source=gmail&ust=1688458718678000&usg=AOvVaw054SXfqSyaQroRlhvtBq62">https://groups.google.com/d/msgid/cp2k/f8c23295-cf28-4f92-a65a-aff594ea3c03n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/f8c23295-cf28-4f92-a65a-aff594ea3c03n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/f8c23295-cf28-4f92-a65a-aff594ea3c03n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1688458718678000&usg=AOvVaw0uwimSgEw9SMJizGKEG8pc">https://groups.google.com/d/msgid/cp2k/f8c23295-cf28-4f92-a65a-aff594ea3c03n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></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/9f3303b1-7a80-49a0-b2fc-cb06f697d0e6n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/9f3303b1-7a80-49a0-b2fc-cb06f697d0e6n%40googlegroups.com</a>.<br />