Unphysical results you got are caused by a quasilinear dependence
in the atomic basis set you are using, as indicated in the output
files. Perturbation theory calculations are typically very sensitive
to such dependences in ill-defined basis sets. To get correct results,
you need to increase the overall precision of all stages involved
into xmcqdpt2 computations. Attached is rar archive with corrected
input files and corresponding outputs from both Firefly v. 8.1.1 and
Firefly v. 8.2.0.
Hope this helps.
On Tue Oct 25 '16 12:07pm, Ole_ wrote
>I'm trying to calculate XMCQDPT2 energies of Pentafluorobenzene along one of its normal modes. The CAS consists of 7 electrons in 6 orbitals, and I'm averaging the lowest 7 states. The resulting energies (MCSCF and XMCQDPT2) in the Franck-Condon structure are very reasonable, and convergence seems good to me.
>When I'm beginning to distort the molecule, the MCSCF potential energy curves of the states are smooth and reasonable. The XMCQDPT2 energies, however, start to jump around and show what I would call unphysical values. At larger distortions, the XMCQDPT2 energies eventually become smooth.
>I included a zip file which has selected input and output files to keep the file size small, as well as plots of the MCSCF and XMCQDPT2 energies to make this more clear (plotted is energy in eV over normal coordinate Q). From the plots, the problem is obvious, I hope. The output of the XMCQDPT2 program seems to be ok (eigenvectors etc), I checked all the points given in the documentation.
>I cannot figure out whats causing this behaviour. I guess it is not a "physical" problem with the MOs or anything, but rather an unfortunate combination of incompatible keywords, because the MCSCF energies are fine. But the XMCQDPT2 step does not require that many stettings... any help would be really appreciated here!
>I'm using the latest Firefly 8.2.0.
>Thank you so much for your help!
>Best regards, Ole
|This message contains the 3118 kb attachment|
[ pfb.rar ]