Support and discussions for Molcas and OpenMolcas users and developers
You are not logged in.
Please note: The forum's URL has changed. The new URL is: https://molcasforum.univie.ac.at. Please update your bookmarks!
You can choose an avatar and change the default style by going to "Profile" → "Personality" or "Display".I want to use Z-matrix input for calculation. symmetry still exists in the xyz coordinate obtained by Z-matrix transformation in log file, but molcas cannot recognize the set point group keywords. How can I use symmetry keywords while using Z-matrix input?
It's not going as follows
&GATEWAY
Title= uo2cl2 in C2V symmetry
symmetry
xy y
XBAS=ANO-RCC-VTZP
zmat
U1
Cl2 1 2.5049656
Cl3 1 2.5049656 2 107.795782
O4 1 1.7488885 3 94.2185517 2 -95.8050332
O5 1 1.7488885 4 165.6565511 3 125.8881242
end of z-matrix
END OF INPUT
not setting symmetry, normal operation and output
&GATEWAY
Title= uo2cl2 in C2V symmetry
XBAS=ANO-RCC-VTZP
zmat
U1
Cl2 1 2.5049656
Cl3 1 2.5049656 2 107.795782
O4 1 1.7488885 3 94.2185517 2 -95.8050332
O5 1 1.7488885 4 165.6565511 3 125.8881242
end of z-matrix
END OF INPUT
&SEWARD
&SCF &End
ITERATIONS
300
PRORbitals
2 5.0
charge
0
&GRID_IT
end of input
Last edited by zz (2022-10-27 13:20:09)
Offline
You have to build your zmat with symmetry-unique atoms only, and make sure it complies with the symmetry requirements (i.e., atoms at symmetry axes or planes are exactly there, otherwise they'll be replicated). You may need to add dummy atoms
symmetry
xy y
zmat
U1
X 1 1.0
Cl1 1 2.5049656 2 53.897891
O2 1 1.7488885 2 97.17138 3 90.0
end of z-matrix
Offline
Small systems can be set manually, but for symmetric systems with more than 50 atoms, is there a program that can help me write the Z-matrix input?
Offline
and when do we need a ghost atom?
Offline
The first atom will be at the origin, the 2nd in the Z axis, and the 3rd in the XZ plane. If that doesn't fit your needs (as in this case, where we don't want any atom in the Z axis), you can use ghost atoms.
Offline
Such an input requires extremely high symmetry, once the atom is a little bit not on the coordinate axis or plane, the symmetry operation will be performed. Is there any keyword to relax the limit value of symmetry operation
Offline
There's SYMTHRESHOLD, but I think it only works with XYZ input.
Offline
Okay, thank you very much!
Offline