In the frame of the Born-Oppenheimer approximation, nuclear motions in crystals can be simulated rather accurately using a harmonic model. In turn, the electronic first-order density matrix (DM) can be expressed as the statistically weighted average over all its determinations each resulting from an instantaneous nuclear configuration. This model has been implemented in a computational scheme which adopts an ab initio one-electron (Hartree-Fock or Kohn-Sham) Hamiltonian in the CRYSTAL program. After selecting a supercell of reasonable size and solving the corresponding vibrational problem in the harmonic approximation, a Metropolis algorithm is adopted for generating a sample of nuclear configurations which reflects their probability distribution at a given temperature. For each configuration in the sample the "instantaneous" DM is calculated, and its contribution to the observables of interest is extracted. Translational and point symmetry of the crystal as reflected in its average DM are fully exploited. The influence of zero-point and thermal motion of nuclei on such important first-order observables as x-ray structure factors and Compton profiles can thus be estimated.

Nuclear motion effects on the density matrix of crystals: An ab initio Monte Carlo harmonic approach

PISANI, Cesare;ERBA, ALESSANDRO;FERRABONE, MATTEO;DOVESI, Roberto
2012-01-01

Abstract

In the frame of the Born-Oppenheimer approximation, nuclear motions in crystals can be simulated rather accurately using a harmonic model. In turn, the electronic first-order density matrix (DM) can be expressed as the statistically weighted average over all its determinations each resulting from an instantaneous nuclear configuration. This model has been implemented in a computational scheme which adopts an ab initio one-electron (Hartree-Fock or Kohn-Sham) Hamiltonian in the CRYSTAL program. After selecting a supercell of reasonable size and solving the corresponding vibrational problem in the harmonic approximation, a Metropolis algorithm is adopted for generating a sample of nuclear configurations which reflects their probability distribution at a given temperature. For each configuration in the sample the "instantaneous" DM is calculated, and its contribution to the observables of interest is extracted. Translational and point symmetry of the crystal as reflected in its average DM are fully exploited. The influence of zero-point and thermal motion of nuclei on such important first-order observables as x-ray structure factors and Compton profiles can thus be estimated.
2012
137
4
044114-1
044114-11
http://dx.doi.org/10.1063/1.4737419
Cesare Pisani; Alessandro Erba; Matteo Ferrabone; Roberto Dovesi
File in questo prodotto:
File Dimensione Formato  
JChemPhys_137_044114.pdf

Accesso riservato

Tipo di file: PDF EDITORIALE
Dimensione 766.48 kB
Formato Adobe PDF
766.48 kB Adobe PDF   Visualizza/Apri   Richiedi una copia
pap.pdf

Accesso aperto

Tipo di file: POSTPRINT (VERSIONE FINALE DELL’AUTORE)
Dimensione 475.96 kB
Formato Adobe PDF
475.96 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2318/133938
Citazioni
  • ???jsp.display-item.citation.pmc??? 0
  • Scopus 6
  • ???jsp.display-item.citation.isi??? 6
social impact