RPA calculation hangs at dipoles

You can post here problems arising when using the devel release of Yambo. Issues as parallelization strategy, performance issues and other technical aspects realted to the new release.

Moderators: Daniele Varsano, andrea.ferretti, andrea marini, Conor Hogan, myrta gruning

RPA calculation hangs at dipoles

Postby vjhalani » Wed Apr 17, 2019 1:38 am

Hi all,

Using Yambo 4.4, I am trying to do an RPA calculation (yambo -o c) on GaAs on pretty fine grids, and after a certain density, I'm finding that the calculation gets stuck after reading the kb_pp* files and is unable to write the dipole databases. This happens at a 24^3 grid and 30^3 grids, but curiously, this problem does not happen when doing BSE (yambo -b -o b -k sex) and I've done full BSE calculations up to 30^3. We can't do BSE on any finer grids, but we'd like to check the RPA spectrum (which should be much less expensive) on finer grids to estimate convergence of the excitons. I've attached the input file and one of the log files (as you can see it ends at about 2 minutes, but the calculation will go until the walltime of one hour). The report file is too large to attach in its entirety, but I attach a version with the k-points and occupations removed.

Thanks you!
You do not have the required permissions to view the files attached to this post.
Vatsal A. Jhalani
Postdoctoral Scholar | Department of Applied Physics
California Institute of Technology
vjhalani
 
Posts: 17
Joined: Mon Jan 22, 2018 8:23 pm

Re: RPA calculation hangs at dipoles

Postby Davide Sangalli » Wed Apr 17, 2019 7:42 am

Thank you Vatsal for reporting this issue.

Can you also attach the reports and logs of the succesfull BSE run ?
Recently I've seen similar problems and were due to missing or incomplete ns.kb_pwscf_fragment_* files.

Moreover, if the problem are the dipoles you can try one of the alternative approaches to compute them as described in our recent yambo paper:
https://iopscience.iop.org/article/10.1 ... 48X/ab15d0

I'd suggest the covariant dipoles.
Just add to the input file
Code: Select all
DipApproach="Covariant"


Best,
D.
Davide Sangalli, PhD
CNR-ISM, Division of Ultrafast Processes in Materials (FLASHit) and MaX Centre
http://www.ism.cnr.it/en/davide-sangalli-cv/
http://www.max-centre.eu/
User avatar
Davide Sangalli
 
Posts: 316
Joined: Tue May 29, 2012 4:49 pm
Location: Via Salaria Km 29.3, CP 10, 00016, Monterotondo Stazione, Italy

Re: RPA calculation hangs at dipoles

Postby vjhalani » Thu Apr 18, 2019 5:56 pm

Ciao Davide,

Sure, I've attached them (I had to compress them since the l-* file is too large). It seems that it gets stuck after kb_pp_pwscf_fragment_412, while the last one in the BSE is *_413. I checked the SAVE/ and ns.kb_pp_pwscf_fragment_413 is there.

I don't know if this could have an affect, but at one point while troubleshooting this, I switched to the official 4.3 release and tried an RPA calculation, and had the same issue. Is it possible that it affected the SAVE/ and there is some compatibility problem?

Grazie,
Vatsal
You do not have the required permissions to view the files attached to this post.
Vatsal A. Jhalani
Postdoctoral Scholar | Department of Applied Physics
California Institute of Technology
vjhalani
 
Posts: 17
Joined: Mon Jan 22, 2018 8:23 pm

Re: RPA calculation hangs at dipoles

Postby vjhalani » Wed Apr 24, 2019 12:28 am

Ciao,

As an added note, the covariant dipole approach resolves the issue!

Grazie,
Vatsal
Vatsal A. Jhalani
Postdoctoral Scholar | Department of Applied Physics
California Institute of Technology
vjhalani
 
Posts: 17
Joined: Mon Jan 22, 2018 8:23 pm


Return to Yambo version 4.0

Who is online

Users browsing this forum: No registered users and 1 guest