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
|