Search found 34 matches

by fhssn1
Sun Sep 15, 2019 6:55 am
Forum: Ground state
Topic: relaxation of supercell different from primitive cell
Replies: 0
Views: 2775

relaxation of supercell different from primitive cell

I have carried out a highly converged primitive cell calculation (optcell 0, ionmov 2, followed by optcell 2) for a wurtzite/hex crystal. The input file is as follows (sorry for using mixed units, Bohr for acell, and Angstrom for xangst): autoparal 1 ecut 25 pawecutdg 25 optcell 2 dilatmx 1.05 ecuts...
by fhssn1
Thu Jul 04, 2019 9:23 pm
Forum: Ground state
Topic: State of hybrid functionals? (B3LYP and HSE06)
Replies: 1
Views: 2155

Re: State of hybrid functionals? (B3LYP and HSE06)

On further inspection of the input files shared in 'Hybrids' page on website, featuring C diamond, N2, and CH calculations, I have a couple more questions: In N2 and CH calcs specifically, it's a 3 step process (1) ionmov with GGA-PBE (2) ionmov with PBE0 (3) ionmov with HSE06. Why was PBE0 used? A ...
by fhssn1
Thu Jul 04, 2019 5:22 am
Forum: Configuring and compiling ABINIT
Topic: autoparal file writing gives MPI ABORT (cluster install)
Replies: 0
Views: 2508

autoparal file writing gives MPI ABORT (cluster install)

I compiled abinit-8.10.2 with mpi and with mpi-io on one of the XSEDE clusters (comet to be specific). I used spack for compilation, and the install was in my local directory (no system wide install). The package appears to be working fine in serial mode. However, when I have a multi-dataset run (nd...
by fhssn1
Thu Jul 04, 2019 12:48 am
Forum: Ground state
Topic: State of hybrid functionals? (B3LYP and HSE06)
Replies: 1
Views: 2155

State of hybrid functionals? (B3LYP and HSE06)

Can abinit 8.10.2 or 8.10.3 be used for hybrid functional calculations? specifically B3LYP or HSE06. I need to do a dimer calculation for two different atom types. The Hybrid page on abinit website says I can do it in two steps: (1) do a ground-state calc in GGA, (2) do a hybrid calc. For step 1, do...
by fhssn1
Tue May 01, 2018 8:33 am
Forum: Ground state
Topic: How to setup a vacancy in a charged state?
Replies: 1
Views: 2459

How to setup a vacancy in a charged state?

I need to create a vacany in a supercell by removing one of the atoms. But that vacancy needs to be in a charged state. I know I can specify 'charge value' in the input file but that is for the whole cell in general, not specific to the vacancy. Is there a way to specify a charge for a vacancy (or e...
by fhssn1
Sun Apr 22, 2018 12:54 am
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Re: Unable to reproduce magnetic moment for YbN

Eric, I have good news. I conducted a '+U' run with the U value reported in the paper and, with a few other changes in the input file (increase of ecut from 20 Ha to 25 Ha and pawecutdg from 40 Ha to 60 Ha, setting N spinat to zero), I was able to get a magnetic moment of 0.183 uB, which is between ...
by fhssn1
Sat Apr 21, 2018 11:03 pm
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Re: Unable to reproduce magnetic moment for YbN

Eric, I just realized I missed a line in PAW paper that mentioned that the calculations were carried out using '+U' formalism (I originally thought +U was used only to generate U fitting values but not for non-AE result validation against AE). I hope this turns out to be a key point and fixes my res...
by fhssn1
Thu Apr 19, 2018 9:42 am
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Re: Unable to reproduce magnetic moment for YbN

Eric, thanks for your two responses. Between FeO, MnO, higher ecut, and turning N spinat off, I think I have a handful of direction to try. I will work on it and see if I get more updates.
by fhssn1
Thu Apr 19, 2018 9:36 am
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Re: Unable to reproduce magnetic moment for YbN

I increased ngkpt to 8 8 8, nstep to 200, ecut is 20.0, pawecutdg 40.0. And used very high values for z component of Yb's spinat (from 10 to 32). And following are the results (spinat vs magnetic-moments): 1.0000000000E+01 1.75317468E+00 1.2000000000E+01 2.33364378E-06 1.4000000000E+01 5.18211156E-0...
by fhssn1
Thu Apr 19, 2018 8:29 am
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Re: Unable to reproduce magnetic moment for YbN

To elaborate further, I think the randomness may have been due to using too little k-points (ngkpt 2 2 2). I ran it again with ngktp 6 6 6 and the magnetization goes even lower. The values so far are as follows: Magnetization (Bohr magneton)= 5.29398019E-08 Magnetization (Bohr magneton)= 2.28816879E...
by fhssn1
Thu Apr 19, 2018 7:25 am
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Re: Unable to reproduce magnetic moment for YbN

I have done ndtset studies for both unconstrainted (spinat ranging from [0,0,1.0] going up to [0,0,12.0]) as well as constrained (magconon_lambda at 1e-5, 1e4, 1e-3, 1e-2, etc, at fixed spinat [0,0,4]). I'm unable to reach any conclusions: - I am using spinat for Yb and N like this: [0,0,value], [0,...
by fhssn1
Wed Apr 18, 2018 10:36 pm
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Re: Unable to reproduce magnetic moment for YbN

Thanks so much Eric for the quick responses today. I will take your advice into account.
by fhssn1
Wed Apr 18, 2018 9:52 pm
Forum: Ground state
Topic: DoS resolution in conduction bands, and prtdos
Replies: 8
Views: 7053

Re: DoS resolution in conduction bands, and prtdos

Thanks Eric. I will take your advise to carry out further computations. In the meantime, would you be kind enough to answer my other questions in Ground State section of this forum? (especially the YbN magnetic moment one) I haven't received any answers in the past almost two weeks and I'm under a b...
by fhssn1
Tue Apr 17, 2018 10:57 am
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Re: Unable to reproduce magnetic moment for YbN

Ok, after a lot of trial and error, I've been able to do a calculation that can at least give 'Total magnetization (exact up - dn)' that is not in the 1e-5 range, but higher. This is as a result of discovering two keywords magconon and magcon_lambda, which have absolutely no mention in the tutorial ...
by fhssn1
Thu Apr 12, 2018 6:30 am
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Re: Unable to reproduce magnetic moment for YbN

Just to add, this is the VLab PAW RE datasets paper, which is also recommended on this forum by the ABINIT team. ( http://www.vlab.msi.umn.edu/resources/repaw/index.shtml ) So is it possible any members of the ABINIT team have already run and verified at least part of that paper, and if yes, are the...
by fhssn1
Thu Apr 12, 2018 6:27 am
Forum: Ground state
Topic: Unable to reproduce magnetic moment for YbN
Replies: 13
Views: 11042

Unable to reproduce magnetic moment for YbN

I'm trying to reproduce the results of one REN (YbN) out of the many REN (rare-earth nitrides) that the following paper has published: Topsakal, Wentzcovich - Accurate projected augmented wave (PAW) datasets for rare-earth elements (RE = La–Lu) I'm using their PAW dataset, with GGA-PBE for Yb, while...
by fhssn1
Wed Apr 11, 2018 7:38 pm
Forum: Ground state
Topic: spin polarized convergence studies?
Replies: 1
Views: 1944

spin polarized convergence studies?

If we need a spin polarized system, do we have to carry out convergence-studies/relaxation in spin polarized configuration as well? or does it suffice to do convergence/relaxation in non-spin polarized and only use spin polarized for main calculation and post-processing calculations (e.g., BS, DoS, ...
by fhssn1
Tue Apr 10, 2018 6:14 pm
Forum: Ground state
Topic: DoS resolution in conduction bands, and prtdos
Replies: 8
Views: 7053

Re: DoS resolution in conduction bands, and prtdos

Well I generated DoS with dosdeltae of 0.001 Ha and it was a bit of an improvement but not too much. I ran with 0.004 Ha and the plot looks like this: noisydos2.png The noise is much reduced but the result is not what I'm looking for. It looks like there are unnecessary oscillations in various regio...
by fhssn1
Tue Apr 10, 2018 8:10 am
Forum: Ground state
Topic: DoS resolution in conduction bands, and prtdos
Replies: 8
Views: 7053

Re: DoS resolution in conduction bands, and prtdos

I made a 2x2x2 supercell out of my system and computed the DoS, and unfortunately, I now have the opposite problem, my DoS is too high-resolution and noisy (I used prtdos 2). noisydos.png High-resolution is okay, but I don't want the noise. In general I know a supercell DoS becomes more noisy becaus...
by fhssn1
Tue Apr 10, 2018 7:47 am
Forum: Ground state
Topic: toldfe for a primitive cell vs a supercell?
Replies: 1
Views: 1668

toldfe for a primitive cell vs a supercell?

Do we relax the toldfe value (for SCF run) when we use a supercell compared to when we use a primitive cell? E.g., if we are using toldfe of 1e-12 for a 4 atom primitive cell, and now we need to do a 2x2x2 supercell calculation (32 atoms). The system is 8 times bigger, so we can use 8e-12 instead an...
by fhssn1
Tue Apr 10, 2018 2:20 am
Forum: Ground state
Topic: structural relaxation for supercell with impurity
Replies: 9
Views: 5396

Re: structural relaxation for supercell with impurity

I was told it may have something to do with a not well-tested PAW dataset of the impurity (because I have difficulties only when I add the impurity). My impurity is one of rare earth elements, and I recently came across VLab PAW datasets (and the paper by Topsakal) of rare earth elements, and I'm g...
by fhssn1
Tue Mar 27, 2018 7:55 pm
Forum: Ground state
Topic: DoS resolution in conduction bands, and prtdos
Replies: 8
Views: 7053

Re: DoS resolution in conduction bands, and prtdos

Thanks so much Eric. I guess I have no more questions on this topic at the moment and I will take your advice to further my calculations.
by fhssn1
Tue Mar 27, 2018 7:53 pm
Forum: Ground state
Topic: structural relaxation for supercell with impurity
Replies: 9
Views: 5396

Re: structural relaxation for supercell with impurity

PS: is your SCF well converged for each relaxation step? Ok I guess my understanding of optcell 0 was incorrect (especially when the log file calls it 'max grad (forces/stress)' not just 'max (forces/stress)'. I guess if it's a grad (gradient) then it makes more sense to me. As for your question, t...
by fhssn1
Tue Mar 27, 2018 12:09 am
Forum: Ground state
Topic: structural relaxation for supercell with impurity
Replies: 9
Views: 5396

Re: structural relaxation for supercell with impurity

The advise to first relax the internal coordinates and then relax the cell is because in many systems doing so will help the algo to converge the relaxation by reducing the number of degrees of freedom. Starting the cell relaxation after will be less far from the equilibrium than from scratch, thou...
by fhssn1
Mon Mar 26, 2018 9:08 pm
Forum: Ground state
Topic: DoS resolution in conduction bands, and prtdos
Replies: 8
Views: 7053

Re: DoS resolution in conduction bands, and prtdos

Dear fhssn1, Your plot looks good to me, what do think is wrong actually? If you want the DOS with more unoccupied states, you have to generate a density file (or WFK) that contains enough bands, otherwise it will indeed not show them (the complain with 16 is because 16 band is not enough band to p...