Matter Modeling Asked by Wychh on August 19, 2021

I am interested in calculating the Fock matrix in the molecular orbital basis with PySCF, though I am not clear on the methodology behind this task.

In my attempt, I use the following script (for the example H$_{2}$ molecule):

```
from pyscf import gto, scf
geometry = '''
H 0.000 0.000 0.000
H 0.000 0.000 0.740
'''
mol = gto.Mole()
mol.atom = geometry
mol.basis = '3-21g'
mol.build()
mf = scf.RHF(mol)
mf.scf()
Fao = mf.get_fock()
Fmo = mf.mo_coeff.T @ Fao @ mf.mo_coeff
print('F_mo')
print(Fmo)
```

In this method, I first calculate the molecular mean-field. I then do matrix multiplication with the molecular coefficient transpose matrix (`mf.mo_coeff.T`

), the Fock matrix in the atomic basis (`Fao`

) and the molecular orbital coefficients (`mf.mo_coeff`

).

The resulting off-diagonal matrix elements are essentially zero for the H$_{2}$ molecule and other larger systems taken to 10 decimal places (CH$_{4}$, NH$_{3}$, H$_{2}$O). This has confused me: I have seen other Fock matrices in the molecular orbital basis with off-diagonal elements present.

I am therefore looking for confirmation of my method, and if there is a better way of doing this task?

For **canonical** MOs, it is expected that the off-diagonal elements of the described product with the Fock matrix are zero (within reasonable accuracy) because they are the result of the diagonalization of the Fock matrix.

However, canonical MOs are not the only possible choice. One can "rotate" the MOs to minimize certain metrics in order to obtain, for instance, **localized** MOs (LMO). Such rotations will keep the total energy and the electron density intact, but will disallow the interpretation of the diagonal values of the Fock matrix (represented in the LMO basis) as orbital energies.

Correct answer by TAR86 on August 19, 2021

1 Asked on January 5, 2022 by quantumx

computational chemistry coupled cluster gaussian quantum chemistry

1 Asked on January 5, 2022

ab initio calculations density functional theory functional many body perturbation theory quantum chemistry

1 Asked on January 3, 2022

density functional theory mathematical modeling quantum chemistry spectroscopy td dft

3 Asked on January 1, 2022 by anyon

1 Asked on January 1, 2022 by bnd

1 Asked on December 22, 2021

lammps molecular dynamics molecular interactions molecular modeling software

2 Asked on December 22, 2021

1 Asked on December 20, 2021

electronic structure fcidump integrals one topic per answer software

1 Asked on December 20, 2021 by c-alexander

density functional theory from phys se quantum chemistry spectroscopy td dft

1 Asked on December 20, 2021

1 Asked on December 18, 2021 by kfc

1 Asked on December 16, 2021

density functional theory from phys se mathematical modeling statistical thermodynamics variational principle

1 Asked on December 16, 2021 by luphys

alloys density functional theory from phys se relativistic software

1 Asked on December 16, 2021

1 Asked on December 13, 2021 by roman-korol

atomic physics full configuration interaction high precision molecular modeling periodic system

1 Asked on December 13, 2021

band structure density functional theory from phys se quantum espresso

9 Asked on December 11, 2021

high performance computing matter modeling specific technology one topic per answer open source software

4 Asked on December 9, 2021 by georgios-karaiskakis

computational chemistry density functional theory functional history

2 Asked on December 9, 2021 by koroma

band structure lattice dynamics perovskites phonons symmetry

Get help from others!

Recent Answers

- Jon Church on Why fry rice before boiling?
- Joshua Engel on Why fry rice before boiling?
- Lex on Does Google Analytics track 404 page responses as valid page views?
- haakon.io on Why fry rice before boiling?
- Peter Machado on Why fry rice before boiling?

Recent Questions

- How Do I Get The Ifruit App Off Of Gta 5 / Grand Theft Auto 5
- Iv’e designed a space elevator using a series of lasers. do you know anybody i could submit the designs too that could manufacture the concept and put it to use
- Need help finding a book. Female OP protagonist, magic
- Why is the WWF pending games (“Your turn”) area replaced w/ a column of “Bonus & Reward”gift boxes?
- Does Google Analytics track 404 page responses as valid page views?

© 2023 AnswerBun.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP