1 Introduction

Nowadays, there has been considerable interest in the isomorphous replacement of Al atoms by other trivalent cations due to the increased interest in zeolites. These micro-porous species possess different physical and chemical properties from their aluminium analogues and offer possibilities for new applications of molecular sieves [1,2,3,4,5]. Gallium is directly below aluminium in the periodic table; therefore, it forms chemically similar, analogous tetrahedra to AlO4. Isomorphous substitution of Ga for Al in alumosilicate natrolite (Na16Al16Si24O80·16H2O) gives gallosilicate natrolite (Ga–natrolite) [1].

NMR spectroscopy is one of the main methods for studying the microstructure of materials, including zeolites [6]. In present paper, the structure of orthorhombic synthetic gallosilicate natrolite (Na16Ga16Si24O80·16H2O) was studied by means of NMR MAS of 29Si and 71Ga nuclei. The temperature dependence of the spin–lattice relaxation time T1 of 71Ga nuclei has also been studied using solid-state NMR. Obtained results were compared with available references [7,8,9,10,11,12,13,14,15].

2 Experimental Procedure

Gallium form of natrolite was hydrothermally synthesized as described in [12]. Water solution of gallium oxide (Ga2O3, 99.99 + %, Aldrich) and sodium hydroxide (50% aqueous solution, Aldrich) was heated up to 100 °C and stirred over night. After cooling to room temperature, colloidal silica (Ludox AS-40) was slowly added, with stirring. Result mixture with oxide composition of 6.0Na2O·1.0Ga2O3·10.0SiO2·150H2O was charged into Teflon-lined 100 mL autoclave and heated up to 170 °C. After 12 days, the autoclave was cooled down. Products were washed three times with water and dried over night in the oven.

Powder X-ray data (XRD) were collected on PANalytical XPERT-PRO diffractometer using Cu–Kα radiation source (λ = 1.541874 Å). Samples were analyzed over the 2θ range of 5–70° with the step size of 0.017° (Fig. 1). Phases were identified by comparing diffraction peaks with the data reported in Inorganic Crystal Structure Database (ICSD) using X’Pert HighScore program. According to ICSD, obtained XRD pattern matches NAT topology and it is in a good agreement with material TNU-4 as in [12]. This suggests that experimental material is gallium form of natrolite.

Fig. 1
figure 1

Powder X-ray diffraction pattern of gallium natrolite prepared in this study

Chemical analyses of synthesized gallium natrolite and natural natrolite were performed using DSM 982 Gemini FE-SEM scanning electron microscope (SEM) ZEISS. The detector consisted of the inlens and lateral electron (SE) detector and the backscatter electron (BSE) detector. BSE detector was used most frequently for individual SEM images. Respective excitation voltage was 0.265 keV. Details of chemical analyses can be found in Table 1. SEM picture of small crystals produced by hydrothermal synthesis can be found in Fig. 2.

Table 1 Data of chemical analyses of Ga natrolites
Fig. 2
figure 2

Scanning electron micrographs of the gallium natrolite

NMR MAS spectra were obtained on polycrystalline sample of Ga–natrolite on a Bruker Avance-400 NMR spectrometer. 29Si MAS NMR spectra were measured at a spinning rate of 10 kHz using 4 mm rotors at a 29Si frequency of 79.490 MHz with π/2 pulse length of 3 μs. Typically, 10,000 scans were accumulated and the 29Si chemical shifts were referenced to TMS. 71Ga MAS NMR spectra were obtained at a 71Ga frequency of 122.0564 MHz in 4 mm rotors at a spinning rate of 14.0 kHz with an acquisition of 1000 pulse transients, which was repeated with RF pulse length of 2.4 μs. The Dmfit program [16] was used to simulate 71Ga spectra to extract isotropic chemical shifts (\(\delta_{iso}\)), quadrupolar coupling constants (CQ), and asymmetry parameters (η). The spin–lattice relaxation time T1 for 71Ga nuclei was measured by the saturation-recovery method.

3 Results and Discussion

Figure 3 shows 29Si MAS NMR spectra of Al–and Ga–natrolites. From these figure, it follows that 29Si MAS NMR spectra of Al– and Ga–natrolite are similar and contain two resonant lines. In natural Al–natrolite, two resonance lines of 29Si MAS NMR spectrum at– 84.9 ppm and – 92.7 ppm with an intensity ratio of ≈ 2:1 were observed in our investigations (Fig. 3).

Fig. 3
figure 3

29Si MAS NMR spectra of Al–natrolite [13] (a) and Ga–natrolite (b)

The number of resonance NMR lines of the nuclei 29Si and 71Ga determined by the degree of ordered of the natrolite structure [6, 14, 17]. It has been demonstrated at first by Lipmaa et al. [17] that chemical shifts of 29Si MAS NMR spectra of zeolite and other aluminoosilicates are sensitive to the number of AlO4 tetrahedra linked to the SiO4 tetrahedron, also described as the second coordination sphere of silicon. The number of aluminium tetrahedra sharing oxygens with the SiO4 tetrahedron defines five different structural units in the alumosilicate framework, namely, Si(nAl), where index n = 0, 1, 2, 3, 4 specify the number of aluminium tetrahedra connected with SiO4 tetrahedron. 29Si chemical shift ranges by about 5 ppm low-field shifts for each additional connected substituent [17]. This allows to establish structures of various zeolites and get information about Si/Al ordering.

The orthorhombic Al–natrolite (space group Fdd2) has three different central sites (T-sites) in TO4 tetrahedra with multiplicities 8:16:16. According to the X-ray and NMR structure determination study [17, 18], the structure of natural Al–natrolite is well-ordered with Si/Al ratio of 3:2 [6, 17]. There are two types of species with silicon atoms containing Si/Al ratio of 1.5, i.e., those connected to three AlO4 tetrahedra and one SiO4 tetrahedron (Si(AlO4)3(SiO4)) and those connected to two AlO4 tetrahedra and two SiO4 tetrahedra (Si(AlO4)2(SiO4)2) in 2:1 ratio [17].

Two resonance lines at – 81.1 ppm and – 91.1 ppm observed in NMR MAS 29Si spectrum (Fig. 3) of gallosilicate natrolite may indicate that the natrolite has sufficiently well-ordered structure and mainly contains Si(GaO4)3(SiO4) and Si(GaO4)2(SiO4)2) structural units. However, the intensity ratio of these two lines (≈ 1.7) is somewhat unusual and has an Si/Ga ratio different than 1.5. This suggests that in the structure of gallosilicate natrolite, there are other Si(GaO4)n(SiO4)4 − n species with n = 0, 1, 4 structural units [11]. Intensities of NMR MAS 29Si lines of these units may be too weak to be observed experimentally [11]. The difference of chemical shifts between two resonance lines of 29Si atoms in Ga–natrolite is 10.0 ppm.

Gallium atoms have four nearest-neighbour oxygen atoms in tetrahedral coordination in Ga–natrolite. The nucleus of gallium atom has the spin of I = 3/2 and it is a quadrupolar nucleus. The quadrupole moment of the nucleus is a sensitive probe of the existence of a nonuniform electric field (EFG) at the location of the nucleus. The experimental and simulated 71Ga MAS NMR spectra of Ga–natrolite are shown in Fig. 4 which clearly shows a line shape that is determined by the second-order quadrupolar interaction [19, 20].

Fig. 4
figure 4

71Ga MAS NMR spectra of Ga–natrolite. a Experimental spectrum at T = 300 K. b Theoretical spectrum with parameters: the quadrupole coupling constant CQ = e2qQ/h = 5 MHz; the asymmetry parameter η = 0.5; chemical shift δiso = 190.54 ppm. Spinning sidebands are marked by asterisks

The experimental temperature dependence of the spin–lattice relaxation time T1 of 71Ga nuclei in Ga–natrolite is shown in Fig. 5. It is important to note that the relaxation process for 71Ga nuclei in Ga–natrolite was well described by a single exponential. It can be assumed that the main physical mechanism responsible for spin–lattice relaxation of 71Ga nuclei in Ga–natrolite is identical with the mechanism of spin–lattice relaxation of 27Al nuclei in Al–natrolite and the spin–lattice relaxation of 71Ga nuclei in Ga–natrolite is concerned with modulation of electric quadrupolar interactions at the sites of the 71Ga nuclei by thermal motion of water molecules [14]. In this case, the rate of spin–lattice relaxation of 71Ga nuclei may be described by the equation as [14, 20,21,22,23]

$$R_{{1{\text{Q}}}} \equiv T_{{1{\text{Q}}}}^{ - 1} \cong \alpha \cdot \overline{{\left( {\delta C_{\text{Q}} } \right)^{2} }} \left[ {\frac{{\tau_{\text{c}} }}{{1 + \omega_{0}^{2} \tau_{\text{c}}^{2} }} + \frac{{4\tau_{\text{c}} }}{{1 + 4\omega_{0}^{2} \tau_{\text{c}}^{2} }}} \right] ,$$
(1)

where parameter α depends on the asymmetry parameter η of the EFG (from the estimations α ≈ 0.13–0.17 [23]); ω0 is the Larmor frequency of the quadrupole nucleus and

$$\delta C_{\text{Q}} = \frac{{e^{2} (\delta q)Q}}{\hbar }$$
(2)

describes fluctuations of the quadrupole coupling constant [20] at the sites of the 71Ga nuclei, \(\tau_{\text{c}}\) is the correlation time, which describes activated translational and reorientational jumps of electric dipoles in water molecules. In Eq. (1), \(\overline{\left( \cdots \right)}\) represents an ensemble average.

Fig. 5
figure 5

Temperature dependences of the spin–lattice relaxation rate of the 71Ga nuclei in Ga–natrolite. Experimental (circles) and theoretical (continuous line)

From Eq. (1) ,it follows that \(R_{{1{\text{Q}}}} \sim \tau_{\text{c}}\) for \(\omega_{0} \tau_{\text{c}} < < 1\), while \(R_{{1{\text{Q}}}} \sim 1/\tau_{\text{c}}\) for \(\omega_{0} \tau_{\text{c}} > > 1\), with a maximum of \(R_{{1{\text{Q}}}}\) (and minimum of \(T_{{1{\text{Q}}}}\)) in the intermediate region. We assume that the correlation time \(\tau_{\text{c}}\) caused by the molecular motion follows the Arrhenius-type temperature dependence [20]:

$$\tau_{\text{c}} = \tau_{0} \exp \left( {E_{\text{a}} /RT} \right) .$$
(3)

Translational and reorientational jumps of electric dipoles of water molecules were assumed to be dynamically heterogeneous (i.e., described by different activation energies \(E_{\text{a}}\)), and as a result, these were characterized by a normal distribution of the activation energies [24]:

$$p(E_{a} ) = \frac{1}{{\sqrt {2\pi \sigma_{E} } }}\exp \left\{ {\frac{{E_{a} - \overline{E}_{a} }}{{2\sigma_{E}^{2} }}} \right\} .$$
(4)

Therefore, Eq. (1) was averaged over distribution function of Eq. (4). Experimental temperature dependence of the spin–lattice relaxation time and results of their approximation by function of Eq. (1) are presented in Fig. 5. Obtained adjusting parameters are given in Table 2.

Table 2 Calculated parameters that describe the contribution \(R_{1Q} \equiv T_{1Q}^{ - 1}\) to the spin–lattice relaxation

Obtained activation energy 7.3 kcal/mol strongly suggests that translational and reorientational jumps of electric dipoles of water molecules are responsible for the relaxation process of the 71Ga nuclei.

According to the data presented in Table 1, averaged fluctuations of quadrupole coupling constant at the sites of the 71Ga nuclei are

$$\sqrt {\overline{{\left( {\delta C_{\text{Q}} } \right)^{2} }} } \approx 0.3\;{\text{MHz}} .$$
(5)

Full constant of quadrupolar interaction for the 71Ga nuclei in the natrolite \(C_{\text{Q}} = {\text{eqQ}}/h = 5\;{\text{MHz}}\), and therefore, the contribution of the electric dipolar moments of the water molecules to the full EFG at the 71Ga sites is 6%.

4 Conclusions

Based on the analysis of obtained NMR MAS spectra of 29Si and 71Ga nuclei in Ga–natrolites, synthetic gallosilicate natrolite has a well-ordered structure with Si/Ga ratio (3:2) and it contains (Si(GaO4)3(SiO4)) and (Si(GaO4)2(SiO4)2) tetrahedra in the 2:1 ratio. From experimental and simulated 71Ga MAS NMR spectrum of Ga–natrolite, it follows that NMR line shape is determined by the second-order quadrupolar interaction. The spin–lattice relaxation of the 71Ga is governed by electric quadrupole interaction with crystal electric field gradients modulated by translational motion of H2O molecules in natrolite pores.