Effective core potential (ECP): for atoms or for elements?


Just Got Here
Hello Admins,

I have a Copper dimer, and I want to use ECP for one atom (to replace 10 core electrons) and all-electron basis for the other one.
So I tagged the two coppers as "Cu1" and "Cu2" as follows. However, it turned out that nwchem 6.5 replaced 10 electrons for both coppers, and I obtained only 38 electrons (expected 29+19=48). It seems that the program applied the ECP for all atoms of the same element, though different tags are used.

Is there any workaround to treat these atoms differently? Thanks.



Input


geometry
Cu1    0.0 0.0 0.0
Cu2 0.0 0.0 2.0
symmetry C1
end

ecp
Cu2 library Cu "Stuttgart RSC 1997 ECP"
end

basis "ao basis" spherical PRINT
Cu1 library Cu "Def2-TZVP" 
Cu2 library Cu "Stuttgart RSC 1997 ECP"
end

dft
direct
xc b3lyp
print "final vectors analysis"
end
task dft



Output


         No. of atoms     :     2
No. of electrons : 38
Alpha electrons : 19
Beta electrons : 19


P.S.
I had tried to explicitly set ECP also for "Cu1" and set all coefficients as zero, but failed with

  int_ecp_build_ecp_ptrs:too many coefficients                 911

PS. 2
I had also tried to treat "Cu1" as a ghost atom with charge and mass specified. It worked ok for the SCF task. But when I turned to call DFT task, I met with the error
  forrtl: severe (174): SIGSEGV, segmentation fault occurred


Forum >> NWChem's corner >> General Topics