Re: [SIESTA-L] << two set of charge calculations >>

[Emilio Artacho]

oops, sorry Camps, my lecturing you was not pertinent!

E

On Jun 20, 2022, at 10:12 PM, I. Camps 
mailto:ica...@gmail.com>> wrote:

I completely agree with Tamas and Emilio BUT my question is not related to 
which charge calculation scheme is "better".

My question is that in my calculations, two different sets of data of the same 
type of charges are appearing in the output file, instead only one for each. I 
have two outputs for Hirshfeld and two outputs for Voronoi charges.

[]'s,

Camps


On Mon, Jun 20, 2022 at 5:02 PM Emilio Artacho 
mailto:e.arta...@nanogune.eu>> wrote:
Tamas’s reply is correct, I just want to add a reminder of the fact
that atomic charges have a fundamental definition problem and none of
the proposals gives the ‘good’ answer. This is a direct consequence
of its responding to an ill-posed question: how many electrons ‘belong’
to a given atom (or can be assigned to it). It is perfectly defined if the atoms
are infinitely separated from each other, but not otherwise.

It is clear, however, that concepts like charge transfer etc are useful
in chemistry and very much support chemical analysis and intuition.
Atomic charges schemes (when used sensibly) are valuable. Just remember
to use them with care (qualitatively, trends etc). There are good comparative
studies assessing their reliability in various chemistry situations.

There are situations for which the question can be rephrased
into something physically well defined (see e,g, the Born effective
charges, or other questions relating to dielectric polarisation).

One can also find claims in the literature for a particular scheme to be
the ‘right’ one. To my mind they all rely on arbitrary choices, which can
be more or less sensible or well motivated, but still arbitrary (as Tamas
says, some depend on the basis set choice while other do not, for
instance).

best

Emilio

On Jun 19, 2022, at 2:47 PM, Tamas Karpati 
mailto:tkarp...@gmail.com>> wrote:

Dear Camps,

Please note that an argument is going on for decades about how to
calculate atomic charges. Different methods/schemes give different
results, each is giving better/worse results for different
applications. It is recommended to check how well each performs at
your actual problem and choose which one is to be used. Also
remarkable is the basis set dependence of atomic charges, consider
this a parameter to be calibrated.

Regards,
 t

On Fri, Jun 17, 2022 at 10:02 PM I. Camps 
mailto:ica...@gmail.com>> wrote:

Hello Alberto,

Here it is the info about the SIESTA version:

Siesta Version  : siesta-max-R3--710-676-597
Architecture: unknown
Compiler version: ifort (IFORT) 19.1.1.217 20200306
Compiler flags  : mpifort -fPIC -O2 -march=core-avx2 -axCore-AVX512 -fp-model 
precise
PP flags: -DFC_HAVE_ABORT -DF2003 -DMPI -DCDF -DNCDF -DNCDF_4 
-DNCDF_PARALLEL 
-I/cvmfs//soft.computecanada.ca/easybuild/software/2020/avx2/MPI/intel2020/openmpi4/netcdf-fortran-mpi/4.5.2/include
Libraries   : libncdf.a libfdict.a -Wl,-Bstatic -Wl,--start-group 
-lmkl_scalapack_lp64 -lmkkl_blacs_openmpi_lp64 -lmkl_intel_lp64 
-lmkl_sequential -lmkl_core -Wl,--end-group -Wl,-Bdynamic -lnetcdff
PARALLEL version
NetCDF support
NetCDF-4 support
NetCDF-4 MPI-IO support

And here is the output section:

siesta: Final energy (eV):
siesta:  Band Struct. =   -8272.290139
siesta:   Kinetic =   19960.524774
siesta:   Hartree =  151423.860682
siesta:   Eldau   =   0.00
siesta:   Eso =   0.00
siesta:Ext. field =   0.00
siesta:   Enegf   =   0.00
siesta:   Exch.-corr. =  -11180.064205
siesta:  Ion-electron = -320401.282309
siesta:   Ion-ion =  129282.468462
siesta:   Ekinion =   0.00
siesta: Total =  -30914.492596
siesta: Fermi =  -4.212218

siesta: Stress tensor (static) (eV/Ang**3):
siesta: 0.0001260.00   -0.00
siesta: 0.000.000101   -0.49
siesta:-0.00   -0.49   -0.016465

siesta: Cell volume =   7672.635004 Ang**3

siesta: Pressure (static):
siesta:SolidMolecule  Units
siesta:   0.5895  0.5941  Ry/Bohr**3
siesta:   0.00541292  0.00545494  eV/Ang**3
siesta:   8.67254766  8.73987328  kBar
(Free)E+ p_basis*V_orbitals  =  -30859.763440
(Free)Eharris+ p_basis*V_orbitals  =  -30859.763491
spin moment: S , {S} =0.0   0.0   0.0   0.0

siesta: Electric dipole (a.u.)  =0.000.0432460.00
siesta: Electric dipole (Debye) =0.010.1099190.00

Hirshfeld Net Atomic Populations:
Atom #Qatom  Species
10.149  B
20.149  B
30.149  B
40.149  B
5   -0.149  N
...
  155   -0.149  N
  156   -0.149  N
  1570.149  B
  1580.149  B
  1590.149  B
  1600.149  B

Voronoi Net 

Re: [SIESTA-L] << two set of charge calculations >>

[Nick Papior]

The 2nd data output relates to your ldos range. So they are effectively not
equivalent. The first is the final scf, the 2nd ldos.

I think this is the case, but I am not 100 sure (not at my computer to
check).

On Tue, 21 Jun 2022, 22:00 I. Camps,  wrote:

> I completely agree with Tamas and Emilio BUT my question *is not* related
> to which charge calculation scheme is "better".
>
> My question is that in my calculations, two different sets of data of the
> same type of charges are appearing in the output file, instead only one for
> each. I have two outputs for Hirshfeld and two outputs for Voronoi charges.
>
> []'s,
>
> Camps
>
>
> On Mon, Jun 20, 2022 at 5:02 PM Emilio Artacho 
> wrote:
>
>> Tamas’s reply is correct, I just want to add a reminder of the fact
>> that atomic charges have a fundamental definition problem and none of
>> the proposals gives the ‘good’ answer. This is a direct consequence
>> of its responding to an ill-posed question: how many electrons ‘belong’
>> to a given atom (or can be assigned to it). It is perfectly defined if
>> the atoms
>> are infinitely separated from each other, but not otherwise.
>>
>> It is clear, however, that concepts like charge transfer etc are useful
>> in chemistry and very much support chemical analysis and intuition.
>> Atomic charges schemes (when used sensibly) are valuable. Just remember
>> to use them with care (qualitatively, trends etc). There are good
>> comparative
>> studies assessing their reliability in various chemistry situations.
>>
>> There are situations for which the question can be rephrased
>> into something physically well defined (see e,g, the Born effective
>> charges, or other questions relating to dielectric polarisation).
>>
>> One can also find claims in the literature for a particular scheme to be
>> the ‘right’ one. To my mind they all rely on arbitrary choices, which can
>> be more or less sensible or well motivated, but still arbitrary (as Tamas
>> says, some depend on the basis set choice while other do not, for
>> instance).
>>
>> best
>>
>> Emilio
>>
>> On Jun 19, 2022, at 2:47 PM, Tamas Karpati  wrote:
>>
>> Dear Camps,
>>
>> Please note that an argument is going on for decades about how to
>> calculate atomic charges. Different methods/schemes give different
>> results, each is giving better/worse results for different
>> applications. It is recommended to check how well each performs at
>> your actual problem and choose which one is to be used. Also
>> remarkable is the basis set dependence of atomic charges, consider
>> this a parameter to be calibrated.
>>
>> Regards,
>>  t
>>
>> On Fri, Jun 17, 2022 at 10:02 PM I. Camps  wrote:
>>
>>
>> Hello Alberto,
>>
>> Here it is the info about the SIESTA version:
>>
>> Siesta Version  : siesta-max-R3--710-676-597
>> Architecture: unknown
>> Compiler version: ifort (IFORT) 19.1.1.217 20200306
>> Compiler flags  : mpifort -fPIC -O2 -march=core-avx2 -axCore-AVX512
>> -fp-model precise
>> PP flags: -DFC_HAVE_ABORT -DF2003 -DMPI -DCDF -DNCDF -DNCDF_4
>> -DNCDF_PARALLEL -I/cvmfs//
>> soft.computecanada.ca/easybuild/software/2020/avx2/MPI/intel2020/openmpi4/netcdf-fortran-mpi/4.5.2/include
>> Libraries   : libncdf.a libfdict.a -Wl,-Bstatic -Wl,--start-group
>> -lmkl_scalapack_lp64 -lmkkl_blacs_openmpi_lp64 -lmkl_intel_lp64
>> -lmkl_sequential -lmkl_core -Wl,--end-group -Wl,-Bdynamic -lnetcdff
>> PARALLEL version
>> NetCDF support
>> NetCDF-4 support
>> NetCDF-4 MPI-IO support
>>
>> And here is the output section:
>>
>> siesta: Final energy (eV):
>> siesta:  Band Struct. =   -8272.290139
>> siesta:   Kinetic =   19960.524774
>> siesta:   Hartree =  151423.860682
>> siesta:   Eldau   =   0.00
>> siesta:   Eso =   0.00
>> siesta:Ext. field =   0.00
>> siesta:   Enegf   =   0.00
>> siesta:   Exch.-corr. =  -11180.064205
>> siesta:  Ion-electron = -320401.282309
>> siesta:   Ion-ion =  129282.468462
>> siesta:   Ekinion =   0.00
>> siesta: Total =  -30914.492596
>> siesta: Fermi =  -4.212218
>>
>> siesta: Stress tensor (static) (eV/Ang**3):
>> siesta: 0.0001260.00   -0.00
>> siesta: 0.000.000101   -0.49
>> siesta:-0.00   -0.49   -0.016465
>>
>> siesta: Cell volume =   7672.635004 Ang**3
>>
>> siesta: Pressure (static):
>> siesta:SolidMolecule  Units
>> siesta:   0.5895  0.5941  Ry/Bohr**3
>> siesta:   0.00541292  0.00545494  eV/Ang**3
>> siesta:   8.67254766  8.73987328  kBar
>> (Free)E+ p_basis*V_orbitals  =  -30859.763440
>> (Free)Eharris+ p_basis*V_orbitals  =  -30859.763491
>> spin moment: S , {S} =0.0   0.0   0.0   0.0
>>
>> siesta: Electric dipole (a.u.)  =0.000.0432460.00
>> siesta: Electric dipole (Debye) =0.010.1099190.00
>>
>> Hirshfeld Net Atomic Populations:
>> 

Re: [SIESTA-L] << two set of charge calculations >>

[Tamas Karpati]

Dear Camps, Emilio,

True, I overlooked the point in the original question (my apologies!)
but thought that the doubling of charge prints is simply due to a spin
polarized calculation.

I was wrong, as that doubling also appears in case of closed shell
singlets' outputs. What is more, the two lists are quite different.

Conclusion: I don't know the reason (only I can think of the initial
and final data of an iterative charge calculation but it's just an
idea).

Question: Could someone clarify this?

Best regards,
  t


On Tue, Jun 21, 2022 at 10:03 PM I. Camps  wrote:
>
> I completely agree with Tamas and Emilio BUT my question is not related to 
> which charge calculation scheme is "better".
>
> My question is that in my calculations, two different sets of data of the 
> same type of charges are appearing in the output file, instead only one for 
> each. I have two outputs for Hirshfeld and two outputs for Voronoi charges.
>
> []'s,
>
> Camps
>
>
> On Mon, Jun 20, 2022 at 5:02 PM Emilio Artacho  wrote:
>>
>> Tamas’s reply is correct, I just want to add a reminder of the fact
>> that atomic charges have a fundamental definition problem and none of
>> the proposals gives the ‘good’ answer. This is a direct consequence
>> of its responding to an ill-posed question: how many electrons ‘belong’
>> to a given atom (or can be assigned to it). It is perfectly defined if the 
>> atoms
>> are infinitely separated from each other, but not otherwise.
>>
>> It is clear, however, that concepts like charge transfer etc are useful
>> in chemistry and very much support chemical analysis and intuition.
>> Atomic charges schemes (when used sensibly) are valuable. Just remember
>> to use them with care (qualitatively, trends etc). There are good comparative
>> studies assessing their reliability in various chemistry situations.
>>
>> There are situations for which the question can be rephrased
>> into something physically well defined (see e,g, the Born effective
>> charges, or other questions relating to dielectric polarisation).
>>
>> One can also find claims in the literature for a particular scheme to be
>> the ‘right’ one. To my mind they all rely on arbitrary choices, which can
>> be more or less sensible or well motivated, but still arbitrary (as Tamas
>> says, some depend on the basis set choice while other do not, for
>> instance).
>>
>> best
>>
>> Emilio
>>
>> On Jun 19, 2022, at 2:47 PM, Tamas Karpati  wrote:
>>
>> Dear Camps,
>>
>> Please note that an argument is going on for decades about how to
>> calculate atomic charges. Different methods/schemes give different
>> results, each is giving better/worse results for different
>> applications. It is recommended to check how well each performs at
>> your actual problem and choose which one is to be used. Also
>> remarkable is the basis set dependence of atomic charges, consider
>> this a parameter to be calibrated.
>>
>> Regards,
>>  t
>>
>> On Fri, Jun 17, 2022 at 10:02 PM I. Camps  wrote:
>>
>>
>> Hello Alberto,
>>
>> Here it is the info about the SIESTA version:
>>
>> Siesta Version  : siesta-max-R3--710-676-597
>> Architecture: unknown
>> Compiler version: ifort (IFORT) 19.1.1.217 20200306
>> Compiler flags  : mpifort -fPIC -O2 -march=core-avx2 -axCore-AVX512 
>> -fp-model precise
>> PP flags: -DFC_HAVE_ABORT -DF2003 -DMPI -DCDF -DNCDF -DNCDF_4 
>> -DNCDF_PARALLEL 
>> -I/cvmfs//soft.computecanada.ca/easybuild/software/2020/avx2/MPI/intel2020/openmpi4/netcdf-fortran-mpi/4.5.2/include
>> Libraries   : libncdf.a libfdict.a -Wl,-Bstatic -Wl,--start-group 
>> -lmkl_scalapack_lp64 -lmkkl_blacs_openmpi_lp64 -lmkl_intel_lp64 
>> -lmkl_sequential -lmkl_core -Wl,--end-group -Wl,-Bdynamic -lnetcdff
>> PARALLEL version
>> NetCDF support
>> NetCDF-4 support
>> NetCDF-4 MPI-IO support
>>
>> And here is the output section:
>>
>> siesta: Final energy (eV):
>> siesta:  Band Struct. =   -8272.290139
>> siesta:   Kinetic =   19960.524774
>> siesta:   Hartree =  151423.860682
>> siesta:   Eldau   =   0.00
>> siesta:   Eso =   0.00
>> siesta:Ext. field =   0.00
>> siesta:   Enegf   =   0.00
>> siesta:   Exch.-corr. =  -11180.064205
>> siesta:  Ion-electron = -320401.282309
>> siesta:   Ion-ion =  129282.468462
>> siesta:   Ekinion =   0.00
>> siesta: Total =  -30914.492596
>> siesta: Fermi =  -4.212218
>>
>> siesta: Stress tensor (static) (eV/Ang**3):
>> siesta: 0.0001260.00   -0.00
>> siesta: 0.000.000101   -0.49
>> siesta:-0.00   -0.49   -0.016465
>>
>> siesta: Cell volume =   7672.635004 Ang**3
>>
>> siesta: Pressure (static):
>> siesta:SolidMolecule  Units
>> siesta:   0.5895  0.5941  Ry/Bohr**3
>> siesta:   0.00541292  0.00545494  eV/Ang**3
>> siesta:   8.67254766  8.73987328  kBar
>> (Free)E+ p_basis*V_orbitals  =  -30859.763440
>> 

Re: [SIESTA-L] << two set of charge calculations >>

[I. Camps]

I completely agree with Tamas and Emilio BUT my question *is not* related
to which charge calculation scheme is "better".

My question is that in my calculations, two different sets of data of the
same type of charges are appearing in the output file, instead only one for
each. I have two outputs for Hirshfeld and two outputs for Voronoi charges.

[]'s,

Camps


On Mon, Jun 20, 2022 at 5:02 PM Emilio Artacho 
wrote:

> Tamas’s reply is correct, I just want to add a reminder of the fact
> that atomic charges have a fundamental definition problem and none of
> the proposals gives the ‘good’ answer. This is a direct consequence
> of its responding to an ill-posed question: how many electrons ‘belong’
> to a given atom (or can be assigned to it). It is perfectly defined if the
> atoms
> are infinitely separated from each other, but not otherwise.
>
> It is clear, however, that concepts like charge transfer etc are useful
> in chemistry and very much support chemical analysis and intuition.
> Atomic charges schemes (when used sensibly) are valuable. Just remember
> to use them with care (qualitatively, trends etc). There are good
> comparative
> studies assessing their reliability in various chemistry situations.
>
> There are situations for which the question can be rephrased
> into something physically well defined (see e,g, the Born effective
> charges, or other questions relating to dielectric polarisation).
>
> One can also find claims in the literature for a particular scheme to be
> the ‘right’ one. To my mind they all rely on arbitrary choices, which can
> be more or less sensible or well motivated, but still arbitrary (as Tamas
> says, some depend on the basis set choice while other do not, for
> instance).
>
> best
>
> Emilio
>
> On Jun 19, 2022, at 2:47 PM, Tamas Karpati  wrote:
>
> Dear Camps,
>
> Please note that an argument is going on for decades about how to
> calculate atomic charges. Different methods/schemes give different
> results, each is giving better/worse results for different
> applications. It is recommended to check how well each performs at
> your actual problem and choose which one is to be used. Also
> remarkable is the basis set dependence of atomic charges, consider
> this a parameter to be calibrated.
>
> Regards,
>  t
>
> On Fri, Jun 17, 2022 at 10:02 PM I. Camps  wrote:
>
>
> Hello Alberto,
>
> Here it is the info about the SIESTA version:
>
> Siesta Version  : siesta-max-R3--710-676-597
> Architecture: unknown
> Compiler version: ifort (IFORT) 19.1.1.217 20200306
> Compiler flags  : mpifort -fPIC -O2 -march=core-avx2 -axCore-AVX512
> -fp-model precise
> PP flags: -DFC_HAVE_ABORT -DF2003 -DMPI -DCDF -DNCDF -DNCDF_4
> -DNCDF_PARALLEL -I/cvmfs//
> soft.computecanada.ca/easybuild/software/2020/avx2/MPI/intel2020/openmpi4/netcdf-fortran-mpi/4.5.2/include
> Libraries   : libncdf.a libfdict.a -Wl,-Bstatic -Wl,--start-group
> -lmkl_scalapack_lp64 -lmkkl_blacs_openmpi_lp64 -lmkl_intel_lp64
> -lmkl_sequential -lmkl_core -Wl,--end-group -Wl,-Bdynamic -lnetcdff
> PARALLEL version
> NetCDF support
> NetCDF-4 support
> NetCDF-4 MPI-IO support
>
> And here is the output section:
>
> siesta: Final energy (eV):
> siesta:  Band Struct. =   -8272.290139
> siesta:   Kinetic =   19960.524774
> siesta:   Hartree =  151423.860682
> siesta:   Eldau   =   0.00
> siesta:   Eso =   0.00
> siesta:Ext. field =   0.00
> siesta:   Enegf   =   0.00
> siesta:   Exch.-corr. =  -11180.064205
> siesta:  Ion-electron = -320401.282309
> siesta:   Ion-ion =  129282.468462
> siesta:   Ekinion =   0.00
> siesta: Total =  -30914.492596
> siesta: Fermi =  -4.212218
>
> siesta: Stress tensor (static) (eV/Ang**3):
> siesta: 0.0001260.00   -0.00
> siesta: 0.000.000101   -0.49
> siesta:-0.00   -0.49   -0.016465
>
> siesta: Cell volume =   7672.635004 Ang**3
>
> siesta: Pressure (static):
> siesta:SolidMolecule  Units
> siesta:   0.5895  0.5941  Ry/Bohr**3
> siesta:   0.00541292  0.00545494  eV/Ang**3
> siesta:   8.67254766  8.73987328  kBar
> (Free)E+ p_basis*V_orbitals  =  -30859.763440
> (Free)Eharris+ p_basis*V_orbitals  =  -30859.763491
> spin moment: S , {S} =0.0   0.0   0.0   0.0
>
> siesta: Electric dipole (a.u.)  =0.000.0432460.00
> siesta: Electric dipole (Debye) =0.010.1099190.00
>
> Hirshfeld Net Atomic Populations:
> Atom #Qatom  Species
> 10.149  B
> 20.149  B
> 30.149  B
> 40.149  B
> 5   -0.149  N
> ...
>   155   -0.149  N
>   156   -0.149  N
>   1570.149  B
>   1580.149  B
>   1590.149  B
>   1600.149  B
>
> Voronoi Net Atomic Populations:
> Atom #Qatom  Species
> 10.167  B
> 20.167  B
> 30.167  B
> 4

Re: [SIESTA-L] << two set of charge calculations >>

[Emilio Artacho]

Tamas’s reply is correct, I just want to add a reminder of the fact
that atomic charges have a fundamental definition problem and none of
the proposals gives the ‘good’ answer. This is a direct consequence
of its responding to an ill-posed question: how many electrons ‘belong’
to a given atom (or can be assigned to it). It is perfectly defined if the atoms
are infinitely separated from each other, but not otherwise.

It is clear, however, that concepts like charge transfer etc are useful
in chemistry and very much support chemical analysis and intuition.
Atomic charges schemes (when used sensibly) are valuable. Just remember
to use them with care (qualitatively, trends etc). There are good comparative
studies assessing their reliability in various chemistry situations.

There are situations for which the question can be rephrased
into something physically well defined (see e,g, the Born effective
charges, or other questions relating to dielectric polarisation).

One can also find claims in the literature for a particular scheme to be
the ‘right’ one. To my mind they all rely on arbitrary choices, which can
be more or less sensible or well motivated, but still arbitrary (as Tamas
says, some depend on the basis set choice while other do not, for
instance).

best

Emilio

On Jun 19, 2022, at 2:47 PM, Tamas Karpati 
mailto:tkarp...@gmail.com>> wrote:

Dear Camps,

Please note that an argument is going on for decades about how to
calculate atomic charges. Different methods/schemes give different
results, each is giving better/worse results for different
applications. It is recommended to check how well each performs at
your actual problem and choose which one is to be used. Also
remarkable is the basis set dependence of atomic charges, consider
this a parameter to be calibrated.

Regards,
 t

On Fri, Jun 17, 2022 at 10:02 PM I. Camps 
mailto:ica...@gmail.com>> wrote:

Hello Alberto,

Here it is the info about the SIESTA version:

Siesta Version  : siesta-max-R3--710-676-597
Architecture: unknown
Compiler version: ifort (IFORT) 19.1.1.217 20200306
Compiler flags  : mpifort -fPIC -O2 -march=core-avx2 -axCore-AVX512 -fp-model 
precise
PP flags: -DFC_HAVE_ABORT -DF2003 -DMPI -DCDF -DNCDF -DNCDF_4 
-DNCDF_PARALLEL 
-I/cvmfs//soft.computecanada.ca/easybuild/software/2020/avx2/MPI/intel2020/openmpi4/netcdf-fortran-mpi/4.5.2/include
Libraries   : libncdf.a libfdict.a -Wl,-Bstatic -Wl,--start-group 
-lmkl_scalapack_lp64 -lmkkl_blacs_openmpi_lp64 -lmkl_intel_lp64 
-lmkl_sequential -lmkl_core -Wl,--end-group -Wl,-Bdynamic -lnetcdff
PARALLEL version
NetCDF support
NetCDF-4 support
NetCDF-4 MPI-IO support

And here is the output section:

siesta: Final energy (eV):
siesta:  Band Struct. =   -8272.290139
siesta:   Kinetic =   19960.524774
siesta:   Hartree =  151423.860682
siesta:   Eldau   =   0.00
siesta:   Eso =   0.00
siesta:Ext. field =   0.00
siesta:   Enegf   =   0.00
siesta:   Exch.-corr. =  -11180.064205
siesta:  Ion-electron = -320401.282309
siesta:   Ion-ion =  129282.468462
siesta:   Ekinion =   0.00
siesta: Total =  -30914.492596
siesta: Fermi =  -4.212218

siesta: Stress tensor (static) (eV/Ang**3):
siesta: 0.0001260.00   -0.00
siesta: 0.000.000101   -0.49
siesta:-0.00   -0.49   -0.016465

siesta: Cell volume =   7672.635004 Ang**3

siesta: Pressure (static):
siesta:SolidMolecule  Units
siesta:   0.5895  0.5941  Ry/Bohr**3
siesta:   0.00541292  0.00545494  eV/Ang**3
siesta:   8.67254766  8.73987328  kBar
(Free)E+ p_basis*V_orbitals  =  -30859.763440
(Free)Eharris+ p_basis*V_orbitals  =  -30859.763491
spin moment: S , {S} =0.0   0.0   0.0   0.0

siesta: Electric dipole (a.u.)  =0.000.0432460.00
siesta: Electric dipole (Debye) =0.010.1099190.00

Hirshfeld Net Atomic Populations:
Atom #Qatom  Species
10.149  B
20.149  B
30.149  B
40.149  B
5   -0.149  N
...
  155   -0.149  N
  156   -0.149  N
  1570.149  B
  1580.149  B
  1590.149  B
  1600.149  B

Voronoi Net Atomic Populations:
Atom #Qatom  Species
10.167  B
20.167  B
30.167  B
40.167  B
5   -0.168  N
...
  155   -0.168  N
  156   -0.168  N
  1570.168  B
  1580.168  B
  1590.168  B
  1600.168  B
Bader Analysis core-charge setup. Radii (standard, H):  1.000 0.600

dhscf: Vacuum level (max, mean) =   -0.038479   -0.112800 eV

siesta: LDOS info
siesta: E1 -- E2 [eV]:  -20.000 --0.000

Hirshfeld Net Atomic Populations:
Atom #Qatom  Species
10.227  B
20.227  B
30.227  B
40.227  B
50.888 

Re: [SIESTA-L] << two set of charge calculations >>

[Tamas Karpati]

Dear Camps,

Please note that an argument is going on for decades about how to
calculate atomic charges. Different methods/schemes give different
results, each is giving better/worse results for different
applications. It is recommended to check how well each performs at
your actual problem and choose which one is to be used. Also
remarkable is the basis set dependence of atomic charges, consider
this a parameter to be calibrated.

Regards,
  t

On Fri, Jun 17, 2022 at 10:02 PM I. Camps  wrote:
>
> Hello Alberto,
>
> Here it is the info about the SIESTA version:
>
> Siesta Version  : siesta-max-R3--710-676-597
> Architecture: unknown
> Compiler version: ifort (IFORT) 19.1.1.217 20200306
> Compiler flags  : mpifort -fPIC -O2 -march=core-avx2 -axCore-AVX512 -fp-model 
> precise
> PP flags: -DFC_HAVE_ABORT -DF2003 -DMPI -DCDF -DNCDF -DNCDF_4 
> -DNCDF_PARALLEL 
> -I/cvmfs//soft.computecanada.ca/easybuild/software/2020/avx2/MPI/intel2020/openmpi4/netcdf-fortran-mpi/4.5.2/include
> Libraries   : libncdf.a libfdict.a -Wl,-Bstatic -Wl,--start-group 
> -lmkl_scalapack_lp64 -lmkkl_blacs_openmpi_lp64 -lmkl_intel_lp64 
> -lmkl_sequential -lmkl_core -Wl,--end-group -Wl,-Bdynamic -lnetcdff
> PARALLEL version
> NetCDF support
> NetCDF-4 support
> NetCDF-4 MPI-IO support
>
> And here is the output section:
>
> siesta: Final energy (eV):
> siesta:  Band Struct. =   -8272.290139
> siesta:   Kinetic =   19960.524774
> siesta:   Hartree =  151423.860682
> siesta:   Eldau   =   0.00
> siesta:   Eso =   0.00
> siesta:Ext. field =   0.00
> siesta:   Enegf   =   0.00
> siesta:   Exch.-corr. =  -11180.064205
> siesta:  Ion-electron = -320401.282309
> siesta:   Ion-ion =  129282.468462
> siesta:   Ekinion =   0.00
> siesta: Total =  -30914.492596
> siesta: Fermi =  -4.212218
>
> siesta: Stress tensor (static) (eV/Ang**3):
> siesta: 0.0001260.00   -0.00
> siesta: 0.000.000101   -0.49
> siesta:-0.00   -0.49   -0.016465
>
> siesta: Cell volume =   7672.635004 Ang**3
>
> siesta: Pressure (static):
> siesta:SolidMolecule  Units
> siesta:   0.5895  0.5941  Ry/Bohr**3
> siesta:   0.00541292  0.00545494  eV/Ang**3
> siesta:   8.67254766  8.73987328  kBar
> (Free)E+ p_basis*V_orbitals  =  -30859.763440
> (Free)Eharris+ p_basis*V_orbitals  =  -30859.763491
>  spin moment: S , {S} =0.0   0.0   0.0   0.0
>
> siesta: Electric dipole (a.u.)  =0.000.0432460.00
> siesta: Electric dipole (Debye) =0.010.1099190.00
>
> Hirshfeld Net Atomic Populations:
> Atom #Qatom  Species
>  10.149  B
>  20.149  B
>  30.149  B
>  40.149  B
>  5   -0.149  N
> ...
>155   -0.149  N
>156   -0.149  N
>1570.149  B
>1580.149  B
>1590.149  B
>1600.149  B
>
> Voronoi Net Atomic Populations:
> Atom #Qatom  Species
>  10.167  B
>  20.167  B
>  30.167  B
>  40.167  B
>  5   -0.168  N
> ...
>155   -0.168  N
>156   -0.168  N
>1570.168  B
>1580.168  B
>1590.168  B
>1600.168  B
> Bader Analysis core-charge setup. Radii (standard, H):  1.000 0.600
>
> dhscf: Vacuum level (max, mean) =   -0.038479   -0.112800 eV
>
> siesta: LDOS info
> siesta: E1 -- E2 [eV]:  -20.000 --0.000
>
> Hirshfeld Net Atomic Populations:
> Atom #Qatom  Species
>  10.227  B
>  20.227  B
>  30.227  B
>  40.227  B
>  50.888  N
> ...
>1550.886  N
>1560.885  N
>1570.227  B
>1580.227  B
>1590.227  B
>1600.227  B
>
> Voronoi Net Atomic Populations:
> Atom #Qatom  Species
>  10.119  B
>  20.120  B
>  30.120  B
>  40.120  B
>  50.996  N
> ...
>1550.993  N
>1560.993  N
>1570.119  B
>1580.119  B
>1590.119  B
>1600.119  B
>
> >> End of run:  10-NOV-2021  11:48:50
> Job completed
>
>
> []'s,
>
> Camps
>
>
> On Thu, Jun 16, 2022 at 5:02 PM Alberto Garcia  wrote:
>>
>> Hi,
>>
>> I cannot reproduce your results. Which version of Siesta are you using? Can 
>> you show your output?
>>
>> The expected behavior is something like this (obtained with the 4.1 branch 
>> version):
>>
>> [...]
>> siesta: Electric dipole (a.u.)  =   -0.000.558297   -0.00
>> siesta: Electric dipole (Debye) =   -0.001.419050   -0.00
>>
>> Hirshfeld Net Atomic Populations:
>> Atom #Qatom  Species
>>  1   -0.224  O
>>  20.113  H
>>  30.113  H
>>
>> Voronoi Net Atomic Populations:
>> Atom #Qatom  Species
>>  1   -0.164  O
>>  20.082  H
>>  30.082  H
>> Bader Analysis core-charge setup. Radii (standard, H):  1.000 0.600
>>
>> 

Re: [SIESTA-L] << two set of charge calculations >>

[I. Camps]

Hello Alberto,

Here it is the info about the SIESTA version:










*Siesta Version  : siesta-max-R3--710-676-597Architecture:
unknownCompiler version: ifort (IFORT) 19.1.1.217 20200306Compiler flags  :
mpifort -fPIC -O2 -march=core-avx2 -axCore-AVX512 -fp-model precisePP flags
   : -DFC_HAVE_ABORT -DF2003 -DMPI -DCDF -DNCDF -DNCDF_4
-DNCDF_PARALLEL
-I/cvmfs//soft.computecanada.ca/easybuild/software/2020/avx2/MPI/intel2020/openmpi4/netcdf-fortran-mpi/4.5.2/include
Libraries
  : libncdf.a libfdict.a -Wl,-Bstatic -Wl,--start-group
-lmkl_scalapack_lp64 -lmkkl_blacs_openmpi_lp64 -lmkl_intel_lp64
-lmkl_sequential -lmkl_core -Wl,--end-group -Wl,-Bdynamic -lnetcdffPARALLEL
versionNetCDF supportNetCDF-4 supportNetCDF-4 MPI-IO support*

And here is the output section:






































































































*siesta: Final energy (eV):siesta:  Band Struct. =   -8272.290139siesta:
Kinetic =   19960.524774siesta:   Hartree =  151423.860682siesta:
Eldau   =   0.00siesta:   Eso =   0.00siesta:
 Ext. field =   0.00siesta:   Enegf   =   0.00siesta:
Exch.-corr. =  -11180.064205siesta:  Ion-electron = -320401.282309siesta:
Ion-ion =  129282.468462siesta:   Ekinion =   0.00siesta:
  Total =  -30914.492596siesta: Fermi =  -4.212218siesta:
Stress tensor (static) (eV/Ang**3):siesta: 0.0001260.00
-0.00siesta: 0.000.000101   -0.49siesta:-0.00
-0.49   -0.016465siesta: Cell volume =   7672.635004 Ang**3siesta:
Pressure (static):siesta:SolidMolecule
 Unitssiesta:   0.5895  0.5941  Ry/Bohr**3siesta:
0.00541292  0.00545494  eV/Ang**3siesta:
8.67254766  8.73987328  kBar(Free)E+ p_basis*V_orbitals  =
 -30859.763440(Free)Eharris+ p_basis*V_orbitals  =  -30859.763491
 spin moment: S , {S} =0.0   0.0   0.0   0.0siesta:
Electric dipole (a.u.)  =0.000.0432460.00siesta:
Electric dipole (Debye) =0.010.1099190.00Hirshfeld Net
Atomic Populations:Atom #Qatom  Species 10.149  B
20.149  B   30.149  B
  40.149  B   5   -0.149  N
 ...   155   -0.149  N 156   -0.149  N
1570.149  B 1580.149  B
 1590.149  B 1600.149  B
 Voronoi Net Atomic Populations:Atom #Qatom  Species 10.167  B
  20.167  B   30.167  B
  40.167  B   5   -0.168  N
 ...   155   -0.168  N 156   -0.168  N
1570.168  B 1580.168  B
1590.168  B 1600.168  B
 Bader Analysis core-charge setup. Radii (standard, H):  1.000 0.600dhscf:
Vacuum level (max, mean) =   -0.038479   -0.112800 eVsiesta: LDOS
infosiesta: E1 -- E2 [eV]:  -20.000 --0.000Hirshfeld Net Atomic
Populations:Atom #Qatom  Species 10.227  B
 20.227  B   30.227  B   4
   0.227  B   50.888  N  ...   155
   0.886  N 1560.885  N 157
 0.227  B 1580.227  B 159
 0.227  B 1600.227  B  Voronoi Net
Atomic Populations:Atom #Qatom  Species 10.119  B
20.120  B   30.120  B
  40.120  B   50.996  N
 ...   1550.993  N 1560.993  N
1570.119  B 1580.119  B
 1590.119  B 1600.119  B   >>
End of run:  10-NOV-2021  11:48:50Job completed*


[]'s,

Camps


On Thu, Jun 16, 2022 at 5:02 PM Alberto Garcia  wrote:

> Hi,
>
> I cannot reproduce your results. Which version of Siesta are you using?
> Can you show your output?
>
> The expected behavior is something like this (obtained with the 4.1 branch
> version):
>
> [...]
> siesta: Electric dipole (a.u.)  =   -0.000.558297   -0.00
> siesta: Electric dipole (Debye) =   -0.001.419050   -0.00
>
> Hirshfeld Net Atomic Populations:
> Atom #Qatom  Species
>  1   -0.224  O
>  20.113  H
>  30.113  H
>
> Voronoi Net Atomic Populations:
> Atom #Qatom  Species
>  1   -0.164  O
>  20.082  H
>  30.082  H
> Bader Analysis core-charge setup. Radii (standard, H):  1.000 0.600
>
> dhscf: Vacuum level (max, mean) =0.636991   -0.068255 eV
>
> cite: Please see "h2o.bib" for an exhaustive BiBTeX file.
> [...]
>
> in which one gets two 

Re: [SIESTA-L] << two set of charge calculations >>

[Alberto Garcia]

Hi,

I cannot reproduce your results. Which version of Siesta are you using? Can you 
show your output?

The expected behavior is something like this (obtained with the 4.1 branch 
version):

[...]
siesta: Electric dipole (a.u.)  =   -0.000.558297   -0.00
siesta: Electric dipole (Debye) =   -0.001.419050   -0.00

Hirshfeld Net Atomic Populations:
Atom #Qatom  Species
 1   -0.224  O
 20.113  H
 30.113  H

Voronoi Net Atomic Populations:
Atom #Qatom  Species
 1   -0.164  O
 20.082  H
 30.082  H
Bader Analysis core-charge setup. Radii (standard, H):  1.000 0.600

dhscf: Vacuum level (max, mean) =0.636991   -0.068255 eV

cite: Please see "h2o.bib" for an exhaustive BiBTeX file.
[...]

in which one gets two blocks, one for Voronoi and another one for Hirshfeld 
populations.

  Alberto


- El 14 de Junio de 2022, a las 22:18, I. Camps ica...@gmail.com escribió:

| Hello,
| 
| I set my input to calculate and export the charges using Voronoi, Bader and
| Hirshfeld approaches.
| 
| My output has at the end two sets, one after the energy decomposition/final
| energy/etc. section, and then after some info about Bader/Vacuum level/LDOS
| info.
| 
| Both sets return different charges.
| 
| My questions are:
| - Why two sets of charges?
| - Which one is the "good" one?
| 
| []'s,
| 
| Camps
| 
| 
| --
| SIESTA is supported by the Spanish Research Agency (AEI) and by the European
| H2020 MaX Centre of Excellence (http://www.max-centre.eu/)

-- 
SIESTA is supported by the Spanish Research Agency (AEI) and by the European 
H2020 MaX Centre of Excellence (http://www.max-centre.eu/)