QCISD keywords

QCI methods 117 QCISD keywords 9,114 E4T option 186 quadratic Cl methods 114 quadrupole moment 21 quotation li... 

We can compute all of the results except those in the first row by running just three jobs QCISD(T,E4T] calculations on HF and fluorine and a Hartree-Fock calculation on hydrogen (with only one electron, the electron correlation energy is zero). Note that the E4T option to the QCISDfT) keyword requests that the triples computation be included in the component MP4 calculation as well as in the QCISD calculation (they are not needed or computed by default). 

Optimize the structure of acetyl radical using the 6-31G(d) basis set at the HF, MP2, B3LYP and QCISD levels of theory. We chose to perform an Opt Freq calculation at the Flartree-Fock level in order to produce initial force constants for the later optimizations (retrieved from the checkpoint file via OptsReadFC). Compare the predicted spin polarizations (listed as part of the population analysis output) for the carbon and oxygen atoms for the various methods to one another and to the experimental values of 0.7 for the C2 carbon atom and 0.2 for the oxygen atom. Note that for the MP2 and QCISD calculations you will need to include the keyword Density=Current in the job s route section, which specifies that the population analysis be performed using the electron density computed by the current theoretical method (the default is to use the Hartree-Fock density). 

Compute the isotropic hyperfine coupling constant for each of the atoms in HNCN with the HF, MP2, MP4(SDQ) and QCISD methods, using the D95(d,p) basis set Make sure that the population analysis for each job uses the proper electron density by including the Density=Current keyword in the route section. Also, include the 5D keyword in each job s route sectionfas was done in the original study). 

Quadratic Cl (QCI) and coupled cluster (CC) exemplify more complex methods that are not strictly variational in character, but include physical corrections similar to those of higher-order perturbation theory. Keywords for these methods also include a specification of substitutions from the reference FIF configuration, such as QCISD or CCSD, respectively, for QCI or CC methods with all single and double substitutions. More complete descriptions of these methods are beyond the scope of this appendix. 

The E4T option to the UQCISD method keyword tells Gaussian to run fhe component MP4 calcdations at the MP4(SDTQ) level, rather than the default of MP4(SDQ). When you set up the QCISD(T) calculations for methane, include the additional option IOP(2/16=1) in the route section (which says to ignore any symmetry changes during the scan) and also include Gue s=(Always b() for the unrestricted casr. Mix requests that the HOMO and LUMO be mixed so as to destroy alpha-beta and spatial spin symmetries (this option is also useful for producing unrestricted wave functions for singlet systems), and Always says to recompute a new guess wavefunction at each point. 

---