Disabling Symmetry II
seth carlton olsen
solsen at students.uiuc.edu
Tue Dec 7 20:01:30 GMT 1999
I have included a job below that does not seem recognize the
'nosym' command (it doesn't mention the disabling of symmetry in the
output):
***anion
memory,15,m
file,1,a2ex.int,new
file,2,a2ex.wfn,new
basis=dzp
geomtyp=xyz
geometry={nosym;
21
anion
H,4.418131843,-1.138721853, 0.051255551
H, 4.264339139, 1.404894735,-0.104088897
C,-0.396193228, 0.911431206,-0.160974690
C, 0.541691376,-2.776197552, 0.078609747
N,-1.144786098,-4.549697884, 0.153618288
O,-4.280321237,-1.559377677,-0.085235444
C,-2.576851015, 0.190394947,-0.158690858
N, 1.424893952,-0.026350050,-0.070674000
C, 2.157917393, 1.218944224,-0.132087971
H, 3.489164154,-0.585947040, 0.000941898
C, 0.042676502,-0.080605725,-0.093160601
C,-0.868413811,-1.135962410,-0.046520025
H,-0.510039182,-2.526047438, 0.044173600
C,-1.437419513,-3.507525212, 0.085632956
H,-2.875535278,-3.227568058, 0.042556322
C,-3.225943491,-1.809636335,-0.049792587
C,-2.273666019,-0.851292361,-0.089844735
H, 3.477582231, 0.790900650,-0.082259476
C, 2.329983714,-1.100181572, 0.010763198
H, 1.817605296, 2.387956735,-0.208581960
O,-3.723512970,-4.118806158, 0.081341065
}
gdirect
hf
sym,-1
wf,98,1,0
casscf;occ,50;closed,48;state,2;noextra
casscf;occ,51;closed,47;state,2;noextra
casscf;occ,52;closed,46;state,2;noextra
---
Curiously, the following job, which differs from the former only
in the labeling of the atoms and their
coordinates, DOES seem to recognize the command. The input deck was
generated by taking the first 21 atoms of the cartesian caffeine structure
(an example in the Molpro Manual), and making minor adjustments to the
coordinate input (insertion of commas, etc.) so as to best mimic the
above input deck.
***anion,test deck
memory,15,m
file,1,a2ex.int,new
file,2,a2ex.wfn,new
basis=dzp
geomtyp=xyz
geometry={nosym;
21
anion
C, 0.8423320060,-0.3654865620, 0.0000000000
C,-0.2841017540,-1.1961236000, 0.0000000000
N, 2.0294818880,-1.1042264700, 0.0000000000
N, 0.0774743850,-2.5357317920, 0.0000000000
N,-1.6472646000,-0.6177952290, 0.0000000000
C, 1.4531962870,-2.3678913120, 0.0000000000
C, 0.6373131870, 1.1735112670, 0.0000000000
C,-1.7812691930, 0.7688916330, 0.0000000000
N,-0.6771444680, 1.6306355000, 0.0000000000
O, 1.6106752160, 1.9349693060, 0.0000000000
O,-2.9202890400, 1.2510058880, 0.0000000000
C,-0.9202462430, 3.1094501020, 0.0000000000
C,-2.8623938560,-1.4824503660, 0.0000000000
C, 3.4552156930,-0.6811094280, 0.0000000000
H, 2.0878150460,-3.2451913360, 0.0000000000
H,-1.4989252090, 3.4222116470,-0.8897886280
H,-1.4989252090, 3.4222116470, 0.8897886280
H, 0.0071905670, 3.7148499490, 0.0000000000
H,-3.4903070930,-1.2888938190,-0.8907763360
H,-3.4903070930,-1.2888938190, 0.8907763360
H,-2.6289534570,-2.5638654230, 0.0000000000
}
gdirect
hf
sym,-1
wf,98,1,0
casscf;occ,50;closed,48;state,2;noextra
casscf;occ,51;closed,47;state,2;noextra
casscf;occ,52;closed,46;state,2;noextra
---
I've been trying to figure out at what point the command stops
working as you change the test deck into the first deck, and all that I
can determine is that as soon as you change the atom names and
coordinates, the command is no longer recognized. I'll write again if I
determine anything more. Any help or insight into the problem will be
greatly appreciated.
Seth Olsen
More information about the Molpro-user
mailing list