V (with respect to Li
cathode) which matches perfectly with the experimental report.
Looking forward to your further advice Sir.
with regards,
On Fri, 30 Dec 2022 at 21:01, xavier rocquefelte
wrote:
One more point the negative electrode is usually graphite and not
graphene.
30/12/2022 16:26, xavier rocquefelte wrote:
Dear Shamik
I certainly did not follow all the discussion but it seems to me that
your objective is reachable if you properly define your parameters.
You mention that the band structure is finally not the essential point
and you want to have a
IST
at:http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
--
--------
Professeur des Universités de Rennes 1
Institut des Sciences Chimiques de Rennes (ISCR)
Univ Rennes - CNRS - UMR6226, France
https://iscr.univ-rennes1.fr/xavier-rocquefelte
ac.at/index.html
--
Professeur des Universités de Rennes 1
Institut des Sciences Chimiques de Rennes (ISCR)
Univ Rennes - CNRS - UMR6226, France
https://iscr.univ-rennes1.fr/xavier-rocquefelte
___
Wien m
e Rennes 1
Institut des Sciences Chimiques de Rennes (ISCR)
Univ Rennes - CNRS - UMR6226, France
https://iscr.univ-rennes1.fr/xavier-rocquefelte
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.
ding various more or less empirical "tricks".
Eventually, you may try a meta-GGA ., but this is probably also
not very good.
Best regards
Peter Blaha
Am 2/11/22 um 08:19 schrieb xavier rocquefelte:
Dear Shamik,
To my point of view using the strategy (1) is not correct. I
un
Dear Atefe Marasi,
Infinite epsilon means that you extrapolate the epsilon value for the
zero of energy.
You must plot the real part of the dielectric function to properly
estimate this value.
You must be careful because if you have a band gap and a bad description
of the gap value ... the
nobody else has thought", Albert Szent-Gyorgi
www.numis.northwestern.edu <http://www.numis.northwestern.edu>
On Tue, Feb 16, 2021, 01:03 Xavier Rocquefelte
<mailto:xavier.rocquefe...@univ-rennes1.fr>> wrote:
Dear Colleagues,
We are doing EECE calculations
Dear Colleagues,
We are doing EECE calculations and we want to constrain the electronic
configuration of Eu.
It is possible to do it using DFT+U and -orbc option.
Do you know if we can do similar calculations using EECE? We look at
runsp and we found the use of .forcedmat and .forceorb.
Co
1. You cannot simulate paramagnetism using DFT. Indeed, you need to go
beyond in order to take into account the effect of temperature on the
dynamical disorder of the magnetic moments. Using DFT you will have
access to a static picture, which is OK if you want to simulate:
- a non-magnetic sta
e Peter knows.
From: Wien mailto:wien-boun...@zeus.theochem.tuwien.ac.at>> on behalf of
Xavier Rocquefelte mailto:xavier.rocquefe...@univ-rennes1.fr>>
Sent: Thursday, June 4, 2020 9:28 AM
To: wien@zeus.theochem.tuwien.ac.at
<mailto:wien@zeus.theochem.t
s with your cif file?
*From:* Wien on behalf of
Xavier Rocquefelte
*Sent:* Thursday, June 4, 2020 9:13 AM
*To:* wien@zeus.theochem.tuwien.ac.at
*Subject:* Re: [Wien] Structure optimization of Beta -Sn
Dear Fabien
This is strange... Perhaps a conversion problem. Yes I have used the
Indeed, I used the one from VESTA! Yesteray night I was working in
parallel ;)
Here is the file I previously sent with modification by hand ... the cif
format was not totally supported by cif2struct.
Regards
Xavier
) 2020 by FIZ Karlsruhe - Leibniz Institute for Information Infrastructure.
ruct, I get "unknown space
group name: I41/amds" and no struct file is generated.
Did cif2struct work for you?
*From:* Wien on behalf of
Xavier Rocquefelte
*Sent:* Wednesday, June 3, 2020 11:54
Dear Shamik
I will do as I do with my colleagues who are learning how to use WIEN2k ;)
Use the cif file I previously sent. If you directory is named Beta-Sn,
rename the cif file Beta-Sn.cif.
Then use cif2struct in this way: cif2struct Beta-Sn.cif. It will
generate Beta-Sn.struct
Then initi
/06/2020 à 23:54, Xavier Rocquefelte a écrit :
Dear Fabien,
I explained to Shamik that the structure he was using was not correct.
I also sent cif and struct file.
See below a proper case.struct file.
Shamik could you please send email only to the wienlist to avoid
multiple answers from the list
Dear Fabien,
I explained to Shamik that the structure he was using was not correct.
I also sent cif and struct file.
See below a proper case.struct file.
Shamik could you please send email only to the wienlist to avoid
multiple answers from the list and many people trying to help you
without
Your structure file seems to be not correct!
You should have only one inequivalent position and the space group is I
41/a m d.
Here is the structure:
http://som.web.cmu.edu/structures/S018-beta-Sn.html
Le 03/06/2020 à 20:40, shamik chakrabarti a écrit :
Dear Dr. Tran,
I
And another comment. Applying the magnetization along a high symmetry
direction does not mean that you will apply SOC in the direction for
which the effect will be maximum... It is why, people study the
magnetocrystalline anisotropy energy (MAE).
Best regards
Xavier
Le 20/05/2020 à 08:07, Pe
B 86,
075149 (2012)
Thanks
Wasim
On Tue, May 12, 2020 at 2:02 AM Xavier Rocquefelte
<mailto:xavier.rocquefe...@univ-rennes1.fr>> wrote:
Dear Wasim,
In VO2, vanadium is 4+ and thus we expect one electron in V(3d)
states, thus a gap based on d-d transitions.
MBJ alo
:
Dear All,
Sorry for my email. MBJ works for this material also.
Please see the reference if anybody interested: PHYSICAL REVIEW B 86,
075149 (2012)
Thanks
Wasim
On Tue, May 12, 2020 at 2:02 AM Xavier Rocquefelte
<mailto:xavier.rocquefe...@univ-rennes1.fr>> wrote:
Dear Wasim,
In VO2, vanadium is 4+ and thus we expect one electron in V(3d) states,
thus a gap based on d-d transitions.
MBJ alone is not expected to open the gap. You will need to put a
Hubbard term at least.
I recommend you to read the following article which deals with DFT+DMFT
and M1 p
LDA for U values I have suggested.
शुक्र, 10 अप्रैल 2020, 15:37 को Xavier Rocquefelte
<mailto:xavier.rocquefe...@univ-rennes1.fr>> ने लिखा:
One more comment, LDA+U (+SOC) will be better than everything in a
DFT level, because of iridium.
Le 10/04/2020 à 12:02, Xavier Rocq
One more comment, LDA+U (+SOC) will be better than everything in a DFT
level, because of iridium.
Le 10/04/2020 à 12:02, Xavier Rocquefelte a écrit :
I really recommend to look carefully at the litterature of this system.
Using mbj or hybrid will work but for what?
In addition, the states
atest/
http://hauleweb.rutgers.edu/tutorials/
On Fri, Apr 10, 2020 at 5:46 AM Xavier Rocquefelte
<mailto:xavier.rocquefe...@univ-rennes1.fr>> wrote:
For such a system DFT is not sufficient. You must use DMFT.
In addition, you should include spin-orbit coupling.
It
For such a system DFT is not sufficient. You must use DMFT.
In addition, you should include spin-orbit coupling.
It is thus a very difficult situation because DMFT+SOC is not trivial at
all.
Cheers
Xavier
Le 10/04/2020 à 11:39, Amit Chauhan ph17d008 a écrit :
Dear All:
I am working on a
Dear Shamik,
We can answer your question in many ways.
First of all, DFT is a theory of ground state properties. We thus do not
expect good simulations of excited states properties!
However, we are commonly using DFT to reproduce excited states
properties (optical, magnetic ...), while we mu
I am working on similar compounds but not with Mn.
I would recommend to use PBE+U or PBE0 on-site hybrid.
In the first case (PBE+U) you must define an appropriate value for Ueff
of Mn-3d states.
In the second case, use alpha = 0.25 for only M-3d states.
It should be a not to bad approximatio
e then 200 iterations never get convegence ???
>
> Le jeudi 21 novembre 2019 à 18:27:12 UTC+1, Xavier Rocquefelte
> <mailto:xavier.rocquefe...@univ-rennes1.fr>> a écrit :
>
>
> If you want our help we need to have details about your system.
>
> Then we can advice you
ybrid
functionel but the process did not converged after more than 200
iteration ??
can you Suggests me solutio
Le jeudi 21 novembre 2019 à 16:51:45 UTC+1, Xavier Rocquefelte
a écrit :
Hummm ... Here are many options ... choose the one you like ;)
1/ You consider the importance to have hybrid function
Hummm ... Here are many options ... choose the one you like ;)
1/ You consider the importance to have hybrid functional to properly
describe the atomic structure. If not, optimize using regular DFT and
then estimate the property of interest with hybrid functional using to
optimized structure i
I recommend you the following article:
https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.108.047201
Bulk gold is diamagnetic and it exhibits a tiny Pauli and Orbital
paramagnetic state.
Pauli paramagnetism is due to the non-zero density at the Fermi level
(metal) and the orbital para
Dear Anup,
The method you want to use is an "old" one which is now not anymore
used, except if you really have no other choice.
It is preferred to use an Hubbard (or an on-site hybrid) correction
which will allow to correct the DFT error concerning the treatment of
the 4f states of Sm.
I r
See the WIEN2k userguide at pages 177, 178 ...
One can constrain individual positions incase.inMor define linear
constrains for several po-sitions usingcase.constraint(thanks to
B.Yanchitsky (Kiev, y...@imag.kiev.ua); for detailssee comments in
the SRCtemplates/template.constraint file).
You can use a trick consisting to replace part of the magnetic atoms by
non-magnetic atoms (having similar radii).
For instance, Zinc for copper ... It works really nicely and a way to
check the consistency is to do the sum of the individual contributions
and see if you recover the total signa
Dear Pablo,
My understanding is that you have a mixed basis set, with plane wave
outside the atomic spheres which are augmented by atomic functions
inside the sphere. You have a good explanation in the David Singh book
but also in the WIEN2k userguide I believe.
Best Regards,
Xavier
Le 22
fudge -- the mixer
only works when the math is right!
On Thu, Nov 1, 2018 at 9:34 AM Xavier Rocquefelte
<mailto:xavier.rocquefe...@univ-rennes1.fr>> wrote:
I have exactly the same experience. It is really a matter of
cooking to
my point of view and strongly depends on t
I have exactly the same experience. It is really a matter of cooking to
my point of view and strongly depends on the system and what is
happening at the Fermi level.
In my case I am using DFT+U+SOC for such systems. The hubbard term
indeed allows to separate the 4f states. However, the SOC is
I have exactly the same experience. It is really a matter of cooking to
my point of view and strongly depends on the system and what is
happening at the Fermi level.
In my case I am using DFT+U+SOC for such systems. The hubbard term
indeed allows to separate the 4f states. However, the SOC is
Dear Xiangyan,
10 meV is more than very small!
However, accuracy and parameters to adjust (or tune) will really depend
on the material, the size of the unit cell, the nature of the bonding...
If you want a more accurate answer, we need a more accurate question ;)
Cheers
Xavier
Le 12/09/2
As you Laurence, I was thinking about the effect of shifting or not the
kmesh!
Peter, do you think it will lead to a better convergence?
Cheers
Xavier
Le 03/09/2018 à 13:59, Laurence Marks a écrit :
What you are doing "should" work -- I have done similar things myself.
I have also managed t
Dear Sylwia
I will give you a chemist answer, which appears to work properly in many
cases.
If you have 3d elements, the crystal field (CF) is usually larger than
SOC, thus I will do first orb, and then include SO.
If you have 4f elements, the SOC is usually larger than the CF, thus I
will
Dear Pablo
It really depends on the property you are considering:
1. In terms of chemical bond, the effect of SO on 3d and 4f elements
will be negligeable.
2. In terms of spectroscopy it can be significant if you have access to
d3/2 and d5/2 states, but also p1/2 and p3/2 ... and similar for
It seems you forget the option -c.
More generally speaking to help users we need more details.
Best Regards
Xavier
Le 08/03/2018 à 04:45, Lalit Mohan a écrit :
My version is 7.1
During execution of band structure task
at command x lapw1 -band
i am getting following message
Commandline: *x
Your cif file contains 3 settings explaining why cif2struct cannot work.
Here is the first setting of the cif file.
Cheers
Xavier
Le 06/03/2018 à 12:27, Sherif Yehia a écrit :
Dear Users
The attached .cif file data have been compiled from
the crystallographic data sheet for
"HoCo3 Crys
for your system.
Cheers
Xavier
Le 27/01/2018 à 16:01, chin Sabsu a écrit :
Sir, it is A2BX6, x is halogen.
Sent from Yahoo Mail on Android
<https://overview.mail.yahoo.com/mobile/?.src=Android>
On Sat, Jan 27, 2018 at 7:23 PM, Xavier Rocquefelte
bands with the structure having forces ~2-3
mRy/Bohr?
Thank you very much!
Chin S.
On Saturday 27 January 2018, 6:31:35 PM IST, Xavier Rocquefelte
wrote:
Dear Chin S.
Your system is metallic from the point of view of DFT which is not a
surprise.
In your last calculation, which i
Dear Chin S.
Your system is metallic from the point of view of DFT which is not a
surprise.
In your last calculation, which is the more accurate you simply catch
this fact ... In other words in your previous calculations you had 0.1
eV gap due to less kpoints and/or less accurate calculation
Dear Marcelo
Another option would be to do a chemical substitution of one element far
from the probed atom, for which you want to estimate EFG, in such a way
to compensate the charge.
This will allow you to use hybrid. Indeed, EFG is a local probe and thus
using a supercell with a reasonnabl
impact of such approximation ? Considering your last comments, it
seems not due to the execution of the orb program only in lapw1 and
not during the perturbative procedure.
Thank you by advance.
All the best,
William Lafargue-dit-Hauret
Le 16/01/2018 à 17:05, Xavier Rocquefelte a écrit :
I was
in case.inso)
Please check the presence of case.vorbud. It must not be there for EECE.
You could also test Laurence suggestion, running:
x lapw1 -up/dn -orb
x lapwso -up (no -orb !!!)
and see of it makes a difference.
On 01/16/2018 04:52 PM, Xavier Rocquefelte wrote:
Dear Peter
You are
/2018 04:52 PM, Xavier Rocquefelte wrote:
Dear Peter
You are totally correct. We are doing SO non-selfconsistent by using
a standard procedure for EECE calculations:
runsp_lapw -eece -p -ec 0.1 -NI
and then we estimate the MAE using this non-SCF procedure :
Increase EMAX in case.in1c
that for -eece -so, case.vorbud is NOT present (from
previous LDA+U).
Peter
On 01/16/2018 02:50 PM, Xavier Rocquefelte wrote:
Here is a document showing the results graphically.
https://filesender.renater.fr/?s=download&token=8ac3a214-edfa-4894-fa1f-27aba5a5522f
It really looks like
additional soft modes which only poorly converge.
On Tue, Jan 16, 2018 at 8:09 AM, Xavier Rocquefelte
wrote:
Dear Laurence
Here is the point. Our results show that in version 17 a problem occurs
related to the time-reversal, which appears only if we do GGA+U or EECE.
a) The calculations we
this should not matter.
b) Hack the runsp_lapw script so that -orb is applied to lapw1, not in lapwso
c) Run without the up/dn component (-noorbud)
On Tue, Jan 16, 2018 at 7:50 AM, Xavier Rocquefelte
wrote:
Here is a document showing the results graphically.
https://filesender.renater.fr/?s=dow
ture and not WIEN2k_17.
We are now comparing the two versions of the code.
Regards
Xavier
Le 16/01/2018 à 12:10, Xavier Rocquefelte a écrit :
Dear All
Finally the problem is not completely solved.
More precisely, when we are doing GGA+SO calculations and using a
correct kmesh (no temp
recent modifications in SRC_ORB
and SRC_LAPW2 related to the manipulation of case.vorbup, case.vorbdn
and case.vorbud files.
Surprisingly, the EECE+SO calculations in WIEN2k_16 are symmetric, while
not in WIEN2k_17.
Next soon ... I hope.
Xavier
Le 10/01/2018 à 15:10, Xavier Rocquefelte a
d to my question. It was
essential to identify such a mistake which has a huge impact on the
results.
Best wishes
Xavier
Le 10/01/2018 à 10:47, Xavier Rocquefelte a écrit :
Dear Lyudmila
The fact we have a small angle with axes is expected (also observed
experimentally). It is related to the
Dear Lyudmila
The fact we have a small angle with axes is expected (also observed
experimentally). It is related to the monoclinic symmetry of the system
which permits it. However, you gave me an idea that I will test now and
comment soon ;)
Cheers
Xavier
Le 10/01/2018 à 10:40, Lyudmila a
th my best regards
Jaroslav
On 09/01/18 09:44, Xavier Rocquefelte wrote:
Dear Colleagues
I recently obtained a surprising result concerning the calculation of
the magnetocrystalline anisotropy energy (MAE) of SeCuO3.
This compound has a monoclinic symmetry (SG. P21/n) and is known to
be ant
Dear Colleagues
I recently obtained a surprising result concerning the calculation of
the magnetocrystalline anisotropy energy (MAE) of SeCuO3.
This compound has a monoclinic symmetry (SG. P21/n) and is known to be
antiferromagnetically ordered at low temperature.
Here I provide the results
Indeed, it seems that there is an error in the read_input.f file of SRC_hf.
If SCREENED if TRUE everything is fine, but not when SCREENED is FALSE,
because in this case the program does not read SMU.
Here I did a simple modification (in red) of read_input.f :
...
read(4,*) screened
o dig
through the code to see if the vector values being output are time
forward or backwards."
What you observe is certainly due to this aspect because point group of
the iron atoms in this cell is 4mm.
Cheers
Xavier
Le 03/01/2018 à 15:41, Xavier Rocquefelte a écrit :
Dear Gerhard
One cla
Dear Gerhard
One clarification is needed I think. The discussion was about applying
an external ELECTRIC field (not a magnetic one).
Thus one part of your answer concerns something else which is also
interesting :) Indeed, my PhD student has written a modification of
WIEN2k to take into acco
‘outward’ for the other layer, is
not correct. Yet I don’t see why.
Thanks!
Stefaan
*Van:*Wien [mailto:wien-boun...@zeus.theochem.tuwien.ac.at] *Namens
*Xavier Rocquefelte
*Verzonden:* dinsdag 2 januari 2018 15:38
*Aan:* wien@zeus.theochem.tuwien.ac.at
*Onderwerp:* Re: [Wien] zigzag potential
To be honest I also feel that something is missing in my last arguments.
What is the electronic configuration of Fe at the surface? The orbital
occupancy could play a role in the understanding of the present
observation.
Le 02/01/2018 à 15:37, Xavier Rocquefelte a écrit :
Dear Stefaan
As
Dear Stefaan
As always it is very nice to read your posts :)
I will only react on your "Thought 3". What will happen if you do the
same calculation along 00-1? To my point of view, you will obtain the
same result. Indeed, the magnetic anisotropy (MAE) of bulk-Fe must be
symmetric. Here you br
Dear Chin
Your question is a little aside from the WIENLIST. Indeed, it is
thermodynamics.
In equation 1, you evaluate the enthalpie of formation of the compound
ABO3 from the reactants AO and BO2 using internal energies. Indeed, you
neglect the volumetric effect (P deltaV). From your DFT ca
Dear Bingrui Peng
What are the eigenvalues and how many bands do you have in the valence
states.
When you did "x lapw1 -band" what was the energy range of your calculation.
Look at the case.output1 file, it contains the eigenvalues for each
k-points. You will then see the number of bands and
Dear Timothy
The main problem here is that in such intermetallic compound you must
have both itinerant and localized magnetisms which is not so trivial to
treat from DFT and adding an Hubbard correction will not help. To my
point of view the problem can be solved by looking at the experimental
Dear Timothy
The main problem here is that in such intermetallic compound you must
have both itinerant and localized magnetisms which is not so trivial to
treat from DFT and adding an Hubbard correction will not help. To my
point of view the problem can be solved by looking at the experimental
Dear Shamik
The userid and password are generated when you start for the first time
the w2web interface by typing in the unix terminal:
w2web
If it has been already generated, you have then created a directory
named .w2web which contains configuration files.
One easy solution is to remove
One of the best approach is to do a SCF calculation for a small number
of k-points and then increase gradually this number.
Then you plot the total energy as a function of the number of k-points
and you will choose the smallest number giving rise to a good convergency.
Ideally you should do s
It seems that you did a mistake in the generation of the case.struct
file. If you have the structure in VESTA you can export in CIF format
and then use cif2struct.
Best Regards
Xavier
Le 01/02/2017 à 08:21, Rajneesh Chaurasiya a écrit :
Dear Wien2k user,
I want to compute the electronic pr
7;t think of it immediately.
On Fri, Jan 20, 2017 at 2:03 PM, Xavier Rocquefelte
mailto:xavier.rocquefe...@univ-rennes1.fr>> wrote:
Dear Colleagues
I did recently a calculation which has been published long time ago
using a old WIEN2k version (in 2008).
It correspon
n many cases
these have small eigenvalues in the mixing Jacobian which are removed
when symmetry is imposed. All one can do is use MSEC3 or some of the
additional flags (see the mixer README) such as "SLOW".
Fifth...probably exists, but I can't think of it immed
Dear Colleagues
I did recently a calculation which has been published long time ago
using a old WIEN2k version (in 2008).
It corresponds to a spin-polarized calculation for the compound CuO. The
symmetry is removed and the idea is to estimate the total energies for
different magnetic orders
niversité Sétif1 (Algeria)/
/http://laboratoires.univ-setif.dz/L.E.S.I.M.S//
Le Vendredi 13 janvier 2017 20h04, Xavier Rocquefelte
a écrit :
I imagine that you have done a spin-polarized calculations. If so, you
must sum the spin-up and spin-down contributions when you are doing
the epsilon2 spe
I imagine that you have done a spin-polarized calculations. If so, you
must sum the spin-up and spin-down contributions when you are doing the
epsilon2 spectrum without SO.
Cheers
Xavier
Le 13/01/2017 à 19:06, Nacir GUECHI a écrit :
Dear Professor Peter Blaha and Wien2k users.
I calculated
Dear Bhamu
The Fermi level when expressed in Rydberg is given with respect to the
printed value in the scf file.
In contrast, by convention the Fermi level is fixed at 0 when given in eV.
In your case, if you have the following data :FER = +0.134 Ry, it means
that you must supply this value
Dear Bhamu
In the present case, you are dealing with Ag+ ions and thus a d10
electronic configuration for silver.
In such a case applying a Hubbard correction will mainly lead to correct
the position of the Ag(4d) states which are all occupied and below the
O(2p) band. It should not affect s
Just to add one more point to this funny discussion, the term
"paramagnetic" is sometimes used in the DFT litterature in an improper way.
It could clearly lead to misunderstanding for researchers who do not
know so much on how magnetic properties could evolve with temperature
and applied magne
Dear Bhamu
First of all, your band gap to be accurate need many kpoints when taken
from the SCF file.
The best is to plot the band structure. Then you will have the
fundamental band gap, not the optical band gap (usually larger).
In your case, you must also include the spin-orbit coupling t
Dear Mengxi Wu
If I clearly understand you, you are considering the GAMMA point where
the bands are degenerate.
At the GAMMA point these bands have the same symmetry to my point of
view, explaining why the dipole matrix is non-zero. Thus your result is
expected to my point of view.
Best Re
Could you clarify your question. What do you mean by "each range"?
- Energy range
- Density value
Le 13/11/2016 à 16:54, boudiaf khadidja a écrit :
Dear Wien2k user
I use VIETA to plot the charge density ;But I need to know a way to
plot or get the values of each range for the charge densit
Dear R. Chouhan
Just to complete the very nice answer of Peter. I have used the force
theorem as explained by Peter using GGA+U to estimate MCA and in the
cases I have considered it works amazingly nicely.
As Peter said, I was using P1 symmetry and checking carefully the
convergency before i
Dear Colleague
The option "-o" does not exist.
It seems that you wanted to put "-so" ... Is it the case ?
Cheers
Xavier
Le 02/11/2016 à 14:24, Abderrahmane Reggad a écrit :
Thanks Dr Assmann for your quick reply
I repetad the calculation and I same problem and I want to note that I
modifi
Dear Colleague
The notation is with three indices.
It will be along the directions defined in your case.struct file.
In other words, if you put 0 0 1, it will be along the c-axis of your
case.struct file.
Best Regards
Xavier
Le 09/10/2016 à 12:52, Abderrahmane Reggad a écrit :
Dear Wien
i manually created the case.indm
file, the scf runs only a single cycle . Thereafter showing a Stop message
Thanks
On Tue, Aug 2, 2016 at 1:45 PM, Xavier Rocquefelte
<mailto:xavier.rocquefe...@univ-rennes1.fr>> wrote:
Dear Shakeel
You must create two files when initializin
Dear Shakeel
You must create two files when initializing a GGA+U (or LDA+U)
calculation, i.e. case.inorb and case.indm (or case.indmc if no
inversion center in the space group).
Normally when you start such a calculation from the w2web interface, the
interface open a template file to allow y
Indeed I never use the scf file when I restart a calculation.
Have a nice week
Xavier
Le 25/04/2016 08:38, Stefaan Cottenier a écrit :
You must use the command
restore_lapw -d directory
Apparently, your computer works better than ours. I also don't see the
case.scf restored unless the scri
Dear Dr. Bhamu
you forget to put "-d" which explain why you did not restore from the
files contained in the directory.
You must use the command
restore_lapw -d directory
More details, here:
http://100.36.165.205/usersguide/5_Shell_scripts.html#SECTION05222000
Regards
Xavier
Le 2
Dear Brik
What you found is usual and is simply the signature of a high absorption in
this high energy region.
Indeed, a simple way to understand the negative value of eps1 is to look at the
following relation:
eps1 = n^2 - k^2
with n the real part and k the imaginary part of the complex refr
Dear Fhokrul
Your structure looks strange to me. The angle is 60° and not 120° and it seems
that you do not take benefit of the symetry.
Here is a cif file of graphene based on the graphite structure in which I have
simply increase the c parameter.
The structure you are using may explain why yo
About the Néel temperature, I know two main strategies:
- If you know an analytical formula which relates TN with J values (J: magnetic
exchange parameters of your system) then you can estimate TN from DFT
calculations. It will require the calculation of all the J's values of your
system. Such
You must also read this document produced by Stefaan Cottenier.
http://www.wien2k.at/reg_user/textbooks/DFT_and_LAPW-2_cottenier.pdf
In particular p43 and the following pages.
Regards
Xavier
sikander Azam a écrit :
> Dear all
> When we plot the band structure so what the brillion zone symmet
Dear Azam
Many excellent books could help you to understand how a band structure works.
I recommend you to look at the following web sites:
http://www.chem.uci.edu/~lawm/Hoffmann.pdf
http://www.chemeddl.org/alfresco/service/api/node/content/workspace/SpacesStore/addbf469-059f-44ff-beb5-3417c58e
Dear Rishi
Indeed, WIEN2K provides numbers with 4 digits however it does not mean that the
value is accurate up to 4 digits.
The accuracy will depend on:
- the convergence criteria you are using.
- the RMT you are using. Indeed the magnetic moment for each atom is estimated
inside the muffin-ti
Dear Farshad,
As mentionned by Pascal you must use a supercell.
For instance, if you use a small supercell with 2a x 2b x c, you will
have 8 positions for the 2a site. In such a sitiuation you will be able
to do the following occupation:
- 1Li and 7 Fe, i.e. occupations of 0.125/0.875
- 2Li an
for your quick and detailed reply. I am
> calculating the valance optical properties, where the hole is
> positioned in the valance states; JDOS, Energy loss function etc.
>
> Thanks
>
> - Original Message -
> From: "Xavier Rocquefelte"
> To: "A Mail
1 - 100 of 105 matches
Mail list logo
