[DFTB-Plus-User] Re: Failure in diagonalisation routine ZHEEVR for Li-S systems

Tue Nov 8 15:41:40 CET 2022

Dear Dr. Arabi,

I got confused when I was replying to the message below. I think I should
have edited the subject line before sending it out. Sorry for the
confusion. Here is the new reply with the correct subject line.

Best regards,
Qisheng Wu

On Tue, Nov 8, 2022 at 9:37 AM <
dftb-plus-user-request at mailman.zfn.uni-bremen.de> wrote:

> 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: DFTB-Plus-User Digest, Vol 99, Issue 6 (Wu, Qisheng)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 8 Nov 2022 09:37:01 -0500
> From: "Wu, Qisheng" <qisheng_wu at brown.edu>
> To: dftb-plus-user at mailman.zfn.uni-bremen.de
> Subject: Re: [DFTB-Plus-User] DFTB-Plus-User Digest, Vol 99, Issue 6
> Message-ID:
>         <
> CAPhmeq2ie6R+Ux+nLPQbiPBkW872Cf-Jc86QJvPkDepJioJ7xg at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Dear Dr. Bálint Aradi,
>
> Many thanks for the quick reply. You are right that the S is supposed to be
> with "d" as the maximal angular shell in the mio-set. I made the change and
> re-run the geometry optimizations. Unfortunately, this does not solve the
> problem. For the Li2S molecule, It turns out that it makes very little
> difference in terms of Li-S bond length (2.005 angstroms vs 2.002
> angstroms). But of course, "d" still should be used anyway.
>
> I have also looked at the input file generated within MS and attached it to
> this email (Li2S_tetra.hsd). The version of DFTB+ in MS 2020 might be
> different (I do not know what version it is), and there are some keywords
> that I did not previously use in my input. The Driver section is also
> empty, whose possible reason is that MS uses its own driver to conduct
> geometry optimizations (it's just my guess). After tailoring my own inputs
> based on what I found in the input from MS (dftb_in.hsd as attached), I was
> still not able to solve the problem.
>
> What I do not understand is that the calculation works for the Li2S
> molecule but not for the Li2S bulk. It's really puzzling. I hope this issue
> could be fixed in the end.
>
> By the way, I have tested using different versions of DFTB+ (e.g., 17.1),
> and the same error showed up for geometry optimization of the Li2S bulk,
> but not for the Li2S molecule.
>
> Thanks again.
>
> Best regards,
> Qisheng Wu
>
>
> On Tue, Nov 8, 2022 at 6:00 AM <
> dftb-plus-user-request at mailman.zfn.uni-bremen.de> wrote:
>
> > 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: Failure in diagonalisation routine ZHEEVR for Li-S
> >       systems (Bálint Aradi)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Mon, 7 Nov 2022 16:12:28 +0100
> > From: Bálint Aradi <aradi at uni-bremen.de>
> > To: dftb-plus-user at mailman.zfn.uni-bremen.de
> > Subject: Re: [DFTB-Plus-User] Failure in diagonalisation routine
> >         ZHEEVR for Li-S systems
> > Message-ID: <8d8e2f4f-6368-a719-f281-d263b107c054 at uni-bremen.de>
> > Content-Type: text/plain; charset="utf-8"; Format="flowed"
> >
> > Dear Qisheng Wu,
> >
> > The DFTB+ in MS is identical to the open source version, therefore, I am
> > a little bit puzzled, why you get different behavior for the bulk
> > calculation.
> >
> > Check, whether the dftb+ input file generated by MS is equivalent to
> > yours. My guess would be the maximal angular momenta. Actually, in the
> > mio-set the S is supposed to be with "d" as maximal angular shell (and
> > you use "p" instead). Also probably MS uses only "s" for the Li and not
> > "p" as in your input?
> >
> > Best regards,
> >
> > Bálint
> >
> >
> > On 07.11.22 15:47, Wu, Qisheng wrote:
> > > Dear DFTB+ developers,
> > >
> > > I am currently working on a Li-S system and fitted Li-Li and L-S
> > > parameters as an extension of the /mio/ set using the DFTB+ module and
> > > parameterization tools as implemented in the Materials Studio software.
> > > S-S is directly taken from the /mio/ set. To verify the parameters, I
> > > first conducted geometry optimizations of Li2S bulk and Li2S molecule,
> > > which ran very well through the DFTB+ module in the Materials Studio
> > > software.
> > >
> > > Then I switched to DFTB+ 22.1 code (precompiled executable for x86_64),
> > > which I will need to use for large-scale simulations. The geometry
> > > optimization of the Li2S molecule worked very well, but that of the L2S
> > > bulk failed, with the error information:
> > >
> > > "-> Failure in diagonalisation routine ZHEEVR, unable to complete
> > > Cholesky factorization of B     85"
> > >
> > > I checked some previous messages in this email list, and am aware that
> > > this could be due to unphysical geometry. But I do not think this is
> the
> > > case for my simulation since it worked well at least for the Li2S
> > molecule.
> > >
> > > I appreciate it if anyone can help address this issue, which will be
> > > very helpful for me to continue my research project. I have put all the
> > > inputs (dftb_in.hsd), outputs (out), and parameter files (*skf) in my
> > > google drive. Feel free to download that using the link below:
> > >
> > >
> >
> https://drive.google.com/file/d/1LL8GHzfI8pH4Nlr2ZJ8IN4q3xTVpQmy1/view?usp=share_link
> > <
> >
> https://drive.google.com/file/d/1LL8GHzfI8pH4Nlr2ZJ8IN4q3xTVpQmy1/view?usp=share_link
> > >
> > >
> > > Best regards,
> > > Qisheng Wu
> > >
> > > _______________________________________________
> > > 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/20221107/2202d405/attachment-0001.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 99, Issue 6
> > *********************************************
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20221108/9a800232/attachment.htm
> >
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: Li2S_tetra.hsd
> Type: application/octet-stream
> Size: 3842 bytes
> Desc: not available
> URL: <
> http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20221108/9a800232/attachment.obj
> >
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: dftb_in.hsd
> Type: application/octet-stream
> Size: 3997 bytes
> Desc: not available
> URL: <
> http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20221108/9a800232/attachment-0001.obj
> >
>
> ------------------------------
>
> 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 99, Issue 7
> *********************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20221108/e820d739/attachment.htm>

[DFTB-Plus-User] ​Re: Failure in diagonalisation routine ZHEEVR for Li-S systems

[DFTB-Plus-User] Re: Failure in diagonalisation routine ZHEEVR for Li-S systems