Structural and Adsorptive Characteristics of 2D Multilayer Nanoflakes of NiCo Phosphates for Chromium(VI) Removal: Experimental and Monte Carlo Simulations

Metal phosphates are efficient adsorbent materials for heavy elements present in industrial effluents because of their promising properties. Hexachromium ions are among the most dangerous contaminants owing to their harmful properties and non-degradability. Accordingly, the present work offers a simplified study of the preparation of bimetallic phosphate materials from nickel cobalt phosphate (NiCo-Ph) based on the sol–gel method in an equimolar ratio. Characterization of the bulk, crystal phase, texture profile, and nanosize of NiCo-Ph was carried out using various techniques such as Fourier transform infrared spectroscopy, X-ray diffraction, X-ray photoelectron spectroscopy, nitrogen adsorption–desorption isotherm measurements, field emission scanning electron microscopy, transmission electron microscopy, and Raman spectroscopy. In this regard, the adsorption performance of NiCo-Ph was exemplified through six batch experiments, elucidating the impacts of the sorbent dose, initial concentration of pollutants, sorption time, temperature, pH, and shaking rate. According to UV/vis spectrophotometry measurements and their related calculations of NiCo-Ph, the maximum removal efficiency (RE %) of 92% and adsorption capacity (qm) of 37 mg/g were achieved at pH = 6, a dose of 5.0 g/L, 100 mg/L of [Cr(VI)], 300 rpm, adsorption time of 60 min, and 298 K. Monte Carlo simulations were also carried out to correlate the experimental data with theoretical calculations that provided a higher negative value (−911.62 kcal mol–1) for the adsorption energy of Cr(VI) in acidic medium. The adsorbent NiCo-Ph prepared by this direct method is therefore recommended for the quantification of Cr(VI) under slightly acidic solutions and at room temperature, which can maintain its efficiency even up to six cycles.


INTRODUCTION
One of the common toxic contaminants is heavy metals. They are not only non-biodegradable but also toxic due to their bioaccumulation in living cells. Among the documented toxic elements, Cr, As, Cd, Pb, Ni, Cu, Zn, and Hg 1 are not only dangerous for human health but also affect aquatic and environmental ecosystems. Consequently, it is necessary to purify the polluted waters from these heavy metals before being discharged into water bodies. Accordingly, several approaches were investigated for the removal of heavy metals including chemical precipitation, 2 electrodialysis, 3 electrocoagulation, 4 membrane filtration, 5 adsorption, 6 and photocatalysis. 7 Recently, adsorption has been a significant interest area with a special focus on recognizing low-cost adsorbents. 8−10 Chromium exists in the environment generally in hexavalent and trivalent forms. Chromium in the hexavalent form is very poisonous and is located in the first assembly of carcinogenic materials. Chromium is one of the most commonly used heavy metals in different industrial applications such as dyeing, electroplating, tannery, and fertilizer production. When the polluted water of these manufacturing industries enters the environment, the Cr metal spreads through the environment. 11 The U.S. Environmental Protection Agency has stated hexavalent chromium as one of the dangerous toxic contaminants, even in very low doses. Extreme chromium consumption by humans leads to renal and hepatic injury, gastrointestinal irritation, capillary damage, and central nervous system disorders. Chromium could lead to cancer of the digestive tract and lungs. 5 In wastewaters, Cr(VI) is predominantly present in the form of extremely soluble oxyanions, for instance, as bichromates, dichromates, and chromates. In order to decrease the severe harm caused by them, in particular, to ecology, it becomes the responsibility of these industries to diminish the Cr(VI) dose in their wastes to a suitable level prior to releasing them into the environment. Also, the recovery of chromium present in wastewaters is an attractive choice for economic reasons. 4 The utilization of the above diverse approaches for removal of chromium appears to be restricted owing to the low removal capacity, comparatively high costs, release of wastes formed by treatment, the need for expending a great energy, the need for fixed system monitoring, production of toxic sludges, and generally the need for a high capital investment and high operating costs. 5−11 Ultimately, the adsorption technique should be an efficient, facile, less-expensive, and eco-friendly approach. Additionally, it should be among the most operative and public-appropriate technology that is regularly utilized over other methods in the worldwide ecological protection zones. 12,13 Adsorption has numerous benefits over the orthodox approaches, which comprises sorbent reusability, low selectivity for precise metals, less operational cost, no chemical sludges, and a small operation time. Different scientists have investigated the efficiency of several other adsorbents for Cr(VI) removal from wastewater using fixed-bed and batch methods. Some of the adsorbents are Leucaena leucocephala seed shell-activated carbon, 14 activated carbon based on luffa sponge, 15 polyaniline doped with sulfuric acid, 16 magnetite nanoparticles, 17 and palm trunk charcoal. 18 Eventually, these sorbents are non-eco-friendly and hardly reflected in their costs. On the other hand, bimetal phosphates are an alternative that are naturally approachable and have been extensively explored as adsorbents. 19 The well-structured morphology of metal phosphates has much advantage and displayed considerable success in numerous applications as well as in technological advance. 20 They have been used in lithium batteries, 21 as catalysts for chemical conversion, 22 in electrocatalytic conversion processes, 23 and as supercapacitors. 24 Surprisingly, bimetal phosphate mesoporous adsorbents were scarcely reported in the literature. 20 In this report, novel nanomesoporous NiCo-Ph was synthesized with a good surface area via a facile and costeffective method. The introduced bimetal NiCo-Ph was characterized using field emission scanning electron microscopy (FE-SEM), transmission electron microscopy (TEM), Xray diffraction (XRD), Fourier transform infrared (FTIR) spectroscopy, Raman spectroscopy, and X-ray photoelectron spectroscopy (XPS), and the pore size and BET surface area were also examined. Meanwhile, the nanomesoporous NiCo-Ph was utilized in the removal of hexavalent chromium ions from aqueous media. The batch technique was utilized; parameters including the contact time, pH, and initial chromium dose were investigated. Furthermore, Monte Carlo (MC) simulations were performed to confirm the experimental adsorption findings.

EXPERIMENTAL SECTION
2.1. Synthesis of NiCo-Ph. All the chemicals utilized in this approach were of analytical grade and used directly without any additional refinement. All the reagents were obtained from Sigma-Aldrich and Chemical Reagent Co., Ltd. Figure 1 illustrates the sol−gel synthesis pathway to prepare the sample. In a typical blend of equimolar Ni and Co, precursors were hydrolyzed in a mixture of bidistilled water and EtOH (100:50). This is followed by adding a hexadecyl pyridinium chloride surfactant (HDPCl) with regular addition of urea to the previous blend under stirring for 15 min. Then, this mixture was dropwise transferred to an appropriate volume of ammonium monohydrogen phosphate (NH 4 ) 2 HPO 4 by adjusting the pH to 6−7, and this mixture was subjected to stirring for further 3 h at a constant temperature of 40°C. The resultant blend, which appeared as a sol, was sonicated for half an hour, then digested at room temperature to the next day, followed by filtration, washing, and drying at 90°C for a day. The attained sample was abbreviated as NiCo-Ph, and it was characterized using various techniques.
2.2. Material Characterization. FTIR investigation was carried out in the range 400−4000 cm −1 using a Bruker FTIR spectrometer. Characterization of the bulk and the crystal phase of the sample as powder was done using an X-ray diffractometer at 2θ from 10 to 80°. XPS data were obtained as reported previously. 25 The specific surface area (S BET ) and pore diameter distribution were also determined via studying the nitrogen adsorption−desorption isotherms (Micromeritics ASAP2010) at 77 K using the Brunauer−Emmett−Teller method. The morphology of the powder sample was explored using a field emission scanning electron microscope equipped with an energy-dispersive X-ray spectrometer (JSM-5410 model JEOL, Japan). The particle profile and nanosize were estimated by TEM using a Jeol-1230 electron microscope working at 200 keV. Raman spectroscopy was utilized to study the stability of the chemical structure of the prepared adsorbent using a Raman spectrometer (Horiba Scientific, SUA) before and after the adsorption of Cr(VI).
2.3. Batch Adsorption Experimentations. All experiments were performed in the 298−328 K temperature range. The adsorption tests were completed in batches. The experiments of batch adsorption were accomplished by shaking 0.5 g of NiCo-Ph with 100 mL of chromium(VI) solutions (5.0 g/L). The bottles were covered and stirred in a shaking water bath (THERMOLAB) at an agitated speed of 200−400 rpm. The concentration of chromium(VI) in the experimental medium was spectrophotometrically measured using a UV/vis spectrophotometer (model Shimadzu UV-1800 Japan) in the wavelength range from 180 to 800 nm and in a 1.0 cm quartz cell. After the adsorption experiment, the absorbance values were measured at λ max = 360.0 nm. All the tests were repeated three times, and the average of the results was recorded. The standard deviation (±SD) was calculated between 2.0 and 4.0% for all experiments. The influence of the NiCo-Ph sorbent mass (0.05−1.0 g/100 mL) on Cr(VI) removal was examined. The influence of the contact time on Cr(VI) adsorption was investigated from 0 to 200.0 min. The pH values of the tested experiments were adjusted from 1.0 to 8.0 using 0.1 M KOH or 0.01 M HNO 3 . The removal efficacy (RE/%) of Cr(VI) and the amount of adsorbed Cr(VI) Q e (mg/g) at equilibrium were calculated using the following equations i k j j j j j y where W (g) represents the sorbent weight, V (L) symbolizes the solution volume, And C e and C i (mg/L) are the equilibrium and initial Cr(VI) concentrations, respectively.  26 The adsorption of Cr(VI) ions in acidic (hydronium ions + water molecules), basic (hydroxide ions + water molecules), and neutral (water molecules only) media and on the surface of NiCo-Ph was achieved in a simulation box (32.27 Å × 32.27 Å × 50.18 Å) with periodic boundary tuning. The Forcite classical simulation engine was employed to improve the energy of Cr(VI) ions. 27 Also, the specifics of computational calculations were illustrated in our former published research. 27,28 3. RESULTS AND DISCUSSION 3.1. Material Characterization. Figure 2A depicts the FTIR spectrum of the NiCo-Ph sample. The bands in the region of 3440−3030 cm −1 show the vibrational mode of crystalline water. A wideband at 3440−3250 cm −1 is due to the interference of OH groups in adsorbed water and the OH groups in ammonium ions. 29 Correspondingly, the bands of 1580 and 1450 cm −1 are ascribed to the bending mode of H 2 O. A group of small bands lying in the region of 1030−970 and 860 cm −1 is assigned to the vibrational mode of the anion of PO 4 3− as symmetric and asymmetric O−P−O, respectively. In the low-frequency region, a sharp band of bending of P−O is detected at 540 cm −1 . Figure 2B illustrates the XRD pattern of the studied sample, pointing out the effective sol−gel route in synthesizing the bimetal Ni and Co phosphates. The pattern shows piercing diffraction peaks at 2θ = 11.48, 13.40, 18.47, 23.30, 28.23, 30.36, and 33.35°. Furthermore, the predominant (hkl) planes of (110), (020), (200), (101), (031), (3̅ 01), and (3̅ 21), respectively, are consistent with the described phase of the monoclinic hydrated nickel phosphate and monoclinic hydrated cobalt phosphate (XRD cards 00-033-0951 and 00-041-0375). 30 The grain size was assessed to be 66 nm using the full width at half-maximum (fwhm) value and the Debye− Scherrer equation. This result is in good agreement with the particle size, as estimated by SEM studies.
To clarify the distribution of the hydrodynamic size of the NiCo-Ph sample, Figure 2C displays its dynamic light scattering (DLS) distribution in the investigated medium of Cr(VI). DLS measurement displays a broad size distribution of 333−2350 nm and is maximized at 1100 nm. On the other hand, this hydrodynamic size is greater than the particle size distribution illustrated by SEM and TEM methods, which is due to the electric double layer on the particle interface.
The nitrogen adsorption/desorption isotherm plot ( Figure  2D) displays the isotherm type, which is indexed to the type IV isotherm and characterized to be mesoporous nanomaterials. The specific surface area S BET and pore width distributions were estimated to be 8 m 2 /g and 8.4 nm, respectively. This good value in terms of pore width distribution that is in the mesoporous range is the reason behind the efficiency of the NiCo-Ph sample in Cr(VI) removal. This in turn is needed to achieve a reasonable Cr(VI) adsorption capacity in a quite short contact time.
XPS investigation of the NiCo-Ph sample was carried out to show the valence or electronic surface states of the contents. The XPS spectra ( Figure 3A) and their data ( Table 1) show good agreement with the XRD data. Figure 3A shows all regions of the XPS spectrum that is characterized for Ni 2p (859.8 eV), Co 2p (784.1 eV), O 1s (530.6 eV), and P 2p (136.3 eV). The short XPS spectrum points to the oxidation states of NiCo-Ph as Ni (II), PO 4 3− , and Co 3+ . Minutely, the high-resolution XPS spectrum ( Figure 3B) in the range of binding energy of Ni 2p displays typical peaks at 879.8 and 856.4 eV as well as satellite peaks, which are attributed to Ni 2p 1/2 and Ni 2p 3/2 , respectively, 31 that confirmed the presence of Ni as Ni 2+ in the prepared NiCo-Ph sample. Correspondingly, peaks at 798 and 782 eV ( Figure 3C) are attributed to Co 3+ and Co 2+ ions and are distinguished as the binding energy peaks of Co 2p 1/2 and Co 2p 3/2 , respectively. 32 The two main peaks go along with two satellite peaks centered at 788.3 and 784.88 eV due to the high spin of the Co 2+ ion. The peak at 782 eV was due to the Co oxidation state being linked with Co−PO x , which is probably caused by the surface oxidation. 33−35 The O 1s spectrum is shown in Figure 3D, the deconvoluted or fitted O 1s peaks at 530.1, 531.7, and 532.6 eV are attributed to M−O or hydroxyl groups. Figure 3E represents a peak at 134.2 eV, which is attributed to P−O of (PO 4 ). 36 SEM micrographs at different magnifications are depicted in Figure 4A,B. As is clearly shown, the examined sample showed a nanoflake morphology in two dimensions, which accumulated above each other forming a multilayer with a nanoflower shape. The lengths of nanoflakes range from 100 to 200 nm. Analysis of SEM results was done with ImageJ software. The histogram of nanoparticles ( Figure 4F) shows a reasonable size distribution from 3 to 40 nm with a mean size distribution of 20 ± 2 nm.
A further study of the morphology and topology of the NiCo-Ph sample was also carried out using TEM ( Figure  4C,D). Typically, the profile of the sample, where the nanoflakes appear on one side and full nanoflakes mixed with semispherical particles with diameters ranging from 20 to 80 nm. These results are in full agreement with those obtained using SEM.
The corresponding SAED image is presented in Figure 4E, which illustrates all rings of different planes with a lattice gap (d-spacing) of 6.5 Å matching with the high relative intensity (100%) of the (020) plane of Co 2p and Ni 2p, which is proportionate to the XRD findings.
3.2. Adsorption Study. Toward estimating the ability of the prepared bimetal phosphate NiCo-Ph as a potential adsorbent for Cr(VI) ions from aqueous phases, six analogous control experiments were carried out to study and identify the optimum parameters required for the adsorption of chromium ions using the adsorbent prepared.
At first, Figure 5 depicts the variations between the absorption spectra at time intervals ranging from 0 to 100 min for 100 mL sample of Cr(VI) ions (100 mg/L) containing 0.5 g/100 mL of NiCo-Ph at 298 K; there are two distinct peaks, an obvious sharp peak at λ = 360 nm and a slight broad peak at 425 nm. As noticeable, there is a diminution in the intensity of spectral curves with increasing time, where the sharp peak at 360 nm slightly faded after 100 min of contact time. The difference in intensities of the UV−vis spectra confirms the superior performance of the NiCo-Ph adsorbent, which is due to the low concentrations of Cr(VI) species. Therefore, Figure 6A shows the effect of the sorbent dose within the range of 0.1−1 g/100 mL in a control experiment at pH = 6, with an initial concentration of Cr(VI) of 100 mg/L, a T of 298 K, and a contact time of 60 min. The removal efficiency (%) distinctly increases with the increase of the adsorbent weight and reaches 92% by consuming 5.0 g/L of the examined sorbent. After that, the upturn is slight and does not clearly intrude on the efficiency of the adsorbent. The enhancement of the removal efficiency by increasing the adsorbent dose can be attributed to the availability of the free active sites on the surface of our adsorbent. The binding may be due to the electrostatic attraction between the charged positive surfaces of NiCo-Ph and the anionic chemistry of solutions of chromium ions in an adsorption process at pH = 6. The optimum mass of the sorbent was identified to be 5.0 g/L of Cr(VI) ions (100 mg/L).
Another control experiment disclosed the sufficient contact time for the adsorption of Cr(VI) ions on the surface of the bimetallic phosphate sample ( Figure 6B). This control experiment was carried out at pH = 6, with the adsorbent dose of 5.0 g/L, and [Cr(VI)] of 100 mg/L at a contact time interval of 15−200 min. This behavior facilitated enhancing the removal efficiency of up to 92% at 60 min accompanied by a flat plateau even after undergoing further sorption. This asserted the consumption of free active sites accessible for the Cr(VI) sorption. The initial Cr(VI) dose is shown in Figure  6C for another control experiment using an adsorbent dose of 5.0 g/L, pH = 6, a sorption time of 60 min, and a concentration range of 20−250 mg/L. As detectable from the Cr(VI) dose graph, there is a slight change observed in the Cr(VI) sorption on the NiCo-Ph surface upon increasing the Cr(VI) dose up to 250 mg/L, especially after 100 mg/L, which corresponds to the removal of 92%. This reasonable removal % value of 92% is sufficiently obtainable for a high Cr(VI) dose of 100 mg/L. Furthermore, the removal % decreased from 92 to 75% with increasing Cr(VI) concentrations, which caused the active centers of the used adsorbent to become saturated and clogged. Figure 6D exemplifies the influence of the shaking speed on Cr(VI) sorption, which was scrutinized under a controlled experiment at pH = 6, 25°C, 0.5 g NiCo-Ph/100 mL, an adsorption time of 60 min, and [Cr(VI)] = 100 mg/L, within speed limits of 200−400 rpm. To clarify the role of the shaking rate, the shaking rate was increased, and a drastic change in the removal % value was noted from 75% at 200 rpm until reaching 97.5% at 300 rpm. However, increase of the shaking rate to higher than 300 rpm has no significant change in the removal rate, declaring the adequacy of this speed rate (300 rpm) for the outer transfer from the bulk solution of Cr(VI) ions to the frontier surface of the used adsorbent. Previous analogous effects reported the adsorption of Cr(VI) solutions. 37,38 Another control experiment was also conducted for studying the impact of alteration of temperature. The test was conducted in the temperature range of 298, 308, 318, and 328 K, keeping the other control parameters of pH = 6, 0.5 g NiCo-Ph/100 mL, adsorption time of 60 min, and [Cr(VI)] = 100 mg/L the same ( Figure 7A). Unexpectedly, the increase in temperature causes a slight diminution in the rate of Cr(VI) removal (%) from 97% at 298 K to 85% at 328 K. This outcome may be elaborated in terms of the nature of the pores of the investigated adsorbent. In addition, a competitive desorption process has occurred instead of the adsorption process due to the high thermal energy. This distinguishes our absorbent NiCo-Ph for its efficiency under ambient conditions. A comparable approach also illuminated this effect for the adsorptive removal of Cr(VI) from an aqueous solution. 39 Figure 7B depicts the effect of varying pH from 4 to 9on the removal rate %, with the same control parameters of 0.5 g NiCo-Ph/100 mL, an adsorption time of 60 min, 25°C, and [Cr(VI)] = 100 mg/L. As is clearly shown, an increase in the removal rate % to 92% has been observed upon increasing the pH value to reach a maximum of 6.0. Then, a decline in the removal % value was observed at pH values >6, and the removal % was as low as 30% at pH = 9.
From these insights, it can be said that the adsorption of hexavalent chromium ions is fundamentally contingent on the pH value of a solution. This can be due to the active groups on the adsorbent surface as well as the structure or valence of chromium ions present in the solution. Cr(VI) removal using NiCo-Ph as an adsorbent is preferred at slightly acidic pH. Previous contributions reported that at acidic pH values, the metal ions in bimetal phosphates exist as positive ions until pH 6, which are then converted to negatively charged metal phosphate ions at pH greater than 7. 40 Therefore, as described by XPS outcomes, bimetal phosphates (NiCo-Ph) are present in Ni 2+ and Co 3+ species. Besides, the prevailing ionic forms of Cr ions in the aqueous solution are HCr 2 O 7 − or Cr 2 O 7 2− at slightly acidic pH. 41 Under these conditions, an electrostatic attraction arises between the anionic species of chromium and the positively charged surface of NiCo-Ph through an anionic adsorption process, which is prominently improved. While with the increment of pH values, the surface layer of the adsorbent NiCo-Ph became negatively charged as MPO 4 − . In addition, CrO 4 2− and OH − ions were initiated in an intense rivalry, hence a mutual repulsion takes place resulting in poor sorption of Cr(VI). Moreover, this effect can be elucidated in terms of pH-response characteristics, which enables the examined sorbent to behave as a hydrophilic adsorbent under acidic conditions and as an amphiphilic adsorbent under slightly neutral conditions. This can be attributed to the  use of HDPCl as a cationic surfactant in the preparation route of the synthesized adsorbent. The amine surfactants have higher amphiphilic property at low pH values once the NH 2 group is protonated. Consequently, these amine groups in the used surfactant likely act as pH-responsive amphiphiles. 42 Thus, they can competently enhance the degradation of adsorbed Cr(VI) ions with outstanding recyclability; a similar previous study was reported by Shi et al. 43 It can also be interpreted in terms of adsorption energy and its correlation with the removal ability of Cr(VI), as will be elaborated in the following section of MC Simulations. Figure 8A discloses the regeneration of the NiCo-Ph sorbent through seven repetitive cycles under optimum conditions, which were selected as illustrated previously (100 mL of 100 mg/L Cr(VI), 5.0 g/L sorbent, pH = 6, and 25°C). The probability of the used adsorbent due to its activity preserved the outstanding adsorption behavior over several cycles. The outcomes revealed a minor lessening in the removal efficiency of hexachromium ions from 96% up to ∼80% after seven cycles. Figure 8B shows the Raman spectra of NiCo-Ph before and after 60 min of Cr(VI) removal. This is to approve the establishment of NiCo-Ph and its stabilization even after 60 min as optimum adsorption time for Cr(VI) removal. A pronounced sharp peak was observed at 560 cm −1 , which is attributed to ν antisymmetric bending vibrations of the phosphate bond P−O−P. 44 A slight shift of the peak at 780 cm −1 (before the adsorption process) at a lower wavenumber was noted that may be pronounced for phosphate network modifiers or a system of transition metals doped with P 2 O 5 . 45 Moreover, this lower shift was also observed at 865 cm −1 for the NiCo-Ph sample after the adsorption process at 900 cm −1 , which is ascribed to the vibrational stretching terminal O atom containing the PO 4 group, 44 46 Furthermore, the sample of NiCo-Ph after association with chromium(VI) disclosed a typical peak in the range of 630 cm −1 which is associated with the existence of Cr(III), declaring the efficiency of the prepared NiCo-Ph in the reduction of Cr(VI) to Cr(III) ions and thereby its potential for remediation of water from Cr(VI) ions. 47,48 3.3. Adsorption Isotherm Considerations. In order to estimate the adsorption efficiency of NiCo-Ph nanoflakes, chromium(VI) adsorption onto 0.5 g nano-NiCo-Ph/100 mL Cr(VI) at equilibrium is investigated using the models of Langmuir and Freundlich isotherms. 49 The Freundlich model was commonly used to describe the inhomogeneity of the surface. The model of Langmuir assumes that adsorption takes place on the adsorbent's external interface. The Langmuir where C e is the Cr(VI) equilibrium concentration in medium (mg/L), n is the Freundlich exponent interrelated to adsorption power, K f represents the Freundlich constant (mg/g), 1/n characterizes the adsorption density, q m represents the maximum adsorption efficacy (mg/g), K L is the Langmuir constant (L/mg), and q e symbolizes the adsorption effectiveness at equilibrium (mg/g). The aboverevealed factors can be obtained from the relation of C e /q e versus C e and log q e versus log C e plots. The value of n was found to be 2.21 (greater than 1), which designates that the adsorption process has promising ability. The linear plots of Langmuir and Freundlich isotherm models for Cr(VI) at pH = 6 and a sorbent weight of 0.5 g/100 mL Cr(VI) are presented in Figure 9A,B. The computed parameters of the investigated adsorption models are recorded in Table 2. Based on the correlation coefficient values (R 2 ), the data fit well to the model of Langmuir (R 2 = 0.986) than to the Freundlich model (R 2 = 0.965). The characteristics and feasibility of the Langmuir model according to the factor of dimensional separation, R L , is given using the following formula 51 where C i is the initial concentration of Cr(VI). The values of R L confirm if the adsorption is promising (0 < R L < 1), irreversible (R L = 0), linear (R L = 1), or unfavorable (R L > 1). The R L value for the investigated sorbent NiCo-Ph nanoflakes was found to be 0.023 (in the range of 0 < R L < 1). This value specifies a suitable adsorption isotherm type of Cr(VI) ions onto the NiCo-Ph surface in the examined concentration range. The n value changes with the sorbent heterogeneity and indicates the favorability of the adsorption route. The n value should be 10 < n > 1 under promising adsorption conditions. The value of n in the present report was found to be 2.20, which displays that the Cr(VI) sorption onto NiCo-Ph nanoflakes is favorable.

MC Simulations.
MC simulations were carried out to study the adsorption between the Cr(VI) ions and the NiCo-Ph surface. Figure 10 therefore shows the highest proper adsorption arrangements for the Cr(VI) ions on the NiCo-Ph surface achieved using the adsorption locator module. 52 Furthermore, the attained outcomes from the MC simulations are shown in Table 3 such as the adsorption energy for relaxed adsorbate molecules, rigid adsorption energy for unrelaxed adsorbate molecules, and deformation energy for relaxed adsorbate molecules. 53 Table 3 divulges that Cr(VI) (−911.62 kcal mol −1 ) has a higher negative value of adsorption energy in acidic medium in comparison with the value of Cr (VI) in basic or neutral medium, which confirms the powerful adsorption of Cr(VI) ions on the NiCo-Ph surface in acidic medium. Additionally, the values listed in Table 3 manifest that the adsorption energy values of Cr(VI) in acidic medium preand postgeometry optimization steps (i.e., unrelaxed and relaxed values of −822.88 and −88.7379 kcal mol −1 , respectively) are further negative than those in basic or neutral medium, asserting the higher adsorption of Cr(VI) ions on the NiCo-Ph surface in acidic medium.
The adsorption energy of an adsorbate when one of the adsorbates has been omitted is elucidated through the dE ads / dN i values. 54 The dE ads /dN i values for Cr(VI) in acidic medium (−162.59 kcal mol −1 ) are higher than those in basic or neutral medium, which declares the excellent adsorption of Cr(VI) ions on the NiCo-Ph surface in acidic medium. Moreover, the dE ads /dN i values of water molecules, hydronium ions, and hydroxide ions are around −13.87, −32.11, and −25.42 kcal mol −1 , respectively, which are little compared to the values of Cr(VI) ions on the NiCo-Ph surface in different media, revealing the stronger adsorption of Cr(VI) ions than other ions, which supports the exchanging of water molecules, hydronium ions, and hydroxide ions for Cr(VI) ions. Therefore, the Cr(VI) ions are decisively adsorbed on the NiCo-Ph surface; this is affirmed by experimental and theoretical studies collectively.   Table 4 displays the findings of a comparison of the adsorption capacity and pH of solution for Cr(VI) removal for some reported adsorbents. 55−62 The adsorption capacity, q m (mg/g), was found to be 37.0 mg/g for the fabricated NiCo-Ph. This outcome indicates that this material displays a very significant adsorption capability compared to other kinds of composites and materials, for instance, nano-Al 2 O 3 has a q m of 17.7 mg/g 55 and hydrous stannic oxide has a q m of 17.7 mg/g, 56 as recorded in Table 4.

CONCLUSIONS
In the present study, the structural and adsorptive characteristics of 2D multilayer nanoflakes of NiCo-Ph for chromium-(VI) removal were thoroughly scrutinized. The bimetal NiCo-Ph was simply manufactured by a sol−gel route in an aqueous solution of HDPCl with the regular addition of a ureastabilizing agent. In addition, the synthesized NiCo-Ph was  fully characterized via various methods including XRD, FTIR, DLS, XPS, SEM, TEM, and studying the nitrogen adsorption− desorption isotherms. SEM and TEM characterizations asserted the self-assembly of the nanoflake structure resulting in the formation of NiCo-Ph flowers with diameters ranging from 20 to 80 nm. XRD and SAED outcomes emphasized the formation of the (020) plane of Co 2p and Ni 2p with a lattice gap (d-spacing) of 6.5 Å related to Ni 3 (PO 4 ) 2 ·xH 2 O and Co 3 (PO 4 ) 2 ·xH 2 O, respectively. Besides, XPS results proclaim the following oxidation states for NiCo-Ph: Ni 2+ , PO 4 3− , and Co 3+ . To correlate NiCo-Ph and its adsorbtion potential for the Cr(VI) removal, six control experiments were carried out to identify the optimum parameters. Batch adsorption tests were carried out at room temperature, pH = 6, a sorption time of 60 min, 0.5 g of sorbent per 100 mL of Cr(VI) ions (100 mg/L), and a rate of 300 rpm, which afforded a removal ratio of 92% and an adsorption capacity (q m ) of 37 mg/g. Significantly, NiCo-Ph displayed pronounced adsorption of Cr(VI) throughout the regeneration tests by maintaining ∼80% of its intrinsic efficacy even after the 7th adsorption cycle. Raman spectroscopy analysis established the stability of the synthesized adsorbent NiCo-Ph and further asserted its effectiveness for the remediation of water from Cr(VI) ions and their perfect reduction to the Cr(III) state. Likewise, MC simulations also affirmed the facile action of our adsorbent based on the highest proper adsorption configurations for the Cr(VI) ions on the NiCo-Ph surface. Furthermore, the higher negative adsorption energy values of Cr(VI) in the acidic medium of unrelaxed and relaxed adsorbate molecules of −822.88 and −88.7379 kcal mol −1 , respectively, emphasize the reasonable adsorption of Cr(VI) ions in acidic medium. Given the above, the bimetal phosphate NiCo-Ph nanoflakes are considered among the most promising adsorbents for Cr(VI) species uptake in a slightly acidic medium and at ambient temperature.