[DFTB-Plus-User] DFTB-Plus-User Digest, Vol 106, Issue 11
Saleh Md Abu
saleh.md.abu.314 at s.kyushu-u.ac.jp
Mon Jun 12 16:58:11 CEST 2023
Dear Bálint,
Thanks for your crystal-clear explanation. Now, I understand the issue of the dipole moment.
have a good day!
Thanks and best regards,
Saleh
________________________________
From: DFTB-Plus-User <dftb-plus-user-bounces at mailman.zfn.uni-bremen.de> on behalf of dftb-plus-user-request at mailman.zfn.uni-bremen.de <dftb-plus-user-request at mailman.zfn.uni-bremen.de>
Sent: Monday, June 12, 2023 9:55 PM
To: dftb-plus-user at mailman.zfn.uni-bremen.de <dftb-plus-user at mailman.zfn.uni-bremen.de>
Subject: DFTB-Plus-User Digest, Vol 106, Issue 11
Send DFTB-Plus-User mailing list submissions to
dftb-plus-user at mailman.zfn.uni-bremen.de
To subscribe or unsubscribe via the World Wide Web, visit
https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user
or, via email, send a message with subject or body 'help' to
dftb-plus-user-request at mailman.zfn.uni-bremen.de
You can reach the person managing the list at
dftb-plus-user-owner at mailman.zfn.uni-bremen.de
When replying, please edit your Subject line so it is more specific
than "Re: Contents of DFTB-Plus-User digest..."
Today's Topics:
1. Re: Printing of Dipole Moment (Bálint Aradi)
2. Re: Question about dptools (Bálint Aradi)
----------------------------------------------------------------------
Message: 1
Date: Mon, 12 Jun 2023 14:48:29 +0200
From: Bálint Aradi <aradi at uni-bremen.de>
To: dftb-plus-user at mailman.zfn.uni-bremen.de
Subject: Re: [DFTB-Plus-User] Printing of Dipole Moment
Message-ID: <bf59b3a7-fe1c-e085-75ee-2248ef6064ca at uni-bremen.de>
Content-Type: text/plain; charset="utf-8"; Format="flowed"
Dear Saleh,
in periodic system, the absolute value of the dipole is generally ill
defined, since you can shift the cell boundaries arbitrarily around. If
you think about a 1D example of periodic alternating charges, depending
how you lay your periodic boundaries, the dipole vector points points
either to the right or to the left (as you calculate if from the charges
within your unit cell):
... | - + | - + | - + | ...
... - | + - | + - | + ...
The warning should just ensure, that you do not use the printed dipole
vector for something, you are not supposed to.
Best regards,
Bálint
On 08.06.23 02:40, Saleh Md Abu wrote:
> Dear concern,
> I have recently updated DFTB+ 20.1 to DFTB+ 22.2. After that, I ran an
> NPT simulation of 50 water molecules and got a warning message;*Dipole
> moment is not defined absolutely!*. I ran the same simulation in the
> previous version but didn't see any warning message of dipole moment. I
> do not understand the meaning of this warning. I am afraid I have
> prepared poor systems, or is there something wrong with my systems?
>
> Can anyone help me in this regard?
> Thanks and best regards,
> Saleh
>
> *Input file as below*
> Geometry = genFormat {
> <<< "geometry.gen"
> }
>
> Driver = VelocityVerlet{
> TimeStep [fs] = 1.0
> Thermostat = NoseHoover {
> Temperature [Kelvin] = 300
> CouplingStrength [cm^-1] = 3200
> }
> Barostat = {
> Pressure = 3.44396E-09
> # Coupling
> Timescale = 10.0
> }
> Steps = 1000000
> MovedAtoms = 1:-1
> MDRestartFrequency = 100
> }
>
> Hamiltonian = DFTB {
> Scc = Yes
> MaxSCCIterations = 2000
> MaxAngularMomentum {
> O = "p"
> H = "s"
> # C = "p"
> # N = "p"
> # Cl = "d"
> # S = "d"
> # P = "d"
> }
> # Charge = 1.0
> SlaterKosterFiles = Type2FileNames { # File names with two atom
> type names
> Prefix = "~/3ob-3-1/" # Prefix before first type name
> Separator = "-" # Dash between type names
> Suffix = ".skf" # Suffix after second type name
> }
> ThirdOrderFull = Yes
> HubbardDerivs {
> # C = -0.1492
> H = -0.1857
> O = -0.1575
> # N = -0.1535
> # Cl = -0.0697
> # S = -0.11
> # P = -0.14
> }
> HCorrection = H5 {
> RScaling = 0.68
> WScaling = 0.23
> H5Scaling {
> O = 0.045
> # N = 0.0675
> }
> }
> Dispersion = DftD3 {
> Damping = ZeroDamping {
> sr6 = 1.25 #3OB
> alpha6 = 29.61 #3OB
> }
> s6 = 1.0
> s8 = 0.49 # 3OB
> HHRepulsion = Yes
> }
> SCCTolerance = 1e-4
>
> Mixer = Anderson {
> MixingParameter = 0.05
> Generations = 8
> }
>
> Solver = DivideAndConquer {}
>
> Filling = Fermi {
> Temperature [Kelvin] = 1000.0
> }
>
> KPointsAndWeights = {
> 0.000000000000E+000 0.000000000000E+000 0.000000000000E+000
> 1.00000000000
> }
>
> }
>
> Options {
> WriteChargesAsText = Yes
> }
>
>
> Analysis {
> CalculateForces = Yes
> }
>
> ParserOptions {
> ParserVersion = 8
> }
>
> #Parallel {
> # Groups = 4
> # UseOmpThreads = Yes
> #}
>
> *Part of the output as below*
> MD step: 10100
> Nose-Hoover chain variables
> x:
> 0.9709126675E+00 0.5095393583E+01 -0.1201993198E+02
> v:
> 0.1172380185E-03 -0.6143591866E-04 -0.9400160624E-03
> g:
> -0.5122282026E-06 0.7590791869E-06 -0.5381051028E-05
> Lattice vectors (A)
> 0.11511686E+02 0.00000000E+00 0.00000000E+00
> 0.00000000E+00 0.11511686E+02 0.00000000E+00
> 0.00000000E+00 0.00000000E+00 0.11511686E+02
> Volume: 0.102947E+05 au^3 0.152552E+04 A^3
> Pressure: -0.435908E-05 au -0.128249E+09 Pa
> Gibbs free energy: -203.6228907991 H -5540.8608 eV
> Gibbs free energy including KE: -203.4334648600 H -5535.7062 eV
> Potential Energy: -203.6229262537 H -5540.8617 eV
> MD Kinetic Energy: 0.1894259385 H 5.1545 eV
> Total MD Energy: -203.4335003152 H -5535.7072 eV
> MD Temperature: 0.0008475433 au 267.6327 K
> *Warning: dipole moment is not defined absolutely!*
> Dipole moment: -16.09138492 -13.84621712 3.59721433 au
> Dipole moment: -40.90022412 -35.19357631 9.14320756 Debye
> MD step: 10200
> Nose-Hoover chain variables
> x:
> 0.9443649828E+00 0.6825341898E+01 -0.1345300994E+02
> v:
> 0.1130087503E-03 0.5289183249E-03 -0.1920912684E-02
> g:
> -0.2748048030E-06 0.3238016101E-06 -0.5105070806E-05
> Lattice vectors (A)
> 0.11516625E+02 0.00000000E+00 0.00000000E+00
> 0.00000000E+00 0.11516625E+02 0.00000000E+00
> 0.00000000E+00 0.00000000E+00 0.11516625E+02
> Volume: 0.103079E+05 au^3 0.152748E+04 A^3
> Pressure: 0.149716E-06 au 0.440480E+07 Pa
> Gibbs free energy: -203.6199934071 H -5540.7819 eV
> Gibbs free energy including KE: -203.4232467089 H -5535.4282 eV
> Potential Energy: -203.6200289073 H -5540.7829 eV
> MD Kinetic Energy: 0.1967466983 H 5.3538 eV
> Total MD Energy: -203.4232822090 H -5535.4291 eV
> MD Temperature: 0.0008802984 au 277.9759 K
> Warning: dipole moment is not defined absolutely!
> Dipole moment: -9.66148980 -15.98281037 -0.30577511 au
> Dipole moment: -24.55705956 -40.62425512 -0.77720286 Debye
> MD step: 10300
> Nose-Hoover chain variables
> x:
> 0.8823886131E+00 0.8379095854E+01 -0.1387168258E+02
> v:
> 0.1058884319E-03 -0.9580387030E-03 0.9837814733E-03
> g:
> 0.7523473872E-07 -0.3729004755E-06 -0.4466987244E-05
> Lattice vectors (A)
> 0.11528007E+02 0.00000000E+00 0.00000000E+00
> 0.00000000E+00 0.11528007E+02 0.00000000E+00
> 0.00000000E+00 0.00000000E+00 0.11528007E+02
> Volume: 0.103385E+05 au^3 0.153201E+04 A^3
> Pressure: -0.152318E-04 au -0.448135E+09 Pa
> Gibbs free energy: -203.6052263574 H -5540.3801 eV
> Gibbs free energy including KE: -203.3982624316 H -5534.7483 eV
> Potential Energy: -203.6052619629 H -5540.3811 eV
> MD Kinetic Energy: 0.2069639236 H 5.6318 eV
> Total MD Energy: -203.3982980394 H -5534.7493 eV
> MD Temperature: 0.0009260131 au 292.4115 K
> *Warning: dipole moment is not defined absolutely!*
> Dipole moment: -12.68722536 -22.05729297 17.53312324 au
> Dipole moment: -32.24771287 -56.06405106 44.56475768 Debye
> MD step: 10400
> Nose-Hoover chain variables
> x:
> 0.8346498308E+00 0.6569701582E+01 -0.1188304703E+02
> v:
> -0.1442483628E-04 0.7331923441E-03 0.6808164405E-03
> g:
> 0.1603784843E-06 -0.5291815472E-05 -0.4847254387E-05
> Lattice vectors (A)
> 0.11532605E+02 0.00000000E+00 0.00000000E+00
> 0.00000000E+00 0.11532605E+02 0.00000000E+00
> 0.00000000E+00 0.00000000E+00 0.11532605E+02
> Volume: 0.103509E+05 au^3 0.153385E+04 A^3
> Pressure: 0.977452E-05 au 0.287576E+09 Pa
> Gibbs free energy: -203.5923194080 H -5540.0289 eV
> Gibbs free energy including KE: -203.3778966071 H -5534.1941 eV
> Potential Energy: -203.5923550562 H -5540.0299 eV
> MD Kinetic Energy: 0.2144228023 H 5.8347 eV
> Total MD Energy: -203.3779322538 H -5534.1951 eV
> MD Temperature: 0.0009593861 au 302.9498 K
> Warning: dipole moment is not defined absolutely!
> Dipole moment: -0.26787305 0.89270894 16.61393390 au
> Dipole moment: -0.68086544 2.26903997 42.22841124 Debye
> MD step: 10500
> Nose-Hoover chain variables
> x:
> 0.9051690473E+00 0.7504479737E+01 -0.1183789952E+02
> v:
> 0.3205883830E-04 -0.2075510106E-02 -0.9431590510E-03
> g:
> -0.8509819391E-07 -0.4925412607E-05 -0.1077083201E-05
> Lattice vectors (A)
> 0.11540175E+02 0.00000000E+00 0.00000000E+00
> 0.00000000E+00 0.11540175E+02 0.00000000E+00
> 0.00000000E+00 0.00000000E+00 0.11540175E+02
> Volume: 0.103713E+05 au^3 0.153687E+04 A^3
> Pressure: -0.468445E-06 au -0.137821E+08 Pa
> Gibbs free energy: -203.6140395857 H -5540.6199 eV
> Gibbs free energy including KE: -203.4081380876 H -5535.0171 eV
> Potential Energy: -203.6140753041 H -5540.6209 eV
> MD Kinetic Energy: 0.2059014980 H 5.6029 eV
> Total MD Energy: -203.4081738061 H -5535.0180 eV
> MD Temperature: 0.0009212595 au 290.9104 K
> Warning: dipole moment is not defined absolutely!
> Dipole moment: -4.06856578 -13.63431827 14.24332489 au
> Dipole moment: -10.34126354 -34.65498313 36.20292366 Debye
> MD step: 10600
> Nose-Hoover chain variables
> x:
> 0.9056278317E+00 0.6487463254E+01 -0.1048597708E+02
> v:
> 0.8879065586E-04 0.4594279028E-02 0.2246002129E-03
> g:
> -0.4339333373E-06 -0.1860775486E-05 0.1572257438E-04
> Lattice vectors (A)
> 0.11540304E+02 0.00000000E+00 0.00000000E+00
> 0.00000000E+00 0.11540304E+02 0.00000000E+00
> 0.00000000E+00 0.00000000E+00 0.11540304E+02
> Volume: 0.103717E+05 au^3 0.153692E+04 A^3
> Pressure: -0.262813E-05 au -0.773222E+08 Pa
> Gibbs free energy: -203.6008275147 H -5540.2604 eV
> Gibbs free energy including KE: -203.4105803559 H -5535.0835 eV
> Potential Energy: -203.6008632343 H -5540.2614 eV
> MD Kinetic Energy: 0.1902471584 H 5.1769 eV
> Total MD Energy: -203.4106160759 H -5535.0845 eV
> MD Temperature: 0.0008512177 au 268.7930 K
> Warning: dipole moment is not defined absolutely!
> Dipole moment: 7.22136767 -16.73774786 -7.02930762 au
> Dipole moment: 18.35488726 -42.54311496 -17.86671928 Debye
>
>
> _______________________________________________
> DFTB-Plus-User mailing list
> DFTB-Plus-User at mailman.zfn.uni-bremen.de
> https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user
--
Dr. Bálint Aradi
Bremen Center for Computational Materials Science, University of Bremen
http://www.bccms.uni-bremen.de/cms/people/b-aradi/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: OpenPGP_signature
Type: application/pgp-signature
Size: 840 bytes
Desc: OpenPGP digital signature
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20230612/b26d00af/attachment-0001.sig>
------------------------------
Message: 2
Date: Mon, 12 Jun 2023 14:55:11 +0200
From: Bálint Aradi <aradi at uni-bremen.de>
To: mehrdad shiri via DFTB-Plus-User
<dftb-plus-user at mailman.zfn.uni-bremen.de>
Subject: Re: [DFTB-Plus-User] Question about dptools
Message-ID: <e3f250e0-eeb3-3b98-8202-d2a589c1940c at uni-bremen.de>
Content-Type: text/plain; charset="utf-8"; Format="flowed"
Dear Mehrdad,
do you use the most recent version of dptools and dftb+? If yes, please
make a bug report on GitHub:
https://github.com/dftbplus/dftbplus/issues
Also, please note, that there is no support team for dftb+. We answer
user questions on an absolutely voluntarily basis by sparing some of our
research time for helping others with their dftb+ related problems. ;-)
Best regards,
Bálint
On 08.06.23 01:17, mehrdad shiri via DFTB-Plus-User wrote:
> Dear support team
> When I am using the dp_bands command from dp_tools module it always
> gives an error of " nkpt = len(eigvalarrays) / nspin ZeroDivisionError:
> integer division or modulo by zero"
>
> It seems that the parameters are not being extracted correctly from
> band.out file. how can I fix this issue?
>
> Best
>
> Mehrdad
>
> _______________________________________________
> DFTB-Plus-User mailing list
> DFTB-Plus-User at mailman.zfn.uni-bremen.de
> https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user
--
Dr. Bálint Aradi
Bremen Center for Computational Materials Science, University of Bremen
http://www.bccms.uni-bremen.de/cms/people/b-aradi/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: OpenPGP_signature
Type: application/pgp-signature
Size: 840 bytes
Desc: OpenPGP digital signature
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20230612/475f027b/attachment.sig>
------------------------------
Subject: Digest Footer
_______________________________________________
DFTB-Plus-User mailing list
DFTB-Plus-User at mailman.zfn.uni-bremen.de
https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user
------------------------------
End of DFTB-Plus-User Digest, Vol 106, Issue 11
***********************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20230612/5ca64afe/attachment.htm>
More information about the DFTB-Plus-User
mailing list