1. Introduction
The major drawbacks of the industrial development are the quantity and diversity of wastes which are discharged in the environment. One of the most dangerous groups of inorganic pollutants is represented by the heavy metals since these are not susceptible to biological degradation [1,2].
Porous activated carbons, zeolites, bio-adsorbent materials and carbon nanotubes are extensively used as adsorbents for the removal of heavy metals [3]. However, practical applications of these materials are limited by their low adsorption capacities, low efficiencies, or high cost. With the rapid progress in new material development, metal-organic frameworks (MOFs) and coordination networks (CNs) have received increasing attention in recent years [4]. MOFs are mostly constructed based on metal organic carboxylic derivatives from metal ion nodes linked by organic linkers to form a variety of 1D chain, 2D layer and a three-dimensional (3D) crystal structures with micropores. Phosphonate coordination networks are fast gaining an essential position amongst the families of CNs materials. To expand the adsorption capacity of CNs in a quite large range of pH, it is suggested to choose CNs which are stable in water media [5]. Compared with other adsorbent materials, the main advantages of MOFs and CNs in adsorption processes are their large specific surface area, well ordered unique structures, stable and homogeneous pores of specific sizes. Certainly, MOFs demonstrate good absorbance capacities for a high variety of species that include heavy metals [6,7], drugs [8,9] and dyes [10] from wastewaters.
Cadmium contamination results from many sources such as metal plating, iron and steel production, mining operations, phosphate fertilizer manufacture and use [11]. Even at a low dosage, it can be harmful to both human health and the environment [12]. For this reason, a lot of treatment methods like ion exchange, precipitation, filtration, oxidation-reduction, membrane separation, and adsorption have been developed for the treatment of wastewaters with heavy metal content [13,14,15]. The most economical, feasible and selective method for heavy metal removal from aqueous solutions is the adsorption technique [16,17]. Therefore, researchers are focused on the development of new and more efficient adsorbent materials ranging from natural substances to highly selective synthetic systems to be used as hazardous metal adsorbents [2,18,19,20,21,22,23]. Heavy metal ions, even at small concentrations, are extremely toxic to alive organisms, because they are non-biodegradable, and they tend to accumulate in the environment.
To the best of our knowledge, so far only a few studies have reported the removal of cadmium ions from aqueous solutions throughadsorption onto MOF materials, in which some composites such as cyclodextrin metal-organic framework-based nanoporous carbon [24] and sulfonated MOF loaded onto iron oxide nanoparticles (Fe3O4@MOF235(Fe)–OSO3H were used [25]. The preparation of these materials involves the use of greater quantities of reagents and many preparation steps, which lead to an increase of the production costs.
Taking into account the International Union of Pure and Applied Chemistry (IUPAC) nomenclature and terminology recommendation [26] that coordination networks (CNs) are a subdivision of coordination polymers and MOFs a further subset of coordination networks, we will henceforth name our materials as CNs.
In the present paper, the use of two coordination networks based on cobalt and nickel were used as adsorbent materials in the removal process of cadmium ions from aqueous solutions. The structure, morphology, and properties of materials were investigated by Fourier-transform infrared spectroscopy (FTIR), scanning electron microscopy (SEM) and thermal gravimetric analysis (TGA), which were previously described [27,28]. Inspired by these adsorption properties, we performed PM3 semiempirical analyses of structural properties to predict and understand better some special properties of these compounds. Therefore, octahedral models were calculated for networks containing Ni2+ and Co2+ ions and N,N-bis(phosphonomethyl)glycine. Bond lengths/angles, torsion angles and partial charges for the central metal ions Ni and Co coordination networks are compared.
2. Materials and Methods All chemicals were of reagent grade quality achieved from commercial sources and used without further purification. Ni(CH3COO)2·4H2O and Co(NO3)2·6H2O were purchased from Merck (Milipore, Darmstadt, Germany), N,N-bis(phosphonomethyl)-glycine and Sodium hydroxide (Sigma Aldrich Chemie GmbH (München, Germany) and urea from Alfa Aesar (Karlsruhe, Germany). 2.1. Instrumentation The specific surface area together with a pore volume of Co–Gly and Ni–Gly were measured with an ASAP 2020 BET surface area analyzer (Micrometrics, Micrometrics Instrument Corporation, Norcross, GA, USA) by cold nitrogen adsorption. SEM images were registered with a FEG 250 microscope (Quanta, Field Electron and Iron Company (FEI), Hillsboro, OR, USA), equipped with an EDAX/ZAF quantifier. Cadmium ion concentrations were measured via a SpectrAA 280 FS atomic adsorption spectrophotometer (Varian, Melbourne, Australia). Thermal analysis (TG-DTA) data were recorded on an SDT-Q600 analyzer from TA Instruments (New Castle, DE, USA). A Diamond thermogravimetric analyzer (Perkin Elmer, New York, NY, USA) was used applying temperatures between 30 and 680 °C under a N2 flow increasing the heating at a rate of 10 °C/min. The adsorption studies of investigated materials were made in batch mode using a SW23 shaker bath (Julabo Labortechnik GmbH, Sellbach, Germany). 2.2. Materials Synthesis
A 250 mL Erlenmeyer flask was filled with Ni(CH3COO)2·4H2O (50.0 mmol) or Co(NO3)2·6H2O (50.0 mmol) and bidistilled water (50 mL). The materials were stirred with a constant speed of 1000 rpm until a clear (green or violet) solution was formed. In another flask N,N-bis(phosphonomethyl)-glycine, urea (50.0 mmol), and bidistilled water (50 mL) were mixed in the same conditions till a incolor clear solution was formed. Both solutions were mixed in a 250 mL Erlenmeyer flask and the pH was adjusted to 4.5 in the case of Ni containing synthesis and 2.8 in the case of Co containing synthesis with an aqueous solution of NaOH (0.1 M). Then the clear green or violet solution was heated in an oil-bath at 80°C for 75 h, unperturbed. After 75 h heating crystalline green (Ni–Gly) and violet crystals (Co–Gly) materials precipitated and were isolated by filtration and finally air dried (yield: 52–75%) [27,28,29].
2.3. Adsorption Studies All the adsorption studies were conducted in batch mode. In the first step the influence of the pH upon the adsorption capacity of Co–Gly and Ni–Gly was determined. For each experiment 25 mL of a solution containing 30 mL of cadmium ions were treated with 0.05 g of adsorbent material for 1 h at a constant speed of 200 rpm, using a Julabo SW23 shaker bath. After 1h of reaction, the samples were filtered and the residual concentration of Cd(II) ions was analyzed in the filtrate by atomic adsorption spectrophotometer. The pH adjustment of the solution was done using 1.0 M NaOH or 1.0 M HCl and was measured using a pH-meter (Mettler Toledo, Giessen-Germany).
The adsorption capacity of the studied materials in the removal process of Cd(II) was calculated according to the mass balance (Equation (1)):
q=(C0−Ce)·Vm
where q is the amount of Cd(II) adsorbed (mg/g); C0 and Ce represents the initial and equilibrium concentration of Cd(II) in the solutions (mg/L), respectively. V represents the solution volume (L) and m represents the adsorbent mass (g) used in the experiments.
To study the effect of contact time on adsorption, further experiments were carried out using the same S:L ratio, the same concentration in solution of Cd, an initial pH of the solutions equal to 5, but the suspension were kept in contact for different times (15–120 min) at 25 °C. After the contact time had passed, the suspensions were filtered and the liquid was collected for analysis of the residual concentration of cadmium. Pseudo-first and second order kinetic models were applied to estimate the adsorption rate constants and the adsorbent mechanism. The influence of the initial concentration of Cd(II) ions upon the adsorption capacity of Co–Gly and Ni–Gly was measured using the same S:L ratio at different initial concentrations (range: 5–300 mg/L). The non-linearized isotherm models of Langmuir, Freundlich, and Redlich-Peterson were employed to correlate the experimental adsorption data. Their adsorption capacities have been studied as a function of pH, contact time and cadmium initial concentration. The studied CNs were regenerated with HCl solution 0.2 M having an initial pH = 2. For the recovery of Cd ions from the CNs surface, a S:L ratio of 1 g/L was used, and the samples were mixed for 15 min. After the regeneration process the phases were separated the recycled adsorbent was used in other adsorption process and the extracted Cd ions from the solution were determined. The materials were used in five adsorption-desorption process cycles. 2.4. Computational Studies
Two coordination networks—Ni–Gly and Co–Gly—containing the basic units [Ni(HO3PCH2)2N(H)CH2COO)(H2O)2] and [Co(HO3PCH2)2N(H)CH2COO)(H2O)2], respectively, were built and visualized with the aid of Mercury 4.1.3 software (Cambridge Crystallographic Data Centre, Cambridge, UK) [30]. The generated CNs were geometrically optimized by means of the semi-empirical PM3-RHF method implemented in HyperChem version 7.52 (Hypercube, Inc., Gainesville, FL, USA) [31] software. Polak-Ribiere conjugate gradient algorithm and a RMS gradient norm limit of 10−2 kcal/A were used, while the self-consistent field (SCF) convergence criterion was considered 10−5. Maestro version 12.0.012 from the Schrodinger package was used for the computation of the surface areas. [https://www.schrodinger.com/maestro].
3. Results and Discussion 3.1. Materials Characterization
The morphology of the synthesized CNs is presented in Figure 1. It can be observed that the CNs based on Ni ions present a more ordered structure, with particles of well-defined sizes and shapes compared with the CN based on Co ions, which surface is more compact with particle conglomerates of various sizes and shapes.
The specific surface area and the pore volume of the synthesized CNs are presented in Table 1. It can be observed that the Ni–Gly sample presents a higher specific surface area and a higher pore volume compared with Co–Gly sample. In accordance with the results of the characterization studies, due to its structure and morphology, it is expected that Ni–Gly to develop higher adsorption capacity in the removal process of Cd ions from aqueous solutions compared to Co–Gly sample.
TGA data for Ni–Gly (Ni(C4H9O8NP2)·2H2O) and Co–Gly (Co(C4H9O8NP2)·2H2O) shows that the removal of water molecules starts almost immediately upon heating and is lost slowly between 290 °C and 370 °C followed by decomposition at ~400 °C. The total weight loss caused by decomposition of all the organic composition of Ni–Gly and Co–Gly is around 35% for the former and 38% in case of Co–Gly and ocurrs at approximately 700 °C, suggesting an endothermic process (Figure 2). Detailed X-Ray Diffraction (XRD)characterisation studies for Co–Gly and Ni–Gly are under way.
3.2. The pH Influence upon the Adsorption Studies
The solutions’ pH affects the properties and the degree of protonation of the adsorbent surface. Due to the fact that at higher values of pH the Cd(II) ions could precipitate the studies were carried out in the 2–8 pH range. The experimental data regarding the dependence of Cd(II) ions adsorbed by the studied materials as a function of the initial pH of the solutions are presented in Figure 3.
The initial pH of Cd(II)-containing solutions has a significant effect upon the adsorption performance of the studied materials, displaying a maximum adsorption capacity at an initial pH of 5. The adsorption capacity decreases with the change of pH around this value. This behavior can be explained by the surface loading of the adsorbent material and by the competitive adsorption of protons [17,32]. At lower pH values the adsorbent surface is positively charged and therefore there is an electrostatic repulsion between the adsorbent surface and cadmium cations [33,34]. In the same time at higher pH values, Cd ions could precipitate under Cd(OH)2 and then the adsorption is inhbitated [33,35,36]. Further experiments were carried out with Cd(II) solution having an initial pH of 5.
3.3. Kinetics of Adsorption
Figure 4 shows the kinetics of Cd(II) adsorption onto Co–Gly and Ni–Gly, respectively. Adsorption is fast, in both cases, and the equilibrium between the adsorbent and adsorbate was achieved after 60 min. The kinetic data were fitted by non-linear regression using the Lagergren (pseudo-first order kinetic model) and Ho and McKay (pseudo-second order kinetic model) equations [14,17]. Table 2 presents the calculated parameters and the correlation coefficients obtained after fitting.
The regression coefficients R2 showed that the Ho and McKay model fitted the kinetic behaviour of the process when Co–Gly and Ni–Gly were used as adsorbents. For both adsorbents, the adsorption capacities calculated at equilibrium are in agreement with those experimental values obtained. The adsorption of Cd(II) ions onto Co–Gly and Ni–Gly has a chemo-sorption profile. 3.4. Equilibrium of Adsorption
The equilibrium adsorption data of Cd(II) onto Co–Gly and Ni–Gly were analyzed by using the Langmuir, Freundlich and Redlich-Peterson models and non-linear analysis in order to predict the overall adsorption behavior. The isotherm parameters obtained after fitting the experimental data for the adsorption of Cd(II) on the two materials are presented in Table 3. Figure 5 presents the experimental data and the isotherms obtained by simulations of the mathematical models used.
As seen in Table 3, the Langmuir model fitted the data the best over the whole concentration range. The maximum adsorption capacities experimentally obtained were 48.2 mg/g for Cd(II) adsorption onto Co–Gly and 55 mg/g for Cd(II) adsorption onto Ni–Gly, respectively. These values are close to the maximum adsorption capacities obtained when the data are fitted by a Langmuir model (qmax = 51.5 mg/g for Co–Gly and 58.1 mg/g for Ni–Gly). From Table 3 it can be observed that the non-homogeneity factor n in Redlich-Peterson model has values closes to 1, so that the behavior of the samples obeys a Langmuir model. The Langmuir isotherm idea involves a monolayer coverage of adsorbate above a homogeneous adsorbent surface [17]. The essential characteristics of the Langmuir isotherm is communicated in relations of a dimensionless constant separation factor RL that is specified by the following Equation (2):
RL=11+KL C0
where KL is the Langmuir constant and C0 is the initial concentration of Cd(II) ions. The value of the separation parameter RL offers important data about the type of adsorption. The value of RL point out the category of Langmuir isotherm to be irreversible (RL = 0), favorable (0< RL< 1), linear (RL = 1), or unfavorable (RL > 1) [32]. The RL was established to be between 0 and 1 for the entire concentration interval, and for both studied material which indicates the favorable adsorption of cadmium onto the studied materials.
The studied CNs were used in five adsorption-desorption process cycles and it was observed that their adsorption capacity remains constant for four adsorption-desorption cycles, then it decreases by 20% because the recovery of Cd ions from the CNs’ surface decreases (Figure 6). A decreasing Cd ion recovery capacity means that the available sites for adsorption decrease, and for this reason the adsorption capacity decreased after four adsorption-desorption process cycles.
The maximum adsorption capacity achieved by the studied materials in the removal process of Cd(II) ions from aqueous solutions were compared with the maximum adsorption capacities obtained using other adsorbents and reported in the specialty literature. The results are presented in Table 4. It can be observed that the coordination networks present a higher efficiency in the removal process of Cd(II) ions from aqueous solutions than other low cost adsorbent materials. It could be observed that the CN-based materials reported until now in the literature, displayed higher adsorption capacities in the removal process of Cd ions from aqueous solutions, but in these cases they involved some expensive composite materials, not only CNs. Therefore, a synergistic effect of the synthesized CNs and other materials such as nanoporous carbon or iron oxide nanoparticles from the composite structures could be proposed.
3.5. Computational Semiempirical Studies
In order to design and geometrically optimize the following CN models: [Ni2((HO3PCH2)2N(H)CH2COO))3*2H2O]2− and [Co2((HO3PCH2)2N(H)CH2COO))3*2H2O]2− in silico methods were applied. The optimized structures of these networks are presented in Figure 7.
N,N′-bis-phosphonomethylglycine (Gly) in a coordination complex with Mg ion has been synthesized and structurally characterized by Demadis and co-workers [44] and was used as model structure (CCDC Reference Code 729893). N,N′-bis-phosphonomethylglycine (Gly) in coordination complex with Mg present a 2D layered architecture [44] as can be seen in Figure 8.
From Figure 8 it can be seen that the oxygen atoms belonging to the carboxylate and phosphonate groups hold the layers together by H bonding. Each layer is formed by Mg−O (phosphonate) bonds. It has been observed that the slightly distorted octahedral geometry of the Mg central ion is also maintained in the case of Ni/Co networks. The Ni2+/Co2+ ions coordinate with six oxygen atoms, two belonging to water molecules and four pertaining to phosphonate groups of the Gly ligands (Figure 7). The two axially-oriented oxygen atoms belong to phosphonate moieties, while of the oxygen atoms occupying the four equatorial positions of the octahedral geometry, two belong to phosphonate moieties, and the other two appertain to water molecules (Figure 7). As can be seen in Figure 7, the networks have rings (cavities) consisting of eight atoms; one metal ion (Co2+ or Ni2+) and one protonated nitrogen binding between them two -O-PO2H-CH2 radicals. The network cavities may confer various practical application to these materials, such as, gas storage or different atoms/molecules adsorption properties.
In order to measure and compare the geometrical parameters of the CNs, the most important atoms of the networks were numbered (Figure 9). Co–Gly is exemplified and the atom’s numbers are equivalent for all considered networks. Bond lengths/angles, torsion angles and partial charges for the numbered atoms of the Ni–Gly and Co–Gly are presented in Table 5 and Table 6.
Considering the crystallographic structure of Mg complexed with N,N′-bis-phosphonomethyl-glycine (Mg–Gly) (Figure S1) as the most similar network to our models, some of their geometrical features were analyzed and compared. Thus, the bond length values for Me-O are in the range of 1.8–2.0 Å for Ni and Co compared with 2.0–2.1 Å resulted for Mg (Table S1). The O=P bond (from the cycle) has differences from up to 0.3 Å between Co and Mg; and up to 0.26 Å, respectively, for Ni and Mg CNs. The tendency of the length differences remains about the same for O-P bond (from the ring), shrinking the difference between Co–Gly and Mg–Gly to 0.14 Å. On the other hand, The O=P bond (from the outer phosphonate group) presents similar values in a range of 1.48–1.50 Å for all three networks. the calculated N–C and P–C bond lengths of around 1.5 and 1.9 Å, respectively, are also quite similar to the experimentally determined ones.
There is evidence of the importance of the atomic charges for modeling metal-organic frameworks [45]. The electrostatic attraction between the negatively charged CNs (Me–Gly) and the positively charged heavy metal ions (Cd2+) is the main factor that causes adsorbtion [46]. Due to the charge transfer from oxygen atoms to metal ions, Ni and some Co ions tend to have negative partial charges (Table 5 and Table 6). Thus, Ni27 and Ni56 ions have negative partial charge values of −0.555 and −0.594, respectively. As one can see, a negative value of −0.0996 was obtained for Co56, while the central ion Co27 has a positive charge of 0.075 (Table 5 and Table 6).
The metal ions (Me56) which belong to the marginal rings of eight atoms are in particular the ones with negative charges and they present increased chances of giving up electrons. The surface areas for these both marginal cycles (of Ni–Gly and Co–Gly) were computed and the results show a higher value of 39.924 for Ni–Gly compared with 30.547 obtained for Co–Gly. Taking into account these informations and compare it with the absorbance affinity of Co–Gly and Ni–Gly for Cd(II) ions we can assign this effect to the accentuated negativity of the both Ni ions present in the network higher than the central Co ions as well as due to the bigger area surface for the marginal Ni-ring.
A plethora of potential properties such as adsorption, gas storage, heterogeneous catalysis, separation, ion exchange, magnetism, and sensors, can be explained by the negative charge values of the ions in these CN materials [46,47].
The adsorption efficiency of CNs is well corelated with surface area, pore size and distribution. [21]. In order to evaluate the ability of a chemical structure to donate electrons, the investigation of the highest occupied molecular orbitals (HOMO) values and their localization is of great interest.
From the orbitals component analysis (Figure 10) it can be observed that the HOMOs are located over the marginal cycles of eight atoms, mainly on metal ions (Me56). These findings are in accordance with the aforementioned results, which attested that negative charges of the Co56/Ni56 ions present increased chances of giving up electrons to the adsorbants. The Langmuir isotherm concept adopts monolayer coverage of adsorbate above a homogeneous adsorbent surface.
A series of electronic properties (heat of formation, free energy, vibrational zero-point energy, minimum and maximum fundamental vibrations, frontier orbitals and the energy gap between them) resulted from semiempirical PM3 calculations are presented in Table 7. The positive values of fundamental vibration (νmin) regarding both CN models show that these semiempirical calculated geometries are not transition states, confirming their stability.
The energy gap between the frontier orbitals (LUMO and HOMO) can be used to estimate the strength and stability of coordination networks. A significant band gap correlates with structural and kinetic stability [48]. Another measure for chemical stability is the heat of formation (∆Hform), the lower is the value, the more stable is the complex. Analyzing Table 7 it can be observed that these calculated properties advocates for a good stability of the models, in special for Co–Gly network. Like other studies [49], our findings show that the semiempirical methods are valuable tools to be used in adsorption processes involving CNs.
4. Conclusions In this study, CNs were used as adsorbent materials for the removalof Cd(II) ions from aqueous solutions. The present investigations showed that the studied CNs have a good affinity for the removal of Cd(II) ions from aqueous solution compared with other materials reported in the specialty literature. By applying suitable kinetic models to the experimental data it was established that the adsorption of Cd(II) ions on studied material is defined by a pseudo-second-order kinetic model. The equilibrium sorption data were modeled using Langmuir, Freundlich, and Redlich-Peterson isotherms and the first one provided an excellent fit of the experimental data, giving a maximum adsorption capacity of 51.5 mg/g and 58.1 mg/g for Co–Gly and Ni–Gly, respectively. The higher specific surface area and pore volume of Ni–Gly together with the higher negative partial charges of Ni in the network shown by PM3 semiempirical computations increase the electrostatic attraction between the Ni–Gly and the positively charged heavy metal ions (Cd2+). This finding is the key factor that causes higher adsorbtion capacity for Ni–Gly. Thus, Ni–Gly is the adsorbant material with a greater adsorbance capacity, as confirmed both by experimental and theoretical methods.
Enlarge this image.Figure 1. Scanning electron microscopy (SEM) images of the synthesized coordination networks (CNs) (a) Co-Gly; (b) Ni-Gly.
Enlarge this image.Figure 3. pH effect upon the adsorption capacity of the studied materials in the removal process of Cd(II) ions from aqueous solutions.
Enlarge this image.Figure 6. The adsorption performance of the studied CNs in various adsorption-desorption cycles (a) the adsorption capacity, after each cycle (b) Cd recovery, after each cycle.
Enlarge this image.Figure 8. Partial view of two adjacent layers in the crystal structure of Mg-Gly model structure.
Enlarge this image.Figure 9. Co-Gly atoms numbering. For simplification only the numbers of significant atoms around the metal ions were represented.
Enlarge this image.Figure 10. The highest occupied molecular orbital (HOMO) components (a) Ni-Gly; (b) Co-Gly.
| Adsorbent | Specific Surface Area, m2/g | Pore Volume, cm3/g |
|---|---|---|
| Co–Gly | 32 | 0.25 |
| Ni–Gly | 45 | 0.85 |
| Adsorbent | qexp | Lagergren Modelqt=qe(1−e−k1t) | Ho and McKay Modelqt=qe(1−1k2 qet+1) | ||||
|---|---|---|---|---|---|---|---|
| qe (mg/g) | k1 (1/min) | R2 | qe (mg/g) | k2 (g/mg·min) | R2 | ||
| Co–Gly | 13.2 | 8.39 | 0.0302 | 0.8146 | 16.36 | 2.5 × 10-3 | 0.9968 |
| Ni–Gly | 15.2 | 10.76 | 0.0357 | 0.9848 | 18.14 | 1.11 × 10-3 | 0.9968 |
| Adsorbent | Langmuir Modelqe=qmaxKL·Ce1+KL·Ce | Freundlich Modelqe=KF·Ce1/n | ||||||
|---|---|---|---|---|---|---|---|---|
| qmax (mg/g) | KL | R2 | KF | 1/n | R2 | |||
| Co–Gly | 51.5 | 0.0739 | 0.9994 | 4.47 | 0.5143 | 0.9391 | ||
| Ni–Gly | 58.1 | 0.0818 | 0.9968 | 5.54 | 0.4972 | 0.9406 | ||
| Adsorbent | Redlich-Peterson Modelqe=qmaxKRP·Ce1+(KRP·Ce)n | |||||||
| qmax (mg/g) | KRP | n | R2 | |||||
| Co–Gly | 0.485 | 4.47 | 1.02 | 0.9312 | ||||
| Ni–Gly | 0.503 | 5.54 | 0.952 | 0.9419 | ||||
| Adsorbent | qm, mg/g | pH | References |
|---|---|---|---|
| Orange Peels | 4.9 | 5 | [37] |
| Kaolinite | 7.407 | 8 | [38] |
| Metakaolinite | 9.174 | 8 | [39] |
| Cyclodextrin metal-organic framework based nanoporous carbon | 140.85 | 7 | [24] |
| Sulfonated metal organic framework loaded on iron oxide nanoparticles | 163.9 | 3 | [25] |
| Grass char | 115.8 | 6.8 | [39] |
| Kaolin | 1.46 | 6.8 | [40] |
| Sediments | 10.01 | 5.5 | [41] |
| Zerovalent iron particles | 714.3 | - | [42] |
| Natural cheese | 5.12 | 6 | [43] |
| Co–Gly | 51.5 | 5 | Present paper |
| Ni–Gly | 58.1 | 5 |
| Ni–Gly | ||||||||
|---|---|---|---|---|---|---|---|---|
| Atom | ID | Charge | Bond | Distance | Bond Angle | Degree | Torsion Angle | Degree |
| Carbon | C15 | −0.723 | Ni56-O29 | 1.815 | O29-Ni56-O36 | 94.636 | C46-P25-O29-Ni56 | 9.656 |
| Carbon | C18 | −0.697 | Ni56-O36 | 1.792 | Ni56-O36-P26 | 138.283 | P25-O29-Ni56-O36 | −92.872 |
| Carbon | C43 | −0.599 | O36-P26 | 1.853 | O36-P26-C43 | 93.346 | O29-Ni56-O36-P26 | 130.169 |
| Carbon | C46 | −0.807 | P26-C43 | 1.986 | P26-C43-N30 | 124.21 | Ni56-O36-P26-O32 | 41.057 |
| Carbon | C74 | −0.763 | C43-N30 | 1.493 | C43-N30-C46 | 116.868 | O36-P26-O32-Ni27 | 1.422 |
| Hydrogen | H34 | 0.2853 | N30-C46 | 1.502 | N30-C46-P25 | 130.108 | P26-O32-Ni27-O4 | −94.052 |
| Hydrogen | H35 | 0.2188 | C46-P25 | 1.935 | C46-P25-O29 | 113.781 | P26-O32-Ni27-O8 | 170.186 |
| Hydrogen | H54 | 0.2133 | P25-O37 | 1.467 | P25-O29-Ni56 | 95.464 | P26-O32-Ni27-O33 | −2.211 |
| Hydrogen | H55 | 0.2568 | P25-O38 | 1.700 | O4-Ni27-O8 | 97.744 | P26-O32-Ni27-O65 | 69.181 |
| Nitrogen | N5 | 0.6925 | P26-O41 | 1.766 | O8-Ni27-O53 | 89.307 | O32-Ni27-O4-P1 | −150.244 |
| Nitrogen | N30 | 0.7692 | P26-O32 | 1.786 | O53-Ni27-O65 | 85.753 | O32-Ni27-O8-P2 | 28.512 |
| Nickel | Ni27 | −0.555 | O32-Ni27 | 1.842 | O65-Ni27-O33 | 71.3 | Ni27-O8-P2-C15 | 53.962 |
| Nickel | Ni56 | −0.594 | Ni27-O4 | 1.869 | O33-Ni27-O32 | 93.215 | O4-Ni27-O8-P2 | −62.239 |
| Oxygen | O4 | −0.567 | O4-P1 | 1.705 | O32-Ni27-O4 | 90.604 | O32-Ni27-O65-P57 | 70.339 |
| Oxygen | O7 | −0.835 | P1-O9 | 1.492 | O4-Ni27-O53 | 91.125 | ||
| Oxygen | O8 | −0.433 | P1-O10 | 1.634 | O4-Ni27-O65 | 162.982 | ||
| Oxygen | O9 | −0.886 | P1-C18 | 1.912 | O4-Ni27-O33 | 91.804 | ||
| Oxygen | O10 | −0.670 | C18-N5 | 1.499 | O8-Ni27-O65 | 100.934 | ||
| Oxygen | O13 | −0.700 | N5-C15 | 1.494 | O8-Ni27-O33 | 171.094 | ||
| Oxygen | O29 | −0.476 | C15-P2 | 1.898 | O8-Ni27-O32 | 91.428 | ||
| Oxygen | O32 | −0.323 | P2-O7 | 1.456 | O53-Ni27-O33 | 85.817 | ||
| Oxygen | O33 | 0.2309 | P2-O13 | 1.634 | O53-Ni27-O32 | 178.044 | ||
| Oxygen | O36 | −0.415 | P2-O8 | 1.743 | O65-Ni27-O32 | 92.328 | ||
| Oxygen | O37 | −0.873 | Ni27-O8 | 1.836 | ||||
| Oxygen | O38 | −0.699 | Ni27-O53 | 1.925 | ||||
| Oxygen | O41 | −0.650 | O53-H54 | 0.959 | ||||
| Oxygen | O53 | −0.085 | O53-H55 | 0.986 | ||||
| Oxygen | O60 | −0.852 | Ni27-O33 | 1.947 | ||||
| Oxygen | O65 | −0.604 | O33-H34 | 1.012 | ||||
| Oxygen | O66 | −0.692 | O33-H35 | 0.979 | ||||
| Phosphorus | P1 | 2.0639 | Ni27-O65 | 1.853 | ||||
| Phosphorus | P2 | 2.0923 | O65-P57 | 1.733 | ||||
| Phosphorus | P25 | 2.1319 | P57-O60 | 1.450 | ||||
| Phosphorus | P26 | 1.1473 | P57-O66 | 1.683 | ||||
| Phosphorus | P57 | 2.0932 | P57-C74 | 1.848 | ||||
| Co–Gly | ||||||||
|---|---|---|---|---|---|---|---|---|
| Atom | ID | Charge | Bond | Distance | Bond Angle | Degree | Torsion Angle | Degree |
| Carbon | C15 | −0.7165 | Co56-O29 | 1.874 | O29-Co56-O36 | 103.345 | C46-P25-O29-Co56 | −57.484 |
| Carbon | C18 | −0.6848 | Co56-O36 | 2.025 | Co56-O36-P26 | 136.025 | P25-O29-Co56-O36 | −10.336 |
| Carbon | C43 | −0.731 | O36-P26 | 1.560 | O36-P26-C43 | 104.476 | O29-Co56-O36-P26 | 78.858 |
| Carbon | C46 | −0.7151 | P26-C43 | 1.892 | P26-C43-N30 | 121.608 | Co56-O36-P26-O32 | 36.937 |
| Carbon | C74 | −0.7589 | C43-N30 | 1.494 | C43-N30-C46 | 114.412 | O36-P26-O32-Co27 | −16.624 |
| Hydrogen | H34 | 0.209 | N30-C46 | 1.503 | N30-C46-P25 | 127.21 | P26-O32-Co27-O4 | −104.721 |
| Hydrogen | H35 | 0.2187 | C46-P25 | 1.956 | C46-P25-O29 | 102.822 | P26-O32-Co27-O8 | 159.331 |
| Hydrogen | H54 | 0.2528 | P25-O37 | 1.463 | P25-O29-Co56 | 139.07 | P26-O32-Co27-O33 | −18.618 |
| Hydrogen | H55 | 0.2879 | P25-O38 | 1.707 | O4-Co27-O8 | 96.057 | P26-O32-Co27-O65 | 75.019 |
| Nitrogen | N5 | 0.7035 | P26-O41 | 1.641 | O8-Co27-O53 | 83.833 | O32-Co27-O4-P1 | −138.862 |
| Nitrogen | N30 | 0.7478 | P26-O32 | 1.659 | O53-Co27-O65 | 88.608 | O32-Co27-O8-P2 | 17.287 |
| Cobalt | Co27 | 0.075 | O32-Co27 | 1.918 | O65-Co27-O33 | 93.076 | Co27-O8-P2-C15 | 46.317 |
| Cobalt | Co56 | −0.0996 | Co27-O4 | 1.917 | O33-Co27-O32 | 97.175 | O4-Co27-O8-P2 | −69.782 |
| Oxygen | O4 | −0.7808 | O4-P1 | 1.675 | O32-Co27-O4 | 86.851 | O32-Co27-O65-P57 | −105.659 |
| Oxygen | O7 | −0.8435 | P1-O9 | 1.504 | O4-Co27-O53 | 90.374 | ||
| Oxygen | O8 | −0.7288 | P1-O10 | 1.642 | O4-Co27-O65 | 178.93 | ||
| Oxygen | O9 | −0.8995 | P1-C18 | 1.913 | O4-Co27-O33 | 86.533 | ||
| Oxygen | O10 | −0.6691 | C18-N5 | 1.501 | O8-Co27-O65 | 84.167 | ||
| Oxygen | O13 | −0.6707 | N5-C15 | 1.493 | O8-Co27-O33 | 170.364 | ||
| Oxygen | O29 | −0.7278 | C15-P2 | 1.921 | O8-Co27-O32 | 92.241 | ||
| Oxygen | O32 | −0.7066 | P2-O7 | 1.461 | O53-Co27-O33 | 86.876 | ||
| Oxygen | O33 | −0.0615 | P2-O13 | 1.68 | O53-Co27-O32 | 174.936 | ||
| Oxygen | O36 | −0.7897 | P2-O8 | 1.669 | O65-Co27-O32 | 94.188 | ||
| Oxygen | O37 | −0.8422 | Co27-O8 | 1.888 | ||||
| Oxygen | O38 | −0.6909 | Co27-O53 | 1.972 | ||||
| Oxygen | O41 | −0.6614 | O53-H54 | 0.965 | ||||
| Oxygen | O53 | −0.2364 | O53-H55 | 0.995 | ||||
| Oxygen | O60 | −0.8173 | Co27-O33 | 2.000 | ||||
| Oxygen | O65 | −0.7366 | O33-H34 | 0.978 | ||||
| Oxygen | O66 | −0.5273 | O33-H35 | 0.973 | ||||
| Phosphorus | P1 | 2.0692 | Co27-O65 | 1.936 | ||||
| Phosphorus | P2 | 2.0995 | O65-P57 | 1.632 | ||||
| Phosphorus | P25 | 2.0392 | P57-O60 | 1.451 | ||||
| Phosphorus | P26 | 1.9959 | P57-O66 | 1.823 | ||||
| Phosphorus | P57 | 2.0514 | P57-C74 | 1.853 | ||||
| CNs/Electronic Properties | ∆Hform (kcal/mol) | Free Energy (kcal/mol) | Zero-pt Vib Energy (kcal/mol) | νmin (cm−1) | νmax (cm−1) | HOMO (eV) | LUMO (eV) | LUMO—HOMO (eV) |
|---|---|---|---|---|---|---|---|---|
| Co–Gly | −1990.3819 | −284904 | 342.5893 | 16.06 | 3959.68 | −4.7667 | 3.2892 | 8.0559 |
| Ni–Gly | −1828.6042 | −296474 | 345.3689 | 9.67 | 3963.54 | −4.6634 | 2.7312 | 7.3946 |
Supplementary Materials
The following are available online at https://www.mdpi.com/2079-4991/10/5/899/s1, Figure S1: Superposition of the crystal structure Mg-Gly over the optimized structure of Mg-Gly, Table S1: Geometric properties of the Mg-Gly model.
Author Contributions
A.V. and B.M. synthesized and characterized the adsorbent materials; L.L. performed the adsorbance studies and analyzed the results; A.B. performed the computational studies, L.C. accomplished some computational measurements and edited the manuscript. All authors discussed the outcomes and commented on the paper. All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.
Funding
This work was partially supported by Program no 1 and 2, from the "Coriolan Dragulescu" Institute of Chemistry Timisoara, Romania and by a grant of the Romanian National Authority for Scientific Research, CNCS-UEFISCDI, project number PN-III-P1-1.1- TE-2016-2008, within PN III.
Acknowledgments
We thank to Mircea and Maria Mracec for access to the HyperChem 7.52 The authors wish to thank Schrodinger Inc for providing an academic trial license to complete the calculation for this paper.
Conflicts of Interest
The authors declare no conflict of interest.
1. Shafique, U.; Anwar, J.; Anzano, J.M.; Salman, M. Hazardous Metals: Sources, Toxicity, and Control; Nova Science Publisher Inc.: New York, NY, USA, 2011; pp. 157-186.
2. Ciesielczyk, F.; Bartczak, P.; Jesionowski, T. Removal of cadmium(II) and lead(II) ions from model aqueous solutions using sol-gel-derived inorganic oxide adsorbent. Adsorption 2016, 22, 445-458.
3. Fu, F.L.; Wang, Q. Removal of heavy metal ions from wastewaters: A review. J. Environ. Manag. 2011, 92, 407-418.
4. Giliopoulos, D.; Zamboulis, A.; Giannakoudakis, D.; Bikiaris, D.; Triantafyllidis, K. Polymer/Metal Organic Framework (MOF) Nanocomposites for Biomedical Applications. Molecules 2020, 25, 185.
5. Wang, X.; Chen, Y.; Zheng, Y.; Zhang, Z. Study of Adsorption and Desorption Performances of Zr-Based Metal-Organic Frameworks Using Paper Spray Mass Spectrometry. Materials 2017, 10, 769.
6. Zhuang, P.; Zhang, P.; Li, K.; Kumari, B.; Li, D.; Mei, X. Silver Nanoclusters Encapsulated into Metal-Organic Frameworks for Rapid Removal of Heavy Metal Ions from Water. Molecules 2019, 24, 2442.
7. Chowdhury, T.; Zhang, L.; Zhang, J.; Aggarwal, S. Removal of Arsenic (III) from Aqueous Solution Using Metal Organic Framework-Graphene Oxide Nanocomposite. Nanomaterials 2018, 8, 1062.
8. Wang, H.; Zhang, X.; Wang, Y.; Quan, G.; Han, X.; Yan, J. Facile synthesis of magnetic nitrogen-doped porous carbon from bimetallic metal-organic frameworks for efficient norfloxacin removal. Nanomaterials 2018, 8, 664.
9. Luo, Z.; Fan, S.; Liu, J.; Liu, W.; Shen, X.; Wu, C.; Huang, Y.; Huang, G.; Huang, H.; Zheng, M. A 3D stable metal-organic framework for highly efficient adsorption and removal of drug contaminants from water. Polymers 2018, 10, 209.
10. Wei, F.; Chen, D.; Liang, Z.; Zhao, S. Comparison study on the adsorption capacity of Rhodamine B, Congo Red, and Orange II on Fe-MOFs. Nanomaterials 2018, 8, 248.
11. Hutton, M. Sources of cadmium in the environment. Ecotoxicol. Environ. Saf. 1983, 7, 9-24.
12. Anda, A.; Illésand, B.; Soós, G. Effect of cadmium pollution of atmospheric origin on field-grown maize in two consecutive years with diverse weather conditions. Acta Biol. Hung. 2013, 64, 476-489.
13. Mopoung, R.; Kengkhetkit, N. Lead and Cadmium Removal Efficiency from Aqueous Solution by NaOH Treated Pineapple Waste. Int. J. Appl. Chem. 2016, 12, 23-35.
14. Mondal, M.K.; Mishra, G.; Kumar, P. Adsorption of cadmium (II) and chromium (VI) from aqueous solution by waste marigold flowers. J. Sustain. Dev. Energy Water Environ. Syst. 2015, 3, 405-415.
15. Canet, L.; Ilpide, M.; Seta, P. Efficient facilitated transport of lead, cadmium, zinc and silver across a flat sheet-supported liquid membrane mediated by a lasalocid A. Sep. Sci. Technol. 2002, 37, 1851-1860.
16. Cheraghi, E.; Ameri, E.; Moheb, A. Adsorption of cadmium ions from aqueous solutions using sesame as a low-cost biosorbent: Kinetics and equilibrium studies. J. Environ. Sci. Technol. 2015, 12, 2579-2592.
17. Kumar, P.S.; Ramakrishnan, K.; Kirupha, S.D.; Sivanesan, S. Thermodynamic and kinetic studies of cadmium adsorption from aqueous solution onto rice husk. Braz. J. Chem. Eng. 2010, 27, 347-355.
18. Cutillas-Barreiro, L.; Ansias-Manso, L.; Fernández-Calviño, D.; Arias-Estévez, M.; Nóvoa-Muñoz, J.C.; Fernández-Sanjurjo, M.J.; Álvarez-Rodríguez, E.; Núñez-Delgado, A. Pine bark as bio-adsorbent for Cd, Cu, Ni, Pb and Zn: Batch-type and stirred flow chamber experiments. J. Environ. Manag. 2014, 144, 258-264.
19. Santhosh, C.; Kollu, P.; Felix, S.; Velmurugan, V.; Jeong, S.K.; Grace, A.N. CoFe2O4 and NiFe2O4@graphene adsorbents for heavy metal ions-kinetic and thermodynamic analysis. RSC Adv. 2015, 5, 28965-28972.
20. Wysokowski, M.; Klapiszewski, Ł.; Moszyński, D.; Bartczak, P.; Szatkowski, T.; Majchrzak, I.; Siwińska-Stefańska, K.; Bazhenov, V.V.; Jesionowski, T. Modification of chitin with kraft lignin and development of new biosorbents for removal of cadmium(II) and nickel(II) ions. Mar. Drugs 2014, 12, 2245-2268.
21. Li, J.; Wang, X.; Zhao, G.; Chen, C.; Chai, Z.; Alsaedi, A.; Hayat, T.; Wang, X. Metal-organic framework-based materials: Superior adsorbents for the capture of toxic and radioactive metal ions. Chem. Soc. Rev. 2018, 47, 2322-2356.
22. Manousi, N.; Giannakoudakis, D.A.; Rosenberg, E.; Zachariadis, G.A. Extraction of Metal Ions with Metal-Organic Frameworks. Molecules 2019, 24, 4605.
23. Rocío-Bautista, P.; Pacheco-Fernández, I.; Pasán, J.; Pino, V. Are Metal-Organic Frameworks able to provide a new generation of solid-phase microextraction coatings?-A review. Anal. Chim. Acta 2016, 939, 26-41.
24. Liu, C.; Wang, P.; Liu, X.; Yi, X.; Liu, D.; Zhou, Z. Ultrafast removal of Cd(II) by green cyclodextrin metal-organic framework based nanoporous carbon: Adsorption mechanism and application. Chem. Asian J. 2019, 14, 261-268.
25. Moradi, S.E.; Dadfarnia, S.; Emami, S.; Shabani, A.M.H. Sulfonated metal organic framework loaded on iron oxide nanoparticles as a new sorbent for the magnetic solid phase extraction of cadmium from environmental water samples. Anal. Methods 2016, 8, 6337-6346.
26. Batten, S.R.; Champness, N.R.; Chen, X.M.; Garcia-Martinez, J.; Kitagawa, S.; Öhrström, L.; O'Keeffe, M.; Suh, P.M.; Reedijk, J. Therminology of metal-organic frameworks and coordination polymers (IUPAC Recommendations 2013). Pure Appl. Chem. 2013, 85, 1715-1724.
27. Maranescu, B.; Lupa, L.; Visa, A. Heavy metal removal from waste waters by phosphonate metal organic frameworks. Pure Appl. Chem. 2018, 90, 35-47.
28. Maranescu, B.; Lupa, L.; Visa, A. Synthesis, characterizations and Pb(II) sorption properties of cobalt phosphonate materials. Pure Appl. Chem. 2016, 88, 979-992.
29. Mateescu, A.; Gabriel, C.; Raptis, R.G.; Baran, P.; Salifoglou, A. pH-Specific synthesis, spectroscopic, and structural characterization of an assembly of species between Co(II) and N,N-bis(phosphonomethyl)glycine. Gaining insight into metal-ion phosphonate interactions in aqueous Co(II)-organophosphonate systems. Inorg. Chim. Acta 2007, 360, 638-648.
30. Macrae, C.F.; Sovago, I.; Cottrell, S.J.; Galek, P.T.A.; McCabe, P.; Pidcock, E.; Platings, M.; Shields, G.P.; Stevens, J.S.; Towler, M.; et al. Mercury 4.0: From visualization to analysis, design and prediction. J. Appl. Cryst. 2020, 53, 226-235.
31. Hyperchem 7.52 Release for Windows; HyperCube, Inc.: Gainesville, FL, USA, 2002.
32. Appel, C.; Ma, L. Concentration, pH, and Surface Charge Effects on Cadmium and Lead Sorption in Three Tropical Soils. J. Environ. Qual. 2002, 31, 581-589.
33. He, Z.; Ren, B.; Hursthouse, A.; Wang, Z. Efficient removal of Cd(II) using SiO2-Mg(OH)2 nanocomposites derived from sepiolite. Int. J. Environ. Res. Public Health 2020, 17, 2223.
34. Naeem, M.A.; Imran, M.; Amjad, M.; Abbas, G.; Tahir, M.; Murtaza, B.; Zakir, A.; Shahid, M.; Bulgariu, L.; Ahmad, I. Batch and column scale removal of cadmium from water using raw and acid activated wheat straw biochar. Water 2019, 11, 1438.
35. Klapiszewski, Ł.; Siwinska-Stefanska, K.; Kołodynska, D. Development of lignin based multifunctional hybrid materials for Cu(II) and Cd(II) removal from the aqueous system. Chem. Eng. J. 2017, 330, 518-530.
36. Kuchekar, S.R.; Patil, M.P.; Aher, H.R.; Gaikwad, V.B.; Han, S.H. Adsorptive removal of cadmium (II) ion from industrial wastewater by natural adsorbent. J. Mater. Environ. Sci. 2019, 10, 1117-1122.
37. Mahmoud, M.A.; El-Halwany, M.M. Adsorption of cadmium onto orange peels: Isotherms, kinetics, and thermodynamics. J. Chromatogr. Sep. Tech. 2014, 5, 1-6.
38. Essomba, J.S.; Ndi Nsami, J.; Belibi Belibi, P.D.; Tagne, G.M.; Ketcha Mbadcam, J. Adsorption of cadmium(II) ions from aqueous solution onto kaolinite and metakaolinite. Pure Appl. Chem. Sci. 2014, 2, 11-30.
39. Li, Z.; Katsumi, T.; Imaizumi, S.; Tang, X.; Inui, T. Cd(II) adsorption on various adsorbents obtained from charred biomaterials. J. Hazard. Mater. 2010, 183, 410-420.
40. Asci, Y.; Nurbas, M.; Acikel, Y.S. Sorption of Cd(II) onto kaolin as a soil component and desorption of Cd(II) from kaolin using rhamnolipid biosurfactant. J. Hazard. Mater. 2007, 139, 50-56.
41. Oh, S.; Kwak, M.Y.; Shin, W.S. Competitive sorption of lead and cadmium onto sediments. Chem. Eng. J. 2009, 152, 376-388.
42. Boparai, H.K.; Joseph, M.; O'Carroll, D.M. Kinetics and thermodynamics of cadmium ion removal by adsorption onto nano zerovalent iron particles. J. Hazard. Mater. 2011, 186, 458-465.
43. Leyva-Ramos, R.; Bernal-Jacome, L.A.; Acosta-Rodriguez, I. Adsorption of cadmium(II) from aqueous solution on natural and oxidized corncob. Sep. Purif. Technol. 2005, 45, 41-49.
44. Demadis, K.D.; Famelis, N.; Cabeza, A.; Aranda, M.A.; Colodrero, R.M.; Infantes-Molina, A. 2D corrugated magnesium carboxyphosphonate materials: Topotactic transformations and interlayer "decoration" with ammonia. Inorg. Chem. 2012, 51, 7889-7896.
45. Hamad, S.; Balestra, S.R.G.; Bueno-Pérez, R.; Calero, S.; Ruiz-Salvador, A.R. Atomic charges for modeling metal-organic frameworks: Why and How. J. Solid State Chem. 2015, 223, 144-151.
46. Zhu, H.; Yuan, J.; Tan, X.; Zhang, W.; Fang, M.; Wang, X. Efficient removal of Pb2+ by Tb-MOFs: Identifying the adsorption mechanism through experimental and theoretical investigations. Environ. Sci. Nano 2019, 6, 261-272.
47. Maranescu, B.; Visa, A.; Mracec, M.; Ilia, G.; Maranescu, V.; Simon, Z.; Mracec, M. Lamellar Co2+ vinylphosphonate metal organic framework. PM3 semi-empirical analysis of structural properties. Rev. Roum. Chim. 2011, 56, 473-482.
48. Burdett, J.K.; Coddens, B.A.; Kulkarni, G.V. Band gap stability of solids. Inorg. Chem. 1988, 27, 3259-3261.
49. Daniel, C.R.A.; Rodrigues, N.M.; da Costa, N.B., Jr.; Freire, R.O. Are Quantum Chemistry Semiempirical Methods Effective to Predict Solid State Structure and Adsorption in Metal Organic Frameworks? J. Phys. Chem. C 2015, 119, 23398-23406.
Aurelia Visa1,†, Bianca Maranescu1,†, Lavinia Lupa2,†, Luminita Crisan1,† and Ana Borota1,*,†
1“Coriolan Dragulescu” Institute of Chemistry, 24 M. Viteazul Ave, 300223 Timişoara, Romania
2Faculty of Industrial Chemistry and Enviromental Engineering, University Politehnica Timisoara, 2 Piata Victoriei, 300006 Timişoara, Romania
*Author to whom correspondence should be addressed.
†All authors contributed equally to this work.
© 2020. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Crisan, Luminita; Borota, Ana