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.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.