Application of GISAXS in the Investigation of Three-Dimensional Lattices of Nanostructures

: The application of the grazing-incidence small-angle X-ray scattering (GISAXS) technique for the investigation of three-dimensional lattices of nanostructures is demonstrated. A successful analysis of three-dimensionally ordered nanostructures requires applying a suitable model for the description of the nanostructure ordering. Otherwise, it is possible to get a good agreement between the experimental and the simulated data, but the parameters obtained by ﬁtting may be completely incorrect. In this paper, we theoretically examine systems having di ﬀ erent types of nanostructure ordering, and we show how the choice of the correct model for the description of ordering inﬂuences the analysis results. Several theoretical models are compared in order to show how to use GISAXS in the investigation of self-assembled arrays of nanoparticles, and also in arrays of nanostructures obtained by ion-beam treatment of thin ﬁlms or surfaces. All models are supported by experimental data, and the possibilities and limitations of GISAXS for the determination of material structure are discussed. and I.M.; analysis, M.M., and L.B.; data curation, K.S., M.M., L.B., and S.B.; writing—original draft preparation, L.B., and M.M.; review and editing, all authors.


Introduction
Grazing-incidence small-angle X-ray scattering (GISAXS) is a powerful method for fast and accurate structural analysis of thin-film materials containing nano-objects (NOs), like quantum dots, nanoparticles, nanowires, and similar nanostructures. It is essentially a diffraction method performed at grazing incidence angle to achieve relevant statistics when probing thin films, measured at small diffraction angles to capture objects of a few nanometers in size. The method is especially efficient for the analysis of systems that have some regularity in the NO ordering. Due to the regularity, distinctive Bragg peaks appear in GISAXS maps, because GISAXS reflects the structural properties in the reciprocal space [1][2][3][4][5][6][7][8][9][10][11][12]. Thus, a visual inspection of a GISAXS map is usually enough to roughly determine the type of the NO arrangement and the separations between the NOs. Of course, a detailed analysis and fitting of a GISAXS map is needed to get the precise values of the structural parameters including NO shape, size, arrangement properties, and their statistical distributions. Several excellent papers describe the GISAXS technique and the general theory behind it [1][2][3]13]. In addition, there are several very useful software platforms for the analysis of GISAXS data [14][15][16][17][18].
The main problem in the accurate analysis of GISAXS data is that each material type needs a specific model that describes the type of NO ordering inside it [3]. If the model is not properly chosen, it is possible to get good agreement between measured and simulated GISAXS maps, but with wrong values for the structural parameters. However, in some cases, the experimental data cannot be fitted by the NOs (this contribution is called shape factor), and by the arrangement of the NOs (called structure factor). The shape factor is the Fourier transform of the NO shape; thus, it is usually not difficult to define it. However, the structure factor depends on the NO ordering and different approximations have to be applied to get the function that properly describes the NO ordering. Therefore, we firstly concentrate on the description of the structure factor.
A GISAXS intensity distribution in general is given by the following expression [3]: where A is a constant, δρ is the difference in the electron densities of the NO material and the matrix, R and R' are the position vectors of the NOs, Q is the scattering vector, q is the complex scattering vector corrected to refraction at the vacuum-substrate interface, and Ω R (q) is the shape factor of a NO occurring atposition R; the shape factor is derived from the function that is unity in the NO volume and zero outside it; t i and t f are the Fresnel transmittivities [20] of the substrate surface corresponding to the primary and scattered waves, respectively ( t i t f 2 gives the Yoneda wing in a GISAXS map).
The brackets indicate averaging over the positions and shapes of the NOs. In this paper, we consider the case in which the NO sizes are not correlated with their positions (decoupling approximation); for details about different approximations, please see Reference [1]. As visible from Equation (1), a GISAXS map reflects the properties of the NO shape and ordering in the reciprocal (Q) space. The key parameter for the successful application of Equation (1) is to make good models for the NO shape (shape factor) and NO ordering (structure factor), as mentioned previously. The shape factor Ω(q) is the Fourier transform of the NO shape; therefore, if we know the shape precisely, it is usually not difficult to calculate it. However, for the determination of the structure factor, the summation should be performed over all pairs of the NO positions R and R' within the probed volume. This is the most important task in the determination of the structure factor, and, for that purpose, we must use approximations for the definitions of the NO positions, which make the calculation of the sum and its averaging in Equation (1) easier.
In this paper, we use a paracrystal model to describe the NO positions [3,21]. The paracrystal model for a one-dimensional (1D) case assumes that the position of the n-th NO (at the position R n ) can be expressed as a sum of basis vector a of an ideal (undisturbed) lattice and deviation vector: R n = na + D n , where D n denotes the total deviation of a NO with index n from its ideal position (na). The description of D n depends on the type of NO ordering, as described in detail in the next section.
Sometimes, extra attention should be given to the choice of the shape factor because it may be different than the assumed one. For example, some NOs make a chemical reaction with the surrounding material; thus, a kind of a shell appears around the NO and we have a different type of shape factor (core/shell instead of the NO alone). We consider the influence of the proper choice of the shape factor in Section 4.3.

Ordering Types in Nanostructured Materials
As explained in the text above using a 1D example (R n ), the positions of the 3D regularly ordered NOs deviate from the ideal positions that they would have in a perfect crystal lattice defined by basis vectors a 1 -a 3 . While the ideal positions define the position of peaks in the GISAXS map of a certain nanomaterial, the deviations are very important because they are closely related to the intensity distribution of each peak in the map. Therefore, the statistical properties of the NO deviations from the ideal positions (vectors D n ) are crucial for the successful analysis of a GISAXS map. The mathematical description of the statistical properties of the deviations is sometimes complicated. As mentioned before, we use a 3D paracrystal model to describe the NO positions and their deviations. This model treats each component of the basis vectors (a 1 -a 3 ) and the related deviations (D n 1 -D n 3 ) as independent variables. While the ideal positions are precisely defined with the basis vectors, the deviations are usually described by long-range ordering (LRO) or short-range ordering (SRO) models [3]. Additionally, we consider the possibility of a complete absence of regularity in ordering, which leads to the random ordering (RO) model. Figure 1 shows the main differences in these ordering types. We consider a 1D case to keep simplicity, but each component of a basis vector shows the same behavior. Thus, in random ordering, there is no correlation in the NO positions; thus, the structure factor for this type of ordering is a constant. On the other hand, the LRO and SRO assume the existence of correlation in the QD positions (regularity of the NO ordering) with the main difference that the NO positions fluctuate around a predefined separation in the SRO model, while they fluctuate around predefined positions in the LRO model. Therefore, if we fix the position of the first NO (indicated by number 1 in Figure 1b, the possible positions of the second NO have a distribution characterized by the width parameter σ around the ideal position (separation a). We assume it is a normal distribution with the standard deviation σ. The possible positions of the third NO in the SRO fluctuate around the ideal position (2a) but with the distribution having a deviation of 2σ (see Figure 1b). Actually, the separations between neighboring NOs fluctuate with the distribution described by σ, but the width of the distribution describing the position of each NO from the first NO (origin) consequently increases. The position of NO in such a 1D SRO model is given by ( , , ) ( ) = . (13) We assume that the basis vectors a1 and a2 are in the x, y plane, while the third basis vector may have all three components, depending on the ordering type. Therefore, only the equations related to the basis vector a3 ( ( ) (q);) have the refraction effects included (the others contain only the x and y components; thus, there is no need). Thus, for the successful analysis of a GISAXS map from a particular sample we have to construct the function G3D (Equation (14)) using 1D functions G(q) (Equations (7-13)). We must be careful to choose the correct functions for the description, because a GISAXS map strongly depends on the ordering type. Actually, we must know the main ordering properties of the NOs before starting the GISAXS analysis. These properties usually can be determined from the production process of the nanomaterial. If the material is produced by some kind of self-assembly process, then usually the separations fluctuate; therefore, we have to use the SRO model. However, if the positions are predefined from some reason (pre-pattering of the substrate, predefined layer thicknesses ( vertical positions of NOs), ion-beam-induced nucleation, ion tracks, etc.), the LRO model should be applied. In some cases, we have a combination of LRO and SRO related to one basis vector. Most usually, this appears in films composed on multilayers with predefined thicknesses, i.e., vertical positions, that contain NOs. If these NOs are ordered by some self-assembly process, then their lateral (x and y) components show SRO, while the z component has LRO. In that case, we use Equation (11), i.e., ( ) ( ), for the description of ordering along basis vector a3. We sometimes have to combine several models in different special directions, as illustrated in the examples given later. We can estimate the ordering type from the measured GISAXS map, if it shows more than one order of peaks related to the NO ordering. In the Q-space, the SRO and LRO differ by the behavior of the correlation peaks: the width of the peaks increases with the order of the peak for the SRO, while In the 3D case, we have basis vectors a 1 -a 3 and we assume that the x, y, z components (δ x;y;z ) of the random deviation δ n are normally distributed with zero mean and root-mean-square (r.m.s.) dispersion σ x , σ y , and σ z , respectively. Hence, to describe the system of NOs, we need to know the component of the basis vectors that characterize their ordering a 1−3 x,y,z (we interchanged the upper and lower indices of the standard notation a 1 -a 3 to include the space x, y, z components and to have the same notation as in Reference [3]), number of periods in each space direction N 1 -N 3 , the disorder parameters in a particular space direction σ 1−3 x,y,z , and, of course, the NO shape parameters. After insertion of Equations (2) or (3) into Equation (1), treating each component independently, the calculation of the sums and averaging is easy to perform analytically; thus, we come to the final equation for the GISAXS intensity distribution in the 3D case as follows: where f 1 and f 2 are the components of the shape factor, calculated as follows: where Ω FT (q) is the Fourier transformation of the shape function Ω R (q). F(q z ) and Γ(q z ) are defined as follows: where N 1 -N 3 are the numbers of the QD position (R n1,n2,n3 ) indices (number of periods) in the x, y and z directions, respectively.
G 3D (q) is a product of three 1D correlation functions G (i) ot (q); i = 1, 2, 3 and they depend on the ordering type in the particular space direction.
We assume that the basis vectors a 1 and a 2 are in the x, y plane, while the third basis vector may have all three components, depending on the ordering type. Therefore, only the equations related to the basis vector a 3 (G Thus, for the successful analysis of a GISAXS map from a particular sample we have to construct the function G 3D (Equation (14)) using 1D functions G(q) (Equations (7)-(13)). We must be careful to choose the correct functions for the description, because a GISAXS map strongly depends on the ordering type. Actually, we must know the main ordering properties of the NOs before starting the GISAXS analysis. These properties usually can be determined from the production process of Crystals 2019, 9, 479 6 of 13 the nanomaterial. If the material is produced by some kind of self-assembly process, then usually the separations fluctuate; therefore, we have to use the SRO model. However, if the positions are predefined from some reason (pre-pattering of the substrate, predefined layer thicknesses (vertical positions of NOs), ion-beam-induced nucleation, ion tracks, etc.), the LRO model should be applied. In some cases, we have a combination of LRO and SRO related to one basis vector. Most usually, this appears in films composed on multilayers with predefined thicknesses, i.e., vertical positions, that contain NOs. If these NOs are ordered by some self-assembly process, then their lateral (x and y) components show SRO, while the z component has LRO. In that case, we use Equation (11) MIX (q), for the description of ordering along basis vector a 3 . We sometimes have to combine several models in different special directions, as illustrated in the examples given later.
We can estimate the ordering type from the measured GISAXS map, if it shows more than one order of peaks related to the NO ordering. In the Q-space, the SRO and LRO differ by the behavior of the correlation peaks: the width of the peaks increases with the order of the peak for the SRO, while the peak width is constant for the LRO model and its intensity decreases with the order of the peak (for details, please see the Reference [3]). Hence, the peaks in the GISAXS map show the same basic behavior as the distribution of the NO positions when measured from the origin (Equations (2) and (3)). In Section 4, we examine three examples to illustrate the effects of ordering of NOs and their shape.

3D Ge Quantum Dot Lattices
In this first example, we consider a 3D lattice of spheroidal NOs (in this example, quantum dots (QDs)) prepared by multilayer deposition. We firstly consider the production process of that material. The material is prepared by computer-controlled deposition of a multilayer (in this example, alternating layers of Ge and Al 2 O 3 ); thus, the z position of each layer should be predefined (assuming a constant deposition speed). However, the QDs in this material form by diffusion-mediated nucleation, and their positions within the layers (x, y plane) are defined by a self-assembly process (see Reference [4]). Therefore, for this particular case, we should use a combination of LRO in the z direction and SRO in the x and y directions. The correct correlation function is The ordering type of QDs for this case is body-centered tetragonal (BCT), which can be concluded from the arrangement and shape of spots and also from the literature data about similar samples [10][11][12].
The measured GISAXS map from this material and its simulation are shown in Figure 2a,b. From the visual inspection of the measured map, we see only one lateral peak in the Q y direction and two peaks in the Q z direction. Here, we must be careful with the conclusion that there is no more than one peak in the Q y direction, and that the ordering is worse in that direction, because the intensity visible in GISAXS maps also depends on the shape factor. Where the shape factor is zero, no intensity will be visible in the GISAXS map; thus, the shape factor can "make" higher orders not visible.
The second important step in a GISAXS analysis is the inspection of the peak profiles. They are shown in Figure 2c,d. The profiles are taken along the lines shown in the inset of Figure 2d. They demonstrate the width of the profiles in the z direction. For both cases, the width of both peaks is nearly the same, showing that the chosen model was correct.
Therefore, when we are sure to have the suitable model of ordering, we can continue with the analysis of the intensity map. The third important step is to determine the lattice type of the NO ordering. This is sometimes difficult, and a complemental microscopy analysis should be applied to be sure that the ordering type is correct, because different ordering types can produce very similar GISAXS maps when they are averaged over azimuthal orientations (due to the appearance of random domains). Another possibility is to have some pre-knowledge about the nucleation properties and the Crystals 2019, 9, 479 7 of 13 production process of the nanomaterial and to make detailed simulations of different structures, to see the differences. One good example is given in Reference [6]. Therefore, when we are sure to have the suitable model of ordering, we can continue with the analysis of the intensity map. The third important step is to determine the lattice type of the NO ordering. This is sometimes difficult, and a complemental microscopy analysis should be applied to be sure that the ordering type is correct, because different ordering types can produce very similar GISAXS maps when they are averaged over azimuthal orientations (due to the appearance of random domains). Another possibility is to have some pre-knowledge about the nucleation properties and the production process of the nanomaterial and to make detailed simulations of different structures, to see the differences. One good example is given in Reference [6].
In this particular case, the ordering is body-centered tetragonal (BCT), as illustrated in Figure 3. Thus, for the description of an ideal lattice, we use the basis vectors a1-a3; a1 and a2 are placed in the x, y plane and have the same length |a1| = |a2| = a, while the vertical stacking of the QDs is defined by the vector a3 (see Figure 3a). As explained before, the multilayer has a predefined position of the single layers; thus, the z components of the QD positions (a3 z ; |a3 z |= c) have LRO, as illustrated in Figure 3b. The QDs within the layers show SRO (see Figure 3c), as they grow by a self-assembly process. It is also important to note that the regular ordering appears in domains randomly rotated around the z axis; thus, azimuthal averaging must be taken into account. The parameters of the analyzed QD lattice obtained by fitting of the experimental map from Figure 2a are given in Table 1, and the simulation is shown in Figure 2b. The shape of the QDs is assumed to be spheroidal with the lateral and vertical radii Rx,y and Rz, respectively, with the distribution having the standard deviation σR. We chose a gamma size distribution because it is very easy to get the analytical expressions for the averaged sizes using it. In this particular case, the ordering is body-centered tetragonal (BCT), as illustrated in Figure 3. Thus, for the description of an ideal lattice, we use the basis vectors a 1 -a 3 ; a 1 and a 2 are placed in the x, y plane and have the same length |a 1 | = |a 2 | = a, while the vertical stacking of the QDs is defined by the vector a 3 (see Figure 3a). As explained before, the multilayer has a predefined position of the single layers; thus, the z components of the QD positions (a 3 z ; |a 3 z |= c) have LRO, as illustrated in Figure 3b.
The QDs within the layers show SRO (see Figure 3c), as they grow by a self-assembly process. It is also important to note that the regular ordering appears in domains randomly rotated around the z axis; thus, azimuthal averaging must be taken into account. The parameters of the analyzed QD lattice obtained by fitting of the experimental map from Figure 2a are given in Table 1, and the simulation is shown in Figure 2b. The shape of the QDs is assumed to be spheroidal with the lateral and vertical radii R x,y and R z , respectively, with the distribution having the standard deviation σ R . We chose a gamma size distribution because it is very easy to get the analytical expressions for the averaged sizes using it.

3D networks of Ge Nanowires
The second example shows a 3D network of Ge nanowires grown also by the self-assembly growth presented in the previous example, but with no multilayer deposition. The deposition consisted of continuous Ge + Al2O3 co-deposition. In Figure 4, we present GISAXS maps and their simulations for two films differing by the NW length and the parameters of the NW lattice. The parameters of the simulations are given in Table 2. Details of their preparation are given in Section 2. For this case, there are no predefined positions for the network nodes; thus, all components of the G3D function are SRO, as shown below.

3D Networks of Ge Nanowires
The second example shows a 3D network of Ge nanowires grown also by the self-assembly growth presented in the previous example, but with no multilayer deposition. The deposition consisted of continuous Ge + Al 2 O 3 co-deposition. In Figure 4, we present GISAXS maps and their simulations for two films differing by the NW length and the parameters of the NW lattice. The parameters of the simulations are given in Table 2. Details of their preparation are given in Section 2. For this case, there are no predefined positions for the network nodes; thus, all components of the G 3D function are SRO, as shown below.
G 3D (q) = G  However, the complication in this example is the shape of the NWs. The network is continuous, and its geometry is shown in Figure 5. In an ideal case, this structure consists of cylinders connected in the nodes of a BCT lattice, as shown in Figure 5a. In the real case, there are deviations in the network node positions, as visible from the cross-sectional TEM image of this structure (Figure 5b).  However, the complication in this example is the shape of the NWs. The network is continuous, and its geometry is shown in Figure 5. In an ideal case, this structure consists of cylinders connected in the nodes of a BCT lattice, as shown in Figure 5a. In the real case, there are deviations in the network node positions, as visible from the cross-sectional TEM image of this structure (Figure 5b). The form factor of this shape has a complicated form, and we prefer to have simple analytical formulas for GISAXS fitting, because we wish to have a fast fitting program. Therefore, we use an approximation choosing spheroidal QDs arranged in a BCT lattice (the same as in the previous example), but with the QDs tilted along the basis vector a 3 as shown in the Figure 5c. R x,y denote the radii of the NWs, and R z is the half of their length. Here, we have to perform averaging over all four equivalent configurations of the basis vectors a 3 (connection of the center of the BCT lattice to the four neighboring nodes in the a 1 -a 2 plane), as well as azimuthal averaging. However, the complication in this example is the shape of the NWs. The network is continuous, and its geometry is shown in Figure 5. In an ideal case, this structure consists of cylinders connected in the nodes of a BCT lattice, as shown in Figure 5a. In the real case, there are deviations in the network node positions, as visible from the cross-sectional TEM image of this structure (Figure 5b). The form factor of this shape has a complicated form, and we prefer to have simple analytical formulas for GISAXS fitting, because we wish to have a fast fitting program. Therefore, we use an approximation choosing spheroidal QDs arranged in a BCT lattice (the same as in the previous example), but with the QDs tilted along the basis vector a3 as shown in the Figure 5c. Rx,y denote the radii of the NWs, and Rz is the half of their length. Here, we have to perform averaging over all four equivalent configurations of the basis vectors a3 (connection of the center of the BCT lattice to the four neighboring nodes in the a1-a2 plane), as well as azimuthal averaging.  Using the presented model, we get a good agreement between the experimental and the simulated intensity distributions. The parameters obtained are also in good agreement with the microscopic measurements [22].

Lattices of Core/Shell Quantum Dots
In this example, we demonstrate how the internal structure of QDs can influence the analysis of a GISAXS map. The production process of the core/shell QDs presented in this example is very similar to that in Section 4.1, except for the fact that a very thin Si layer is deposited after each Ge layer, i.e., it is a Ge/Si/Al 2 O 3 multilayer deposition (for details, please see Section 2). Due to such a deposition process, each Ge QD is covered by a thin (about 1 nm thick) Si shell, and this shell is not symmetric around the core. This process is analyzed in detail in References [10,11]. Hence, the shape and ordering type of the QDs is the same as in Section 4.1, but the internal structure of the QDs is slightly different. It seems reasonable that a different internal structure of the QDs leads to different GISAXS intensity distributions, and no additional story about it is necessary. However, here, we want to point out that a very tiny difference in the QD structure, like a shell with the thickness of a few angstroms that is often invisible by microscopy, can cause very strong effects and cause problems in the GISAXS analysis (please see Reference [10] for an example). In many cases, we are not even aware that such a shell exists around NOs that we wish to investigate, because a chemical reaction easily occurs at the interface between the NO and the matrix that surrounds it to form a shell. Therefore, this example illustrates how a thin shell affects a GISAXS intensity distribution.
The correct choice of the shape factor is a Fourier transform of the core-shell nanoparticle with a shifted origin (details about the structure can be found in References [10,11]). The experimentally measured map and the simulation obtained by fitting using that shape factor are given in Figure 6a,b, respectively. However, if we apply the same model as in Section 4.1 (BCT ordering of spheroidal QDs), we can also get a reasonably good fit and good agreement between the experiment and its simulation (Figure 6c). The parameters for both simulations are given in the Table 3. We can see that the obtained sizes for the QDs, especially their vertical sizes, are completely wrong. Only after microscopy analysis of the material cross-section, shown in Figure 6d, the correct model was identified. The simulation with the core-shell shape factor gives a very good agreement with the microscopy measurements (Figure 6e). However, the sizes obtained by the simple model are very different from the real values (Figure 6f). With a careful inspection of the simulation Sim 2, we can notice that the low intensities (at the higher Q z ) are different from in the experiment. Here, we see it because of the logarithmic scale for the intensity distribution; for the linear case, it could be easily overlooked. If we use larger radius values, then we cannot get any good agreement for the peak shape and intensities. Thus, a very careful analysis is needed. The effect depends strongly on the material of the core and shell, i.e., on their electron density contrast. Generally, the strongest difference should have a stronger effect on the GISAXS intensity distribution.  Table 3. Parameters of the two simulations shown in Figure 6b,c for the film GeSi1; t is the shell thickness, and d is the shift of the core with respect to the center of the shell.    Figure 6b,c for the film GeSi1; t is the shell thickness, and d is the shift of the core with respect to the center of the shell.

Discussion
In the first example (given in Section 4.1), it was illustrated how to distinguish LRO and SRO in the particular space directions and how the ordering type influences the peak profiles. To be able to distinguish the SRO and LRO types of ordering, we should have at least two orders of peaks visible in the GISAXS maps. Of course, we can also determine it from the preparation process of the material if we know it. For example, if the entire material grows only by a self-assembly process, like for the Ge NWs given in the second example (Section 4.2), the SRO model should be applied in all directions. Often, we have a combination of different ordering types in different directions, as also shown in Section 4.1, and several nice examples showing combinations of all three ordering types (SRO, LRO, and RO) are given in Reference [2]. In that paper the analysis of GISAXS maps is considered from the nanostructures produced by ion modification of material surfaces. Ordered NOs form along the trace of each ion on the material surface. These NOs have SRO along the ion track, and LRO in the direction perpendicular to it, due to the predefined positions along the ion track. If the irradiation fluence (number of used ions per surface unit) is low, then these traces have RO; however, if the fluence is high, the traces overlap and a superstructure of ion traces is formed which has SRO type.
From the third example (Section 4.3), we can conclude that a very wrong result can be obtained even if the fit results are reasonably in good agreement with the experiment. The problem is of course the shape factor. This is one of the main disadvantages of the GISAXS method: different structures can produce the same or very similar GISAXS maps, and if we do not know enough details about the material, and if we do not have complemental measurements, it is very easy to get completely wrong results.
Finally, we wish to discuss the central region of the GISAXS maps, close to the Q y = 0 axis. This region contains a contribution of the coherent scattering that is not taken into account in the used simulations. Additionally, it also contains the contribution from the interface roughness in the multilayer samples, (not taken into account in the simulations). If we take into account all these contributions, the number of fitting parameters will be huge, and the fit would be completely unreliable. Therefore, in the GISAXS analysis, we do not take this central part into consideration. In addition, the paracrystal-based models that we use for the simulations do not work well in this region because the shape effects of the entire QD lattice become significant there. Therefore, cross-like structures appear in this part in the simulations, instead of sheets (please see the GISAXS simulations in the given examples). Only the central position (at the Q y = 0 plane) of this cross-structure is correct and we should take it into account. An illustration of how to select properly the region of interest for GISAXS is shown in the Figure 7.
As mentioned in the introduction, a new computer program called GisaxStudio [19] is under development. It enables the selection of the region of interest of the desired shape for the fitting. It contains the different models mentioned here (3dLattice for 3D lattices of nanoparticles) and also in Reference [2] (iBeam for ion-beam-modified materials). We hope that it will not be useful only for researchers for whom the structural characterization of nanomaterials is the primary research interest, but also for everyone who is just "passing by". The software is also designed to grow; thus, it can be used for different types of materials in the future. The latest version is available in Reference [19], and a paper describing it precisely is in the preparation process [23]. addition, the paracrystal-based models that we use for the simulations do not work well in this region because the shape effects of the entire QD lattice become significant there. Therefore, cross-like structures appear in this part in the simulations, instead of sheets (please see the GISAXS simulations in the given examples). Only the central position (at the Qy = 0 plane) of this cross-structure is correct and we should take it into account. An illustration of how to select properly the region of interest for GISAXS is shown in the Figure 7. As mentioned in the introduction, a new computer program called GisaxStudio [19] is under development. It enables the selection of the region of interest of the desired shape for the fitting. It contains the different models mentioned here (3dLattice for 3D lattices of nanoparticles) and also in Reference [2] (iBeam for ion-beam-modified materials). We hope that it will not be useful only for researchers for whom the structural characterization of nanomaterials is the primary research interest, but also for everyone who is just "passing by". The software is also designed to grow; thus, it can be used for different types of materials in the future. The latest version is available in Reference [19], and a paper describing it precisely is in the preparation process [23].

Conclusions
In this paper, we showed some basic principles of GISAXS data analysis. It shows the importance of the proper choice of models for the description of the arrangement and shape of nano-objects within a matrix. The correct model type can be determined by careful analysis of the nano-object formation process, and from the experimental GISAXS maps from the materials containing nanoobjects. We also pointed out some main steps in the analysis of GISAXS data. It is especially suitable for beginners in the usage of the GISAXS technique.

Conclusions
In this paper, we showed some basic principles of GISAXS data analysis. It shows the importance of the proper choice of models for the description of the arrangement and shape of nano-objects within a matrix. The correct model type can be determined by careful analysis of the nano-object formation process, and from the experimental GISAXS maps from the materials containing nano-objects. We also pointed out some main steps in the analysis of GISAXS data. It is especially suitable for beginners in the usage of the GISAXS technique.