| Both sides previous revision Previous revision Next revision | Previous revision |
| pes_generators [2023/08/16 09:12] – [POTENTIAL ENERGY SURFACES (XSURF)] rauhutmoschneide | pes_generators [2024/11/13 08:27] (current) – rauhutmoschneide |
|---|
| {xsurf,sym=auto} | {xsurf,sym=auto} |
| </code> | </code> |
| The ''XSURF'' program is based on an iterative algorithm, i.e. grid points will be added automatically to the grid representation of the potential until a convergence threshold will be met. This guarantees a well-balanced description of the different terms in the expansion of the potential and simultaneously minimizes the number of //ab initio// calculations for a representation of the potential. For further details see:\\ | As the elongations along different coordinates within the build-up of the PES will lead to different molecular point groups for the individual electronic structure calculations, any symmetry-specific information in that part of the input needs to be avoided. The ''XSURF'' program is based on an iterative algorithm, i.e. grid points will be added automatically to the grid representation of the potential until a convergence threshold will be met. This guarantees a well-balanced description of the different terms in the expansion of the potential and simultaneously minimizes the number of //ab initio// calculations for a representation of the potential. For further details see:\\ |
| |
| G. Rauhut, //Efficient Calculation of Potential Energy Surfaces for the Generation of Vibrational Wave Functions//, [[https://dx.doi.org/10.1063/1.1804174|J. Chem. Phys.]] **121**, 9313 (2004).\\ | G. Rauhut, //Efficient Calculation of Potential Energy Surfaces for the Generation of Vibrational Wave Functions//, [[https://dx.doi.org/10.1063/1.1804174|J. Chem. Phys.]] **121**, 9313 (2004).\\ |
| * **''MINxD''=//n//** The minimum number of coarse grid points can be controlled by the keywords ''%%MIN1D, MIN2D, MIN3D, MIN4D%%''. These 4 keywords determine the minimum number of //ab initio// calculations in one dimension for each 1D, 2D, 3D and 4D surface. The defaults are currently ''%%MIN1D=4, MIN2D=4, MIN3D=4, MIN4D=4%%''. Presently, values larger than 24 are not supported. | * **''MINxD''=//n//** The minimum number of coarse grid points can be controlled by the keywords ''%%MIN1D, MIN2D, MIN3D, MIN4D%%''. These 4 keywords determine the minimum number of //ab initio// calculations in one dimension for each 1D, 2D, 3D and 4D surface. The defaults are currently ''%%MIN1D=4, MIN2D=4, MIN3D=4, MIN4D=4%%''. Presently, values larger than 24 are not supported. |
| * **''MINxD_CRIT1''=//n//** Minimum number of ab initio calculations per dimension, which are needed before a subsurfaces can be declared to be converged according to criterion 1. For dimension 1-3 the default is 6, while it is reduced to 4 for all higher order subsurfaces. | * **''MINxD_CRIT1''=//n//** Minimum number of ab initio calculations per dimension, which are needed before a subsurfaces can be declared to be converged according to criterion 1. For dimension 1-3 the default is 6, while it is reduced to 4 for all higher order subsurfaces. |
| * **''MPG''=//n//** Symmetry of the normal modes is recognized by the program automatically. Only Abelian point groups can be handled at the moment. Symmetry of the modes will be determined even if the ''NOSYM'' keyword is used in the electronic structure calculations. In certain cases numerical noise can be very high and thus prohibits a correct determination of the symmetry labels. Symmetry can be switched off by using ''MPG=1''. | * **''MPG''=//n//** Symmetry of the normal modes is recognized by the program automatically. Symmetry of the modes will be determined even if the ''NOSYM'' keyword is used in the electronic structure calculations. In certain cases numerical noise can be very high and thus prohibits a correct determination of the symmetry labels. Symmetry can be switched off by using ''MPG=1''. |
| * **''NDIM''=//n//** The keyword ''NDIM=n'' terminates the expansion of the PES after the $n$-body term. The default is set to 3. Please note, when you use ''NDIM=4'' as a keyword for the ''XSURF'' program, you need to pass this information to the ''VSCF'' and ''VCI'' programs also. Otherwise these programs will neglect the 4-body terms. | * **''NDIM''=//n//** The keyword ''NDIM=n'' terminates the expansion of the PES after the $n$-body term. The default is set to 3. Please note, when you use ''NDIM=4'' as a keyword for the ''XSURF'' program, you need to pass this information to the ''VSCF'' and ''VCI'' programs also. Otherwise these programs will neglect the 4-body terms. |
| * **''NGRID''=//n//** Based on a coarse grid of //ab initio// points a fine grid will be generated from automated interpolation techniques. The keyword ''NGRID=n'' determines the number of equidistant grid points in one dimension. ''NGRID=n'' has to be an even number. The default is currently set to 16. | * **''NGRID''=//n//** Based on a coarse grid of //ab initio// points a fine grid will be generated from automated interpolation techniques. The keyword ''NGRID=n'' determines the number of equidistant grid points in one dimension. ''NGRID=n'' has to be an even number. The default is currently set to 16. |
| The ''INTENSITY'' directive of the ''SURF'' program provides the option to alter the electronic structure methods for calculating the dipole surfaces. It also allows to define the VARDIP//n//D[X,Y,Z] variables separately. $n$ describes the dimension of the coupling surface and can be chosen to be 1 - 4. | The ''INTENSITY'' directive of the ''SURF'' program provides the option to alter the electronic structure methods for calculating the dipole surfaces. It also allows to define the VARDIP//n//D[X,Y,Z] variables separately. $n$ describes the dimension of the coupling surface and can be chosen to be 1 - 4. |
| |
| Dipole surfaces can be computed for all those methods for which analytical gradients are available in Molpro. For all methods except Hartree-Fock this requires the keyword ''%%CPHF,1%%'' after the keyword for the electronic structure method. In multi-level schemes for which the variables ''VAR1D'', ''VAR2D'' and ''VAR3D'' are set individually, the VARDIP//n//D[X,Y,Z] variables have to be set accordingly. Symmetry is currently only implemented for the 1D, 2D and 3D dipole surfaces. For 4D terms symmetry will automatically switched off at the moment. The determination of dipole surfaces beyond Hartree-Fock quality effectively doubles the computation time for surface calculations. | Dipole surfaces can be computed for all those methods for which analytical gradients are available in Molpro. For all methods except Hartree-Fock this requires the keyword ''%%CPHF,1%%'' after the keyword for the electronic structure method. In multi-level schemes for which the variables ''VAR1D'', ''VAR2D'' and ''VAR3D'' are set individually, the VARDIP//n//D[X,Y,Z] variables have to be set accordingly. The determination of dipole surfaces beyond Hartree-Fock quality effectively doubles the computation time for surface calculations. |
| |
| === Options === | === Options === |
| === Options === | === Options === |
| |
| * **''ANGLE''=//value//** Rotation angle in degree. | * **''ANGLE''=//value//** Rotation angle in degree. This keyword requires the specification of two normal coordinates provided by ''NM1'' and ''NM2''. |
| * **''BLx''=//n//** This option assigns mode $n$ to block $x$. For example, the line ''%%LINCOMB,BL1=1,BL1=2,BL2=5,BL2=6%%'' for a 4-atomic molecule assigns modes 1 and 2 to a 1st block and thus these two modes will be localized afterwards. Modes 3 and 4 will not be affected and thus refer to standard normal coordinates, while modes 5 and 6 constitute a 2nd block | * **''BLx''=//n//** This option assigns mode $n$ to block $x$. For example, the line ''%%LINCOMB,BL1=1,BL1=2,BL2=5,BL2=6%%'' for a 4-atomic molecule assigns modes 1 and 2 to a 1st block and thus these two modes will be localized afterwards. Modes 3 and 4 will not be affected and thus refer to standard normal coordinates, while modes 5 and 6 constitute a 2nd block |
| * **''LOCAL''=//n//** ''LOCAL''=1 localizes the normal coordinates of the CH-stretchings. Note that this destroys symmetry of these modes. Usually localization has strong impact on subsequent ''VSCF'' calculations. ''LOCAL''=2 localizes the normal coordinates of a molecular cluster to the contributing entities. This localization scheme localizes within the individual irreps, which usually leads to a very faint localization. Switching symmetry off by ''MPG''=1 in the ''XSURF'' program leads to a much stronger localization. | * **''LOCAL''=//n//** ''LOCAL''=1 localizes the normal coordinates of the CH-stretchings. Note that this destroys symmetry of these modes. Usually localization has strong impact on subsequent ''VSCF'' calculations. ''LOCAL''=2 localizes the normal coordinates of a molecular cluster to the contributing entities. This localization scheme localizes within the individual irreps, which usually leads to a very faint localization. Switching symmetry off by ''MPG''=1 in the ''XSURF'' program leads to a much stronger localization. |
| === Options === | === Options === |
| |
| * **''TYPE''=//n//** Type of internal Coordinate. **n=** | * **''TYPE''=//variable//** Type of internal Coordinate. **variable=** |
| * 1: Zmatrix coordinates (need to be provided via the Input geometry) | * zmat: Zmatrix coordinates (need to be provided via the Input geometry) |
| * 2: Delocalised internal coordinates | * ozmat: Zmatrix coordinates (internally generated) |
| * 3: Natural internal coordinates | * pzmat: projected Zmatrix coordinates (Zmatrix needs to be provided via the Input geometry) |
| | * opzmat: projected Zmatrix coordinates (Zmatrix is internally generated) |
| | * nzmat: NOMODECO Zmatrix coordinates (Zmatrix needs to be provided via the Input geometry) |
| | * onzmat: NOMODECO Zmatrix coordinates (Zmatrix is internally generated) |
| | * deloc: Delocalised internal coordinates |
| | * natint: Natural internal coordinates |
| | * pnatint: projected Natural internal coordinates |
| | * poly: polyspherical internal coordiantes |
| * **''ECKART''=//n//** If $n=1$, the geometry is rotated according to the eckart conditions which is recommended for the use of other properties than the energy | * **''ECKART''=//n//** If $n=1$, the geometry is rotated according to the eckart conditions which is recommended for the use of other properties than the energy |
| * **''NDIMGMAT''=//n//** Determines the order of the expansion of the Wilson G matrix. Higher than dimension 2 is not recommended | * **''NDIMGMAT''=//n//** Determines the order of the expansion of the Wilson G matrix. Higher than dimension 2 is not recommended |
| * **''NDIMVG''=//n//**. Determines the order of the expansion of the PSeudopotential. Higher than the dimension 2 is not recommended | * **''NDIMVG''=//n//**. Determines the order of the expansion of the Pseudopotential. Higher than the dimension 2 is not recommended |
| * **''RADII''=//x//** $x$ is the threshold for detecting a bon between two atoms | * **''RADII''=//x//** $x$ is the threshold for detecting a bon between two atoms |
| * **''CUTOFF''=//x//** $x$ is the threshold for detecting clusters (not useable currently) | * **''CUTOFF''=//x//** $x$ is the threshold for detecting clusters (not useable currently) |
| * **''ITMAX''=//n//** Defines the maximum number of iteration when elongate the geometries | * **''ITMAX''=//n//** Defines the maximum number of iteration when elongate the geometries |
| * **''PSCRLAB''=//label//** Defines a label for precomputing the surface to determine a first scaling | * **''PSCRLAB''=//label//** Defines a label for precomputing the surface to determine a first scaling |
| | * **''ELONG''=//n//** Defines different ways on how to compute the elongations within the backtransformation from internal to Cartesian coordinates |
| | * For Z-Matrix based coordinates it can be either a simple construction($n=0$) or the HOPE algorithm ($n=1$) |
| | * For NatInt coordinates it can either be a standard B matrix algorithm ($n=0$), virtuel geometries ($n=1$) or subsequent elongation of single coordinates ($n=3$) |
| * **''VAM''=//n//** If $n=1$, rovibrational coupling is considered when the G-matrix is computed--> | * **''VAM''=//n//** If $n=1$, rovibrational coupling is considered when the G-matrix is computed--> |
| ==== Recommendations ==== | ==== Recommendations ==== |
| ''SURF'',//options// | ''SURF'',//options// |
| |
| The old ''SURF'' program for generating potential energy surfaces is still implemented, but will no longer be supported in the future. The new and much improved ''XSURF'' code follows the same philosophy and many options are the same as in the old ''SURF'' code, but ''XSURF'' cannot read any potential files being generated by ''SURF''. If you need to use the old ''SURF'' program, see older versions of the manual for the corresponding keywords. | The old ''SURF'' program for generating potential energy surfaces is still implemented, but is no longer supported. Several keywords are specific to the ''XSURF'' core and cannot be used together with ''SURF''. The new and much improved ''XSURF'' code follows the same philosophy and many options are the same as in the old ''SURF'' code, but ''XSURF'' cannot read any potential files being generated by ''SURF''. If you need to use the old ''SURF'' program, see older versions of the manual for the corresponding keywords. |
| |
| |