Chemistry - A European Journal
10.1002/chem.201800623
FULL PAPER
each polymer by adjusting the angle α and the height c (Figure S12) so
that layers of sheets directly on top of one another (AA) or offset by half
of the repeating unit (AB) were created. Specific parameters used in
each unit cell can be found in Table S1 and unit cells visualized in MAPS
can be found in Figure S11. The net charge of COP1 and COP2 unit cell
was +4. To neutralize the polymer, we added 4 Cl– counter ions in
random positions into the simulation box. Geometry optimization and MD
simulations were conducted using LAMMPS[53] program with
DreidingX6[54] force field. Equations of states were integrated using a
velocity-Verlet[55] integrator with a timestep of 2 fs in the isothermal-
isobaric (NPT) ensemble for each model. A temperature of 298.15 K and
a pressure of 101.325 kPa (1 atm) were maintained using a Nose-Hoover
thermostat[56] and barostat,[57] respectively. Coulomb interactions were
Congo red. The data were fitted to Ho and McKay’s pseudo second-order
adsorption model t/q
t
= t/q
e
+ 1/(kobs·q
e
2), where q
t
is the quantity of
adsorbed Congo red at time t (min), q
e
is the adsorbed amount at
equilibrium, and kobs is the second-order rate constant (g mg-1 min-1).[59]
To determine maximum adsorption capacity, adsorption isotherms were
obtained by varying the initial concentration of the dye from 0.1 mM to
.5 mM. The adsorption isotherms were fitted (correlation coefficients, R2
1
>
0.96) by using the Langmuir adsorption model: q
e
= q
(mg L-
(mg g ) is the amount of dye adsorbed at equilibrium, C
e
m e e
·b·C /(1+b·C ),
-1
where q
1)
e
is the equilibrium solute concentration remaining in solution when q
e
is
achieved, q
m
is the maximum adsorption capacity corresponding to
-1 [14]
complete monolayer coverage, and b is a constant (L mg ).
[
58]
evaluated using an implementation of the Ewald summation algorithm
-
4
with a relative force accuracy of 1.0 × 10 . Lennard-Jones and Coulomb
interactions were cut off at a distance of 12 Å. Tail corrections were
employed for accurate computation of the long-range interactions. Each
structure was left to sample conformations for 100 ns to ascertain its
stability and to visualize and analyze the equilibrium populations. Data
were recorded with 1 ps intervals with the first 10 ns considered
equilibration and discarded. All observables reported here were obtained
by time series averaging of the snapshots.
Acknowledgements
This work was supported by New York University Abu Dhabi.
We thank NYUAD and the Core Technology Platforms for their
generous support for the research program. This research was
carried out on the High Performance Computing resources at
New York University Abu Dhabi. The authors also thank Dr.
Sajeewa Walimuni Degawe for his help with the MAPS software
and Dr. Thirumurugan Prakasam for his help with EPR
spectroscopy.
To compare the stability of alternative sheet configurations we computed
the average enthalpy ∆H =(1/푡) ∑ tt =1 ∆H , where ∆Hi is the instantaneous
i
enthalpy of snapshot i, and t is the total number of snapshots after
equilibrium is reached. The entropy of inter-layer state is estimated from
the trajectory using Gibbs Entropy ∆S = - k ∫ P(E)ln(P(E))dE , where kB
B
is the Boltzmann constant, and P(E) is the probability density of total
energy, E sampled during the simulations. Using these two entities we
compute the Gibbs free energy as ∆G = ∆H - T∆S.
Keywords: viologen • calix[4]arene • nanosheet • iodine capture
•
Congo red
Iodine capture experiments. For vapor phase experiments, ~10 mg of
the polymer was placed on a small crystallization dish, which was placed
[
1]
H. Li, D.-H. Qu, Sci. China Chem. 2015, 58, 916–921.
in
a
closed chamber containing solid I
2
and incubated at 70 °C.
[2]
[3]
C. L. Bird, A. T. Kuhn, Chem. Soc. Rev. 1981, 10, 49–82.
C. Hua, B. Chan, A. Rawal, F. Tuna, D. Collison, J. M. Hook, D. M.
D’Alessandro, J. Mater. Chem. C 2016, 4, 2535–2544.
O. Buyukcakir, S. H. Je, S. N. Talapaneni, D. Kim, A. Coskun, ACS
Appl. Mater. Interfaces 2017, 9, 7209–7216.
Gravimetric measurements of the polymer were recorded at regular time
intervals up to 24 h, and the % mass increase was calculated as follows:
%
mass increase at [t=t ]= (mass at [t=t ]-mass at [t=0])⁄mass at [t=0].
i i
[
4]
2
For solution phase experiments, 5 mL of 1 mM stock solution of I was
prepared in cyclohexane. 5.0 mg of polymer was added to the stock
solution, and the mixture was continuously stirred. The progress of iodine
capture was monitored by UV-Vis spectrophotometry at regular time
intervals up to 210 min. % iodine uptake was calculated by dividing the
absorbance at λmax = 523 nm at a given time by the absorbance at λmax of
[5]
G. Das, T. Skorjanc, S. K. Sharma, F. Gándara, M. Lusi, D. S.
Shankar Rao, S. Vimala, S. Krishna Prasad, J. Raya, D. S. Han, et
al., J. Am. Chem. Soc. 2017, 139, 9558–9565.
[
[
6]
7]
H. Yen, C. Tsai, S. Chen, G. Liou, Macromol. Rapid Commun. 2017,
2
the starting I solution.
38, 1–7.
M. Burgess, E. Che, R. S. Assary, J. Hui, S. Moore, Chem. Mater.
2016, 28, 7362–7374.
Congo red adsorption experiments. To determine the effect of pH on
Congo red uptake, 25 μM aqueous dye solutions of various pH were
prepared. Solution pH was adjusted with 1 M NaOH or 1 M HCl solutions.
UV-Vis spectra of the starting solutions were measured. Then, 5.0 mg of
polymer was added to 5 mL of each solution, and the mixture was
continuously stirred. The progress of dye uptake was monitored by UV-
Vis spectrophotometry at regular time intervals up to 375 min. Percent
Congo red uptake was calculated by dividing the absorbance of the dye
solution at λmax after 375 min by the absorbance at λmax before polymer
addition and multiplying by 100. Uptake was calculated for each polymer,
and for each of the three redox states of its constituent viologens. To
compare the different redox states of COPs in their ability to adsorb
Congo red, the % uptake after 15 min was considered.
[8]
H. Ma, B. Liu, B. Li, L. Zhang, Y.-G. Li, H.-Q. Tan, H.-Y. Zang, G.
Zhu, J. Am. Chem. Soc. 2016, 138, 5897–5903.
[
9]
G. Das, T. Skorjanc, T. Prakasam, S. Nuryyeva, J.-C. Olsen, A.
Trabolsi, RSC Adv. 2017, 7, 3594–3598.
[10]
T. Skorjanc, F. Benyettou, J.-C. Olsen, A. Trabolsi, Chem. – A Eur.
J. 2017, 23, 8333–8347.
[
[
11]
12]
L. Wang, X. Shi, P. Jia, Y. Yang, J. Polym. Sci. Part A Polym.
Chem. 2004, 42, 6259–6266.
S. Memon, O. Oguz, A. Yilmaz, M. Tabakci, M. Yilmaz, Ş. Ertul, J.
Polym. Environ. 2001, 9, 97–101.
For investigation of the kinetics, the amount of pollutant adsorbed at each
[13]
[14]
Y. Yang, T. M. Swager, Macromolecules 2007, 40, 7437–7440.
D. Shetty, I. Jahovic, J. Raya, F. Ravaux, M. Jouiad, J.-C. Olsen, A.
Trabolsi, J. Mater. Chem. A 2017, 5, 62–66.
time point was determined using the equation q
where q (mg/g) is the quantity of adsorbed Congo red per 1 g of
adsorbent at a given time t (min), C (mmol/L) is the initial concentration
before polymer treatment, C (mmol/L) is the residual concentration of the
dye after treatment with polymer, V is the volume of dye solution, m (g) is
the amount of adsorbent used, and M (g/mol) is the molecular weight of
t 0 t w
= [(C – C )·M ·V]/m,
t
0
[
15]
G. Das, T. Prakasam, S. Nuryyeva, D. S. Han, A. Abdel-Wahab, J.-
C. Olsen, K. Polychronopoulou, C. Platas-Iglesias, F. Ravaux, M.
Jouiad, et al., J. Mater. Chem. A 2016, 4, 15361–15369.
t
w
This article is protected by copyright. All rights reserved.