VSCF - cubic and quartic anharmonicities

From NWChem

Viewed 99 times, With a total of 2 Posts
Jump to: navigation, search

Just Got Here
Threads 1
Posts 2
Dear all,

I have carried out a vscf calculation and I could find the VSCF anharmonic frequencies in the output. However, I am also looking for the cubic and quatric anharmonicities and I could not find them in the output or the ecce.out file.
Does anyone know if and how those can be extracted?

There are however several output files like the *.coup file but it is not clear to me what the numbers in this file actually represent.

Best wishes
Nicole

Gets Around
Threads 1
Posts 144
VSCF doesn't have cubic and quartic anharmonicities because you are not expanding the nuclear potential in a power series. For the implementation in NWChem, you have the option to calculate the effective potential for the vibrational wave functions as a sum of contributions from individual modes, pairs of modes, or triplets of modes, with the code calculating a 1, 2, or 3-dimensional grid of the energy for displacements of the geometry along each relevant mode. The total vibrational wave function (approximated as a product of single mode wave functions) is then calculated using this effective potential. VSCF is a bit like Hartree-Fock for the nuclear vibrations.

The *.diag file contains the 1-dimensional grid for each mode. The *.coup file it appears was intended to contain the 2-dimensional grid for each pair of modes, but it seems it overwrites itself all the time so only the last result for the last pair of modes ends up in the file at the end of the calculation. For the 3-dimensional grid for triplets of modes, it looks like the grid would get written to the normal output file, but I'm not sure because I've never run one of those calculations since they are very expensive for anything but very small systems.

Just Got Here
Threads 1
Posts 2
Thanks a lot for the explanation, Sean!
That was very helpful to me.

Best wishes,
Nicole


Forum >> NWChem's corner >> NWChem functionality



Who's here now Members 0 Guests 0 Bots/Crawler 1


AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC