Cálculo / visualização da densidade eletrônica a partir de informações orbitais escritas no formato Molden

Question

Estou tentando obter uma densidade de elétrons a partir de informações orbitais escritas no formato Molden pela fórmula, $$ \rho(\mathbf{r}) = \sum_{i}^{occ} |\phi_{i}(\mathbf{r})|^{2}. \tag{1} $$

Existe algum programa de pós-processamento para lidar com este trabalho?

Acredito que posso fazer um script sozinho, pois o arquivo do Molden parece fornecer coeficientes de GTO, mas quero verificar se já existem exemplos funcionando.

Seria desejável que os dados da função densidade final estivessem no formato compatível com softwares de visualização usuais (como Molden, VMD, etc.), mas acho que pode ser visualizado se for representado em uma grade 3D.

Answer 1

Eu recomendo fortemente o cclib, que pode ler Molden e calcular a densidade.

Embora eu não os use tanto, tenho certeza de que você pode fazer isso com:

• Orbkit
• Horton

Sem dúvida, existem outras ferramentas - Molden é um formato bastante padrão para intercâmbio químico quântico.