Search found 43 matches

by andyamygto
Fri Dec 04, 2020 7:02 am
Forum: Response calculations
Topic: Strange divergence in phonon calculation
Replies: 15
Views: 16553

Re: Strange divergence in phonon calculation

Dear Eric,

Thanks for helping report the issue. I look forward the upcoming update.

Best wishes,
Andy
by andyamygto
Fri Nov 20, 2020 9:31 pm
Forum: Response calculations
Topic: Strange divergence in phonon calculation
Replies: 15
Views: 16553

Re: Strange divergence in phonon calculation

Dear Eric, I have tested the nbands from 35 to 95 with 10 nbands increase every time. Seems like the issue still remains. Also, I have tested the relativity high energy cutoff with the PAW dataset (up to 45 Ha with pawecutdg=90) The problem still exists. And there is something that might help. I jus...
by andyamygto
Sat Nov 14, 2020 2:46 pm
Forum: Response calculations
Topic: Strange divergence in phonon calculation
Replies: 15
Views: 16553

Re: Strange divergence in phonon calculation

Dear Eric, Some update.. Same divergence happened even when I use the Atomic PAW from GPAW https://wiki.fysik.dtu.dk/gpaw/setups/setups.html So I think this might be the bug for the PAW+DFPT calculation? However, In another case SnO2 rutile, the PAW + DFPT working well. So, form my perspective, eith...
by andyamygto
Fri Nov 13, 2020 5:58 pm
Forum: Response calculations
Topic: Strange divergence in phonon calculation
Replies: 15
Views: 16553

Re: Strange divergence in phonon calculation

Dear Eric, No matter I use high or low tsmear (0.0001 to 0.1). the divergence will happened. Also, I have tried the GBRV PAW data set before, but the GRBV PAW data set doesn't include the spin polarize info in their dataset, so I wonder if that may cause the problem? All the ground state converged h...
by andyamygto
Wed Nov 04, 2020 9:14 pm
Forum: Response calculations
Topic: Strange divergence in phonon calculation
Replies: 15
Views: 16553

Re: Strange divergence in phonon calculation

Dear Eric, Sorry for the late response. Thanks for suggesting that the tsmear might be large for my calculation. However, the divergence also happened when I use the occopt=1 with spinmagntarget = 0 to define the AFM polarize for my system. I wonder if it is the PAW's problem or there are some sort ...
by andyamygto
Wed Nov 04, 2020 9:08 pm
Forum: Many-body perturbation theory
Topic: "ug for (band, ik_ibz, spin): not stored in memory" ?
Replies: 4
Views: 10089

Re: "ug for (band, ik_ibz, spin): not stored in memory" ?

Dear Users,

I also met this problem in 9.2.1 version.
Is this still an issue?

Regards,
Andy
by andyamygto
Mon Nov 02, 2020 9:39 am
Forum: Response calculations
Topic: Strange divergence in phonon calculation
Replies: 15
Views: 16553

Re: Strange divergence in phonon calculation

Dear Eric, The system is CrI3 monolayer but the error is the same, which is the weird divergence when doing perturbation. I have test for both GGA and LDA with JTH-LDA-atomicdata-1.1 I have no idea what is causing the problem. Here is my input (I was still using the files file in 9.2.1 LUL) In my in...
by andyamygto
Mon Nov 02, 2020 3:59 am
Forum: Response calculations
Topic: Strange divergence in phonon calculation
Replies: 15
Views: 16553

Re: Strange divergence in phonon calculation

Hi,

I met the same problem in the latest product version 9.2.1
Is there any suggestion or information about this error?
I was using the newest JTH PSP from the Abinit site.

Regards,
Andy
by andyamygto
Sun Sep 13, 2020 6:03 pm
Forum: Response calculations
Topic: Negative frequency with denser q-points
Replies: 16
Views: 16451

Re: Negative frequency with denser q-points

Dear Lily you can grep the eigenvector with the corresponding mode and use xcrysden or VESTA to check the vibration mode. In my experiment, LDA usually perform better with DFPT method (I have no idea why is that a case...) So it's not surprising that you said changing to LDA will fix the negative fr...
by andyamygto
Sat Sep 05, 2020 9:10 pm
Forum: Response calculations
Topic: Negative frequency with denser q-points
Replies: 16
Views: 16451

Re: Negative frequency with denser q-points

Dear Lily, If the negative frequency is happened at gamma point and you have tested most the converge variables, It might mean that the system is not stable under GGA. The force I mean is that the force acting on each atom. In particular case, reading the force may find something interesting. Have y...
by andyamygto
Fri Sep 04, 2020 11:20 pm
Forum: Response calculations
Topic: Negative frequency with denser q-points
Replies: 16
Views: 16451

Re: Negative frequency with denser q-points

Dear Lily, Here is few comment and question... 1, which q-point you get the negative frequency? One thing you can do is to print out the phonon eigenvector to see which kind of vibration cause the negative frequency (with q-vector, it usually hard to distinguish the vibration symmetry, but still wor...
by andyamygto
Wed Sep 02, 2020 11:39 am
Forum: Response calculations
Topic: Negative frequency with denser q-points
Replies: 16
Views: 16451

Re: Negative frequency with denser q-points

Dear Lily, You should grep the frequency from the DFPT calculation at each q-point (the abinit output, not the anaddb output) If you are on the linux you can use command to get the frequency calculated form the DFPT at each point: $grep -A 3 cm 'your_result'.out The converge of the vacuum layer shou...
by andyamygto
Sun Aug 30, 2020 8:24 pm
Forum: Response calculations
Topic: Negative frequency with denser q-points
Replies: 16
Views: 16451

Re: Negative frequency with denser q-points

Dear lilymali03, Here is my suggestion. in anaddb files, test the phonon dispersion with the variable rifcsph. The variable decide the Radius of the Interatomic Force Constant sphere. If you use the denser q-point, the short range force range is larger. Which might cause some sort of mistake. But th...
by andyamygto
Fri Jul 17, 2020 3:57 pm
Forum: Ground state
Topic: Non Linear calculation input for ferromagnetic system
Replies: 9
Views: 5625

Re: Non Linear calculation input for ferromagnetic system

Dear all,

Another question I wonder,
Is that means it is impossible to calculate the non-linear properites within DFPT method?
or I can just use occopt = 1 for metal? I think it makes wrong since i can not define the fermi level well.

If there are some advise I will be thankful.

Regards,
Andy
by andyamygto
Sun Apr 26, 2020 9:06 pm
Forum: Ground state
Topic: Non Linear calculation input for ferromagnetic system
Replies: 9
Views: 5625

Re: Non Linear calculation input for ferromagnetic system

Dear Lucas,

Thanks for the info. I'll wait for the good news.
You guys are doing amazing job.

Best wishes,
Andy
by andyamygto
Thu Apr 16, 2020 5:08 pm
Forum: Ground state
Topic: Non Linear calculation input for ferromagnetic system
Replies: 9
Views: 5625

Re: Non Linear calculation input for ferromagnetic system

Dear Eric,

I see. it's unfortunate to know that FM case cannot be done...
Hope that spinmagntarget correction will be include in the future. (or the metallic occupation)
I'll wait the comment form Lucas.

Regards,
Andy
by andyamygto
Thu Apr 09, 2020 3:31 am
Forum: Ground state
Topic: Non Linear calculation input for ferromagnetic system
Replies: 9
Views: 5625

Re: Non Linear calculation input for ferromagnetic system

Ok. I think i will just wait for the good news
Hope that everything is alright.

Cheers
Andy
by andyamygto
Wed Apr 08, 2020 3:53 am
Forum: Ground state
Topic: Non Linear calculation input for ferromagnetic system
Replies: 9
Views: 5625

Re: Non Linear calculation input for ferromagnetic system

Dear Eric,
Thanks for the reply.
I am looking forward it!
I wonder if there is any other way to set the spin without using spinmagntarget?

Cheers,
Andy
by andyamygto
Tue Apr 07, 2020 7:35 pm
Forum: Response calculations
Topic: Phonon DOS [SOLVED]
Replies: 19
Views: 11413

Re: Phonon DOS [SOLVED]

Hi...just some comment and different view I cannot tell if your idea is correct since i am not familiar with the PDOS calculation. But i think you can do the simple test by using the small grid to prove your point. Actually I think that the much lower frequency states is not causing by the missing N...
by andyamygto
Tue Apr 07, 2020 7:11 pm
Forum: Ground state
Topic: Non Linear calculation input for ferromagnetic system
Replies: 9
Views: 5625

Non Linear calculation input for ferromagnetic system

Dear Abinit users The new implementation of DFPT for non-linear response (usepead = 0) works for collinear spin polarized system. but in the test file, all the test of this implementation use the variable spinmagtarget=0. I think for that case, the system must be Anti-ferromagnetic, right? -------- ...
by andyamygto
Tue Mar 31, 2020 4:35 am
Forum: Response calculations
Topic: The Raman(NLO) calculation for the ferromagnetic system
Replies: 5
Views: 4578

Re: The Raman(NLO) calculation for the ferromagnetic system

Use spinmagtarget, it's easiest. Do you confirm your system is insulating? Otherwise there is no point... Hi I just noticed that abinit release the new version. Just wonderIng is there any other sight or solution to the problem I met? I haven’t try the new version(9.0) yet But I want to make sure i...
by andyamygto
Mon Mar 30, 2020 2:19 pm
Forum: Link to other codes
Topic: Pseudopotential file pspxc not equal to input ixc [SOLVED]
Replies: 2
Views: 11674

Re: Pseudopotential file pspxc not equal to input ixc [SOLVED]

Dear Eric

Thanks for the info.The calculation successes and work perfectly,.

Cheers,
Andy Hsu
by andyamygto
Wed Mar 25, 2020 3:34 pm
Forum: Response calculations
Topic: Phonon DOS [SOLVED]
Replies: 19
Views: 11413

Re: Phonon DOS [SOLVED]

Hi From my experience, the ngkpt in the main run is okay if it is commensurate to the q-grids. if you want to use the fine k-grid with only 2x2x2 q-grid, I suggest to use integral multiple ngkpt of the q-grids or it might cause the numerically error.(Image that the q-point is not on the k-grid, so t...
by andyamygto
Wed Mar 25, 2020 1:48 pm
Forum: Link to other codes
Topic: Pseudopotential file pspxc not equal to input ixc [SOLVED]
Replies: 2
Views: 11674

Pseudopotential file pspxc not equal to input ixc [SOLVED]

Hi I was doing some simple test with ABINIT8.10.3 and using the ONCVPSP-PW-PDv0.4 pseudo potentials https://github.com/abinit/pseudo_dojo/tree/master/pseudo_dojo/pseudos And in my calculation, the code says that --- !WARNING src_file: m_pspini.F90 src_line: 1088 message: | Pseudopotential file pspxc...
by andyamygto
Wed Mar 25, 2020 12:40 am
Forum: Response calculations
Topic: Phonon DOS [SOLVED]
Replies: 19
Views: 11413

Re: Phonon DOS [SOLVED]

Hi Yes the qph1l list of qpts doesn't need to match the the main run, but the anaddb input ngqpt need to be the q-grid you use. Once you have the information at the q-grid, the rest of the phonon dispersion can be form by the Fourier transform. so you can choose any path in the Brillouin zone to for...