<div dir="ltr"><div>Dear CP2K users,</div><div><br></div>
<div>I'm trying to set up a QS calculation of EPR magnetic properties of
 NaCl type crystals into which doped transition metal ions. However we 
encountered a convergence problem as soon as the "current calculation " 
job starts. </div><div><br></div>
<div>The message is followed:</div><div><br></div><div>------------------------------<wbr>------------------------------<wbr>--------------------<br></div>
<div>*** Start current Calculation ***<br>
Inizialization of the current environment<br>
To get CURRENT parameters within PBC you need localized zero order orbitals <br>
CURRENT| Gauge used  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â  R_AND_STEP_FUNCTION<br>
CURRENT| Use old gauge code  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â Â T<br>
CURRENT| Compute chi for PBC calculation  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â  T<br>
CURRENT| Orbital center used  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â  WANNIER<br>
CURRENT| Compute  Â 260 selected response functions out of  Â 260 for spin  1<br>
CURRENT| There is a total of  Â  260 (clustered) center(s) for spin  1<br>
CURRENT| Compute  Â 259 selected response functions out of  Â 259 for spin  2<br>
CURRENT| There is a total of  Â  259 (clustered) center(s) for spin  2<br>
CURRENT| Calculation of the p and (r-d)xp operators applied to psi0<br>
CURRENT| current_operators: CheckSum L_x = 0.1075579152983971E+08<br>
CURRENT| current_operators: CheckSum L_y = 0.1075613503806711E+08<br>
CURRENT| current_operators: CheckSum L_z = 0.1074174864868970E+08<br>
CURRENT| current_operators: CheckSum P_x = 0.2095122596008357E+03<br>
CURRENT| current_operators: CheckSum P_y = 0.2095122597965540E+03<br>
CURRENT| current_operators: CheckSum P_z = 0.2094514062936749E+03<br>
*** Self consistent optimization of the response wavefunctions ***</div>
<div>Read response wavefunction from restart file naclrh-RESTART-nmr_p.lr<br>
Read response wavefunction from restart file naclrh-RESTART-nmr_rxp.lr<br>
Read response wavefunction from restart file naclrh-RESTART-nmr_p.lr<br>
Read response wavefunction from restart file naclrh-RESTART-nmr_rxp.lr<br>
Read response wavefunction from restart file naclrh-RESTART-nmr_p.lr<br>
Read response wavefunction from restart file naclrh-RESTART-nmr_rxp.lr<br>
Response to the perturbation operator P_x<br>
Iteration  Â  Method  Â Restart  Â Â  Â Stepsize  Â Â  Â Convergence  Â Â  Â Â  Â Time<br>
------------------------------<wbr>------------------------------<wbr>--------------------<br>
1  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.00E+00  Â Â  Â 0.0005645548  Â Â  Â Â Â 1.88<br>
2  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.65E-02  Â Â  Â 0.0004104050  Â Â  Â Â Â 2.83<br>
3  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.95E-02  Â Â  Â 0.0005342940  Â Â  Â Â Â 3.78<br>
4  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.65E-02  Â Â  Â 0.0005775678  Â Â  Â Â Â 4.81<br>
5  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.12E-01  Â Â  Â 0.0002175021  Â Â  Â Â Â 5.76<br>
6  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.15E-01  Â Â  Â 0.0003167685  Â Â  Â Â Â 6.71<br>
7  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.11E-01  Â Â  Â 0.0002438622  Â Â  Â Â Â 7.85<br>
8  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.17E-01  Â Â  Â 0.0001200927  Â Â  Â Â Â 8.80<br>
9  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.44E-01  Â Â  Â 0.0003625334  Â Â  Â Â Â 9.88<br>
10  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.16E-01  Â Â  Â 0.0005849186  Â Â  Â  10.89<br>
11  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.77E-02  Â Â  Â 0.0002077820  Â Â  Â  11.85<br>
12  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.12E-01  Â Â  Â 0.0001416913  Â Â  Â  12.86<br>
13  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.27E-01  Â Â  Â 0.0001917244  Â Â  Â  13.92<br>
14  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.13E-01  Â Â  Â 0.0001881651  Â Â  Â  14.88<br>
15  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.79E-02  Â Â  Â 0.0001655726  Â Â  Â  15.97<br>
16  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.37E-01  Â Â  Â 0.0001248280  Â Â  Â  16.97<br>
17  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.25E-01  Â Â  Â 0.0001401467  Â Â  Â  17.92<br>
18  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.10E-01  Â Â  Â 0.0002603074  Â Â  Â  18.96<br>
19  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.33E-01  Â Â  Â 0.0001861282  Â Â  Â  19.91<br>
20  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.13E-01  Â Â  Â 0.0001875534  Â Â  Â  20.86<br>
21  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.32E-01  Â Â  Â 0.0002338661  Â Â  Â  21.81<br>
22  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.15E-01  Â Â  Â 0.0002797592  Â Â  Â  22.88<br>
23  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.75E-02  Â Â  Â 0.0001564825  Â Â  Â  23.97<br>
24  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.44E-01  Â Â  Â 0.0004296774  Â Â  Â  25.36<br>
25  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.43E-01  Â Â  Â 0.0003670127  Â Â  Â  26.34<br>
26  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.12E-01  Â Â  Â 0.0001827137  Â Â  Â  27.47<br>
27  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.13E-01  Â Â  Â 0.0002811690  Â Â  Â  28.43<br>
28  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.18E-01  Â Â  Â 0.0002879992  Â Â  Â  29.40<br>
29  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.20E-01  Â Â  Â 0.0002451334  Â Â  Â  30.36<br>
30  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.83E-02  Â Â  Â 0.0002721858  Â Â  Â  31.48<br>
31  Â Â  Â Â Â PCG  Â Â  Â  T  Â Â  Â Â  Â 0.40E-01  Â Â  Â 0.0005828498  Â Â  Â  32.75<br>
32  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.92E-02  Â Â  Â 0.0006058013  Â Â  Â  34.29<br>
33  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.34E-01  Â Â  Â 0.0001909972  Â Â  Â  35.23<br>
34  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.15E-01  Â Â  Â 0.0004781541  Â Â  Â  36.34<br>
35  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.11E-01  Â Â  Â 0.0002734769  Â Â  Â  37.51<br>
36  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.25E-01  Â Â  Â 0.0011724429  Â Â  Â  38.89<br>
37  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.28E-01  Â Â  Â 0.0009580690  Â Â  Â  40.04<br>
38  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.25E-01  Â Â  Â 0.0002626740  Â Â  Â  41.10<br>
39  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.13E-01  Â Â  Â 0.0004089794  Â Â  Â  42.10<br>
40  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.11E-01  Â Â  Â 0.0001815872  Â Â  Â  43.13<br>
41  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.23E-01  Â Â  Â 0.0005746123  Â Â  Â  44.21<br>
42  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.63E-02  Â Â  Â 0.0004449090  Â Â  Â  45.22<br>
43  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.34E-01  Â Â  Â 0.0007450151  Â Â  Â  46.34<br>
44  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.59E-02  Â Â  Â 0.0003882099  Â Â  Â  47.29<br>
45  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.67E-01  Â Â  Â 0.0003010523  Â Â  Â  48.61<br>
46  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.23E-01  Â Â  Â 0.0009782303  Â Â  Â  49.56<br>
47  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.13E-01  Â Â  Â 0.0009607723  Â Â  Â  51.17<br>
48  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.38E-01  Â Â  Â 0.0010698243  Â Â  Â  52.50<br>
49  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.11E-01  Â Â  Â 0.0002619895  Â Â  Â  53.82<br>
50  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.29E-01  Â Â  Â 0.0006729204  Â Â  Â  54.83<br>
The linear solver didnt converge! Maximum number of iterations reached.</div>
<div>------------------------------<wbr>------------------------------<wbr>--------------------<br><br><br>
I was trying to decrease the keyword "ENERGY_GAP 2.0E-004" and increase "MAX_ITER", but the program was aborted afterward.</div><div><br>
The error message is followed:</div><div><br>------------------------------<wbr>------------------------------<wbr>--------------------<br>
 *** Start current Calculation ***  Â Â  Â Â  Â Inizialization of the current environment<br>
 To get CURRENT parameters within PBC you need localized zero order orbitals <br>
 CURRENT| Gauge used  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â  R_AND_STEP_FUNCTION<br>
 CURRENT| Use old gauge code  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â Â T<br>
 CURRENT| Compute chi for PBC calculation  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â  T<br>
 CURRENT| Orbital center used  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â  WANNIER<br>
 CURRENT| Compute  Â 260 selected response functions out of  Â 260 for spin  1<br>
 CURRENT| There is a total of  Â  260 (clustered) center(s) for spin  1<br>
 CURRENT| Compute  Â 259 selected response functions out of  Â 259 for spin  2<br>
 CURRENT| There is a total of  Â  259 (clustered) center(s) for spin  2<br>
 CURRENT| Calculation of the p and (r-d)xp operators applied to psi0<br>
 CURRENT| current_operators: CheckSum L_x = 0.1075579154261990E+08<br>
 CURRENT| current_operators: CheckSum L_y = 0.1075613501246297E+08<br>
 CURRENT| current_operators: CheckSum L_z = 0.1074174867173965E+08<br>
 CURRENT| current_operators: CheckSum P_x = 0.2095122596008357E+03<br>
 CURRENT| current_operators: CheckSum P_y = 0.2095122597965540E+03<br>
 CURRENT| current_operators: CheckSum P_z = 0.2094514062936749E+03<br>
 *** Self consistent optimization of the response wavefunctions ***<br>
 <br>
 Read response wavefunction from restart file naclrh-RESTART-nmr_p.lr<br>
 Read response wavefunction from restart file naclrh-RESTART-nmr_rxp.lr<br>
 Read response wavefunction from restart file naclrh-RESTART-nmr_p.lr<br>
 Read response wavefunction from restart file naclrh-RESTART-nmr_rxp.lr<br>
 Read response wavefunction from restart file naclrh-RESTART-nmr_p.lr<br>
 Read response wavefunction from restart file naclrh-RESTART-nmr_rxp.lr<br>
 Response to the perturbation operator P_x<br>
 Iteration  Â  Method  Â Restart  Â Â  Â Stepsize  Â Â  Â Convergence  Â Â  Â Â  Â Time<br>
 -----------------------------<wbr>------------------------------<wbr>---------------------<br>
 1  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.00E+00  Â Â  Â 0.0111934742  Â Â  Â Â Â 5.20<br>
 2  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.26E+00  Â Â  Â 0.0168502715  Â Â  Â Â Â 6.15<br>
 3  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.80E-01  Â Â  Â 0.0211182672  Â Â  Â Â Â 7.11<br>
 4  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.42E-01  Â Â  Â 0.0141167009  Â Â  Â Â Â 8.05<br>
 5  Â Â  Â Â Â PCG  Â Â  Â  F  Â Â  Â Â  Â 0.66E-01  Â Â  Â 0.0211785481  Â Â  Â Â Â 8.98<br>
 *****************************<wbr>******************************<wbr>********************<br>
 *  Â ___  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â Â *<br>
 *  /  Â \  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â  *<br>
 * [ABORT]  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â *<br>
 *  \___/  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â  tr(Ap_j*p_j) < 0  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â Â *<br>
 *  Â  |  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â *<br>
 *  O/|  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â *<br>
 * /| |  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â *<br>
 * / \  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â Â  Â  qs_linres_methods.F:520 *<br>
 *****************************<wbr>******************************<wbr>********************<br>
 <br>
 ===== Routine Calling Stack ===== <br>
 5 linres_solver<br>
 4 current_response<br>
 3 linres_calculation_low<br>
 2 linres_calculation<br>
 1 CP2K<br> 
------------------------------<wbr>------------------------------<wbr>--------------------<br><br><br>
 <br>
  The question is that are there any parameters need to be changed to 
make the calculation work? Or simply i just missed some important steps?<br>
  <br></div><div>  Input Â files are attached.</div><div><br></div></div>