[molpro-user] orbitals to domains attribution problem on df-lccsd(t)-f12 with interact=1

[Edrisse Chermak]

Dear Molpro developers & users,

I'm trying to run a df-lccsd(t)-f12 single point on a four molecules
system, by taking in account the excitations for intermolecular
interactions (interact=1). After experimenting several options of the
program, I always face two issues:

1) after localizing the orbitals  the code finds 5 molecules (the 5th is
an hydrogen atom).

  - I tried thus to localize either with Pipek-Mezey (2 or 3)
with/without the delete option, or using NPA, but this attribution remains.

2) after defining the molecules, the code defines the domains then stops
because at least one orbital domain belongs to two different molecules.

  - I tried thus to change the CHGFRAC parameter (THRBP in the manual)
from 0.965 to 0.995 when localizing with Pipek-Mezey. I also tried to
change the NPASEL parameter from 0.03 to 0.14, when localizing with NPA,
but still always one orbital domain at least belongs to two molecules
(including the 5th)

  - I also tried to merge the orbital domains, but in this case domains
spread over all the system and thus the code stops.

Please find attached outputs for those cases. I'm running Molpro 2012.1
Thanks in advance for reading and proposing ideas,
Edrisse
--

Edrisse Chermak
Post-Doctoral Fellow
Office : +966 2 8080705
Catalysis center - Modeling group, building 3, 4th floor
KAUST, Thuwal, KSA
http://kcc.kaust.edu.sa/Pages/Home.aspx


________________________________

This message and its contents including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PM_5molecules_orbitalproblem.out
Type: chemical/x-gulp
Size: 37747 bytes
Desc: not available
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20130724/71de5c39/attachment.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NPA_5molecules_orbitalproblem.out
Type: chemical/x-gulp
Size: 40001 bytes
Desc: not available
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20130724/71de5c39/attachment-0001.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NPA_5molecules_mergeproblem.out
Type: chemical/x-gulp
Size: 87550 bytes
Desc: not available
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20130724/71de5c39/attachment-0002.bin>