Energy per unit cell in NWChem differs for different supercells

54 views
Skip to first unread message

Dee Mee

unread,
Feb 10, 2021, 8:53:06 AMFeb 10
to NWChem Forum

I try to compute the unit cell energy of the Ti2CO2 MXene with NWchem.

Computing the 1x1 supercell with 4x4x1 K-points grid I got -1875 Ha.

Computing the 2x2 supercell with 2x2x1 K-points grid I got -7502 Ha, that means -1875.5 Ha per unit cell.

I understand, that for useful results convergence should be reached with K-points and the basis set. But here I just want to understand the source of the 0.5 Ha discrepancy. Previously I thought, that NxM supercell is equivalent to NxM K-points grid. Apparently it is not true. What is the source of such a significant difference?

Inputs for both cases are attached.

I also have asked this question in StackExchange Materials Modelling group here. But for a while there are no decisive answers.

supercell-test.tar.gz

Bylaska, Eric J

unread,
Feb 10, 2021, 1:14:44 PMFeb 10
to nwchem...@googlegroups.com

 

It looks like you’re trying to use the Gaussian DFT code in NWChem.  This code does not contain periodic boundary conditions.  Only the plane-wave band code in the nwpw module can do periodic dft calculations with k-points.

--
You received this message because you are subscribed to the Google Groups "NWChem Forum" group.
To unsubscribe from this group and stop receiving emails from it, send an email to nwchem-forum...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/nwchem-forum/ac16b13b-2f70-4636-88c3-02e157742d4an%40googlegroups.com.

Dee Mee

unread,
Feb 12, 2021, 2:43:22 AMFeb 12
to NWChem Forum
Thank you, Eric! You helped a lot!

Indeed, NWchem seems to supports k-points only in its plane-wave module. I should have noticed that the monkhorst-pack option is in the nwpw section and it is suspicious.
But I still think, that periodic boundary conditions are supported in the Gaussian DFT code. At least a see difference in energy between system crystal and without this section. Please, correct me if I am wrong.

Edoardo Aprà

unread,
Feb 12, 2021, 1:18:50 PMFeb 12
to NWChem Forum
No, periodic boundary conditions are not supported in the Gaussian DFT code.
I have opened a github issue about the code not stopping when these conditions appear in the input file

-

unread,
Feb 12, 2021, 10:34:10 PMFeb 12
to nwchem...@googlegroups.com
Is it just about 3D PBC (system crystal)? Or all periodic systems are not implemented for Gaussian DFT? The documentation here tells, that 2D (surface) and 1D (polymer) symmetries are "not currently available with NWPW module".  But where are they available?

You received this message because you are subscribed to a topic in the Google Groups "NWChem Forum" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/nwchem-forum/CYQeUO-zZv0/unsubscribe.
To unsubscribe from this group and all its topics, send an email to nwchem-forum...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/nwchem-forum/a11e0685-5dab-46ac-8e6e-94256b18106cn%40googlegroups.com.

Dee Mee

unread,
Feb 12, 2021, 10:34:18 PMFeb 12
to NWChem Forum
Is it just about 3D PBC (system crystal)? Or all periodic systems are not implemented for Gaussian DFT? The documentation here tells, that 2D (surface) and 1D (polymer) symmetries are "not currently available with NWPW module".  But where are they available?

Eric Bylaska

unread,
Feb 12, 2021, 10:36:48 PMFeb 12
to nwchem...@googlegroups.com
The plane wave code automatically places 2d and 1d periodic in a padded 3D cell


From: nwchem...@googlegroups.com <nwchem...@googlegroups.com> on behalf of Dee Mee <maz...@gmail.com>
Sent: Friday, February 12, 2021 7:23:59 PM
To: NWChem Forum <nwchem...@googlegroups.com>
Subject: Re: [nwchem-forum] Energy per unit cell in NWChem differs for different supercells
 

Dee Mee

unread,
Feb 12, 2021, 11:00:55 PMFeb 12
to NWChem Forum
So it means, that for Gaussian basis there are neither 3D, nor 2D, nor 1D periodic boundaries. Right?

Edoardo Aprà

unread,
Feb 12, 2021, 11:52:57 PMFeb 12
to NWChem Forum
Reply all
Reply to author
Forward
0 new messages