Humulus scandens-Derived Biochars for the Effective Removal of Heavy Metal Ions: Isotherm/Kinetic Study, Column Adsorption and Mechanism Investigation

Humulus scandens was first adopted as a biomass precursor to prepare biochars by means of a facile molten salt method. The optimized biochar exhibits a high specific surface area of ~450 m2/g, a rich porous structure and abundant oxygen functional groups, which demonstrate excellent adsorption performance for heavy metal ions. The isotherm curves fit well with the Langmuir models, indicating that the process is governed by the chemical adsorption, and that the maximum adsorption capacity can reach 748 and 221 mg/g for Pb2+ and Cu2+, respectively. In addition, the optimized biochar demonstrates good anti-interference ability and outstanding removal efficiency for Cu2+ and Pb2+ in simulated wastewater. The mechanism investigation and DFT calculation suggest that the oxygen functional groups play dominant roles in the adsorption process by enhancing the binding energy towards the heavy metal ions. Meanwhile, ion exchange also serves as the main reason for the effective removal.


Introduction
With the rapid development of China's economy, a large volume of wastewater containing heavy metal ions is being discharged by the mining, metallurgy, chemical, electronics and other production industries [1,2]. Heavy metal ion pollution remaining in different water sources seriously threatens people's safety and health [3,4]. Harmful heavy metal ions such as cadmium, chromium, lead, arsenic, copper, mercury and so on are commonly detected in different water bodies [5,6]. At present, chemical, physical and biological methods including electrochemical treatment, ion exchange, chemical precipitation and adsorption are the mostly commonly employed strategies to remove heavy metal ions from wastewater [7]. Among these methods, adsorption is considered an effective and economic method because of its facile operation and low cost [8,9]. At present, biomass-derived activated carbons with a large specific surface area and resource abundance have been widely used in wastewater treatment [10][11][12]. A large number of studies have proven that the outstanding adsorption performance of biochar originates from the large surface area, stable porous structure, abundant surface oxygen functional groups such as carboxyl and hydroxyl groups, and rich mineral species which can remove the harmful ions through the ion exchange process [13]. Therefore, it is of great significance to develop high-efficiency and low-cost biochar materials to remove heavy metals from polluted water environments [14].
Hitherto, many biomasses, such as agricultural wastes, woody materials, algae and so on [15][16][17][18], are widely employed to prepare biochar for the removal of heavy metal ions from the water system. For example, Undaria was adopted as a raw material to prepare biochar through rapid pyrolysis and physical steam activation. The as-obtained biochar was able to effectively remove Cu 2+ from the aqueous solution [19]. In order to achieve the fast removal of Pb 2+ , biomasses such as palm, soybean meal, shells and so on were also adopted to prepare biochars with more adsorption sites as well as rich ion diffusion channels [20][21][22][23][24]. Based on these studies, it was found that although they are prepared through similar treatment methods, the adsorption performance of biochars varies widely, which implies that the unique tissue structures as well as the composition of the biomass exerts great influence on the adsorption performance of the biochar. Therefore, it is of great significance to explore the family of biomass for high-performance biochars.
In addition, in order to improve the adsorption capacity of biochars, various approaches have been developed to achieve the chemical modification of amino, alkali, sulfur and phosphoric function groups on the biochar surface. For example, Zhang et al. modified the surface of a rice straw-derived biochar with an amino group through nitrification and amination to prepare high-efficiency Pb 2+ adsorbent. The as-introduced amino group shows strong complexation with Pb 2+ to improve the adsorption strength of biochars [23]. Ming et al. adopted EDTA to modify the activated carbon derived from shells, and the mechanism studies indicate that the electrostatic action and complexation of surface active groups play important roles in enhancing the Pb 2+ adsorption capacity [24]. Zhu et al. prepared surface-oxidized biochar from porous biochar with a high Cd 2+ -adsorption capacity, and they found that the adsorption capacity is mainly affected by the surface functional groups [25]. In addition to the oxygen and amino function groups, Cl active sites and sulfur groups were also deliberately introduced on the surface of the biochars to enhance the removal efficiency of Hg + and Cu 2+ , taking advantage of the enhanced binding energy between the active sites and the target heavy metal ions [26,27]. Although many schemes have been developed to achieve the surface modification of biochars, these strategies still suffer from many drawbacks, such as a tedious preparation procedure, high cost, toxic modification reagents and potential secondary contamination risk. In addition to the modification, the current production of biochars usually requires inert gas protection and high-temperature treatment with corrosive activators, which often have many disadvantages such as a high cost and a complicated preparation process. Therefore, it is necessary to develop a facile preparation process to produce biochar with rich surface function groups through a one-step pyrolysis process.
The molten salt method represents a commonly adopted route to prepare porous carbon materials with inert salt as the sealing and activating bi-functional medium. The biochars prepared through the molten salt method are widely used in the fields of environmental adsorption, energy storage devices, air purification and so on [28]. Due to the etching effect of molten salt and the template effect of salt crystals, mesopores and macropores can be introduced into the products [29]. At the same time, oxygen in the atmosphere can further react with highly active carbon atoms, which could introduce a large number of micropores into the products [30]. Therefore, the biochar prepared by the molten salt method has a high specific surface area and rich oxygen functional groups, which is conducive to the adsorption performance [31].
Humulus scandens (HS) is a perennial climbing herb which is widely distributed in the northern region of China. However, the massive growth of HS causes negative impacts on agricultural production and human activities. Moreover, the tissues of HS have regular porous structure and rich cellulose, which allows it to be promising precursor to prepare porous carbon materials. In this study, HS was firstly used as a precursor to prepare biochar by means of the NaOH modified molten salt method. The removal efficiency of Cu 2+ and Pb 2+ ions from an aqueous solution was greatly enhanced due to the welldeveloped porous structure, rich oxygen function groups as well as the embedded mineral species. The maximum Q m can reach 748 and 221 mg/g for Pb 2+ and Cu 2+ , respectively.
In addition, the effects of adsorption conditions including initial pH, contact time, and initial concentration on the adsorption capacities were systematically investigated. The adsorption mechanism was studied, which indicates that oxygen function groups and the imbedded mineral species play dominant roles in the adsorption process by enhancing the binding energy towards the heavy metal ions and facilitating the ion exchange process. In addition, the optimized biochar shows outstanding anti-interference ability and good column adsorption performance, permitting its potential application in industrial-scale utilization. Considering that HS is a harmful biomass for agricultural production, the effective conversion of HS into biochar with excellent adsorption capacity can not only enrich the biomass family to produce carbon materials, but also reduce the adverse effects of the massive growth of HS in the north of China.  Figure S1), Henan Province, China, was washed and dried in an oven at 50 • C for 12 h. The dried and pulverized HS was screened through a 100-mesh sieve for later use. Pb(NO 3 ) 2 , Cu(NO 3 ) 2 ·5H 2 O, NaOH, NaCl and KCl were purchased from Sinopharm Chemical Reagent Co. Ltd. (Beijing, China), were of analytic grade and were used without further purification.

Preparation of the Biochar
NaCl and KCl were mixed in an agate mortar at a mass ratio of 1:1. Then, 3 g of HS powder was mixed with 6 g of mixed salt and placed in a 50 mL ceramic crucible, covered with about 1 cm of mixed salt. The mixture was heated to 400, 500, 600, 700 and 800 • C in a muffle furnace at a rate of 10 • C/min, respectively. After continuous pyrolysis for 3 h, it was naturally cooled and cleaned with DI water to neutral. After drying in a 60 • C vacuum oven for 12 h, a series of biochars were obtained, and the biochar material obtained at a pyrolysis temperature of 600 • C was named HSC-M.
Then, 1, 2 and 3 g of NaOH were added to the mixture of 3 g HS powder and 6 g salt. After grinding evenly, the mixture was transferred into a ceramic pot, and covered with salt mixture with a thickness of about 1 cm. It was heated in a muffle oven at a rate of 10 • C/min to 600 • C and maintained for 3 h. The as-prepared biochars were named HSC-MA-1, HSC-MA-2 and HSC-MA-3 according to the NaOH dosage.
In addition, the washed, dried and crushed HS powder was placed in a ceramic crucible and placed in a tubular furnace. Under N 2 protection, the powder was sintered for 3 h at 600 • C with a heating rate of 10 • C/min. After having cooled naturally, the biochar was demoted as HSC-N 2 .

Characterization of the Biochar
The pyrolysis process of the biomass in the molten salt medium was studied using TG-DSC in the air. The surface morphology and element distribution of biochar were analyzed by means of field emission scanning electron microscopy equipped with an energy dispersion spectrometer (SEM, Hitachi SU8010, Tokyo, Japan). The internal structure of the samples was observed using transmission electron microscopy (TEM, TF20, JEOL 2100F, Tokyo, Japan). X-ray photoelectron spectroscopy (XPS, Escalab250XI, Thermo Fisher Scientific, Waltham, MA, USA) was used to record the surface elemental composition and chemical bonds of the samples. X-ray powder diffraction (XRD, Rigaku DMAX-2000, Rigaku, Tokyo, Japan) was used to demonstrate the lattice space structure of the sample. The contents of Ca, Mg, Fe and Cu in the samples were analyzed by inductively coupled plasma mass spectrometry (ICP-MS, ELAN DRC-e, PerkinElmer Inc, LAS, USA). A Fourier transform infrared spectrometer (FTIR, Nexus 470, GMI Inc., Lebanon, OH, USA) and Raman spectrometer (Labramhr EVO, Horiba France SAS, Palaiseau, France) were used to analyze the surface functional groups of the biochar in the wave number ranges of 4000~400 cm −1 and 3000~50 cm −1 , respectively. The specific surface area and pore diameter distribution of the samples were determined by means of the N 2 adsorption method using a surface area analyzer (BET, MicrotracBEL, Osaka, Japan).

Adsorption Experiment
Cu(NO 3 ) 2 ·5H 2 O and Pb(NO 3 ) 2 were added into two 500 mL volumetric flasks, respectively, and were prepared with 1000 mg/L Cu 2+ and Pb 2+ reserve solutions. Ten milligrams of the as-prepared samples were added into a 250 mL conical flask containing 100 mL, 100 mg/L Pb 2+ or Cu 2+ , respectively. The suspension was oscillated in a constant temperature water bath oscillator at 180 r/min for 4 h at 25 • C. After the reaction, the suspension was filtered using a 0.45 µm cellulose acetate membrane. Finally, the concentration of Pb 2+ and Cu 2+ after adsorption equilibrium was determined using a flame atomic absorption spectrometer.

Influence of Initial pH Value
The initial pH value of 100 mg/L Cu 2+ and Pb 2+ solution was adjusted to 2, 3, 4, and 5 with 0.1 mol/L HNO 3 or 0.1 mol/L NaOH solution. The effect of pH value on the adsorption of Pb 2+ and Cu 2+ by biochar in 100 mL of the Pb 2+ and Cu 2+ solution was determined by using the above method with a HSC-MA-2 dosage of 10 mg and a concentration of 100 mg/L.

Adsorption Kinetics and Isothermal
For the kinetics tests, the initial concentration of Pb 2+ was set as 200 mg/L and the biochar amount was set as 0.2 g/L. Meanwhile, the initial solubility of Cu 2+ was 100 mg/L, and the amount of biochar was 0.5 g/L. The contact time between the biochars and the metal ion solution was 5-300 min, and the adsorption capacities of Pb 2+ and Cu 2+ on biochar were measured at 5, 10, 20, 30, 60, 120, 180, 240 and 300 min, respectively.
For the isothermal tests, the initial concentration of Pb 2+ was set as 20, 50, 100, 150, 200, 250, and 300 mg/L, and the addition amount of biochar was 0.2 g L −1 , and the initial concentration of Cu 2+ was 20, 40, 80, 120, 160, 200, 240 mg/L, and the addition amount of biochar was 0.5 g/L. At a constant temperature of 25 • C, the water bath was oscillated at a speed of 180 RPM for 4 h, and the samples were sampled at adsorption equilibrium. After filtration, the concentrations of Pb 2+ and Cu 2+ in the filtrate were determined, and the adsorption capacities of the biochar HSC-MA-2 on Pb 2+ and Cu 2+ at corresponding initial concentrations were calculated.
The amount of heavy metal ions adsorbed per unit mass of the adsorbent (Q e ) was calculated from Equation (1): The formula for calculating the removal rate of metal ions is: c 0 is the initial concentration of the Cu 2+ and Pb 2+ ions (mg L −1 ); c e is the concentration of the Cu 2+ and Pb 2+ ions at adsorption equilibrium (mg L −1 ); Q e is the amount of heavy metal ions adsorbed per unit mass of the adsorbent (mg g −1 ); V is the total volume of the solution (L); and W is the mass of the biochar (g).
In order to study adsorption mechanism and kinetic parameters, quasi-first order, quasi-second order and intra-molecular diffusion models were employed to fit the adsorption curves, respectively.
Quasi first order dynamic equation: Quasi second order dynamic equation: The equation of the intramedullary diffusion model: In the equation, q e and q t are the adsorption capacity at adsorption equilibrium and the adsorption capacity at time t, respectively (mg/g); t is the adsorption time (min); k 1 is the primary kinetic adsorption rate constant (min −1 ), and k 2 is the secondary kinetic adsorption constant (g (mg/min) −1 ). The correlation coefficient obtained by the fitting calculation was used to determine the adsorption process.
Adsorption isotherms of adsorbents describe the relationship between solid and liquid phases. Langmuir and Freundlich models are two commonly used adsorption isotherm models. Langmuir and Freundlich adsorption isotherm models are defined as follows: Freundlich equation: ln q e = 1 n ln C + ln k f Langmuir equation: In the equation, C is the adsorption concentration (mg/L), q e is the equilibrium adsorption capacity (mg/g), q m is the maximum adsorption capacity (mg/g), k l is the Langmuir constant (L/mg), k f is the characteristic constant proportional to the equilibrium constant (L/mg), and n is the characteristic constant representing the properties of adsorption force.

The Effect of Background Ions
The effect of background ions on the adsorption of Pb 2+ and Cu 2+ ions by biochar was studied; 10 and 20 mg of HSC-MA-2 were added to 100 mL of 100 mg/L Pb 2+ and 100 mg/L Cu 2+ containing sodium, magnesium and calcium ions (with concentration of 1, 3, 5, 7 and 10 mmol/L), respectively, and then oscillated for 4 h at pH 5.0 and 25 • C.

Multi-Metal Adsorption
The interaction of coexisting heavy metals in the adsorption process has more practical significance than the interaction of a single system. A mixture of Pb 2+ , Cu 2+ , and Cr 3+ was prepared at pH 5.0, and the initial concentrations of all metal ions were set as 50, 100, 150 and 200 mg/L. Different doses (0.5-3.0 g/L) of HSC-MS-2 were added into 100 mL of mixed heavy metal solution and reacted at 25 • C for 4 h.

Column Sorption
The actual wastewater simulation was estimated by the column sorption trials. The column adsorptions were carried out in a glass column (11.0 mm in diameter and 150 mm in length). The bottom of column was firstly packed with quartz sand (100 meshes). Then, 1.0 g of HSC-MA-2 was inserted on top of the quartz sand layer. The actual wastewater  Figure S2). The Cu 2+ , Pb 2+ and Cr 3+ concentration were controlled at 100 mg/L with a pH of 5 to determine the effect of wastewater treatment on column sorption.

Gaussian Computations
All computations were performed by the Gaussian 09 Program (Revision D.01, Gaussian, Inc., Wallingford, CT, USA). All the geometries of the carbon domains with different functional groups were optimized at the B3LYP/6-31+G(d) level. The electrostatic potential (ESP) maps describe the electron density distribution, and the red area represents electron-rich areas, while the green area represents electron-deficient areas. According to the ESP maps, the geometries were further optimized at the B3LYP/6-31+G(d)/SDD level. The interaction energy between M and xxx were calculated by considering basis set superposition error (BSSE) correction.

Material Characterization
The pyrolysis process of HS biomass in the air or in a nitrogen atmosphere, and the pyrolysis process of HS in mixed salt with and without NaOH in the air were studied by using a synchronous thermal analyzer. Figure 1 shows the corresponding TG, DSC and DTG curves. For the pyrolysis process of HS in the air, Figure 1a shows that the first~10% of weight loss ends around~170 • C, which corresponds to the evaporation of moisture and decomposition of organic matter in HS. The second weight loss stage occurs in the range of 170-420 • C, with a total weight loss of~55% due to the pyrolysis of organic compounds and the burning of carbonized products. When the temperature rose to~600 • C, the weight loss of HS reached the maximum of~76.8%. When the temperature was elevated to 730 • C, the weight of HS did not change, which indicates that the organic components in HS are completely removed at high temperature without a salt seal. The residue corresponds tõ 23.2% of the initial biomass. Based on Figure 1b, the weight loss under N 2 is controlled at 67.6%, indicating that~32.4% of the biomass was retained after the calcination. Compared with that of the sample test in the air, it could be concluded that many mineral species are present in the biochars. For the mixture of HS and salt at a mass ratio of 1:2, it is found in Figure 1c that the weight loss of HS in the first three stages is basically the same as that in Figure 1a. Based on the calculation, the biomass weight loss at 600 • C is at~60%, which is much lower than that without the salt protection (~75.0%), indicating the salt sealing partially reduces the hydrolysis of the biomass. As depicted in Figure 1d, after the addition of NaOH, it can be observed that the weight loss experiences a slight increase from~60% to~64%, which indicates that the NaOH helps to cave the carbon structure during the high-temperature reaction. Based on these observations, a possible formation mechanism is proposed: due to the etching effect of the molten salt ions as well as the penetrated O 2 atoms at high temperature, marco-and meso-porous structures form in the biochar. In addition, the NaOH not only further increases the pore diameters, but also introduces abundant oxygen functional groups on the biochar surface.  Figure 2a shows an illustration of the preparation process of the HS-derived biochars. After the biomass was ground into a green powder, it was mixed with NaOH and the KCl/NaCl salt mix, and then sealed with KCl/NaCl in a crucible. Afterwards, the crucible was covered and heated to a given temperature for the pyrolysis. Finally, the product was washed with water to remove the excessive NaOH and salt. The biochars was dried for the adsorption of the heavy metal ions. The morphology and microstructure of the HS-derived biochars were studied using a field emission scanning electron microscope (FE-SEM). Figure 2b illustrates the excessive growth of HS in the north of China, which causes negative impacts on both agricultural production as well as the human activities. Figure 2c shows the FE-SEM image of the original biomass, indicating that the surface is smooth, with almost no pore structure detected. As depicted in Figure  2d, the surface of the biochar prepared in a nitrogen atmosphere is rougher, demonstrating that carbonation resulted in the formation of a porous structure in the as-yielded biochars (Supplementary Materials Figure S3). However, the porosities are closely packed together. However, as seen in Figure 2e, under the protection of molten salt, the HSC-MS exhibits an irregular sheet-like structure with a rougher surface. As shown in Supplementary Materials Figure S4, after the dual activation of molten salt and NaOH, the surfaces of the HSC-MAs all become rougher. In addition, the overall particle size of the HSC-MA-2 becomes smaller and the surface shows greater roughness, forming a rich porous structure (Figure 2f,g). Such porous structures allow the biochar to have large specific surface area, which can facilitate the mass diffusion in the pore channels of the biochar. The above results show that the formation of the porous structure of biomass results from the etching effects of both the molten salt and NaOH at high temperature. Meanwhile, the abundant porous structure and large specific surface area can promote the adsorption of heavy metal pollutants on the HS-derived biochars. Figure Figure 2a shows an illustration of the preparation process of the HS-derived biochars. After the biomass was ground into a green powder, it was mixed with NaOH and the KCl/NaCl salt mix, and then sealed with KCl/NaCl in a crucible. Afterwards, the crucible was covered and heated to a given temperature for the pyrolysis. Finally, the product was washed with water to remove the excessive NaOH and salt. The biochars was dried for the adsorption of the heavy metal ions. The morphology and microstructure of the HS-derived biochars were studied using a field emission scanning electron microscope (FE-SEM). Figure 2b illustrates the excessive growth of HS in the north of China, which causes negative impacts on both agricultural production as well as the human activities. Figure 2c shows the FE-SEM image of the original biomass, indicating that the surface is smooth, with almost no pore structure detected. As depicted in Figure 2d, the surface of the biochar prepared in a nitrogen atmosphere is rougher, demonstrating that carbonation resulted in the formation of a porous structure in the as-yielded biochars (Supplementary Materials Figure S3). However, the porosities are closely packed together. However, as seen in Figure 2e, under the protection of molten salt, the HSC-MS exhibits an irregular sheet-like structure with a rougher surface. As shown in Supplementary Materials Figure  S4, after the dual activation of molten salt and NaOH, the surfaces of the HSC-MAs all become rougher. In addition, the overall particle size of the HSC-MA-2 becomes smaller and the surface shows greater roughness, forming a rich porous structure (Figure 2f,g). Such porous structures allow the biochar to have large specific surface area, which can facilitate the mass diffusion in the pore channels of the biochar. The above results show that the formation of the porous structure of biomass results from the etching effects of both the molten salt and NaOH at high temperature. Meanwhile, the abundant porous structure and large specific surface area can promote the adsorption of heavy metal pollutants on the HS-derived biochars. Figure 2h depicts the TEM image of HSC-MA-2. The sample is composed of thin sheets with a large number of macropores with diameters of 50-100 nm. The TEM image with higher magnification (Figure 2i) shows a great number of mesoporous and micropores, due to the NaOH-assisted pyrolysis. The high-resolution TEM (HRTEM) image shown in Figure 2j indicates that many wrinkling lattice fringes can be detected at the surface of the carbon, indicating that the surface of the biochar was highly graphitized. indicates that many wrinkling lattice fringes can be detected at the surface of the carbon, indicating that the surface of the biochar was highly graphitized. X-ray photoelectron spectroscopy (XPS) was used to analyze the chemical elements of the biochars. Figure 3a shows the survey analysis of the HSC-N2, HSC-M and HSC-MA-2. There are two distinct peaks corresponding to C1s and O1s at the binding energies of 284.6 and 532.4 eV in the spectrum. In addition, HSC-MA-2 shows a much enhanced content of mineral species such as Na, Mg and Ca. The detailed element compositions are summarized in Table 1. X-ray photoelectron spectroscopy (XPS) was used to analyze the chemical elements of the biochars. Figure 3a shows the survey analysis of the HSC-N 2 , HSC-M and HSC-MA-2. There are two distinct peaks corresponding to C1s and O1s at the binding energies of 284.6 and 532.4 eV in the spectrum. In addition, HSC-MA-2 shows a much enhanced content of mineral species such as Na, Mg and Ca. The detailed element compositions are summarized in Table 1.   The C1s peak shown in Figure 3b indicates that HSC-N2 and HSC-M have a similar relative carbon content, while HSC-MA-2 has a relatively low carbon content, suggesting that the etching of the molten salt and NaOH could substantially reduce the carbon sp 2 domain and introduce abundant oxygen motifs. Based on the fitting of the C1s signals shown in Supplementary Materials Figure S5, it can be also observed that the peaks belonging to C-C and C=C (at 284.8 eV and 286.5 eV) are reduced. Meanwhile, the peaks ascribed to C=O and COOH (at 288.7 eV and 290.9 eV) are slightly enhanced, and the increased binding energy suggests that carbon motifs with high valence forms are incorporated in the HSC-MA-2 [32,33].
Moreover, this change can be also detected in the O1s spectrum displayed in Figure  3c. It is obvious that the content of oxygen in HSC-MA-2 is the highest, which indicates that the oxygen functional groups are introduced into the biochar by the NaOH modification reagent. The deconvolution of O 1s spectra shown in Supplementary Materials Figure S6 can be fit into four types: C=O (531.3 eV), OH/RCOOR (532.6 eV), O-C=O (533.7 eV) and COOH (534.9 eV). Compared with that of HSC-N2 and HSC-M, the peak ascribing to COOH is greatly enhanced in HSC-MA-2, indicating that carboxyl groups are successfully introduced into the biochars. This result is in good agreement with the C1s  The C1s peak shown in Figure 3b indicates that HSC-N 2 and HSC-M have a similar relative carbon content, while HSC-MA-2 has a relatively low carbon content, suggesting that the etching of the molten salt and NaOH could substantially reduce the carbon sp 2 domain and introduce abundant oxygen motifs. Based on the fitting of the C1s signals shown in Supplementary Materials Figure S5, it can be also observed that the peaks belonging to C-C and C=C (at 284.8 eV and 286.5 eV) are reduced. Meanwhile, the peaks ascribed to C=O and COOH (at 288.7 eV and 290.9 eV) are slightly enhanced, and the increased binding energy suggests that carbon motifs with high valence forms are incorporated in the HSC-MA-2 [32,33].
Moreover, this change can be also detected in the O1s spectrum displayed in Figure 3c. It is obvious that the content of oxygen in HSC-MA-2 is the highest, which indicates that the oxygen functional groups are introduced into the biochar by the NaOH modification reagent. The deconvolution of O 1s spectra shown in Supplementary Materials Figure  S6 can be fit into four types: C=O (531.3 eV), OH/RCOOR (532.6 eV), O-C=O (533.7 eV) and COOH (534.9 eV). Compared with that of HSC-N 2 and HSC-M, the peak ascribing to COOH is greatly enhanced in HSC-MA-2, indicating that carboxyl groups are successfully introduced into the biochars. This result is in good agreement with the C1s spectrum, and the increased carboxyl group population could in turn enhance the surface coordination effects on the heavy metal ions.
Figure 3d-f shows the contents of N, S and P in biochar. It can be seen that the biochars prepared from HS all contain a certain amount of N, P and S. Considering that no N-, Pand S-containing regents were used in the preparation, these heteroatoms originate from the pristine tissue of the HS biomass. Based on the previous reports, N, P and S with a lone pair of electrons show high affinity with the positively charged heavy metal ions [34]. Although HSC-MA-2 experience rigid MS and NaOH etching, the N and P could still be preserved in the biochar, which could further enhance the surface affinity towards the Pb and Cu ions. As depicted in Figure 3g-i, the modification of NaOH can greatly increase the content of Na, Mg and Ca in the biochar (the detail composition values are summarized in Table 1). These mineral ions are beneficial to the adsorption of heavy metals via the ion exchange process. For these biochars, no signal belonging to Pb and Cu could be detected, as shown in Supplementary Materials Figures S7 and S8.
In order to better study the element content in the biochars, we used inductively coupled plasma mass spectrometry (ICP-MS) to better reveal the metallic elements of the biochars. The detail results are shown in Table 2, and it is clear that the as-measured composition is much higher compared with that measured using XPS. Considering that XPS mainly detects the elements distributed on the surface of the biochars, the different composition results of ICP-MS demonstrate that the majority of the mineral species are embedded inside the biochars. Furthermore, HSC-MA-2 exhibits the highest contents in Na, Mg and Ca, which is in accordance with the XPS measurement. In order to analyze the graphitization degree and crystal phase of the biochars, the biochars were analyzed using laser Raman spectra, as shown in Figure 4a. There are characteristic diffraction peaks at~1320 and~1570 cm −1 , corresponding to the D band (disordered band) and G band (graphitized band) of the biochars [11]. The D-band is due to the defects and disorder induced by carbon atoms hybridizing, and the G-band is due to the stretching of carbon atoms. In general, the intensity ratio (I D /I G ) of the D and G bands is used to describe the disorder degree of biochars. The calculated D and G band intensity ratios of HSC-  (C=O, stretching vibration)). Therefore, the three peaks centered at~950,~1000 and~1300 cm −1 could be attributed to the C-O, OH and COOH or CHO groups, respectively. It can be seen that the intensity of oxygen functional group peaks, such as COOH or CHO, are greatly increased in HSC-MA-2 compared with other samples. The increased oxygen functional groups could provide abundant coordinating sites for the heavy metal ions on the biochars, which could greatly facilitate the surface adsorption process. [12,35]   Therefore, the three peaks centered at ~950, ~1000 and ~1300 cm −1 could be attributed to the C-O, OH and COOH or CHO groups, respectively. It can be seen that the intensity of oxygen functional group peaks, such as COOH or CHO, are greatly increased in HSC-MA-2 compared with other samples. The increased oxygen functional groups could provide abundant coordinating sites for the heavy metal ions on the biochars, which could greatly facilitate the surface adsorption process. [12,35] Surface Area and Pore Structure The pore structures of HSC-N2, HSC-M and HSC-MA-2 were studied by means of the nitrogen adsorption/desorption method. Figure 5 shows the nitrogen adsorption/desorption isotherms of the biochars, which belongs to the mixture of the type I and type IV isotherms, indicating the co-existence of the micro-and meso-porous structures. The micropore adsorption initially occurs at a low pressure. The amount of nitrogen adsorbed is between 0.4 and 0.9, corresponding to the capillary condensation of meso-pores, which indicates that HSC-MA-2 has a hierarchical porous structure [36]. In the range of relative pressure greater than 0.9, the nitrogen desorption isotherm of the HSC-MA-2 shows a continuous rising trend compared with that of the other two curves, indicating that there are still large pores in the biochar. The results of nitrogen adsorption and desorption show that a lot of micro-and meso-pores with different size distributions can be detected in the biochars prepared by molten salt with NaOH. As listed in Table 3, the BET specific surface area of the biochar experiences a slight reduction from 182.9 to 161.4 and 171.9 m 2 /g −1 . However, the overall pore volume and pore diameter increase with the addition of NaOH. The results indicate that the optimal HSC-MA-2 exhibits well-developed pore structures with a greater pore diameter and larger pore volume, which could benefit the mass diffusion and facilitate the adsorption process on the exposed inner surfaces.

Surface Area and Pore Structure
The pore structures of HSC-N 2 , HSC-M and HSC-MA-2 were studied by means of the nitrogen adsorption/desorption method. Figure 5 shows the nitrogen adsorption/desorption isotherms of the biochars, which belongs to the mixture of the type I and type IV isotherms, indicating the co-existence of the micro-and meso-porous structures. The micropore adsorption initially occurs at a low pressure. The amount of nitrogen adsorbed is between 0.4 and 0.9, corresponding to the capillary condensation of meso-pores, which indicates that HSC-MA-2 has a hierarchical porous structure [36]. In the range of relative pressure greater than 0.9, the nitrogen desorption isotherm of the HSC-MA-2 shows a continuous rising trend compared with that of the other two curves, indicating that there are still large pores in the biochar. The results of nitrogen adsorption and desorption show that a lot of micro-and meso-pores with different size distributions can be detected in the biochars prepared by molten salt with NaOH. As listed in Table 3, the BET specific surface area of the biochar experiences a slight reduction from 182.9 to 161.4 and 171.9 m 2 /g. However, the overall pore volume and pore diameter increase with the addition of NaOH. The results indicate that the optimal HSC-MA-2 exhibits well-developed pore structures with a greater pore diameter and larger pore volume, which could benefit the mass diffusion and facilitate the adsorption process on the exposed inner surfaces.

Adsorption Tests
The adsorption of Pb 2+ and Cu 2+ on different biochars were studied. In general, the Pb 2+ adsorption capacity of biochar prepared with HS ( Figure 6a) is much higher than that of Cu 2+ , which indicates that the as-prepared biochar has better affinity to Pb 2+ . In addition, the adsorption capacity increases with the NaOH amount, and the maximum value was reached when the mass ratio of biomass to NaOH equaled 3/2. As shown in Figure 6b, the pyrolysis temperature was controlled at 400, 500, 600, 700 and 800 • C to study the temperature influences. Obviously, the biochar prepared at 600 • C demonstrated high adsorption capacity, indicating that 600 • C is the optimal pyrolysis temperature for HS biomass. It is generally believed that the pH value affects the states of metal ions, the ionization degree of the solution as well as the surface charge of biochars. It is of great significance to study the effects of pH on adsorption performances [37]. To avoid precipitation, the pH's effects on the Pb 2+ and Cu 2+ ions adsorption capacity of the biochars were investigated at the pH range of 2.0-5.0 (Figure 6c). The adsorption capacity of the Cu 2+ and Pb 2+ increased with increasing pH values. At a low pH, the competition adsorption of H + on the surface sites greatly limit the adsorption capacity of Pb 2+ and Cu 2+ of the biochar. To make it worse, some minerals in the biochar may also dissolve and release large amounts of cations, such as Ca 2+ and Mg 2+ . These released ions may compete with Pb 2+ and Cu 2+ for the adsorption sites, hindering the adsorption capacity of heavy metal ions. Thus, the sorption kinetics and sorption isotherm experiments were carried out at a pH of 5.
The effects of HSC-MA-2 with different dosages on the adsorption capacities ( Figure  6d) and removal rates (Figure 6e) were studied to determine the optimal dosage. As expected, the adsorption capacity decreased with the increase in the biochar dosage. For Pb 2+ , the removal efficiency was ~75% when the dosage was 0.2 g/L −1 , and reached above ~90% when the dosage was elevated to 0.5 g/L −1 . For Cu 2+ , the removal efficiency reached over ~75% when the dosage was 0.5 g/L −1 , and ~90% removal efficiency was achieved when the dosage was 1.0 g/L −1 . Based on these analyses, the best dosages of Pb 2+ and Cu 2+ adsorbed by HSC-MA-2 are 0.2 g/L −1 and 0.5 g/L −1 , respectively. Considering that HSC- It is generally believed that the pH value affects the states of metal ions, the ionization degree of the solution as well as the surface charge of biochars. It is of great significance to study the effects of pH on adsorption performances [37]. To avoid precipitation, the pH's effects on the Pb 2+ and Cu 2+ ions adsorption capacity of the biochars were investigated at the pH range of 2.0-5.0 (Figure 6c). The adsorption capacity of the Cu 2+ and Pb 2+ increased with increasing pH values. At a low pH, the competition adsorption of H + on the surface sites greatly limit the adsorption capacity of Pb 2+ and Cu 2+ of the biochar. To make it worse, some minerals in the biochar may also dissolve and release large amounts of cations, such as Ca 2+ and Mg 2+ . These released ions may compete with Pb 2+ and Cu 2+ for the adsorption sites, hindering the adsorption capacity of heavy metal ions. Thus, the sorption kinetics and sorption isotherm experiments were carried out at a pH of 5.
The effects of HSC-MA-2 with different dosages on the adsorption capacities ( Figure 6d) and removal rates (Figure 6e) were studied to determine the optimal dosage. As expected, the adsorption capacity decreased with the increase in the biochar dosage. For Pb 2+ , the removal efficiency was~75% when the dosage was 0.2 g/L, and reached above~90% when the dosage was elevated to 0.5 g/L. For Cu 2+ , the removal efficiency reached over~75% when the dosage was 0.5 g/L, and~90% removal efficiency was achieved when the dosage was 1.0 g/L. Based on these analyses, the best dosages of Pb 2+ and Cu 2+ adsorbed by HSC-MA-2 are 0.2 g/L and 0.5 g/L, respectively. Considering that HSC-MA-2 has the best adsorption effect on Pb 2+ and Cu 2+ , the subsequent adsorption mechanism studies will focus on HSC-MA-2, adopting HSC-N 2 and HSC-M as references.

Adsorption Kinetics
The adsorption kinetics of Pb 2+ and Cu 2+ adsorption on different biochars are shown in Figure 7a-e. Table 4 depicts the kinetic parameters of the adsorption of Pb 2+ and Cu 2+ . The correlation coefficient (R 2 ) derived from the PSO model is much higher than that of the PFO model, which indicates that the adsorption of heavy metal ions is mainly dominated by the chemical adsorption, such as complexation and precipitation [4,5]. The adsorption rate is very fast at the beginning of the adsorption process, and then gradually reduces with the increase in contact time, until the adsorption equilibrium is reached after~120 min.  Table 4 depicts the kinetic parameters of the adsorption of Pb 2+ and Cu 2+ . The correlation coefficient (R 2 ) derived from the PSO model is much higher than that of the PFO model, which indicates that the adsorption of heavy metal ions is mainly dominated by the chemical adsorption, such as complexation and precipitation [4,5]. The adsorption rate is very fast at the beginning of the adsorption process, and then gradually reduces with the increase in contact time, until the adsorption equilibrium is reached after ~120 min.  The fitted curves of the intra-particle diffusion model are shown in Figure 7c,f. The intra-particle diffusion analysis shows that the adsorption of Pb 2+ and Cu 2+ by HSC-MA-2, HSC-M and HSC-N2 can be divided into three stages. Rapid adsorption rates are observed at the first stage, which results from the presence of many adsorption sites on the biochar surface [24,[37][38][39]. At the second stage, the adsorption is controlled by the intra-particle diffusion rate. Due to the decrease in heavy metal ion concentration and the very low mass transfer rate, the intra-particle diffusion process is slower than the first  The fitted curves of the intra-particle diffusion model are shown in Figure 7c,f. The intra-particle diffusion analysis shows that the adsorption of Pb 2+ and Cu 2+ by HSC-MA-2, HSC-M and HSC-N 2 can be divided into three stages. Rapid adsorption rates are observed at the first stage, which results from the presence of many adsorption sites on the biochar surface [24,[37][38][39]. At the second stage, the adsorption is controlled by the intra-particle diffusion rate. Due to the decrease in heavy metal ion concentration and the very low mass transfer rate, the intra-particle diffusion process is slower than the first stage. The third stage is the equilibrium stage, and the adsorption rate is greatly reduced [38]. Based on these observations, the HSC-MA-2 shows the fastest adsorption rate among the samples, indicating that the well-developed porous structure could benefit the intra-particle diffusion, thus enhancing the adsorption kinetics.

Adsorption Isotherm
The effects of the initial concentration of Pb 2+ and Cu 2+ on adsorption capacity were investigated. Due to the concentration gradient between the solution and the adsorbent surface, the adsorption capacity increases with the initial concentration of Pb 2+ and Cu 2+ . It is generally believed that a higher initial concentration of Pb 2+ and Cu 2+ provides a greater driving force to overcome the mass transfer resistance in the solution, which results in a higher adsorption capacity. When the initial concentrations of Pb 2+ and Cu 2+ are higher, the adsorption capacity is prone to reach a balance due to the limitation of the maximum adsorption capacities [39]. Freundlich and Langmuir models were used to describe the relationship between adsorption capacity and equilibrium concentration, and the fitting curves are shown in Figure 8. The as-calculated Q max , correlation coefficient (R 2 ), and other relevant constants are displayed in Table 5. The Langmuir isotherm has higher R 2 than the Freundlich model, which further proves that the chemical adsorption governs the removal process. The Q max of HSC-MA-2 is 748.1 and 221.1 mg/g for Pb 2+ and Cu 2+ , respectively. These results are in good agreement with the above experimental results, indicating that the adsorption capacity of the modified biochar for Pb 2+ and Cu 2+ is significantly improved compared to the biochar prepared directly by the molten salt sealing method without the modification reagent (594.5 and 155.0 mg/g for Pb 2+ and Cu 2+ ) or under N 2 protection (310.6 and 110.0 mg/g for Pb 2+ and Cu 2+ ). These results further indicate that under the conditions of an air atmosphere and high temperature, the etching effects of molten salt and NaOH on carbon products can endow the biochar with high specific surface area and rich oxygen functional groups, which is conducive to the adsorption capacity of heavy metal ions. In addition, as shown in Table 6, the optimized HSC-MA-2 demonstrates much enhanced adsorption capacity compared with that of the previously reported biochars prepared from different biomasses [40][41][42][43][44][45][46][47][48]. Table 5. Langmuir and Freundlich constants and correlation coefficients (R 2 ) for Pb 2+ and Cu 2+ adsorption on HSC-N 2 , HSC-M and HSC-MA-2.

Sample
Langmuir Freundlich under the conditions of an air atmosphere and high temperature, the etching effects of molten salt and NaOH on carbon products can endow the biochar with high specific surface area and rich oxygen functional groups, which is conducive to the adsorption capacity of heavy metal ions. In addition, as shown in Table 6, the optimized HSC-MA-2 demonstrates much enhanced adsorption capacity compared with that of the previously reported biochars prepared from different biomasses [40][41][42][43][44][45][46][47][48].

Effects of Competition Ions
In the adsorption process, the background ions (Na + , Mg 2+ and Ca 2+ ) affect the adsorption of target ions on the biochar. As shown in Figure 9, the results indicate that Na + has little effect on the removal of Pb 2+ and Cu 2+ , while the presence of Mg 2+ and Ca 2+ has a greater influence on the adsorption capacity of HSC-MA-2 compared with Na + (monovalent cation). As the concentration of background ions is increased, the removal efficiencies of Pb 2+ and Cu 2+ decrease due to the enhanced ion strength of Mg 2+ and Ca 2+ , which could be more competitive for the active sites on HSC-MA-2. In addition, the influence of Ca 2+ on the adsorption of Pb 2+ and Cu 2+ is greater than that of Mg 2+ . The higher valence sates of Mg 2+ and Ca 2+ may lead to stronger hydration, which makes Mg-OH more stable than Ca 2+ in aqueous solutions and thus reduces its interaction with the surface active sites [39]. Therefore, the effect of Mg 2+ on Pb 2+ and Cu 2+ adsorption by HSC-MA-2 is lower than that of Ca 2+ , but greater than that of Na + .

Multi-Metal Adsorption and Column Test
The interaction of coexisting heavy metal ions in the adsorption process has more practical significance than in the single heavy metal ion system. Therefore, the adsorption capacity in the poly-metal ion system was studied (Figure 10a-c). As expected, the removal efficiencies of Pb 2+ , Cu 2+ , and Cr 3+ increase with the increase in the biochar dosage. For HSC-MA-2, when the biochar dosage is increased from 50 to 300 mg, the removal efficiency of all heavy metal ions is enhanced significantly. The removal efficiency of all heavy metal ions decreases with the increase in the concentration of metal ions. As the dosage of HSC-MA-2 is increased to 300 mg, the removal efficiency of all heavy metal ions reaches ~100%. These results suggest that when the metal ions compete for the same adsorption sites, the metal ions with higher affinity replace the other ions with

Multi-Metal Adsorption and Column Test
The interaction of coexisting heavy metal ions in the adsorption process has more practical significance than in the single heavy metal ion system. Therefore, the adsorption capacity in the poly-metal ion system was studied (Figure 10a-c). As expected, the removal efficiencies of Pb 2+ , Cu 2+ , and Cr 3+ increase with the increase in the biochar dosage. For HSC-MA-2, when the biochar dosage is increased from 50 to 300 mg, the removal efficiency of all heavy metal ions is enhanced significantly. The removal efficiency of all heavy metal ions decreases with the increase in the concentration of metal ions. As the dosage of HSC-MA-2 is increased to 300 mg, the removal efficiency of all heavy metal ions reaches~100%. These results suggest that when the metal ions compete for the same adsorption sites, the metal ions with higher affinity replace the other ions with lower affinity, and the most intense competition takes place when there is a high concentration of heavy metal ions. For HSC-MA-2 (200 mg), the removal efficiencies of Pb 2+ , Cu 2+ and Cr 3+ are 98.05%, 94.57% and 96.18%, respectively, when the concentration of all metal ions is set at 150 mg/L. It is evidenced that the presence of Cr 3+ has a significant effect on the adsorption of Pb 2+ and Cu 2+ on the biochar, due to the strong affinity for Cr 3+ . The experimental results prove that the HSC-MA-2 can effectively take up many types of heavy metal ions at the same time. Based on the abovementioned batch experiments, column adsorption experiments are also carried out in order to investigate the possible application of HSC-MA-2 in large-scale operations for the practical wastewater treatment. The column adsorption experiments were conducted using HSC-MA-2 as the filler (Figure 10d). The simulated wastewater was prepared from the Yellow River; 100 mg/L Cu 2+ , Pb 2+ and Cr 3+ were added, and the average discharge of the column was adjusted to 0.4 mL/min. As the operation time increases, heavy metal ions adsorption gradually decreases. This phenomenon indicates that HSC-MA-2 is gradually saturated regarding heavy metal ions' adsorption capacity. Cu 2+ and Cr 3+ become saturated within 600 min, while the transition of Pb 2+ occurs at 500 min and reaches saturation at 720 min. Generally, the longer the saturation time, the better the adsorption effect of the metal ions. Therefore, in the simulated column adsorption of practical wastewater, the affinity of HSC-MA-2 to different heavy metal ions follows the sequence of Pb 2+ > Cu 2+ > Cr 3+ . Considering the outstanding removal effects for the heavy metal ions, the HSC-MA-2 provides a highly efficient removal rate for the practical wastewater treatment. Generally, the longer the saturation time, the better the adsorption effect of the metal ions. Therefore, in the simulated column adsorption of practical wastewater, the affinity of HSC-MA-2 to different heavy metal ions follows the sequence of Pb 2+ > Cu 2+ > Cr 3+ . Considering the outstanding removal effects for the heavy metal ions, the HSC-MA-2 provides a highly efficient removal rate for the practical wastewater treatment.

Adsorption Mechanism
As indicated in previous reports, the sorption mechanism probably involved surface complexation, metal precipitation, ion exchange, π-π interactions, and physical sorption [39,42]. In order to investigate the high sorption capacity of HSC-MA-2, the biochars be-

Adsorption Mechanism
As indicated in previous reports, the sorption mechanism probably involved surface complexation, metal precipitation, ion exchange, π-π interactions, and physical sorption [39,42]. In order to investigate the high sorption capacity of HSC-MA-2, the biochars before and after metal ion adsorption were investigated in detail. As shown in Figure 11a, the Raman signal is greatly subtracted after the adsorption of the heavy metal ions, especially for the case of the Pb 2+ adsorption, which indicates that the signal is greatly affected by the interaction between the biochar and the positively charged heavy metal ions. The reduced I G signal reflects the deviation of delocalized π electrons in the carbon domain, suggesting the strong π-π interaction between the biochar and the metal ions. In particular, the I D intensity is much more reduced than that of the I G , which indicates that the defects remaining in the carbon domain serve as the most effective sites for the adsorption. This phenomenon can be further proved by the FTIR spectra. As depicted in Figure 11b, the absorption peaks centered at~950,~1000 and~1300 cm −1 are greatly abated, which is also attributed to the coordination of the oxygen group with the metal ions. Furthermore, the high-resolution XPS spectra on the C1s and O1s are disclosed in Figure 11c,d. The reduced C1s peaks can be explained by the strong interaction on the surface, which results from the heavy adsorption loading of the metal ions on the biochar. In addition, the O1s can be delineated into three peaks, corresponding to the OH, C=O and COOH, respectively. It is obvious that after the Pb 2+ adsorption, the peaks ascribed to C=O and COOH are greatly reduced, which indicates that the C=O and COOH groups play significance roles in the heavy metal adsorption. Based on the Gauss simulation, COOH and C=O more effectively lead to the charge polarization on the carbon surface, which gives rise to the higher surface affinity towards the metal ions (362.8 kcal/mol) (Figure 11e    In order to investigate the ion exchange during the adsorption process, the release of Ca 2+ and Mg 2+ ions was monitored during the adsorption process (Figure 12a,b). The Ca 2+ and Mg 2+ concentrations experience continuous enhancement with prolonged adsorption In order to investigate the ion exchange during the adsorption process, the release of Ca 2+ and Mg 2+ ions was monitored during the adsorption process (Figure 12a,b). The Ca 2+ and Mg 2+ concentrations experience continuous enhancement with prolonged adsorption times, which suggests that the release of Ca 2+ and Mg 2+ is accompanied by the heavy metal ion adsorption. In addition, the amount of Ca 2+ released is much higher than that of the Mg 2+ , implying that Ca 2+ is the most effective mineral species helping the adsorption process. Moreover, the high-resolution XPS spectra shown in Figure 12c,d indicate that after the Pb 2+ and Cu 2+ adsorption, the Ca 2+ and Mg 2+ content in the spectra are greatly suppressed, especially for the Pb 2+ adsorption case. Taking into account the reduced Ca 2+ and Mg 2+ remaining in the biochar, as well as the enhanced concentration of Ca 2+ and Mg 2+ ions in the final adsorption system, it can be concluded that ion exchange plays an important role in facilitating heavy metal ion adsorption.

Conclusions
The NaOH-modified molten salt method was adopted to prepare a series of biochars using Humulus scandens as a biomass precursor. The optimized biochar HSC-MA-2 exhibits excellent Pb 2+ and Cu 2+ adsorption performance, and the Qm can reach 748 and 221 mg/g −1 , respectively. In addition, HSC-MA-2 shows good anti-interference ability and high removal efficiency in the simulated wastewater using Yellow River water for practical application. The enhanced performance of the HSC-MA-2 can be ascribed to the following qualities of the proposed NaOH-modified molten salt process. Firstly, the synthetic etching effect of the molten salt and the NaOH media introduce a large surface area and a large porous structure to offer abundant adsorption sites as well as mass diffusion channels. Secondly, the penetrating O2 molecules and the NaOH endow the carbon surface with a large population of the oxygen functional groups such as C=O and COOH, which substantially increase the coordinating strength with the positively charged metal ions. Finally, the NaOH helps to preserve the pristine mineral species from the biomass in the biochar, which serve as the active centers to adsorb the foreign heavy metal ions through ion exchange. Therefore, this method could serve as a general protocol to produce high-performance biochar for water treatment.

Conclusions
The NaOH-modified molten salt method was adopted to prepare a series of biochars using Humulus scandens as a biomass precursor. The optimized biochar HSC-MA-2 exhibits excellent Pb 2+ and Cu 2+ adsorption performance, and the Q m can reach 748 and 221 mg/g, respectively. In addition, HSC-MA-2 shows good anti-interference ability and high removal efficiency in the simulated wastewater using Yellow River water for practical application. The enhanced performance of the HSC-MA-2 can be ascribed to the following qualities of the proposed NaOH-modified molten salt process. Firstly, the synthetic etching effect of the molten salt and the NaOH media introduce a large surface area and a large porous structure to offer abundant adsorption sites as well as mass diffusion channels. Secondly, the penetrating O 2 molecules and the NaOH endow the carbon surface with a large population of the oxygen functional groups such as C=O and COOH, which substantially increase the coordinating strength with the positively charged metal ions. Finally, the NaOH helps to preserve the pristine mineral species from the biomass in the biochar, which serve as the active centers to adsorb the foreign heavy metal ions through ion exchange. Therefore, this method could serve as a general protocol to produce high-performance biochar for water treatment.