Support and discussions for Molcas and OpenMolcas users and developers
You are not logged in.
Please note: The forum's URL has changed. The new URL is: https://molcasforum.univie.ac.at. Please update your bookmarks!
You can choose an avatar and change the default style by going to "Profile" → "Personality" or "Display".Dear all,
I am calculating UV and X-ray absorption spectrum with OpenMolcas. I would like to know how to visualize the excitation. For instance, in a valence excitation, one electron is excited from HOMO to LUMO level. In LR-TDDFT calculations, Q-Chem or NWchem will output the related Kohn-Sham orbitals of one excitation. How to find out similar information from MOLCAS output?
Is the Natural orbitals the same for different roots of RASSCF(CASPT2) calculations?
Thanks!
Offline
Natural orbitals are different for each root, so you cannot just compare the occupations; but if you compare the orbitals and find out they are almost the same, the difference in occupations can give you a hint. The "right" way is to examine the CSF coefficients, which refer to the average orbitals, and these are the same for each root. Of course, analyzing the CSF coefficients can be hard, so you can use the WFA tool to (https://plassergroup.wordpress.com/2018 … in-molcas/) to help you get a more visual picture.
Offline
Thanks a lot!
Offline
Dear Ignacio,
Do you have more tutorial or manual on libwfa or any similar software? I haven't figured out how to visualize the excitation or the average orbitals you mentioned. It seems that libwfa write some visualization information to the hdf5 it reads. But pegamoid.py crashes when I tried to visualize it.
Thanks a lot!
Offline
Pegamoid hasn't been tried with WFA files yet, so there could be some incompatibilities. Maybe you can ask in https://github.com/libwfa/libwfa, if Felix doesn't read this.
If you stay within RASSCF (no RASSI or CASPT2), you can get the transition orbitals with the TDM option, and that should work with Pegamoid.
Offline