1
|
Tao YH, Dai X, Moggach SA, Clode PL, Fitzgerald AJ, Hodgetts SI, Harvey AR, Wallace VP. The spectrum of Ih ice using terahertz time-domain spectroscopy. J Chem Phys 2024; 160:214503. [PMID: 38828818 DOI: 10.1063/5.0193458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 05/06/2024] [Indexed: 06/05/2024] Open
Abstract
Here, we report the frequency-dependent spectrum of ice Ih in the range of 0.2-2 THz. We confirm the presence of a feature that blue-shifts from around 1.55-1.65 THz with a decreasing temperature from 260 to 160 K. There is also a change in the trend of the refractive index of ice corresponding to a dispersion, which is also around 1.6 THz. The features are reproduced in data acquired with three commercial terahertz time-domain spectrometers. Computer-simulated spectra assign the feature to lattice translations perpendicular to the 110 and 1̄10 planes of the ice Ih crystal. The feature's existence should be recognized in the terahertz measurements of frozen aqueous solution samples to avoid false interpretations.
Collapse
Affiliation(s)
- Yu Heng Tao
- Department of Physics, The University of Western Australia, Crawley, WA 6009, Australia
| | - Xiangyu Dai
- College of Physics and Optoelectronic Engineering, Shenzhen University, Shenzhen 518060, Guangdong, China
| | - Stephen A Moggach
- School of Molecular Sciences, The University of Western Australia, Crawley, WA 6009, Australia
| | - Peta L Clode
- Centre for Microscopy, Characterisation, and Analysis, The University of Western Australia, Crawley, Western Australia, Australia
| | - Anthony J Fitzgerald
- Department of Physics, The University of Western Australia, Crawley, WA 6009, Australia
| | - Stuart I Hodgetts
- School of Human Sciences, The University of Western Australia, Crawley, WA 6009, Australia
- Perron Institute for Neurological and Translational Science, Nedlands, WA, Australia
| | - Alan R Harvey
- School of Human Sciences, The University of Western Australia, Crawley, WA 6009, Australia
- Perron Institute for Neurological and Translational Science, Nedlands, WA, Australia
| | - Vincent P Wallace
- Department of Physics, The University of Western Australia, Crawley, WA 6009, Australia
| |
Collapse
|
2
|
Mitoli D, Maul J, Erba A. First-Principles Anharmonic Infrared and Raman Vibrational Spectra of Materials: Fermi Resonance in Dry Ice. J Phys Chem Lett 2024; 15:888-894. [PMID: 38241167 DOI: 10.1021/acs.jpclett.3c03372] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2024]
Abstract
We introduce a computational tool for the quantum-mechanical simulation of anharmonic infrared and Raman vibrational spectra of materials. The approach, implemented in the CRYSTAL software, stems from Taylor's expansion of the potential energy surface (PES) on the basis of normal modes up to cubic and quartic terms. The PES can be sampled with four different numerical schemes at the level of density functional theory (DFT), with local, generalized-gradient, and hybrid density functional approximations. Anharmonic states are obtained by solving Shrödinger's nuclear equation with either the vibrational self-consistent field (VSCF) or vibrational configuration interaction (VCI) methods. Nuclear quantum effects (NQEs) are thus fully accounted for. Infrared intensities are computed numerically through a Berry phase approach or analytically through a coupled-perturbed (CP) approach. Raman intensities are computed analytically via the CP approach. A variety of anharmonic features of vibrational spectra of materials can be simulated, including band shifts, combination bands, overtones, resonances (first-order Fermi, second-order Darling-Dennison), and hot bands. We showcase the effectiveness of the approach on the description of a first-order Fermi resonance (FR) in CO2 dry ice: a challenging test-case given that the FR occurs in the Raman spectrum, requires NQEs, and involves two- and three-mode couplings. Fundamental mechanistic differences with respect to the well-known FR in molecular CO2 are addressed. This application represents the first quantum-mechanical, periodic description of FR in dry ice.
Collapse
Affiliation(s)
- Davide Mitoli
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jefferson Maul
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| |
Collapse
|
3
|
Erba A, Desmarais JK, Casassa S, Civalleri B, Donà L, Bush IJ, Searle B, Maschio L, Edith-Daga L, Cossard A, Ribaldone C, Ascrizzi E, Marana NL, Flament JP, Kirtman B. CRYSTAL23: A Program for Computational Solid State Physics and Chemistry. J Chem Theory Comput 2023; 19:6891-6932. [PMID: 36502394 PMCID: PMC10601489 DOI: 10.1021/acs.jctc.2c00958] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Indexed: 12/14/2022]
Abstract
The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.
Collapse
Affiliation(s)
- Alessandro Erba
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jacques K. Desmarais
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Silvia Casassa
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Bartolomeo Civalleri
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Lorenzo Donà
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Ian J. Bush
- STFC
Rutherford Appleton Laboratory, Chilton Didcot, Oxfordshire OX11 0QX, United Kingdom
| | - Barry Searle
- SFTC
Daresbury Laboratory, Daresbury, Cheshire WA4 4AD, United Kingdom
| | - Lorenzo Maschio
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Loredana Edith-Daga
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Alessandro Cossard
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Chiara Ribaldone
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Eleonora Ascrizzi
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Naiara L. Marana
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université
de Lille, CNRS, UMR 8523 — PhLAM — Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department
of Chemistry and Biochemistry, University
of California, Santa
Barbara, California 93106, United States
| |
Collapse
|
4
|
Banks PA, Kleist EM, Ruggiero MT. Investigating the function and design of molecular materials through terahertz vibrational spectroscopy. Nat Rev Chem 2023; 7:480-495. [PMID: 37414981 DOI: 10.1038/s41570-023-00487-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/13/2023] [Indexed: 07/08/2023]
Abstract
Terahertz spectroscopy has proved to be an essential tool for the study of condensed phase materials. Terahertz spectroscopy probes the low-frequency vibrational dynamics of atoms and molecules, usually in the condensed phase. These nuclear dynamics, which typically involve displacements of entire molecules, have been linked to bulk phenomena ranging from phase transformations to semiconducting efficiency. The terahertz region of the electromagnetic spectrum has historically been referred to as the 'terahertz gap', but this is a misnomer, as there exist a multitude of methods for accessing terahertz frequencies, and now there are cost-effective instruments that have made terahertz studies much more user-friendly. This Review highlights some of the most exciting applications of terahertz vibrational spectroscopy so far, and provides an in-depth overview of the methods of this technique and its utility to the study of the chemical sciences.
Collapse
Affiliation(s)
- Peter A Banks
- Department of Chemistry, University of Vermont, Burlington, VT, USA
| | - Elyse M Kleist
- Department of Chemistry, University of Vermont, Burlington, VT, USA
| | | |
Collapse
|
5
|
Desmarais JK, De Frenza A, Erba A. Efficient calculation of derivatives of integrals in a basis of non-separable Gaussians. J Chem Phys 2023; 158:2882252. [PMID: 37094000 DOI: 10.1063/5.0144841] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2023] [Accepted: 04/03/2023] [Indexed: 04/26/2023] Open
Abstract
A computational procedure is developed for the efficient calculation of derivatives of integrals over non-separable Gaussian-type basis functions, used for the evaluation of gradients of the total energy in quantum-mechanical simulations. The approach, based on symbolic computation with computer algebra systems and automated generation of optimized subroutines, takes full advantage of sparsity and is here applied to first energy derivatives with respect to nuclear displacements and lattice parameters of molecules and materials. The implementation in the Crystal code is presented, and the considerably improved computational efficiency over the previous implementation is illustrated. For this purpose, three different tasks involving the use of analytical forces are considered: (i) geometry optimization; (ii) harmonic frequency calculation; and (iii) elastic tensor calculation. Three test case materials are selected as representatives of different classes: (i) a metallic 2D model of the Cu(111) surface; (ii) a wide-gap semiconductor ZnO crystal, with a wurtzite-type structure; and (iii) a porous metal-organic crystal, namely the ZIF-8 zinc-imidazolate framework. Finally, it is argued that the present symbolic approach is particularly amenable to generalizations, and its potential application to other derivatives is sketched.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy
| | - Alessandro De Frenza
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy
| |
Collapse
|