forum@abinit.org
Subject: The ABINIT Users Mailing List ( CLOSED )
List archive
- 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
--
- How to deal with the output files?, Huiyong Deng, 09/06/2006
- Re: [abinit-forum] How to deal with the output files?, Anglade Pierre-Matthieu, 09/06/2006
- Re: [abinit-forum] How to deal with the output files?, Huiyong Deng, 09/06/2006
- Re: [abinit-forum] How to deal with the output files?, Anglade Pierre-Matthieu, 09/06/2006
- Re: [abinit-forum] How to deal with the output files?, Nuno A. G. Bandeira, 09/12/2006
- Re: [abinit-forum] How to deal with the output files?, Matthieu Verstraete, 09/12/2006
- Re: [abinit-forum] How to deal with the output files?, Huiyong Deng, 09/13/2006
- Re: [abinit-forum] How to deal with the output files?, Nuno A. G. Bandeira, 09/13/2006
- Re: [abinit-forum] How to deal with the output files?, Nuno A. G. Bandeira, 09/12/2006
- Re: [abinit-forum] How to deal with the output files?, Nuno A. G. Bandeira, 09/06/2006
- Re: [abinit-forum] How to deal with the output files?, Anglade Pierre-Matthieu, 09/06/2006
- Re: [abinit-forum] How to deal with the output files?, Huiyong Deng, 09/06/2006
- Re: [abinit-forum] How to deal with the output files?, Anglade Pierre-Matthieu, 09/06/2006
Archive powered by MHonArc 2.6.16.