In this case the “scan” option is the one we want to use. Rigid PES are much easier to perform although they can be less inforative in terms of a real dynamic process. Common mistakes and error messages during geometry scans are associated to lack of memory, maximum number of steps reached or energy divergence due to inconsistent geometries when analyzing for an error is always best to take a look at the optimization progress in order to find out at which step the molecule started to break apart. This same procedure holds for scanning any other kind of geometrical variable. for our hypothetical calculation we could write the second line as 2 1 5 8 S 4 45.0 with which the dihedral angle would be scanned from 0.0 degrees to 180.0 degrees in 45.0 degrees steps) In this way, the second step of our hypothetical calculation would look like the molecule on figure 3, and the convergence criteria would be reached much faster. The first line will make all other dihedral angles around the C2-C3 bond to behave in the same way as D, steps is the number of increments in the variable to be performed from the starting geometry and increment is the step size to be taken in the variable (Note this values are always in integer format so a decimal point is always to be used! e.g. Hence, before defining the variable to be scanned we must select carefully, through the use of wildcards all other variables associated with it. This would result in an unefficient method for computing the desired barrier. Since the value of D is kept constant until convergence then the algorithm (provided a decent level of theory is indicated when it comes to molecules more interesting than ethane) will increase the values of dihedral angles D and D until reaching the convergence criteria at a regular structure. If we define the dihedral angle 2-1-5-8 (fig.1) to be scanned from 180.0 deg to 0.0 by 45.0 degrees increments then the second step of the scan would look distorted like the molecule on figure 2. Lets say we want to calculate the energy barrier associated with the rotation of a dihedral angle on ethane. Therefore the use of wildcards is compulsory and it is illustrated below. The worst case scenario would consist of the molecule’s complete distortion from any physically achievable structure, making the calculation end with an error or by providing meaningless results. This means that functional groups can be destroyed on each step and in the best case scenario the optimization algorithm will put it back together again at, of course, some unnecessary computational cost. A common mistake is to define the variable to scan without taking into consideration all the other variables which depend from the first. Using internal coordinates becomes compulsory and a well-defined Z-Matrix is preferable. For the most tasks, where z-matrices are the superior input, I recommend building it by hand.The use of Internal Redundant coordinates (through the Opt=ModRedundant option) must not be overlooked! This option performes a geometry optimization at each step while maintaining the scanned variable constant, which is referred to as a Relaxed Potential Energy Surface (PES) Scan. Whether or not this z-matrix does actually fits the purpose you are needing it for, is a different question. The following file is produced (it is a valid Gaussian input file): # HF/6-31G* Test The complete command is: newzmat -ixyz -ozmat -round -rebuildzmat If you allow rounding ( -round) it will work. Input z-matrix variables are not compatible with final structure. This will likely fail on your input with an error: Therefore you also need to tell Gaussian to rebuild it ( -rebuildzmat). In Gaussian terms, a block of cartesian coordinates is also a z-matrix. This means the first line denotes the number of atoms, the second line denotes a comment, and both need to be skipped. The format you have posted is actually (simple) xmol. According to the manual of the newzmat utility of Gaussian 16, it accepts unadorned Cartesian coordinates ( G16 online manual).
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |