An extension of the CRYSTAL program is presented allowing for calculations of anharmonic infrared (IR) intensities and Raman activities for periodic systems. This work is a follow-up of two papers devoted to the computation of anharmonic vibrational states of solids from DFT (density functional theory) calculations (Erba et al. J. Chem. Theory Comput. 2019, 15, 3755-3765 and Erba et al. J. Chem. Theory Comput. 2019, 15, 3766-3777). The approach presented here relies on the evaluation of integrals of the dipole moment and polarizability operators over anharmonic wave functions obtained from either VSCF (vibrational self-consistent field) or VCI (vibrational configuration interaction) calculations. With this extension, the program now allows for a more complete characterization of the vibrational spectroscopic features of solids within the density functional theory. In particular, it is able (i) to provide reliable positions and intensities for the most intense spectral features and (ii) to check whether a first overtone or a combination band has a nonvanishing IR intensity or Raman activity. Therefore, it becomes possible to assign the transition(s) corresponding to satellite peak(s) around a fundamental transition or the overtones or combination bands that may be as intense as their corresponding fundamental transitions through the strongest mode-mode couplings, as in so-called Fermi resonances. The present method is assessed on two molecular systems, H2O and H2CO, as well as on two solid state cases, boron hydrides BH4 and their deuterated species BD4 in a crystalline environment of alkali metals (M = Na, K). The solid state cases are particularly insightful as, in the B-H (or B-D) stretching region here considered, they exhibit many spectral features entirely due to anharmonic effects: two out of three in the IR spectrum and four out of six in the Raman spectrum. All IR and Raman active overtones and combination bands experimentally observed are correctly predicted with our approach. The effect of the adopted quantum-chemical model (DFT exchange-correlation functional/basis set) for the electronic structure calculations on the computed spectra is discussed and found to be significant, which suggests some special care is needed for the analysis of subtle spectral features.

Calculation of Anharmonic IR and Raman Intensities for Periodic Systems from DFT Calculations: Implementation and Validation

Carbonniere P.;Erba A.;Dovesi R.;Rerat M.
2020-01-01

Abstract

An extension of the CRYSTAL program is presented allowing for calculations of anharmonic infrared (IR) intensities and Raman activities for periodic systems. This work is a follow-up of two papers devoted to the computation of anharmonic vibrational states of solids from DFT (density functional theory) calculations (Erba et al. J. Chem. Theory Comput. 2019, 15, 3755-3765 and Erba et al. J. Chem. Theory Comput. 2019, 15, 3766-3777). The approach presented here relies on the evaluation of integrals of the dipole moment and polarizability operators over anharmonic wave functions obtained from either VSCF (vibrational self-consistent field) or VCI (vibrational configuration interaction) calculations. With this extension, the program now allows for a more complete characterization of the vibrational spectroscopic features of solids within the density functional theory. In particular, it is able (i) to provide reliable positions and intensities for the most intense spectral features and (ii) to check whether a first overtone or a combination band has a nonvanishing IR intensity or Raman activity. Therefore, it becomes possible to assign the transition(s) corresponding to satellite peak(s) around a fundamental transition or the overtones or combination bands that may be as intense as their corresponding fundamental transitions through the strongest mode-mode couplings, as in so-called Fermi resonances. The present method is assessed on two molecular systems, H2O and H2CO, as well as on two solid state cases, boron hydrides BH4 and their deuterated species BD4 in a crystalline environment of alkali metals (M = Na, K). The solid state cases are particularly insightful as, in the B-H (or B-D) stretching region here considered, they exhibit many spectral features entirely due to anharmonic effects: two out of three in the IR spectrum and four out of six in the Raman spectrum. All IR and Raman active overtones and combination bands experimentally observed are correctly predicted with our approach. The effect of the adopted quantum-chemical model (DFT exchange-correlation functional/basis set) for the electronic structure calculations on the computed spectra is discussed and found to be significant, which suggests some special care is needed for the analysis of subtle spectral features.
2020
16
5
3343
3351
Carbonniere P.; Erba A.; Richter F.; Dovesi R.; Rerat M.
File in questo prodotto:
File Dimensione Formato  
ANHARM_INTENS_JCTC.pdf

Accesso riservato

Tipo di file: PDF EDITORIALE
Dimensione 2.93 MB
Formato Adobe PDF
2.93 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
MANUSCRIPT_INT_ANHARM_CRYSTAL_PC_AE_FR_RD_MR_REVISED_V3.docx

Accesso riservato

Tipo di file: POSTPRINT (VERSIONE FINALE DELL’AUTORE)
Dimensione 743.25 kB
Formato Microsoft Word XML
743.25 kB Microsoft Word XML   Visualizza/Apri   Richiedi una copia
MANUSCRIPT_INT_ANHARM_CRYSTAL_PC_AE_FR_RD_MR_REVISED_V3.pdf

Open Access dal 29/04/2021

Tipo di file: POSTPRINT (VERSIONE FINALE DELL’AUTORE)
Dimensione 492.84 kB
Formato Adobe PDF
492.84 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/1751676
Citazioni
  • ???jsp.display-item.citation.pmc??? 0
  • Scopus 8
  • ???jsp.display-item.citation.isi??? 8
social impact