Size and Shape-Dependent Solubility of CuO Nanostructures.

In our theoretical study, the enhanced solubility of CuO nanoparticles in water saturated by air is predicted based on a simple thermodynamic model. CuO is considered in the form of nanoparticles with various shapes. The interfacial energy of a solid CuO/dilute aqueous solution interface was assessed by applying the average CuO surface energy and contact angle of a sessile drop of water. The equilibrium CuO solubility was calculated using Gibbs energy minimization technique. For the smallest spherical nanoparticles considered in this work (r = 2 nm), the solubility is significantly higher than the solubility of bulk material. In the case of cylindrical nanoparticles, the solubility increase is even more considerable. The CuO spherical nanoparticles solubility was also calculated using the Ostwald-Freundlich equation which is known to overestimate the solubility as discussed in this contribution.

1. Introduction

In recent years, an increasing interest in the research of transition metal oxides nanostructures emerged due to their unique properties and behavior differing from the bulk counterparts [1,2,3,4]. In this research area, it is very challenging to measure exactly particle sizes or particle-size distribution curves, as well as to experimentally determine the solubility of such nanoparticles [5,6]. Such results can then be compared with theoretical calculations based on thermodynamic modeling revealing the equilibrium background of these nanosizing effects.

CuO nanostructures (CuO-ns) are potentially toxic substances to many different organisms, such as plants, algae, bacteria, fish, and mammalian cells [7,8,9]. The mechanism of their toxicity is not definitely explained but strong evidence has been put forward that the CuO dissolution in an aqueous media plays a dominant role [10,11]. On the other hand, some observations indicate that nanoparticles themselves might be a source of toxicity due to enhanced formation of reactive oxygen species (ROS) [12,13,14]. It is well known that the solubility is size and shape-dependent [15,16], which has been observed for CuO [17,18,19] as well as for other oxide nanoparticles, e.g., ZnO [9,20] and CeO2 [21].

Enhanced solubility of nanoparticles is one of the numerous consequences of reducing the size of matter. Such behavior of small (nano-)particles has been predicted more than one century ago. Following Gibbs treatment of interfaces and their influence on equilibrium between solid particles and a liquid solution, the Ostwald–Freundlich equation could be derived in the form: where Vm is molar volume, T is absolute temperature, R is gas constant, c and c,∞ stand for equilibrium solubility of a substance i in the form of spherical particles of radius r and bulk material (r → ∞), respectively, and γsl represents the interfacial energy between solid particle and solution. Equation (1) holds for the ideal solution of i in a solvent and is based on the assumption that the interface energy γsl is not size-dependent. It was shown in our previous work [22], that this equation gives reliable results for spherical nanoparticles with a diameter greater than approx. 20 nm. For smaller nanoparticles, it predicts more enhanced solubility compared to the more advanced approach based on complex equilibrium calculation. This is due to the fact that the Ostwald–Freundlich equation—Equation (1), predicts the enhanced solubility as a consequence of the lower nanoparticles stability in comparison with bulk materials, and no response to changes in solution composition is considered. Moreover, the influence of the nanoparticles shape is also not considered in Equation (1). Based on a study published by Mistra et al. [19], CuO solubility is highly influenced by the shape of nanostructured CuO (spheres, rods, etc.). The size-dependence has also been observed in the case of bactericidal activity of CuO nanoparticles [23].

The aim of this work is to predict the size and shape dependence of CuO solubility in the form of spherical nanoparticles, nanorods, and nanoplates. The equilibrium solubility (composition of saturated solution in equilibrium with CuO) was calculated using Gibbs energy minimization technique considering thirteen Cu containing aqueous species and six solid Cu based compounds. The essential parameters for the calculations, namely interfacial energy between solid CuO and dilute aqueous solution, were estimated on the basis of the literature values of CuO surface energy and the contact angle of a water drop on the CuO surface.

2. Thermodynamic Description of CuO Dissolution

2.1. Dissolution of Bulk CuO

The solubility of bulk CuO in pure water has been investigated many times (for review see example [24]), and the equilibrium total Cu content (molality) in solution in the range 10−7–10−6 mol kg−1 was obtained (≈8 × 10−6–8 × 10−5 g CuO in one liter of solution) [25,26]. The solubility is significantly dependent on the pH of the solution achieving a minimum at around pH ≈ 10 and increases with both acidity and alkalinity increase of the solution. Cu2+(aq) is the dominant aqueous copper species in a neutral and acidic solution, while a variety of hydrated Cu(OH)(2− (aq) species (products of Cu2+(aq) hydrolysis) are present in alkaline solutions. The solubility is also affected by the presence of some other complexing species, for example, Cl−, CO32−, or PO43−. Due to the carbonate ion or dissolved carbon dioxide, solid alkaline carbonates, namely malachite Cu2(OH)2CO3 or azurite Cu3(OH)2(CO3)2, can precipitate from aqueous solutions [26,27]. Moreover, some complex aqueous species can be formed which changes the total copper content as well as the copper speciation in the solution.

To calculate the equilibrium solubility of CuO in pure water and aqueous solutions, the non-stoichiometric method based on minimization of the total Gibbs energy of the system on the set of points satisfying the material balance conditions was used. The calculation algorithm and the computer program CHEMEQ were described elsewhere [28]. The liquid phase was considered as a dilute solution of 19 species in water (see Table 1). Chemical potentials of the solvent and the solutes were expressed with respect to Raoultian (pure substance) and Henrian (hypothetical ideal solution of the relevant solute in water at unit molality) standard states, respectively. To express the non-ideal behavior of the dilute aqueous solution, the Davies equation for solutes activity coefficients γ was used [29]: where z is the ion charge and the Debye–Hückel constant 0.509 applied for T = 298.15 K and the ionic strength on the molal basis Im is defined as: where m is molality.

Table 1

Substances considered in equilibrium calculation s and input thermodynamic data [30,31] (standard Gibbs energies of formation at T = 298.15 K and po = 100 kPa).

Substance	∆fGo (kJ/mol)	Substance	∆fGo (kJ/mol)
H2O(l)	−237.129	Cu2(OH)22+(aq)	−285.1
H+(aq)	0	Cu3(OH)42+(aq)	−633.0
OH−(aq)	−157.244	CuCO3(aq)	−501.5
O2(aq)	−16.4	Cu(CO3)22−(aq)	−1048.98
CO2(aq)	−385.98	CuHCO3+(aq)	−532.08
CO32−(aq)	−527.81	N2(g)	0
HCO3−(aq)	−586.77	O2(g)	0
Cu+(aq)	48.87	CO2(g)	−394.359
CuOH(aq)	−122.32	Cu2O(s)	−147.90
Cu(OH)2−(aq)	−333.05	CuO(s)	−128.29
Cu2+ (aq)	65.04	Cu(OH)2(s)	−359.92
CuOH+(aq)	−126.66	CuCO3(s)	−528.20
Cu(OH)2(aq)	−316.54	Cu2(OH)2CO3(s)	−902.35
Cu(OH)3−(aq)	−493.98	Cu3(OH)2(CO3)2(s)	−1431.43
Cu(OH)42−(aq)	−657.48

The gas-phase containing N2, O2, and CO2 were considered to be in equilibrium with the aqueous solution. The oxygen partial pressure p(O2)/po was fixed at 0.21 (standard pressure po = 100 kPa) and partial carbon dioxide partial pressure p(CO2)/po at 4 × 10−4. Apart from CuO, five other solid substances, namely Cu2O, Cu(OH)2, CuCO3, Cu2(OH)2CO3, and Cu3(OH)2(CO3)2 were included in the calculations.

The values of standard Gibbs energies of formation (ΔfGo) were used as input thermodynamic data. The values for Cu aqueous species as well as for solid Cu compounds were adopted from the recent assessment by Puigdomenech and Taxén [30], while all other aqueous species, liquid water, and gaseous O2 and CO2 were adopted from the NBS tables [31].

2.2. Dissolution of CuO Nanostructures

Enhanced dissolution of nanostructured solids is due to the interface effect, which results in lower thermodynamic stability of materials with large specific interface area in comparison with bulk materials. As the first approximation, the enhanced solubility of spherical nanoparticles can be simply calculated using the Ostwald–Freundlich equation—Equation (1). As discussed later, this equation predicts a higher solubility for small nanoparticles (diameter under approx. 20 nm) than a more advanced approach based on complex equilibrium calculation taking into account the mass exchange between the solid and solution phase.

A simple procedure for equilibrium calculation in nanosystems was presented recently [32] and this approach was used for the demonstration of size/shape-dependent solubility of ZnO nanoparticles [22]. The interface effect can be simply treated by redefinition of the chemical potential of a solid substance in the form of nanoparticles surrounded by a liquid (fluid in general) phase. It is applicable for nanoparticles of general shape whose chemical potential μnp can be expressed as: where the so-called differential shape factor α′ is defined as [32]: and reqv = (3Vnp/4π)1/3 is the radius of a spherical nanoparticle with the same volume as a non-spherical one. The interfacial energy γsl is dependent on the composition of the surrounding solution but in the case of dilute solution, a constant value referring to the solid/solvent interface can be used.

The above-described calculation of bulk CuO solubility in pure water was extended to CuO nanoparticles of various shapes (spherical as well as cylindrical with different aspect ratios). The chemical potential of CuO is given as: with α′ = 1 for spherical nanoparticles. In the case of cylindrical nanoparticles with a height h and radius ρ with a constant aspect ratio x = h/ρ, the differential shape factor α′ is given as [32]: and reqv = ρ(3x/4)1/3. The shape factor values α′ = 1.145 and 1.726 are obtained for x = 2 and 20, respectively. The value of Vm,CuO = 12.207 × 10−6 m3 mol−1 was used for solid CuO [33]. As for the interfacial energy, the average value γsl = 690 mJ m−2 was estimated as described in the next section.

Equilibrium calculations are more complicated in systems where two (or more) solid substances (phases) can coexist. In such a case, various geometrical models are available [34]—two (or more) single-component particles or one two (or more) component particles with core-shell or Janus geometry. Some other parameters such as interfacial energies on solid/solid interfaces are indispensable for such calculations.

3. Estimation of the Interfacial Energy CuO/Aqueous Solution

The interfacial energy γsl at the CuO/aqueous solution interface can be estimated using the Young equation and the values of surface energy γsg of solid CuO and contact angle θ of the water drop on the CuO surface (pure water can be used as an approximation of a highly diluted aqueous solution):

The surface energy of liquid water γlg = 72.0 mJ m−2 was used for calculation [35]. Surface energies for various crystallographic planes of CuO with monoclinic crystal lattice have been calculated using an ab-initio approach [36,37,38]. The individual values (see Table 2) were averaged according to the formula [39]: yielding the resulting surface energy 1106.7 mJ m−2. This value refers to vacuum surroundings, but for real conditions, one must consider surface energy lowering due to the adsorption of some species from the surrounding atmosphere on the solid surface (primarily water vapor during contact angle measurements of the water drop). Analyzing literature data on calorimetrically determined surface energies for hydrated and anhydrous surfaces for a number of oxides [40], one can estimate that this reduction is approx. 38% and so γsg for a hydrated CuO surface is ~690 mJ m−2.

Table 2

Calculated values of surface energies γ( for solid CuO [36,37,38].

(hkl)	γ(hkl) (mJ m−2)
	Ref. [36]	Ref. [37]	Ref. [38]
(111)	740	720	750
(−111)	860		890
(011)	930	910	940
(101)	1160		1170
(110)	1290	1180	1185
(010)	1370	1680	1485
(100)	2280	2240	1755
Average-Equation (9)	1094.6	1142.4	1083.2

As for the contact angle of a water drop on a copper oxide surface, a lot of measurements have been performed exclusively on thin nanostructured films and coatings [41,42,43,44]. The obtained values cover quite a large range from 10° to 165° according to the actual composition of those films (presence of Cu, Cu2O, or Cu(OH)2) and their condition (as prepared, aged, further treated, and so on). Moreover, contact angles show film thickness-dependence in some cases. In such a situation, we decided to use a neutral value of contact angle 90° (cosθ = 0) and thus γsl ≈ γsg = 690 mJ m−2.

4. Results and Discussion

The solubility of bulk CuO was calculated at T = 298.15 K and p = 101.325 kPa. The aqueous solution was considered as saturated with O2 and CO2 by contact with gaseous-phase with fixed partial pressures p(O2)/po = 0.21 and p(CO2)/po = 4 × 10−4. The results are summarized in Table 3.

Table 3

Calculated solubility of CuO in water at T = 298.15 K and p = 101.325 kPa (initial conditions: no(CuO) = 1 mol, no (H2O) = 1 kg (55.5084 mol), p(O2)/po = 0.21).

p(CO2)/po	m(Cu)tot(mol kg−1)	pH	Im(mol kg−1)	Dominant Aqueous Cu Apecies #
4 × 10−4	7.74 × 10−6	6.41	2.29 × 10−5	Cu2+(95.8), Cu(OH)+(2.6), CuCO3 (1.0)
0	1.10 × 10−7	7.37	3.21 × 10−7	Cu2+(77.3), Cu(OH)+(19.9), Cu(OH)2 (2.5)

# Number in parentheses means the percentage of the total Cu content in solution.

Under the above-mentioned conditions and at p(CO2)/po = 4 × 10−4 (average concentration of carbon dioxide in air), the calculated CuO solubility is m(Cu)tot = 7.74 × 10−6 mol kg−1 (0.616 mg CuO) at pH = 6.41 and no other solid Cu substances is stable. Malachite can be formed by Equation (10) in a carbonized aqueous solution

At 298.15 K, the equilibrium value of CO2(aq) molality is m(CO2(aq)) = 1.07 × 10−4 mol kg−1 which corresponds to CO2 partial pressure p(CO2)/po = 3.15 × 10−3. Hence, since the CO2 activity considered in our calculation is one order of magnitude lower, malachite is not formed at these conditions. The calculated CuO solubility in the limiting case p(CO2)/po → 0, namely m(Cu)tot = 1.10 × 10−7 mol kg−1 (8.75 × 10−3 mg CuO) at pH = 7.37, is in good agreement with the data from the literature [24,25,26].

Enhanced solubility of CuO was calculated first for spherical nanoparticles of radius r according to the Ostwald−Freundlich equation—Equation (1) which, after substituting the numerical values, assumes the form:

The size-dependence of solubility is shown in Figure 1. For the smallest particles considered in the present study (r = 2 nm), such calculation gives almost 30-times higher solubility than for bulk CuO.

Figure 1

Calculated solubility of CuO−np in water at T = 298.15 K and p = 101.325 kPa.

The results based on complex equilibrium calculations are summarized in Table 4 and plotted in Figure 1. It is obvious that there are serious discrepancies between the Ostwald−Freundlich prediction and our calculations. It is due to the fact that the calculation according to the Ostwald−Freundlich equation does not take into account the response of the solution to the reduced stability of dissolving nanoparticles. While for bulk CuO the resulting pH is 6.41 for the most soluble nanoparticles of radius r = 2 nm, pH = 6.91. Such a change in the final pH of the solution decreases the solubility of bulk CuO in approx. 6-times [24]. This effect is indeed not considered within the Ostwald−Freundilch equation.

Table 4

Calculated solubility of CuO−np in water at T = 298.15 K and p = 101.325 kPa (initial conditions: no(CuO) = 1 mol, no (H2O) = 1 kg (55.5084 mol), p(O2)/po = 0.21, p(CO2)/po = 4 × 10−4).

r or rekv (nm)	m(Cu)tot (mol kg−1)
	Sphere	Cylinder (x = 2)	Cylinder (x = 20)
2	2.959 × 10−5	3.779 × 10−5	1.433 × 10−4
3	1.811 × 10−5	2.072 × 10−5	3.818 × 10−5
4	1.450 × 10−5	1.595 × 10−5	2.395 × 10−5
5	1.275 × 10−5	1.373 × 10−5	1.871 × 10−5
6	1.171 × 10−5	1.245 × 10−5	1.601 × 10−5
8	1.054 × 10−5	1.104 × 10−5	1.327 × 10−5
10	9.911 × 10−6	1.028 × 10−5	1.189 × 10−5
12	9.509 × 10−6	9.798 × 10−6	1.106 × 10−5
14	9.232 × 10−6	9.473 × 10−6	1.050 × 10−5
16	9.031 × 10−6	9.236 × 10−6	1.011 × 10−5
18	8.877 × 10−6	9.056 × 10−6	9.893 × 10−6
20	8.755 × 10−6	8.914 × 10−6	9.579 × 10−6
bulk	7.741 × 10−6

As the shape-dependence of solubility is concerned, it is obvious that the larger the shape factor α′ the higher is the CuO solubility. It should also be noted, that the shape factor (and thus the solubility) increases for x < 2 as it follows from Figure 2.

Figure 2

Dependence of differential shape factor for a cylinder on the aspect ratio x = h/ρ according to Equation (7).

Any quantitative comparison with experimentally determined solubility of CuO nanoparticles should be questionable due to various media in which the dissolution experiments have been carried out. The enhanced solubility of nanoparticles in comparison with bulk materials (microparticles) has been observed in a number of studies [17,18,19]. Wang et al. [11] have studied the dissolution of CuO nanoparticles (r = 20 nm) and have observed that approx. 0.25 mg L−1 of Cu released into the suspension of nanoparticles in 3% NaCl solution (320 mg L−1 CuO nanoparticles after 72 h exposition). Our calculated value 0.56 mg L−1 is almost twice greater but there is no evidence of the saturation of solution in the work of Wang et al. [11]. Misra et al. [19] have investigated dissolution of CuO nanospheres (diameter 7 nm) in 1mM NaNO3 solution (pH = 6.7). They have determined equilibrium solubility 1 mg Cu L−1 which is the same as our calculated value 1.71 = mg L−1 at pH = 6.69. They have also examined rod-shape CuO nanoparticles with length h = (40 ± 10) nm and diameter 2ρ = (7 ± 1) nm (x = 11.43, α′ = 1.484). The observed ratio of Cu content in the solutions was approx. 2.2 for nano-sphere and nano-rods. Our calculations for the same conditions give a lower value 1.3. It should be noted, that the same value of interfacial energy for CuO/solution interface was considered, but it is reasonable to suppose that the increase of surface-to-volume ratio for a cylinder might be partly compensated by enhanced exposition of lower energy surfaces in real nanoparticles. This would increase the stability of the nano-rods, reduce their solubility, and increase the calculated ratio of Cu content in solutions for nano-spheres and nano-rods.

5. Conclusions

In this work, we have used the theoretical approach to determine various CuO nanostructures (spheres or cylinders) solubility in an aqueous environment. The equilibrium CuO solubility was calculated using Gibbs energy minimization technique considering thirteen Cu containing aqueous species, out of which only Cu2+, CuOH+, Cu(OH)2, and CuCO3 are significant in neutral or slightly acidic aqueous solutions. For the smallest spherical nanoparticles considered in this work (r = 2 nm), the solubility was 3.8-times the solubility of bulk CuO, suggesting significantly increased solubility due to nanosizing. The increase of solubility was even more considerable in the case of cylindrical nanoparticles. It is a well-known fact that copper oxide nanostructures are potentially toxic to many different organisms. The mechanism of their toxicity is not definitely explained but the strong evidence that the CuO dissolution plays the dominant role has been brought. Our study of CuO nanoparticles dissolution can be thus used for environmental issues, especially for water environment.