1 Introduction

Vanadyl(IV) pyrophosphate [that is, (VIVO)2P2O7] has long been used industrially for the selective oxidation of n-butane to maleic anhydride (MAN) [1]. The catalyst is often referred to as vanadium phosphorus oxide (VPO) due to the presence of multiple vanadium(V) orthophosphate by-phases. Structure types of such by-phases or single stable phases are αI-VOPO4, αII-VOPO4, β-VOPO4, γ-VOPO4, and δ-VOPO4 [2, 3]. These were subject of various studies into their structures and their catalytic properties [2, 3]. It was found that the catalytic reaction mechanism may involve vanadium in an intermediate oxidation state [4]. This, in turn, has encouraged experimental efforts to introduce mixed-valency in vanadyl phosphates, hoping for new catalyst materials [5,6,7].

Electron paramagnetic resonance spectroscopy (EPR) is an established powerful method for investigating systems with unpaired electrons [8]. The dependency of the electron and nuclear Zeeman effects on the external magnetic field, in contrast to the field independence of other terms of the EPR Hamiltonian, such as hyperfine and quadruple interactions, renders acquiring EPR spectra at multiple microwave frequencies helpful and often essential for an unequivocal determination of the plethora of involved magnetic parameters [8, 9]. When a single crystal of the system under investigation is present, the orientation of the magnetic tensors (such as the g and A) and consequently the relative orientations of molecular units and chemical bonds with respect to the crystal axes can be determined by recording EPR spectra at different orientations of the crystal [10].

While vanadium(V) has no unpaired electrons (S = 0) and thus is EPR silent, vanadium(IV) has one unpaired electron (S = ½). Thus, continuous wave (cw) EPR spectroscopy has been invoked previously to study the solid solution (VIVxVV1–2xWVIx)OPO4, which adopts the β‐VOPO4 structure type for x ≤ 0.01 [5]. However, some features of the spectra, and a small splitting of ca. 0.7 mT in particular, could not be explained [5]. Here, this system is revisited by recording cw X-band EPR spectra of a single crystal of such a β-VOPO4 mixed valence species in addition to S-, X- and Q-band EPR spectra of its powder. The resulting consistent set of magnetic parameters from the simulation of all the data sets is reported, and reveals a significant rotation of the g- and A-tensors with respect to each other.

2 Experimental Details

2.1 Synthesis and Crystallization

The previously reported literature procedures were used for the synthesis of β‐VOPO4, WVOPO4, (VIVxVV1–2xWVIx)OPO4, and the making of the single crystals of (V1–xWx)OPO4 [5].

2.2 EPR Measurement

The cw EPR spectra of all samples were acquired on a Bruker Elexsys E580 spectrometer at ambient temperature. Powdered samples were measured at S-, X- and Q-band frequencies (Figs. 1,  2, left), with MgO:Mn2+ as magnetic field calibration standard [11]. For S-, X- and Q-band frequencies (corresponding to ca. 3.5 GHz, ca. 9.7 GHz and ca. 33.7 GHz), resonators ER4118MS5, ER4119HS and EN5107D2 from Bruker were utilized, respectively. Single-crystal EPR measurements (Fig. 4) were performed at X-band, using a programmable E218-1001 goniometer also from Bruker. The crystals were mounted on a homemade mounting rod, which consisted of a quartz tube equipped with a flat-faced Teflon tip. All spectral simulations were performed using EasySpin, a MATLAB® toolbox for simulating and fitting EPR spectra [12].

Fig. 1
figure 1

EPR spectra of the powder sample of (VIV0.01VV0.98WVI0.01)OPO4 a at S-band and b at Q-band microwave frequencies. The simulations are overlaid as red lines. All spectra were acquired at ambient temperature. The simulation parameters are listed in Table 1 (colour figure online)

Fig. 2
figure 2

Shortcoming of the single-frequency EPR. a The experimental X-band EPR spectrum of (VIV0.01VV0.98WVI0.01)OPO4 (black, identical to that of [5]), the simulation with multifrequency parameters (red), and the previously reported simulation [5] (magenta). b The top inset focuses on a small but significant mismatch at the central region, if g and A are considered to be collinear. c The two insets at the bottom zoom in on the most pronounced set of smaller peaks, which are suspected to stem from hyperfine coupling to one or more of the phosphorus atoms. d Q-band comparison of the multifrequency simulation (red), and the simulation with previously reported X-band-extracted parameters in [5] (magenta) (colour figure online)

3 Results and Discussion

The solid solution β-(V0.99W0.01)OPO4, which assumes the β structure type (β-(VIV0.01VV0.98WVI0.01)OPO4) [13], gives rise to an EPR spectrum due to the presence of paramagnetic V(IV) in its structure [5]. The cw S-, X- and Q-band EPR spectra of powdered (V0.99W0.01)OPO4, are shown in Fig. 1. These spectra are dominated by the octuplet hyperfine peaks brought about by the 51 V (I = 7/2) nuclei. A previously reported effort in the simulation of the X-band EPR spectrum of the compound resulted in a presentable spectral agreement albeit failing to reproduce some small but significant features of the spectrum [5]. Using here the previously reported values of the g (g||= 1.935, g = 1.965) and A(51 V) tensor (A||= 504 MHz, A = 159 MHz) drastically fails to reproduce the spectrum at Q-band (Fig. 2, right). A close inspection of the Q-band spectrum, where the g resolution is larger than at S- or X-band, reveals that g is actually orthorhombic and not as previously reported axial. In addition, initial attempts to simultaneous simulate the EPR spectra from all three frequencies indicated that g and A are not collinear.

The orientation of the two tensors with respect to each other was then established via single crystal X-band EPR measurements (Fig. 3).

Fig. 3
figure 3

Single crystal EPR measurement of β-(V0.99W0.01)OPO4 at X-band. The angles refer to the rotations around the crystal axis a (right), and a perpendicular rotation axis in the b–c plane (left) (Fig. 4)

Fig. 4
figure 4

Relative orientation of the molecular units and crystal axes in β-(V0.99W0.01)OPO4 (adopted from [13]). The atoms are colour coded as follows: phosphorus yellow, oxygen in vanadyl purple, other oxygen atoms white, vanadium grey. gz // a and Az tilted away from a by 18° (whether Az is in the ac or a–b plane could not be determined) (colour figure online)

Since the g tensor is practically axial at X-band, it was sufficient to record two sets of EPR spectra, with incremental rotation of the crystal around two perpendicular rotation axes, one parallel to a and one in the b and c plane (Fig. 4) [13]. Figure 3 shows the resulting spectra and their simulations. The agreement is very good and the parameters are consistent with the powder EPR data at S-, X-, and Q-band microwave frequencies. Successful simulation of the single crystal and Q-band powder measurements were only possible with Euler rotation angles α, β, γ of 0°, 18° and 0°, respectively, that is, the g and A tensors are off by an 18° rotation around the y axis (Fig. 4). A comparison with the previously reported parameters (Table 1) shows also that indeed a small g and A orthorhombicity has to be included to fully describe the spectra as well as a non-collinearity of the tensors. Figure 4 also shows the zigzag of the octahedra, and how gz closely aligns with the a axis of the crystal, while Az aligns with the V≡O bonds.

Table 1 Simulation parameters for the powder sample of (VIV0.01VV0.98WVI0.01)OPO4 at S-, X-, and Q-band microwave frequencies (Figs. 1, 2, left) as well as its single crystal EPR spectra at X-band (Fig. 4)

Although the presented simulation model demonstrates a consistent and good agreement over three frequency bands and for powder as well as single crystal, it is not perfect. Of particular interest are the two sets of line splittings at X-band, each about 0.7 mT, centered around 328 mT and 360 mT, respectively (Fig. 2, left—bottom insets). We suspect that these doublets may rise from unaccounted-for hyperfine interactions of the phosphorus atoms, as opposed to the presence of two virtually identical species at a 1:1 ratio. There are three types of P atoms, connected to each V through O bridges, at distances ranging from 3.1 to 3.3 Å [13]. These phosphorus atoms were not included in the simulations to avoid overfitting the data.

4 Conclusion

Here, we have shown that multi-frequency EPR in combination with single-crystal EPR was necessary to obtain a full picture of the g and A(51 V) tensors of β-(V0.99W0.01)OPO4. However, the improvement with respect to the previous reported tensors is fairly subtle and does not change the conclusions drawn in the earlier paper [5].