Practical mineralogical quantification of bentonites supported for a PXRD calibrated hkl model

In the bentonite industry — exploration and beneficiation — ore quality control has been traditionally carried out using the swelling parameter, denominating those that swell as naturally sodium bentonites and those that do not as naturally non-sodium bentonites. However, the properties of the bentonites can vary due to the quantity and cationic type of montmorillonite. The variation in those properties cannot always be predicted, because there is still no consolidated and convenient method for clays mineralogical quantification. A quantification via PXRD using the Rietveld method is only reliable when all the crystalline structures of the mineral phases are known. This is not the case for bentonites, since the montmorillonite shows turbostratic disorder that leads to extremely wide and asymmetric non-basal reflections that are not considered in the available structural models. Thus, in this work, is generate a calibrated hkl model for a Brazilian montmorillonite applying the combined Rietveld-Le Bail-Internal Standard method developed by Paz et al . (2018). The study concluded that the combined method is a good choice for the mineralogical quantification of the bentonites, handling turbostratic disorder of montmorillonite, since the develop hkl phase model showed good results in mixtures with >50% montmorillonite content (greater precision and reproducibility).


INTRODUCTION
The term bentonite was first used in 1898 by W.C. Knight to describe a plastic clay located in Benton shale, Rock Creed Region, Wyoming, USA, with the capacity to adsorb large amounts of water and forming thixotropic suspensions (Eisenhour and Brown 2009). Geologically, bentonites are clay rocks composed essentially by Na-montmorillonite, mainly formed from the chemical alteration of volcanic glass of pyroclastic and/or volcanoclastic rocks (Christidis and Huff 2009). However, due to the lack of natural occurrences of bentonites like the North American, the term is generally used for a clay with similar properties consisting essentially of 60-80% of montmorillonite, which are commercialized regardless of occurrence or geological origin (Grim 1973, Elzea andMurray 1994).
Montmorillonites are dioctahedral 2:1 clay minerals of the smectite group. They are end-members of the montmorillonite-beidellite series with mainly octahedral charge generated by substitution of Al 3+ or Fe 3+ by Mg 2+ and layer charge of 0.2 to 0.6 mol/f.u. ). The features of the smectite group such as the large variability in chemical composition, variable layer charge, octahedral and tetrahedral charges and the Fe content ) can lead to an unambiguously classification of the dioctahedral types.
According to these features, the Formosa Mg-montmorillonite, a relatively new occurrence from northern Brazil (Paz et al. 2011(Paz et al. , 2012a, which is of interest in this work, is classified as a beidellitic montmorillonite with octahedral layer charge of > 50%, Fe content of < 0.30 mol/f.u., and layer charge between 0.426-0.6 mol/f.u. Bentonites are widely used as industrial raw material in a large number of technical applications, due to the physical and chemical properties of montmorillonites that underlie its commercial value. Their major uses are as a thixotropic agent in drilling fluids, iron pelletizing, foundry processes (metal casting binder), filtering, clarifying, decolorizing, odor absorbent, cat litter, paper industry, civil engineering, waterproofing, and many other applications (Odom 1984).
Since the end of 20 th century, the uses for bentonite and its consumption have increased dramatically, requiring a special attention to ore quality control. In the bentonite industryexploration and beneficiation -ore quality control has been traditionally carried out using the swelling parameter, denominating those that swell as naturally sodium bentonites and those that do not as naturally non-sodium bentonites. However, the properties of the bentonites can vary widely due to the quantity and cationic type of montmorillonite, as well as the amount of accessory minerals present. These variations cannot always be predicted, because there is no consolidated and convenient method for the mineralogical quantification of clays.
Mineralogical quantification methods are still under development, including imaging analytical techniques and Powder X-ray diffractometry (PXRD) (Scarlett andMadsen 2006, Paz et al. 2018). The latter is considered as the most appropriate routine quantitative analysis method compared to other analytical techniques, such as Fourier Transform Infrared spectroscopy (FTIR), chemical analysis, and electron microscopy (Bish andPost 1993, Zhou et al. 2018). However, there are still many uncertainties related to the PXRD technique, even for simple geological matrices (Paz et al. 2018), having many ways to quantify without a reliable method, obtaining results with great variability.
The Rietveld refinement (Rietveld 1969) is one of the most commonly used methods to quantify polycrystalline mixtures via PXRD. This method allows dealing with a wide variety of sample-related effects by using the full diffraction pattern as a whole (Hill and Howard 1987, Scarlett and Mandsen 2006, Paz et al. 2018. Nevertheless, this is only reliable when all the crystalline structures of the mineral phases within the sample are known (Paz et al. 2018, Taylor andRui 1992). In addition, an important requirement of the method is that the diffraction patterns are formed by well-defined Bragg reflections (Wang et al. 2018). This is not the case for bentonites, since the montmorillonite shows turbostratic disorder: random rotation and/ or translation of individual layers to each other (Biscoe and Warren 1942, Ufer et al. 2008, Wang et al. 2018, Taylor and Matulis 1994, which causes a broad diffraction band as a result of the union of several neighboring reflections (Paz et al. 2018). This disorder leads to extremely wide and asymmetric non-basal reflections (Ufer et al. 2008) that are not considered in the available Crystallographic Information File of the Inorganic Crystal Structure Database (CIF-ICSD) structural models, having only partial crystallographic information. Additionally, the variety of interchangeable cations in the interlayer space and the isomorphous substitutions occurring in the octahedral and/or tetrahedral sheets affect the position of characteristic reflections during quantification (Zhou et al. 2018).
Due to these characteristics, the diffraction patterns of montmorillonites are more complex than the patterns of crystalline materials (Zhou et al. 2018), making it difficult to quantify the phases using the conventional Rietveld method, which makes the mineralogical quantification of bentonites a major challenge.
An alternative for cases where the crystallographic information is partial has been to combine three widely known and reliable methods used for mineralogical quantification: Rietveld, Le Bail and Internal Standard, prioritizing speed, convenience and quality. Thus, the objective of this work was to generate a calibrated hkl model for the Formosa Mg-montmorillonite applying the combined Rietveld-Le Bail-Internal Standard method developed by Paz et al. (2018) for the quantification of the main mineral components of bentonite clays.

Materials
High purity fluorite and concentrated Mg-montmorillonite samples (Tab. 1) were used to prepare the standard mixtures presented in Table 2. The Mg-montmorillonite sample was obtained by particle size separation and concentration of the clay fraction (< 2 μm) from the Formosa bentonite, described for the first time by Paz et al. (2012a).
The occurrence of this bentonite is associated with altered Mesozoic volcanism basalts of the Mosquito formation, from the Parnaíba Sedimentary Basin, in Formosa da Serra Negra County, southern Maranhão State, northern Brazil (Paz et al. 2012a).
The clay fraction (< 2 μm) was obtained according to the following procedure: 50 g of the material was dispersed in 200 mL of deionized water in ultrasonic immersion for 30 min. Then, the material was wet sieved manually to separate the sediment sand fraction (> 75 μm, 250 mesh). Afterward, the suspension was centrifuged for 5 min at 1,400 rotations per minute (RPM) in a Novatecnica NT820 centrifuge to separate the sedimented silt fraction (75-2 μm) from the clay fraction (suspended in the overflow). The overflow was further centrifuged for 10 min at 2,800 RPM to concentrate the clay fraction (< 2 μm) and then was dried in an oven at 60°C. Finally, the concentrated clay fraction was manually pulverized in an agate mortar and to ensure the homogeneity of the sample, quartering was conducted according to the elongated heap method.
Fluorite sample is a reference material produced and certified by the Sigma Aldrich laboratory, Germany. It was used as an internal standard in the calibration of the hkl model and in the standard mixtures for its validation.

Preparation of standard mixtures
A series of 7 binary standard mixtures of montmorillonite-fluorite was prepared, weighted in accordance with the proportions in Table 2. For the calibration of the hkl model, three replicates of the MF8020 standard mixture were prepared. A quantity of 3.0000 ± 0.0006 g of each replicate was prepared using the composition of two parts of 1.5000 ± 0.0006 g. For validation, 1.5000 ± 0.0006 g of each standard mixture were prepared. The mixtures were manually pulverized and homogenized in an agate mortar.
Proposed methodology: the combined Rietveld-Le Bail-Internal Standard Powder X-Ray Diffraction Method This method was initially proposed by Paz et al. (2018), which combines three widely known and trustworthy methods used for mineralogical quantification of PXRD results: Rietveld, Le Bail and Internal Standard, prioritizing speed, convenience, and quality in the quantification.
According to the proposed methodology, the first step is to generate a calibrated hkl phase model for montmorillonite, using the Le Bail fitting profile method. Then, the calibrated hkl model is validated using the Rietveld refinement method, inserting simultaneously the calibrated hkl phase model for montmorillonite and the CIF-ICSD file for fluorite (internal standard).

Equipment
PXRD measurements were performed using a PANalytical Empyrean divergent beam diffractometer, with a θ-θ goniometer, a sealed ceramic Co X-ray tube (Kα1 = 1.78901 Å), with a line focus of 1,800 W, a Fe kβ filter, and a PIXcel 3D area detector, operating in a linear scanning mode (1D), with active length of 3.3473° 2θ (255 active channels).

Instrumental conditions for PXRD measurements
The samples were measured using the following instrumental conditions: 40 kV, 35 mA, soller slits of 0.04°rad (incident and diffracted beams), scan range from 2 to 110° 2θ, step size of 0.0066° 2θ with 19.266 s of time/step in the continuous scanning mode, divergence slit of ¼ o , anti-scatter slit of ½ o , irradiated sample size of 10 mm, anti-scatter slit for the diffracted beam of 7.5 mm (PIXcel), and sample spinning 2 rotations/s. The total time of the analysis was approximately 34 min. Data acquisition was carried out with the PANalytical X'Pert Data Collector software version 5.1a. Instrumental resolution was determined using the LaB 6 NIST/SRM 660b standard.

Software
The mineralogical identification was carried out with the PANalytical X'Pert HighScore Plus software, version 3.0, using the Crystallography Open Database (COD) (Gražulis et al. 2009).
Mineralogical quantification by applying the combined method was performed using the FullProf software, version 2011 (Rodríguez-Carvajal 1993), with the specific graphical interface Full version 2012 (Paz et al. 2012b). The CIF files for montmorillonite, fluorite, and LaB 6 are referenced in Table 3.

Chemical analysis
The X-ray fluorescence (XRF) spectrometer employed was a wavelength dispersive (WDS) sequential Malvern PANanalytical Axios Minerals model, with a Rhodium (Rh) ceramic X-ray tube operating at a maximum power potential level of 2.4 kW. The specimen was fused with lithium tetraborate-Li 6 B 4 O 7 (1:6 g sample/flux). Data acquisition was performed using the PANalytical SuperQ Manager software. For the loss on ignition (LOI) analysis, samples were dried at 105°C, followed by calcination at 1,020°C for 2 hours in a muffle furnace.

Mössbauer spectroscopy
The spectra were recorded at room temperature (298 K) in a previously pulverized sample, using a conventional spectrometer operating in a constant acceleration mode with triangular reference signals, and 57 Co source diffused in a Rh layer. Data accumulation was performed with a multichannel analyzer with 1,024 channels and a speed range from -11 to +11 mm/s. Velocity range was calibrated with respect to the α-Fe metallic spectrum collected.

Scanning electron microscopy
A Zeiss SIGMA-VP scanning electron microscope was used for SEM analysis, operating with a constant accelerating voltage of 20 kV. The sample was sprayed over a holder with a double-sided carbon tape and then coated with gold using an Emitech K550X sputter coater.

Fourier transformed infrared spectroscopy
The Fourier transformed infrared spectroscopy (FTIR) spectra were recorded in a Thermo Scientific Nicolet Is50 FTIR spectrophotometer, in the spectral range of 4,000-400 cm -1 and 100 scans were recorded with a frequency resolution of 4 cm -1 . The sample was prepared as pressed KBr discs (1:200 mg sample/ KBr).

Particle size distribution
The measurements were conducted in a Malvern MASTERSIZER Hydro 2000MU particle size analyzer, whit measured size range of 0.020 to 2,000.000 μm. The parameters used were: spherical particle model, rotation of 2,500 RPM, dispersant volume of 800 mL, dispersant refraction index (water) of 1.330, particle refraction index (montmorillonite) of 1.560 according to Mie's theory, ultrasonic immersion time of 1min, and obscurity range of 5-10%.

Chemical analysis
The chemical composition of the Formosa montmorillonite is shown in Table 4. The main components present are SiO 2 and Al 2 O 3 , dominant constituents of all clay minerals. The low percentage of TiO 2 corresponds to Anatase, identified by PXRD analysis. It is important to emphasize that the high percentage of Fe 2 O 3 , typical for most Brazilian bentonites, corresponds to the total Fe content. However, Fe is related to at least two different minerals identified by PXRD: montmorillonite and hematite. Typically, in montmorillonite, Fe 2+ and Fe 3+ can substitute Al 3+ in the octahedral sheet ). Thus, the Mössbauer spectroscopy was used to distinguish between the Fe related to those two minerals.

Mössbauer spectroscopy
The Mössbauer spectrum consists of three spectral components: two overlapping Fe 3+ doublets and one sextet (Fig. 2).
The inner broad doublet with an isomer shift of 0.38 mm/s and a quadrupole splitting of 1.06 mm/s is relative to the structural Fe 3+ in octahedral position of the montmorillonite. The outer doublet that is much more intense than the inner doublet is related to impurities of amorphous character (with no longrange magnetic ordering) (Cao et al. 1997). Such impurities are present at quantities too low for detection by PXRD analysis. A narrow sextet is attributed to hematite, with a hyperfine magnetic field of ~51.1 T (Murad and Wagner 1998). The quadrupole splitting of -0.22 mm/s suggests a weakly ferromagnetic state, isomorphic substitutions of < 10% and/or particle sizes of > 20 nm (Murad 2010). This mineral was identified at low quantities by PXRD analysis. The Isomer shifts (δ), quadrupole splitting (D), and peak widths together with the assignment of different components and its percentages of relative spectral area (RA) are summarized in Table 5.

Structural formula calculation
To ensure that the chemical analysis represents only the montmorillonite's composition, XRF data were corrected using the calibration of the Fe 2 O 3 result obtained from the Mössbauer analysis. The structural formula calculated by the XRD: X-ray diffractometry; Mnt: montmorillonite; Kln: kaolinite; Ant: anatase; Kfs: K-feldspar; Hem: hematite.      (Schultz 1969). The total layer charge is -0.5 mol/f.u., which is balanced by the charge on the interlayer cations (0.5 mol/f.u.). Substitutions of 0.30 Mg 2+ and 0.06 Fe 3+ for Al 3+ in the octahedral sheet accounts for 60.16% of the total layer charge. The 39.84% of the total layer charge is due to 0.20 Al 3+ for Si 4+ substitution in the tetrahedral sheet. Fe 3+ < 0.3 mol/f.u. and Al > 1.4 mol/f.u. indicate that the sample is cis-vacant ). According to Emmerich et al. (2009), smectite can be classified as a highly charged beidellitic montmorillonite (Tab. 6). Figure 3 shows the SEM micrographs of the Formosa montmorillonite with typical smectite morphology: very extensive stacking of the sheets, evidenced by large undulations and tactoids, tactoids inside the aggregates, and arrangement of  aggregates inside the clusters (Souza Santos 1989, Neaman et al. 2003, Paz et al. 2011.

Particle size distribution
Particle Size Distribution (PSD) curves (Fig. 5) show a particle size range from ~0.3 to ~100 μm and ~60% (by number) of the particles are < 10 μm (ideal particle size). The frequency distribution curve exhibits a monomodal lognormal distribution with the local maxima at 8.9 μm and ~54% particles are  finer than this value. This distribution curve together with the coefficient of uniformity (> 4) and the coefficient of curvature (~1) indicate a well-graded sample with a homogeneous particle size distribution (ASTM 2011).
The average particle size of 13.08 μm and the D 90 > 28.34 μm, suggest an acceptable particle size distribution for the mineralogical quantification, considering that all particles are < 50 μm (Payzant 2008). These values show that the preparation method (manual pulverization-homogenization) was excellent, allowing to reduce errors associated to particle size. Table 7 summarizes the PSD parameters. Calibrated hkl phase model for montmorillonite Figure 6 shows the diffraction pattern of the Formosa Mg-montmorillonite with the lines of a similar montmorillonite indexed standard (ICSD-CIF 159274).
Significant discrepancies were observed in the relative intensities and in the peak positions when comparing both diffraction patterns. These differences reveal a lack of crystalline structure information related to structural defects and divergences in unit cell dimensions and spatial groups (Paz et al. 2018).
The Formosa Mg-montmorillonite exhibit a turbostratic stacking disorder which causes a broad diffraction band in the region 22-37 °2θ CoKα (Highlighted region in Fig. 6). It is the most severe kind of layer stacking disorder, defined as a random rotation and/or translation of the individual layers tetrahedral-octahedral-tetrahedral (TOT) to each other in the a-b plane (Fig. 7) (Biscoe and Warren 1942, Ufer et al. 2009, Wang et al. 2018. The broad diffraction band region shows a typically broadened peak with strong asymmetry (Ufer et al. 2009, Wang et al. 2018, without the complete design of each of them up to the background line (Angélica et al. 2018, Paz et al. 2018. This effect difficults the use of the conventional Rietveld method, as the method allows dealing with phases which contain only minor disorders. However, it presents severe problems of adjustment when dealing with severe disorders such  as turbostratic disorder, preventing the correct determination of the phase content (Ufer et al. 2009, Ufer et al. 2008. For a successful refinement by the Rietveld method, it is necessary that the structural model must be as similar as possible to the measured diffraction pattern. This is not the case of montmorillonite, because the complex region of the diffraction band is not adequately considered in the available structural models. Moreover, the montmorillonite is characterized by a small tactoid size, which combined with the turbostratic disorder make it difficult to develop a proper structural model in terms of atomic coordinates (x, y, z), with reliable thermal and occupational factors (Dermatas and Dadachov 2003). Nevertheless, a successful quantitative analysis for mineralogical phases with partial crystallographic information is possible by combining the Rietveld, Le Bail, and Internal Standard methods (Paz et al. 2018).
Concerning the Le Bail profile fitting (Le Bail et al. 1988), an hkl model is developed from a prepared mixture with known amounts of the phase of interest plus a well-characterized standard material (Paz et al. 2018). In this method, peak positions of the phase of interest are limited according to the space group, while allowing variations in the peak individual intensities and the scale factor that represents an empirical structure factor related to the concentration of the phase is derived (Scarlett andMadsen 2006, Paz et al. 2018). This hkl model is used in the Rietveld refinement for quantitative phase analysis.
The ICSD-CIF 159274 (Gournis et al. 2008) structural model was the most compatible with the Formosa Mg-montmorillonite among the models available in the ICSD database. Thus, unit cell dimension data and the space group were used to generate the calibrated hkl model. Accordingly, using the Le Bail method, an hkl model for montmorillonite was generated and calibrated from the standard binary mixture MF8020 (montmorillonite 73.7%, fluorite 19.8%, impurities 6.5%). An average scale factor of 2.96x10 -4 for montmorillonite was manually adjusted at the end of the refinement, and the average scale factor of 4.07x10 -4 for fluorite was calibrated from the known mass fraction. The resultant hkl file contains a list of 195 effective reflections, with the information of the hkl planes, multiplicity, peak position, and absolute intensities (Fig. 8).
A good fit of the model (Fig. 9) was obtained when the initial interval of 22.4-23.1 °2θ CoKα of the diffraction band was excluded. The background was adjusted by the adjustable points method, with the manual insertion of 15 points in free positions, and the profile was adjusted using the Thompson Source: Gournis et al. (2008). XRD: X-ray diffractometry; CIF: crystallographic information file; ICSD: Inorganic Crystal Structure Database.  Cox Hastings function (TCH ), without considering the structural information (atomic coordinates), which significantly reduces the number of parameters to be refined.
The method was reproducible for three replicates of the MF8020 mixture (A, B and C) (Tab. 8), obtaining statistical indices of refinement χ 2 = (R wp /R exp ) 2 and R Bragg lower than 5, which were considered satisfactory (Tab. 9). The stability and greater precision along with the speed and convenience of the results show that the combined method is a good choice to be considered in the industry for the mineralogical quantification of bentonites.

Summary of refined parameters in generating the calibrated hkl phase model for montmorillonite
For the construction of the hkl phase model, the Pseudo-Voigt and TCH peak profile functions were tested. According to the graphical adjustment, the visual inspection of the residue and the statistical errors (R Bragg , R wp , R exp and χ 2 ), the TCH function was the best option to refine the peak shape of the montmorillonite and fluorite. It gave a better matching to the experimental diffraction lines, without extending the convergence time of this multi-parametric system (Gournis et al. 2008). This function allows analyzing the tactoid-size and strain broadening effects (Thompson et al. 1987).
To ensure greater stability of the method a parameter was included at each step. In each of the refined parameters, the iterative adjustment calculation was executed for several cycles until full convergence.
The hkl model was generated according the following best refinement sequence: starting with the Rietveld method for the two phases (montmorillonite and fluorite). The first parameter refined was the sample displacement. After refining that,  the relaxation factors for the profile fitting, and atomic and global parameters were reduced from 0.99 to 0.2. The background was initially adjusted by the polynomial function with five coefficients.
Subsequently, the small region at the beginning of the montmorillonite diffraction band was excluded (22.4 to 23.1° 2θ CoKα). Removal this interval means that only the upper part of the intensities is considered in the hkl model, and the lower part is considered as a background (Paz et al. 2018).
Afterward, the polynomial function of the background was changed to a linear interpolation with given points, with the simultaneous manual insertion of 15 points in the diffraction pattern at peak-free positions, which were not refined initially. The insertion of the adjustable points allows to consider the background as a discontinuous function (two continuities). This step is necessary due to the removal of the small region at the beginning of the diffraction band.
The first refined phase is the one with the highest concentration, in this case, montmorillonite. Due to the small tactoid size and the turbostratic disorder, the montmorillonite phase exhibits cell unitary defects and asymmetric and broad peaks.
Therefore, beginning with the Rietveld method, unit cell parameters (a, b, c, β) were the first parameters to be refined. Afterward, the Y factor in the TCH function that corresponds to the tactoid size was refined. The third largest effect of the montmorillonite phase on the diffraction pattern was the peak asymmetry that was refined in the TCH function (two asymmetry factors). Consequently, the montmorillonite refinement was changed from the Rietveld refinement to fitting profile Le Bail method (profile matching S = Cte, in FULL interface) and then, the unitary cell content was informed (calculated structural formula).
As regards fluorite, no significant discrepancies were observed between the diffractometric profile observed and the structural model used (CIF-ICSD 28730). However, a characteristic feature of fluorite is the broad peaks caused by microstrain. Therefore, this was the parameter refined, adjusted by the X and U parameters of the TCH function and two asymmetry factors.
Subsequently, the background was refined, leaving the extreme points of the two continuities fixed. Finally, with the best-fitted diffraction model, the montmorillonite scale factor was manually modified until a standard mass value of montmorillonite equal to 73.7% was obtained. Then, the *.hkl file (calibrated hkl phase model for montmorillonite) was created (Fig. 8).  Validation of the calibrated hkl phase model for montmorillonite The method was reproducible for three replicates of the binary mixtures MF8020. The statistical indices of the refinement were satisfactory since the global index χ 2 for each mixture was ~1 and the R Bragg of the two phases was ≤ 5 (Tab. 10). Figure 10 represents the residual plot for each quantified mixture (residual difference between the standard mass and estimated mass values). A random distribution and the absence of trend behavior observed in the mixtures with > 50% of montmorillonite, show statistical independence of the residues. The difference between the calculated mass fraction and the standard mass fraction in these mixtures was below ± 2.5%. The global indices χ 2 for each mixture were considered satisfactory since the values were < 5. As regards R Bragg indices for each phase, the values tend to be high when the phase concentration is low, not only because of the lack of fit (Tab. 11) (Paz et al. 2018).
Regarding the binary mixtures with < 50% of montmorillonite, an underestimation of the expected value was observed (Tab. 10). This is due to the fact that the number of effective reflections decreases as the phase concentration decreases, affecting the success and stability of the refinement, without having a full convergence in each of the refined parameters.

Summary of the refined parameters in the quantification of the montmorillonite and fluorite binary mixtures
The refinement sequence for the quantification of the binary mixtures was the same used in the generation of the hkl model, but, this time, already using the calibrated hkl phase model for montmorillonite. It is important to emphasize that the removal of the initial region of the broad diffraction band is necessary along with the adjustment of the background using a linear interpolation between the given points. Mean and SD 74.9 ± 0.1 18.5 ± 0.1 SD: standard deviation. Figure 10. Residual plot of quantification of seven standard binary mixtures using the calibrated hkl file generated for Mg-montmorillonite.

CONCLUSIONS
The combined Rietveld-Le Bail-Internal Standard method of Paz et al. (2018) -for the purpose of mineralogical quantification -was successfully reproduced in this work. The calibrated hkl phase model developed can be used for the mineralogical quantification of bentonites, due to the reproducibility and greater precision in the results. The exclusion of a poorly resolved montmorillonite peak helped to solve the lack of fit in the complex region of the diffraction band, handling turbostratic disorder of montmorillonite, without hampering the generation of the hkl model. These results make the combined method a good choice for the bentonite industry, being convenient and fast.
According to the chemical analysis, the Formosa Mg-montmorillonite was classified as highly charged beidellitic montmorillonite. Knowing the degree of purity of the sample was a fundamental requirement for the development of the calibrated hkl model. However, the chemical composition did not influence the validation of the calibrated hkl model using the combined method.

ACKNOWLEDGMENTS
This study was partly funded by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) -Finance Code 001, along with a scholarship to the first author. We also thank the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for its financial support: Edital Universal MCTI/CNPq No. 01/2016(Grant 402.427/2016 Table 11. Quantitative result of the binary mixtures and statistic error indices of the Rietveld refinement using the calibrated hkl phase model for Mg-montmorillonite.