Qing Guo,Xiao-Yan Liu,Si-Cheng Liu,Yi-Ning Li,Yi Yin,Peng Zhang
A theoretical analysis of the vibrational modes of ammonium metavanadate† † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ra02053c
Vanadium( v ) is an extremely rare and precious metal, mainly used in aerospace equipment and new energy construction. However, an efficient, simple, and environmentally friendly method for separating V from its compounds is still lacking. In this study, we used first-principles density functional theory to analyse the vibrational phonon density of states of ammonium metavanadate and simulated its infrared absorption and Raman scattering spectra. By analysing the normal modes, we found that the V-related vibration has a strong infrared absorption peak at 711 cm−1, while other significant peaks above 2800 cm−1 are from N–H stretching vibrations. Therefore, we propose that providing high-power terahertz laser radiation at 711 cm−1 may facilitate the separation of V from its compounds through phonon–photon resonance absorption. With the continuous progress of terahertz laser technology, this technique is expected to be developed in the future, and it may offer new technological possibilities.
A theoretical analysis of the vibrational modes of ammonium metavanadate† † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ra02053c

Vanadium(v) is an extremely rare and precious metal, mainly used in aerospace equipment and new energy construction. However, an efficient, simple, and environmentally friendly method for separating V from its compounds is still lacking. In this study, we used first-principles density functional theory to analyse the vibrational phonon density of states of ammonium metavanadate and simulated its infrared absorption and Raman scattering spectra. By analysing the normal modes, we found that the V-related vibration has a strong infrared absorption peak at 711 cm−1, while other significant peaks above 2800 cm−1 are from N–H stretching vibrations. Therefore, we propose that providing high-power terahertz laser radiation at 711 cm−1 may facilitate the separation of V from its compounds through phonon–photon resonance absorption. With the continuous progress of terahertz laser technology, this technique is expected to be developed in the future, and it may offer new technological possibilities.

Vanadium(v) is an extremely rare metal with high melting and boiling points. It is resistant to both hydrochloric acid and sulfuric acid and outperforms most stainless steels in some ways. Because of its extensive utility, it is known as the metal ‘vitamin’. V is mainly used in the production of high-strength low-alloy steels, specialised steels, and aerospace alloys, as well as in other applications. In the field of aerospace, the excellent improvement effect of V in titanium alloys has been discovered, making it useful in body structures and jet engines. The global market for V redox flow batteries is also expected to grow at an annual rate of 59.7% between 2018 and 2025. Considering the production and demand data of V for the energy technology industry specifically, the demand for V is expected to increase by up to 73% by 2050. This dramatic increase will be due to the growing demand for new materials, making mining and recycling of renewable resources strategically important for many countries. Currently, there are various methods of extracting V directly from ores and secondary raw materials, including sodium roasting–water leaching V extraction and calcification roasting–acid leaching V extraction. However, these methods have disadvantages, such as causing pollution due to the use of various additives in the roasting and leaching processes and high costs. It is worth noting that China is one of the countries with the richest V mineral resources. However, an efficient and environmentally friendly method of extracting V has not yet been developed, and there is a lack of research on the recovery of secondary V resources, such as titanium–ferromagnetic slag. Since the discovery of the crystal structure of ammonium metavanadate (NH4VO3) in 1950, many studies using infrared (IR) and Raman spectroscopy have focused on V in NH4VO3. Among these studies, Waal and Twu et al. investigated the pyrolysis process of NH4VO3 under the action of N2 and NH3 + H2O using Raman spectroscopy at the molecular level. However, there has been a lack of theoretical study of the lattice dynamic processes based on vibrational spectroscopy. In this work, we simulated the vibrational spectrum of NH4VO3 and analysed the normal modes. Through our analysis, we were able to assign V-related vibrational peaks and determine the IR-active modes of V in the NH4VO3 spectrum. Based on these findings, we propose a new method to assist V separation using photon–phonon resonance absorption (PPRA).

We performed geometry optimisation and phonon calculations using the CASTEP code, which implements the first-principles density functional theory method. We adopted the generalised gradient approximation in the form of the revised Perdew–Burke–Ernzerhof (RPBE) exchange–correlation functional because the gradient of electron density in NH4VO3 varies widely. The convergence tolerance for energy and self-consistent field (SCF)was set to 1 × 10−9 eV per atom to eliminate virtual frequencies. The energy cut-off was set to 750 eV and the K-point mesh to 3 × 1 × 2 to calculate phonons using a norm-conserving pseudopotential and the linear response method. The property of polarisability, IR and Raman spectra were calculated. The simulated spectra can be compared with experimental data and the vibrational normal modes from the phonon calculation can be used for assignments. NH4VO3 has a pyroxene (Si2O6) structure, and its primitive cell contains 36 atoms with space group Pbcm. The crystal consists of four NH4+ clusters and two VO3− ribbons oriented in one direction. N3− is bonded in tetrahedral geometry to four H1+ atoms in NH4 clusters, while V5+ is bonded to four O2− atoms to form corner-sharing VO4 tetrahedra in each VO3 ribbon. In the analysis of vibrational modes, intermolecular vibrations were classified as translation and rotation modes, while intramolecular vibrations were classified as bending and stretching modes. The study focused on the V-related vibrational modes for further analysis and assignments.

Fig. 1 shows the simulated infrared (IR) spectrum, Raman spectrum, and phonon vibrational density of states (VDOS) of NH4VO3. Considering the normal vibrations at the gamma point, there are 108 normal modes for a 36-atom primitive cell, corresponding to 36 × 3 = 108 dispersion curves that make up the VDOS. Excluding three acoustic branches, there are 105 optical branches that can be detected by IR absorption or Raman scattering through phonon–photon coupling. Due to the unique structure of NH4VO3, its IR- and Raman-active phonons are fully complementary. Table 1 presents a comparison between the calculated normal modes and the experimental IR-active and Raman-active vibrational modes. Park et al. reported the Raman spectrum of NH4VO3 in 1989 and assigned the peaks at 85.5, 149, 170, 195, and 209 cm−1 to translational modes, and the peaks at 269 and 311 cm−1 to rotational modes. Similarly, Du et al. reported the Raman and IR spectra and assigned the Raman-active peaks at 260 and 210 cm−1 to V–O bending vibrations, while Onodera et al. assigned the peaks at 205 and 255 cm−1 (Raman) and the peak at 223 cm−1 (IR) to V–O–V bending vibrations. In this work, we identified corresponding Raman-active vibrational modes at 44, 84, 137, 152, 195, 215, 261, and 312 cm−1, while the IR-active vibrational modes were at 99, 126, 172, 207, and 221 cm−1. In the low-frequency band, the vibrational modes are intermolecular collective vibrations. Our vibrational analysis confirmed that the modes at 44 and 84 cm−1 belong to NH4+ rotations, the modes at 99, 126, and 137 to VO3− rotations, the modes at 152, 172, 195, and 207 to NH4+ translations, and the modes at 215, 221, 261, and 312 cm−1 to skeletal vibrations, where the two ionic groups vibrate together. Fig. 2 presents examples to illustrate these four types of vibrational modes. In Fig. 3, three Raman-active normal modes are evident at 317, 342, and 367 cm−1, which are bending vibrations of the V–O bond. Onodera et al. compared NH4VO3 with KVO3 and assigned a peak at 375 cm−1 to the hydrogen bond vibration between NH4+ and VO3−. Du et al. attributed the mode at 385 cm−1 to the V PBM data was replaced with SVG by xgml2pxml: <glyph-data id="z.dbd" format="PBM" resolution="300" x-size="8" y-size="10" xml:space="preserve"> 00000000 00000000 00000000 00000000 11111111 00000000 11111111 00000000 00000000 00000000 </glyph-data> <svg xmlns="http://www.w3.org/2000/svg" version="1.0" width="13.200000pt" height="16.000000pt" viewBox="0 0 13.200000 16.000000" preserveAspectRatio="xMidYMid meet"><metadata> Created by potrace 1.16, written by Peter Selinger 2001-2019 </metadata><g transform="translate(1.000000,15.000000) scale(0.017500,-0.017500)" fill="currentColor" stroke="none"><path d="M0 440 l0 -40 320 0 320 0 0 40 0 40 -320 0 -320 0 0 -40z M0 280 l0 -40 320 0 320 0 0 40 0 40 -320 0 -320 0 0 -40z"/></g></svg> O bending vibration. Park et al. discovered two Raman peaks at 326 and 344 cm−1, one of medium intensity and the other a shoulder peak. According to Fig. 1, the Raman peaks in this region are very weak. In the IR spectrum, Du et al. assigned various combined vibrations of V–O bonds between 500 and 850 cm−1 and a symmetric vibration of a Raman-active peak at 496 cm−1. Similarly, Onodera and Park assigned the Raman peak at 497 cm−1 to the symmetric vibration. They tentatively assigned the peaks using group theory. Table 1 shows that the simulated normal mode at 489 cm−1 corresponds to a skeletal rotation, as shown in Fig. 3. There are 10 normal modes ranging from 700 to 1000 cm−1, which are all related to V–O stretching. Onodera and Park assigned three Raman peaks at 643, 895, and 925 cm−1, and three at 646, 897, and 936 cm−1, to asymmetric and symmetric V–O stretching vibrations, respectively. Onodera also observed IR peaks at 690, 850, and 935 cm−1. Du assigned two IR peaks at 895 and 935 cm−1 to V–O stretching. Our simulations are in good agreement with the experiments. We found five IR-active modes at 711, 862, 924, 975, and 994 cm−1, and five Raman-active modes at 712, 898, 948, 959, and 994 cm−1. Note that the IR-active mode at 994 cm−1 is not the same as the Raman-active mode at 994 cm−1. These two wavenumbers are not degenerate energies; they appear equal only due to number rounding. As shown in Table 1, these modes all correspond to V–O stretching. Two examples are shown in Fig. 4. The IR-active modes at 711 and 862 cm−1 correspond to the experimental IR peaks at 690, 850 cm−1 (ref. 34), and 895 cm−1 (ref. 30). The dynamic process of V–O asymmetric stretching at 711 cm−1 is presented in the ESI files.† Compared with the other V-related IR-active modes at 862, 924, 975, and 994 cm−1, the intensity ratios to 711 cm−1 are 45.55%, 17.23%, 17.02%, and 18.63%, respectively. It is clear that the IR peak at 711 cm−1 is very strong, which means that the PPRA effect at this peak should be very efficient. If a terahertz laser at this frequency is applied to NH4VO3, the V–O bonds will absorb the radiation energy efficiently, which may potentially break the chemical bonds and facilitate the separation of V from NH4VO3. In the higher-frequency region ranging from 1418 to 3289 cm−1, we found that all vibrations are related to the NH4+ group due to the low mass of the ions. Waal et al. attributed 1414 and 1422 cm−1 to the triply degenerate bending vibrations of NH4+ and pointed out a shoulder peak at 1406 cm−1. Onodera et al. suggested that the Raman peaks at 1420, 1440, and 1650 cm−1 represent NH4+ bending, while the peak at 1412 cm−1 represents NH4+ bending in the IR spectrum. Du et al. reported that the IR peak at 1416 cm−1 was related to the bending vibration of the N–H bond.Table 1 shows that the nine vibration modes from 1418 to 1673 cm−1 all represent NH4+ bending. In the higher-frequency range above 2800 cm−1, the IR and Raman spectra exhibit several distinct characteristic peaks corresponding to N–H stretching vibrations. Waal et al. assigned vibrational peaks at 2839, 2926, and 3019 cm−1 to symmetric stretching, and peaks at 3122 and 3207 cm−1 to triply degenerate asymmetric stretching. Park et al. assigned Raman peaks at 2920 and 3050 cm−1 to NH4+ stretching. Onodera et al. observed IR peaks at 2800, 2980, and 3190 cm−1, while Du et al. assigned IR peaks at 2790, 2950 and 3200 cm−1 to NH4+ stretching. In this work, we identified that the modes at 2874, 2875, 2876, 3285, 3285, and 3289 cm−1 correspond to asymmetric stretching, while those at 2910, 2911, 2912, 2956, 3251, and 3253 cm−1 correspond to symmetric stretching. Although there are some high-intensity IR peaks in this range corresponding to N–H stretching modes, the PPRA effect does not directly facilitate V–O breaking. Two examples of normal modes at 1420 and 2876 cm−1 are shown in Fig. 4.

Based on density functional theory simulations of the VDOS of NH4VO3, we analysed the dynamic processes of the normal modes. The results show that the IR-active modes and Raman-active modes are fully complementary. Each vibrational normal mode is either IR-active or Raman-active. In particular, we confirmed that the normal modes from 711 to 994 cm−1 represent the V–O stretching vibrations. The highest-intensity peak in the IR spectrum is at 711 cm−1, indicating that the PPRA effect of IR radiation at this frequency is very strong. Although there are still some high-intensity peaks in the region above 2800 cm−1, they do not stimulate the PPRA effect of V-related vibrations. V is typically obtained from V-bearing titanomagnetite and ilmenite ore through metallurgical processing, where it is produced as a by-product. The V in the ore is usually in the form of powdered V2O3, which is then dissolved in water or an acidic or alkaline solution to form V-containing ion clusters. The two main chemical methods currently used for the industrial extraction of V are the sodium roasting–water leaching process and the calcium roasting–acid leaching process. Based on our mode analysis, we propose the use of a high-power terahertz laser radiation at 711 cm−1 to assist in breaking the V–O bonds and separating V from NH4VO3. With the continuous progress of terahertz laser technology, this PPRA method could offer new application prospects. By utilising this PPRA physical method, it may be possible to achieve an environmentally friendly and efficient extraction of V from ores.

There are no conflicts to declare.