Page 1 of 2

DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Tue Jul 30, 2019 10:33 am
by onion2440
Hi everyone,
I want to calculate the phonon dispersion of Nb with spin orbit coupling. The Psp is GBRV. However, when abinit 8.10.3 calculate the q point : qpt2 0.00000000000000D+00 0.00000000000000D+00 0.00000000000000D+00, the abinit crashed and give me the informations below:

--- !BUG
src_file: m_paw_occupancies.F90
src_line: 693
mpi_rank: 8
message: |
pawrhoij(:)%rhoijim must be allocated!
...

anyone can help me?

Best
yang

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Wed Jul 31, 2019 3:04 am
by onion2440
I changed Psp to hgh, and the problems solved.

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Wed Jul 31, 2019 5:38 pm
by ebousquet
Dear Yang,
Are you using GGA or LDA in your calculation with PAW potential?
Otherwise, shifting to HGH norm-conserving means that the ecut must be much larger than for PAW. If it is for production I would advise you to use the pseudodojo norm-conserving potentials:
http://www.pseudo-dojo.org
In "type" when you have FR means fully relativistic with SOC included.
Best wishes,
Eric

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Fri Aug 02, 2019 2:41 pm
by onion2440
ebousquet wrote:Dear Yang,
Are you using GGA or LDA in your calculation with PAW potential?
Otherwise, shifting to HGH norm-conserving means that the ecut must be much larger than for PAW. If it is for production I would advise you to use the pseudodojo norm-conserving potentials:
http://www.pseudo-dojo.org
In "type" when you have FR means fully relativistic with SOC included.
Best wishes,
Eric


Dear Eric,
I'm using LDA with PAW potential, since it seems DFPT method only support LDA PAW now.
Thank you for your suggestion, I will change to NC potetial.

Best,
yang

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Mon Aug 05, 2019 10:20 am
by ebousquet
Dear Yang,
DFPT works for both NC and PAW (if using the latest version) the only special requirement is to use pawxcdev=0 if you want to use PAW+GGA DFPT.
Best wishes,
Eric

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Thu Aug 08, 2019 5:42 am
by onion2440
ebousquet wrote:Dear Yang,
DFPT works for both NC and PAW (if using the latest version) the only special requirement is to use pawxcdev=0 if you want to use PAW+GGA DFPT.
Best wishes,
Eric


Dear Eric,
I'm confusing with the special requirement of set pawxcdev=0 when we use PAW+GGA DFPT.
The document said (https://docs.abinit.org/variables/paw/#pawxcdev):
Be careful: GGA requires pawxcdev > 0

Is there any special case here?

Best wishes,
Yang

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Thu Aug 08, 2019 9:21 pm
by ebousquet
onion2440 wrote:
ebousquet wrote:Dear Yang,
I'm confusing with the special requirement of set pawxcdev=0 when we use PAW+GGA DFPT.
The document said (https://docs.abinit.org/variables/paw/#pawxcdev):
Be careful: GGA requires pawxcdev > 0

Hum, strange, it should be pawxcdev=0 to do DFPT+PAW+GGA...

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Fri Aug 09, 2019 11:40 am
by ebousquet
I confirm that it should be pawxcdev=0! It'll be corrected, thanks for reporting it.
Eric

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Mon Oct 07, 2019 3:17 pm
by Hola
Hello Eric,

I want to do a DFPT calculation with a PAW pseudopotential (for calculating phonon spectra).
With the latest version of abinit (8.10.3), I should mentionne in the input file
pawxcdev 0, otherwise the calculation crash.

With an older version of abinit (8.8.3), installed in a supercomputer, the calculation is running
well, whatever I mentionne pawxcdev 0 or pawxcdev 1, in my input file. Is it normal ?

Yours sincerely,

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Sun Oct 20, 2019 8:44 pm
by NPike
Hello Hola,

As we release new versions of ABINIT it is possible that the requirements on input variables change. The various help menus and variable documentation always corresponds to the most recent release of ABINIT.

-Nick

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Tue Nov 26, 2019 4:49 pm
by andyamygto
Hello Yang,

I met the same problem when i doing the FULL DFPT calculation with PAW-PP
I wonder if anyone have an idea how to fix this problem without change the PP-file(since i want to do the DFT+U+DFPT,so i need to use the PAW-PP for my calculation)

Best wishes,
Andy

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Wed Dec 04, 2019 3:11 pm
by ebousquet
Dear Andy,
What is the error message you get?
With DFP+U there were a small bug in the code and it was frozen in the 8.10.3. This has been solved since a while but not sure if this is fixed in 8.10.3...
Best wishes,
Eric

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Sun Dec 08, 2019 1:52 pm
by andyamygto
Dear Eric,
Here is the WARNING massage i got

pawio_print_ij: WARNING -
The matrix seems to have high value(s) !
( 64 components have a value greater than 50.0).
It can cause instabilities during SCF convergence.
Action: you should check your atomic dataset (psp file)
and look for "high" projector functions...

I haven't apply the Hubbard U for my calculation yet.
So I think the problem is cause by the FULL DFPT + PAW
Since it cause the divergence in my 2DEtotal during the perturbation,so i just quit the job for now
and change to the NC-PP.

Best wishes,
Andy

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Fri Dec 20, 2019 10:28 am
by torrent
Dear all,

Some comments...
* DFPT+PAW+U is not guaranteed, especially in 8.10.3.
This is an on-going development... that will be better in v9.0 (january).

* The warning "The matrix seems to have high value!" is just a warning. It means -- as explained -- that the PAW datasets is a bit hard to converge because abinit has difficulties toorthogonalize projectors for your particular system.
Universal atomic datatsets do not exists although we are doing our best to build as transferable as possible datasets.
And, if the dataset is from GBRV table, this is another story.
Anyway, as soon as your calculation converges, you don't have to worry about this warning.

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Fri Dec 20, 2019 6:36 pm
by andyamygto
torrent wrote:
Fri Dec 20, 2019 10:28 am
Dear all,

Some comments...
* DFPT+PAW+U is not guaranteed, especially in 8.10.3.
This is an on-going development... that will be better in v9.0 (january).

* The warning "The matrix seems to have high value!" is just a warning. It means -- as explained -- that the PAW datasets is a bit hard to converge because abinit has difficulties toorthogonalize projectors for your particular system.
Universal atomic datatsets do not exists although we are doing our best to build as transferable as possible datasets.
And, if the dataset is from GBRV table, this is another story.
Anyway, as soon as your calculation converges, you don't have to worry about this warning.
Dear torrent

Thanks a lot for the comment.
I wonder is there any way to deal with the divergence total energy (2DEtotal) in perturbation without changing the PAW to NCPP?
I try both the general Atomic dataset and the one from GBRV table . They both leads to the weird divergence during the perturbation. (2DEtotal up to 10^30 Hatree)
BTW, I am looking forward the v9.0. You guys are doing amazing job!

Regards,
Andy

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Sat Jan 18, 2020 4:01 pm
by ebousquet
Dear Andy,
Regarding the divergence of the SCF, it can be specific for each case. The first thing would be to play with the mixing parameters (diemix, diemixmag if magnetism, nline, etc) see if reducing them could improve the situation.
Bets wishes,
Eric

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Sun Feb 02, 2020 10:45 am
by andyamygto
Dear Eric,

Thanks for the advise , i was trying to tune the mixing parameters and see the results.
It doesn't seems to be better while doing perturbation. I will keep trying on this.

BTW, I was wonder is there any tip for setting these kind of parameters?

Regards,
Andy

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Sun Feb 09, 2020 12:09 am
by ebousquet
Dear Andy,
How to tune mixing parameter is something coming from experience but in general the best to have them as large as possible such that you reduce them step by step to see when this can improve your case.
By the way, at this stage it'll be hard to help you without knowing which system are you running on, I've made several tests of DFPT+PAW (with and without U) and most of them go smoothly without much tuning mixing. Is your system converging well during ground state or relaxation?
Best wishes,
Eric

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Sun Feb 09, 2020 12:54 pm
by andyamygto
Dear Eric,
Thanks for the reply.
Sorry, I was busy from another project recently so i haven't do many tests, but there is still some progress and will be take on later.
The system i am running is a monolayer CrI3 , my system is well converge during ground state and relaxation.
The only thing might cause the problem is that for the monolayer CrI3, the ground state energy with AFM and FM state is extreme close( 0.001 Ha difference).
But there's another topic i was trying to solve is that i cannot do the FM calculation in the FULL DFPT
viewtopic.php?p=12726#p12726

I don't know if that cause the problem or not.
I will try to do more test about mixing parameter.

Regards,
Andy

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Tue Feb 11, 2020 3:17 am
by ebousquet
andyamygto wrote:
Sun Feb 09, 2020 12:54 pm
The only thing might cause the problem is that for the monolayer CrI3, the ground state energy with AFM and FM state is extreme close( 0.001 Ha difference).
Yes this can be a problem, when you are close to phase transitions, the energy landscape flattens and things have tendency to diverge... diemix and diemixmag could do the job here, but they might need to be small and requiring more SCF... You can also increase nline to 10 or 12 to help.

For the FM, you could switch to metallic occupation (and putting a small tsemar if insulating) as stated in the post you mention, did you try it?

Best wishes,
Eric

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Tue Feb 11, 2020 4:29 am
by andyamygto
ebousquet wrote:
Tue Feb 11, 2020 3:17 am
For the FM, you could switch to metallic occupation (and putting a small tsemar if insulating) as stated in the post you mention, did you try it?
Yes i have tried this way at first. but seems like the code won't let me do the calculation.
It's sending the ERROR message:

chkint_prt: ERROR -
Context: the value of the variable usepead is 0.
Context: the value of the variable optdriver is 5.
The value of the input variable occopt is 4, while it must be
equal to 1
Action: you should change one of the input variables occopt,
usepead or optdriver.

So i tried the spinmagntarget later on, but both way i can't get the calculation down. Did i miss something?
I will try the diemix and diemixmag etc... first. Hope i can get a good result.

Regards,
Andy

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Sat Feb 15, 2020 8:23 pm
by ebousquet
andyamygto wrote:
Tue Feb 11, 2020 4:29 am
chkint_prt: ERROR -
Context: the value of the variable usepead is 0.
Context: the value of the variable optdriver is 5.
The value of the input variable occopt is 4, while it must be
equal to 1
Action: you should change one of the input variables occopt,
usepead or optdriver.
Hummm, looks strange in my case it is running without problem with the 8.10.3. Sounds like there are flags that might be conflicting in your case. Do you mind to show your input?

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Tue Feb 18, 2020 1:00 pm
by andyamygto
ebousquet wrote:
Sat Feb 15, 2020 8:23 pm
Hummm, looks strange in my case it is running without problem with the 8.10.3. Sounds like there are flags that might be conflicting in your case. Do you mind to show your input?
Hi, here is my input.
I am running this job in abinit 8.10.3
I was using LDA psp8 for the calculation, but I think the psedopotential will not effect this error.

regards,
Andy

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Thu Feb 20, 2020 2:59 pm
by ebousquet
Dear Andy,
OK, I see that you are doing more than the phonons since in your input you also calculate the non-linear responses. I guess up to dataset #4 things are running but you get the crash at one of the following dataset, could you confirm this and tell me at which dataset you have the error message?
Not sure if the non-linear response is implemented with SOC, I'm going to ask to the dev of this part.
Best wishes,
Eric

Re: DFPT phonon with soc problem:rhoijim must be allocated!

Posted: Fri Feb 21, 2020 1:22 pm
by andyamygto
Dear Eric,

Yes I got the error message for dataset 5,6,7 which is the non-linear response part.

I have done a lot research and seems like the DFPT doesn't support the SOC.
But for my calculation, I haven't turn on the SOC yet (nspinor).
I just want to do the ferromagnetic calculation with DFPT but it won't work.
It should be able to do only the spin polarized calculations without the SOC, right?(Since I was getting all the information from the tutorial. please correct me if I am wrong.)

I would appreciate it if you can give me more information.

Best wishes,
Andy