Differences

This shows you the differences between two versions of the page.

--- geometry_optimization_optg [2023/08/01 09:25] – [Options for the PQS optimizer] werner
+++ geometry_optimization_optg [2024/01/08 13:24] (current) – external edit 127.0.0.1
@@ Line 77: / Line 77: @@
-There are many further options, but only the following may sometimes be relevant:
+There are many further options, but only the following may be sometimes relevant:
 ''OPTC:''
@@ Line 97: / Line 97: @@
   * ''ITYPE=0'' - standard molecular optimization (default)
   * ''ITYPE=1'' - optimization of molecule reacting/adsorbed on model surface
-  * ''ITYPE=2'' - optimization of molecular clusters (requires ''MOLECULES'' option)
+  * ''ITYPE=2'' - optimization of molecular clusters (see below)
 ''ICONS:''
   * ''ICONS=0'' - full optimization (no constraints, default)
   * ''ICONS=1'' - constrained optimization using Lagrange multipliers (default if constraints are given see [[geometry optimization (OPTG)#Constraint optimisations using the PQS optimizer|constraint optimisations using the PQS optimizer]])
-  * ''ICONS=2'' - constrained optimization using Penalty functions
+  * ''ICONS=2'' - constrained optimization using penalty functions
-''MOLECULES:''
+''Cluster optimizations:''
-  * ''MOLECULES=[n1,n2,n3,...]'' - For cluster optimisations. The molecules must come one after the other in the geometry input. ''n1'', ''n2'',... are the number of atoms in each molecule. The sum of these must be equal to the total number of atoms. This option sets automatically ''ITYPE=2''
-Unfortunately, the cluster optimization in pqsopt is buggy and currently not working.
+Cluster optimizations use inverse distance coordinates for intermolecular degrees of freedom.
+  * ''MOLECULES=[//n1,n2,n3,...//]'' - For cluster optimisations. The molecules must come one after the other in the geometry input. //''n1'', ''n2'',...// are the number of atoms in each molecule. The sum of these must be equal to the total number of atoms. This option sets automatically ''ITYPE=2''
+  * ''CUTOFF=//value//'' - Distance cutoff for intermolecular coordinates in Angstrom. The program will generate inverse distance coordinates for between all atoms of different molecules, unless the distance is larger than this value. Default: 5 Angstrom.
+  * ''SCALE=//value//'' - Scaling factor for inverse distance coordinates. Typical values are between 1 and 10. Default is 1.0
+Unfortunately, according to our current limited experience, the cluster optimization in pqsopt is neither well working in the original PQS program nor in Molpro and should therefore be used with care.
+''Optimizations with z-matrix input:''
+By default, optimizations with z-matrix geometry input are done with the old RF optimizer, and only the active variables in the z-matrix are optimised (cf. [[geometry optimization (OPTG)#Defining active geometry parameters (ACTIVE)|Defining active geometry parameters]]). It is possible, however, to force pqs optimization by giving the optg option ''PQS'' or ''METHOD=PQS''. In this case, all cartesian coordinates will be optimised, and in each iteration the new structure is back-transformed to the z-matrix. If the z-matrix contains constants, these may change and will be overwritten. If the same variable is used in several places of the z-matrix and the corresponding new parameters are not the same, the variable is only kept in the first place, and the further ones are overwritten by the new fixed parameters. The updated z-matrix is printed in the log-file in each iteration.
 ==== Constraint optimisations using the PQS optimizer ====