As nanotechnology advances, nanosafety and the risks of toxicity and ecotoxicity associated with nanomaterials have become a matter of great concern due to their related contributions to the human health and the environment; since the nanomaterials could appear in air, water, soils, plants and consequently in human and animal bodies [1, 2].
QDs are inorganic semiconductor nanocrystals composed of a heavy-metal core and usually a shell, which have grasped a significant attention in recent years . Their unique optical and electronic properties are captivating, leading to a variety of research and commercial applications including bioimaging, solar cells, LEDs, diode lasers, and transistors [4-7]. Until recently, cadmium-based QDs have been assuredly the most progressed, since the synthesis is straight-forward and their band gaps lie in the visible region of the spectrum, allowing for simple characterization. The inherent (eco)toxicity of these QDs has hindered their applicability, motivating researches into alternative, less toxic QDs [8-12]. I−III−VI CuInS2 QDs are environmentally friendly and biocompatible which have emerged as particularly exciting materials for the synthesis of a new class of QDs, not only they contain no heavy metal ions, but also due to the unique structural and electronic properties that arise from the composition and structure of ternary semiconductor compounds in general [13, 14].
Various approaches have been reported for CuInS2 syntheses such as solvothermal synthesis, single-source precursor routes and hot injection techniques which are based on organic solvents so their applications have been faced with several problems [15-18]. Direct hydrothermal synthesis of CuInS2 QDs offers advantages such as lower reaction temperature with comparable photoluminescence quantum yield (PLQY), not using toxic and expensive organometallic reagents, does not need any surface functionalization during synthesis, no longer producing harmful bi-products, comparatively smaller particle size and tunable nanoparticles᾽ size and morphology by controlling the precursors᾽ concentration, temperature and time of reaction [19, 20].
There are few methods for direct hydrothermal synthesis of CuInS2 QDs, while to tackle this problem, developing a proper synthesis method seems to be a necessity. According to an extensive survey on the literature, a very efficient facile direct synthesis method producing small size CuInS2 QDs has been discovered in a work that Liu et al. do in order to prepare water-soluble high-quality ternary CuInS2 QDs with mercaptopropionic acid (MPA) as the stabilizer by a novel hydrothermal synthesis route[8, 21, 22]. Many factors affect the CuInS2 QDs synthesis such as reaction temperature, reaction time, precursors concentration ratio, pH and ligand/precursor concentration ratio. These parameters change the PL intensity and the size of CuInS2 QDs. So it is necessary to optimize them for producing appropriate nanoparticles [19, 23, 24].
In order to achieve the optimum synthesis condition, previously one parameter has been changed while the others have been kept unchanged during the experiments. By this way, the influence of each parameter could be understood solely [25-30]. These traditional optimization methods have some disadvantageous. They require a considerable amount of work and time due to the high number of experiments. Also, the influence of parameters interaction couldn’t be investigated . In order to optimize the experimental synthesis conditions, multivariate statistical techniques can be employed. Among these optimization techniques, response surface methodology is the most popular tool in process optimization, consequently, the most applied . Response surface methodology is a mathematical and statistical method that widely used for experiential modeling in order to evaluate the influence of independent variables on dependent ones [33, 34]. Optimal designs initiate with a pseudo-random set of model points that enable the model fitting. The first choice can usually be upgraded by replacing a subset of the points with improved selections . The optimal design uses several criteria to decide which replacements are suitable for the studied system, which the D-optimal algorithm is the most used criteria, which chooses runs that minimizes the determinant of the variance-covariance matrix .
To the best knowledge of authors, no study has been done on the optimization of the synthesis condition of CuInS2 QDs using the response surface methodology approach. In a study, Pian Wuet al. optimized the formulation variables and process of acid-modified ZnSe/ZnS core/shell QDs using the response surface methodology . In another study, Box–Behnken design (BBD) and response surface methodology were adopted to optimize the synthesis condition for ZnSe/ZnS core/shell QDs via microwave irradiation . So, employing of RSM technique for CuInS2QDs synthesis condition can be considered as a new aspect of this work. At the present research, RSM coupled with D-optimal design was used to study the relationship between the independent (reaction time, reaction temperature, pH, precursors concentration ratio and ligand/precursor concentration ratio) and dependent (PL intensity and the desirable production output based on maximum PL intensity) variables in the optimized condition.
MATERIALS AND METHODS
All reagents were analytical grade and used directly without any purification. Copper(II) chloride dehydrate (CuCl2.2H2O), Indium(III) chloride tetrahydrate (InCl3.4H2O), Sulfourea (CS (NH2)2), Mercaptopropionic acid (MPA), and Sodium hydroxide (NaOH) were purchased from Sigma-Aldrich Corporation (St. Louis, MO, USA).
Synthesis of MPA-capped CuInS2 QDs
The hydrothermal synthesis of hydrophilic MPA-capped CuInS2 QDs was achieved by adopting a literature method in which the surfactant MPA acts not only as both stabilizing ligand and source of Sulphur for the nanoparticles but also as the reaction solvent [8, 22, 24]. The underlying principle behind this method is that excess of thiol promotes complete surface ligand coverage and therefore, good colloidal stability.
In this work, Indium(III) chloride and Copper(II) chloride dehydrate were used, because of their demonstrated reactivity and solubility. In a typical experiment, CuCl2.2H2O (0.15 mmol) and InCl3.4H2O (0.15 mmol) were dissolved in distilled water (10.5 ml). Then MPA (1.8 mmol) was injected into the solution, producing opaque yellow granules immediately. The pH value of the mixture solution was adjusted to 11.3 by the drop-wise addition of 2 mol/L NaOH solution with gentle magnetic stirring. During this process, the solution changed from turbid to clear pink (Fig.1, step 3). After stirring for 10 min, CS (NH2)2 (0.30 mmol) was dissolved in the solution. The Cu–to–In–to–S and Cu–to–MPA precursor ratios were 1:1:2 and 1:12, respectively. All the above mentioned experimental procedures were performed at room temperature, and then the solution was transferred into a Teflon-lined stainless steel autoclave with a volume of 15 ml. The autoclave was maintained at 150 ̊C for 21 h after which cooled down to room temperature, quenching the reaction by a hydrocooling process. A purification method based on solvent extraction was used to separate unreacted precursors and reaction by-products from the as-synthesized hydrophilic QDs suspended in a nonorganic solvent. Ethanol was added to the solution to obtain CuInS2 QDs precipitate, and the process was repeated three times. The unreacted residues were removed by the cycled washing. The CuInS2 QDs was dried at 60 ̊C for 4–6 h (Fig. 1, step 5). The obtained powder was used for further measurements [8, 22, 24].
Field emission Scanning electron microscopy (FESEM) experiments were performed on a TESCAN-Mira 3-XMU (Kohoutovice, Czech Republic) - operating at 3-30 KV voltage. FTIR spectra were recorded with a Beijing Rayleigh Analytical Instrument Corporation (Beijing, China). BRAIC-WQF-510 FTIR Spectrometer equipped with a DTGS detector (16 scans). UV-VIS absorption spectra of CIS QDs were obtained using UV-VIS spectrophotometer, LAMBDA™ 25, Wavelength range 190-1100 nm, PerkinElmer Co. (Waltham, MA., USA). The Photoluminescence spectra were measured by spectrofluorometer, Avaspec 2048 TEC, Avantes Co. (Apeldoorn, The Netherlands).
Design of experiment using response surface methodology and D-optimal approach
In order to perform a regression model, The Design Expert Software, version 10, based on RSM and D-optimal design was used. According to the primary studies, MPA/Cu (A), time (B), temperature (C), pH (D), and In/Cu (E) were selected as independent variables (Table 1) and PL intensity was chosen as the response to this research with the purpose of achieving the highest PL intensity which is desired in QDs applications.
The complete experimental runs and the corresponding responses are illustrated in Table 2.
The response behaviors can be defined using a quadratic polynomial model  by the following equation:
Where y is the response (dependent variable),
xi and xj are independent variables, ε is the residual, β0 is the constant coefficient, βi, βii and βij are the coefficients for the linear, quadratic, and the interaction terms, respectively .
RESULTS AND DISCUSSION
Response surface methodology and D-optimal approach
The PL intensity response can be described by a reduced quadratic equation based on coded factors as follow:
PL intensity= -0.020+0.16 A+0.098 B+0.73 C
+0.39 D+0.12 E-0.25 A2+0.55 C2+0.16 E2(2)
The interaction and some quadratic coefficients of PL intensity model were insignificant and excluded from the model. The equation demonstrated that the PL intensity model was sensitive to all factors. Model assessment was done by the coefficient of determination (R-squared). According to (Fig. 2), the R-squared value of 0.9201 for PL intensity and model shows that 92.01% of the variability is matched by the experimental data (Fig. 2), which a reasonable agreement between predicted and actual values is understood.
The analysis of variance (ANOVA) was applied to verify the model accuracy [33, 40]. The ANOVA for the predicted model is listed in Table 3.
The model significance can be described using the p-value and F value, which the smaller p-value (<0.1) and larger F value show the more significant model coefficient [34, 41]. As shown in Table 3, the F-value of 20.14 for PL intensity confirms that the model is highly significant. In addition, the model adequacy was further verified by the p-value less than 0.0001.
The influence of replicate points and errors dependence of model can be surveyed using the lack of fit test [33, 40]. As insignificant lack of fit is necessary, the model fitting was confirmed by the p-values of 0.8097 for PL intensity .
Moreover, the adequacy of the model is also should be studied using plots of normal probability and residuals vs. predicted [33, 38, 40] (Fig. 3). As shown in (Fig. 3), as expected the straight line of residual points on the normal probability plot for PL intensity (Fig. 3a) is revealed, confirming the normal distribution of errors [38, 40]. In addition, the assumption of constant variance is tested by residual versus the predicted responses. The plot should be a random scatter, which observed for PL intensity (Fig. 3b) .
The perturbation plots of the PL intensity as a response is shown in (Fig. 4). This plot helps to compare the impacts of all the factors at a particular point in the design space, which the response is plotted by changing only one factor over its range at a constant value of other factors . As shown in (Fig. 4) the most effective parameters on PL intensity was temperature, while the other factors were less sensitive to PL intensity. As the reaction temperature increases the PL intensity increases due to better ligand passivation, while it has been decreased severely at higher temperatures. This phenomenon happens because of the pyrolysis of MPA and poor ligand passivation.  Increase of nanoparticles size resulting in agglomeration is another reason for the reduction in PL intensity [44, 28].
The 3D plots at (Fig. 5) show the simultaneous effect of temperature (as the most important parameter) and the other factors on PL intensity, which increased significantly by temperature enhancement. A linear and nonlinear growing trend in PL intensity was observed by increasing pH (Fig. 5c) and MPA/Cu (Fig. 5a) while the time and In/Cu did not show a considerable effect on PL intensity (Fig. 5b and d). MPA/ Cu and pH can affect the decomposition of Cu-In-MPA complex. In the low concentration of MPA as the temperature and pH increases, PL intensity of nanoparticles decreases due to the formation of Cu(OH)2 precipitates. Excess MPA/Cu distort the surface to form new non-radiative defects which quench the PL intensity of as-prepared QDs . Shortage of surface ligands causes instability and agglomeration of QDs resulted in diminishing PL intensity.
Model verification and optimization
As the maximum values of PL intensity were crucial for optimal condition, the multi-response optimization was applied using desirability function. Table 4 shows optimized new experimental test as optimum run without constrains. The result verified that the experimental data and predicted values were in good agreement and the relative error was around 5%, confirming that the model had enough accuracy to predict removal percentage for this system.
Preparation of CuInS2 QDs under the optimum condition
Structural characterization: The surface structural analysis of the as-synthesized CuInS2 QDs were undertaken using FESEM and FTIR. FESEM observation for CuInS2 nanoparticles produced under the optimum condition is shown in (Fig. 6). By analysis of the FESEM image, we determined the particle size of most of CuInS2 QDs obtained was approximately 7 nm.
Because of their small dimensions, QDs have an inherently large surface to volume ratio, therefore, surface properties play a crucial role in their conjugation and also photoluminescence properties. As such, it is important to investigate the functionalization of the surface, determining the extent of surface coverage, the nature of the capping ligands, including their tendency to associate with surface atoms and surface defects [44-47]. To further characterize the as-synthesized CuInS2 QDs, FTIR has been carried out (Fig. 7).
The FTIR spectra pointed that most functional groups of the as-synthesized QDs could be clearly found through the characteristic peaks of O–H (3448 cm-1 stretching vibration), –COOH (2370 cm-1 asymmetric stretching vibration), C=O (1705 cm-1 stretching vibration), C-O(1480 cm-1 stretching vibration). The absence of the characteristic peak of S– H between 2550 to 2680 cm-1 indicated that the final CuInS2 nanoparticles contain MPA on their surface, which might be caused by the covalent bonds between thiols and metal atoms [19, 24] of the ternary QDs.
Optical Characterization: The emissive properties of CuInS2 QDs were explored with fluorescence spectroscopy; a typical PL spectrum for CuInS2 QDs under the optimum condition is shown in (Fig. 8a). It can be seen that the nanoparticles exhibit excellent fluorescence emission spectra with the emission peak around 675 nm. As obvious red shift of emission for as-synthesized QDs at the optimum condition compared to the synthesis in other conditions is clear. The redshift of emission peak demonstrates that the nanoparticles are in the near infrared region which is more suitable for biological imaging and detection . The emission band is narrow and symmetrical. A typical absorption spectrum is also shown in (Fig. 8b); The UV-VIS absorption spectra of the as-prepared CuInS2 QDs have been measured at room temperature.
At the present research, CuInS2 quantum dots have been successfully synthesized by hydrothermal method. The influence of independent variables including MPA/Cu, time, temperature, pH, and In/Cu on the PL intensity as a response was studied using RSM and D-optimal design and the optimum condition for producing nanoparticles with high PL intensity was achieved. UV-VIS, FESEM, and FTIR analysis on the produced QDs in the optimum condition showed the nanoparticles have a near-infrared narrow PL spectra with high PL intensity. A reduced quadratic equation has been proposed to describe PL intensity. The statistical analysis and reduced quadratic models showed that the temperature variables were the most effective parameters of PL intensity. A good agreement between the predicted results and the experimental data was observed in model analysis and optimization step.
CONFLICT OF INTEREST
The authors declare that there are no conflicts of interest regarding the publication of this manuscript.