Electrophilicity of Aldehydes and Ketones (Reactor context)

I am using Reactor 5.7.0.  I obtained unexpected results while using a Borch reductive amination reaction with this selectivity rule:  -electrophilicity(ratom(1));  -nucleophilicity(ratom(2))     and a tolerance of 0.01.  The reaction is defined as shown below. 

Aryl aldehydes are more electrophilic than aryl ketones, yet this is not what is observed.

Expected:
CC(=O)C1=CC=CC(C=O)=C1.C1CCNCC1>>CC(=O)C1=CC=CC(CN2CCCCC2)=C1 |c:5,9,22,32,t:3,20|
Calculated:
CC(=O)C1=CC=CC(C=O)=C1.C1CCNCC1>>CC(=O)C1=CC=CC(CN2CCCCC2)=C1.CC(N1CCCCC1)C1=CC=CC(C=O)=C1 |c:5,9,22,32,45,49,t:3,20,43|

The Calculated result is clearly wrong.    I would like the electrophilicity calculation to be fixed so that it makes the correct predictions.

Interestingly, for a related example, the expected and calculated results were the same, namely:

CC(=O)C1=CC=CC=C1C1=CC(C=O)=CC=C1.NCC1=CC(=CC=C1)C1=C(N)C=CC=C1>>CC(=O)C1=C(C=CC=C1)C1=CC(CNCC2=CC(=CC=C2)C2=C(N)C=CC=C2)=CC=C1

 

I suppose, that currently only aromatic electropholicity is calculated on Huckel-basis, but my colleague will correct me if I am  wrong.
I would suggest to use something like this instead:

hCount(ratom(1))

where the carbonyl carbon is mapped with number 1. Not too sophisticated, but it is fast and always prefers aldehydes to ketones.

Hi,

Yes ,  the simple Huckel molecular orbital (SHMO) calculation can  not  make difference between the  "keto" and the "aldehyde" groups. Because only  delocalizable atoms are included in the scope of the SHMO calculation.The  "Methyl" side of the "keto" group and  the "H" atom of the "aldehyd" groups are ineffectives in the SHMO, this is why the calculated  electrophilicity of the two functional groups will be the same.

We are going to discuss inhouse,  what we can do  with this problem.

 

Jozsi