Skip to content

Commit dbb1a02

Browse files
committed
Corrections to add symmetry check for ORCA
1 parent 94dc459 commit dbb1a02

File tree

1 file changed

+9
-3
lines changed

1 file changed

+9
-3
lines changed

GSASII/exports/G2export_PDB.py

Lines changed: 9 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -276,15 +276,21 @@ def Exporter(self,event=None):
276276
Cell = General['Cell'][1:7]
277277
A,B = G2lat.cell2AB(Cell)
278278
self.Write('# GSAS-II generated ORCA input file\n# Basic Mode')
279+
self.Write('#%Sym SymThresh 0.01 End')
279280
self.Write('# after xyz is charge & multiplicity; change as needed')
280281
self.Write('* xyz 0 1')
281282
fmt = '{:4s}'+3*'{:12.4f}'
282283
for atom in Atoms:
283284
xyz = np.inner(A,np.array(atom[cx:cx+3]))
284285
self.Write(fmt.format(atom[ct],*xyz))
285-
self.Write('*\n%%Compound\nGeometry %s;\nVariable Optimize = 0;\nVariable CC;'%Oname)
286-
self.Write('if (Optimize=0) then\n New_Step\n ! RHF SP def2-SVP\n Step_End')
287-
self.Write('else\n New_Step\n ! BLYP SVP Opt\n Step_End\nEndIf')
286+
self.Write('*\n%%Compound\nGeometry %s;\nVariable Optimize = 0;Variable CC;'%Oname)
287+
self.Write('if (Optimize=0) then\n New_Step\n ! RHF SP def2-SVP')
288+
self.Write(' #%Sym SymThresh 0.01 End\n # uncomment to do symmetry check')
289+
self.Write(' Step_End')
290+
self.Write('else\n New_Step\n ! BLYP SVP Opt')
291+
self.Write(' #%Sym SymThresh 0.01 End\n # uncomment to do symmetry check')
292+
self.Write(' Step_End')
293+
self.Write('EndIf')
288294
self.Write('%s.Read();'%Oname)
289295
self.Write('%s.WriteXYZFile(filename="%s_ORCA.xyz");'%(Oname,Oname))
290296
self.Write('End')

0 commit comments

Comments
 (0)