DFPT_CGWF WARNING : d2edt2 is zero, skipping update

Hey everyone, I’m getting a warning when calculating the nonlinear optic susceptibility that wasn’t mentioned in the tutorial(https://docs.abinit.org/tutorial/nlo/#1-ground-state-properties-of-alp-and-general-parameters). Have any of you encountered this issue before? And how can I get rid of this warning?Thanks!

Dear Jiang,

Essentially WARNING messages in Abinit should be almost harmless to the physical properties you are calculating, otherwise you have the ERROR messages that put a hard stop to your calculation when an incongruent set of parameters is provided.

In your case, if you want to make sure about what’s going on, you can increase the amount of printed information with prtvol (a value of 11 or higher should do the job). You will be able to see for which band the d2edt2 value is zero.
My suspicion is data the number of bands or nbdbuf have been changed from one dataset to another.
Let me know!
Bogdan

Dear Bogdan

You are right. I tested nbdbuf and the warning disappeared when nbdbuf was set to 0. Whether the warning is present or not, the output results are the same, and the war
ning does not affect the final result. Unfortunately, even when I increased the accuracy to 30 Ha, 161616 k points, the nonlinear optical susceptibility still had a significant difference from the literature(https://journals.aps.org/prb/pdf/10.1103/PhysRevB.71.125107). I have attached the input file for AlAs, and I would be grateful if someone could tell me where the problem is.
tnlo_2.abi (5.7 KB)
tnlo_2.abo (179.2 KB)
tnlo_4.abi (1.3 KB)
tnlo_4.abo (17.0 KB)

Dear Jiang,
You should (always) test the convergence w.r.t. kpoints and ecut, here kpoints should be quite enough (it might be more than necessary) but ecut=30Ha might be a low end value for norm-conserving pseudopotentials. Non-linear responses are quite “delicate” to obtain and you have to be sure that you have enough plane waves and k-points.
In any case, convergence test should always be done for whatever property you calculate.
Best wishes,
Eric

Dear Everyone,
I’m sorry.My previous statement was incorrect. After tests, it is necessary to strictly limit nbdbuf to 0 in the calculation. Unfortunately, even when I set the ecut to 70 Ha and the ngkpt to 161616, I still cannot obtain reasonable results.I hope someone can give me advice.Thanks!
Best wishes,
Jiang
tnlo_4.abi (1.3 KB)
tnlo_4.abo (17.0 KB)
tnlo_6.abi (7.3 KB)
tnlo_6.abo (256.9 KB)