Mistakes and inconsistencies regarding adsorption of contaminants from aqueous solutions: A critical review

Mistakes and inconsistencies regarding adsorption of contaminants from aqueous solutions: A critical review

Accepted Manuscript Mistakes and inconsistencies regarding adsorption of contaminants from aqueous solutions: A critical review Hai Nguyen Tran, Sheng...

13MB Sizes 2 Downloads 126 Views

Accepted Manuscript Mistakes and inconsistencies regarding adsorption of contaminants from aqueous solutions: A critical review Hai Nguyen Tran, Sheng-Jie You, Ahmad Hosseini-Bandegharaei, Huan-Ping Chao PII:

S0043-1354(17)30269-5

DOI:

10.1016/j.watres.2017.04.014

Reference:

WR 12811

To appear in:

Water Research

Received Date: 29 October 2016 Revised Date:

29 March 2017

Accepted Date: 6 April 2017

Please cite this article as: Tran, H.N., You, S.-J., Hosseini-Bandegharaei, A., Chao, H.-P., Mistakes and inconsistencies regarding adsorption of contaminants from aqueous solutions: A critical review, Water Research (2017), doi: 10.1016/j.watres.2017.04.014. This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

ACCEPTED MANUSCRIPT

Mistakes and inconsistencies regarding adsorption of contaminants from aqueous

2

solutions: A critical review

3 4

Hai Nguyen Trana,b*, Sheng-Jie Youb, Ahmad Hosseini-Bandegharaeic,d, Huan-Ping Chaob*

RI PT

1

5 6 a

8 9

b

10 11 12 13 14 15

c

Department of Civil Engineering, Chung Yuan Christian University, Chungli 320, Taiwan

SC

Department of Environmental Engineering, Chung Yuan Christian University, Chungli 320, Taiwan Wastewater Division, Faculty of Health, Sabzevar University of Medical Sciences, PO Box 319, Sabzevar, Iran d Department of Engineering, Kashmar Branch, Islamic Azad University, PO Box 161, Kashmar, Iran

16

M AN U

7

Corresponding authors:

18

H. N. Tran ([email protected]) and H.- P. Chao ([email protected])

19 20 21

Department of Environmental Engineering, Chung Yuan Christian University, Chungli 320, Taiwan

EP AC C

22

TE D

17

1

ACCEPTED MANUSCRIPT

Abstract

24 25

In recent years, adsorption science and technology for water and wastewater treatment has attracted substantial attention from the scientific community. However, the number of

26 27

publications containing inconsistent concepts is increasing. Many publications either reiterate previously discussed mistakes or create new mistakes. The inconsistencies are reflected by the

28 29

increasing publication of certain types of article in this field, including “short communications”, “discussions”, “critical reviews”, “comments”, “letters to the editor”, and

30

“correspondence (comment/rebuttal)”. This article aims to discuss (1) the inaccurate use of

31 32 33 34 35 36 37 38 39

technical terms, (2) the problem associated with quantities for measuring adsorption performance, (3) the important roles of the adsorbate and adsorbent pKa, (4) mistakes related to the study of adsorption kinetics, isotherms, and thermodynamics, (5) several problems related to adsorption mechanisms, (6) inconsistent data points in experimental data and model fitting, (7) mistakes in measuring the specific surface area of an adsorbent, and (8) other mistakes found in the literature. Furthermore, correct expressions and original citations of the relevant models (i.e., adsorption kinetics and isotherms) are provided. The authors hope that this work will be helpful for readers, researchers, reviewers, and editors who are interested in the field of adsorption studies.

40 41 42

Keywords: Adsorption; mistake; comment; inconsistency; critical review

AC C

EP

TE D

M AN U

SC

RI PT

23

2

ACCEPTED MANUSCRIPT

43

Table of Contents 1. Introduction ........................................................................................................................5 2. Technical terms used in the study of adsorption ..................................................................6

47 48

3. Quantities for measuring adsorption performance ...............................................................8 4. Incorrect assumptions regarding pKa ...................................................................................9

49

5. Adsorption kinetics ...........................................................................................................11

50

5.1. The important role of initial time............................................................................11

51 52

5.2. Adsorption equilibrium time for porous materials ..................................................13 5.3. Pseudo-first-order (PFO) equation..........................................................................14

53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79

5.3.1. Derivation of the PFO equation ...................................................................14 5.3.2. Problems in the application of the PFO equation..........................................15 5.4. Pseudo-second-order (PSO) equation .....................................................................17 5.4.1. Derivation of the PSO equation ...................................................................17 5.4.2. Problems in the application of the PSO equation..........................................20 5.5. Elovich equation ....................................................................................................23 5.5.1. Derivation of the Elovich equation ..............................................................23 5.5.2. Incorrect expression of the Elovich equation................................................24 5.6. Intra-particle diffusion model .................................................................................24 6. Adsorption isotherms ........................................................................................................25 6.1. Adsorption equilibrium ..........................................................................................25 6.2. Langmuir equation .................................................................................................29 6.2.1. Derivation of the Langmuir equation ...........................................................29 6.2.2. Derivation of the separation factor ...............................................................31 6.2.3. Mistakes concerning the Langmuir model and the separation factor.............31 6.3. Freundlich equation ...............................................................................................33 6.3.1. Derivation of the Freundlich equation ..........................................................33 6.3.2. Mistakes concerning the Freundlich equation ..............................................34 6.4. Redlich–Peterson equation .....................................................................................35 6.4.1. Derivation of the Redlich–Peterson equation ...............................................35 6.4.2 Mistakes concerning the Redlich–Peterson equation .....................................36 6.5. Dubinin–Radushkevich equation ............................................................................36 6.5.1. Derivation of the Dubinin–Radushkevich equation ......................................36 6.5.2. Mistakes concerning the Dubinin–Radushkevich equation ...........................37 7. Adsorption thermodynamics .............................................................................................38 7.1. Principles of adsorption thermodynamics ...............................................................38 7.2. Equilibrium constant derived from the Langmuir constant (KL) ..............................40 3

AC C

EP

TE D

M AN U

SC

RI PT

44 45 46

ACCEPTED MANUSCRIPT

80

7.3. Equilibrium constant derived from the Freundlich constant (KF) ............................44 7.4. Equilibrium constant derived from the partition coefficient (Kp).............................45 7.5. Equilibrium constant derived from the distribution coefficient (Kd) ........................47 8. Others mistakes ................................................................................................................49

84 85

8.1. Inconsistent data points in experimental data and model fitting ..............................49 8.2. Oxidation state of chromium ..................................................................................51

86

8.3. Incorrect labels.......................................................................................................53

87

8.4. BET specific surface area of an adsorbent ..............................................................54

88

8.5. Maximum absorption wavelength in dye adsorption studies ...................................57

89 90 91 92 93 94

8.6. Cπ-cation and π-π interactions................................................................................57 8.6.1. Cπ-cation interactions ..................................................................................57 8.6.2. π-π interactions ............................................................................................59 8.7. Other miscellaneous errors .....................................................................................62 9. Nonlinear-optimization technique .....................................................................................66 10. Conclusions ....................................................................................................................68

95 96 97

Acknowledgements ..............................................................................................................69

M AN U

SC

RI PT

81 82 83

AC C

EP

TE D

98

4

ACCEPTED MANUSCRIPT

99

1. Introduction Nowadays, the number of publications in international scientific journals has increasingly

101

become a standard for the assessment of scientists (Geckeis and Rabung, 2004). As the number of

102

scientific publications increase, mistakes and misconceptions enter scientific literature and some

103

can also be repeated in subsequent publications. Once errors and mistakes enter the literature, it is

104

difficult to eradicate them. The entrance of mistakes in scientific publications, as well as their

105

propagation, can result from the subjectivity of the authors and/or objective reasons.

SC

RI PT

100

A typical study on the adsorption of organic and inorganic contaminants comprises several

107

sections, such as the effects of solution pH and ionic strength, studies of adsorption kinetics,

108

isotherms, thermodynamics, desorption, and regeneration. Adsorption processes can be conducted

109

using either batch or column techniques. It has been found that many scientific articles in the field

110

of adsorption science and technology contain several mistakes and misconceptions including the

111

use of technical terms, application of quantities, determining the roles of some constants like pKa,

112

performing kinetic studies and modeling kinetic adsorption data, finding a suitable isotherm model

113

to fit equilibrium data, calculating thermodynamic parameters, correct interpretation of adsorption

114

mechanisms, and measuring the specific surface area of an adsorbent. Moreover, a survey of the

115

scientific publications related to the adsorption of contaminants in aqueous solutions reveals that

116

some less obvious inconsistencies and mistakes have unavoidably slipped the attention of authors

117

and were repeated in subsequent publications.

TE D

EP

AC C

118

M AN U

106

The observation of these inconsistencies and mistakes, regardless of their origin, along

119

with comments and open discussions can prevent their propagation in scientific knowledge

120

transfer. Therefore, the presentation of a comprehensive review on common mistakes in

121

adsorption studies, besides introducing the most correct approaches to use and giving a source of 5

ACCEPTED MANUSCRIPT

122

up-to-date literature on this issue, seems to be valuable to both readers and researchers in this

123

field. The main purpose of this review is to identify mistakes in past publications and prevent the

125

propagation of such mistakes in future scientific literature. Herein, we present some common

126

mistakes in various aspects of adsorption studies and analyze the reasons for their entrance into

127

the literature. This was accomplished by discussing the published papers in the field of adsorption

128

of contaminants from aqueous solutions, commonly including “comments on”, “reply for

129

comments”, “response to”, “authors’ response to comments on”, “comments on the authors’

130

response to comments on”, “comments on comment on”, “in reaction to”, “critical review”, and “a

131

note on the comments” published works. In addition to our comprehensive review on the mistakes

132

existing in publications, some correct approaches are suggested where necessary to avoid the

133

propagation of such mistakes in the future scientific literature. This paper primarily focuses on

134

three basic constituents of adsorption theory: adsorption equilibria, kinetics, and thermodynamics.

135

2. Technical terms used in the study of adsorption

TE D

M AN U

SC

RI PT

124

Adsorption processes have their own vernacular. Therefore, correctly understandings the

137

technical terms used in adsorption technology can prevent the introduction of several unexpected

138

ambiguities and discrepancies. Some basic adsorption terms are summarized in Figure 1.

140

AC C

139

EP

136

Figure 1

To some extent, a thorough understanding of the properties of adsorbents might prevent

141

common mistakes made in interpreting adsorption processes and adsorption mechanisms.

142

Numerous techniques can be used to characterize an adsorbent (Figure 2). More information on the

143

basic properties of adsorbents (i.e., biosorbent, biochar, hydrochar, and activated carbons) and the

144

relations between these properties has been published in our recent work (Tran et al. 2017c). 6

ACCEPTED MANUSCRIPT

145

Figure 2 The appropriate use of technical words regarding the biosorption mechanism of heavy

147

metals has been highlighted by Robalds et al. (2016). To avoid reader confusion, they

148

recommended that “surface precipitation” be used instead of “microprecipitation” (also written as

149

“micro-precipitation”). Similarly, it is necessary to distinguish between “chelation”,

150

“complexation”, and “coordination”. The word “complexation” includes two meanings:

151

“chelation” and/or “coordination”. For example, the incorrect sentence “in another sense, it can

152

also be defined as a collective term for a number of passive accumulation processes which in any

153

particular case may include ion exchange, coordination, complexation, chelation, adsorption and

154

microprecipitation” should be revised as “biosorption may include ion exchange, complexation

155

(including coordination and/or chelation), physisorption, or microprecipitation” (Robalds et al.,

156

2016). In addition, the use of the term “electrostatic adsorption” may cause misunderstanding and

157

confusion; therefore, depending on the intended meaning, the appropriate terms (i.e., “electrostatic

158

attraction”, “electrostatic force”, or “electrostatic interaction”) are recommended. Notably, the

159

authors also suggested a new classification system for the (bio)sorption mechanisms of heavy

160

metals (Figure 3).

SC

M AN U

TE D

EP

162

Figure 3

The electrical state of an adsorbent’s surface in solution is usually characterized by either

AC C

161

RI PT

146

163

the point of zero charge (PZC) or the isoelectric point (IEP). Many authors have complained about

164

the misuse of the terms IEP and PZC outside their normal meanings (Somasundaran, 1968;

165

Kosmulski, 2009). The PZC is defined as the solution conditions under which the surface charge

166

density equals zero. Potentiometric titration, or related methods, is used to determine the PZC by

167

finding the point at which the apparent surface charge density in the presence of an inert electrolyte 7

ACCEPTED MANUSCRIPT

is not dependent on ionic strength (Figure 4a). The IEP occurs when the electrokinetic (ζ) potential

169

at the shear plane equals zero; ζ is determined from measuring electrokinetic (electrophoresis and

170

streaming potential) and electroacoustic phenomena (Figure 4b) (Kosmulski, 2009). IEP values

171

clearly represent only the external surface charge of particles in solution, whereas the PZC varies

172

in response to the net total (external and internal) surface charge of the adsorbent (as the

173

characteristic of amphoteric surfaces). Therefore, the difference between the PZC and IEP values

174

for a porous carbon system can be interpreted as a measure of the surface charge distribution. A

175

difference (PZC − IEP) of greater than zero indicates that the external particle surface is more

176

negatively charged than the internal particle surface, and a difference of close to zero corresponds

177

to a more homogeneous distribution of surface charge. The PZC is only equal to the IEP in the

178

absence of specific or counter ions. This means that equality between the PZC and IEP exists only

179

if no specific adsorption of counter ions from the solution occurs (Radovic, 1999; Onjia and

180

Milonjić, 2002).

TE D

M AN U

SC

RI PT

168

181 182

Figure 4

3. Quantities for measuring adsorption performance In the field of adsorption, adsorption performance can be expressed as the amount of

184

adsorbate adsorbed at equilibrium or the percentage of removed adsorbate. The amount of

185

adsorbate adsorbed at equilibrium (qe; mg/g) is often calculated using the material balance of an

186

adsorption system; the adsorbate, which has disappeared from the solution, must be in the

187

adsorbent (Volesky, 2007).

AC C

EP

183

qe =

(C o − C e ) V1 m1

(1)

8

ACCEPTED MANUSCRIPT

188

where C0 (mg/L) and Ce (mg/L) are the initial and equilibrium adsorbate concentration in solution,

189

respectively; m (g) is the dried mass of used adsorbent; and V (L) is the volume of the adsorbate

190

solution. Depending on the purpose of the study, the parameter qe may be expressed in different units

192

as follows: (1) mg/g for evaluating practical and engineering processes, (2) mmol/g or meq/g for

193

examining the stoichiometry of a process and for studying functional groups and metal-binding

194

mechanisms, and (3) mmol/g for comparing the selective adsorption performance of various

195

adsorbates (Volesky, 2007). Of these units, mg/g (milligrams of adsorbate adsorbed per gram of

196

dried adsorbent) is used most commonly for qe in adsorption studies.

M AN U

SC

RI PT

191

In addition to qe, the adsorption performance can also be expressed as the percentage of

198

removed adsorbate (%removal). However, the unit of %removal needs to be used cautiously, as it

199

is very approximate and can cause misleading conclusions about relative adsorption performance.

200

This unit is only appropriate for the purpose of crude measurements and perhaps for quick and

201

very approximate screening of adsorbent materials (Kratochvil and Volesky, 1998; Hai, 2017) (C o − C e ) × 100 Co

(2)

EP

% removal =

TE D

197

To some extent, the use of %removal in the study of adsorption equilibria can cause

203

incorrect observations and, as a result, inaccurate conclusions (a detailed discussion is provided in

204

Section 6.1).

205

4. Incorrect assumptions regarding pKa

206

AC C

202

A mistake regarding the fundamentals of chemistry was commented on by Rayne (2013).

207

The commented paper reported the adsorption of n-perfluorooctanesulfonamide (PFOSA) onto

208

multiwalled carbon nanotubes with differing oxygen content. However, PFOSA was mistakenly 9

ACCEPTED MANUSCRIPT

assumed to be an “organic base”, resulting in the following invalid interpretation: “for PFOSA

210

(pKa = 6.52), when pH < pKa, protonation occurs on the amino group, and the decreased

211

protonation leads to the increased adsorption, but when pH > pKa, PFOSA exists as neutral

212

molecule in water.” In fact, PFOSA is an acid and, thus, is a neutral compound below its pKa value

213

(estimated at between 6.0 and 7.0) and a dissociated anion at pH > pKa (Rayne, 2013).

215

Recently, the effect of pH on the interaction of copper with the surface functional groups of raw pomegranate-peel biosorbent was reported as follows:

R–OH ↔ R–O− + H+ R–O− + Cu2+ ↔ R–OCu+ R–O− + Cu(OH)+ ↔ R–OCu(OH)

M AN U

R–OH2+ ↔ R–OH + H+

SC

214

RI PT

209

(3) (4) (5) (6)

where –R represents the surface adsorbent; R–OH2+, R–OH, and R–O− represent protonated,

217

neutral, and ionized surface hydroxyl functional groups, respectively; and R–OCu+ and R–

218

OCu(OH) are the formed bonding complexes.

TE D

216

The authors (Ben Ali et al., 2016) explained that, “the effect of pH could be also explained

220

by considering the point of zero charge of the adsorbent (pHpzc). The adsorbent surface is

221

positively charged for pH< 4.7 and it becomes negatively charged at pH value above 4.7.

222

Therefore, for pH values < 4.7, the adsorption is unfavorable because of repulsive electrostatic

223

interactions between metal ions and positively charged functional groups. The maximum

224

adsorption of Cu(II) occurs at pH above pHpzc value when the adsorbent surface is negatively

225

highly charged.” However, the important role of the pKa values of carboxylic acid groups (–

226

COOH; pKa 1.7–4.7) and hydroxyl groups (–OH; pKa 9.5–13) was ignored (Volesky, 2007).

227

Correspondingly, Eqs. 3–6 indicate a misunderstanding of the surface chemistry of the adsorbent.

AC C

EP

219

10

ACCEPTED MANUSCRIPT

The authors investigated the effect of solution pH on the adsorption process of Cu(II) in solutions

229

with pH values from 2.0 to 6.0. Therefore, dissociation of the –OH groups to –O− did not occur at

230

these solutions pH values. Instead, the –COOH groups dissociated, forming negatively charged

231

carboxylate groups (–COO‒) in the aforementioned pH range, and were the dominant groups that

232

interacted with copper ions (Tran, 2017).

RI PT

228

Moreover, Bui and Choi (2010) made a critical comment that Oleszczuk et al (2010)

234

mistyped the properties of two pharmaceuticals: oxytetracycline (logKow = 2.45, pKa = 7) and

235

carbamazepine (logKow = −1.22, pKa = 3.27, 7.32, 9.11). The correct expression should be

236

oxytetracycline (Sw = 121 mg/L, logKow = −1.22, pKa = 3.27, 7.32, 9.11; density = 1.63 g/cm3) and

237

carbamazepine (Sw = 112 mg/L, logKow = 2.45, pKa = 7; density =1.15 g/cm3). Bui and Choi

238

(2010) also emphasized that the pKa value of carbamazepine reported by Oleszczuk et al (2010) is

239

still debated. Carbamazepine (CBZ) is commonly characterized with two pKa values: pKa1 (= 1 or

240

2.3) for the equilibrium RCONH3+ ⇌ RCONH2 + H+ and pKa2 (= 13.9) for the equilibrium

241

RCONH2 ⇌ RCONH− + H+. The misreporting of the CBZ pKa can cause inaccurate conclusions.

242

For example, Oleszczuk et al. (2010) suggested that an electrostatic interaction occurred between

243

CBZ and multiwalled carbon nanotubes (MWCNT) at pH < 7.0 because they opined that CBZ

244

existed in its cationic form at pH < 7.0 (pKa = 7.0). However, CBZ with pKa1 < 2.3 and pKa2 ≥ 13.9

245

actually exists entirely as a neutral species in the pH range of 1.9–11.9, as calculated using

246

ACDLab. As a result, the electrostatic interaction was ruled out as a contributing interaction in the

247

CBZ-adsorption mechanisms (Bui and Choi, 2010).

248

5. Adsorption kinetics

249

5.1. The important role of initial time

250

AC C

EP

TE D

M AN U

SC

233

In the study of adsorption kinetics, Azizian (2006) proposed that experiments should be 11

ACCEPTED MANUSCRIPT

251

started at t < 2 min (considered as the initial time). Adsorption kinetics data in the initial time

252

periods play an important role in performing accurate modeling and drawing conclusions. Depending on the affinity between the adsorbent and adsorbate, the initial adsorption rate

254

may be very fast. For the adsorption of heavy metals, Tran et al. (2015) investigated the adsorption

255

of Cd2+ onto orange-peel-generated biochar, and found that equilibrium can be reached rapidly in

256

kinetic experiments with a removal rate of 80.6–96.9% (within 1 min). An analogous performance

257

was reported by Guo et al. (2015), with approximately 70.0–96.6% of the total Cd2+ in solution

258

removed within the first minute of contact. Additionally, Müller et al. (2012) noted that

259

approximately 90% of Zn(II) adsorption onto Lewatit S1468 occurs within 1 min, and this

260

adsorption process reaches equilibrium in approximately 5 min.

M AN U

SC

RI PT

253

For the adsorption of organic contaminants, kinetic studies indicated that the adsorption

262

equilibrium of methylene green 5 was quickly established, with a low activation energy required

263

for adsorption (Ea; 3.30–27.8 kJ/mol), and activated carbons removed 50–73% of the dye from

264

solution within 1 min (Tran et al., 2017d; Tran et al., 2017e). Similarly, Canzano et al. (2012)

265

noted that the adsorption process of anionic Congo-red dye onto raw and acid-treated pine cone

266

powder reached equilibrium in the first few minutes (< 5 min), during which no measurements

267

were carried out.

EP

As time goes on, the amount of adsorbate uptake per mass of adsorbent at any time t (qt)

AC C

268

TE D

261

269

approaches that of at equilibrium (qe) and, therefore, plotting the kinetic data is bound to produce a

270

straight line, independent of the real kinetic order followed by the system. In other words,

271

equilibrium data does not describe the way by which the system reaches equilibrium. As a result,

272

the use of data recorded over a significant time interval after the attainment of equilibrium, or very

273

close to equilibrium, is likely to lead to erroneous conclusions regarding adsorption kinetics. 12

ACCEPTED MANUSCRIPT

274

5.2. Adsorption equilibrium time for porous materials The adsorption processes for inorganic and organic contaminants onto porous adsorbents

276

(i.e., activated carbon, biochar, activated carbon spheres, zeolite, and macroreticular resin) can

277

take longer (i.e., several days or weeks) to approach true equilibrium than nonporous adsorbents

278

(i.e., biosorbent and hydrochar). This difference is attributed to the different adsorption

279

mechanisms of these adsorbents. Pore filling is the most common adsorption mechanism for

280

porous adsorbents, along with several other interactions (i.e., electrostatic attraction, hydrogen

281

bonding, surface precipitation, π-π interaction, n-π interaction, and cation exchange) that depend

282

on the characteristics of the porous adsorbents and adsorbates. Many studies have opted for a 24 h

283

equilibrating time for isotherm adsorption experiments, without providing evidence that

284

adsorption in their system attained true equilibrium within the appropriate timeframe.

M AN U

SC

RI PT

275

Therefore, it is necessary to provide the specified timeframe in which the adsorption

286

process approaches a true equilibrium. Figures 5a–c represent the influence of contact time on the

287

adsorption processes of various contaminants onto porous adsorbents (Hung and Lin, 2006;

288

Masson et al.; 2016; Kim et al., 2016). Clearly, the amount of contaminant adsorbed onto the

289

porous adsorbents continued to increase after 24 h. For example, Kim et al. (2016) speculated on

290

the process of methyl violet (MV) adsorption onto synthesized granular mesoporous carbon

291

(GMC) and commercial granular activated carbon (GAC). The results of their kinetic study

292

demonstrated that MV adsorption reached an equilibrium state within 20 d for GMC-polyvinyl

293

alcohol and 26 d for GMC-carboxymethyl cellulose, while adsorption onto GAC had not reached

294

equilibrium even after 75 d. Furthermore, a phenomenon involving desorption from the adsorbent

295

surface after reaching adsorption equilibrium (Figure 5d) has been reported in the literature

296

(Sharma et al., 2013, Halder; Khan et al., 2016). Therefore, studies of the adsorption kinetics play

AC C

EP

TE D

285

13

ACCEPTED MANUSCRIPT

297

an important role in identifying the required equilibration time and the optimal contact time and for

298

an adsorption process.

299

Figure 5 In aqueous-phase adsorption, two kinetic reaction models (i.e., pseudo-first-order and

301

pseudo-second-order equations), the Elovich or Roginsky–Zeldovich model, and the intraparticle

302

diffusion model have been widely applied to mathematically describe the intrinsic kinetic

303

adsorption constant. Furthermore, most mistakes found and discussed in the literature correspond

304

to the use of kinetic models. Therefore, the following sections (5.3–5.6) include a thorough

305

discussion of these models.

306

5.3. Pseudo-first-order (PFO) equation

307

5.3.1. Derivation of the PFO equation

M AN U

SC

RI PT

300

Recently, Ho published a series of critical papers (Ho and McKay, 1998b; Ho, 2005, Ho,

309

2014a; Ho, 2015; Ho, 2016e) discussing two main points. First, Ho (Ho and McKay, 1998b)

310

claimed to offer the correct reference style for citing Lagergren’s paper, with the original

311

translated from “Lagergren, S. (1898), Zur theorie der sogenannten adsorption gelöster stoffe.

312

Kungliga Svenska Vetenskapsakademiens Handlingar, Band 24, No. 4, 1–39” into English as

313

“Lagergren, S., About the theory of so-called adsorption of soluble substances, Kungliga Svenska

314

Vetenskapsakademiens. Handlingar, Band 24, No. 4, 1898, pp. 1–39.” Second, Ho stated that

315

Lagergren’s first-order rate equation has been called pseudo-first-order since 1998 (Ho, 2005) in

316

order to distinguish kinetic equations based on the adsorption capacity of solids from those based

317

on the concentration of a solution. However, the responses of many authors (Özacar; 2005, Kumar,

318

2006d; Utomo et al., 2010) state that Lagergren’s first-order rate expression had been referred to as

319

a pseudo-first-order equation before Ho’s publication in 1998. For example, in 1990, Sharma et al.

AC C

EP

TE D

308

14

ACCEPTED MANUSCRIPT

320

(1990) studied the removal of Ni(II) from an aquatic environment by wollastonite. In the section

321

on adsorption kinetics, they stated: “the adsorption of Ni (II) on wollastonite follows the first order

322

adsorption rate expression of Lagergren (1898).” The great contributions of Ho toward correcting the citation style format for the PFO

324

equation cannot be denied. However, the PFO equation was not originally proposed or initially

325

expressed by Ho. In the scientific community, the authors who first propose a theoretical model

326

should always be given credit for their contributions via correct citation (Kumar and Rattanaphani,

327

2007; Wu, 2009; Utomo et al., 2010). Therefore, the original work by Lagergren (1898), who first

328

presented the first-order-rate equation for the adsorption of oxalic acid and malonic acid onto

329

charcoal, should be cited. The PFO equation can be expressed correctly in nonlinear (Eq. 7) and

330

linear forms (Eq. 8 or Eq. 9).

log(qe − qt ) = log(qe ) −

ln(qe − qt ) = ln(qe ) − k1t

k1 t 2.303

TE D

qt = qe (1 − e −k1t )

M AN U

SC

RI PT

323

(7) (8)

(9)

where qe and qt are the amounts of adsorbate uptake per mass of adsorbent at equilibrium and at

332

any time t (min), respectively; and k1 (1/min) is the rate constant of the PFO equation.

333

5.3.2. Problems in the application of the PFO equation

AC C

334

EP

331

It cannot be denied that the PFO equation is widely applied in adsorption kinetics.

335

However, several problems related to the application of the PFO equation have been discussed in

336

the literature. First, many incorrect expressions for its linear forms have been published. Various

337

incorrect forms are shown in Eq. 10 (Alagumuthu et al., 2010), Eq. 11 (Fu et al., 2007; Jing et al.,

338

2014; Nekouei and Nekouei, 2017), Eq. 12 (Ho, 2014e; Ho, 2014a; Jing et al., 2014; Ho, 2016e; 15

ACCEPTED MANUSCRIPT

Ho, 2016b; Ho, 2016c; Ho, 2016a; Ho, 2016d; Ho, 2017), Eq. 13 (Ho, 2015), Eqs. 14 and 15 (Jing

340

et al., 2014), Eq. 16 (Kumar, 2006a), and Eq. 17 (Kumar and Porkodi, 2007b), which should be

341

avoided in future studies. As a result of these incorrect expressions, recalculation of the PFO

342

equation parameters was necessary to obtain the correct parameters (Wu et al., 2017). log(qe ) − k1 t 2.303

log(qe − qt ) =

log(qe ) − k1 2.303t

(10)

(11)

SC

log(qe − qt ) =

RI PT

339

log(qe − qt ) = log(qe )

Ci kt = 1 Ct 2.303

log(1 −

C i − Ct kt = 1 Ci − Ce 2.303

qe k1 = qe − qt 2.303t

EP

log

k1t 2.303

TE D

log

M AN U

1 k1 1 = + qt qet qe

log(qe − qt ) = log(qe ) −

k1 2.303

(12)

(13)

(14)

(15)

(16)

(17)

where Ci (mg/L), Ct (mg/L), and Ce (mg/L) are the concentrations of adsorbate at the initial time (t

344

= 0), any time t, and equilibrium, respectively.

345

AC C

343

The second problem is the two unknown parameters (qe and k1). Furthermore, in most

346

cases, the PFO equation is only appropriate for the initial 20 to 30 min of contact time, not for the

347

whole range (Ho and McKay, 1998a). Thus, plots of Eq. 8 or Eq. 9 are only linear over

348

approximately the first 30 min; beyond this initial period, the experimental and theoretical data

16

ACCEPTED MANUSCRIPT

will not fit adequately (McKay et al., 1999). Another important problem is the selection of an

350

appropriate qe value. Notably, the adjusted qe value cannot be lower than the maximum measured

351

value of qt. This is because mathematical errors associated with taking the logarithm of a negative

352

number will occur when using Eq. 8 or Eq. 9 (Plazinski et al., 2009). Finally, the qe value

353

calculated using the PFO equation is not equal to the qe value obtained from experiments, which

354

further indicates the inability of the PFO equation to fit kinetic adsorption data (Febrianto et al.,

355

2009). The presence of a boundary layer or external resistance controlling the beginning of the

356

sorption process was argued to be responsible for this discrepancy (McKay et al., 1999).

SC

RI PT

349

Therefore, two methods have been recommended to obtain an accurate estimation of the

358

kinetic parameters in the PFO equation: (1) a trial and error process to obtain the optimal qe value

359

(Ho and McKay, 1998a), or (2) application of a nonlinear optimization technique. A simple guide

360

to using the nonlinear method is introduced in Section 9.

361

5.4. Pseudo-second-order (PSO) equation

362

5.4.1. Derivation of the PSO equation

TE D

M AN U

357

363

dn = K (no − n)2 dt

EP

removal of heavy metals from water using natural zeolites, as follows: (18)

AC C

364

In 1984, Blanchard et al. (1984) initially proposed a second-order rate equation for the

Integration of Eq. 18 gives: 1 − α = Kt no − n

(19)

After rearranging, Eq. 19 becomes: n=

Ktno + αno −1 Kt + α

(20)

17

ACCEPTED MANUSCRIPT

365

where n is the amount of M2+ fixed or the amount of NH4+ released at each instant; no is an

366

exchange capacity; and K is a rate constant. Plotting 1/(n0 − n) as a function of time (t) must produce a straight line, the slope of which

368

gives the rate constant K and the intercept leads to the exchange capacity. Considering the

369

boundary condition of n = 0 for t = 0, it follows that α = 1/n, and Eq. 20 can be rearranged to: (21)

If qt = no, qe = n, and k2 = K, Eq. 21 becomes: qt =

qe2k2t 1+ k2qet

M AN U

370

n2 + Kt 1 + Knt

SC

no =

RI PT

367

(22)

371

where qe (mg/g) and qt (mg/g) are the amount of adsorbate adsorbed at equilibrium and at any t

372

(min), respectively; and k2 (g/mg × min) is the rate constant of the PSO equation. The nonlinear form of the PSO equation (Eq. 22) has been applied elsewhere (Ho, 1995;

374

Ho; Wase et al., 1996). Recently, a series of comments made by Ho implied that the derivation and

375

first application of the PSO equation belongs to his work. Any application of the PSO equation in

376

the field of adsorption studies without citation of the original work or with a misquotation (i.e., the

377

citation of a secondary reference) might suffer from his comments. However, in reaction to Ho’s

378

comments, a series of notes were published by other scholars (Kumar, 2006d; Kumar, 2006c;

379

Kumar and Fávere, 2006; Kumar and Guha, 2006; Kumar and Rattanaphani, 2007) to clarify the

380

correct citation of the PSO equation. Thus, the nonlinear form of the PSO equation for solid-liquid

381

adsorption systems was not originally reported by Ho (1995) or first applied by Ho et al. (1996),

382

but initially proposed by Blanchard et al. (1984).

383

AC C

EP

TE D

373

In the literature, the PSO equation can be correctly expressed in four linear forms (Eqs. 23–

18

ACCEPTED MANUSCRIPT

384

26), and the initial adsorption rate (h) can be described by Eq. 29.

385

Expression for the linear form (Type 1) with a plot of t/qt versus t: t 1 1 = ( )t + qt qe k2qe2

RI PT

386

(23)

Expression of the linear form (Type 2) with a plot of 1/qt versus 1/t: 1 1 1 1 =( 2) + qt k2qe t qe

qt = − ( 388

SC

Expression of the linear form (Type 3) with a plot of qt versus qt/t: 1 qt ) + qe k 2 qe t

M AN U

387

(24)

Expression of the linear form (Type 4) with a plot of qt/t versus qt:

qt = −(k 2 qe )qt + k 2 qe2 t

rearranged to:

t 1 + 2 k 2 qe qe

(27)

Subsequently, when t approaches 0, Eqs. 23 and 27 can be expressed as Eq. 28. Setting h = qt/t, the initial adsorption rate (h; mg/g × min) can then be determined by Eq. 29.

AC C

392

1

EP

qt =

391

(26)

The initial adsorption rate (h) was proposed by Ho et al. (1996). First, Eq. 22 can be

TE D

389 390

(25)

t 1 = qt k 2 qe2



qt = k 2 qe2 t

(28)

h = k 2 qe2

(29)

393

Clearly, only the linear form (Type 1, Eq. 23) might have been reported by Ho et al., (1996)

394

for the adsorption of dyes from waste streams by peat (Özacar, 2005; Kumar, 2006d; Kumar and

19

ACCEPTED MANUSCRIPT

Fávere, 2006; Wu, 2009). Kumar and Fávere (2006) hypothesized that, in adsorption research, the

396

Langmuir equation (Langmuir, 1918) has been the most widely used isotherm to describe

397

adsorption processes under equilibrium conditions. In the literature, four linear types of the

398

Langmuir isotherm have been reported. Irrespective of the linear expression reported, these

399

equations have been widely called the Langmuir isotherm. It is not correct to claim that the

400

Langmuir model was derived or proposed by other authors who transformed the original Langmuir

401

expression into a new linear expression (Kumar and Guha, 2006). In scientific literature, the names

402

of scientists (i.e., Langmuir or Freundlich) have always been presented, even after modification of

403

their equations (Arica, 2003). Similarly, Utomo et al. (2010) suggested that it might not be

404

necessary to cite papers that deal with mathematical modification and different mathematical

405

expressions of the models (i.e., PFO and PSO, and Langmuir). However, from the scientific point

406

of view, the authors who first proposed a theoretical model should be credited because of their

407

great contributions to the advancement of science. Therefore, the original work by Blanchard et al.

408

(1984) should be cited for the expression of the PSO equation (Özacar, 2005; Kumar, 2006d;

409

Kumar and Fávere, 2006; Plazinski et al., 2009; Wu, 2009; Utomo et al., 2010).

410

5.4.2. Problems in the application of the PSO equation

EP

TE D

M AN U

SC

RI PT

395

First, many incorrect expressions of the linear form of PSO equation are found in the

412

literature, as shown in Eq. 30 (Fu et al., 2007; Kumar and Porkodi, 2007b), Eq. 31 (Ho, 2014f),

413

Eqs. 32 and 33 (Ho, 2014g), Eq. 34 (Ho, 2013), Eq. 35 (Lin and Wang, 2009), Eq. 36 (Ho, 2014c),

414

Eq. 37 (Ho, 2014d), Eq. 38 (Ho, 2014), and Eq. 39 (Ho, 2015).

AC C

411

t 1 1 = + 2 qt k 2 qe qe t

(30)

20

ACCEPTED MANUSCRIPT

(31)

t 1 t = + qt k 2 q e q e

(32)

RI PT

t 1 1 = + qt 2k 2 qe

t 1 1 = ( qe2 ) + t qt k2 qe

(33)

t 1 1 = + t 2 qt k2 − qe qe

SC

(34)

M AN U

t 1 1 = ( )t 2 qt (k2 qe ) qe t 1 t = k2 qe2 + qt 2 qe

t 1 1 = + 2 qt k2 qe qe

(36)

(37)

(38)

(39)

EP

t kt t = 22 + qt k2 qe qe

TE D

1 1 1 = + 2 qt k2 qe qe

(35)

Second, based on the best fit of kinetic experimental data using the PSO model, many

416

authors have drawn controversial conclusions. For example, “the adsorption process is the

417

chemisorption, which involve valence forces through sharing electrons between adsorbate and

418

adsorbent”, “the best fit of experimental kinetic data in pseudo-second-order kinetics suggests the

419

chemisorption, which may involve valency forces through sharing of electrons between dye anion

420

and adsorbent”, and other similar conclusions. The problem with such conclusions has been

421

discussed by several researchers (Kumar, 2006a; Lima et al., 2015; Lima et al., 2016), who

AC C

415

21

ACCEPTED MANUSCRIPT

422

claimed that adsorption mechanisms cannot be directly assigned based on observing simple kinetic

423

experiments or by fitting kinetic models (i.e., the PFO and PSO models). Adsorption mechanisms can only be established by (1) using several analytical techniques

425

(i.e., FTIR, SEM, nitrogen adsorption-desorption isotherms, Raman spectroscopy, TGA/DTA,

426

DSC,

427

titration, and solution calorimetry), and (2) having a good sense of the chemical nature of the

428

adsorbate and adsorbent, adsorbent’s surface, and chemical or physical interactions between the

429

adsorbent and adsorbate (Volesky, 2007; Lima et al., 2015; Lima et al., 2016). The use of

430

analytical techniques together with adsorptive thermodynamic data (i.e., changes in enthalpy and

431

entropy) and activation and adsorption energies, are necessary to confirm whether the adsorption

432

of contaminants in aqueous solution is a chemical or physical process (Lima et al., 2015; Tran et

433

al., 2016).

Si and

13

C solid-state NMR, XRD, XPS, pHPZC, pHIEP, CHN element analysis, Boehm

M AN U

SC

29

RI PT

424

Finally, in most cases, experimental data for adsorption kinetics are well fitted by the linear

435

form of the PSO equation. However, caution should be taken in drawing conclusions based on

436

fitting to the linear form of the PSO equation. For example, Figure 6 presents a plot of t/qt versus t.

437

Clearly, the R2 values for the linear form of the PSO equation are very high (R2 > 0.99); however,

438

the corresponding R2 values for the nonlinear form of the PSO equation are significantly lower (R2

439

= 0.53–0.68). This finding means that the adsorption process of methylene green 5 (MG5) onto

440

synthesized activated carbon and commercial activated carbon (porous materials) is not

441

adequately described by the PSO equation. In fact, Tran et al. (2017a) concluded that the process

442

of MG5 adsorption onto these porous materials involved π-π interactions and pore filling, and was

443

not related to chemical adsorption. Obviously, the nonlinear method can be applied to obtain

444

kinetic model parameters that are more accurate than those obtained using the linear method.

AC C

EP

TE D

434

22

ACCEPTED MANUSCRIPT

Figure 6

445 446

5.5. Elovich equation

447

5.5.1. Derivation of the Elovich equation An empirical equation was firstly proposed by Roginsky and Zeldovich (1934) for the

449

adsorption of carbon monoxide onto manganese dioxide. However, this equation is now generally

450

known as the Elovich equation and has been extensively applied to chemisorption data

451

(McLintock, 1967). This equation can be expressed mathematically as follows:

SC

RI PT

448

(40)

M AN U

dqt = α exp(− βqt ) dt 452

where qe and qt are the amounts of adsorbate uptake per mass of adsorbent at equilibrium and at

453

any time t (min), respectively; α (mg/g × min) is the initial rate constant because dqt/dt → α when

454

qt → 0; and β (mg/g) is the desorption constant during any one experiment.

become Eq. 41 (nonlinear). qt =

460

ln(1 + αβ t )

(41)

EP

linear form (Eq. 42) can be obtained: qt =

459

β

To simplify the Elovich equation, Chien and Clayton (1980) assumed αβt >> 1. Thus, a

457 458

1

AC C

456

By applying the boundary conditions of qt = 0 at t = 0, the integrated form of Eq. 40 will

TE D

455

1

β

ln(αβ ) +

1

β

ln(t )

(42)

and a plot of qt versus lnt should give a linear relationship with a slope of (1/β) and an

intercept of (1/β)ln(αβ).

23

ACCEPTED MANUSCRIPT

461

5.5.2. Incorrect expression of the Elovich equation Errors in the presentation of the nonlinear form of the Elovich model are common in the

463

literature. These incorrect equations (Eqs. 43–45) can be found in several previous papers (Ho and

464

McKay, 1998a; Lin and Wang, 2009; Alagumuthu et al., 2010; Ben Ali et al., 2016).

RI PT

462

qt = β ln(αβ ) + ln(t ) 1

β

ln(t )

(44)

SC

qt = β ln(αβ ) +

(43)

465 466 467

5.6. Intra-particle diffusion model

M AN U

qt = β ln(αβ ) + β ln(t )

(45)

The linearized transformation of the intra-particle diffusion model (Weber and Morris, 1963) is presented as follows: qt = k p t + C

(46)

where kp (mg/g × min) is the rate constant of the intra-particle diffusion model and C (mg/g) is a

469

constant associated with the thickness of the boundary layer, where a higher value of C

470

corresponds to a greater effect on the limiting boundary layer.

TE D

468

Generally, although the PSO model can adequately describe adsorption kinetic

472

experimental data, this model does not reveal the adsorption mechanisms. Similarly, the Elovich

473

equation seems to describe various reaction mechanisms, such as bulk diffusion and surface

474

diffusion. In contrast, the intra-particle diffusion model can be useful for identifying the reaction

475

pathways and adsorption mechanisms and predicting the rate-controlling step. In a solid-liquid

476

sorption process, adsorbate transfer is often characterized by film diffusion (also known as

477

external diffusion), surface diffusion, and pore diffusion, or combined surface and pore diffusion.

478

In short, if a plot of qt against t0.5 is linear and passes through the origin, the adsorption is entirely 24

AC C

EP

471

ACCEPTED MANUSCRIPT

479

governed by intra-particle diffusion. In contrast, if the intra-particle diffusion plot gives multiple

480

linear regions, then the adsorption process is controlled by a multistep mechanism. Four steps associated with transport processes during adsorption by porous adsorbents

482

were originally proposed in reference (Walter, 1984). The first stage is transport in the solution

483

phase (known as “bulk transport”; occurs quickly), which can occur instantaneously after the

484

adsorbent is transferred into the adsorbate solution; therefore, it does not control engineering

485

design. In most cases, this stage occurs too rapidly and its contribution is considered negligible.

486

The second stage is “film diffusion” (occurs slowly). In this stage, the adsorbate molecules are

487

transported from the bulk liquid phase to the adsorbent’s external surface through a hydrodynamic

488

boundary layer or film. The third stage involves diffusion of the adsorbate molecules from the

489

exterior of the adsorbent into the pores of the adsorbent, along pore-wall surfaces, or both (known

490

as “intraparticle diffusion”; occurs slowly). The last stage, adsorptive attachment, often occurs

491

very quickly; therefore, it is also not significant for design. These four steps are summarized in

492

Figure 7 (Weber and Smith, 1987).

TE D

M AN U

SC

RI PT

481

Figure 7

493

6. Adsorption isotherms

495

6.1. Adsorption equilibrium

EP

494

Milonjić (2009a) and Milonjić (2010) suggested that the amount of adsorbate adsorbed onto

497

an adsorbent depended on the equilibrium concentrations of metal ions, equilibrium solution pH,

498

and temperature. An adsorption isotherm should be given as the equilibrium adsorbed amounts

499

versus the equilibrium ion concentrations for a constant equilibrium solution pH and temperature.

500

Thus, environmental parameters in the sorption system (especially solution pH) must be carefully

501

controlled at the given value over the entire contact period until the sorption equilibrium is

AC C

496

25

ACCEPTED MANUSCRIPT

reached. These comments were in good agreement with the suggestion that metal-ion adsorption is

503

pH dependent (Volesky, 2003; Tien, 2007; Volesky, 2007; Tien, 2008). However, the pH effect

504

has invariably been examined in terms of the initial pH of the aqueous solution. Not a single case

505

has reported data demonstrating the changes in the solution pH during the course of an adsorption

506

experiment.

RI PT

502

In the scientific literature, adsorption equilibrium studies in batch experiments have been

508

conducted in different ways as follows: (1) solution pH values were not mentioned in the

509

manuscript, (2) initial pH values were neither measured nor adjusted and final pH values were

510

measured, (3) initial and final pH values were measured, (4) initial pH values were adjusted and

511

final pH values were not mentioned, (5) initial pH values were adjusted and final pH values were

512

measured, (6) initial and final pH values were adjusted, (7) pH was controlled several times during

513

the course of the reaction, (8) a constant pH was maintained throughout the equilibration time by

514

adding acidic or alkaline solutions, and (9) experiments were conducted at a constant pH using a

515

buffered solution (Smičiklas et al., 2009).

TE D

M AN U

SC

507

Kumar (2006a) recommended that presenting a plot of qe versus Ce for the complete

517

adsorption isotherm in an adsorption study plays a key role in identifying the regions (e.g., Henry,

518

Freundlich, Langmuir, and BET) in which the experimental data relating to adsorption equilibrium

519

are actually located. In contrast, plots of the Freundlich, Langmuir, or other adsorption isotherm

520

models do not help to identify these regions. Figures 8a and 8b provide typical complete

521

adsorption isotherms that should be presented in the investigation of liquid-phase adsorption.

522

Kumar and Sivanesan (2006) also recommended that using equilibrium data covering the

523

complete isotherm was the best way to obtain the parameters in isotherm expressions; equilibrium

524

data with a partial isotherm was not sufficient, instead equilibrium data that covers the complete

AC C

EP

516

26

ACCEPTED MANUSCRIPT

525

isotherm was required. Furthermore, an inconsistent definition of the adsorption isotherm plot in the literature was

527

recognized by Li and Pan (2007); the plot of qe against Ce/Wo (with Wo being a liquid/solid ratio) is

528

invisible with the traditional definition of the complete adsorption isotherm (the plot of qe versus

529

Ce).

RI PT

526

Recently, Senthil Kumar et al. (2014) presented a plot of %removal versus C0 in an

531

adsorption equilibrium study instead of presenting the complete adsorption isotherm (Figure 8a).

532

However, Hai (2017) made a critical note that it is clear that at every initial concentration of

533

methylene blue (MB; 50, 100, 150, 200, 250 mg/L), the removal percentage of MB by sulfuric

534

acid-treated orange peel (STOP) at 30 °C was always overwhelmingly higher than the

535

corresponding values at 40, 50, and 60 °C. This result revealed that increasing the temperature

536

from 30 to 60 °C resulted in a decrease in the adsorption capacity (qe) of MB dye onto STOP

537

(Figure 8a). However, the maximum adsorption capacities of STOP from the Langmuir model

538

(Qomax) exhibited the following order: 60 °C (83.333 mg/g) > 50 °C (71.428 mg/g) > 40 °C (62.5

539

mg/g) > 30 °C (50 mg/g). This logic problem might be caused by the adsorption process of MB

540

onto STOP not reaching equilibrium at the different temperatures (see adsorption isotherms in

541

Figure 8a) owing to the low initial MB concentrations (50–250 mg/L) used at a high solid/liquid

542

ratio of 4.0 g/L (Hai, 2017). Therefore, in adsorption equilibrium studies, presenting plots of

543

%removal versus Co, %removal versus Ce, or qe versus Co should be avoided, and a plot of qe

544

versus Ce for the complete adsorption isotherm must be presented.

M AN U

TE D

EP

AC C

545

SC

530

Figure 8b provides a typical example of the comparison of adsorption performance under

546

the same operating conditions (i.e., temperature and solution pH). Clearly, “adsorbent B” exhibits

547

a higher affinity than “adsorbent A” at low initial adsorbate concentrations (Co) and low

27

ACCEPTED MANUSCRIPT

548

equilibrium adsorbate concentrations (Ce), whereas the opposite order (“adsorbent A” >

549

“adsorbent B”) was found at higher Co and Ce values (Volesky, 2007). It was recently identified that many researchers have made questionable conclusions when

551

studying the effects of agitation rate or shaking speed on the adsorption capacity (qe) or maximum

552

adsorption capacity (Qomax) of an adsorbent toward an adsorbate. For example, Chu et al. (2005)

553

concluded that the adsorption capacity of vitamin E onto silica varies as a function of the agitation

554

rate (120–180 rpm). The Qomax values exhibited the following order: 43.71 mg/g at 180 rpm >

555

24.20 mg/g at 160 rpm > 24.20 mg/g at 140 rpm > 17.62 mg/g at 120 rpm. However, Choong and

556

Chuah (2005) made a critical remark that the agitation rate only impacts on the speed at which a

557

system reaches equilibrium, but not on the equilibrium itself. This means that when the adsorption

558

process obtains a true equilibrium (sufficient contact time between adsorbent and adsorbate), the

559

Qomax values at different agitation rates will be insignificantly different. This is because the

560

decrease of the boundary layer resistance a higher agitation rate will facilitate the transport of the

561

adsorbate to the adsorbent.

SC

M AN U

TE D

562

RI PT

550

Figure 8

Although adsorption isotherms can insignificantly contribute to elucidating adsorption

564

mechanisms, they are less helpful in this regard than data on adsorption kinetics and

565

thermodynamics (Volesky, 2003; Tien, 2007; Tien, 2008). However, collecting adsorption

566

isotherms is a useful strategy to both describe the relationship between the adsorbate concentration

567

in solution (liquid phase) and the adsorbent (solid phase) at a constant temperature and design

568

adsorption systems (Tran et al., 2016).

AC C

EP

563

569

A wide variety of models of adsorption isotherms have been applied in the literature. These

570

models can be classified as follows: (1) irreversible isotherms and one-parameter isotherms (i.e.,

28

ACCEPTED MANUSCRIPT

Henry isotherm); (2) two-parameter isotherms (i.e., Langmuir, Freundlich, Dubinin–

572

Radushkevich, Temkin, Flory–Huggins, and Hill); (3) three-parameter isotherms (i.e., Redlich–

573

Peterson, Sips, Toth, Koble–Corrigan, Khan, Fritz–Schluender, Vieth–Sladek, and Radke–

574

Prausnitz); and (4) more than three-parameter isotherms (i.e., Weber–van Vliet, Fritz–Schlunder,

575

and Baudu) (Kumar and Porkodi, 2006; Hamdaoui and Naffrechoux, 2007; Foo and Hameed,

576

2010). Of these, the Langmuir and Freundlich models are the most commonly used, followed by

577

the Dubinin–Radushkevich and Redlich–Peterson models, because of the usefulness of their

578

model parameters, their simplicity, and their easy interpretability.

579

6.2. Langmuir equation

580

6.2.1. Derivation of the Langmuir equation

M AN U

SC

RI PT

571

The theoretical Langmuir equation (Langmuir, 1918), which was originally applied to the

582

adsorption of gases on a solid surface, was developed using the following assumptions: (1) a fixed

583

number of accessible sites are available on the adsorbent surface and all active sites have the same

584

energy; (2) adsorption is reversible; (3) once an adsorbate occupies a site, no further adsorption

585

can occur on that site; and (4) there is no interaction between adsorbate species. The nonlinear

586

form of the Langmuir model is described in Eq. 47 and its four linearized forms are shown in Eqs.

587

48–51.

588

o Qmax + K L Ce 1 + K LCe

AC C

qe =

EP

TE D

581

(47)

Hanes–Woolf linearization of the Langmuir model (Type 1): Ce 1 1 = ( o )Ce + o qe Qmax Qmax K L

589

(48)

Lineweaver–Burk linearization of the Langmuir model (Type 2):

29

ACCEPTED MANUSCRIPT

1 1 1 1 =( o ) + o qe Qmax K L Ce Qmax Eadie–Hoffsiee linearization of the Langmuir model (Type 3): qe = ( 591

− 1 qe o ) + Qmax K L Ce

(50)

Scatchard linearization of the Langmuir model (Type 4): qe o = − K L qe + Qmax KL Ce

RI PT

590

(49)

SC

(51)

where Qomax (mg/g) is the maximum saturated monolayer adsorption capacity of an adsorbent, Ce

593

(mg/L) is the adsorbate concentration at equilibrium, qe (mg/g) is the amount of adsorbate uptake

594

at equilibrium, and KL (L/mg) is a constant related to the affinity between an adsorbent and

595

adsorbate. For a good adsorbent, a high theoretical adsorption capacity Qomax and a steep initial

596

sorption isotherm slope (i.e., high KL) are generally desirable (Kratochvil and Volesky, 1998;

597

Wang, 2009).

TE D

M AN U

592

The limitations of using the four linear forms of the Langmuir model have been

599

highlighted by Bolster and Hornberger (2007). The transformation of data for linearization can

600

result in modifications of error structure, introduction of error into the independent variable, and

601

alteration of the weight placed on each data point, which often leads to differences in the fitted

602

parameter values between linear and nonlinear versions of the Langmuir model. For the Type 1

603

linearization, because Ce and Ce/qe are not independent, the correlation between Ce and Ce/qe is

604

overestimated, i.e., this equation may provide good fits to data that do not conform to the

605

Langmuir model. For the Type 2 linearization, transformation leads to clustering of data points

606

near the origin and is extremely sensitive to variability at low values of qe (high values of 1/qe).

607

For the Type 3 linearization, the abscissa is not error free, and qe/Ce and qe are not independent. 30

AC C

EP

598

ACCEPTED MANUSCRIPT

Thus, the correlation between qe/Ce and qe is underestimated, i.e., this equation may provide poor

609

fits to data that do conform to the Langmuir model. Finally, the variables qe and qe/Ce are also

610

not independent in the Type 4 linearization. In this case, the correlation between qe and qe/ Ce is

611

underestimated, i.e., this equation may provide a poor fit to data that do conform to the Langmuir

612

model.

613

6.2.2. Derivation of the separation factor

RI PT

608

If the experimental data are adequately described by the Langmuir model, it is essential to

615

calculate the separation factor. Hall et al. (1966) originally proposed that the essential

616

characteristics of the Langmuir isotherm model can be expressed in terms of a dimensionless

617

constant called the separation factor or equilibrium parameter RL, which is defined as follows: RL =

M AN U

SC

614

1 1 + K LCo

(52)

where RL is a constant separation factor (dimensionless) of the solid-liquid adsorption system, KL

619

is the Langmuir equilibrium constant, and Co (mg/L) is the initial adsorbate concentration.

TE D

618

Notably, the isotherm shape can be predicted using the separation factor (RL) (Weber and

621

Chakravorti, 1974) and the Freundlich exponent n (Worch, 2012). Table 1 summarizes the various

622

isotherm shapes.

624 625

Table 1

AC C

623

EP

620

6.2.3. Mistakes concerning the Langmuir model and the separation factor Wan, et al. (2014) published a corrigendum because of a misapplication of the nonlinear

626

form of the Langmuir equation (Eq. 53). A similar mistake in presenting the nonlinear form of the

627

Langmuir equation (Eq. 54) was observed by Kumar and Porkodi, (2007b).

31

ACCEPTED MANUSCRIPT

(53)

Ce 1 1 = + o o qe K LQmax Qmax Ce

(54)

RI PT

Ce C K = oe + o L qe Qmax Qmax

Furthermore, an unexpected typographical error was made in reference (Egirani, 2004). As

629

a result, the sentence “a plot of Ce/qe versus Ce should indicate a straight line of slope 1/KLQomax

630

and an intercept of 1/Qomax” was corrected to “a plot of Ce/qe versus Ce should indicate a straight

631

line of slope 1/Qomax and an intercept of 1/ KLQomax” by Ho (2004a).

SC

628

The incorrect expression and inaccurate calculation of RL were recognized and discussed

633

by Kumar (2006a) and Kumar (2006b). First, the separation factor RL is not a Langmuir constant,

634

and is only employed for isotherms that obey the Langmuir model. Second, RL values will vary

635

with the initial adsorbate concentration, irrespective of the shape of the isotherm, and RL values at

636

different initial adsorbate concentrations cannot be calculated directly from the Langmuir

637

isotherm, but can be calculated from the Langmuir equilibrium constant KL. These

638

misunderstandings and mistakes are summarized in Table 2.

TE D

M AN U

632

Furthermore, a mistake related to the definition of RL compared with the original definition

639

was stated by Fu et al. (2007). The incorrect presentation is shown in Eq. 55.

641

1 o 1 + K LQmax Co

(55)

AC C

RL =

EP

640

Recently, Hai (2017) commented that “the authors’ conclusion about separation factor

642

(RL) was not valid. It is noted that there is no information or value of RL mentioned in the paper’s

643

results and discussion. It is impossible for the authors to conclude that the values of RL were found

644

to be between 0 and 1, indicating the favorable adsorption of MB dye onto STOP.”

645

Clearly, the separation factor is related to the Langmuir equation in the adsorption 32

ACCEPTED MANUSCRIPT

isotherm. Nevertheless, Goswami et al. (2016) recently described an incorrect relationship

647

between the separation factor and adsorption kinetic models (i.e., PFO, PSO, and intra-particle

648

diffusion) at different initial cadmium concentrations. This is a fundamental misunderstanding that

649

should be avoided.

651

6.3. Freundlich equation

652

6.3.1. Derivation of the Freundlich equation

SC

Table 2

650

RI PT

646

The Freundlich equation is one of the earliest empirical equations used to describe

654

equilibrium data and adsorption characteristics for a heterogeneous surface (Freundlich, 1906).

655

Unlike the Langmuir equation, the Freundlich isotherm can describe neither the (arithmetic)

656

linearity range at very low concentrations nor the saturation effect at very high concentrations.

657

Hence, the Freundlich isotherm does not describe the saturation behavior of an adsorbent (Wang,

658

2009). The nonlinear and linear forms of the Freundlich equation can be expressed as shown in

659

Eqs. 56 and 57, respectively.

qe = K F Cen

TE D

M AN U

653

log qe = n log Ce + log K F

(56)

EP

(57)

where qe (mg/g) is the amount of adsorbate uptake at equilibrium, Ce (mg/L) is the adsorbate

661

concentration at equilibrium, KF (mg/g)/(mg/L)n is the Freundlich constant, and n (dimensionless)

662

is the Freundlich intensity parameter, which indicates the magnitude of the adsorption driving

663

force or the surface heterogeneity.

AC C

660

33

ACCEPTED MANUSCRIPT

6.3.2. Mistakes concerning the Freundlich equation

665 666

Several investigators have identified an incorrect expression (Eq. 58) of the linear form of the Freundlich equation (Kumar and Porkodi, 2007b; Alagumuthu et al., 2010).

log qe =

1 + log K F n log Ce

RI PT

664

(58)

Furthermore, several discussions that are inconsistent with the original concepts of the

668

Freundlich theory have been identified in the literature (Lu, 2008; Hai, 2017). Hai (2017)

669

commented that the exponent n of the Freundlich equation must be dimensionless. Therefore, the

670

original authors incorrectly reported the units of the exponent n (g/L); and there is no validation for

671

the statement of the significance of the n value for the adsorption of MB dye onto STOP is as

672

follows: when n = 1 the adsorption is linear, when n < 1 the adsorption is a chemical process, and

673

when n > 1 the adsorption is a physical process. According to the Freundlich theory, the adsorption

674

isotherm becomes linear when n = 1, favorable when n < 1, and unfavorable when n > 1. Similarly,

675

Harsha et al. (2015) also made a basic mistake when they reported the units of the exponent n as

676

g/L.

TE D

M AN U

SC

667

Lu (2008) highlighted that in many cases KF is not the maximum adsorption capacity and

678

the KF value is only equal to the maximum adsorption capacity when n approaches infinity.

679

Although KF is not defined as the maximum adsorption capacity (Qomax), KF values and Qomax

680

values should be of the same order. Such a mistake has been discussed by Tran (2017), where the

681

reported KF values followed the order: 4.471 at 303 K > 3.439 at 313 and 323 K > 2.877 at 333 K,

682

while the Qomax values were 30.12 mg/g at 313 K > 21.367 mg/g at 323 K > 21.276 mg/g at 333 K

683

> 20.492 mg/g at 303 K. This mistake might be attributable to a miscalculation.

684

AC C

EP

677

Notably, Kumar (2006a) highlighted that values for the exponent n (Eq. 59) in the range of

34

ACCEPTED MANUSCRIPT

0–10 suggest favorable adsorption. If the experimental equilibrium data do not lie in the Henry

686

region, then the linear regression method will just check the hypothesis instead of verifying the

687

theory behind the Freundlich model; therefore, it is not possible in practice to obtain n > 10

688

(Kumar, 2006a). A corresponding error was also found in a published paper (Huang et al., 2014),

689

in which the n values reported for the removal of aniline and Cr(VI)/aniline + Cr(VI) by an

690

activated carbon/chitosan composite were 10.20 and 11.34, respectively.

RI PT

685

SC

qe = K F Ce1 / n 6.4. Redlich–Peterson equation

692

6.4.1. Derivation of the Redlich–Peterson equation

M AN U

691

(59)

A three-parameter Redlich–Peterson isotherm was proposed upon considering the

694

limitations of the Freundlich and Langmuir isotherms (Redlich and Peterson, 1959). This model

695

incorporates the features of the Freundlich and Langmuir models and might be applicable for

696

demonstrating adsorption equilibrium over a wide range of adsorbate concentrations. The

697

nonlinear form of this empirical model is given as follows:

K RPCe 1 + a RP Ceg

(60)

EP

qe =

TE D

693

where KRP (L/g) and aRP (mg/L)−g are the Redlich–Peterson constants and g (dimensionless) is an

699

exponent whose value must lie between 0 and 1.

700

AC C

698

Eq. 60 becomes a linear isotherm (Henry’s law equation) at low surface coverage (g = 0),

701

reduces to the Langmuir isotherm when g = 1, and transforms into the Freundlich isotherm when

702

KRP and aKP >> 1 and g = 1. Therefore, if the g value is outside the range of 0–1, the data is not

703

adequately explained by the Redlich–Person equation. Accurate calculation of g can be helpful in

704

expatiation of where the isotherm presentations, neither in the Freundlich or Langmuir (Kumar

35

ACCEPTED MANUSCRIPT

705

and Sivanesan, 2006).

706

6.4.2 Mistakes concerning the Redlich–Peterson equation A series of comments related to mistakes in calculating the exponent g have been reported

708

elsewhere (Inbaraj, 2006; Kumar and Porkodi, 2007a; Kumar et al., 2007; Kumar et al., 2007;

709

Kumar and Porkodi, 2008). These mistakes involved values of the exponent g that were outside the

710

range of 0–1 (Table 3). It is not easy to predict the parameters of the Redlich–Peterson equation

711

because there are three unknown parameters. Although a new linear form of the Redlich–Peterson

712

equation was published by Wu et al. (2010), the nonlinear form might be more appropriate for

713

accurately calculating the parameters of adsorption models with more than two unknown

714

parameters, such as the Redlich–Peterson model.

M AN U

SC

RI PT

707

Table 3

715

6.5. Dubinin–Radushkevich equation

717

6.5.1. Derivation of the Dubinin–Radushkevich equation

structure of an adsorbent (Dubinin and Radushkevich, 1947), and is expressed as follows:

qe = q DR e − K RDε

2

EP

719

The Dubinin–Radushkevich equation was developed to account for the effect of the porous

(61)

The linear form of the Dubinin–Radushkevich equation is:

AC C

718

TE D

716

ln q e = − K DR ε

ε = RT ln(1 +

2

+ ln q DR

(62)

1 ) Ce

(63)

By inserting Eq. 63 into Eq. 61, Eq. 64 can be obtained:

36

ACCEPTED MANUSCRIPT

2

2

2

ln q e = − K DR R T ln (1 +

1 ) + ln q DR Ce

(64)

The parameters qDR and KDR in Eq. 64 can be obtained as follows: (1) a plot of lnqe against

721

ln2(1 + 1/Ce) has a slope = −KDRR2T2 and an intercept = lnqRD, and the E value can be obtained

722

using Eq. 65, and (2) a plot of lnqe against R2T2ln2(1 + 1/Ce) has a slope = −KDR and an intercept =

723

lnqRD, and the E value can be obtained using Eq. 66. Notably, the E values obtained from Eqs. 65

724

and 66 are the same (Tran et al., 2016).

E=

1 1 = 2 K DR − 2slope

SC

1 RT = 2 K DR − 2 slope

M AN U

E=

RI PT

720

(65)

(66)

where qRD (mg/g) is the adsorption capacity, KRD (mol2/kJ2) is a constant related to the sorption

726

energy, ɛ is the Polanyi potential, E (kJ/mol) is the mean adsorption energy, R is the gas constant, T

727

is the temperature in Kelvin, and qe and Ce are obtained from Eq. 1.

728

6.5.2. Mistakes concerning the Dubinin–Radushkevich equation

731

equation (Cavas, 2008; Fu et al., 2008). These erroneous forms are shown in Eqs. 67–70.

EP

730

Some authors have highlighted erroneous expressions of the Dubinin–Radushkevich

Incorrect nonlinear form of the Dubinin–Radushkevich equation:

AC C

729

TE D

725

2

q e = q DR exp( − K DR ( RT (ln 1 + 1 / C e ) ) 732

Incorrect linear form of the Dubinin–Radushkevich equation:

ln q e = −2 K DR RT ln(1 + 733

(67)

1 ) + ln q DR Ce

(68)

Incorrect expressions for the Polanyi potential (ɛ):

37

ACCEPTED MANUSCRIPT

ε = RTe1/ C

(69)

1 ) Ce

(70)

e

RI PT

ε = RT ln(

Additionally, the magnitude of E may give useful information about the type of adsorption

735

process (physical or chemical) (Fu et al., 2008; Tran et al., 2016). Fu et al. (2008) commented that

736

the E values for the removal of phenolic compounds by organophilic bentonite can be categorized

737

as corresponding to ion exchange (E = 8–16 kJ/mol), whereas the original authors made an

738

inconsistent conclusion that the adsorption of two phenolic compounds could be regarded as

739

physical adsorption.

M AN U

SC

734

Recently, Tran (2017) made the following comment on the conclusions of the original

741

paper: “the obtained value of E in this work varies from 32.596 kJ/mol to 40.572 kJ/mol for studied

742

temperature which is higher than the adsorption energy values previously enumerated.” However,

743

these adsorption energy (E) values might not have been calculated correctly. Assuming the β

744

values at 303 K (β = 1.93 × 10−8 mol2/kJ2), 313 K (2.99 × 10−8), 323 K (2.00 × 10−8), and 333 K

745

(2.35 × 10−8) reported by these authors are correct, the recalculated adsorption energy (E) values

746

are 5090 kJ/mol at 303 K, 4089 kJ/mol at 313 K, 5000 kJ/mol at 323 K, and 4613 kJ/mol at 333 K.

747

Thus, the E values in this study ranged from 4089 to 5091 kJ/mol, which is impossible for a

748

heavy-metal biosorption process.

749

7. Adsorption thermodynamics

750

7.1. Principles of adsorption thermodynamics

AC C

EP

TE D

740

751

Thermodynamic studies are an indispensable component of predicting adsorption

752

mechanisms (e.g., physical and chemical). The key distinctions between physical and chemical

38

ACCEPTED MANUSCRIPT

753

adsorption were summarized in detail in a recent report (Tran et al., 2016). The thermodynamic

754

parameters can be computed according to the laws of thermodynamics using the following

755

equations: ∆G o = − RT ln K C

RI PT

756

(71)

The relationship between ∆G° and ∆H° and ∆S° is described as follows: ∆Go = ∆H 0 − T∆S 0

The well-known van’t Hoff equation is obtained by substituting Eq. 71 into Eq. 72 − ∆H o 1 ∆S o x + R T R

(73)

M AN U

ln KC =

SC

757

(72)

758

where R is the universal gas constant (8.3144 J/(mol × K)) and T is the absolute temperature in

759

Kelvin.

The Gibbs energy change (∆G°) is directly calculated from Eq. 71, whereas the enthalpy

761

change (∆H°) and entropy change (∆S°) are determined from the slope and intercept, respectively,

762

of a plot of lnKC versus 1/T (Eq. 73). The common units of ∆G°, the gas constant, and temperature

763

are J/mol, J/(mol × K), and K, respectively; therefore, the equilibrium constant KC in Eq. 73 must

764

be dimensionless (Milonjić, 2007; Milonjić, 2009b; Canzano et al., 2012; Dawood and Sen, 2012;

765

Zhou et al., 2012; Zhou and Zhou, 2014; Tran et al., 2016; Ghosal and Gupta, 2017; Hai, 2017;

766

Rahmani-Sani et al., 2017; Zhou, 2017).

EP

AC C

767

TE D

760

Eq. 72 describes the relationship of ∆G° with ∆H° and ∆S°, and it is useful for checking the

768

logic of calculated values for these thermodynamic parameters. Unfortunately, associated

769

mistakes have been found in the literature. Fu et al. (2008) substituted the values of ∆H° and ∆S°

770

into Eq. 72 to obtain ∆G°, but the obtained ∆G° values were significantly different from the ∆G°

771

values obtained using Eq. 71. Fu et al. (2008) highlighted that this mistake might result from

39

ACCEPTED MANUSCRIPT

772

incorrect application of the method for calculating the equilibrium constant. Clearly, accurate estimation of thermodynamic parameters is directly dependent on

774

accurate determination of the equilibrium constant between two phases (KC; dimensionless). In the

775

literature, the thermodynamic parameters can be calculated from KC values derived from

776

adsorption-isotherm constants (i.e., Langmuir, Freundlich, Frumkin, Flory–Huggins, and Henry)

777

or the partition coefficient (Liu, 2009; Doke and Khan, 2013; Tran et al., 2016). However, to

778

obtain an appropriate calculation method, several factors need to be considered thoroughly. First,

779

the equilibrium constant (KC) must be dimensionless. Second, the linear regression coefficient (R2)

780

of the van’t Hoff equation must be high. Third, the adsorbate in solution must have a high or low

781

concentration. In addition, the temperatures used for calculating the thermodynamic parameters

782

must have units of Kelvin (K), not degrees Celsius (°C). The incorrect use of temperature units led

783

to the publication of an erratum by Ho (2004a), with accurate recalculations.

M AN U

SC

RI PT

773

As the equilibrium constant KC can be derived using a number of approaches, considerable

785

variation in the thermodynamic parameters can be obtained; therefore, the most appropriate

786

approach should be determined. Here, we thoroughly discuss the derivation of the equilibrium

787

constant KC from adsorption isotherms (i.e., Langmuir, Freundlich, and Henry constants) and the

788

partition coefficient because of their popularity in the literature.

789

7.2. Equilibrium constant derived from the Langmuir constant (KL)

EP

AC C

790

TE D

784

The Langmuir equation was initially derived from a kinetic study and then subsequently

791

from a thermodynamic study. The derivation of the Langmuir equation from the thermodynamic

792

perspective can be described as follows (Crittenden et al., 2012; Tran et al., 2016). The

793

relationship between vacant surface sites on the surface of an adsorbent (Sv; mmol/m2), adsorbate

794

species in solution (A; mmol), and adsorbate species bound to surface sites (SA; mmol/m2) can be

40

ACCEPTED MANUSCRIPT

795

described by the following reaction:

SV + A ↔ SA

(74)

On the basis of the Langmuir expression, it is assumed that the reaction has a constant

797

Gibbs energy change (∆G°; J/mol) for all sites, so the thermodynamic equilibrium constant (KC;

798

dimensionless) can be expressed as:

RI PT

796

− ∆G o

⇔ KC = e

∆G = − RT ln K C

RT

(75)

SC

o

Another Langmuir assumption is that each site is capable of binding at most one adsorbate

800

molecule (monolayer). According to the equilibrium condition, the thermodynamic equilibrium

801

constant may be written as:

SA KC = =e SV [ A]

M AN U

799

− ∆G o

RT

(76)

where [A] is the concentration of adsorbate A in solution at equilibrium (mg/L), R is the universal

803

gas constant (8.314 J/mol × K), and T is absolute temperature K (273 + °C).

TE D

802

The critical problem with Eq. 76 is that there are two unknown parameters (i.e., SV and

805

[A]). However, this problem can be solved if the total sites available or monolayer coverage (ST;

806

mol/m2) are fixed:

807

SA + SA KC [ A]

AC C

ST = SV + SA =

EP

804



SA =

K C [ A]ST 1 + K C [ A]

(77)

The expression of SA in units of mmol/m2 is not useful in mass balance, and units of mass

808

of adsorbate adsorbed per unit mass of adsorbent (qe; mg/g) are much more useful. The mass

809

loading (qe) can be obtained in units of mg/g (Eq. 1) by multiplying both sides of Eq. 71 by the

810

surface area (SBET; m2/g) of the adsorbent and the molecular weight of the adsorbate (Mw; g/mol).

41

ACCEPTED MANUSCRIPT

811

Thus, Eq. 77 becomes:

K [ A]ST S BET M w SA.S BET .M w = C 1 + K C [ A]



o Qmax K L Ce qe = 1 + K L Ce

(78)

where qe (mg/g) = SA × SBET × Mw is described by Eq. 1; Ce (mg/L) = [A] is the concentration of

813

adsorbate in solution at equilibrium; Qomax (mg/g) = ST × SBET × Mw is the maximum monolayer

814

adsorptive capacity of the adsorbent when the surface sites are saturated with adsorbate; the

815

relationship between KL and KC is described in the following discussion.

SC

RI PT

812

However, the main problem is that the Langmuir constant KL is dimensional with common

817

units of L/mmol or L/mg, while the equilibrium constant KC is dimensionless (without units).

818

Thus, the direct application of KL (L/mmol or L/mg) in the calculation of thermodynamic

819

parameters produces incorrect results, as discussed by many authors (Milonjić, 2007; Milonjić,

820

2009b; Canzano et al., 2012; Dawood and Sen, 2012; Zhou et al., 2012; Zhou and Zhou, 2014;

821

Anastopoulos and Kyzas, 2016; Tran et al., 2016; Ghosal and Gupta, 2017; Hai, 2017,

822

Rahmani-Sani et al., 2017). To solve this unit problem, several methods have been recommended

823

(Milonjić, 2007; Zhou and Zhou, 2014; Tran et al., 2016). Depending on the units of KL, the

824

equilibrium constant KC can easily be obtained as a dimensionless constant.

EP

TE D

M AN U

816

When an adsorption study is conducted in aqueous solution and KL has units of L/mmol, KC

826

can be easily obtained as a dimensionless parameter by multiplying KL by 55.5 and then by 1,000

827

(Eq. 79). This method was originally proposed by Milonjić (2007) and then developed by Zhou

828

and Zhou (2014) and Tran et al. (2016). The values of the parameters ∆G°, ∆H°, and ∆S° can be

829

calculated using the following equations:

AC C

825

KC = 55.5×1,000× KL

(79)

42

ACCEPTED MANUSCRIPT

∆Go = −RT ln(55.5×1,000× KL ) ln(55 .5 × 1,000 × K L ) =

(80)

− ∆H o 1 ∆S o × + R T R

(81)

where the factor 55.5 is the number of moles of pure water per liter (1,000 g/L divided by 18

831

g/mol) and the term 55.5 × 1,000 × KL is dimensionless.

RI PT

830

In the case of KL with units expressed in L/mg, Milonjić (2007 and 2009b) stated that KL

833

could be obtained as a dimensionless parameter by multiplying KL by 106 (Eq. 82). However, Zhou

834

and Zhou (2014) recommended that KC could be obtained as a dimensionless parameter by

835

multiplying KL by the molecular weight of the adsorbate (Mw; g/mol), by 1000, and then by 55.5

836

(Eq. 85).

M AN U

SC

832

KC = 106 K L ∆G o = − RT ln(106 K L )

TE D

− ∆H o 1 ∆S o ln(10 K L ) = x + R T R 6

(82) (83) (84)

where the factor 106 is the solution density (assuming the density of pure water is 1.0 g/mL) and

838

the term 106 × K is dimensionless.

EP

837

K C = M w × 55.5 ×1,000 × K L

AC C

(85)

∆Go = −RT ln(Mw × 55.5×1,000× KL )

(86)

− ∆ H o 1 ∆S o ln( M w × 55.5 ×1,000 × K L ) = × + R T R

(87)

839

where the factor 55.5 is the number of moles of pure water per liter and the term Mw × 55.5 × 1,000

840

× KL is dimensionless.

43

ACCEPTED MANUSCRIPT

In a recent paper, Tran et al. (2016) applied Eqs. 79–87 to compare thermodynamic

842

parameters calculated from the Langmuir constants derived from four linear forms of the

843

Langmuir model. They concluded that (1) both the recommendations of Milonjić (2007 and

844

2009b) and Zhou and Zhou (2014) provide the same methods for calculating the thermodynamic

845

parameters, and (2) there is some deviation between these methods for the results from Types 1, 3,

846

and 4; only Type 2 produces identical results.

847

7.3. Equilibrium constant derived from the Freundlich constant (KF)

SC

RI PT

841

The Freundlich equation is consistent with the thermodynamics of heterogeneous

849

adsorption (Crittenden, et al., 2012). The Freundlich constant KF can be obtained as a

850

dimensionless value using Eq. 88 (Ghosal and Gupta, 2015; Tran et al., 2016).

M AN U

848

K ρ 106 (1− ) KC = F ( ) n 1000 ρ 1

K ρ 106 (1− ) ∆G o = − RT ln( F ( ) n ) 1000 ρ

K F ρ 10 6 (1− n ) − ∆ H o 1 ∆S o x + ( ) )= R T R 1000 ρ 1

ln(

(89)

(90)

where ρ is the density of pure water (assumed as ~1.0 g/mL).

EP

851

TE D

1

(88)

It is noted that the units of KF rely on the units used for the liquid-phase concentration (C)

853

and solid-phase concentration (q). Units of mg/L or mmol/L for C and mg/g or mmol/g for q are

854

used most frequently to demonstrate adsorption from water solutions. The differing units of KF can

855

be reciprocally converted using Eq. 91 (Worch, 2012).

AC C

852

K F = K F' ( M w )1−1 / n 856

(91)

where the units of the Freundlich constant KF are (mg/g)/(mg/L)1/n, with Co and Ce (mg/L) and qe

44

ACCEPTED MANUSCRIPT

857

(mg/g); the units of KF′ are (mmol/g)/(mmol/L)1/n, with Co and Ce (mmol/L) and qe (mmol/g); and

858

Mw is the molecular weight of the adsorbate. The exponent 1/n is not affected by the units of the

859

liquid-phase concentration (C) and solid-phase concentration (q) Based on a comparison of the thermodynamic parameters calculated from the Freundlich

861

constant (dimensionless) and Langmuir constant (dimensionless), Tran et al., (2016) concluded

862

that the signs and magnitudes of ∆G°, ∆H°, and ∆S° calculated from KF were consistent with those

863

calculated from KL. Of course, the experimental data in the adsorption isotherms must be fitted

864

well by the Langmuir and Freundlich models, and the R2 value of the van’t Hoff equation must be

865

higher than 0.90 (Tran et al., 2016). It is also evident that, to some extent, the application of KF to

866

estimate the thermodynamic parameters should be approached with caution.

867

7.4. Equilibrium constant derived from the partition coefficient (Kp)

M AN U

SC

RI PT

860

868

reported by Biggar and Cheung (1973). The equilibrium constant can be defined as follows:

Kp =

as γ s C s = ae γ eCe

TE D

869

Changes in the thermodynamic partition coefficient with changes in temperature were first

(92)

where as is the activity of the adsorbate adsorbed onto the adsorbent, ae is the activity of the

871

adsorbate in solution at equilibrium, γs is the activity coefficient of the adsorbate adsorbed onto the

872

adsorbent, γe is the activity coefficient of the adsorbate in solution at equilibrium, Cs is the

873

concentration of adsorbate adsorbed onto the adsorbent at equilibrium (mg/L), and Ce is the

874

concentration of adsorbate in solution at equilibrium (mg/L). Cs is defined by the mass balance of

875

adsorbate that disappears from the solution, which should appear on the adsorbent.

876 877

AC C

EP

870

When the concentration of adsorbate in the solution approaches zero, which results in Cs → 0 and Ce → 0, the activity of coefficient γ approaches unity, and Eq. 92 can be written as:

45

ACCEPTED MANUSCRIPT

Cs

Lim C Cs →0

e

=

as = Kp ae

(93)

Kp values can be obtained plotting ln(Cs/Ce) versus Cs and extrapolating Cs to zero. If a

879

straight line fits the data with a high regression coefficient (R2) and its intersection with the vertical

880

axis provides the value of Kp, the partition coefficient will be in unison with the equilibrium

881

constant. The ∆G° value can be directly calculated using Eq. 71, while the values of ∆H° and ∆S°

882

are determined from the slope and intercept, respectively, of Eq. 73.

SC

RI PT

878

This method has been widely applied in the literature. A summary of the applications of

884

this method can found in the critical review in reference (Doke and Khan, 2013). However,

885

incorrect applications of Biggar and Cheung’s definition have also been identified in the literature.

886

In a typical experiment, Senthil Kumar et al. (2014) made a mistake in the application of Kp during

887

calculation of the thermodynamic parameters for the adsorption of MB onto STOP (Table 4). Hai

888

(2017) highlighted that it is impossible to obtain an equilibrium constant, Gibbs energy change

889

(∆G°), enthalpy change (∆H°), and entropy change (∆S°) for the adsorption process at every initial

890

MB concentration, as shown in Table 4. Assuming that the experimental data reported by Senthil

891

Kuma et al. (2014) are correct, Hai (2017) recalculated the thermodynamic parameters based on

892

the method of Biggar and Cheung (1973), and the corrected thermodynamic parameters are also

893

listed in Table 4. Likewise, Lima et al. (2015) also determined the value of the adsorption

894

equilibrium constant (KC) from the value of the best fit nonlinear isotherm equilibrium model

895

instead of using the equilibrium constant (Kp) calculated from the initial and equilibrium

896

concentrations of the adsorbate.

AC C

EP

TE D

M AN U

883

897

Table 4

898

Importantly, according to the United States Environmental Protection Agency (USEPA,

46

ACCEPTED MANUSCRIPT

1999), the partition coefficient is only appropriate for calculating the thermodynamic parameters

900

if the initial adsorbate concentration is low. Analogous conclusions can be found in the literature

901

(Liu, 2009; Salvestrini et al., 2014; Lima et al., 2015; Tran et al., 2016; Hai, 2017). It has been

902

highlighted that the partition coefficient is only equal to the thermodynamic equilibrium constant

903

when the adsorbed concentrations in the solutions are very low. In this situation, the partition

904

coefficient can be used to calculate the thermodynamic parameters. Recently, Tran et al. (2016),

905

who compared the thermodynamic parameters calculated from the Langmuir constant and the

906

partition coefficient, concluded that although the signs of the calculated ∆G°, ∆H°, and ∆S°

907

parameters were the same, the partition coefficient might not be appropriate for calculating the

908

thermodynamic parameters owing to the low regression coefficient of the van’t Hoff equation (R2

909

= 0.53).

910

7.5. Equilibrium constant derived from the distribution coefficient (Kd)

M AN U

SC

RI PT

899

In 1987, Khan and Singh (1987) proposed a method for calculating thermodynamic

912

parameters (Eq. 94) based on the method of Biggar and Cheung (1973). The distribution

913

coefficient (Kd) can be defined as follows:

(94)

In this method, Kd values are obtained by plotting ln(qe/Ce) against Ce and extrapolating to

AC C

914

qe Ce

EP

Kd =

TE D

911

915

zero Ce. If a straight line fits the data with a high regression coefficient (R2), then its intersection

916

with the vertical axis provides the value of Kd. Fundamentally, Eq. 94 can be derived from the

917

Freundlich and Langmuir equations. If n = 1, the Freundlich equation (Eq. 56) will become

918

Henry’s equation for a linearized isotherm (Eq. 95), and the Freundlich constant (KF) will be equal

919

to the Henry constant (KH). For very low concentrations of adsorbate (KL and Ce << 1), the

47

ACCEPTED MANUSCRIPT

920

Langmuir equation (Eq. 47) will become Henry’s equation (Eq. 95), and the Henry constant (KH)

921

will become QomaxKL. The Henry constant (KH) is also referred to as the distribution coefficient

922

(Kd).

qe = K H Ce

RI PT

(95)

This method has been widely applied to calculating thermodynamic parameters in the

924

literature. Unfortunately, the units of Kd in Eq. 94 are L/g as Ce has units of mg/L and qe has units

925

of mg/g. Therefore, it is impossible to directly employ Kd as the equilibrium constant for the

926

calculation of ∆G°, ∆H°, and ∆S°. This problem associated with the units has been discussed by

927

several scientists (Patrickios and Yamasaki, 1997; Milonjić, 2007; Canzano et al., 2012; Dawood

928

and Sen, 2012; Tran et al., 2016; Zhou, 2017). Canzano, Iovino et al. (2012) suggested that the

929

partition coefficient Kd (L/g) could be converted to KC (dimensionless) by multiplying Kd by a

930

factor of 1,000, as originally proposed by Milonjić (2007 and 2009b).

M AN U

SC

923

Similar to the partition coefficient, the distribution coefficient may be appropriate for

932

estimating the thermodynamic parameters when the initial concentration of the adsorbate is low

933

(Salvestrini et al., 2014; Tran et al., 2016). Tran et al. (2016) concluded that the signs of ∆G°, ∆H°,

934

and ∆S° estimated from the distribution coefficient and the Langmuir constant are the same.

935

However, the units of Kd (L/g) are different from those of KL (L/mg), and the adsorption

936

equilibrium data do not fit Henry’s adsorption isotherm. As a result, application of the distribution

937

coefficient was not appropriate for calculating the thermodynamic parameters in this study; the

938

low regression coefficient of the van’t Hoff equation (R2 = 0.42) is further evidence of the

939

inappropriateness of this approach.

AC C

EP

TE D

931

940

In summary, to determine accurate thermodynamic parameters, several points should be

941

thoroughly considered. First, the equilibrium constant (KC) must be dimensionless. Second, if the 48

ACCEPTED MANUSCRIPT

concentrations of adsorbate used to obtain the adsorption isotherm are very low (diluted solutions),

943

the partition (Kp) or distribution coefficients (Kd) are appropriate for calculation of the

944

thermodynamic parameters; otherwise, the Langmuir (KL) or Freundlich (KF) constant models may

945

be more appropriate. Third, the plot of lnKC against 1/T corresponding to the van’t Hoff equation

946

must be linear with a high regression coefficient (R2). Additionally, it is necessary to consider

947

whether the adsorption process can reach equilibrium and consideration of the adsorption isotherm

948

shapes and the adsorption model fit are also recommended. Thus, presentation of the complete

949

adsorption isotherm (plot of qe versus Ce) is strongly recommended. Finally, the determined

950

thermodynamic parameters must have a logical relationship with the experimental data.

M AN U

SC

RI PT

942

Notably, numerous recent studies have reported methods of calculating thermodynamic

952

parameters from the partition coefficient (Kp) or distribution coefficient (Kd) that differed from the

953

traditional concepts defined in Sections 7.4 and 7.5, respectively. A typical mistake is present in

954

the work of Nekouei et al. (2016). They studied the adsorption thermodynamics at only one initial

955

adsorbate concentration and different temperatures instead of using a range of adsorbate

956

concentrations (i.e., from 100 mg/L to 1,000 mg/L) and different temperatures (i.e., 10 °C, 20 °C,

957

30 °C, and 50 °C). This mistake led to invalid conclusions regarding the adsorption

958

thermodynamics.

959

8. Others mistakes

960

8.1. Inconsistent data points in experimental data and model fitting

EP

AC C

961

TE D

951

Significant differences between the number of data points in an experiment and those used

962

for model fitting should be avoided. For example, Azizian (2008a) made the following critical

963

comments about the paper of Karadag et al. (2007): there were six experimental data points

964

regarding the adsorption kinetics (plots of qt versus t), but only five experimental data points were

49

ACCEPTED MANUSCRIPT

used for fitting of the experimental data to the PSO model (plot of t/qt versus t). This omission

966

resulted in inaccurate calculations and questionable conclusions. The R2 values of the PSO model

967

were very high (0.992–1.000), but the difference between the experimental (qe,exp) and calculated

968

(qe,cal) values was more than 100% (for example qe,exp = 6.720 mg/g and qe,cal = 15.526 mg/g).

969

Clearly, the qe values calculated from the experimental data and the model were significantly

970

different, but the authors (Karadag et al., 2007) concluded that “The pseudo-second-order kinetic

971

model agrees very well with the dynamic behavior for the adsorption of dyes RR239 and RB5 onto

972

CTAB-zeolite under several different dye concentrations, temperatures, and pH values”; the

973

validity of this conclusion is very questionable.

M AN U

SC

RI PT

965

Lima et al. (2015) observed a change in kinetic parameters when five experimental points

975

were deleted to achieve a good R2 value. They concluded that the strategy of deleting some data

976

points to improve the R2 value caused serious errors in the values of qe and k1 in the PFO equation

977

or k2 in the PSO equation.

TE D

974

Azizian (2008) pointed out that in reference (Karadag et al., 2007), a comparison was made

979

between the results of fitting four experimental equilibrium data points to the Langmuir model and

980

fitting five experimental equilibrium data points to the Freundlich model. Therefore, it is

981

impossible for the authors (Karadag et al., 2007) to conclude that “the Freundlich model exhibited

982

a slightly better fit than the Langmuir model.” To compare the fits of two different models, the

983

same experimental data (including the same number of experimental data points) should be used

984

(Azizian, 2008a). Likewise, Harsha et al. (2015) examined adsorption isotherms using seven

985

experimental data points. However, in the plots for calculating the Langmuir, Freundlich, and

986

Dubinin–Radushkevich models, only three experimental points were used. It is impossible to

AC C

EP

978

50

ACCEPTED MANUSCRIPT

987

obtain the Langmuir, Freundlich, and Dubinin–Radushkevich parameters from only three

988

experimental points. A similar mistake was found in reference (Goswami et al., 2016).

989

8.2. Oxidation state of chromium Many investigators (Aoyama, 2003; Park et al., 2006a; Park et al., 2006b; Park et al.,

991

2006c) have commented that mistakes in analyzing chromium species in aqueous solution have

992

resulted in the incorrect elucidation of hexavalent chromium adsorption, attributing Cr(VI)

993

removal from an aqueous solution to “anionic adsorption”. Generally, when Cr(VI) comes into

994

contact with organic substances or reducing agents, especially in an acidic medium, it is easily or

995

spontaneously reduced to Cr(III), as Cr(VI) has a high redox potential (above +1.3 V under

996

standard conditions). Therefore, there is a strong possibility that the mechanism of Cr(VI) removal

997

by biomaterials or biomaterial-based activated carbon is not “anionic adsorption” but

998

“adsorption-coupled reduction”. Various authors (Aoyama, 2003; Park et al., 2006a; Park et al.,

999

2006b; Park et al., 2006c) also suggested that in Cr(VI) adsorption studies, it is necessary to

1000

analyze both the Cr(VI) and total Cr in aqueous solution, as well as determining the oxidation state

1001

of chromium bound to the biomaterial or activated carbon.

TE D

M AN U

SC

RI PT

990

The authors (Aoyama, 2003; Park et al., 2006a; Park et al., 2006b; Park et al., 2006c)

1003

pointed out that in most papers, only Cr(VI) in aqueous solutions was analyzed by the standard

1004

colorimetrical method; the pink colored complex formed from between 1,5-diphenycarbazide and

1005

Cr(VI) in an acidic solution can be spectrophotometrically analyzed at 540 nm. Similarly, the total

1006

Cr in aqueous solutions has been analyzed using atomic absorption spectroscopy (AAS) or

1007

inductively coupled plasma-atomic emission spectroscopy (ICP-AES), which cannot distinguish

1008

Cr(VI) from the total Cr. Additionally, Aoyama (2003) underlined that the total Cr present in a

1009

solution can be determined by oxidizing any Cr(III) formed with KMnO4, followed by the same

AC C

EP

1002

51

ACCEPTED MANUSCRIPT

procedures described for the determination of Cr(VI). The adsorption density can be calculated

1011

from the difference between the initial concentration of Cr(VI) and the final concentration of total

1012

Cr. In contrast, the Cr(III) concentration is calculated from the difference between the final total Cr

1013

and final Cr(VI) concentrations.

RI PT

1010

Several techniques have been applied in the literature to ascertain the oxidation state of

1015

chromium on an adsorbent, such as X-ray photoelectron spectroscopy (XPS) and X-ray absorption

1016

spectroscopy (XAS). For example, Dambies et al. (2001) investigated the characterization of

1017

Cr(VI) interactions with chitosan using XPS. On the basis of their XPS data, they concluded that

1018

Cr(VI) was entirely reduced to Cr(III) on glutaraldehyde crosslinked beads, while only 60% of

1019

Cr(VI) was found in its reduced form on raw beads. Additionally, Cr(VI) binding and reduction to

1020

Cr(III) by agricultural byproducts of Avena monida (oat) biomass was characterized by

1021

Gardea-Torresdey et al. (2000). Their XAS studies further corroborated that, although Cr(VI)

1022

could bind to oat biomass, it was easily reduced to Cr(III) by positively charged functional groups,

1023

and Cr(III) was subsequently adsorbed by available carboxyl groups.

TE D

M AN U

SC

1014

A study on using fermentation waste for detoxification of Cr(VI) contaminated aqueous

1025

solutions was published by Park et al. (2008). In this study, the colorimetric method was used in

1026

combination with excess potassium permanganate to analyze both the Cr(VI) and total Cr in

1027

aqueous solution, and XPS was used to ascertain the oxidation state of the chromium bound on the

1028

biomass. Park et al. (2008) concluded that the mechanism of Cr(VI) removal by the fermentation

1029

waste involved reduction of Cr(VI) to Cr(III) (redox reaction). As a result, kinetic and equilibrium

1030

models based on the “reduction” mechanism were used to describe the behavior of Cr(VI) and total

1031

Cr in aqueous solution.

1032

AC C

EP

1024

Two main mechanisms for removal of Cr(VI) from aqueous solution by nonliving

52

ACCEPTED MANUSCRIPT

biomass were proposed by Park et al. (2005). In the first mechanism (direct reduction), Cr(VI)

1034

was directly reduced to Cr(III) in the aqueous phase by contact with electron-donor groups of the

1035

biomass, i.e., groups with lower reduction potentials than that of Cr(VI) (+1.3 V). However, the

1036

second mechanism comprised three steps: (1) binding of anionic Cr(VI) ion species to positively

1037

charged groups on the biomass surface, (2) reduction of Cr(VI) to Cr(III) by adjacent

1038

electron-donor groups, and (3) release of the Cr(III) ions into the aqueous phase, owing to

1039

electronic repulsion between the positively charged groups and the Cr(III) ions, or complexation

1040

of Cr(III) with adjacent groups capable of binding chromium. If there are a small number of

1041

electron-donor groups in the biomass or a low concentration of protons in the aqueous phase,

1042

chromium bound on the biomass can remain in the hexavalent state. Therefore, the extent to

1043

which mechanisms I and II operate depends on the nature of the biosorption system, including

1044

solution pH, temperature, functional groups on the biomass, and biomass and Cr(VI)

1045

concentrations (Park et al., 2006b; Park et al., 2006a). These two main mechanisms are

1046

summarized in Figure 9.

Figure 9

1047

8.3. Incorrect labels

EP

1048

TE D

M AN U

SC

RI PT

1033

Incorrect labeling of compounds should also be avoided. A publication entitled “Selenite

1050

adsorption and desorption in main Chinese soils with their characteristics and physicochemical

1051

properties” included such a mistake, as pointed out by Goldberg (2016). The purpose of this work

1052

was to study the adsorption of Se(IV) (selenite), not Se(VI) (selenate); however, the paper reported

1053

that the authors prepared Se solutions using Na2SeO4 (hexavalent Se(VI)). This mistake was

1054

explained by the incorrect labeling of Na2SeO3 as Na2SeO4 in the Materials and Methods section.

AC C

1049

53

ACCEPTED MANUSCRIPT

Similarly, the author of the paper entitled “Adsorption sequence of toxic inorganic anions

1056

on a soil” prepared solutions from Na2HAsO4·7H2O as As(V) species. However, they mislabeled

1057

Na2HAsO4·7H2O as NaAsO2, a trivalent As(III) species (Goldberg, 2009).

1058

8.4. BET specific surface area of an adsorbent

RI PT

1055

Using the Brunauer–Emmett–Teller (BET) surface area of pahokee peat (0.88 m2/g) and

1060

assuming that the atrazine molecule is a sphere with a radius (r) of 4.16 × 10−10 m, Borisover and

1061

Graber (1997) estimated the maximum adsorption capacity of pahokee peat for atrazine. First, the

1062

cross-section area (A) of atrazine can be obtained as 3.14 × (4.16 × 10−10)2 = 54.4 × 10−20 m2;

1063

subsequently, the surface covered by 1 mol of atrazine (MA) can be calculated as (54.4 × 10−20 m2)

1064

× (6.023 × 1023/mol) = 3.3 × 105 m2/mol. Second, the maximum possible concentration of atrazine

1065

on the surface area (Qomax,covered; maximum adsorption capacity) was obtained by dividing the BET

1066

surface area (0.88 m2/g) by the surface covered by 1 mol of atrazine (3.3 × 105 m2/mol) to give 2.7

1067

× 10−6 mol/g. Finally, Qomax,covered was converted into weight units using (2.7 × 10−6 mol/g) ×

1068

(215.7 g/mol) to give 582 µg/g. This value (582 µg/g) is substantially higher than the maximum

1069

adsorption capacity obtained from the Langmuir equation (161 µg/g). Although these calculations

1070

are correct, these inconsistent values might originate from a misunderstanding by the authors

1071

(Borisover and Graber, 1997) who state that “N2 BET measured external surface area”, and thus

1072

attributed the value of 582 µg/g to the maximal possible atrazine concentration on the external

1073

surface of pahokee peat.

M AN U

TE D

EP

AC C

1074

SC

1059

According to the 2011 guide for the accelerated surface area and porosimetry (ASAP)

1075

system, the definition of the BET specific surface area (m2/g) includes both the t-Plot external

1076

surface area (m2/g) and the t-Plot micropore area (m2/g). Table 5 gives a typical example of correct

1077

and incorrect use of the Micromeritics report (experimental data measured using a Micromeritics

54

ACCEPTED MANUSCRIPT

ASAP 2020 sorptometer at 77 K). Although there were no problems found in the N2

1079

adsorption/desorption isotherms, the report results for BET analysis of the prepared activated

1080

carbon in Table 5 contain discrepancies, as the BET specific surface area value is less than the

1081

external surface area. This error resulted from calculation of the specific surface area (SBET) in the

1082

range 0.38 < p/p0 < 0.57, instead of 0.05 < p/p0 < 0.3.

Table 5

1083

RI PT

1078

Recently, Ben Ali et al. (2016) used the “iodine number” to determine the surface area of a

1085

biosorbent without any treatment (pomegranate peel; PGP) and concluded that “the specific

1086

surface area obtained is equal to 598.78 m2/g. Iodine number is generally used as an

1087

approximation for surface area and microporosity of active carbons with good precision”. Tran

1088

(2017) commented that there are two serious misconceptions in this work that need to be

1089

discussed.

M AN U

SC

1084

First, the iodine number of activated carbon is often determined following the internal

1091

method, ASTM D4607-14 (D4607-14 2014). According to the ASTM, the iodine number (mg/g;

1092

amount of iodine adsorbed (mg) by 1.0 g of activated carbon) is a relative indicator of porosity in

1093

activated carbon. Iodine molecules (≈0.27 nm) can be adsorbed into the micropores (pore width >

1094

1 nm) of porous materials. However, the iodine number does not necessarily provide a measure of

1095

the ability of activated carbon to adsorb other species. Although the iodine number may be used to

1096

approximate the surface area for several types of activated carbon, it must be realized that the

1097

relationship between surface area and iodine number cannot be generalized, as it varies with

1098

changes in the raw carbon material, processing conditions, and pore volume distribution.

1099

Biosorbents that have not undergone any treatment have never been defined as porous materials;

1100

thus, the porosity of a biosorbent is negligible.

AC C

EP

TE D

1090

55

ACCEPTED MANUSCRIPT

Therefore, the iodine number cannot be applied for the determination of porosity for any

1102

biosorbent. Moreover, it is impossible to obtain the very high iodine number reported for PGP

1103

(602 mg/g). The high iodine number for PGP can be attributed to the high sulfur content (0.89%)

1104

and prolonged contact time (>30 s) (Tran, 2017). According to the ASTM (D4607-14 2014), after

1105

transferring activated carbon (g) into hydrochloric acid solution (5 wt%), the mixture should be

1106

boiled gently for 30 s to remove any sulfur that may interfere with the test results. Therefore, it is

1107

clear that the iodine number is greatly affected by high sulfur contents. Notably, the contact time

1108

between activated carbon and iodine solution (0.1 N) strongly impacts the determined iodine

1109

number. The contact time recommended by the ASTM was approximately 30 s, while the contact

1110

time used by Ben Ali et al. was approximately 4 min. Recently, Tran et al. (2017c) prepared

1111

activated carbons from golden shower using different chemical activation methods with K2CO3

1112

(GSAC, GSBAC, and GSHAC). Their results indicated that the iodine numbers (mg/g) of the

1113

prepared activated carbons at a contact time of 30 s were significantly lower than those at 5 min;

1114

GSAC (2,604 mg/g < 2,883 mg/g), GSBAC (1,568 mg/g < 2,296 mg/g), GSHAC (2,695 mg/g <

1115

4,842 mg/g).

TE D

M AN U

SC

RI PT

1101

Second, the textural properties of an adsorbent (i.e., specific surface area, total pore

1117

volume, and micropore volume) can be obtained from conventional analysis of nitrogen

1118

adsorption-desorption isotherms, which are measured at 77 K using a sorptometer (e.g.,

1119

Micromeritics ASAP 2020). The BET method is widely employed to compute specific surface

1120

area (SBET) (Marsh and Reinoso, 2006). Therefore, it is impossible to measure the surface area of

1121

an adsorbent using the iodine number method. This mistaken assumption will lead to incorrect

1122

result, such as the extremely high surface area of PGP (598.78 m2/g). The BET surface areas of

1123

various biosorbents are approximately 40 m2/g for yellow passion-fruit shell, 20 m2/g for orange

AC C

EP

1116

56

ACCEPTED MANUSCRIPT

peel, 8.17 m2/g for wheat straw, 8.13 m2/g for sargassum, 4.01 m2/g for Moringa oleifera lamarck

1125

seed powder, 1.31 m2/g Spirogyra species, 1.21 m2/g for waste pomace from an olive oil factory,

1126

0.76 m2/g for soy meal shells, 0.48 m2/g for rubber tree leaves, and 0.48 m2/g for rice bran (Farooq

1127

et al., 2010; Tran et al., 2016). Thus, the discussions and conclusions regarding the surface area of

1128

PGP made by Ben Ali et al. (2016) are not valid, as they are inconsistent with the ASTM

1129

definitions and fundamental understanding of the porosities of materials.

1130

8.5. Maximum absorption wavelength in dye adsorption studies

SC

RI PT

1124

A comment on the method of using the maximum absorption wavelength (λmax) to calculate

1132

Congo red (CR) solution concentrations was published by Zhou et al. (2011). They pointed out that

1133

CR is very sensitive to pH and changes from red to blue owing to a π-π* transition of the azo group

1134

following protonation (Figure 10a). At lower pH values, CR is protonated and cationic CR shows

1135

two tautomeric forms: ammonium rich variety and azonium variety. Therefore, the maximum

1136

absorption wavelength used to calculate the concentration of CR in solution is strongly dependent

1137

on the solution pH (both initial and final) (Figure 10b); for example, λmax was 576 nm at pH 2.18–

1138

3.16, 567 nm at pH 3.86, and 496 nm at pH ≥ 4.71.

TE D

Figure 10

EP

1139

M AN U

1131

Tien (2007) and Tien (2008) highlighted that commercial dyes are often mixtures of active

1141

ingredients and filler materials are often not recognized. Therefore, the possibility that

1142

experiments may involve bi-solute adsorption is often not considered.

1143

8.6. Cπ-cation and π-π interactions

1144

8.6.1. Cπ-cation interactions

AC C

1140

1145

Recently, Morosanu et al. (2016) investigated the biosorption of lead ions onto rapeseed

1146

biomass that was used as a biosorbent without any treatments. On the basis of FTIR measurements

57

ACCEPTED MANUSCRIPT

before and after adsorption, they concluded that the disappearance of the peak at 1710 cm−1, which

1148

is characteristic of C=O stretching, indicated surface complexion through Cπ-cation interactions.

1149

Similarly, Medellin-Castillo et al. (2017) examined the biosorption of Pb2+ and Cd2+ onto a

1150

biosorbent derived from industrial chili seeds. They also found changes in the intensity of the

1151

absorbance peak at 1657 cm−1 in heavy-metal-loaded chili seeds, suggesting that the adsorption of

1152

Cd2+ and Pb2+ involved π-cation interactions between the aromatic rings of the lignin and the

1153

Cd2+ and Pb2+ cations in the solution.

SC

RI PT

1147

In general, Cπ-cation interactions are attributed to electrostatic interactions between the

1155

aromatic rings of basic carbonaceous materials (i.e., biochar, carbon foam, carbon nanotubes,

1156

activated carbon, and graphene) and metallic cations. For example, Tran et al. (2015) used

1157

orange-peel-derived biochar, which was produced through a carbonization process under

1158

limited-oxygen conditions (also known as pyrolysis) at a high temperature (>400 °C). Biochar is

1159

known as a carbon-enriched porous material, similar to activated carbon. Thus, biochar also

1160

possesses a graphitic structure (C=C bonds; π-electrons). As a result, Cπ-cation interactions played

1161

a primary role in the adsorption of heavy metals onto biochar. Analogous conclusions were

1162

reached by other scholars (Swiatkowski et al., 2004; Uchimiya et al., 2010; Rivera-Utrilla and

1163

Sánchez-Polo, 2011).

TE D

EP

The contribution of Cπ-cation interactions to the mechanism of heavy-metal adsorption

AC C

1164

M AN U

1154

1165

onto carbonaceous material cannot be denied. However, unlike biochar and activated carbon,

1166

biosorbents (commonly derived from lignocellulose materials) do not possess a graphitic

1167

structure. In addition, although hydorchar prepared through hydrothermal carbonization (i.e., 150–

1168

250 °C controlled temperature for 2–72 h at a specific pressure) exhibits an aromatic carbon

1169

network, its properties are similar to those of biosorbents (Tran et al., 2017c). Therefore, the

58

ACCEPTED MANUSCRIPT

1170

presence of Cπ-cation interactions in heavy-metal adsorption onto biosorbents and hydrochar is

1171

still controversial and needs to be further investigated. Bui and Choi (2010) considered the role of cation-π bonding between the protonated

1173

amino group of oxytetracycline (OTC) and the graphitic π electrons of multiwalled carbon

1174

nanotubes (MWCNTs). Essentially, the cation-π bonding is dominated by the electrostatic forces

1175

between the cation and the permanent quadrupole of the π-electron-rich aromatic structure and

1176

cation-induced polarization (Ji et al., 2009). Like π-π interactions, cation-π bonding would be

1177

suppressed with increasing pH. This is because at a higher pH, deprotonation of the charged amino

1178

group and protonation of the enol groups will be enhanced and so the electron-acceptor ability of

1179

the oxytetracycline molecule is weaker (Ji et al., 2009). Bui and Choi (2010) opined that

1180

Oleszczuk et al. (2010) ignored the important roles of Ca2+ in the solution; for example,

1181

Ca2+ possibly forms a complex with OTC in the solution or Ca2+ can act as a cation bridge between

1182

the negatively charged OTC and the negative charges of the MWCNTs. Therefore, Bui and Choi

1183

(2010) suggested that Ca2+ may simultaneously bind with the negatively charged OTC and interact

1184

either with the negative charges or the graphitic π electrons of the MWCNTs via cation-π bonding.

1185

As pH is increased, both MWCNTs and OTC become more negatively charged and, consequently,

1186

probably interact more strongly with Ca2+, leading to lower desorption of OTC. This hypothesis

1187

could be further verified by studying the desorption of OTC in the presence of NaCl, instead of

1188

CaCl2.

1189

8.6.2. π-π interactions

AC C

EP

TE D

M AN U

SC

RI PT

1172

1190

π-π interactions (also known as π-π electron donor-acceptor interactions) between

1191

the π-electrons in a carbonaceous adsorbent and the π-electron in the aromatic ring of an adsorbate

1192

were initially proposed by Coughlin and Ezra (1968). In short, the addition of

59

ACCEPTED MANUSCRIPT

electron-withdrawing groups (i.e., oxygen-containing and nitrogen-containing functional groups)

1194

at the edges of the individual graphene layers within a carbonaceous solid causes a considerable

1195

drop in the π-electron density. Positive holes are consequently created in the conduction π-band of

1196

the individual graphene layers, and the interactions between the π-electrons of a carbonaceous

1197

adsorbent and the π-electrons of the adsorbate aromatic rings become weaker.

RI PT

1193

Recently, Tran et al. (2017b), and Tran et al. (2017d) investigated the mechanisms of

1199

methylene green 5 (MG5) adsorption by commercial activated charcoal, synthesized activated

1200

carbons, and prepared biochars. They proposed that the primary mechanisms in MG5 adsorption

1201

were π-π interactions and pore filling, while hydrogen bonding and n-π interactions were minor

1202

contributors (Figure 11). To identify the existence of π-π interactions, they used two pieces of

1203

experimental evidence. First, FTIR analysis showed that a peak corresponding to the skeletal

1204

vibration of aromatic C=C bonds decreased in intensity and upshifted after MG5

1205

adsorption. Secondly, oxygenation of the surface of the carbonaceous solids (i.e., biochar and

1206

activated carbon) through a hydrothermal process with acrylic acid resulted in a decrease in MG5

1207

adsorption and indicated the importance of π-π interactions to the adsorption process.

TE D

M AN U

SC

1198

As discussed in Section 8.6.1, hydrochar also possesses an aromatic carbon network;

1209

therefore, π-π interactions might be expected between the π-electrons of the aromatic C=C bonds

1210

in hydrochar and the π-electrons of an adsorbate. However, a recent study (Tran et al., 2017e)

1211

demonstrated that oxygenation of a hydrochar surface through a hydrothermal process with acrylic

1212

acid contributed to increasing MG5 adsorption and indicated the negligible role of π-π interaction

1213

in the adsorption process. FTIR analysis demonstrated that the aromatic C=C peak did not

1214

significantly decrease in intensity or shift toward higher/lower wavenumbers after adsorption,

1215

which further confirms the insignificant contribution of π-π interactions. Electrostatic attraction

AC C

EP

1208

60

ACCEPTED MANUSCRIPT

played a major role in the adsorption mechanism, while hydrogen bonding and n-π interactions

1217

were minor contributors. The primary adsorption mechanisms of MG5 onto hydrochar were

1218

similar to those on biosorbents, but dissimilar to biochar and activated carbon (i.e., π-π interaction

1219

and pore filling). An identical result was obtained in the investigation of phenol, MG5, and acid

1220

read 1 adsorption onto a commercial glucose-derived spherical hydrochar functionalized with

1221

triethylenetetramine (Tran, et al., 2017a). Therefore, it can be concluded that the addition of

1222

oxygen- and nitrogen-containing functional groups (electron-withdrawing groups) on the surface

1223

of hydrochar does not cause a considerable drop in the π-electron density, but this process provides

1224

abundant adsorption sites on the hydrochar surface.

M AN U

SC

RI PT

1216

Jiang et al. (2009) and Chen et al. (2009) argued over the primary adsorption mechanisms

1226

of hydroxyl- and amino-substituted aromatics on carbon nanotubes and graphite (CNTs/graphite).

1227

Chen et al. (2009) proposed that π-π interactions played a major role in the adsorption mechanism.

1228

However, Jiang et al. (2009) remarked that the experiments of Chen et al. (2009) were

1229

insufficient to support the presence of π-π electron-donor−acceptor (EDA) interactions between

1230

hydroxyl-substituted aromatics and CNTs/graphite. According to Jiang et al. (2009), phenolic

1231

compounds were not appropriate candidates to assess the presence of π-π EDA interactions

1232

between π-electron-donating aromatics and CNTs because of the significant effects of oxygen.

1233

Therefore, the oxygen effect should be considered in more detail in order to gain in-depth insights

1234

into the adsorption mechanisms of environmentally relevant phenolic compounds onto

1235

carbonaceous materials.

AC C

EP

TE D

1225

61

ACCEPTED MANUSCRIPT

Figure 11

1236 1237

8.7. Other miscellaneous errors Azizian (2007) commented that the theoretical development of empirical PFO and PSO

1239

kinetic adsorption models using statistical rate theory (SRT) was first reported by Azizian (2004),

1240

not by Rudzinski and Plazinski (2007). Furthermore, the idea that the PFO model can only be

1241

derived theoretically at a nearly constant bulk concentration was proposed by Azizian (2004). In

1242

addition, the conclusion that the PFO and PSO models are simplified forms of a more general

1243

equation was also proposed by Azizian (2004), not Rudzinski and Plazinski (2007). Similarly,

1244

mistakes related to either no citation of the original paper or incorrect citation were pointed out by

1245

Azizian (2008b).

M AN U

SC

RI PT

1238

A series of mistakes related to incorrectly digitized data, typos, inaccurate calculations,

1247

incorrect use of conditional equilibrium constants, incorrect units, incorrect references, incorrect

1248

application of modeling procedures, and other miscellaneous errors have been highlighted by

1249

Gustafsson and Lumsdon (2014).

TE D

1246

The distinction between heat of adsorption, adsorption energy, and activation energy in

1251

adsorption was reported elsewhere (Inglezakis and Zorpas, 2012). The authors concluded that it is

1252

important to distinguish the differential enthalpy at zero coverage (∆Hdo) from the differential

1253

enthalpy or isosteric heat (∆Hd) and the standard heat of complete coverage or integral enthalpy

1254

(∆H°), which expresses the total heat generation for the adsorption and is related to the

1255

thermodynamic equilibrium constant.

AC C

1256

EP

1250

Li et al. (2010) pointed out that (1) although the cited documents were solid and sound in

1257

their scientific aspects, the reviewers seemed to fail to take full advantage of the intrinsically

1258

contained information; (2) there were some critical errors in the main text; and (3) the writing style

62

ACCEPTED MANUSCRIPT

1259

of the paper was misleading and, to some extent, lacked cautiousness. They focused on 11

1260

comments and concluded that the original paper’s conclusions were not optimistic. However, most

1261

of these comments were rejected by the original authors (Mohan and Pittman, 2011). Regarding the correct citation for determining pHPZC using the batch equilibration method,

1263

Milonjić (2009a) commented that this method was originally proposed by Milonjić et al. (1975).

1264

Furthermore, Milonjić (2009a) noted that using correct and updated citations were very important

1265

for researchers to find relevant information, pioneer ideas, and make progress in a particular

1266

subject. Additionally, from a scientific point of view, it is always necessary to give credit to the

1267

authors who first proposed a method or theoretical model. Recently, Tran et al. (2017b) applied the

1268

“drift method” to determine the pHPZC of commercial activated charcoal (CAC). The effects of

1269

various operation conditions (i.e., different degassing times with N2, background electrolytes,

1270

concentrations of an electrolyte, solid/liquid ratios, and contact times) on the pHPZC were

1271

investigated. The results demonstrated that the pHPZC (9.81 ± 0.07) of CAC was insignificantly

1272

dependent on the operation conditions.

TE D

M AN U

SC

RI PT

1262

Geckeis and Rabung (2004) pointed out that photon correlation spectroscopy (PCS) is a

1274

method for colloid size determination, not for zeta-potential analysis. Methods that have been

1275

correctly applied for zeta-potential analysis include electrophoretic light scattering and laser

1276

Doppler velocimetry (LDV), and the units for the measured zeta potentials are mV, not eV.

AC C

1277

EP

1273

Lee (2017) considered the potential effects of flow velocity on the concentration and

1278

distribution of triphenyltin chloride (TPT) on nano zinc oxide (nZnO). The author remarked that

1279

TPT was possibly re-suspended by the hydraulic flow that occurs during hydrant flushing in a

1280

reactor because TPT can be loosely deposited on the nZnO surface. As reported in the literature,

1281

the flow velocity in a wastewater treatment system is one of the factors affecting the removal

63

ACCEPTED MANUSCRIPT

efficiency for a tested pollutant. The distribution of TPT on nZnO might vary during the

1283

adsorption process according to the time of hydrant flushing; therefore, the process of TPT

1284

adsorption onto nZnO might be affected by the operating conditions of a reactor. However, the

1285

authors of the commented paper replied that the application of nZnO for the removal of TPT from

1286

dockyard wastewater was examined using a batch adsorption technique. As a result, the effects of

1287

flow velocity on the adsorption process conducted in the batch experiments were negligible

1288

compared to that in pilot studies and/or column experiments.

SC

RI PT

1282

Kamagaté et al. (2016) and Komárek et al. (2016) considered the influence of Si species

1290

originating from the partial dissolution of quartz on the competitive adsorption of Cd(II), Cr(VI),

1291

and Pb(II) onto nanomaghemite- and nanomaghemite-coated quartz. They recommended that to

1292

accurately assess the interfacial reaction mechanisms occurring at the Fe-oxide/water interface in

1293

the presence of Si-complex mineral assemblages or Fe-coated sand systems, more attention should

1294

be paid to the possible release of silicate from Si-bearing minerals and its subsequent adsorption on

1295

reactive phases.

TE D

M AN U

1289

A critical paper entitled “comments on “zirconium-carbon hybrid sorbent for removal of

1297

fluoride from water: oxalic acid mediated Zr(IV) assembly and adsorption mechanism”” was

1298

published by Zhao (2015). The authors of the original paper concluded that the fluoride adsorption

1299

occurs on the surface sites of ZrOx-AC (commercial activated carbon modified with Zr together

1300

with oxalic acid) with –OH displacement and/or interacts with the positive charge of zirconium to

1301

form –COOH groups in the oxalate acid; however, Zhao (2015) suggested that if the fluoride

1302

concentration is high enough, it would completely desorb the oxalate acid from the ZrOx-AC

1303

surface. According to the results of previous experiments, Zhao (2015) suggested that it is

AC C

EP

1296

64

ACCEPTED MANUSCRIPT

1304

necessary to further consider or add the exchangeable role of chloride and oxalic acid in the

1305

proposed fluoride-adsorption mechanism. Fan et al. (2007) additionally discussed the mechanisms of fluoride adsorption by calcined

1307

Mg-Al-CO3 layered double hydroxides (CLDH). Based on the experimental data published in the

1308

original paper and analysis of other literature results, they concluded that XRD measurements

1309

indicated that mixed Mg-Al oxides, as well as partially restored LDH, was present in

1310

fluoride-loaded CLDH after freeze-drying treatment. Furthermore, fluoride can be effectively

1311

adsorbed onto MgO, demonstrating that fluoride adsorption over the mixed oxide along with a

1312

memory effect accounts for the effective removal of fluoride by CLDH.

M AN U

SC

RI PT

1306

An analytical article on the magnitude of estimated maximum surface-area-normalized

1314

adsorption capacities (Qmax*) for the adsorption of polycyclic aromatic hydrocarbons (PAHs) and

1315

polychlorobiphenyls (PCBs) on soot and soot-like carbon materials was published by Werner and

1316

Karapanagioti (2005). The original authors (van Noort et al., 2005) used the Langmuir isotherm to

1317

extrapolate Qmax* at the solubility limit from solid-water distribution coefficients (Ks) measured in

1318

the pg/L to ng/L range along with estimated Langmuir affinities for adsorption (b) on a

1319

carbonaceous surface. After comparison with estimated surface-area-normalized monolayer

1320

adsorption capacities (Qmono*) and with empirical sorption data from the literature, Werner and

1321

Karapanagioti (2005) concluded that the Langmuir isotherm cannot be used to extrapolate

1322

maximum sorption capacities (Qmax*) near the aqueous-sorbate saturation limit from distribution

1323

coefficients measured at extremely low aqueous-sorbate concentrations. Adsorption sites in

1324

carbon materials are not uniform and other processes, such as multilayer adsorption, condensation

1325

in capillary pores, and absorption into the polymeric matrix, may be relevant near the solubility

1326

limit.

AC C

EP

TE D

1313

65

ACCEPTED MANUSCRIPT

Müller et al. (2017) analyzed strong conflicts between the XPS and XRD data and their

1328

interpretation by the original author in a study on the synthesis of Yb doped CuFe2O4 nanoferrite

1329

adsorbents (i.e., CuYb0.5Fe1.5O4) and the application of such adsorbents in the removal of different

1330

pollutants (e.g. methyl orange, safranin, Cr3 +, and Pb2 +). After a polemical discussion, Müller et

1331

al. (2017) concluded that because any experimental study stands and falls on the reliability of the

1332

investigated specimen the lack of any significance therefore holds for this study as a whole.

RI PT

1327

Interestingly, a series of recommendations and suggestions were written by the editor of

1334

Carbon for improving manuscript quality and avoiding mistakes (Thrower, 2007; Thrower,

1335

2008b; Thrower, 2008a; Thrower, 2010; Thrower, 2011). The author focused on five main topics:

1336

titles and abstracts, introduction and references, experimental, results and discussion, and

1337

language.

1338

9. Nonlinear-optimization technique

M AN U

SC

1333

To calculate the parameters of kinetic and isotherm models accurately in both batch and

1340

column experiments, application of the nonlinear method instead of the linear method has been

1341

recommended by many researchers (Porter, 1985; Ho, 2004b; Ho et al., 2005; Bolster and

1342

Hornberger, 2007; Han et al., 2007; Kumar, 2007; El-Khaiary et al., 2010; Chowdhury and Das

1343

Saha, 2011; Lima et al., 2015; Tran et al., 2015; Tran et al., 2017d).

EP

Furthermore, to identify the best-fit model, calculation of the chi-squared (χ2, Eq. 96) value

AC C

1344

TE D

1339

1345

is recommended in addition to calculating the coefficient of determination (R2, Eq. 97) for the

1346

nonlinear method. In the chi-squared test, the squares of the differences between the experimental

1347

data and data calculated using the models are divided by the corresponding data obtained and then

1348

summed. If the data obtained using a model are similar to the experimental data, χ2 is close to zero.

1349

High χ2 values indicate high bias between the experiment and model. Therefore, analyzing the data

66

ACCEPTED MANUSCRIPT

set using the chi-squared test is necessary to confirm the best-fit isotherm for a given sorption

1351

system (Ho, 2004b; Ho et al., 2005; Chowdhury and Das Saha, 2011; Lima et al., 2015; Tran et al.,

1352

2015; Tran et al., 2017b; Tran et al., 2017d; Tran et al., 2017e).

χ =∑ 2

R

2

(qe,exp − qe,cal ) 2 qe,cal

∑(q =1 − ∑ (q

e, exp

− qe, mean)2

=

∑ (q ∑ (q

e, cal

-q

e, cal 2 e, mean

- qe,mean)2

) + ∑(qe, cal - qe,exp )2

(96)

(97)

SC

e,exp

− qe,cal )2

RI PT

1350

where qe,exp (mg/g) is the amount of adsorbate uptake at equilibrium obtained from Eq. 1, qe,cal

1354

(mg/g) is the amount of adsorbate uptake achieved from the model using the Solver add-in, and

1355

qe,mean (mg/g) is the mean of the qe,exp values.

M AN U

1353

Lima et al. (2015) confirmed that the use of nonlinear equilibrium and kinetic adsorption

1357

models gives values that are reliable and statistically relevant for modeling the isotherm and

1358

kinetics of adsorption. Moreover, they did not suggest using linearized equilibrium and kinetic

1359

adsorption models, even in the cases of the linearized Langmuir Type 1 and linearized PSO Type 1

1360

equations that provide R2 values close to unity. The linearization of equilibrium and kinetic

1361

adsorption models could make the parameters determined from these models meaningless.

EP

TE D

1356

Owing to the wide application of the nonlinear method in the study of adsorption

1363

processes, we provide here a short introduction to the application of this method using the “Solver

1364

add-in” in Microsoft Office Excel (2013 version). First, it is necessary to load the Solver add-in in

1365

Microsoft Excel: File → Excel options → Add-ins (select Excel Add-ins in Manage box) → select

1366

the Solver Add-in check box. Second, all Ce and qe values obtained experimentally are considered

1367

as input data, and the qe values from the Langmuir equation are calculated based on two variables

1368

(Qomax and KL ≠ 0). Using the Solver Add-in (Figure 12), the Qomax, KL, and R2 values can be

AC C

1362

67

ACCEPTED MANUSCRIPT

1369

obtained. Notably, other programs (i.e., Origin) also provide similar results and can accurately

1370

calculate the parameters of kinetic and isotherm models.

1371

Figure 12 10. Conclusions

RI PT

1372

In scientific publications, the correct use of technical terms and accurate calculations are

1374

essential. According to the viewpoints discussed in this review, we put forward the following key

1375

conclusions and perspectives:

1378 1379 1380 1381 1382

Adsorption performance should be expressed as qe (mg/g), while the use of

the unit of %removal for qe should be avoided or used cautiously. ii.

M AN U

1377

i.

The use of accurate citations and correct mathematical expressions in

original works is also recommended. iii.

Measurements of adsorption kinetics should be started at an initial time of

less than 2 min and finish when the adsorption process reaches true equilibrium. iv.

TE D

1376

SC

1373

A complete adsorption isotherm (plotting qe versus Ce) should be presented

1383

when investigating liquid-phase adsorption to identify the regions in which the

1384

experimental data relating to adsorption equilibrium are actually located. v.

Differences between the number of data points in an experiment and those

EP

1385

used for model fitting (i.e., the pseudo-second-order or Langmuir models) should be

1387

avoided.

1388

AC C

1386

vi.

For adsorption thermodynamics, the equilibrium constant (KC) must be

1389

dimensionless. The optimal method for calculating thermodynamic parameters should

1390

be appropriately selected according to the data from the adsorption equilibrium

1391

experiment. Notably, it is impossible to calculate thermodynamic parameters from an

1392

experiment using one initial adsorbate concentration at different temperatures. 68

ACCEPTED MANUSCRIPT

1393

vii.

To accurately estimate the parameters of adsorption kinetic and isotherm

models, the nonlinear optimization technique should be applied to decrease the bias

1395

between the qe values calculated from the experimental data and those estimated from

1396

the models. The chi-squared (χ2) test should be used alongside the nonlinear

1397

determination coefficient (R2) to obtain the best-fit models for adsorption kinetics and

1398

isotherms.

1399

viii.

RI PT

1394

An in-depth understanding of the models applied to adsorption kinetics and

isotherms is necessary to avoid misapplication of these models, as well as inaccurate

1401

discussions, calculations, and conclusions. ix.

More attention should be paid to various common misunderstandings about

M AN U

1402

SC

1400

1403

fundamental adsorption and chemistry concepts. Understanding the unique properties

1404

of adsorbents and target adsorbates can help avoid mistakes made when explaining

1405

adsorption mechanisms.

Although the comments and recommendations of some authors are not always correct, their

1407

observations and suggestions should be acknowledged because of their great contributions in

1408

transferring scientific knowledge.

1409

Acknowledgements

1410

This current work was financially supported by Chung Yuan Christian University (CYCU) in

1411

Taiwan. The first author would like to thank CYCU for the Distinguished International Graduate

1412

Students (DIGS) scholarship to pursue his doctoral studies. The authors gratefully acknowledge

1413

the editor and anonymous reviewer for their invaluable insight and helpful suggestion.

1414

Reference

1415

Alagumuthu, G., V. Veeraputhiran and M. Rajan (2010). "Comments on “Fluoride removal from

1416

water using activated and MnO2-coated Tamarind Fruit (Tamarindus indica) shell:

1417

Batch and column studies”." Journal of Hazardous Materials 183(1–3): 956-957. 69

AC C

EP

TE D

1406

ACCEPTED MANUSCRIPT

1418

Anastopoulos, I. and G. Z. Kyzas (2016). "Are the thermodynamic parameters correctly

1419

estimated in liquid-phase adsorption phenomena?" Journal of Molecular Liquids 218:

1420

174-185. Aoyama, M. (2003). "Comment on “Biosorption of chromium(VI) from aqueous solution by

1422

cone biomass of Pinus sylvestris”." Bioresource Technology 89(3): 317-318.

1423

Arica, M. Y. (2003). "In reaction to the comment by Dr. Y.-S. Ho on our publication “Affinity dye-ligand

composite

membrane

for

1425

adsorption lysozyme and kinetic properties, Biochemical Engineering Journal 13

1426

(2003) 35–45”." Biochemical Engineering Journal 15(1): 79-80.

SC

1424

Azizian, S. (2004). "Kinetic models of sorption: a theoretical analysis." Journal of Colloid and Interface Science 276(1): 47-52.

1428

M AN U

1427

poly(hydroxyethylmethacrylate/chitosan

RI PT

1421

1429

Azizian, S. (2006). "A novel and simple method for finding the heterogeneity of adsorbents on

1430

the basis of adsorption kinetic data." Journal of Colloid and Interface Science 302(1):

1431

76-81.

Azizian, S. (2007). "Comment on “Kinetics of Solute Adsorption at Solid/Solution Interfaces:  A

1433

Theoretical Development of the Empirical Pseudo-First and Pseudo-Second Order

1434

Kinetic Rate Equations, Based on Applying the Statistical Rate Theory of Interfacial

1435

Transport”." The Journal of Physical Chemistry B 111(1): 318-318.

TE D

1432

Azizian, S. (2008a). "Comments on “Adsorption Equilibrium and Kinetics of Reactive Black 5

1437

and Reactive Red 239 in Aqueous Solution onto Surfactant-Modified Zeolite”

1438

(Karadag, D.; Turan, M.; Akgul, E.; Tok, S.; Faki, A. J. Chem. Eng. Data 2007, 52,

1439

1615−1620)." Journal of Chemical & Engineering Data 53(1): 322-323.

1441

AC C

1440

EP

1436

Azizian, S. (2008b). "Comments on “Biosorption isotherms, kinetics and thermodynamics” review." Separation and Purification Technology 63(2): 249-250.

1442

Ben Ali, S., I. Jaouali, S. Souissi-Najar and A. Ouederni (2016). "Characterization and adsorption

1443

capacity of raw pomegranate peel biosorbent for copper removal." Journal of Cleaner

1444

Production.

1445

Biggar,

J.

W.

and

M.

W.

Cheung 70

(1973).

"Adsorption

of

Picloram

ACCEPTED MANUSCRIPT

1446

(4-Amino-3,5,6-Trichloropicolinic Acid) on Panoche, Ephrata, and Palouse Soils: A

1447

Thermodynamic Approach to the Adsorption Mechanism1." Soil Science Society of

1448

America Journal 37(6).

1450 1451 1452

Blanchard, G., M. Maunaye and G. Martin (1984). "Removal of heavy metals from waters by

RI PT

1449

means of natural zeolites." Water Research 18(12): 1501-1507.

Bolster, C. H. and G. M. Hornberger (2007). "On the Use of Linearized Langmuir equation." Soil Science Society of America Journal 71(6): 1796-1806.

Borisover, M. D. and E. R. Graber (1997). "Comment on “Competitive Sorption between

1454

Atrazine and Other Organic Compounds in Soils and Model Sorbents”."

1455

Environmental Science & Technology 31(5): 1577-1577.

M AN U

SC

1453

1456

Bui, T. X. and H. Choi (2010). "Comment on “Adsorption and Desorption of Oxytetracycline

1457

and Carbamazepine by Multiwalled Carbon Nanotubes”." Environmental Science &

1458

Technology 44(12): 4828-4828.

Canzano, S., P. Iovino, S. Salvestrini and S. Capasso (2012). "Comment on “Removal of anionic

1460

dye Congo red from aqueous solution by raw pine and acid-treated pine cone powder

1461

as adsorbent: Equilibrium, thermodynamic, kinetics, mechanism and process design”."

1462

Water Research 46(13): 4314-4315.

1463

TE D

1459

Cavas, L. (2008). Comment on Equilibrium sorption isotherm studies of Cd(II), Pb(II) and Zn (II) ions detoxification from waste water using unmodified and EDTA-modified maize

1465

husk.

1467 1468

Chen, W., L. Duan, L. Wang and D. Zhu (2009). "Response to Comment on “Adsorption of

AC C

1466

EP

1464

Hydroxyl- and Amino-Substituted Aromatics to Carbon Nanotubes”." Environmental

Science & Technology 43(9): 3400-3401.

1469

Chien, S. H. and W. R. Clayton (1980). "Application of Elovich Equation to the Kinetics of

1470

Phosphate Release and Sorption in Soils." Soil Science Society of America Journal

1471

44(2): 265-268.

1472

Choong, T. S. Y. and T. G. Chuah (2005). "Comment on “Separation of vitamin E from palm

1473

fatty acid distillate using silica: I. Equilibrium of batch adsorption by B.S. Chu et al. 71

ACCEPTED MANUSCRIPT

1474

[Journal of Food Engineering 62 (2004) 97–103]”." Journal of Food Engineering 67(3):

1475

379. Chowdhury, S. and P. Das Saha (2011). "Comparative Analysis of Linear and Nonlinear Methods

1477

of Estimating the Pseudo-Second-Order Kinetic Parameters for Sorption of Malachite

1478

Green onto Pretreated Rice Husk." Bioremediation Journal 15(4): 181-188.

RI PT

1476

Chu, B. S., B. S. Baharin, Y. B. C. Man and S. Y. Quek (2005). "Reply from the authors to the

1480

comments by Choong and Chuah in a recent Letter to the Editor on “Separation of

1481

vitamin E from palm fatty acid distillate using silica: I. Equilibrium of batch

1482

adsorption by B.S. Chu et al. [Journal of Food Engineering 62 (2004) 97–103]”."

1483

Journal of Food Engineering 69(1): 131-132.

SC

1479

Coughlin, R. W. and F. S. Ezra (1968). "Role of surface acidity in the adsorption of organic

1485

pollutants on the surface of carbon." Environmental Science & Technology 2(4):

1486

291-297.

M AN U

1484

Crittenden, J. C., R. R. Trussell, D. W. Hand, K. J. Howe and G. Tchobanoglous (2012).

1488

Adsorption. MWH's Water Treatment: Principles and Design, Third Edition, John

1489

Wiley & Sons, Inc.: 1117-1262.

1490 1491

TE D

1487

D4607-14, A. (2014). "Standard Test Method for Determination of Iodine Number of Activated Carbon." American Society for Testing and Materials, Philadelphia, U.S.A. Dambies, L., C. Guimon, S. Yiacoumi and E. Guibal (2001). "Characterization of metal ion

1493

interactions with chitosan by X-ray photoelectron spectroscopy." Colloids and

1494

Surfaces A: Physicochemical and Engineering Aspects 177(2–3): 203-214.

AC C

EP

1492

1495

Dawood, S. and T. K. Sen (2012). "Author's Responses to the comment by Canzano et al and

1496

also corrigendum to “Removal of anionic dye Congo red from aqueous solution by

1497 1498 1499 1500 1501

raw pine and acid-treated pine cone powder as adsorbent: Equilibrium, thermodynamic,

kinetics, mechanism and process design” published in Water Research, Vol. 46, pp.

1933–1946, 2012." Water Research 46(13): 4316-4317. Doke, K. M. and E. M. Khan (2013). "Adsorption thermodynamics to clean up wastewater; critical review." Reviews in Environmental Science and Bio/Technology 12(1): 25-44.

72

ACCEPTED MANUSCRIPT

1503 1504 1505

Dubinin, M. and L. Radushkevich (1947). "Equation of the characteristic curve of activated charcoal." Chem. Zentr 1(1): 875. Egirani, D. (2004). "Comments on: Removal of copper ions from aqueous solution by tree ferns." Water Research 38(20): 4535.

RI PT

1502

1506

El-Khaiary, M. I., G. F. Malash and Y.-S. Ho (2010). "On the use of linearized

1507

pseudo-second-order kinetic equations for modeling adsorption systems." Desalination

1508

257(1–3): 93-101.

Fan, J., Z. Xu and S. Zheng (2007). "Comment on “Factors influencing the removal of fluoride

1510

from aqueous solution by calcined Mg–Al–CO3 layered double hydroxides”." Journal

1511

of Hazardous Materials 139(1): 175-177.

M AN U

SC

1509

1512

Farooq, U., J. A. Kozinski, M. A. Khan and M. Athar (2010). "Biosorption of heavy metal ions

1513

using wheat based biosorbents – A review of the recent literature." Bioresource

1514

Technology 101(14): 5043-5053.

Febrianto, J., A. N. Kosasih, J. Sunarso, Y.-H. Ju, N. Indraswati and S. Ismadji (2009).

1516

"Equilibrium and kinetic studies in adsorption of heavy metals using biosorbent: A

1517

summary of recent studies." Journal of Hazardous Materials 162(2–3): 616-645.

1518 1519

TE D

1515

Foo, K. Y. and B. H. Hameed (2010). "Insights into the modeling of adsorption isotherm systems." Chemical Engineering Journal 156(1): 2-10. Freundlich, H. (1906). "Uber die adsorption in lo sungen." Z Physics Chemistry 57: 385-471.

1521

Fu, L., X. Huang and J. Wang (2007). "Comment on “Removal of Supranol Yellow 4GL by

1522

adsorption onto Cr-intercalated montmorillonite”." Journal of Hazardous Materials

AC C

1523

EP

1520

148(3): 771-772.

1524

Fu, L., J. Wang, H. Lu, Y. Su and A. Ren (2008). "Comment on “The removal of phenolic

1525

compounds from aqueous solutions by organophilic bentonite”." Journal of Hazardous

1526

Materials 151(2–3): 851-854.

1527

Gardea-Torresdey, J. L., K. J. Tiemann, V. Armendariz, L. Bess-Oberto, R. R. Chianelli, J. Rios,

1528

J. G. Parsons and G. Gamez (2000). "Characterization of Cr(VI) binding and reduction

73

ACCEPTED MANUSCRIPT

1529

to Cr(III) by the agricultural byproducts of Avena monida (Oat) biomass." Journal of

1530

Hazardous Materials 80(1–3): 175-188. Geckeis, H. and T. Rabung (2004). "Comment on the paper: X. Wang, Th. Rabung, H. Geckeis,

1532

Effect of pH and humic acid on the adsorption of cesium onto γ-Al2O3." Journal of

1533

Radioanalytical and Nuclear Chemistry 260(2): 437-438.

RI PT

1531

Ghosal, P. S. and A. K. Gupta (2015). "An insight into thermodynamics of adsorptive removal of

1535

fluoride by calcined Ca-Al-(NO3) layered double hydroxide." RSC Advances 5(128):

1536

105889-105900.

1538

Ghosal, P. S. and A. K. Gupta (2017). "Determination of thermodynamic parameters from Langmuir isotherm constant-revisited." Journal of Molecular Liquids 225: 137-146.

M AN U

1537

SC

1534

1539

Goldberg, S. (2009). "Comments on “Adsorption Sequence of Toxic Inorganic Anions on a Soil”

1540

by K. Saeki Bull Environ Contam Toxicol (2008) 81:508–512." Bulletin of

1541

Environmental Contamination and Toxicology 83(1): 1-2.

Goldberg, S. (2016). "A comment on “Selenite adsorption and desorption in main Chinese soils

1543

with their characteristics and physicochemical properties.” J Soils Sediments

1544

15(5):1150–1158 (2015), doi:10.1007/s11368-015-1085-7." Journal of Soils and

1545

Sediments 16(1): 324-324.

TE D

1542

Goswami, R., J. Shim, S. Deka, D. Kumari, R. Kataki and M. Kumar (2016). "Characterization

1547

of cadmium removal from aqueous solution by biochar produced from Ipomoea

1548

fistulosa at different pyrolytic temperatures." Ecological Engineering 97: 444-451.

1550 1551

Guo, H., S. Zhang, Z. Kou, S. Zhai, W. Ma and Y. Yang (2015). "Removal of cadmium(II) from

AC C

1549

EP

1546

aqueous solutions by chemically modified maize straw." Carbohydrate Polymers

115(0): 177-185.

1552

Gustafsson, J. P. and D. G. Lumsdon (2014). "Comment on “Citrate adsorption can decrease

1553

soluble phosphate concentration in soils: Results of theoretical modelling” by Marek

1554

Duputel, Nicolas Devau, Michel Brossard, Benoît Jaillard, Davey L. Jones, Philippe

1555

Hinsinger and Frédéric Gérard (2013)." Applied Geochemistry 46: 85-89.

1556

Hai, T. N. (2017). "Comments on “Effect of Temperature on the Adsorption of Methylene Blue 74

ACCEPTED MANUSCRIPT

1557

Dye

onto

Sulfuric

Acid–Treated

1558

Communications 204(1): 134-139.

Orange

Peel”."

Chemical

Engineering

Halder, G., A. A. Khan and S. Dhawane (2016). "Fluoride Sorption Onto a Steam-Activated

1560

Biochar Derived From Cocos nucifera Shell." CLEAN – Soil, Air, Water 44(2):

1561

124-133.

RI PT

1559

Hall, K. R., L. C. Eagleton, A. Acrivos and T. Vermeulen (1966). "Pore- and Solid-Diffusion

1563

Kinetics in Fixed-Bed Adsorption under Constant-Pattern Conditions." Industrial &

1564

Engineering Chemistry Fundamentals 5(2): 212-223.

SC

1562

Hamdaoui, O. and E. Naffrechoux (2007). "Modeling of adsorption isotherms of phenol and

1566

chlorophenols onto granular activated carbon: Part II. Models with more than two

1567

parameters." Journal of Hazardous Materials 147(1–2): 401-411.

M AN U

1565

1568

Han, R., Y. Wang, W. Zou, Y. Wang and J. Shi (2007). "Comparison of linear and nonlinear

1569

analysis in estimating the Thomas model parameters for methylene blue adsorption

1570

onto natural zeolite in fixed-bed column." Journal of Hazardous Materials 145(1–2):

1571

331-335.

Harsha, N., K. V. S. Krishna, N. K. Renuka and S. Shukla (2015). "Facile synthesis of

1573

[gamma]-Fe2O3 nanoparticles integrated H2Ti3O7 nanotubes structure as a

1574

magnetically recyclable dye-removal catalyst." RSC Advances 5(38): 30354-30362.

1577 1578

EP

1576

Ho, Y.-S. (1995). "Adsorption of heavy metals from waste streams by peat." Ph.D. Thesis, The University of Birmingham, Birmingham, U.K. Ho, Y.-S. (2004a). "Erratum to “Removal of copper ions from aqueous solution by tree fern”:

AC C

1575

TE D

1572

[Water Res. 37 (2003) 2323–2330]." Water Research 38(20): 4536-4537.

1579

Ho, Y.-S. (2004b). "Selection of optimum sorption isotherm." Carbon 42(10): 2115-2116.

1580

Ho, Y.-S. (2005). "Comment on “Two-stage batch sorber design using second-order kinetic

1581

model for the sorption of metal complex dyes onto pine sawdust” by Özacar, M. and

1582 1583

Şengýl, İ.A." Biochemical Engineering Journal 23(3): 291-292. Ho, Y.-S. (2013). "Comments on “Removal of zirconium(IV) from aqueous solution by Coriolus

75

ACCEPTED MANUSCRIPT

1584

versicolor: Equilibrium and thermodynamic study”." Ecological Engineering 58: 1-2. Ho, Y.-S. (2014a). "Comments on “Adsorption of 2-mercaptobenzothiazole from aqueous

1586

solution by organo-bentonite” by P. Jing, M.H. Hou, P. Zhao, X.Y. Tang, H.F. Wan."

1587

Journal of Environmental Sciences 26(12): 2571-2572.

RI PT

1585

1588

Ho, Y.-S. (2014b). "Comment on “New Calix[4]arene Appended Amberlite XAD-4 Resin with

1589

Versatile Perchlorate Removal Efficiency”." Journal of Chemical & Engineering Data

1590

59(6): 2131-2131.

1592

Ho, Y.-S. (2014c). "Comments on “Adsorption characteristics and behaviors of graphene oxide

SC

1591

for Zn(II) removal from aqueous solution”." Applied Surface Science 301: 584. Ho, Y.-S. (2014d). "Comments on “Simultaneous adsorption of aniline and Cr(VI) Ion by

1594

activated carbon/chitosan composite”." Journal of Applied Polymer Science 131(22):

1595

n/a-n/a.

1597 1598 1599 1600 1601

Ho, Y.-S. (2014e). "Comments on the paper “Removal of Cu (II) by loofah fibers as a natural and low-cost adsorbent from aqueous solutions”." Journal of Molecular Liquids 198: 322. Ho, Y.-S. (2014f). "The real pseudo-second-order rate equation." Industrial Crops and Products

52: 17.

TE D

1596

M AN U

1593

Ho, Y.-S. (2014g). "Using of “pseudo-second-order model” in adsorption." Environmental Science and Pollution Research 21(11): 7234-7235. Ho, Y.-S. (2015). "Comment on “Genetic characterization, nickel tolerance, biosorption, kinetics,

1603

and uptake mechanism of a bacterium isolated from electroplating industrial

1604

effluent”." Canadian Journal of Microbiology 61(11): 881-882.

AC C

EP

1602

1605

Ho, Y.-S. (2016a). "Comments on the paper ‘Synthesis and application of ion-imprinted resin

1606

based on modified melamine-thiourea for selective removal of Hg(II)’." Polymer

1607

International 65(7): 845-846.

1608

Ho, Y.-S. (2016b). "Comments on the paper “Removal of aqueous Hg(II) and Cr(VI) using

1609

phytic acid doped polyaniline/cellulose acetate composite membrane”." Journal of

1610

Hazardous Materials 311: 273-274.

76

ACCEPTED MANUSCRIPT

1611

Ho, Y.-S. (2016c). "Comments on the paper “Simultaneous removal of mixed contaminants by

1612

organoclays — Amoxicillin and Cu(II) from aqueous solution”." Applied Clay Science

1613

126: 297-298.

1615

Ho, Y.-S. (2016d). "Comments on using of “pseudo-first-order model” [J. Taiwan Inst. Chem.

RI PT

1614

Eng. Vol. 59]." Journal of the Taiwan Institute of Chemical Engineers 68: 136-137. Ho, Y.-S. (2016e). "Comments on using of “pseudo-first-order model” in adsorption [Int. J. Biol.

1617

Macromol., vol. 81]." International Journal of Biological Macromolecules 88:

1618

505-506.

SC

1616

Ho, Y.-S. (2017). "Comments on the paper “Application of Mn/MCM-41 as an adsorbent to

1620

remove Methyl Blue from aqueous solution”." Journal of Colloid and Interface

1621

Science.

M AN U

1619

Ho, Y.-S., W.-T. Chiu and C.-C. Wang (2005). "Regression analysis for the sorption isotherms of

1623

basic dyes on sugarcane dust." Bioresource Technology 96(11): 1285-1291.

1624

Ho, Y. S. and G. McKay (1998a). "A Comparison of Chemisorption Kinetic Models Applied to

1625

Pollutant Removal on Various Sorbents." Process Safety and Environmental Protection

1626

76(4): 332-340.

1629 1630 1631 1632 1633

Engineering Journal 70(2): 115-124.

Ho, Y. S., D. A. J. Wase and C. F. Forster (1996). "Kinetic Studies of Competitive Heavy Metal

EP

1628

Ho, Y. S. and G. McKay (1998b). "Sorption of dye from aqueous solution by peat." Chemical

Adsorption by Sphagnum Moss Peat." Environmental Technology 17(1): 71-77. Huang, R., B. Yang, Q. Liu and Y. Liu (2014). "Simultaneous adsorption of aniline and Cr(VI)

AC C

1627

TE D

1622

ion by activated carbon/chitosan composite." Journal of Applied Polymer Science

131(4): n/a-n/a.

1634

Hung, H.-W. and T.-F. Lin (2006). "Adsorption of MTBE from contaminated water by

1635

carbonaceous resins and mordenite zeolite." Journal of Hazardous Materials 135(1–3):

1636 1637

210-217. Inbaraj, B. S. (2006). "Comment on “Adsorption of Reactive Dyes from a Textile Effluent Using

77

ACCEPTED MANUSCRIPT

1638

Sawdust as the Adsorbent”." Industrial & Engineering Chemistry Research 45(21):

1639

7362-7362. Inglezakis, V. J. and A. A. Zorpas (2012). "Heat of adsorption, adsorption energy and activation

1641

energy in adsorption and ion exchange systems." Desalination and Water Treatment

1642

39(1-3): 149-157.

RI PT

1640

Ji, L., W. Chen, L. Duan and D. Zhu (2009). "Mechanisms for strong adsorption of tetracycline

1644

to carbon nanotubes: A comparative study using activated carbon and graphite as

1645

adsorbents." Environmental Science & Technology 43(7): 2322-2327.

SC

1643

Jiang, J., S.-Y. Pang and J. Ma (2009). "Comment on “Adsorption of Hydroxyl- and

1647

Amino-Substituted Aromatics to Carbon Nanotubes”." Environmental Science &

1648

Technology 43(9): 3398-3399.

M AN U

1646

1649

Jing, P., M. Hou, P. Zhao, X. Tang and H. Wan (2014). "Reply to comments on “Adsorption of

1650

2-mercaptobenzothiazole from aqueous solution by organo-bentonite” by Yuhshan

1651

Ho." Journal of Environmental Sciences 26(12): 2573-2574. Kamagaté, M., J. Lützenkirchen, F. Huber and K. Hanna (2016). "Comment on “Competitive

1653

Adsorption of Cd(II), Cr(VI), and Pb(II) onto Nanomaghemite: A Spectroscopic and

1654

Modeling Approach”." Environmental Science & Technology 50(3): 1632-1633.

TE D

1652

Karadag, D., M. Turan, E. Akgul, S. Tok and A. Faki (2007). "Adsorption Equilibrium and

1656

Kinetics of Reactive Black 5 and Reactive Red 239 in Aqueous Solution onto

1657

Surfactant-Modified Zeolite." Journal of Chemical & Engineering Data 52(5):

1658

1615-1620.

1660

AC C

1659

EP

1655

Khan, A. A. and R. P. Singh (1987). "Adsorption thermodynamics of carbofuran on Sn (IV) arsenosilicate in H+, Na+ and Ca2+ forms." Colloids and Surfaces 24(1): 33-42.

1661

Komárek, M., C. M. Koretsky, K. J. Stephen, D. S. Alessi and V. Chrastný (2016). "Response to

1662

Comment on “Competitive Adsorption of Cd(II), Cr(VI), and Pb(II) onto

1663

Nanomaghemite: A Spectroscopic and Modeling Approach”." Environmental Science

1664

& Technology 50(3): 1634-1635.

1665

Kim, Y., J. Bae, H. Park, J.-K. Suh, Y.-W. You and H. Choi (2016). "Adsorption dynamics of 78

ACCEPTED MANUSCRIPT

1666

methyl violet onto granulated mesoporous carbon: Facile synthesis and adsorption

1667

kinetics." Water Research 101: 187-194. Kosmulski, M. (2009). Surface charging and points of zero charge, CRC Press.

1669

Kratochvil, D. and B. Volesky (1998). "Advances in the biosorption of heavy metals." Trends in

1671 1672 1673 1674

Biotechnology 16(7): 291-300.

Kumar, K. V. (2006a). "Comments on “Adsorption of acid dye onto organobentonite”." Journal of hazardous materials 137(1): 638-639.

Kumar, K. V. (2006b). "Comments on “Equilibrium studies for the adsorption of Acid dye onto

SC

1670

RI PT

1668

modified hectorite”." Journal of Hazardous Materials 137(2): 1252-1253. Kumar, K. V. (2006c). "A note on the comments by Dr. Y.S. Ho on “Nitrate removal from

1676

aqueous solution by adsorption onto various materials”." Journal of Hazardous

1677

Materials 136(3): 995-996.

M AN U

1675

Kumar, K. V. (2006d). "A note on the comments by Dr. Y.S. Ho on “Remediation of soil

1679

contaminated with the heavy metal (Cd2+)”." Journal of Hazardous Materials 136(3):

1680

993-994.

1681 1682 1683

TE D

1678

Kumar, K. V. (2007). "Optimum sorption isotherm by linear and non-linear methods for malachite green onto lemon peel." Dyes and Pigments 74(3): 595-597. Kumar, K. V. and V. T. Fávere (2006). "Reply to Comments on ‘Chitosan functionalized with 2[-bis-(pyridylmethyl)

aminomethyl]4-methyl-6-formyl-phenol:

1685

kinetics of copper(II) adsorption’ by Yuh-Shan Ho: Discussion on pseudo second order

1686

kinetic expression." Polymer 47(6): 1772-1773.

Equilibrium

and

AC C

EP

1684

1687

Kumar, K. V. and B. Guha (2006). "Reply to the comments on “Study on biosorption of Cr (VI)

1688

by Mucor hiemalis” by Y.-S. Ho, Biochem. Eng. J. 26 (2005) 82–83." Biochemical

1689

Engineering Journal 30(2): 222-223.

1690

Kumar, K. V. and K. Porkodi (2006). "Relation between some two- and three-parameter isotherm

1691

models for the sorption of methylene blue onto lemon peel." Journal of Hazardous

1692

Materials 138(3): 633-635.

79

ACCEPTED MANUSCRIPT

1693

Kumar, K. V. and K. Porkodi (2007a). "Comments on “adsorption of 4-chlorophenol from

1694

aqueous solutions by xad-4 resin: Isotherm, kinetic, and thermodynamic analysis”."

1695

Journal of Hazardous Materials 143(1–2): 598-599.

1697

Kumar, K. V. and K. Porkodi (2007b). "Comments on “Biosorption of nickel from protonated

RI PT

1696

rice bran”." Journal of Hazardous Materials 147(1–2): 679.

Kumar, K. V. and K. Porkodi (2008). "Comments on “Removal of Congo red from aqueous

1699

solution by anilinepropylsilica xerogel” by Pavan FA, Dias SLP, Lima EC, Benvenutti

1700

EV. Dyes and Pigments 2008;76:64–9." Dyes and Pigments 77(2): 481-482.

1701

Kumar, K. V., K. Porkodi and F. Rocha (2007). "Comments on “Removal of lead from aqueous

1702

solution using Syzygium cumini L.: Equilibrium and kinetic studies”." Journal of

1703

Hazardous Materials 147(1–2): 677-678.

1705

M AN U

1704

SC

1698

Lagergren, S. (1898). "About the theory of so-called adsorption of soluble substances." Kungliga Svenska Vetenskapsakademiens Handlingar 24(4): 1-39.

Langmuir, I. (1918). "THE ADSORPTION OF GASES ON PLANE SURFACES OF GLASS,

1707

MICA AND PLATINUM." Journal of the American Chemical Society 40(9):

1708

1361-1403.

TE D

1706

Lee, D.-G. (2017). "Comment on: Olushola S, Olalekan S, Folahan A, Leslie F and Ximba BJ

1710

(2014) Application of nano zinc oxide (nZnO) for the removal of triphenyltin chloride

1711

(TPT) from dockyard wastewater." Water SA (Online) 43(1): 175.

1713

Li, W. and G. Pan (2007). "Comments on “Ion adsorption components in liquid/solid systems”." Journal of Environmental Sciences 19(6): 768.

AC C

1712

EP

1709

1714

Li, Z., S. Imaizumi, T. Katsumi, X. Tang and Y. Chen (2010). "Comment on JHM 142 (2007) 1–

1715

53 ‘Arsenic removal from water–wastewater using adsorbents—A critical review’ by

1716

D Mohan and CU Pittman Jr." Journal of Hazardous Materials 175(1–3): 1116-1117.

1717

Lima, É. C., M. A. Adebayo and F. M. Machado (2015). Kinetic and equilibrium models of

1718

adsorption. Carbon Nanomaterials as Adsorbents for Environmental and Biological

1719

Applications, Springer: 33-69.

80

ACCEPTED MANUSCRIPT

Lima, E. C., A. R. Cestari and M. A. Adebayo (2016). "Comments on the paper: a critical review

1721

of the applicability of Avrami fractional kinetic equation in adsorption-based water

1722

treatment studies." Desalination and Water Treatment 57(41): 19566-19571.

1723

Lin, J. X. and L. Wang (2009). "Comment on “Adsorption of Supranol Yellow 4 GL from

1724

aqueous solution by surfactant-treated aluminum/chromium-intercalated bentonite”."

1725

Journal of Hazardous Materials 172(1): 516-517.

1727

Liu, Y. (2009). "Is the Free Energy Change of Adsorption Correctly Calculated?" Journal of Chemical & Engineering Data 54(7): 1981-1985.

SC

1726

RI PT

1720

Lu, X. (2008). "Comment on “Thermodynamic and isotherm studies of the biosorption of Cu(II),

1729

Pb(II), and Zn(II) by leaves of saltbush (Atriplex canescens)”." The Journal of

1730

Chemical Thermodynamics 40(4): 739-740.

1731

Marsh, H. and F. R. Reinoso (2006). Activated carbon, Elsevier.

1732

Masson, S., M. Gineys, S. Delpeux-Ouldriane, L. Reinert, S. Guittonneau, F. Béguin and L.

1733

Duclaux (2016). "Single, binary, and mixture adsorption of nine organic contaminants

1734

onto a microporous and a microporous/mesoporous activated carbon cloth."

1735

Microporous and Mesoporous Materials 234: 24-34.

1737

TE D

1736

M AN U

1728

McKay, G., Y. S. Ho and J. C. Y. Ng (1999). "Biosorption of Copper from Waste Waters: A Review." Separation & Purification Reviews 28(1): 87-125. McLintock, I. (1967). "The Elovich equation in chemisorption kinetics." Nature 216: 1204-1205.

1739

Medellin-Castillo, N. A., E. Padilla-Ortega, M. C. Regules-Martínez, R. Leyva-Ramos, R.

1740

Ocampo-Pérez and C. Carranza-Alvarez (2017). "Single and competitive adsorption of

1742 1743 1744 1745 1746

AC C

1741

EP

1738

Cd(II) and Pb(II) ions from aqueous solutions onto industrial chili seeds (Capsicum

annuum) waste." Sustainable Environment Research.

Milonjić, S. K. (2007). "A consideration of the correct calculation of thermodynamic parameters of adsorption." Journal of the Serbian chemical society 72(12): 1363-1367.

Milonjić, S. K. (2009a). "Comments on “factors influencing the removal of divalent cations by hydroxyapatite”." Journal of Hazardous Materials 162(2–3): 1588-1589.

81

ACCEPTED MANUSCRIPT

1747

Milonjić, S. K. (2009b). "Comments on “removal of uranium (VI) from aqueous solution by

1748

adsorption of hematite”, by X. Shuibo, Z. Chun, Z. Xinghuo, Y. Jing, Z. Xiaojian, W.

1749

Jingsong." Journal of Environmental Radioactivity 100(10): 921-922. Milonjić, S. K. (2010). "Comments on the authors’ response to the comments on “Factors

1751

influencing the removal of divalent cations by hydroxyapatite”, by Smiciklas et al."

1752

Journal of Hazardous Materials 176(1–3): 1126-1127.

1754

Milonjić, S. K., A. L. Ruvarac and M. V. Šušić (1975). "The heat of immersion of natural magnetite in aqueous solutions." Thermochimica Acta 11(3): 261-266.

SC

1753

RI PT

1750

Mohan, D. and C. U. Pittman Jr (2011). "Reply to the comments on HAZMAT 142 (2007) 1–53

1756

‘Arsenic removal from water/wastewater using adsorbents – A critical review’ by D.

1757

Mohan and C.U. Pittman Jr. made by Zhenze Li et al. [HAZMAT 175 (2010) 1116–

1758

1117]." Journal of Hazardous Materials 185(2–3): 1614-1617.

M AN U

1755

1759

Müller, G., K. Janošková, T. Bakalár, J. Cakl and H. Jiránková (2012). "Removal of Zn(II) from

1760

aqueous solutions using Lewatit S1468." Desalination and Water Treatment 37(1-3):

1761

146-151.

Müller, F., T. Faidt, S. Grandthyll and H. Natter (2017). "Comment on “CuYb0.5Fe1.5O4

1763

nanoferrite adsorbent structural, morphological and functionalization characteristics

1764

for multiple pollutant removal by response surface methodology” by M.A. Rehman et

1765

al., J. Mol. Liq. 224 (2016) 1256–1265." Journal of Molecular Liquids.

EP

1766

TE D

1762

Nekouei, F. and S. Nekouei (2017). "Comments on the paper “Adsorptive removal of methylene blue by rhamnolipid-functionalized graphene oxide from wastewater”." Water Research

1768

108: 462-463.

AC C

1767

1769

Nekouei, F., H. Kargarzadeh, S. Nekouei, I. Tyagi, S. Agarwal and V. Kumar Gupta (2016).

1770

"Preparation of Nickel hydroxide nanoplates modified activated carbon for Malachite

1771 1772

Green removal from solutions: Kinetic, thermodynamic, isotherm and antibacterial

studies." Process Safety and Environmental Protection 102: 85-97.

1773

Oleszczuk, P., B. Pan and B. Xing (2010). "Response to Comment on “Adsorption and

1774

Desorption of Oxytetracycline and Carbamazepine by Multiwalled Carbon

82

ACCEPTED MANUSCRIPT

1775 1776 1777

Nanotubes”." Environmental Science & Technology 44(12): 4829-4829. Onjia, A. E. and S. K. Milonjić (2002). Influence of the background electrolyte on the point of zero charge of chromium(III)-oxide. Materials Science Forum. 413: 87-92. Özacar, M. (2005). "Response to comment on “Two-stage batch sorber design using

1779

second-order kinetic model for the sorption of metal complex dyes onto pine sawdust”

1780

by Dr. Y.-S. Ho, Biochem. Eng. J. 21 (2004) 39–45." Biochemical Engineering Journal

1781

23(3): 293.

RI PT

1778

Park, D., Y.-S. Yun, J. Y. Kim and J. M. Park (2008). "How to study Cr(VI) biosorption: Use of

1783

fermentation waste for detoxifying Cr(VI) in aqueous solution." Chemical Engineering

1784

Journal 136(2–3): 173-179.

1786

M AN U

1785

SC

1782

Park, D., Y.-S. Yun and J. M. Park (2005). "Studies on hexavalent chromium biosorption by chemically-treated biomass of Ecklonia sp." Chemosphere 60(10): 1356-1364.

1787

Park, D., Y.-S. Yun and J. M. Park (2006a). "Comment on “Chromate ion adsorption by

1788

agricultural by-products modified with dimethyloldihydroxylethylene urea and choline

1789

chloride” by Wartelle and Marshall." Water Research 40(7): 1501-1504. Park, D., Y.-S. Yun and J. M. Park (2006b). "Comment on the Removal Mechanism of

1791

Hexavalent Chromium by Biomaterials or Biomaterial-Based Activated Carbons."

1792

Industrial & Engineering Chemistry Research 45(7): 2405-2407.

TE D

1790

Park, D., Y.-S. Yun and J. M. Park (2006c). "Mechanisms of the removal of hexavalent

1794

chromium by biomaterials or biomaterial-based activated carbons." Journal of

1795

Hazardous Materials 137(2): 1254-1257.

AC C

EP

1793

1796

Patrickios, C. S. and E. N. Yamasaki (1997). "A correction to the calculation of the Gibbs free

1797

energy of adsorption for biomolecules in ion-exchange systems." Biophysical

1798

Chemistry 69(2–3): 219-220.

1799

Plazinski, W., W. Rudzinski and A. Plazinska (2009). "Theoretical models of sorption kinetics

1800

including a surface reaction mechanism: A review." Advances in Colloid and Interface

1801

Science 152(1–2): 2-13.

83

ACCEPTED MANUSCRIPT

1802

Porter, W. R. (1985). "Comment on: "Equilibrium adsorption of polycyclic aromatic

1803

hydrocarbons from water onto activated carbon"." Environmental Science &

1804

Technology 19(9): 869-870.

1806

Radovic, L. R. (1999). "Surface chemistry of activated carbon materials: state of the art and

RI PT

1805

implications for adsorption." Surfactant science series 78: 529-565.

Rahmani-Sani, A., R.-r. Shan, L.-g. Yan and A. Hosseini-Bandegharaei (2017). "Response to

1808

“Letter to Editor: Minor correction to the thermodynamic calculation using the

1809

distribution constant by Shan et al. and Rahmani-Sani et al.”." Journal of Hazardous

1810

Materials 325(5).

1811

SC

1807

Rayne, S. (2013). "Comment on “Adsorption of ionizable organic contaminants on multi-walled carbon nanotubes with different oxygen contents [Li et al., J. Hazard. Mater. 186 (2011)

1813

407–415]”." Journal of Hazardous Materials 252–253: 1.

1814 1815

M AN U

1812

Redlich, O. and D. L. Peterson (1959). "A Useful Adsorption Isotherm." The Journal of Physical Chemistry 63(6): 1024-1024.

Rivera-Utrilla, J. and M. Sánchez-Polo (2011). "Adsorbent-adsorbate interactions in the

1817

adsorption of organic and inorganic species on ozonized activated carbons:

1818

a short review." Adsorption 17(3): 611-620.

1821 1822

biosorption." Journal of Hazardous Materials 304: 553-556.

EP

1820

Robalds, A., G. M. Naja and M. Klavins (2016). "Highlighting inconsistencies regarding metal

Roginsky, S. and Y. B. Zeldovich (1934). "The catalytic oxidation of carbon monoxide on manganese dioxide." Acta Phys. Chem. USSR 1: 554.

AC C

1819

TE D

1816

1823

Rudzinski, W. and W. Plazinski (2007). "Reply to “Comment on “Kinetics of Solute Adsorption

1824

at Solid/Solution Interfaces:  A Theoretical Development of the Empirical Pseudo-First

1825 1826 1827

and Pseudo-Second Order Kinetic Rate Equations, Based on Applying the Statistical

Rate Theory of Interfacial Transport”." The Journal of Physical Chemistry B 111(1): 319-319.

1828

Salvestrini, S., V. Leone, P. Iovino, S. Canzano and S. Capasso (2014). "Considerations about the

1829

correct evaluation of sorption thermodynamic parameters from equilibrium 84

ACCEPTED MANUSCRIPT

1830

isotherms." The Journal of Chemical Thermodynamics 68: 310-316. Senthil Kumar, P., P. S. A. Fernando, R. T. Ahmed, R. Srinath, M. Priyadharshini, A. M. Vignesh

1832

and A. Thanjiappan (2014). "Effect of temperature on the adsorption of methylene blue

1833

dye onto sulfuric acid–treated orange peel " Chemical Engineering Communications

1834

201(11): 1526-1547.

RI PT

1831

Sharma, P., N. Hussain, D. J. Borah and M. R. Das (2013). "Kinetics and Adsorption Behavior of

1836

the Methyl Blue at the Graphene Oxide/Reduced Graphene Oxide Nanosheet–Water

1837

Interface: A Comparative Study." Journal of Chemical & Engineering Data 58(12):

1838

3477-3488.

SC

1835

Sharma, Y. C., G. S. Gupta, G. Prasad and D. C. Rupainwar (1990). "Use of wollastonite in the

1840

removal of Ni(II) from aqueous solutions." Water, Air, and Soil Pollution 49(1): 69-79.

1841

Smičiklas, I., A. Onjia, S. Raičević and D. Janaćković (2009). "Authors’ response to comments

1842

on “Factors influencing the removal of divalent cations by hydroxyapatite”." Journal

1843

of Hazardous Materials 168(1): 560-562.

1845

Somasundaran, P. (1968). "Zeta potential of apatite in aqueous solutions and its change during equilibration." Journal of Colloid and Interface Science 27(4): 659-666.

TE D

1844

M AN U

1839

Swiatkowski, A., M. Pakula, S. Biniak and M. Walczyk (2004). "Influence of the surface

1847

chemistry of modified activated carbon on its electrochemical behaviour in the

1848

presence of lead(II) ions." Carbon 42(15): 3057-3069.

EP

1846

Thrower, P. (2011). "Writing a scientific paper V – Language." Carbon 49(15): 4957-4960.

1850

Thrower, P. A. (2007). "Writing a Scientific Paper: I. Titles and Abstracts." Carbon 45(11):

1851 1852 1853

AC C

1849

2143-2144.

Thrower, P. A. (2008a). "Writing a scientific paper: II. Introduction and references." Carbon

46(2): 183-184.

1854

Thrower, P. A. (2008b). "Writing a scientific paper: III. Experimental." Carbon 46(8): 1113-1114.

1855

Thrower, P. A. (2010). "Writing a scientific paper: IV. Results and discussion." Carbon 48(10):

1856

2675-2676.

85

ACCEPTED MANUSCRIPT

1858 1859 1860

Tien, C. (2007). "Remarks on adsorption manuscripts received and declined: An editorial." Separation and Purification Technology 54(3): 277-278. Tien, C. (2008). "Remarks on adsorption manuscripts revised and declined: An editorial." Journal of Hazardous Materials 150(1): 2-3.

RI PT

1857

1861

Tran, H. N. (2017). "Comments on “Characterization and adsorption capacity of raw

1862

pomegranate peel biosorbent for copper removal”." Journal of Cleaner Production 144:

1863

553-558.

Tran, H. N., F.-C. Huang, C.-K. Lee and H.-P. Chao (2017a). Activated carbon derived from

1865

spherical

functionalized

with

Triethylenetetramine:

synthesis,

1866

characterizations, and adsorption application. Green Processing and Synthesis.

M AN U

hydrochar

SC

1864

1867

Tran, H. N., Y.-F. Wang, S.-J. You and H.-P. Chao (2017b). "Insights into the mechanism of

1868

cationic dye adsorption on activated charcoal: the importance of π-π interactions."

1869

Process Safety and Environmental Protection.

Tran, H. N., S.-J. You and H.-P. Chao (2015). "Effect of pyrolysis temperatures and times on the

1871

adsorption of cadmium onto orange peel derived biochar." Waste Management &

1872

Research.

TE D

1870

Tran, H. N., S.-J. You and H.-P. Chao (2016). "Thermodynamic parameters of cadmium

1874

adsorption onto orange peel calculated from various methods: A comparison study."

1875

Journal of Environmental Chemical Engineering 4(3): 2671-2682.

EP

1873

Tran, H. N., S.-J. You and H.-P. Chao (2017c). "Activated Carbons from Golden Shower upon

1877

Different Chemical Activation Methods: Synthesis and Characterizations." Adsorption

1878

AC C

1876

Science & Technology.

1879

Tran, H. N., S.-J. You and H.-P. Chao (2017d). "Fast and efficient adsorption of methylene green

1880

5 on activated carbon prepared from new chemical activation method." Journal of

1881

Environmental Management 188: 322-336.

1882

Tran, H. N., S.-J. You and H.-P. Chao (2017e). "Insight into adsorption mechanism of cationic

1883

dye onto agricultural residues-derived hyrochars: negligible role of π-π interaction."

1884

Korean Journal of Chemical Engineering. 86

ACCEPTED MANUSCRIPT

Uchimiya, M., I. M. Lima, K. Thomas Klasson, S. Chang, L. H. Wartelle and J. E. Rodgers

1886

(2010). "Immobilization of Heavy Metal Ions (CuII, CdII, NiII, and PbII) by Broiler

1887

Litter-Derived Biochars in Water and Soil." Journal of Agricultural and Food

1888

Chemistry 58(9): 5538-5544.

RI PT

1885

1889

USEPA (1999). Understanding variation in partition coefficient, Kd, values. United States Office

1890

of Air and Radiation. EPA 402-R-99-004A. Washington, USA, Environmental

1891

Protection Agency. 1.

Utomo, R. P., R. Koswojo, Y.-H. Ju, A. Ayucitra, F. E. Soetaredjo, J. Sunarso and S. Ismadji

1893

(2010). "Reply to the comment on "Acid Green 25 removal from wastewater by

1894

organo-bentonite from Pacitan" by R. Koswojo, R. P. Utomo, Y.-H. Ju, A. Ayucitra, F.

1895

E. Soetaredjo, J. Sunarso, S. Ismadji [Applied Clay Science 48 (2010) 81–86]."

1896

Applied Clay Science 50(1): 165-166.

M AN U

SC

1892

1897

van Noort, P. C. M., M. T. O. Jonker and A. A. Koelmans (2005). "Response to Comment on

1898

“Modeling Maximum Adsorption Capacities of Soot and Soot-like Materials for PAHs

1899

and PCBs”." Environmental Science & Technology 39(1): 383-384. Vasanth Kumar, K., K. Porkodi and F. Rocha (2007). "Comments on “Equilibrium and kinetic

1901

studies for the biosorption system of copper(II) ion from aqueous solution using

1902

Tectona grandis L.f. leaves powder”." Journal of Hazardous Materials 146(1–2):

1903

428-429.

TE D

1900

Vasanth Kumar, K. and S. Rattanaphani (2007). "Reply to ‘Comments on “An adsorption and

1905

kinetic study of lac dyeing on silk”’ by Yuh-Shan Ho: Discussion on pseudo second

1906

order kinetic expression." Dyes and Pigments 75(1): 253-254.

AC C

EP

1904

1907

Vasanth Kumar, K. and S. Sivanesan (2006). "Equilibrium data, isotherm parameters and process

1908

design for partial and complete isotherm of methylene blue onto activated carbon."

1909

Journal of Hazardous Materials 134(1–3): 237-244.

1910

Volesky, B. (2003). Sorption and biosorption, BV Sorbex.

1911

Volesky, B. (2007). "Biosorption and me." Water Research 41(18): 4017-4029.

1912

Walter, W. J. (1984). "Evolution of a Technology." Journal of Environmental Engineering 110(5): 87

ACCEPTED MANUSCRIPT

1913

899-917. Wang, H., X. Yuan, Y. Wu, H. Huang, G. Zeng, Y. Liu, X. Wang, N. Lin and Y. Qi (2014).

1915

"Corrigendum to “Adsorption characteristics and behaviors of graphene oxide for

1916

Zn(II) removal from aqueous solution” [Appl. Surf. Sci. 279 (2013) 432–440]."

1917

Applied Surface Science 301: 585.

RI PT

1914

1918

Wang, L. K. (2009). Heavy metals in the environment. Boca Raton, CRC Press.

1919

Weber, T. W. and R. K. Chakravorti (1974). "Pore and solid diffusion models for fixed-bed

1922 1923 1924 1925

SC

1921

adsorbers." AIChE Journal 20(2): 228-238.

Weber, W. J. and J. C. Morris (1963). "Kinetics of adsorption on carbon from solution." Journal of the Sanitary Engineering Division 89(2): 31-60.

M AN U

1920

Weber, W. J. and E. H. Smith (1987). "Simulation and design models for adsorption processes." Environmental Science & Technology 21(11): 1040-1050. Werner, D. and H. K. Karapanagioti (2005). "Comment on “Modeling Maximum Adsorption Capacities of Soot and Soot-like Materials for PAHs and PCBs”." Environmental

1927

Science & Technology 39(1): 381-382.

1928 1929

TE D

1926

Worch, E. (2012). Adsorption technology in water treatment: fundamentals, processes, and modeling, Walter de Gruyter.

Wu, C.-H. (2009). "Response to comment on “Adsorption of direct dyes from aqueous solutions

1931

by carbon nanotubes: Determination of equilibrium, kinetics and thermodynamics

1932

parameter” by Dr. Ho." Journal of Colloid and Interface Science 333(1): 413-414.

1934 1935

Wu, F.-C., B.-L. Liu, K.-T. Wu and R.-L. Tseng (2010). "A new linear form analysis of Redlich–

AC C

1933

EP

1930

Peterson isotherm equation for the adsorptions of dyes." Chemical Engineering

Journal 162(1): 21-27.

1936

Wu, Z., H. Zhong, X. Yuan, H. Wang, L. Wang, X. Chen, G. Zeng and Y. Wu (2017). "Reply for

1937

comment on “Adsorptive removal of methylene blue by rhamnolipid-functionalized

1938 1939

graphene oxide from wastewater”." Water Research 108: 464-465. Zhao, Y. (2015). "Comment on “Zirconium–Carbon Hybrid Sorbent for Removal of Fluoride

88

ACCEPTED MANUSCRIPT

1940

from Water: Oxalic Acid Mediated Zr(IV) Assembly and Adsorption Mechanism”."

1941

Environmental Science & Technology 49(19): 11982-11983. Zhou, Q., C. Xie, W. Gong, N. Xu and W. Zhou (2011). "Comments on the method of using

1943

maximum absorption wavelength to calculate Congo Red solution concentration

1944

published in J. Hazard. Mater." Journal of Hazardous Materials 198: 381-382.

1945

Zhou, X. (2017). "Minor correction to the thermodynamic calculation using the distribution

1946

constant by Shan et al. and Rahmani-Sani et al." Journal of Hazardous Materials 323,

1947

Part B: 735-736.

SC

1949

Zhou, X., H. Liu and J. Hao (2012). "Letters to the Editor." Adsorption Science & Technology

30(7): 647-649.

M AN U

1948

RI PT

1942

1950

Zhou, X. and X. Zhou (2014). "The unit problem in the thermodynamic calculation of adsorption

1951

using the Langmuir eqaution." Chemical Engineering Communications 201(11):

1952

1459-1467.

1953

AC C

EP

TE D

1954

89

ACCEPTED MANUSCRIPT

Table captions Table 1. Relationship between isotherm parameters and isotherm shape

RI PT

Table 2. Mistakes in the presentation and calculation of the separation factor (RL) Table 3. Incorrect values of the exponent g in the Redlich–Peterson equation

Table 4. Comparison of incorrect and correct thermodynamic parameters for methylene blue

SC

adsorption onto sulfuric acid-treated orange peel

Table 5. A typical example of correct and incorrect use of a Micromeritics report (data not

AC C

EP

TE D

M AN U

published)

ACCEPTED MANUSCRIPT

Table 1. Relationship between isotherm parameters and isotherm shape

AC C

EP

TE D

M AN U

SC

RI PT

Freundlich exponent Separation factor Isotherm shapes Remarks n=0 RL = 0 Irreversible Horizontal n<1 RL < 1 Favorable Concave n=1 RL = 1 Linear Linear n>1 RL > 1 Unfavorable Convex Note: Information was adapted from Worch (2012). Copyright (2012), with permission from De Gruyter. RL and n were calculated from equations 52 and 56, respectively.

ACCEPTED MANUSCRIPT

Table 2. Mistakes in the presentation and calculation of the separation factor (RL)

AC C

EP

TE D

M AN U

SC

RI PT

Typical examples Remakes Langmuir parameters The RL is not the Langmuir o Q max KL RL constant, so it is mistake for Orange peel 54.5 0.101 0.01–0.17 presentation. OP Biochar 115 0.019 0.08–0.57 Langmuir parameters The RL values must be calculated RL o Q max KL at different range of initial adsorbate concentration, so it is Orange peel 54.5 0.101 0.17 mistake for calculation. OP Biochar 115 0.019 0.57 Langmuir parameters RL Qomax KL Correction presentation and calculation. Orange peel 54.5 0.101 0.01–0.17 OP Biochar 115 0.019 0.08–0.57 Note: The data were modified from our previous publication (Tran, You et al. 2015).

ACCEPTED MANUSCRIPT

Table 3. Incorrect values of the exponent g in the Redlich–Peterson equation Adsorbate

Sawdust

Reactive black

1.10

(Inbaraj, 2006)

Sawdust

Reactive red

1.04

(Inbaraj, 2006)

Congo red

2.91

(Kumar and Porkodi, 2008)

Xad-4 resin

4-chlorophenol at 303 K

1.07

(Kumar and Porkodi, 2007)

Xad-4 resin

4-chlorophenol at 318 K

1.19

(Kumar and Porkodi, 2007)

2+

–0.35

(Vasanth Kumar et al., 2007)

2+

–1.23

(Kumar et al., 2007

Tectona grandis

Cu

Syzygium cumini

Pb

Activated carbon/ chitoshan composite Activated carbon/ chitoshan composite Activated carbon

Aniline/aniline + Cr(VI)

1.04

(Huang et al., 2014)

Cr(VI)/aniline + Cr(VI)

1.04

(Huang et al. 2014)

Cr(VI)

1.02

(Huang et al. 2014)

AC C

EP

TE D

/chitoshan composite

SC

xerogel

Reference

M AN U

Anilinepropylsilica

g

RI PT

Adsorbent

ACCEPTED MANUSCRIPT

Table 4. Comparison of incorrect and correct thermodynamic parameters for methylene blue adsorption onto sulfuric acid-treated orange peel Incorrect thermodynamic parameters reported by Senthil Kumar et al. (2014)

50 100 150 200 250

Temperature

30 °C –12.5 –10.1 –8.86 –8.03 –7.29

RI PT

40 °C 50 °C 60 °C –53.7 –137 –10.4 –8.87 –8.49 –37.4 –90.3 –9.15 –8.04 –7.46 –32.9 –79.4 –8.19 –7.27 –6.52 –27.5 –64.2 –7.54 –6.64 –6.19 –23.7 –54.4 –6.61 –6.17 –5.64 KC values 30 °C 40 °C 50 °C 60 °C 142 54.9 27.2 21.5 54.8 33.6 19.9 14.8 33.7 23.3 14.9 10.5 24.2 18.1 11.9 9.35 18.1 12.7 9.97 7.66 Correct thermodynamic parameters suggested by Hai (2017) ∆Go ∆Ho ∆So KC van't Hoff equation (kJ/mol) (kJ/mol) (J/mol) 5.25

40 °C

4.33

50 °C

3.54

60 °C

3.27

-4.18 -3.81 -3.40

EP

30 °C

13.65

-31.39 y = 1642 – 3.78 R2 = 0.978

-3.28

The unit of J/mol/K was expressed as the original paper; the correction unit might be kJ/mol.

AC C

a

∆Go (kJ/mol)

SC

Co (mg/L)

∆So (J/mol/K)a

M AN U

50 100 150 200 250

∆Ho (kJ/mol)

TE D

Co (mg/L)

ACCEPTED MANUSCRIPT

Table 5. A typical example of correct and incorrect use of a Micromeritics report Correct data

Prepared activated carbon

Prepared Biochar

BET surface area

613

536

Langmuir surface area

1684

Micropore area

NA

External surface area

698

RI PT

Incorrect data

660

348

188

AC C

EP

TE D

M AN U

SC

NA: The micropore area is not reported because either the micropore volume is negative or the calculation of external surface area is larger than the total surface area (BET surface area). The unit is m2/g. Data is not published.

ACCEPTED MANUSCRIPT

Figure captions Figure 1. Some basic terms used in adsorption science and technology

RI PT

Figure 2. Basic properties of an adsorbent determined by various common techniques Figure 3. Schematic illustration of a new classification system for metal biosorption mechanisms

SC

Figure 4. A typical example of difference between (a) point of zero charge and (b) zeta potential of glucose-derived spherical biochar prepared from 800 °C

M AN U

Figure 5. Effect of contact time on the adsorption process of (a) methyl violet onto granulated mesoporous carbon (GMC) and granular activated carbon (GAC); (b) methyl tert-butyl ether onto carbonaceous resin (Ambersorb 563); (c) ibuprofen (IBP), carbamazepine (CBZ), ofloxacin (OFX), bisphenol-A (BPA), diclofenac (DFN), mecoprop (MCP), pentachlorophenol (PCP), benzotriazol (BZT), and caffeine

TE D

(CAF) onto activated carbon cloths (resin phenolic precursor); and (d) methyl blue onto graphene oxide

Figure 6. Fits of kinetic adsorption data to the linear pseudo-second-order equation (the study

EP

of methylene green 5 adsorption onto commercial activated carbon (CAC) and synthesized activated carbons prepared from golden shower through different

AC C

chemical activation methods (GSAC, GSBAC, GSHAC, and GSHBAC))

Figure 7. Transport processes during adsorption by a porous adsorbent Figure 8. (a) non-equilibrium adsorption isotherms of methylene blue onto sulfuric acid-treated orange peel (STOP) at various temperatures, and (b) comparison of adsorption isotherm curves; adsorbent B performs better (higher qe at q10) than adsorbent A at lower equilibrium concentrations (e.g., Ce = 10 mg/L)

ACCEPTED MANUSCRIPT

Figure 9. Proposed mechanisms for Cr(VI) biosorption by nonliving biomass Figure 10. (a) Chemical structure of Congo Red and (b) UV-vis spectra of Congo Red

RI PT

solutions at different solution pH values Figure 11. An example of a graphene layer and proposed mechanisms of methylene adsorption green 5 onto biochar, synthesized activated carbon, and commercial

SC

activated charcoal

Figure 12. A simple guide to the nonlinear method of calculating the Langmuir parameters

AC C

EP

TE D

M AN U

(RSS: residual sum of squares, TSS: total sum of squares)

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

AC C

EP

Figure 1. Some basic terms used in adsorption science and technology

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

AC C

EP

TE D

Figure 2. Basic properties of an adsorbent determined by various common techniques

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Figure 3. Schematic illustration of a new classification system for metal (bio)sorption mechanisms

AC C

EP

TE D

(Adapted from Robalds et al., 2016. Copyright (2016), with permission from Elsevier)

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Figure 4. A typical example of difference between (a) point of zero charge and (b) zeta

AC C

EP

potential of glucose-derived spherical biochar prepared from 800 °C (Data unpublished)

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Figure 5. Effect of contact time on the adsorption process of (a) methyl violet onto granulated mesoporous carbon (GMC) and granular activated carbon (GAC) (Reprinted from Kim et al.

EP

2016. Copyright (2016), with permission from Elsevier); (b) methyl tert-butyl ether onto carbonaceous resin (Ambersorb 563) (Reprinted from Hung and Lin 2006. Copyright (2006), with permission from Elsevier); (c) ibuprofen (IBP), carbamazepine (CBZ), ofloxacin (OFX),

AC C

bisphenol-A (BPA), diclofenac (DFN), mecoprop (MCP), pentachlorophenol (PCP), benzotriazol (BZT), and caffeine (CAF) onto activated carbon cloths (resin phenolic precursor) (Reprinted from Masson et al. 2016. Copyright (2016) with permission from Elsevier); and (d) methyl blue onto graphene oxide (Reprinted with permission from Sharma et al. 2013. Copyright (2013) American Chemical Society)

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Figure 6. Fits of kinetic adsorption data to the linear pseudo-second-order equation (the study of methylene green 5 adsorption onto commercial activated carbon (CAC) and synthesized

EP

activated carbons prepared from golden shower through different chemical activation methods (GSAC, GSBAC, GSHAC, and GSHBAC)) (Adapetd from Tran et al., 2017d.

AC C

Copyright (2017), with permission from Elsevier)

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Figure 7. Transport processes during adsorption by a porous adsorbent (Adapted with

AC C

EP

permission from Weber and Smith, 1987. Copyright (1987) American Chemical Society)

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Figure 8. (a) non-equilibrium adsorption isotherms of methylene blue onto sulfuric acid-treated orange peel (STOP) at various temperatures (Adapted from Hai 2016. Copyright (2016) Taylor and Francis), and (b) comparison of adsorption isotherm curves: adsorbent B performs better (higher qe at q10) than adsorbent A at lower equilibrium concentrations (e.g., Ce = 10 mg/L) (Adapted from Volesky 2007. Copyright (2007), with permission from Elsevier)

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

TE D

Figure 9. Proposed mechanisms for Cr(VI) biosorption by nonliving biomass

AC C

EP

(Adapted from Park et al., 2005. Copyright (2005), with permission from Elsevier).

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Figure 10. (a) Chemical structure of Congo Red and (b) UV-vis spectra of Congo Red solutions at different solution pH values (Adapted from Zhou et al., 2011. Copyright (2011), with permission from Elsevier)

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Figure 11. An example of a graphene layer and proposed mechanisms of methylene

EP

adsorption green 5 onto biochar, synthesized activated carbon, and commercial activated

AC C

charcoal (Adapetd from Tran et al., 2017 b,d. Copyright (2017), with permission from Elsevier)

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

Figure 12. A simple guide to the nonlinear method of calculating the Langmuir parameters (RSS:

AC C

EP

TE D

residual sum of squares, TSS: total sum of squares)

ACCEPTED MANUSCRIPT

Highlights • This study summarizes common literature mistakes in the field of

RI PT

adsorption • Reiteration and propagation of such mistakes in future publications should be avoided.

SC

• Correct expression and citation of the models used in adsorption studies are provided

AC C

EP

TE D

further investigated

M AN U

• We also highlight some problems that need to be thoroughly discussed or