Some arsenic compounds can show extraordinary polymorphism. Realgar (As4S4) is among several minerals with various crystal forms and is one of the most important sources of arsenic for pharmaceutical use. Currently, realgar is used as an arsenic source in many industries, such as weaponry, publishing, textiles, cosmetics, and health products. In this paper, we used and reported new methods for the purification, nanonization, and structural morphological investigations of As4S4 by using planetary ball mills process for nanonization of the compound. The product was characterized using X-ray powder diffraction analysis, Fourier transform infrared spectrometry spectra, and field emission scanning electron microscope (FESEM) imaging. We investigated the morphological properties of FESEM-imaged realgar nanoparticles by an image-processing technique that calculates fractal dimensions using values on a computer with MATLAB software. We applied the Statistical Package for the Social Sciences software for statistics data extracted from the FESEM image and obtained the statistics results of the fractal dimension and histogram plot for the FESEM image.

1. Introduction

The anticancer action of arsenic compounds was reportedly discovered in Boston Hospital in 1878, where the influence of Fowler’s solution (1% As2O3 in K2CO3) in reducing the number of white blood cells in leukemia patients was clinically defined [1]. Today, arsenic oxide As2O3 (under the brand name Trisenox) is widely applied for treating acute forms of leukemia [24]; however, its solubility is low and a certain type of pretreatment is required to enhance the solubility and/or to prepare nanosized particles with superior bioavailability and improved curative capacity [5].

Realgar is among the numerous minerals used for thousands of years in China as a medicine for treating anthrax, abdominal pains, children’s convulsions, burning, and drying warts and treating abscesses and insect bites [6]; some reports indicate its use in treating CML and APL cancers as well [7]. Realgar has historically been used as a pigment. Realgar has a cage-like structure, with four arsenic atoms and four sulfur atoms, and can be found in veins and sublimation deposits in volcanic rocks, often found with gold, silver, and lead. The related mineral orpiment (As2S3) is also found in association with this mineral. Realgar, orpiment, and arsenopyrite are the principal ore minerals which provide for global arsenic demand [8]. In general, arsenic compounds are considered both toxic and pharmaceutical, depending on their chemical forms [9]. It is interesting to know that two factors, including molecular structure and dosage, determine either the pharmaceutical or toxic characteristics of arsenic [10].

In this paper, α-realgar was investigated as a polymorph of realgar. The difference between α and β As4S4 structures lies in their different molecular arrangement. The preparation of realgar was carried out by both traditional and modern methods, including acid extraction, membrane analysis, calcinations, various forms of milling (dry-, cryo-, and nanomilling), and chemical synthesis of particles by a stoichiometric As-S reaction [11].

Mechanochemistry can also be used to synthesize realgar because it can produce nanosized particles with improved solubility by a simple solid-state method [1216]. As for mechanochemistry’s application in pharmaceutical sciences, features such as bulk disordering, polymorphous transformations, and preparation of nanocrystalline powders via milling techniques have brought new impacts to pharmaceutical science [1721]; here, the improved rate of dissolution and bioavailability plays a crucial role in the application of new drugs [2225].

Image processing of nanostructures involves several preparation steps, separation, and postprocessing that finally lead to the extraction of quantitative data. Recently, computer-based image-processing methods have advanced rapidly, allowing quantification of complex colors, shapes, texture properties, and sizes. Image-processing methods are regularly used in tandem with mechanical and instrumental devices to change human manipulation in the display of an assumed process [26, 27].

Fractals are irregular forms with equal irregularities across all scales. The fractal dimension is a parameter for investigating the degree of data complexity and, in contrast to the naturally numbered Euclidian dimension, can be a real number [28]. The fractal geometry within the data may be calculated; exploration of the fractal dimension can lead to better investigation and prediction of the future behavior of data [29].

Various methods are used for estimating the fractal dimension. The box-counting method [30] is one of the simplest methods, which focuses on counting boxes in a straight line plot, where boxes with different lengths from and the number of boxes (tiles) are counted. If presents the obtained length in the scale, the fractal dimension can be calculated using the slope of the line in the log-log domain.

Our study presents descriptive data for the fractal dimension of As4S4 nanoparticles. This method is basically a function in MATLAB (MathWorks, Natick, Massachusetts, USA), which is a powerful software package for drawing and analyzing data, programming, and performing engineering and research calculations [3134]. Thus, we analyzed data from MATLAB with Statistical Package for the Social Sciences (SPSS; IBM, Armonk, New York, USA) software for discussing statistical aspects and give comparable results of our study.

2. Experimental

2.1. Materials and Methods

The starting materials were obtained from Merck (Berlin, Germany) and were used without further purification. Nanomaterials underwent spectrometry using a Fourier transform infrared (FT-IR) Bruker-Tensor 27 in the 400–4000 cm−1 range. The surface morphology of product was characterized by a field emission scanning electron microscope (FESEM) (Hitachi S 4160, Japan) with an accelerating voltage of 20 kV. X-ray powder diffraction (XRPD) measurements were performed using a Philips X’pert diffractometer (PANalytical, Almelo, the Netherlands) with monochromatized CuKα radiation (0.154 nm). The sizes of the selected samples were estimated using the Scherer method. For identification, the FESEM samples were coated with gold.

2.2. Synthesis of As4S4 Nanoparticles

First, 15 mg of powdered realgar (As4S4) was added to 100 mL distilled water and the solution was boiled for 3 hours to remove soluble impurities. After three stages of filtration, and reduction of the pH to 5 by the addition of a few drops of 2 M HCl, the remaining solid realgar was put in an oven for 10 hours at 80 degrees centigrade to dry. The dried compound was then transferred to one of the polymer jars of the planetary ball mills that contained 15 balls with different sizes for 5 hours of nanotreatment at room temperature (ball mill process).

2.3. Computational Method

Using MATLAB, 20 sections of the FESEM images of As4S4 were selected and their fractal dimensions were calculated. The obtained fractal dimensions were fed into SPSS software to analyze the normal distribution, correlations, standard deviation, mean, cumulative frequency, and variance.

Considering the resulting numbers, we drew a logarithmic diagram in which the vertical and horizontal axes were and , respectively, and the coefficient of the box size was 1.3. We used a program written in MATLAB to calculate the fractal dimension. There are several methods for calculating the fractal definition, including (1) the self-similarly definition [35], (2) the mass fractal dimension [36], (3) the spare-sphere fractal dimension [37], (4) the Euclid dimension [38], (5) the grid style [39], (6) an analytical approach for fractal diminishing [40], and (7) the box-counting fractal dimension [30].

We used the self-similar method for determining the fractal dimension from the FESEM image, which employed the box-counting method for the fractal dimension. Using the box-counting dimension method, we calculated the fractal’s dimension theoretically by applying the following formula:where is the fractal dimension and is the number of squares including a part of the considered fractal dimensions.

For the box-counting fractal dimension, we assumed that is a nonempty and bounded subset of and we supposed that is the fewest number of collections with a maximum dimension of that can cover . The box-counting dimensions under and above (): if they are equal we say they are equal to the box-counting dimension or the box dimension; then we define:

The circle of cubes in coordinates from indicates cubes in the following form:

3. Results and Discussions

3.1. Interpretation of the FT-IR Spectrum of As4S4 Nanoparticles

The FT-IR spectrum was taken within the 400–4000 cm−1 range to evaluate the purity of the synthesized realgar nanoparticles. Realgar is an inorganic compound and so lacks the spectral complexity of many organic materials and is without obvious factor groups. Upon adjusting the pH, some acid and alkali were added to the solution, which caused absorption peaks; however, there were some impurities, and the widened peaks at 3500 cm−1 clearly show the included water (Figure 1). These spectra show the strong peaks at 2924 and 1152 cm−1 for the C-H and C-O vibrations that arose when acids and bases were added for pH adjustment [41].

The main spectrum at 481 cm−1 shows the As-S vibration, which confirms the realgar compound; many special vibrations for realgar are below 400 cm−1 and these vibrations do not appear in the spectra [9, 42].

3.2. Variance X-Ray Powder Diffraction Analysis of As4S4 Realgar Nanoparticles

To obtain the crystallite size distribution of the particles by the half-height method, an FWHM XRPD pattern was used (Figure 2). As observed, the XRPD pattern presents maximum intensity peaks, indicating the small sizes of crystals in the samples synthesized. To obtain the particle size and the peak width by the half-height of FWHM, we applied the Debye-Scherrer equation to get the approximate crystallite sizes of the particles:where is particle size by nm; is constant equivalent to 0.99; is wavelength of CuKα (equivalent to 0.154 nm); is peak width by half height; and is variance angle by degrees:The particle size was calculated by the Debye-Scherrer Equation and revealed that the particle size of the maximum FWHM (0.0984) was around 91 nm.

3.3. Morphological Properties

The image of the surface and morphological structure of the sample indicates that the realgar nanopowder obtained was condensed and highly porous and was comprised of 20 nm spherical particles; larger particles observed in the surface of the image were formed as a result of the aggregation of the smaller particles, which manifested as irregular spherical shapes. FESEM images of the realgar nanoparticles are shown in Figure 3.

3.4. Fractal Properties
3.4.1. General

For analyzing fractals in MATLAB, FESEM images are required. In the present research, FESEM images were used for imaging particle surfaces to determine the fractal dimension. Having used MATLAB for determining the fractal dimension of realgar, SPSS software was used to analyze the MATLAB data. SPSS can determine the frequency of groups in one variable, calculate a simple mean for data, prepare data for testing the relationships between variables, calculate correlations, perform multivariate regressions, and analyze diagnostic functions, linear logarithms, and other multivariate logarithms. In the present research, other characteristics, such as statistical analyses, normality tests, correlation of variables, cumulative frequency, standard deviation, variance, mean, and harmonic mean, were extracted in SPSS.

3.4.2. Analysis of MATLAB and SPSS Data

The descriptive statistics of the fractal dimension are listed in Table 1. This table shows the fractal dimensions of 20 random, separated pieces of the FESEM image with frequency and the valid and cumulative percent of each piece. When dealing with statistical information, we come through numbers, which are sometimes repeated; after the similar value appears more than once, it is recognized as the frequency of a number. The total frequency of all modules less than the higher class boundary of a quantified class is named the cumulative frequency of that class. Cumulative frequency is found by addition of a class interval to the frequency and the frequencies of the former intervals up to that class interval. This is another way of expressing frequency distribution. Cumulative frequency was used to determine the number of explanations that were below or above a specific value. The highest advantage of using cumulative percentage is that it offers an easier method to compare different sets of data.

In statistics and probability theory, skewness reflects the degree of asymmetry of the probability distribution. Skewness is a measure of the presence or absence of distribution function symmetry: for a perfectly symmetrical distribution, the skewness is zero; for an asymmetric distribution with stretching toward higher quantities, the skewness is positive; and for an asymmetric distribution with stretching toward smaller quantities, the skewness is negative. Considering Table 2, we observed the skewness to be −2.178, which indicates an asymmetric distribution; because the obtained figure is negative, the asymmetry is skewed towards smaller values.

A grayscale digital image is collected from separate points of gray tones, or brightness, before incessantly variable tones. A normal image is separated into a number of specific points of brightness, and each of those points is defined via a digital data value. A pixel is the greatest elementary element of a digital image; each brightness point is a pixel of the digital image. Essentially, the image histogram shows the distribution of the pixel intensities in the image and those continuously used as a reference. Interactive thresholding can be very effective and afford fast, precise information [43, 44]. The histogram plot extracted with SPSS software is given in Figure 4.

3.4.3. Model Description

In the obtained sections of As4S4 nanoparticles, the minimum and maximum fractal dimensions were 0.84 and 1.78, respectively; the mean fractal dimension was 1.53 (in statistics, mean signifies central tendency). The normal and detrended normal P-P and Q-Q correlations are plotted in Figure 5. It was observed that the distribution of the fractal dimensions of the As4S4 nanoparticles was aligned in a straight line, and the hypothesis of a normal distribution of errors was true. The drawing normal and detrended normal P-P and Q-Q plots for residuals are fitted here. If the fitted model is appropriate, the figure must be symmetric, and the points must be scattered consistently; this indicates a fixed variance of errors, which is optimal in nature. Summarizing all these data can help us determine the fractal characteristics (Figure 5).

Correlations are associated with value and direction, like a vector: its value indicates the degree of relationship between two variables, and its direction indicates the direction of behavior (either same direction or the opposite) of the variables under investigation. It is a statistic for measuring the power of degree of the linear relationship between several variables and can range from +1 to −1. Values close to +1 and −1 show strong correlation (either positive or negative); however, the more they approach zero, the less powerful the correlation will be. In SPSS, one of the correlation coefficients can be used to determine the degree and type of relationships, based on the type and nature of variables. In the present paper, we used the Pearson correlation coefficient, which is applicable to quantitative variables.

Where continuous variables are considered in normal distributions, this equation is strongly affected by sample size. Therefore, we investigated more points to ensure that we applied a powerful test. The Pearson correlation coefficient is based on the covariance and deviations of the variables, where their estimations can be used for calculating the Pearson correlation coefficient. Also, as Table 3 indicates, the correlation of the fractal dimension of nanoparticles of As4S4 is 0.994; this indicates a strong positive correlation between the variables. The obtained coefficient indicates a strong relationship between the different sections of the image, which all confirm that this mineral is highly fractal. In short, an extremely complex order can be observed in an irregular and disordered appearance.

4. Conclusions

New methods for the purification, nanonization, and structural and morphological investigations of As4S4, using a planetary ball mills for nanonization of the realgar compound, were reported. The realgar nanoparticles experimentally characterized by using XRPD analysis, FT-IR spectra and FESEM image analysis, and the obtained FESEM images were used for investigate the phase composition, morphological, and fractal dimension calculations by MATLAB and SPSS software. The image of the surface and morphological structure of the sample indicates that the obtained nanopowder was condensed. We conclude that the calculation of the fractal dimensions of the FESEM images of 20 randomly selected particles of the nanomaterial and calculation of the Max, Min, and many other data can be used for the analysis of nanoparticles; therefore, homogeneity and uniformity were considered. The correlation of the fractal dimension of the nanoparticles of As4S4  gave a value of 0.994, which indicates a strong positive correlation between the variables.

The realgar nanoparticles are planned for use in patients with gastric cancer, since it can alleviate pain for a while; it will not be used in trials with healthy individuals because of the highly toxic nature of inorganic arsenic compounds. Finally, investigating the use of this compound for therapeutic purposes will be considered by oncologists.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors gratefully acknowledge the financial support from the Research Council of Imam Khomeini International University. Daniel Glossman-Mitnik is a researcher of CONACYT and CIMAV and acknowledges partial support from both institutions.