CCSD(T) fails

175 views
Skip to first unread message

Claire

unread,
Sep 8, 2020, 8:17:53 PM9/8/20
to NWChem Forum
Hello,
I am trying to run a doublet CCSD(T) calculation on a cation.

CCSD(T) fails when it tries to transform the integrals, and I am not sure why.  I was running it on 1 node, 44 cores, just as a test.


Attached is the output, with the input echo'd.
out.txt

Kowalski, Karol

unread,
Sep 8, 2020, 8:41:59 PM9/8/20
to nwchem...@googlegroups.com, clair...@gmail.com

Please use the following input (I have increased the “stack” memory allocation & added tilesize kwyword).

This is pretty big calculations with CCSDT – try to use the symmetry (this is atomic system  - you may want to use d2h symmetry –

it should significantly reduce memory requirements associated with storing triply excited amplitudes) – please try:

geometry units angstrom

symmetry d2h

Yb     0.00000      0.00000     0.00000

end

 

Best,

K

 

 

 

start uo2_sodft

echo

memory stack 1600 mb heap 100 mb global 8000 mb

charge 3

geometry  units angstrom

symmetry d2h

Yb     0.00000      0.00000     0.00000

end

BASIS "ao basis" PRINT

#BASIS SET: (23s,16p,12d,6f) -> [18s,12p,9d,3f]

Yb    S

1358076.401925000            0.0496959645

603589.5119670000          -0.0108376169

268262.0053190000           0.1124431010

119227.5579190000           0.0897513516

  52990.0257420000           0.3207996474

  23551.1225520000           0.5578901152

Yb    S

  10467.1655790000           1.0000000000

Yb    S

   4652.0735910000           1.0000000000

Yb    S

   2067.5882620000           1.0000000000

Yb    S

    918.9281170000           1.0000000000

Yb    S

    408.4124960000           1.0000000000

Yb    S

    181.5166650000           1.0000000000

Yb    S

     80.6740730000           1.0000000000

Yb    S

     35.8551440000           1.0000000000

Yb    S

     15.9356190000           1.0000000000

Yb    S

      7.0824980000           1.0000000000

Yb    S

      3.1477770000           1.0000000000

Yb    S

      1.3990120000           1.0000000000

Yb    S

      0.6217830000           1.0000000000

Yb    S

      0.2763480000           1.0000000000

Yb    S

      0.1228210000           1.0000000000

Yb    S

      0.0545870000           1.0000000000

Yb    S

      0.0242610000           1.0000000000

Yb    P

  18741.9577990000           0.0100228710

   7496.7831200000           0.0158264627

   2998.7132480000           0.0766658728

   1199.4852990000           0.2489842474

    479.7941200000           0.7465576446

Yb    P

    191.9176480000           1.0000000000

Yb    P

     76.7670590000           1.0000000000

Yb    P

     30.7068240000           1.0000000000

Yb    P

     12.2827290000           1.0000000000

Yb    P

      4.9130920000           1.0000000000

Yb    P

      1.9652370000           1.0000000000

Yb    P

      0.7860950000           1.0000000000

Yb    P

      0.3144380000           1.0000000000

Yb    P

      0.1257750000           1.0000000000

Yb    P

      0.0503100000           1.0000000000

Yb    P

      0.0201240000           1.0000000000

Yb    D

   1379.8152340000           0.0075831557

    501.7509940000           0.0443885329

    182.4549070000           0.2483587502

     66.3472390000           0.8033434192

Yb    D

     24.1262690000           1.0000000000

Yb    D

      8.7731890000           1.0000000000

Yb    D

      3.1902500000           1.0000000000

Yb    D

      1.1600910000           1.0000000000

Yb    D

      0.4218510000           1.0000000000

Yb    D

      0.1534000000           1.0000000000

Yb    D

      0.0557820000           1.0000000000

Yb    D

      0.0202840000           1.0000000000

Yb    F

     41.1409390000           0.1076071148

     13.7136460000           0.3540592739

      4.5712150000           0.5193050550

      1.5237380000           0.3584834866

Yb    F

      0.5079130000           1.0000000000

Yb    F

      0.1693040000           1.0000000000

END

SCF

  DOUBLET

  UHF

END

TCE

  CCSDT

   TILESIZE 15

END

TASK TCE ENERGY

--
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/b9ff28f2-3049-46ae-b3ee-3590d7e10224n%40googlegroups.com.

Kowalski, Karol

unread,
Sep 8, 2020, 8:45:19 PM9/8/20
to nwchem...@googlegroups.com, clair...@gmail.com

By the way, the input you have attached is for the  CCSDT run.

If you want to run CCSD(T) please use the “ccsd(t)” keyword in the input.

(please also add “tilesize 15” to the TCE group).

 

Best,

K

 

 

From: <nwchem...@googlegroups.com> on behalf of Claire <clair...@gmail.com>
Reply-To: "nwchem...@googlegroups.com" <nwchem...@googlegroups.com>
Date: Tuesday, September 8, 2020 at 5:18 PM
To: NWChem Forum <nwchem...@googlegroups.com>
Subject: [nwchem-forum] CCSD(T) fails

 

Hello,

--

Claire

unread,
Sep 16, 2020, 11:06:01 AM9/16/20
to NWChem Forum
Thank you!  This works and I have changed it to CCSD(T). 

However I am finding, with eom-CCSD(T), the application exits prematurely (before the wall time would kill the job).  Is there a page which lists all the exit codes for nwchem?

Here is the error termination message:
 Iteration 146 using  371 trial vectors
Application 19744640 network throttled: 8 nodes throttled, 00:12:47 node-seconds
Application 19744640 exit codes: 134
Application 19744640 exit signals: Killed
Application 19744640 resources: utime ~40563s, stime ~2088s, Rss ~2450916, inblocks ~0, outblocks ~1141696

Edoardo Aprà

unread,
Sep 16, 2020, 12:53:36 PM9/16/20
to NWChem Forum
134 is a system exit code
128+6=134
6=SIGABRT
Are you running on a Cray system? If this is the case, This is likely to due the jobs running out of the memory. Could you post the complete output/error files?

Claire

unread,
Sep 16, 2020, 1:34:48 PM9/16/20
to NWChem Forum
Ah yes, this is the problem then.  I am running it on Cray, 8 nodes of 44 cores/each with memory:
memory stack 1600 mb heap 100 mb global 8000 mb

Attached is output.
out.txt

Edoardo Aprà

unread,
Sep 16, 2020, 1:56:47 PM9/16/20
to NWChem Forum
How much physical memory is present on each node?

jeff.science

unread,
Sep 16, 2020, 7:29:54 PM9/16/20
to NWChem Forum
"Iteration 146 using  371 trial vectors" seems like a bad sign, both in terms of EOMCCSD convergence and memory usage.  371 trial vectors where each vector is 20 MB (I think "t2 file size = 2585459" is given in double) is going to consume a lot of memory.

I wonder if Karol knows what the right convergence tricks are here.  Level shifting might help but I haven't run the code in a while.

Jeff

Claire

unread,
Sep 17, 2020, 6:09:36 PM9/17/20
to NWChem Forum
Each node has around 128 GB of memory.

Thanks!  I'll try level-shifting.

Edoardo Aprà

unread,
Sep 17, 2020, 6:38:41 PM9/17/20
to nwchem...@googlegroups.com
128*0.8/44=2.32GB max memory/core
Your input line has
8+1.6+0.1=9.7GB

You should change it to

memory stack 800 mb heap 100 mb global 1400 mb

If you want to stick to the large memory input line, you would need to use fewer processors/node.

128*0.8/9.7=10.55

That mean you should use just 10 processors/node, instead of 44

Kowalski, Karol

unread,
Sep 17, 2020, 8:39:43 PM9/17/20
to nwchem...@googlegroups.com

Yes, there are ways to improve convergence but I would need a copy of input to reproduce the problem.

Claire

unread,
Sep 18, 2020, 1:05:49 AM9/18/20
to NWChem Forum
Hello, attached is the input and pbs script.  Thank you for any insight!

nwchem-pbs.txt
input.txt
Reply all
Reply to author
Forward
0 new messages