Methodology for calculating J-couplings

From NWChem

Viewed 1843 times, With a total of 4 Posts
Jump to: navigation, search

Clicked A Few Times
Threads 4
Posts 6
Hi, I just wanted to confirm whether the implementation of indirect spin-spin couplings (J-couplings) is based on CPDFT in nwchem. If so, the methodology should be the same as in Gaussian and I should get identical results. However, when I calculate the J-coupling constant for H2, the FC contribution comes out slightly different in Gaussian compared to nwchem. I'm using the following input:

start h2_nmr
geometry units angstroms
 H       3.18205 0.0 0.0
H 3.78205 0.0 0.0
 H library cc-pvdz
 xc pw91lda
grid lebedev
tolerances tight
 spinspin 1 1 2
task dft property

Even though the implementations differ a little (nwchem has advanced support for relativistic calculations), I should still get the same results for light elements. Could I be missing something essential in the input?
Edited On 6:25:12 AM PDT - Fri, Apr 22nd 2016 by Raul l

  • Niri Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
    Send PM
Forum Regular
Threads 4
Posts 209
Hi Raul,

Our implementation appears to be working correctly.
Can you try comparing with pure Hartree-Fock ?



Just Got Here
Threads 0
Posts 2
spin spin coupling
Dear NWChem users

I have managed to compute spin spin couplings (and in good agreement with the experiment!) when the atom pair under investigation is H-H for several standard neutral organic molecules. Just an open shell calculation and using the spinspin keyword in the property section as described in Documentation :

SPINSPIN [<integer> number_of_pairs <integer> pair_list]

However, when i tried a C-H coupling instead the calculation always stops during the iterative solution of linear equations and with no error....

I know that C12, the most abundant, is not NMR active... so i define isotope C13 by specifying the keyword mass, followed by 13 after the line of the Cartesian coordinate of this Carbon atom under study.

Same outcome for C12 and C13. The calculation always stops during the itterative solution of linear equations and with no error....

Has anyone been able to run successfully a C-H coupling prediction of any standard organic neutral organic molecule? If so, please share the input how to do so.

Thank you very much in advance


Forum Vet
Threads 10
Posts 1624
Dear Inaki
Could you please provide an input file reproducing the bug you are reporting?

Just Got Here
Threads 0
Posts 2
Thank you very much for answering me.

Firstly, i am using version 6.5. I have requested to our system admin to update the software but still pending...

It may be fixed in a newer version. I cannot test that, sorry.

After lots of trial and error i noticed that i can run the coupling calculation but only 1 by 1. In other words, i can not compute 2 C-H couplings in the same calculation.

Please find below the input of the calculation that failed:

scratch_dir /scratch

geometry units angstroms noautosym noautoz
 C        0.04907        0.82213       -0.65775 mass 13.003355
C -1.15981 1.28999 -1.16070
H -1.20486 2.13123 -1.84724
N -2.30970 0.69954 -0.80166
C -2.24978 -0.34780 0.04068
H -3.19802 -0.80119 0.31094
N -1.15919 -0.90083 0.59434
C -0.04509 -0.28180 0.22092 mass 13.003355
N 1.37198 1.20191 -0.81802
N 1.24761 -0.57303 0.58134
C 1.66725 -1.61610 1.49685
H 1.12804 -2.53639 1.26672
H 1.44762 -1.32294 2.52807
H 2.73910 -1.77789 1.37965
C 2.03623 0.35215 -0.07465
H 3.11035 0.33755 0.05285

* library 6-31G**


spinspin 2 16 1 16 8
task scf property

Please let me know if you get the same unfinished run without error.

If you compute only 1 C-H coupling, (either spinspin 1 16 1 or spinspin 1 16 8) the calculation is fine

Many thanks in advance


Forum >> NWChem's corner >> NWChem functionality

Who's here now Members 0 Guests 1 Bots/Crawler 0

AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC