Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision		 Previous revision
pes_transformations [2025/01/16 11:03] – rauhutpes_transformations [2025/05/22 09:29] (current) – rauhut
Line 23: Line 23:
   * **''NDIMDIP''=//n//** Term after which the $n$-body expansions of the dipole surfaces are truncated. The default is set to 3. Note that ''NDIMDIP'' has to be lower or equal than ''NDIM''.   * **''NDIMDIP''=//n//** Term after which the $n$-body expansions of the dipole surfaces are truncated. The default is set to 3. Note that ''NDIMDIP'' has to be lower or equal than ''NDIM''.
   * **''NDIMPOL''=//n//** Term after which the $n$-body expansions of the polarizability tensor surfaces are truncated. The default is 0. ''NDIMPOL'' has to be lower or equal than ''NDIM'' and must be smaller than 4.   * **''NDIMPOL''=//n//** Term after which the $n$-body expansions of the polarizability tensor surfaces are truncated. The default is 0. ''NDIMPOL'' has to be lower or equal than ''NDIM'' and must be smaller than 4.
 +  * **''NDIMQUAD''=//n//** Term after which the $n$-body expansions of the quadrupole tensor surfaces are truncated. The default is 0. ''NDIMQUAD'' has to be lower or equal than ''NDIM'' and must be smaller than 4.
   * **''NDIMMU''=//n//** An n-mode expansion of the $\mu$-tensor will be generated and transformed to polynomials up to 3rd order (default). The order can be changed by the integer passed to the ''NDIMMU'' keyword.   * **''NDIMMU''=//n//** An n-mode expansion of the $\mu$-tensor will be generated and transformed to polynomials up to 3rd order (default). The order can be changed by the integer passed to the ''NDIMMU'' keyword.
   * **''PSUM''=//n//** Maximum number of basis functions to be used within the fits of the surfaces (sum of exponents).   * **''PSUM''=//n//** Maximum number of basis functions to be used within the fits of the surfaces (sum of exponents).
Line 53: Line 54:
   * **''FITMETHOD''=//n//** (=1 Default) Within the iterative build-up of the individual subsurfaces, intermediate fitting will be used. This can be based on true multidimensional Kronecker product fitting (''FITMETHOD''=1) or on fitting along one-dimensional cuts (''FITMETHOD''=2).   * **''FITMETHOD''=//n//** (=1 Default) Within the iterative build-up of the individual subsurfaces, intermediate fitting will be used. This can be based on true multidimensional Kronecker product fitting (''FITMETHOD''=1) or on fitting along one-dimensional cuts (''FITMETHOD''=2).
   * **''FITxD''=//n//** The maximum order of the polynomials used for fitting within the iterative interpolation scheme can be controlled by the keywords ''%%FIT1D, FIT2D, FIT3D, FIT4D%%''. The default is given by 9. However in certain cases higher values may be necessary, but require an appropriate number of coarse grid points, which can be controlled by ''MIN1D'' etc. (See ''XSURF'' options) -->   * **''FITxD''=//n//** The maximum order of the polynomials used for fitting within the iterative interpolation scheme can be controlled by the keywords ''%%FIT1D, FIT2D, FIT3D, FIT4D%%''. The default is given by 9. However in certain cases higher values may be necessary, but require an appropriate number of coarse grid points, which can be controlled by ''MIN1D'' etc. (See ''XSURF'' options) -->
-  * **''ONLY''=//variable//** Sets one fit function for all coordinates. The possible fit functons are the same as for the option ''COORD''.+  * **''ONLY''=//variable//** Sets one fit function for all coordinates. The possible fit functions are the same as for the option ''COORD''. 
 + 
 +==== Scaling factors ==== 
 + 
 +''PSCAL'',//options// 
 + 
 +The 1D potential functions can be scaled by user-defined factors. This allows to fit potentials to experimentally determined frequencies. 
 + 
 +=== Options === 
 + 
 +  * **''MODE''=//n//** This is a mandatory information, which determines the mode to be scaled. 
 +  * **''SFAC''=//value//** Scaling factor to be used for the respective mode. 
 + 
 +The following example shows the scaling of the three modes of water. 
 + 
 +<code> 
 + {poly 
 +  pscal,mode=1,sfac=0.99924787 
 +  pscal,mode=2,sfac=0.99850975 
 +  pscal,mode=3,sfac=0.99900217} 
 +</code>
  
 <!-- <!--
Line 74: Line 95:
   * **''DVEC''=//n//** (=0 Default) This keyword controls the shift vector within the transformation. ''DVEC=0'' sets this vector to zero.    * **''DVEC''=//n//** (=0 Default) This keyword controls the shift vector within the transformation. ''DVEC=0'' sets this vector to zero. 
   * **''ECKART''=//n//** By default the Eckart transformation matrix needed within the ''PESTRANS'' program will be computed explicitly. ''ECKART=0'' replaces the Eckart transformation matrix by a unit matrix.   * **''ECKART''=//n//** By default the Eckart transformation matrix needed within the ''PESTRANS'' program will be computed explicitly. ''ECKART=0'' replaces the Eckart transformation matrix by a unit matrix.
-  * **''NDIM''=//n//** Order of the $n$-mode potential for the transformed potential. The information of the original potential will be taken from the the ''POLY'' calculation. This option is only available for ''XSURF'' calculations.+  * **''NDIM''=//n//** Order of the $n$-mode potential for the transformed potential. The information of the original potential will be taken from the the ''POLY'' calculation. Note, the ''NDIMDIP'' keyword requests the intensity directive. 
   * **''THRSMAT''=//value//** Threshold controlling the analysis of the **S**-matrix indicating the accuracy of the transformation (default: ''THRSMAT=0.3'').   * **''THRSMAT''=//value//** Threshold controlling the analysis of the **S**-matrix indicating the accuracy of the transformation (default: ''THRSMAT=0.3'').
   * **''THRQ''=//value//** Elements in the displacement vectors below this threshold (default ''%%THRQ=10^{-8}%%'') will be neglected within the transformation.   * **''THRQ''=//value//** Elements in the displacement vectors below this threshold (default ''%%THRQ=10^{-8}%%'') will be neglected within the transformation.
Line 111: Line 132:
  
  
-The following example shows a calculation for water including vibrational-rotational coupling surfaces, the corresponding VSCF and VCI calculations a subsequent transformation of the potential to doubly deuterated water.+The following example shows a standard calculation for water, the VSCF and VCI calculations after the call of the ''PESTRANS'' program refer a to doubly deuterated water.
  
 <code> <code>