Multicomponents of spin-spin relaxation, anisotropy of the echo decay, and nanoporous sample structure

We have experimentally and theoretically investigated multicomponent 1H nuclear magnetic resonance (NMR) echo decays in a-Si:H films containing anisotropic nanopores, in which randomly moving hydrogen molecules are entrapped. The experimental results are interpreted within the framework of the previously developed theory, in which a nanoporous material is represented as a set of nanopores containing liquid or gas, and the relaxation rate is determined by the dipole–dipole spin interaction, considering the restricted motion of molecules inside the pores. Previously, such characteristics of a nanostructure as the average volume of pores and their orientation distribution were determined from the angular dependences of the spin–spin and spin–lattice relaxation times. We propose a new approach to the analysis of the NMR signal, the main advantage of which is the possibility of obtaining nanostructure parameters from a single decay of the echo signal. In this case, there is no need to analyze the anisotropy of the relaxation time T2, the determination of which is a rather complicated problem in multicomponent decays. Despite multicomponent signals, the fitting parameter associated with the size and shape of nanopores is determined quite accurately. This made it possible to determine the size and shape of nanopores in a-Si:H films, herewith our estimates are in good agreement with the results obtained by other methods. The fitting of the decays also provides information about the nanostructure of the sample, such as the standard deviations of the angular distribution of pores and the polar and azimuthal angles of the average direction of the pore axes relative to the sample axis, with reasonable accuracy. The approach makes it possible to quantitatively determine the parameters of the non-spherical nanoporous structure from NMR data in a non-destructive manner.


Introduction
Nuclear magnetic resonance (NMR) and nuclear magnetic resonance imaging (MRI) are the most sensitive and powerful tools for identifying the structure of both pure compounds and mixtures as either solids or liquids [1][2][3][4]. A remarkable feature of these methods is their ability to be used for research on different levels of organization of matter, from atoms [1,2] and molecules [5] to solids [3,4], cells [6,7], membranes [7], and tissues of living organisms [7,8].
The spin-spin interactions play an important role in NMR and MRI due to the fact that coupling constants characterizing the interactions depend on the structural parameters of the tested materials [1][2][3][4][5]. In solids, dipole-dipole interactions (DDI) of nuclear spins are often the main interactions, which are responsible for the shape of NMR spectra and spin-spin and spinlattice relaxation processes [1][2][3][4]. This is especially true for the spins of hydrogen atoms [9] owing to the weak electron-nuclear coupling. Measurements of the relaxation times allow one to study the dynamic and structural properties of materials, including their microstructure and composition [3,5,8].
In recent years, NMR and MRI methods have been used to study the nanostructures of nanoporous materials, in which mobile gas or liquid molecules are trapped in pores less than a few hundred nanometers in size [10][11][12][13]. Such materials can be found commonly in nature, including biological systems and living organisms, and possess a number of unique properties.
One of the main advantages of magnetic resonance methods is their non-invasiveness and nondestructiveness.
In our previous works, we studied the relaxation processes in molecules that are entrapped in nano-sized cavities (nanopores) of the materials [14][15][16][17][18][19][20][21][22]. We have demonstrated that the parameters of nanostructures can be extracted from the angular dependence of the relaxation times T1, T1ρ, and T2. To explain in detail the anisotropy of the relaxation times for liquid and gas entrapped in nanopores that have a characteristic size up to hundreds of nanometers, the model considered the restricted random motion of molecules at averaging the spin Hamiltonian [14][15][16][17][18][19][20][21][22].
It was shown [16][17][18][19][20][21][22] that the averaged Hamiltonian which describes the dipole-dipole interactions in the nanopores can be characterized by an effective coupling constant. Averaging the spatial part of the dipolar Hamiltonian, we developed the models which describe very nontrivial spin dynamics and allow us to obtain analytical expressions for the relaxation times T1, T1ρ, and T2. These expressions establish the relationship between the sample nanostructure and the anisotropy of the relaxation times. However, determining the relaxation times, especially the transverse relaxation time from experiments where the measured signal is multicomponent, is not a well-defined procedure [23][24][25][26]. Usually, a relaxation time is determined using one-or two-exponential functions for fitting experimental data [25,26]. It was shown [14][15][16][17][18][19][20][21][22] that the nanostructure of nanoporous materials can be analyzed using the angular dependence of the relaxation time obtained by a simple fitting of experimental data for a sample, that is put at different angles relative to an applied magnetic field. The fitting by one-or two-exponential functions is more problematic for nanoporous samples where the measured signal is formed by a superposition of signals from a set of nanopores with different orientations, volumes, etc. [14,[24][25][26].
The purpose of the present study was to establish a correlation between the structural features of nanoporous materials (orientation of nanopores, their volume and shape) and the echo decays, bypassing the procedure of determination of the transverse relaxation time. The advantage of the proposed method is the approximation of experimental data by a multicomponent exponential function. The developed approach opens up a way of studying nanostructures by analyzing single echo decay.
This paper is organized as follows: in Section 2 for convenience of the readers, we briefly show that the angular dependence (anisotropy) of the transverse relaxation time results from the residual dipole-dipole interaction in liquid or gas entrapped in nanopores. In Section 3, we describe the experimental details. The method to determine the internal nanostructural parameters in hydrogenated amorphous silicon (a-Si-H) films by using the theoretical model and the experimental data is considered in Section 4. In the last section, we summarize our results.

THEORY
In this section, we first describe the spin-spin relaxation in a single nanopore containing liquid or gas molecules and then consider the effect of nanopore orientation disorder in a sample on the echo decay.

Echo decay in a single nanopore
For gas or liquid entrapped in a nanopore with sizes less than 1 m the Hamiltonian of a spin system can be averaged over the random motion of molecules [14][15][16][17][18][19][20][21][22]. The averaged Hamiltonian describing the dipole-dipole (DD) interaction in the magnetic field 0 H along the zaxis ( 0 Hz ) is given by [27,28]: where  is the nuclear gyromagnetic ratio,  is the reduced Planck constant, V is the nanopore volume and F is the form-factor [28],  denotes the angle between the main C Z -axis of the nanopore and the external magnetic field 0 H ( Fig. 1), zj I is the projections of the spin operator j I onto the Z − axis of the j -th nuclear spin. Eq (1) was obtained for ellipsoidal or cylindrical nanopores. Consideration of nanopores as ellipsoids or cylinders is correct in many cases, including nanopores in a-Si:H films [21,22,[29][30][31][32][33].
. To detect the anisotropy of the relaxation processes, the sample is put at different angles relative to the direction of the applied magnetic field. In this case, the measure of change in orientation of the sample is the angle θS between the normal to the sample surface and the direction of the magnetic field. Therefore, to describe all possible orientations of a nanopore, it is necessary to rewrite Eq. (1), in the form accounting the sample rotation: where ( ) , , cos cos sin sin cos S  is the angle between the external magnetic field 0 H and the normal to the sample surface Naxis ( Fig.1).
 is the spin density of liquid or gas.
According to the phenomenological Bloch equations a transverse relaxation from a single nanopore can be described by an exponential decay: where M0 is the initial signal intensity at t = 0.

Echo decay in a sample with orientational distribution of nanopores
Eqs. primordially, we assumed that the distribution is done by the binary Gaussian function of the polar and azimuthal angles, which yielded good agreement between the experimental data and the theoretical results for various samples [14][15][16][17][18][19][20][21][22], and that the angular distributions of all types of pores in a sample are the same. Replacing pore volumes and form factors by their effective values, we obtain the normalized signal averaged over the whole sample: where ( ) ,   is the density function that is assumed to be described by a binary Gaussian 7 function:

Sample
The studied a-Si:H films were grown at the Laboratoire de Physique des Interfaces. The mm onto the c-Si substrate [34].

NMR measurements
The 1 H NMR measurements were carried out at ambient temperature (293±0.5K) using a Tecmag pulsed solid-state NMR spectrometer, a home-built NMR probe and an Oxford superconducting magnet with the magnetic field B0 = 8.0 T corresponding to the 1 H resonance The echo decay of protons in hydrogen molecules in the a-Si:H films was measured at different angles s  between the normal to the film surface and the applied magnetic field (Fig.1). The angles were established within the accuracy of 0 1  . The spin-spin relaxation times T2 were measured using the Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence [38]. The duration of the π/2 radiofrequency pulse in the sequences was 1.6 μs and the amplitude of the radio frequency field was of H1 = 39 G, which corresponds to the frequency 11 166.6 kHz H  == . All measurements were carried out in a static mode.
All the decays are similar; some of them are very closed to each other due to nonmonotonic decay dependence on the angle θS. As an example, the results of our measurements of the echo decay of the a-Si:H films for the angles θS =40 0 , θS =65 0 , and θS =90 0 are shown in Fig. 2. The results for these angles clear demonstrate the multicomponent character of spin-spin relaxation and its angular dependence.

Results and Discussion
In our previous papers [20,21] it was assumed that the room-temperature 1 H NMR signals from a-Si:H films consist of two components associated with (i) hydrogen atoms belonging to Si-H bonds and (ii) H2 molecules entrapped in ellipsoid-like nanocavities. Therefore, the echo decays were fitted by a superposition of two exponentials. The longer relaxation time T2 was associated with the signal of the H2 molecules included in ellipsoid-like nanocavities. The angular dependence of this time was obtained [21,22] and analysis of this anisotropy allowed one to estimate parameters of the nanostructure [21,22].
At first, we fitted the experimental data by a two-exponential function,  7)). However, the characteristic decay times, which are usually considered as the relaxation times, differ by an order of magnitude, T2s<<T2 (Table 1). At 15 20 t − ms the signal is practically determined by the long-time decay processes which are characterized by T2. This allowed us to separate the slowly decaying signals attributed to mobile hydrogen molecules and further analyze only these processes. As it is seen from the inset in Fig. 2 done in semi-log scale, the obtained the long-time decays (at t>20 ms) differ from a single exponential function.

Analysis of the relaxation and determination of fitting parameters
A signal associated with hydrogen molecules located inside ellipsoid-like nanocavities was fitted by expression (7). The fitting parameters ( , standard deviations   and   , and averaged angles 0  and 0  ) were determined by means of the developed computer program based on Piton package, using the non-linear least squares minimization method.
The typical results of fitting of the long-time decays, associated with mobile hydrogen molecules, and the corresponding fitting parameters are presented in Fig. 3 and in Table 2   For other angles the obtained parameters were close to the values given in Table 2 Table 2). The reason for this may be that the relaxation rates in a single nanocavity possess angle symmetry in 3D space [18]. Another reason for the low accuracy in determining these parameters can be relatively weak dependence of the fitting results on the averaged angles. For example, at 0 S  = the function (7) does not depend on 0  (see Eq. 3). The fitting parameter  and standard deviations   and   agree well with these parameters obtained previously [21,22] by analyzing the angle dependence of the relaxation time T2. The approach developed here allows us to determine their values using a single time dependence of the signal.

Estimations of the structure parameters of nanopore
Let us estimate the shape and characteristics sizes of the nanopores using the results of fitting ( into Eq. (10) and assuming that according to [30][31][32] the spin density for H2 is 3 26 nm  − = we obtain, An experimental study of the nanostructure of a-Si:H films by other methods [27][28][29][30][31][32][33][34][35] showed that nanopores have an ellipsoidal closed to spherical shape. The method of growing our films was similar to that described by Gericke et al. in [35]. Assuming that the nanopores are spherical, the authors of [35] estimated the diameter of nanopores in the films deposited at a low deposition rate as 1.3 nm. Using this estimation and Eq. (11) we estimate the form factor F = 4.7943×10⁻².

CONCLUSION
We studied experimentally and theoretically the multicomponent room-temperature 1 H NMR echo decays in a-Si:H films. The anisotropy of the spin-spin relaxation of protons in H2 molecules entrapped in nanopores of these films was experimentally observed and theoretically explained. Despite the multicomponent signal decays, one of the fitting parameters, η, which is related to the size and shape of the nanocavities, is determined quite accurately. This made it possible to determine the size and shape of nanopores in a-Si:H films; our estimations are in good agreement with the results obtained by other methods [33,35]. Some characteristics of the film nanostructure were also determined with an acceptable accuracy. The experimental results are interpreted within the framework of the previously developed theory [14][15][16][17][18][19][20][21][22] in which a nanoporous material is represented as a set of nanopores containing liquid or gas and the relaxation rate is determined by the dipole-dipole spin interaction, taking into account the restriction of liquid or gas molecule moving inside the pores. Previously, the parameters and nanostructure characteristics such as averaged volume of the pores and their orientation distribution were determined from the angle dependencies of the spin-spin and spin-lattice relaxation times [21,22].
The main difference of the proposed approach is that the above-mentioned parameters can be obtained using a single signal. The method developed by us for the analysis of experimental NMR data is not based on analysis of anisotropy of the relaxation time T2, the determination of which is a rather complicated problem in the case of multicomponent decays.
Good agreement between the estimated pore sizes and the results obtained by other methods was achieved by determining several fitting parameters in expression (7). The fitting parameters also provides information about the nanostructure of the sample, such as the standard deviations for the angular distribution of pores, as well as the polar and azimuth angles of the average direction of the pore axes relative to the sample axis. The approach allows one to use the NMR data for quantitative non-destructive determination of the structural parameters of nanoporous materials.

Ethical Approval
We did not carry out any experiments or investigation involving human participants or living animals.