Author Topic: how to get the data of vibration eignvectors ?  (Read 2137 times)

0 Members and 1 Guest are viewing this topic.

Offline berlin

  • Heavy QuantumATK user
  • ***
  • Posts: 31
  • Country: cn
  • Reputation: 0
    • View Profile
how to get the data of vibration eignvectors ?
« on: June 24, 2016, 19:32 »
We can use "Vibration Visualizer" plugin to see the vibration eignvectors. But how to get or save numerical data of vibration eignvectors ?

What i can do now is to read the nc file, and  get some data from "VibrationalMode_gID00*_eigenvectors_real"  and "VibrationalMode_gID00*_eigenvectors_imag". However, I do not know the units of those data. It sames not be  Angstrom. And more, those have been divided by sqrt(mass_of_atom) or not?

« Last Edit: June 25, 2016, 02:29 by berlin »

Offline Jess Wellendorff

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 933
  • Country: dk
  • Reputation: 29
    • View Profile
Re: how to get the data of vibration eignvectors ?
« Reply #1 on: June 27, 2016, 08:42 »
You should use the "eigenvectors" method of the VibrationalMode analyser: http://www.quantumwise.com/documents/manuals/latest/ReferenceManual/index.html/ref.vibrationalmode.html. The returned eigenvector is the solution to an eigernvalue problem involving the Dynamical Matrix, so please also see the documentation for the DynamicalMatrix analysis object: http://www.quantumwise.com/documents/manuals/latest/ReferenceManual/index.html/ref.dynamicalmatrix.html

Offline berlin

  • Heavy QuantumATK user
  • ***
  • Posts: 31
  • Country: cn
  • Reputation: 0
    • View Profile
Re: how to get the data of vibration eignvectors ?
« Reply #2 on: June 27, 2016, 09:53 »
I got the method from the first url:
Quote
mode(mode_index, temperature): Return the vibrational mode with the specified mode_index. The eigenmode,v, is obtained from the eigenvectors,u, as v = u/sqrt(M)*s, where the scaling factor s = sqrt(kB*T)/omega gives the amplitude of a classical oscillator with frequency omega.

But the scaling factor is defined as
s = sqrt(kB*T)/omega
which may only make sense when kB*T >> omega,

For low tempture, any suggestion?
« Last Edit: June 27, 2016, 10:04 by berlin »

Offline Jess Wellendorff

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 933
  • Country: dk
  • Reputation: 29
    • View Profile
Re: how to get the data of vibration eignvectors ?
« Reply #3 on: June 28, 2016, 11:03 »
"s" is just a scaling factor for all the eigenmodes, which is useful for visualization purposes. You may scale all your eigenmodes by some constant, they will still be the vibrational eigenmodes.