printing out the atom number in spinor analysis

48 views
Skip to first unread message

Valérie Vallet

unread,
Oct 22, 2020, 5:14:41 AM10/22/20
to dirac-users
Hi everyone, 

Is there a specific keyword in the ANALYZE block to print out the atom number attached to a given spinor, to help the analysis of EOM calculations, in particular if one wants to know which atom contributes to a given excitation.

Thanks in advance for any advice.

Valérie

Johann Pototschnig

unread,
Oct 23, 2020, 4:10:08 AM10/23/20
to dirac-users
In the Mulliken population analysis it prints how much different atomic orbitals contribute.
If a spinor mostly consists of orbitals from one atom you can assign it to that atom,
which probably can be automatized by a script.

Of course, there will be spinors assigned to multiple atoms and there is the problem of basis set dependence for Mulliken.
In such cases orbital localization can be useful.

Johann

Haase, Pi A. B.

unread,
Oct 23, 2020, 4:41:51 AM10/23/20
to dirac...@googlegroups.com
Dear Valerie,

With "atom number" you mean which particular atom?

I'm not so familiar with the ANALYZE module but would it be an idea in your mol file to append a number to each of the atoms, such as in this example: http://diracprogram.org/doc/release-19/molecule_and_basis/molecule_using_mol.html. Maybe someone else can correct me, but I think this should be the atom labels which are printed in ANALYZE.

Best regards,
Pi

--
You received this message because you are subscribed to the Google Groups "dirac-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dirac-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/dirac-users/868b031c-2d55-4428-8720-723864fed457n%40googlegroups.com.

Trond SAUE

unread,
Oct 23, 2020, 4:44:04 AM10/23/20
to dirac...@googlegroups.com

Hi Valérie,

do you have an example ? Would this be symmetry-independent or symmetry-dependent atoms ?

All the best,

   Trond

On 10/22/20 11:14 AM, Valérie Vallet wrote:

Is there a specific keyword in the ANALYZE block to print out the atom number attached to a given spinor, to help the analysis of EOM calculations, in particular if one wants to know which atom contributes to a given excitation.
--
Trond Saue
Laboratoire de Chimie et Physique Quantiques
UMR 5626 CNRS --- Université Toulouse III-Paul Sabatier
118 route de Narbonne, F-31062 Toulouse, France

Phone : +33/561556361 Fax: +33/561556065
Mail : trond...@irsamc.ups-tlse.fr
Web : http://dirac.ups-tlse.fr/saue
DIRAC : http://www.diracprogram.org/
ESQC : http://www.esqc.org/
Book: Principles and Practices of Molecular Properties: Theory, Modeling, and Simulations

Valérie Vallet

unread,
Oct 24, 2020, 4:34:18 AM10/24/20
to dirac-users
Thanks for the tip Pi. Indeed I managed to get the labels using a mol coordinate file.
However, ideally, one could imagine appending the atom number (order) taken from the XYZ file, this would be helpful as XYZ input files are much easier to generate automatically from eg. MD snapshots. 
Could this option be added to DIRAC?
In my specific case, I want to monitor core orbitals attached to light elements. There might be several carbon/oxygen atoms in the same molecule, and it is crucial to know who is who?

I attach an example of MOL and XYZ file.

Thanks for your help.
hcooh.aug-cc-pvtz.ccsdt.mol
hcooh.aug-cc-pvtz.ccsdt.xyz
Reply all
Reply to author
Forward
0 new messages