Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] How to deal with the output files?

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] How to deal with the output files?


Chronological Thread 
  • From: "Huiyong Deng" <hydeng2007@gmail.com>
  • To: forum@abinit.org
  • Subject: Re: [abinit-forum] How to deal with the output files?
  • Date: Wed, 13 Sep 2006 16:30:28 +0800
  • Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=beta; d=gmail.com; h=received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:references; b=pE+nDdCzlZ3CQsSoIkWE/dBQFfk2errSNOFWbexJ8lsXwdUwREjStWEPy0J37sojtnrvyJ6hNoaealXuKWwrutrw6KgZekMPLz7F86biZtQJdodcuCAlT9r3l+eS7hX6sn79ZsPyZodl5IHndK5OPqMQ/Bv5/66tjNKPDwApLTQ=

Nuno A. G. Bandeira,
 
Having read your explation about Mulliken population , however, I still don't know whether the Abinit can perform the  Mulliken population calculation. My questions is: 
How to perform the Mulliken population calculation using the Abinit?Thanks you!
 
Thanks Anglade Pierre-Matthieu for his kind help! Maybe I've known how to deal with the result now.
 
These days, I've performed the GW calculation. However,after finishing the KSS calculation (DATASET1 )  and begin to read some data of the DATASET2, the Abinit aborts and has no ERROR  tips. Why? What shall I do?  Is there some detailed literature about GW?
 
 


 
2006/9/13, Nuno A. G. Bandeira <nuno.bandeira@ist.utl.pt>:
Anglade Pierre-Matthieu wrote:

>>     (2) Can the Abinit perform the calculation of the electronic Mulliken
>> population analysis ?
>
> Sorry I don't know.

Say you have an MO with 2 AOs:

PHI=c1.chi1+c2.chi2

When you square it:

PHI^2=c1^2.chi1^2 + c2^2.chi2^2 + 2c1c2.chi1.chi2

when you integrate over all space you have:

1= c1^2 + c2^2 + 2c1c2.S12

2c1c2.S12 is called the "(Mulliken) overlap population" (S12 being the
overlap integral) and it is in a way a measure of the bond strength.
There's a problem here with the cross term 2c1c2.chi1.chi2 because you
can't assign the overlap charge density to any one AO. So what Mulliken
did was to democratically assign the charge density equally between the
2 AOs thus:

chi1.chi2 ~ (chi1^2+chi2^2)/2

so this way you have "Mulliken gross atomic population" in atoms 1 and 2

PHI^2 = (c1^2+c1c2)chi1^2 + (c2^2+c1c2)chi2^2

N1 = c1^2+c1c2
N2 = c2^2+c1c2

With these values you can obviously define a Mulliken charge. With big
basis sets some values may become chemically counterintuitive, for
instance with diffuse gaussian basis functions, but they are nonetheless
widely used as the main method of population analysis.

The Mulliken population has been extended to solids by Hoffmann with the
creation of COOP (crystal orbital overlap population). You can read his
seminal article in JACS,1983,105,3528 and in his book (Solids and
Surfaces: A Chemist's view of bonding in extended structures, Wiley-VCH)


Regards,
--
Nuno A. G. Bandeira, AMRSC
Graduate researcher and molecular sculptor
Inorganic and Theoretical Chemistry Group,
Faculty of Science
University of Lisbon - C8 building, Campo Grande,
1749-016 Lisbon,Portugal
http://cqb.fc.ul.pt/intheochem/nuno.html
Doctoral student @ IST,Lisbon
--




Archive powered by MHonArc 2.6.16.

Top of Page