[Wien] Questions about imposing external magnetic field on no-magnetic system
Fecher, Gerhard
fecher at uni-mainz.de
Wed Jul 19 18:01:57 CEST 2017
An additional thought about the B field effect.
The answer to what a B field is doing to the electronic structure might have also
some more subtle aspects.
One aspect is -- as usual -- the role of symmetry. Similar to the case of spin-orbit
interaction in a ferromagnet, the application of a magnetic field will indeed
change the symmetry (see Koster et al), even though the bare calculation scheme
of Wien2k does not make use of it directly (means: there is no init_bfield).
To have a most easy case, let's assume sodium and le's apply a magnetic field.
What we expect is to find a Zeeman type splitting of the core level, that are 2s and 2p.
If checking the irreps (here for Oh) and energies after a spinpolarized
calculation using a field of 100 T, the result is
(shortened version, longer versions are attached):
from irrep -up
bnd ndg eigval
1 1 -3.887865 =G1+ ==> 2s (a1g)
2 3 -1.822866 =G4- ==> 2p (t1u)
5 1 -0.233852 =G1+ ==> lowest band
from irrep -dn
1 1 -3.886985 =G1+ ==> 2s (a1g)
2 3 -1.821986 =G4- ==> 2p (t1u)
5 1 -0.232453 =G1+ ==> lowest band
that is one has a splitting of the 2s and the 2p states, however the 2p is split
only into two level but we expect 6 !
Performing the same calculation with spin-orbit action respected will reduce the symmetry
(probaly one may fail with running irrep for that situation with
"X not equal for all elements in the class", but more about that in another task).
As a result of the lowered symmetry one finds the irreps and energies for C4h
(I hope the signs at the mj of |l,j,mj> are without typos):
bnd ndg eigval
1 1 -3.887859 =G5+ ==> 2s |0, 1/2, +1/2>
2 1 -3.886980 =G6+ ==> 2s |0, 1/2, -1/2>
3 1 -1.831327 =G6- ==> 2p |1, 1/2, -1/2>
4 1 -1.831034 =G5- ==> 2p |1, 1/2, +1/2>
5 1 -1.818491 =G8- ==> 2p |1, 3/2, +3/2>
6 1 -1.818185 =G5- ==> 2p |1, 3/2, +1/2>
7 1 -1.817893 =G6- ==> 2p |1, 3/2, -1/2>
8 1 -1.817611 =G7- ==> 2p |1, 3/2, -3/2>
9 1 -0.233828 =G5+ ==> lowest s band is magnetically split
10 1 -0.232482 =G6+
that is the result reveals the Zeemann splitting as expected.
Note: Have a look not just on the energies but also on the irreps and spinor wave functions !
If one has more complicated atoms or compounds then pronounced splitting effects
may also appear in the band structure and should not be neglected.
Overall I would conclude that the application of a B field makes only sense if
one uses it together with spin-orbit interaction because otherwise the
calculation will be for a "wrong" symmetry.
Further, I think one should clearly distinguish between the "symmetry of the atomic
positions" and the "symmetry of the field or the (spinor) wave functions".
A certain mirror operation may keep the positions intact but may change the spin
or reverse the direction of magnetisation. The overall symmetry has to keep
all, atomic positions AND field direction.
This concerns also electric fields, indeed.
Ciao
Gerhard
DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:
"I think the problem, to be quite honest with you,
is that you have never actually known what the question is."
====================================
Dr. Gerhard H. Fecher
Institut of Inorganic and Analytical Chemistry
Johannes Gutenberg - University
55099 Mainz
and
Max Planck Institute for Chemical Physics of Solids
01187 Dresden
________________________________________
Von: Wien [wien-bounces at zeus.theochem.tuwien.ac.at] im Auftrag von Peter Blaha [pblaha at theochem.tuwien.ac.at]
Gesendet: Sonntag, 16. Juli 2017 21:25
An: wien at zeus.theochem.tuwien.ac.at
Betreff: Re: [Wien] Questions about imposing external magnetic field on no-magnetic system
Once more: A magnetic field influences the spin and orbital degrees of
freedom.
The spin effects can approximately be taken care of as described in the
UG for NMR in metals. It leads to a trivial (or non-trivial if there is
screening) Zeeman splitting. Since even a large field of 100 T is only 1
mRy splitting, you get in first approximation 2 rigid band structures
shifted by that value. In semiconductors, that shift is probably
everything, however, in metals scf effects may affect this a little bit.
You may get estimates of the induced magnetic moments, or the spin
suszeptibility.
The magnetic field induces also an orbital current. This current is
calculated in the NMR module (you can even plot it) and the orbital
suszeptibility as well as the induced magnetic field is also calculated,
however, only at the position of the nuclei, not in the whole crystal.
In addition, as I mentioned before, this magnetic field breaks
translational symmetry and without that, the concept of "bandstructure"
is in principle not valid anymore.
The "magnetic field effect" in case.inorb as described by Pavel Novak is
a central field (single free atom) approximation and can be used to get
the induced orbital magnetic moment for atoms with localized 4f (or
maybe 3d) electrons. It cannot be used as a first principle method to
obtain all magnetic field effects in every kind of solid.
I'm not an expert in in this kind of physics and thus cannot say much
more about it, but eg. the anomalous Hall effect can be obtained from
the off-diagonal epsilon in spin-orbit optics calculations (Jan Kunes)
and the de Haas-van Alphen measurements give you effectively the ground
state cross sections of the Fermi surface (one does not even need a
magnetic field to calculate this, although experimentally you may even
need very large fields to observe these oszillations).
Am 16.07.2017 um 10:40 schrieb Karel Vyborny:
> As for NMR in not-too-strong B-field, it may indeed not be necessary to
> consider the orbital effects of B. I am not an expert in NMR. What I had
> in mind is related to my quantum Hall effect background and could be
> explained e.g. with bulk GaAs as an example.
>
> For B=0, the conduction band of this semiconductor is (to a good
> approximation) parabolic and centered around Gamma. The corresponding
> DOS is ~\sqrt{E-Eg} (when Ef is at the top of the valence band and Eg is
> the gap). With magnetic field switched on, two things happen. First,
> spin up and spin down bands get split by the Zeeman energy. The less
> trivial effect is the Landau quantisation (which is what I mean by
> "orbital effects"). In principle for any finite B, the smooth DOS breaks
> up into a comb of van Hove singularities ~1/\sqrt{E-Eg-Eorb*(n+1/2)},
> n=0,1,2 (Eorb is the cyclotron energy). In reality (for a real sample),
> any disorder will smear out this structure and DOS~\sqrt{E-Eg} is
> recovered unless B is really strong. Condition for "strong" is something
> like Eorb>>Gamma (disoreder broadening).
>
> The fact that Shubnikov-de Haas and de Haas-van Alphen oscillations can
> be observed in some bulk solids shows that "strong B" is indeed
> achievable. Those 1728 T mentioned below would certainly be strong
> enough for many
> real systems. However, fields of max. several tesla would not - unless
> we deal with a very clean system (like 2DEGs needed for QHE).
> Nevertheless, I'd think twice before showing any "band structure with B
> switched on" as calculated by WIEN.
>
> KV
>
>
> --- x ---
> dr. Karel Vyborny
> Fyzikalni ustav AV CR, v.v.i.
> Cukrovarnicka 10
> Praha 6, CZ-16253
> tel: +420220318459
>
>
> On Sat, 15 Jul 2017, Gavin Abo wrote:
>
>>
>> I looked at the Landau quantization Wikipedia entry [1]. However, it was
>> not clear to me whether this was needed to describe a system with moving
>> spin (e.g., oscillating spins).
>>
>> If so, I think the answer to your question it that your not missing
>> anything
>> and WIEN2k does not have an external magnetic field implementation for
>> Landau quantization.
>>
>> In Chapter 10 Landau Quantization on page 182 of the book titled "Quantum
>> Hall Effects: Recent Theoretical and Experimental Developments" by
>> Zyun F.
>> Ezawa, it mentions that spinless theory is frequently considered when the
>> spin degree of freedom can be ignored, such that a spin frozen system
>> becomes a good approximation under the condition that the Zeeman
>> energy is
>> large.
>>
>> Previously, I didn't understand Dr. Novak's reference to the frozen spin
>> method [2], but it seems now that might be why he mentioned it.
>>
>> The NMR slides [3,4] do show B_ext in the H_NMR equation, but I don't
>> see it
>> described in which input file it is to be included (or if just part of a
>> result in an output file). There is the external magnetic field value
>> that
>> can be entered in case.inorb [5]. Perhaps, the NMR program also uses
>> that
>> too.
>>
>> Of note, it was estimated before that a Bext value of a least 1728 T
>> may be
>> needed to see any noticeable effect in the plots (if the default
>> autoscale-like settings are used) [6].
>>
>> [1] https://en.wikipedia.org/wiki/Landau_quantization
>> [2]
>> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg01508.html
>> [3] http://susi.theochem.tuwien.ac.at/events/ws2015/rolask_nmr.pdf
>> [4]http://susi.theochem.tuwien.ac.at/reg_user/textbooks/WIEN2k_lecture-notes_2
>>
>> 013/nmr-chemical-shift.pdf
>> [5]
>> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg12904.html
>> [6]
>> https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg11093.html
>>
>>
>> On 7/15/2017 4:56 AM, Karel Vyborny wrote:
>> Interesting, I didn't know that WIEN2k can figure out what "band
>> structure with B>0" is... I thought there ought to be some
>> Landau quantisation which is hard to do except for idealised
>> systems. Am I missing something here?
>>
>> KV
>>
>>
>> --- x ---
>> dr. Karel Vyborny
>> Fyzikalni ustav AV CR, v.v.i.
>> Cukrovarnicka 10
>> Praha 6, CZ-16253
>> tel: +420220318459
>>
>>
>> On Sat, 15 Jul 2017, Peng Bingrui wrote:
>>
>> Dear professor Blaha
>>
>> Thank you very much for your suggestions. However,
>> I'm still kind of
>> confused, because my purpose is to see the change of
>> band structure under
>> external magnetic field, and l'm wondering whether
>> NMR calculation can do
>> this ? I'm sorry for my limited knowledge as an
>> undergraduate student.
>>
>> Sincerely yours,
>> Bingrui Peng
>> from the Department of Physics, Nanjing University,
>> China
>>
>> ___________________________________________________________________________
>>
>> _
>> From: Wien <wien-bounces at zeus.theochem.tuwien.ac.at>
>> on behalf of pieper
>> <pieper at ifp.tuwien.ac.at>
>> Sent: Wednesday, July 12, 2017 1:15:41 AM
>> To: A Mailing list for WIEN2k users
>> Subject: Re: [Wien] Questions about imposing
>> external magnetic field on
>> no-magnetic system
>> In case no one has answered this up to now:
>>
>> ad 1) The procedure itself is ok. You might want
>> switch on SO first and
>> converge that without the orbital potential to
>> establish a zero-field
>> base line. Remember to put in LARGE fields - your
>> off-the-shelf lab
>> field of 10 T will not show up at any energy
>> precision you can achieve.
>> Estimate the energy of 1 mu_B in 10 T field in Ry
>> units to see that.
>>
>> Note that your not-so-recent version of Wien2k is
>> not the best for the
>> task. The latest version is 17.1. With 16.1 came the
>> NMR package which
>> should be much better suited to calculate the
>> effects of a magnetic
>> field.
>>
>> ad 2) If you apply a magnetic field experimentally
>> in the lab you do it
>> at all atoms. I suppose you want to model that
>> situation. imho it makes
>> little sense to exempt one or two of your atoms from
>> the field.
>>
>> Good luck
>>
>> ---
>> Dr. Martin Pieper
>> Karl-Franzens University
>> Institute of Physics
>> Universitätsplatz 5
>> A-8010 Graz
>> Austria
>> Tel.: +43-(0)316-380-8564
>>
>>
>> Am 10.07.2017 12:20, schrieb Peng Bingrui:
>> > Dear professor Blaha and WIEN2K users
>> >
>> > I'm running WIEN2K of 14 version on Linux system.
>> I'm going to impose
>> > external magnetic field on LaPtBi, a no-magnetic
>> material. The
>> > procedure that I'm going to use is :
>> >
>> > 1、Do a no-SO calculation : runsp_c_lapw.
>> >
>> > 2、Do a SO calculation : runsp_c_lapw -so -orb,
>> while including
>> > external magnetic field as orbital potential in
>> the same time.
>> >
>> > My questions are:
>> >
>> > 1、Whether this procedure is OK ? If it is not OK,
>> what is the right
>> > one ?
>> >
>> > 2、Which atoms and which orbitals should I treat
>> with orbital
>> > potential ? The electron configurations of these 3
>> atoms are: La (5d1
>> > 6s2) ; Pt (4f14 5d9 6s1); Bi (4f14 5d10 6s2
>> 6p3).
>> >
>> > Thanks very much for your attention.
>> >
>> > Sincerely yours,
>> >
>> > Bingrui Peng
>> >
>> > from the Department of Physics, Nanjing
>> University, China
>>
>>
>>
>
>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
> SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
>
--
--------------------------------------------------------------------------
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-165300 FAX: +43-1-58801-165982
Email: blaha at theochem.tuwien.ac.at WIEN2k: http://www.wien2k.at
WWW:
http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: Na_B100_irreps.txt
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20170719/3531c1df/attachment.txt>
More information about the Wien
mailing list