EPR g-factors / Summary
Hello,
Attached below are all the e-mail I received from the CCL regarding
ab initio evaluations of g-factors. Thanks to all who responded. In
addition to the programs mentioned below it seems that the program DALTON
(http://www.kjemi.uio.no/software/dalton/dalton.html#Capabilities)
also has
this capability.
Best regards,
Jan Jensen
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
Jan H. Jensen Assistant Professor
Department of Chemistry jan-jensen %! at !% uiowa.edu
University of Iowa Phone:(319) 335-1108
Iowa City, IA 52242 FAX: (319) 335-1270
http://www.uiowa.edu/~chemdept/faculty/jensen/
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
From: frisch %! at !% lorentzian.com (Mike Frisch)
Subject: Re: CCL:EPR g-factors
To: uunet!blue.weeg.uiowa.edu!jhjensen %! at !% uunet.uu.net (Jan H. Jensen)
Date: Thu, 24 Sep 1998 22:27:37 -0400 (EDT)
MIME-Version: 1.0
Jan H. Jensen writes:
>
>
> Hello,
>
> I am interested in calculating anisotropic EPR g-factors using ab
> initio electronic structure theory. I would like to know what programs can
> calculate this property, and what kind of accuracy one can expect.
>
> Best regards, Jan Jensen
>
Gaussian 98 can do this. There are references in the G98 manual (and the
web site)
to papers of Barone's which make detailed comparisons with experiment.
Mike Frisch
From: Jussi Eloranta <eloranta %! at !% endor.chem.jyu.fi>
Subject: Re: CCL:EPR g-factors
To: jhjensen %! at !% blue.weeg.uiowa.edu (Jan H. Jensen)
Date: Fri, 25 Sep 1998 11:33:59 +0300 (EEST)
Dear Jan,
> I am interested in calculating anisotropic EPR g-factors using ab
>initio electronic structure theory. I would like to know what programs can
>calculate this property, and what kind of accuracy one can expect.
>
I have been looking for the same information but have not found any program
that can do it (there are some private codes that can - as you can see from
some publications in JCP, for example). The closest I have found is GAMESS-UK
in which the g-tensor code was once in but has been commented out for some
reason. I have been trying to get Dr. M. F. Guest to put it back but I
think there has not been enough requests for it (or he is extremely busy
or something). Currently gamess-uk can calculate isotropic and anisotropic
hfc for various wavefunctions (including MR-CI too!). BTW this was the only
program I found that could do even this for MR-CI wavefunctions!
You can find more info from: http://wserv1.dl.ac.uk/CFS/
(the program is commercial but very cheap, includes very good support,
very good manuals, etc.)
Best regards,
Jussi Eloranta
University of Jyvaskyla
Finland
Hi
Im just sticking my head into this area myself so I'l share what i've found.
there is a good general review on dft with some stuff on ESR at
http://129.43.1.11/tmmec/current/CURRENT_REVIEWS/oscar_ventura/oscar5.html#120
if you can get onto it -I've had some trouble over the past few days.
also
This is something that appeared on the list in about 1996 and refers to using
the gaussian program for computing anisotropic coupling constants.
its from Vicence Barone and his address is at the end of the item.
I have received several requests about the computation of ESR anisotropic
coupling constants by Gaussian94. Since this seems a general question I
directly post this message for the whole CCL community.
These constants are nothing else than field gradients computed with the
spin rather than the total electronic density and not including the nuclear
contribution. Furthermore the resulting tensor must be put in the zero-trace
form. In gaussian94 it is possible to force these options by setting in the
keyword list the following items:
PROP IOP(6/17=2,6/26=4)
Here follows an input and the relevant part of the output for a STO-3G
computation of H2NO
----------------------------------------------------------------------------
#UHF/PROP IOP(6/17=2) IOP(6/26=4)
H2NO
0 2
X
X 1 1.0
N 2 1.0 1 90.0
H 3 NH 2 90.0 1 THETA
H 3 NH 2 90.0 1 -THETA
O 3 NO 2 ALPHA 1 180.0
NH=1.0179
NO=1.2778
THETA=58.9235
ALPHA=110.0
-----------------------------------------------------------------------------
Fermi contact analysis (atomic units).
1
1 N .042923
2 H -.005038
3 H -.005038
4 O .098025
**********************************************************************
Electrostatic Properties Using The SCF Density
**********************************************************************
Warning! Using spin rather than total density!
--- Only the electronic contributions will be computed ---
-----------------------------------------------------
Center ---- Electric Field Gradient ----
XX YY ZZ
-----------------------------------------------------
1 Atom .119043 -.295145 -.363285
2 Atom .002506 .031422 .029384
3 Atom .002506 .031422 .029384
4 Atom 2.130337 -1.663300 -1.698867
-----------------------------------------------------
-----------------------------------------------------
Center ---- Electric Field Gradient ----
( tensor representation )
3XX-RR 3YY-RR 3ZZ-RR
-----------------------------------------------------
1 Atom .298838 -.115349 -.183489
2 Atom -.018598 .010318 .008280
3 Atom -.018598 .010318 .008280
4 Atom 2.540947 -1.252690 -1.288257
-----------------------------------------------------
these are the principal values of anisotropic coupling constants
-------------------------------------------------------------------------------
Since the directions of principal moments are often significant and transforma-
tion to more conventional units can be performed once for ever, I have modified
the links 601 and 602 of gaussian to obtain the following output for the same
input
-------------------------------------------------------------------------------
------------------------------------------------------------------------------
Isotropic Fermi Contact Couplings
------------------------------------------------------------------------------
Atom a.u. MegaHertz Gauss 10(-4) cm-1
1 N(14) .04292 13.86856 4.94865 4.62605
2 H -.00504 -22.51979 -8.03562 -7.51179
3 H -.00504 -22.51979 -8.03562 -7.51179
4 O(17) .09802 -59.42481 -21.20426 -19.82198
------------------------------------------------------------------------------
**********************************************************************
Electrostatic Properties Using The SCF Density
**********************************************************************
Warning! Using spin rather than total density!
--- Only the electronic contributions will be computed ---
Atomic Center 1 is at -.021142 .543231 .000000
Atomic Center 2 is at .158563 1.036966 .871810
Atomic Center 3 is at .158563 1.036966 -.871810
Atomic Center 4 is at -.021142 -.734569 .000000
-----------------------------------------------------
Center ---- Spin Dipole Couplings ----
3XX-RR 3YY-RR 3ZZ-RR
-----------------------------------------------------
1 Atom .298838 -.115349 -.183489
2 Atom -.018598 .010318 .008280
3 Atom -.018598 .010318 .008280
4 Atom 2.540947 -1.252690 -1.288257
-----------------------------------------------------
Center ---- Spin Dipole Couplings ----
XY XZ YZ
-----------------------------------------------------
1 Atom -.074772 .000000 .000000
2 Atom .004800 .007098 .035617
3 Atom .004800 -.007098 -.035617
4 Atom -.099769 .000000 .000000
------------------------------------------------------------------------------
Anisotropic Spin Dipole Couplings in Principal Axis System
------------------------------------------------------------------------------
Atom a.u. MegaHertz Gauss 10(-4) cm-1 Axes
Baa -.1835 -7.077 -2.525 -2.361 .0000 .0000
1.0000
Baa -.1835 -7.077 -2.525 -2.361 .0000 .0000
1.0000
1 N(14) Bbb -.1284 -4.953 -1.767 -1.652 .1724 .9850
.0000
Bcc .3119 12.030 4.293 4.013 .9850 -.1724
.0000
1/R**3 -.5394 -20.803 -7.423 -6.939
Baa -.0268 -14.276 -5.094 -4.762 -.2361 -.6571
.7158
2 H Bbb -.0193 -10.278 -3.667 -3.428 .9631 -.2560
.0828
Bcc .0460 24.554 8.762 8.190 .1288 .7090
.6933
1/R**3 .0633 33.780 12.053 11.268
Baa -.0268 -14.276 -5.094 -4.762 .2361 .6571
.7158
3 H Bbb -.0193 -10.278 -3.667 -3.428 .9631 -.2560
-.0828
Bcc .0460 24.554 8.762 8.190 .1288 .7090
-.6933
1/R**3 .0633 33.780 12.053 11.268
Baa -1.2883 93.221 33.264 31.095 .0000 .0000
1.0000
4 O(17) Bbb -1.2553 90.837 32.413 30.300 .0263 .9997
.0000
Bcc 2.5436 -184.059 -65.677 -61.395 .9997 -.0263
.0000
1/R**3 -1.2318 89.138 31.807 29.733
------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Vincenzo Barone |
Professor of Theoretical Chemistry |
Dipartimento di Chimica | tel. +39-81-5476503
Universita' Federico II | fax +39-81-5527771
via Mezzocannone 4 | e-mail ENZO %! at !%
CHEMNA.DICHI.UNINA.IT
I-80134 Napoli |
Italy
you might also want to read some of his papers i.e.
V. Barone and C. Adamo Theor. Chim. Acta 91 (1995) 129
.. .. Chem. Phys. Letters 224 (1994) 432
.. J. Chem Phys 101 (1994) 10666
or posibly
Leif A Eriksson Internationa Journal of Quantum Chemistry vol52 879-901
(1994)
The impresion I get -so far- is that good acuracy can be got for
anistropic constants
with most ab-initio methods. For isortopic ones you need good basis sets
with split
cores and DFT b3lyp seems to do particularily well. If not using DFT meed big
CC or CI calculations which mean V Small molecules only
hope this all helps
Larry Cuffe
Chemistry Dept
U.C.D
Dublin
From: gaussian.com!fox %! at !% lorentzian.com (Doug Fox)
Subject: Re: CCL:EPR g-factors
To: uunet!blue.weeg.uiowa.edu!jhjensen%lorentzian.com %! at !% uunet.uu.net (Jan
H.
Jensen)
Date: Fri, 25 Sep 1998 15:56:23 -0400 (EDT)
MIME-Version: 1.0
Dr. Jensen,
Gaussian 98 includes anisotropic hyperfine couplings for HF, DFT etc.
Anything which can create a UHF style density.
You can check Chem.Phys. Lett, vol 262. pp 201 (1996) and J.Chem.Phys
vol 105 pp 11060 (1996) for details on implementation and some examples.
Mime-Version: 1.0
Date: Tue, 29 Sep 1998 16:39:26 -0600
To: "Jan H. Jensen" <jhjensen %! at !% blue.weeg.uiowa.edu>
From: schrecke %! at !% t12.lanl.gov (Georg Schreckenbach)
Subject: Re: CCL:EPR g-factors
Dear Jan:
>
> I am interested in calculating anisotropic EPR g-factors using ab
>initio electronic structure theory. I would like to know what programs can
>calculate this property, and what kind of accuracy one can expect.
>
It seems to me that the g-tensor is currently a "hot" field, and
various
groups seem to work on it. This means, on the other hand, that no approach
is well-established yet. Also, I would doubt that any standard program can
do the g-tensor yet.
We have developed a DFT based method that seems to work reasonably well
for first-row radicals with one unpaired electron. Cf. our paper in J.
Phys. Chem. A 1997, 101, 3388; the program is currently not publicly
available. The DFT results were not too bad. One real difficulty turned out
to be the extremely strong influence of the host matrix on experimental
g-tensors.
Another DFT g-tensor paper is J. Chem. Phys. 1997, 107, 2488. For ab
initio, cf. J. Chem. Phys. 1998, 108, 7587 and, in particular, various
papers by Lushington and Grein; they are referenced in the papers cited.
Best regards, Georg
P.S. I am very interested in your summary!
--
==============================================================
Dr. Georg Schreckenbach Tel: (USA)-505-667 7605
Theoretical Chemistry T-12 FAX: (USA)-505-665 3909
M.S. B268, Los Alamos National E-mail: schrecke %! at !% t12.lanl.gov
Laboratory, Los Alamos, New Mexico, 87545, USA
Internet: http://www.t12.lanl.gov/~schrecke/
==============================================================
Your best best for decent results without a facing a huge learning curve is to
use DFT level treatments implemented in conjunction with ADF (Amsterdam Density
Functional). In terms of methodology you have two choices, as outlined in each
of the following:
AUTHOR(s): Schreckenbach, Georg
Ziegler, Tom
TITLE(s): Calculation of the G-Tensor of Electron Paramagnetic
Resonance Spectroscopy Using Gauge-Including Atomic
Orbitals and Density Functional Theory.
In: The journal of physical chemistry. a, molecule
MAY 01 1997 v 101 n 18
Page: 3388
AUTHOR(s): van Lenthe, Erik
Wormer, Paul E.S.
van der Avoird, Ad
TITLE(s): Density functional calculations of molecular g-tensors in
the zero-order regular approximation for relativistic
effects.
In: The journal of chemical physics.
AUG 15 1997 v 107 n 7
Page: 2488
In general one can expect to reproduce experimental trends fairly well for most
systems and get within a 20-30% error in comparing with experimental values for
deviations from the free electron g-value (with the exception of very small
deviations). One usually does better with multi-reference CI calculations
(see: J. Chem. Phys. 106 (1997) 3292 and Int. J. Quantum Chem. 60 (1996) 467)
but these are very labor intensive.
- Gerry
On Wed, 23 Sep 1998, Jan H. Jensen wrote:
>
> Hello,
>
> I am interested in calculating anisotropic EPR g-factors using ab
> initio electronic structure theory. I would like to know what programs can
> calculate this property, and what kind of accuracy one can expect.
>
> Best regards, Jan Jensen
>
> =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
> Jan H. Jensen Assistant Professor
> Department of Chemistry jan-jensen %! at !% uiowa.edu
> University of Iowa Phone:(319) 335-1108
> Iowa City, IA 52242 FAX: (319) 335-1270
> http://www.uiowa.edu/~chemdept/faculty/jensen/
> =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
>
>
>
>
> ---
> Administrivia: This message is automatically appended by the mail exploder:
> CHEMISTRY %! at !% www.ccl.net: Everybody | CHEMISTRY-REQUEST %! at !%
www.ccl.net: Coordinator
> MAILSERV %! at !% www.ccl.net: HELP CHEMISTRY or HELP SEARCH | Gopher:
www.ccl.net 73
> Anon. ftp: www.ccl.net | CHEMISTRY-SEARCH %! at !% www.ccl.net -- archive
search
> Web: http://www.ccl.net/chemistry.html
> ---
>
Gerry Lushington
Ohio Supercomputer Center
1224 Kinnear Rd.
Columbus OH 43212-1163
Ph. 614-292-6036
Fax 614-292-7168