Multiwfn forum

i.s.ger

Member

Registered: 2020-12-01

Posts: 55

Constructing of promolecule density

Dear Tian Lu

I'm trying to construct promolecule density of dimers using atom's densities.

I'm running the following set of commands:

./Multiwfn H0.5.wfn 6 31 0,0,-0.47243153 1 32 0,0,.94486306 1 1 0

It produces new.wfn which is incorrect since when I'm trying to integrate it using:

./Multiwfn new.wfn 100 4 1

it leads to SegFault. It comes from `noGUI/sub.f90:2610`

It seems to me that modifying of wfn by that way does not produce the correct wfn.
At least, the number of electrons should be 2 times more (occupation number).

There is an initial wfn and new.wfn:
H0.5.wfn- initial
new.wfn- wfn produced by MultiWFN

Best regards,
Igor

i.s.ger

Member

Registered: 2020-12-01

Posts: 55

Re: Constructing of promolecule density

Dear Tian Lu,

I understand how to fix that part of the MultiWFN code for being it to work correctly (at least as I want).

I found a mistake with variable usage: center_tmp is not declared but used. The correct variable is ncenter_tmp.
There is a needed patch (works both for ifort and gfortran):
0001-center_tmp-is-not-declared.patch.txt

Best regards,
Igor

i.s.ger

Member

Registered: 2020-12-01

Posts: 55

Re: Constructing of promolecule density

If someone will be interested in such duplication of system, this is the patch for doing that:

0002-Copy-fragments-and-MOs-while-dublication-is.patch.txt

It will produce repeated independent structures with independent molecular orbitals.

Best regards,
Igor

sobereva

Tian Lu (Multiwfn developer)

From: Beijing

Registered: 2017-09-11

Posts: 1,768

Website

Re: Constructing of promolecule density

Dear Igor,

I just updated Multiwfn on Multiwfn website. I significantly extended the subfunction 19 of main function 100, it is highly recommended to use this function to combine fragment wavefunctions to form a whole wavefunction.

The patch in #2 has been considered in this update. The problem you reported in #1 has also been fixed.

Best regards,

Tian

i.s.ger

Member

Registered: 2020-12-01

Posts: 55

Re: Constructing of promolecule density

Dear Tian Lu,

Thank you so much!

Unfortunately, using 100.19, that is not possible to construct promolecular density of different dimers without changing .wfn files.

I used that script for generating dimers with promolecular densities for comparing their densities with SCF-optimized densities (they are at H2p_*.wfn files).

for wfn in H2p_*.wfn; do rm H0.5.wfn_001_0300_0434 coord=`grep "(CENTRE 1)" $wfn | sed "s/ \+/ /g" | cut -d" " -f8` coord2=`echo "${coord} * 2" | bc -l` ADD=`echo "1000 10 48 6 31 0,0,-${coord} 1 32 0,0,${coord2} 1 1 -1 100 -5 1 ${wfn} 0" | tr " " "\n" | ./gMultiwfn H0.5.wfn 2>/dev/null | grep "Final result:" | sed "s/ \+/ /g" | cut -d" " -f4` echo "${wfn};${ADD}" done

Probably, in the general case, I need a reader of XYZ-files and atomwfn directory for constructing promolecular densities for given XYZ's.

Best regards,
Igor