How to calculate polarizability

From NWChem

You are viewing a single post from the thread title above
Jump to: navigation, search

Click here for full thread
Just Got Here
Threads 1
Posts 2
Dear NWChem friends,

In order to obtain refractive index of a compound, I want to calculate polarizability
tensor at a specified frequency (ex. 589.3nm=0.077315 au ) then I made the input file
as shown below. But, an error message has appeared. Could you tell me what is wrong
with this input file.

echo
start C6H6-POL
title "Polarizability of Benzene"

geometry units angstrom
C -1.209387 -0.698222 0.000000
C -1.209427 0.698262 0.000000
H -2.150069 -1.241320 0.000000
C -0.000014 1.396475 0.000000
H -2.150108 1.241358 0.000000
H -0.000012 2.482676 0.000000
C 1.209387 0.698222 0.000000
H 2.150068 1.241320 0.000000
C 1.209426 -0.698262 0.000000
H 2.150107 -1.241358 0.000000
C 0.000013 -1.396474 0.000000
H 0.000011 -2.482675 0.000000
end
DFT
 XC  b3lyp 
MULT 1
END
basis
C  library 6-31G*
H library 6-31G*
end
TDDFT
RPA
end
PROPERTY
response 1 0.077315
velocity
dipole
end
task tddft energy


Error message
------------------------------------------------------------------------
tddft_init: non-Abelian symmetry not permitted 0
------------------------------------------------------------------------
There is an error related to the specified geometry
------------------------------------------------------------------------


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


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