Axel,<div><br></div><div>  Thanks for your reply. Now, I understand why the output quadrupole moment by CP2k is not traceless. But, it seems to me that the six quadrupole tensor components are over complete. In fact, only five variables are needed for uniqueness when spherical coordinate notation is used.</div>
<div><br></div><div>  Anyway, the referred page is very helpful to clarify those commonly used definitions of quadrupole moment. </div><div><br></div><div>  Thanks.</div><div><br></div><div>Hanning<br><br><div class="gmail_quote">
On Sun, Jun 20, 2010 at 11:54 AM, Axel <span dir="ltr"><<a href="mailto:akoh...@gmail.com">akoh...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><br>
<br>
On Jun 20, 1:04 am, Hanning Chen <<a href="mailto:chenh...@gmail.com">chenh...@gmail.com</a>> wrote:<br>
> Dear CP2k users,<br>
><br>
>   I ran a simple SCF calculation on CO2 and found that its output quadrupole<br>
> moment (as shown below) apparently does not satisfy the requirement of<br>
> traceless, i.e., Q(XX)+Q(YY)+Q(ZZ)==0.<br>
<br>
>   Does CP2k have its own definition of quadrupole moment rather than the<br>
> standard one?<br>
<br>
</div>there is not one but two "standard" definitions.<br>
a quick search on google reveals this page, for example:<br>
<br>
<a href="http://cccbdb.nist.gov/quadrupole.asp" target="_blank">http://cccbdb.nist.gov/quadrupole.asp</a><br>
<br>
axel.<br>
<div><div></div><div class="h5"><br>
>   Thanks.<br>
><br>
> Hanning Chen<br>
> Ph.D.<br>
> Department of Chemistry<br>
> Northwestern University<br>
> Evasnton, IL 60201<br>
><br>
> OUTPUT_FILE:<br>
><br>
>  Reference Point [Bohr]             9.44863138      9.44863173<br>
> 9.44499632<br>
>   Charges<br>
>     Electronic=      -16.0000      Core=      16.0000      Total=<br>
> 0.0000<br>
>   Dipoles are based on the traditional operator.<br>
>   Dipole moment [Debye]<br>
>     X=          0.000478093061    Y=         -0.000005701601    Z=<br>
> 0.005608188292                   Total=            0.005628532793<br>
>   Quadrupole moment [Debye*Angstrom]<br>
>                 XX=        -18.812852208659         XY=<br>
> -0.089714841097         XZ=         -0.083871152755<br>
>                 YY=        -14.716948452691         YZ=<br>
> -0.002201287063         ZZ=        -14.716586168813<br>
><br>
> INPUT_FILE:<br>
><br>
> &FORCE_EVAL<br>
>   METHOD Quickstep<br>
>   &DFT<br>
>     &MGRID<br>
>       CUTOFF 300<br>
>     &END MGRID<br>
>     &QS<br>
>       EPS_DEFAULT 1.0E-12<br>
>     &END QS<br>
>     &SCF<br>
>       SCF_GUESS RESTART<br>
>       EPS_SCF 1.0E-10<br>
>       MAX_SCF 1000<br>
>       &OT<br>
>        MINIMIZER DIIS<br>
>        PRECONDITIONER FULL_KINETIC<br>
>       &END OT<br>
>     &END SCF<br>
>     &XC<br>
>       &XC_FUNCTIONAL PBE<br>
>       &END XC_FUNCTIONAL<br>
>     &END XC<br>
>     &POISSON<br>
>      POISSON_SOLVER WAVELET<br>
>      PERIODIC NONE<br>
>     &END POISSON<br>
>     &PRINT<br>
>      &MOMENTS<br>
>       FILENAME =dipole<br>
>       REFERENCE COM<br>
>       PERIODIC = F<br>
>       MAX_MOMENT 2<br>
>      &END MOMENTS<br>
>     &END PRINT<br>
>   &END DFT<br>
>   &SUBSYS<br>
>     &CELL<br>
>       ABC 10.0 10.0 10.0<br>
>       PERIODIC NONE<br>
>     &END CELL<br>
>     &COORD<br>
>   C         5.0000000838        5.0000000000        4.9999999999<br>
>   O         3.8098548376        4.9999999994        4.9999999994<br>
>   O         6.1901450757        5.0000000006        5.0000000006<br>
>     &END COORD<br>
>     &KIND C<br>
>       BASIS_SET DZVP-GTH-PBE<br>
>       POTENTIAL GTH-PBE-q4<br>
>     &END KIND<br>
>     &KIND O<br>
>       BASIS_SET DZVP-GTH-PBE<br>
>       POTENTIAL GTH-PBE-q6<br>
>     &END KIND<br>
>   &END SUBSYS<br>
> &END FORCE_EVAL<br>
> &GLOBAL<br>
>   PROJECT CO2<br>
>   RUN_TYPE ENERGY<br>
>   PRINT_LEVEL MEDIUM<br>
> &END GLOBAL<br>
<br>
</div></div><font color="#888888">--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
For more options, visit this group at <a href="http://groups.google.com/group/cp2k?hl=en" target="_blank">http://groups.google.com/group/cp2k?hl=en</a>.<br>
<br>
</font></blockquote></div><br></div>