Cellulose Nanofiber - Graphene Oxide Biohybrids ... - ACS Publications

3 days ago - This was already apparent from the simulation of the ions in aqueous solution where small mixed aggregates made, at most, of five units w...
0 downloads 0 Views 3MB Size
Subscriber access provided by Kaohsiung Medical University

Article

Cellulose Nanofiber - Graphene Oxide Biohybrids: Disclosing the Self-Assembly and Copper ion Adsorption using Advanced Microscopy and ReaxFF Simulations Chuantao Zhu, Susanna Monti, and Aji P Mathew ACS Nano, Just Accepted Manuscript • DOI: 10.1021/acsnano.8b02734 • Publication Date (Web): 11 Jun 2018 Downloaded from http://pubs.acs.org on June 11, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Cellulose Nanofiber - Graphene Oxide Biohybrids: Disclosing the Self-Assembly and Copper ion Adsorption using Advanced Microscopy and ReaxFF Simulations Chuantao Zhu†, Susanna Monti‡ and Aji P Mathew†*

† Division of Materials and Environmental Chemistry, Stockholm University, Stockholm, 10691, Sweden ‡ CNR - Institute of Chemistry of Organometallic Compounds, Area della Ricerca, via Moruzzi 1, 56124 Pisa, Italy

ABSTRACT: The self-assembly of nanocellulose and graphene oxide into highly porous biohybrid materials has inspired the design and synthesis of multifunctional membranes for removing water pollutants. The mechanisms of self-assembly, metal ion capture and cluster formation on the biohybrids at the nano- and molecular scales are quite complex. Their elucidation requires evidence from the synergistic combination of experimental data and computational models. The AFM-based microscopy studies of (2,2,6,6-tetrame-thylpiperidine-1-oxylradical) (TEMPO)mediated oxidized cellulose nanofibers (TOCNF), graphene oxide (GO) and their biohybrid membranes provide strong, direct evidence of self-assembly; small GO nanoparticles first attach and accumulate along a single TOCNF fiber, while the long, flexible TOCNF filaments wrap around the flat, wide GO planes, thus forming an amorphous and porous biohybrid network. The layered structure of the TOCNF and GO membrane, derived from the self-assembly and its surface properties before and after the adsorption of Cu(II), is investigated by advanced microscopy techniques and is further clarified by the ReaxFF molecular dynamics (MD) simulations. The dynamics of the Cu(II)-ion capture by the TOCNF and GO membranes in solution and the ion cluster formation during drying are confirmed by the MD simulations. The

ACS Paragon Plus Environment

1

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 26

results of this multi-disciplinary investigation move the research one step forward by disclosing specific aspects of the self-assembly behavior of bio-species and suggesting effective design strategies to control the pore size and robust materials for industrial applications. KEYWORDS: TEMPO cellulose nanofibers, graphene oxide, biohybrids, self-assembly, atomic force microscopy, ReaxFF, molecular modeling

Over the past few years, nanocellulose- and graphene-based materials have attracted significant attention in numerous intriguing applications as adsorbents and separation materials.

1-3

The

typical derivatives of cellulose and graphene, namely, (2,2,6,6-tetrame-thylpiperidine-1oxylradical) (TEMPO)-mediated oxidized cellulose nanofibers (TOCNF) 4,5 and graphene oxide (GO6,7, are even considered among the most promising “capturing agents”, and this status is why they have been selected as excellent candidates for fabricating sustainable adsorbents, such as biohybrid aerogels, 8 smart composites9 and eco-friendly membranes. 10,11 As a matter of fact, TOCNF, modified by the position-selective catalytic oxidation of the C6 primary hydroxyls from nanocellulose, possess active-surface characteristics due to the presence of many carboxyl and hydroxyl moieties. These are responsible for the outstanding adsorbent capacity of TOCNF that is actuated through the coordination and capture of heavy metal ions and organic pollutants. 12-14 Similarly, to TOCNF, the GO sheets with hydroxyl, carboxyl and epoxy groups decorating the basal planes and edges15 display large surface areas, high surface activity, good hydrophilicity, partial hydrophobicity in the graphene domains, as well as antifouling properties.

15-17

Practically, the GO sheets are amphiphilic systems whose action,

when adsorbed on surfaces, depends on the morphology of the sheets, the binding strength and the multiscale assembly. The synergistic combination of TOCNF and GO, which is realized through favorable interactions between the abundant hydrogen donors and acceptors, generates an active porous self-assembled biohybrid structure. 18,19 These membranes are very promising materials not only for their adsorption efficiency and high capacity of retention of heavy metals and organic dyes present in solution but also for their remarkable mechanical properties and recyclability.

19

All these properties are strongly related to the conformation, supramolecular

organization, self-assembly and combination of the various components, which are all driven by intermolecular interactions. These can be, for example, intra- and inter-molecular hydrogen bonds, van der Waals and electrostatic forces. Although several experimental investigations have tried to disclose the molecular organization of the various components inside the matrices and

ACS Paragon Plus Environment

2

Page 3 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

the intermolecular interactions responsible for aggregation, adsorption and entrapment, through different techniques, 19,20 the collected data are not detailed enough to give a clear picture of all these behaviors, and many questions still remain unanswered. This status is due to the structural and morphological complexity of the hybrid membranes, which can, for example, consist of crystalline and amorphous domains randomly organized in varying proportions and to the dispersed/attached components. 21 In a previous paper,

19

we have investigated experimentally the self-assembly of TOCNF-,

GO-, and nanoGO-based biohybrids disclosing their unique porous structures before and after adsorption of copper ions. The nanoGO sheets acted as linking agents improving the porosity of the material and the water flow, whereas functionalized cellulose played both roles of crosslinker and adsorbent. We demonstrated that these hybrids combinations are efficient adsorbents that could be tuned by modulating their adsorption capacity, flexibility,

wet stability and

mechanical properties. Further more, surface adsorption of copper ions and formation of nanoparticles of copper on TOCNF was discussed earlier 20 evidencing the increasing capacity of the ions capture when the TOCNF carboxylate content was increased. We have earlier shown through AFM, XPS, and ATR-FTIR spectroscopies, the enhanced adsorption of cellulose with monovalent silver ions cellulose surface by means of sulfate and phosphate groups and computational studies confirmed the possibility of cluster formation and the stability of the silver clusters. As already evidenced in a number of studies, computational methods such as all-atom classical reactive molecular dynamics simulations (RMD) can be very useful to interpret and predict the experimental results. 22-30 Indeed, they can follow the evolution of relatively large portions of the materials for relatively long times, realistically reproduce the system’s behavior in elaborate environments by mimicking the experimental set-up during all stages of metal/pollutant adsorption, and simulate bond breakage and formation under specific circumstances created by local effects. Given these premises, in this work, we combined the experimental characterization and used advanced microscopy techniques with simulations based on force fields to scrutinize the ability of TOCNF and GO at the atomic level to assemble and cooperatively capture copper ions in solution. The modeling was carried out by building a series of representative supramolecular models, inspired by the data available in the literature or kindly provided by the authors (see computational details and related references), and by exploring their dynamics through the state-of-the-art RMD simulations. The RMD runs were based on the ReaxFF approach. 31 The force field employed was appropriately parameterized to describe these systems and was already used effectively in previous investigations dealing with cellulose, graphene,

ACS Paragon Plus Environment

3

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 26

graphene oxide and metal ions in various environmental conditions. 22-30 This study integrates experiments and simulations in each phase of the sample treatment, namely, the preparation of TOCNF- and GO-based biohybrids (100:1, wt %), characterization of self-assemblies, depiction of the membrane structures and evaluations of surface properties due to copper-ion adsorption. The results are satisfactorily consistent and support the importance of a multidisciplinary approach to address the challenging task of tailoring the structural features of bio-filters to enhance the metal ion uptake and organic pollutant trapping.

RESULTS AND DISCUSSION Self-assembly of TOCNF and GO biohybrids. The systematic investigation of morphologies and sizes of TOCNF and GO sheets/nanoparticles is crucial to the studies of the biohybrid selfassembly. They were investigated by AFM imaging and are shown in Figure S1. Dispersed TOCNF are visible in Figure S1a. The diameter of TOCNF is 2-3 nm, with the length from tens of nanometers to several micrometers. The GO sheets and nanoparticles usually have broad sizedistributions due to the random chemical exfoliation process in the GO mass production. 32,33 In the present case, the GO size varies at a wide range, as shown from the AFM imaging in Figure S1b,c (supporting information). Considering that the thickness of one layer of GO sheets is approximately 1 nm, the stacking sheets can reach the thickness of several nanometers (Figure S1b). Figure S1c shows some randomly dispersed GO nanoparticles, with the size varying from 1 nm to 30 nm, indicating that the GO suspension contains a wide size range of both GO sheets and nanoparticles. The self-assembly between TOCNF and GO is demonstrated in Figure 1, which exhibits a homogeneous network structure of TOCNF and GO biohybrids. The GO sheets with a thickness of approximately 1 nm either lie flat or embed in the TOCNF-fiber skeleton (Figure 1a,b). The big GO nanosheets were either found covering the TOCNF network or partially embedded in it (see Figure 1a, circled area). The self-assemblies of both GO sheets and nanoparticles with TOCNF are consistent with our early report.

19

Moreover, Figure 1c provides direct, strong

evidence of the self-assembly between TOCNF and smaller GO nanoparticles (several nanometers), which, as far as we know, has never been detected and reported before. Small GO nanoparticles with the 2-3-nm sizes are attached and accumulated along the TOCNF fibers (Figure 1c). This assembly is not only due to the physical interaction between the components in the hybrids but also due to effective chemical interactions/crosslinking between the various components.

ACS Paragon Plus Environment

4

Page 5 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

The RMD simulations confirm the tendency of GO to self-assemble both in solution and when in contact with the TOCNF fibers (Figure 1d,e). The interactions are essentially the intermolecular hydrogen bonds between the hydroxyl and carboxyl moieties at the edges of the sheets or functionalities of the basal planes, which determine the formation of stacked structures and nanoparticles, depending on the degree of basal plane bending. Furthermore, the simulations suggested that the type of aggregation was mainly due to the local concentration of the molecules and to the interference of the surrounding species. In fact, during the simulations, it was noticed that even though the initial positions of the GO molecules were almost parallel to the fibril surface, some of them were displaced from their flat orientation by the solvent action and by the attraction of the nearby fibers. The resulting reorganization was a rotation of the whole GO sheets and their extension toward the adjacent TOCNF fibril. Essentially, GO acted as a crosslinker and induced the formation of large fibril bundles. The bonding between neighboring GO sheets could be enhanced by TOCNF because of the accessible active groups on the fibril surface that provide extra bonding options. 23 A single TOCNF could also interact with the functional groups and lie on the flat surface of large GO sheets forming layered structures. Thus, the experiments supported by the computational results suggest that the self-assembly starts from the interactions between the small GO nanoparticles/sheets and TOCNF fibers, followed by the interactions with the long flexible TOCNF filaments and the flat, wide GO planes that wrap around each other forming amorphous porous biohybrid networks.

ACS Paragon Plus Environment

5

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 26

Figure 1. AFM and modeling images showing the self-assembly of TOCNF and GO. (a) The GO sheets (square area) and large GO nanoparticles (circle area) are embedded in the biohybrid network. (b) The thickness of the GO sheets and nanoparticles in the biohybrids. (c) The small GO nanoparticles are attached and accumulated along TOCNF from a single fiber scale, with highlights in the square area. The hydrogen bonds connecting the molecules to the TOCNF hydroxyls and carboxyl groups are highlighted (green lines). TOCNF is represented by a solvent-accessible surface, where the pink regions identify the TOCNF oxygens. Water and ions have been removed for clarity. (d) Flat adsorption with separate molecules. (e) Inter-connected sheets adsorbed almost perpendicularly on TOCNF extending toward a periodic image.

Dynamic aspect of the self-assembly in GO and TOCNF membranes. The self-assembled TOCNF and GO biohybrids were converted to a membrane via vacuum filtration and characterized with AFM experiments and molecular dynamic simulations. The surface and crosssection of the membrane are displayed in detail in Figure 2. From the AFM characterization shown in Figure 2a, it is apparent that, as expected, the GO sheets are located on the TOCNF+GO membrane surface. A few definite local arrangements and motifs can be identified. For example, in Figure 2a,b, it is visible that a triangular GO sheet is totally wrapped in TOCNF and thus embedded inside the membrane. Another GO sheet was found closer to the top of the membrane covering other TOCNF, leading to a blurred area shown in the square in Figure 2c, whereas GO nanoparticles surrounded and wrapped by TOCNF fibrils are also observed on the surface of the TOCNF+GO membrane (Figure 2d).

ACS Paragon Plus Environment

6

Page 7 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

The SEM analysis of the membrane cross-section (Figure 2e,f) reveals a layered structure with small micropores and hollow regions between the layers. The micropores were marked with red dashed lines in Figure 2e. The pure TOCNF membrane displays a layered structure as well, but in this case, pores or spaces between the layers are hardly discernible (Figure S2a,b). This outcome can be explained by the tight packing arrangement of TOCNF due to its homogenous nanoscale features. The addition of GO can disrupt the tight packing and create a loose network with increased porosity. In our recent studies on CNF/GO layered membranes, we have also observed that ‘standing inserted’ GO sheets at GO content as low as 1 wt% significantly increased the permeability of the membranes. 34

Figure 2. The surface morphology and cross-section of TOCNF+GO membrane. (a) AFM images show the GO sheets embedded inside and on top of the TOCNF+GO membrane. (b) The magnification area of the GO sheet embedded inside the TOCNF+GO membrane from the top square area. (c) The magnification area of the GO sheet on top of the TOCNF+GO membrane from a, bottom square area. (d) The GO nanoparticles randomly embedded or floated on the membrane. (e) The cross-section of the TOCNF+GO membrane, red dotted line indicating the micro pores. (f) The cross-section of the TOCNF+GO membrane, red square showing the layered structure. (g) (h) The dried structure where the GO sheets adsorbed on TOCNF crosslink the fibrils, self-interact

ACS Paragon Plus Environment

7

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 26

and determine pore size, shape and position inside the nanocellulose matrix. (g) x-y periodicity (replicas due to periodic boundary conditions in all dimensions): the pores are infinite canals running along the z-axis (black regions). (h) y-z periodicity: stripes of stacked GO structures delimit the canals, which could be filled with water and ions.

Interesting results, related to the self-assembling ability of the GO species on TOCNF, emerged from our MD trajectories and were thus discussed in relation to the corresponding experimental counterpart. Indeed, as already mentioned, the attractive interaction between TOCNF and GO was essentially due to the formation of a dense network of intermolecular hydrogen bonds involving mainly the hydroxyl groups of the two species. This network was sometimes reinforced by the involvement of carboxyl groups populating the edges of the sheets and different regions of functionalized TOCNF chains. The dynamic nature of the hydrogen bonds (i.e., bond breaking and formation) was responsible for the relatively low strength of the adhesion of the GO sheets in the sense that they could slight change their orientations on the TOCNF surface. The new positions of the sheets were essentially related to the morphology of the TOCNF facets and to the presence of other GO units on the fiber. Indeed, the self-assembly through inter-sheet hydrogen bonds (Figure 2h) was observed in different cases. Although selforganization could prevent a strong adherence of the sheets to the fiber surface and favor their partial extension in solution, it was noted that such arrangements could modulate the matrix porosity and generate efficient ion-trapping cages. Essentially, the surface roughness increased, and new capturing cavities were created. Starting the simulations from the separate initial arrangements of the GO sheets parallel to the TOCNF surface produced a dominant adsorption mode, where the sheets landed almost flat onto the support. However, among these final arrangements, inclined and partially adsorbed structures were also present. These resulted essentially from self-interactions. In particular, intermolecular interactions between molecules connected to different fibers (periodic images) generated partial adsorbents, inclined and attached to the TOCNF through the edge groups (Figure 2g). Such a scenario effectively illustrates how the nanofibrils could be crosslinked and provides a possible cause of the formation of relatively rich pore structures (Figure 2h). The action of the GO molecules can be further explained by focusing on the fibril separation at the beginning and at the end of the dynamics runs. In the simulated models, the initial separation between the fibrils in the y direction was approximately 22 Å, whereas along the z direction it was even smaller (16 Å). The presence of GO modulated the effective fiber distance, reducing it to less than 10 Å in the narrower region. As a matter of fact, the gap size is also

ACS Paragon Plus Environment

8

Page 9 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

correlated to the water content and can be reduced further when the composite (or membrane) is dried. Indeed, more direct evidence of the effects produced by the GO action on TOCNF packing was obtained when the membranes were consolidated by removing water from the simulation box (as done experimentally - drying process) and by equilibrating the configurations (at constant pressure and temperature) to their bulk geometries (Figure 2g,h). The consequences soon became apparent: the box volume was drastically reduced by approximately 65%, which corresponded to an approximately 30-40% increase in the volume of the isolated TOCNF (to be compared with the fibril packing). This new material was scattered with relatively large cavities, where the size, shape and distribution were controlled by the location and assembly of the GO units (Figure 2g,e). The maximum size of the pores observed in the simulated models was approximately 30 Å, estimated along the x-axis, and 15 Å along the z-axis (Figure S3). These dimensions could be further reduced by drying the system at high pressure and then re-equilibrating it at ambient pressure (Figure S4). In this case, the convergence toward an equilibrium structure was faster, and the resulting material was more compact, with a tighter molecular network and lower porosity due to the deformation of the adsorbate conformations. Metal ions capturing and clustering. With regard to the adsorption capability of the membrane toward metal-ions capture and clustering, the computational data analysis was carried out by comparing the results of different simulations in aqueous solution containing ions, GO+ions, TOCNF+ions, and TOCNF+GO+ions at the low (16 ions) and high (32 ions) concentrations. Even though at the beginning of the simulations, all ions were well-separated from one another, the sampling time was sufficiently long to disclose their binding to the supports and their limited self-aggregation. This was already apparent from the simulation of the ions in aqueous solution where small mixed aggregates made, at most, of five units were identified (Figure 3).

ACS Paragon Plus Environment

9

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 26

. Figure 3. Left: representative ionic clusters extracted from the simulation in aqueous solution. The copper ions (orange spheres) are surrounded by oxygen (red) and sulfur (yellow) atoms. Water molecules and all hydrogens are not displayed for clarity. Right: atomic Cu-X radial distribution functions. (a) 16 CuSO4 ions in water. (b) 6 GO molecules+16 CuSO4 ions. (c) TOCNF+16 CuSO4 ions. (d) TOCNF+6GO+16 CuSO4 ions. (e) TOCNF+6GO+32 CuSO4 ions. X: O(TOCNF) = all oxygens of the fibril; O(SO4 & WAT) and O(WAT & SO4) are water oxygens and the oxygens of the SO4 ions. This distinction has been introduced just because in ReaxFF, the connection table is recalculated at each iteration and, even though the SO4 oxygens were initially “labelled”, in principle, they could be replaced by water oxygens and vice-versa; O(GO) all oxygens of GO.

The overall tendency of copper ions at low concentration was to maintain long separations, be surrounded by water molecules and sulfate ions and be connected/adsorbed to/on the supports while preserving a variable amount of coordinated solvent oxygens (belonging to the water or

ACS Paragon Plus Environment

10

Page 11 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

sulfate coordination spheres). Only in a few clusters, the ions were accommodated at relatively short distances (4.5-5.2 Å), but these types of aggregations were rarely observed (approximately 2%). The copper coordination trend is clearly shown in Figure 3 (right hand side), where the atom-atom Cu-X radial distribution functions are displayed. The blue triangles represent the CuCu distances, whereas the red diamonds depict the coordination of copper to the sulfate counterions. All Cu-Cu plots show a similar behavior except the curve at high concentration, which corresponds to the TOCNF+GO+32 ions model. The presence of relatively sharp peaks at shorter distances suggest that due to the increased population of copper ions, the TOCNG+GO capturing sites become closer, and the ions are found at shorter distances entrapped in the matrix network (Figure S6-S7). A visual examination of the copper-oxygen coordination revealed that this could be realized through distorted tetragonal, octahedral, square pyramidal and trigonal bipyramidal geometries with average Cu−O separations of approximately 2.0 Å (Figure 3 – left hand side). The running coordination number of Cu surrounded by a first solvation shell of water molecules is approximately 5.6, in line with the literature, 35,36 whereas the second hydration layer is located at shorter distances, approximately 3.9-4.3 Å, than those in other studies (i.e., 4.55-4.72 Å)

36

because of the presence of the sulfate counterions. The latter, often coordinated with copper but also with each other, was found connected to the surface hydroxyl and carboxyl groups of the TOCNF through hydrogen bonding interactions. Another common trait shown by all Cu-O RDFs is that an exchange between the molecules of the first and second coordination shells is almost absent within the simulated time (and for the chosen sampling step). More evidence of clustering and metal ion immobilization on the biohybrids is provided by the analysis of the spatial distribution functions (SDFs) of the ions around the simulated models. These models are the three-dimensional maps of iso-surfaces that give the highest probability of finding the selected species at specific locations. Two representative SDFs of the copper and sulfur atoms for TOCNF+GO+ions at low and high concentrations are shown in Figure 4a,b with the density contours 1.5 times larger than the average solvent density. The structure displayed in the SDF plots is an average configuration calculated from the last portion of the trajectory. The maximum density regions are marked with solid orange (Cu) or yellow (S, as representative atom of the counterion) spheres.

ACS Paragon Plus Environment

11

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 26

Figure 4. On the left hand side are shown the SDFs of copper (orange) and sulfur (yellow) atoms for TOCNF+GO+ions at low (a) and high (b) concentration. The orange and yellow regions represent the density contours 2 times larger than the average solvent density. The maximum density is highlighted by spheres. The TOCNF (solvent-accessible surface where the oxygens are light pink) and GO molecules (cyan sticks with red oxygens) are average reconstructed conformations. On the righthand side the Atomic Cu-X radial distribution functions are displayed. (c) TOCNF+GO+16 ions. (d) TOCNF+GO+32 ions. X: O(TOCNF) are all oxygens of the fibril; O(GO) all oxygens of GO. This distinction is made because in ReaxFF the connection table is recalculated at each iteration and the oxygen initially connected to SO4 can become a water oxygen and vice-versa. S is the sulfur atom.

As far as the connection of Cu with the TOCNF and GO oxygens is concerned, the inspection of the RDFs displayed in Figure 4 suggests that Cu is well-anchored on the cellulose surface through the cooperation of the hydroxyl and carboxyl groups of both TOCNF and GO. However, a preference toward carboxyl oxygens is observed. This outcome is highlighted by the first peaks at 1.8-2.0 Å visible in all plots. Considering the coordination number, most of the carboxyl-Cu interactions are the result of TOCNF. At high ion concentrations, the adsorbed ions have the tendency to form clusters. In contrast, in the simulations of six GO molecules in solution, it was found that GO were slightly interconnected but not stacked on top of each other, and the propensity of the copper ions was to remain farther at approximately 4.5 Å from the GO

ACS Paragon Plus Environment

12

Page 13 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

capturing groups but still bind through the hydrogen bonds between water molecules of the first coordination shell and the hydroxyl/carboxyl groups of GO (Figure 4d). As shown in Figure 4a,b, the maximum density zones of Cu (orange contours) are all located on top of the TOCNF surface, and sometimes, the Cu atoms are found entrapped between the TOCNF interface and the GO sheets (Figure 4a). Most of the lower density regions are distributed near the TOCNF surface, as well, and only a few Cu ions remain in solution. This outcome clearly indicates that not all of the copper ions were adsorbed onto the supports within the simulated time, and some remained in solution coordinated to their corresponding counterions (Figure 4a,b). It is interesting to observe that the copper-ion distributions are tighter than the sulfur ones because their tendency is to accommodate stably on the support. The stability depends on the type of oxygens involved in the interaction and on the steric capture mechanism induced by the compressing action of the GO sheets. The most stable and persistent attachments comprise carboxyl oxygens. In contrast, sulfur distributions are sparser and spread, suggesting that the SO4 moieties tend to follow the motion of water and hydrated copper ions. The latter ones are more mobile when adsorbed as CuSO4 clusters, especially if their orientation is toward the solvent (Figure S6b), and they often travel from one site to another one in search of a steady state. The alignment of copper densities crossing the fiber, visible in Figure 4, implies the formation of a chain of adsorbed Cu ions (as shown in Figures S5 and S6), which are not strongly connected to the interface because the interaction does not seem to involve carboxyl oxygens. Essentially, their motion is regulated by the reorientation of the hydroxyl moieties of TOCNF, which keeps them connected to the surface. At high copper concentration, the formation of mixed (containing SO4 ions) clusters on the TOCNF surface was observed in dry conditions (Figure 5) Different views of the final configurations of TOCNF+GO+ions were obtained after drying the solvated models extracted from the dynamics. The drying process was performed by removing all water molecules except those belonging to the first coordination sphere of the copper ions (top) or all water molecules (bottom). In each picture (top and bottom), three sides of the final box, namely, XY (top-left), ZY (top-right) and XZ (bottom), are shown as representing the molecular/ionic arrangements. Apparently, a more packed structure with small pores is obtained when a few water molecules remain entrapped between the fibrils, and GO can adopt more relaxed stacked conformations laid on the fibril surface, whereas in the other case, the formation of larger pores seems induced by a different arrangement adopted by the sheets, which are more bent and inclined in relation to the fibers.

ACS Paragon Plus Environment

13

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 26

Figure 5. Final configurations, dried and packed. In the top image all water molecules were removed except those belonging to the first coordination sphere of the copper ions (see Figure S5 for a better view of the small pores). In the bottom image all water molecules were removed. In each picture the XY (top-left), ZY (top-right) and XZ (bottom) points of view are shown. Apparently, a more packed structure with small pores is obtained when a few water molecules remain entrapped between the fibrils and GO can adopt more relaxed stacked conformations laid on the fibril surface, whereas in the other case the formation of larger pores seems induced by a

ACS Paragon Plus Environment

14

Page 15 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

different arrangement adopted by the sheets, which are more bent and inclined in relation to the fibers. Periodic images are displayed to render more clearly coverage.

Figure 6. Surface morphology and modeling structure of TOCNF+GO membrane after adsorption of Cu(II) under different magnification. (a, b, c) The height images. (d, e, f) The corresponding

ACS Paragon Plus Environment

15

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 26

peak force images. The red arrows display the magnifications of the square areas. The circles note the nanolayer and nanoparticles on the surface of the membrane. (g) Spatial distribution functions of copper (orange) and sulfur (yellow) atoms surrounding the TOCNF+GO model at high concentration. The orange and yellow regions represent density contours 2 times larger than the average solvent density. TOCNF (solvent accessible surface where the oxygens are light pink) and GO molecules (cyan sticks with red oxygens) are average conformations (reconstructed)

The experimental evidence confirms this scenario, suggesting that the stationary phase is characterized by a TOCNF+GO membrane covered with a layer of nanoparticles and a few nanoclusters (Figure 6a,b,c). In particular, in Figure 6d, the circled area comprises a square shape particle that resembles a GO sheet combined with copper nanoclusters supported on the membrane. In this case, a detailed characterization is impossible due to the size of GO that partially envelops the metal ions. Further evidence of this phenomenon is provided by the AFM video images, which were used for monitoring the sample and probe movement (Figure S2d,e). They exhibit bright iridescence from the TOCNF+GO membrane without Cu(II) adsorption, which is mainly due to the reflection of GO on the membrane surface under illumination during AFM measurements (Figure S2c). After Cu(II) adsorption, the iridescence becomes much weaker, which is due to the copper nanolayer and nanoclusters covered on the surface of the TOCNF+GO membrane, shading the GO reflection under illumination (Figure S2d). Moreover, the SEM-EDS results consolidate the conclusion that Cu(II) ions were successfully adsorbed by the biohybrid membrane (Figure 6g). CONCLUSIONS Advanced microscopy techniques revealed the self-assembly behavior and the porous structure in TOCNF and GO networks, as well as their biohybrid membrane. The MD simulations significantly clarified the nature of aggregation and pore formation inside the biohybrid membrane, which seemed related to the degree of GO basal plane bending, to the action of the solvent, to the intermolecular H-bonds between the GO sheets and, finally, to their location in the fiber skeleton. It was further demonstrated that copper ions are captured inside the biohybrid network through cooperative coordination of the hydroxyl and carboxyl groups of both the TOCNF and the GO components. Preferential interactions mainly involve the carboxyl oxygens and determine stronger adsorptions. It was observed that the tendency to stably accommodate on the support also depends on the steric capture mechanism due to the compressing action of the GO sheets on the fiber surface. The formation of copper clusters on the fibers is also reproduced, and their fiber crosslinking action is shown. These outcomes perfectly agree with the ACS Paragon Plus Environment

16

Page 17 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

experimental findings. This combined experimental-theoretical investigation successfully clarified, at the nanoscale, the structural and dynamical aspects of the biohybrids that were not yet well-characterized by the studies present in the literature, thus providing the guidelines to understanding the behavior of cellulosic- and GO-based biohybrids for potential applications as ultralight aerogels,

37

energy and memory storage materials,

38

hydrogels, tissue engineering

scaffolds39 and robust supercapacitors.40 MATERIALS AND METHODS Experimental Section: Materials. The TEMPO-oxidized cellulose nanofibers (TOCNF) (1.0 mmol/g, 1 wt%) were prepared and provided by Liu et al41 Stockholm University, Sweden. Graphene oxide (4 mg/mL), copper(II) sulfate pentahydrate (99.999%), and (3-Aminopropyl) triethoxysilane (APTES) were all purchased from Sigma-Aldrich, Germany and used as received. Distilled water was used as the dispersion medium. ScanAsyst Air probe was purchased from Bruker, Santa Barbara, USA and was applied for AFM characterization. Preparation of TOCNF and GO-Based Biohybrid Membranes, and Cu(II) Adsorption. Suspension of 1 g/L TOCNF+GO (100/1, wt %) in 200 mL was prepared by mixing 20 g TOCNF and 0.5 mL GO suspension and dilution with distilled water. Five mL TOCNF+GO suspension was vacuum-filtered after sonication for 3 minutes, using a 0.45-um pore-size filter membrane (DVPP, Millipore) under filtration pressure of 1 bar. Two TOCNF+GO membranes were fabricated, and one of the membranes was immersed in Cu(II) solution for 30 minutes and then filtered and rinsed by distilled water. The pure TOCNF membrane was prepared by the same method as the control sample. The TOCNF+GO membranes with and without adsorption of Cu(II) were air-dried and investigated by the microscopy and spectroscopy techniques both on the surface and cross-sections. Characterizations. Atomic Force Microscopy (AFM). The topography of the samples with and without copper adsorption was performed with a Dimension Icon AFM (Bruker, Nanoscope controller, Santa Barbara, California, USA). The TOCNF and GO dispersed samples were prepared as in our earlier study. 19 The TOCNF, GO and TOCNF+GO (100:1, wt %) suspension was drop-coated on the APTES-modified mica surface on the AFM metal disc and was used to characterize the structure of the nanoparticles and the biohybrids. The TOCNF+GO (100:1 wt%) membranes with and without Cu(II) were directly attached to the AFM metal disc using AFM double side tape. The probes (Model: ScanAsyst-air, Bruker) were treated with the UV Ozone Cleaner (ProCleaner™ Plus, US) for 20 minutes before scanning. The height, peak-force error

ACS Paragon Plus Environment

17

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 26

images were recorded under ScanAsyst-air mode. The diameter of the TOCNF fiber and the thickness of GO were post processed with NanoScope Analysis 1.5 software (Bruker), using section function by measuring the data perpendicular and laterally towards the fiber length direction and GO, respectively. The measured values were directly marked in Figure S1.

Scanning Electron Microscopy - Energy Dispersive X-ray Spectroscopy (SEM-EDS). The morphologies of the pure TOCNF cross-section

19

and TOCNF+GO (100:1 wt%) membranes

were observed using scanning electron microscopy on a JEOL JSM-7401F microscope (Japan). The membranes were cut and anchored at 0 and 90 degrees to the holder to check the surface and cross-sections. The samples were sputter-coated with gold for 20 seconds before measurements to avoid the burning/damaging to the samples during measurements. Images were taken at 2 kV and approximately 10-mm working distance. The EDS spectrums were obtained at 15 kV and 10-mm working distance with the accumulation time of 2 minutes, which was used to confirm the elements in the sample.

Computational Section: Preliminary Model Building. To characterize in detail the ability of cellulose nanofibrils (CNF) and GO to capture copper ions in solution and their expected cooperation, we built a series of representative supramolecular models inspired by the data available in the literature or kindly provided by the authors. 23,42,43 The infinite length of CNF was rendered through a periodic supercell, with chains covalently bonded across the x border, whereas the CNF packing was controlled by adjusting the dimensions of the simulation box in the y and z directions. The crystal structure of Iβ cellulose determined by X-ray and neutron diffraction

44

was employed to create the sixteen CNF chains made of sixteen glucosyl residues

each, which were aligned along the x direction of the supercell. According to the literature42,43,44, our previous experience on fibril structures45,46and present experimental proportions, this size is a good compromise between the computational cost and a meaningful comparison with the experimental data. A reasonable configuration of GO (28 carbon rings), in line with the CNF model, was downloaded from the Automated Topology Builder (ATB) and Repository site47 and modified by replacing eight of the edge hydrogen atoms with carboxyl groups (Figure S1a). The new molecule was quantum-mechanically energy optimized at the DFT/M062X/6-31G(d) level and the resulting structure was used as a starting conformation in all simulations. The choice of the FF was based on suggestions found in the literature48 where a comparison among three different reactive FFs is reported. Considering all the aspects commented there, we decided that for the type of simulations we were to perform the best choice was the force field

ACS Paragon Plus Environment

18

Page 19 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

developed by Chenoweth et al. 49 and made available in the ADF package under the name of HCO.ff. SO42- parameters were taken from Ref.

50

where the authors carried out a MC

optimization of the MgSO4 hydrates (the force was released with the name HOSMg.ff - ADF). The copper parameters were extracted from CuCl-H2O.ff: (Cu/Cl/H/O) 29 and a direct bond with sulfur atoms was not parametrized based on the hypothesis that a direct interaction was with the oxygens surrounding of the SO42- group (essentially, SO42- ions maintained their structure). Before building the composite systems, the behavior of the CNF model in solution was tested through a series of short (hundreds of picoseconds) all-atom molecular dynamics (MD) simulations in the NPT ensemble, at 300 K and 1 atm. The volume of the system was allowed to change along the y and z directions that are perpendicular to the fibril axis (x direction). The analysis of the MD trajectories confirmed the solution stability of the overall structure, which preserved both its morphology and crystallinity, although the exposed fibril portions (the twelve chains surrounding the core – made of four chains) reoriented their “side groups” (OH and CH2OH moieties), not involved in inter-chain hydrogen bonds, toward the solvent. Then, the equilibrated configuration was subsequently modified to obtain a model functionalized with carboxyl groups (TOCNF). The procedure reflected the experimental TEMPO (2,2,6,6tetramethylpiperidine-1-oxylradical)-mediated oxidation where the exposed C6 primary hydroxyls of the surface (twelve chains in our case) were selectively replaced with C6 carboxyls. 4,51

In the final model, approximately 45% of the primary hydroxyls were changed to the COOH

moieties. The new configuration was energy-minimized and equilibrated in aqueous solution through the NPT MD runs (cell length, Lx, fixed) at ambient temperature and pressure. The final arrangement of the chains formed a slightly twisted fiber, inflated around the edges (due to the interaction with the solvent). This model was used as the starting geometry in all multicomponent simulations. The Complex Models and Simulation Set-Up. All MD simulations were carried out by means of the Amsterdam Density Functional (ADF)/ReaxFF

52

and LAMMPs codes53 resorting also to

the HPC resources offered by CINECA through the Italian Super Computing Resource Allocation (ISCRA) project IscrC_GOHNC. The TOCNF model was surrounded by six GO sheets, which were placed in a proximity of the fibril surface, and the complex was inserted in a box of water containing approximately 4500 molecules and 16 CuSO4 ions, randomly dispersed in the solvent (simulations at high ion concentration - 32 ions - were also carried out). The final size of the simulation box, after

ACS Paragon Plus Environment

19

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 26

equilibration, was approximately 8.3 x 5.1 x 3.9 nm3. The multicomponent systems consisted of TOCNF, GO sheets, water molecules, Cu(II) ions and SO4 counterions (Figure 7).

Figure 7. Multicomponent system consisting of TOCNF (solid surface – light red oxygens and gray carbon and hydrogens), six GO sheets (cyan sticks), Cu(II) ions (orange spheres) and SO4 ions (yellow and red balls). Water molecules have not been displayed for clarity. Periodic images are shown together with the simulation box (white parallelepiped).

The most crowded combination contained these species together, whereas the simplest cases were made of just one species and water. To disclose the synergistic effects on the ion trapping process due to the TOCNF+GO association, the capturing agents (TOCNF or GO) and ions were simulated independently, namely, TOCNF +ions and GO+ions, and together but without the ions (TOCNF+GO). The latter simulations were used to prepare realistic models of the TOCNF+GO complexes for further tests. Several other configurations were created and simulated in different environmental conditions to characterize tentative fibril packing, drying, porosity, etc. However, to be concise, we focus only on the three main models mentioned above (TOCNF +GO+ions, TOCNF +ions and GO+ions) and refer to other studies only to explain specific effects and behaviors, observed experimentally. It is worth mentioning that the reactivity was always on, and no restraints were introduced in the systems. All species could interact with each other and modify their conformation in response to the environment.

ACS Paragon Plus Environment

20

Page 21 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

MDs were carried out first in the NVT ensemble at T=300 K for approximately 200 ps to equilibrate the whole system at the chosen temperature and allow the GO sheets to touch the TOCNF surface. Then, the equilibration was extended for another 500 ps in the NPT ensemble to reach the correct packing of all system components and the appropriate density. Afterwards, the production dynamics was performed in the NVE ensemble for approximately 2 ns, and the system structures were collected every 0.05 ps. The temperature and pressure were controlled through the Berendsen’s thermostat and barostat54 with the relaxation constants of 0.1 ps. The time-step was set to 0.25 fs and the dynamics of ten different initial arrangements was explored in each case. The final portions of these trajectories were combined and analyzed together. We are aware that in such a timescale we cannot provide the whole picture of the behavior of these types of materials as a whole (that would be impossible also extending system size, simulation time and number of starting configurations), but we believe that using a few simplified models, such the ones we have used here, representing specific portions of the complex systems, we could give an idea of possible local features and characteristic effects that can be correlated to the experimental findings. The analysis was focused on descriptors that could be useful to disclose ions clustering and capture, namely, the atom-atom radial distribution functions (RDFs), spatial distribution functions (SDFs), various types of atom-atom distances and hydrogen bond networks. To estimate the influence of the TOCNF+GO complexation on matrix packing and porosity, selected structures were extracted from the trajectories, dried, by removing water and resimulated in the NPT ensemble at ambient temperature until complete equilibration.

ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge at the ACS Publications website at http://pubs.acs.org. The height images and measured size of TOCNF, and the GO sheets/nanoparticles; the SEM image of the cross-section of pure TOCNF membrane, the surface of TOCNF+GO membrane and the video images of TOCNF+GO membrane before and after adsorption under illumination during AFM measurements; different views of the pore structure filled with GO, ions and water; drying and packing at high pressure until equilibrium and then equilibration at ambient pressure

ACS Paragon Plus Environment

21

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 26

of a configuration extracted from the production run of the TOCNF+6GO+16 ions simulation; a comparison of two different drying conditions; final configuration of TOCNF+6GO+32 ions sampled after 2 ns; several typical arrangements of the copper ions adsorbed on TOCNF.

The authors declare no competing financial interest.

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]

ACKNOWLEGEMENTS The authors gratefully acknowledge the financial support of the Swedish research council (VR, grant No: 621-2013-5997 and 2017-04254). We thank Liu, Y. Stockholm University, Sweden, who kindly provided the TOCNF.

REFERENCES (1) Klemm, D.; Kramer, F.; Moritz, S.; Lindström, T.; Ankerfors, M.; Gray, D.; Dorris, A. Nanocelluloses: A New Family of Nature-Based Materials. Angew. Chem. Int. Ed. 2011, 50, 5438-5466. (2) Huang, H.; Ying, Y.; Peng, X. Graphene Oxide Nanosheet: An Emerging Star Material for Novel Separation Membranes. J. Mater. Chem. A 2014, 2, 13772-13782. (3) Dervin, S.; Dionysiou, D. D.; Pillai, S. C. 2D Nanostructures for Water Purification: Graphene and Beyond. Nanoscale 2016, 8, 15115-15131. (4) Isogai, A.; Saito, T.; Fukuzumi, H. TEMPO-Oxidized Cellulose Nanofibers. Nanoscale 2011, 3, 7185. (5) Saito, T.; Kimura, S.; Nishiyama, Y.; Isogai, A. Cellulose Nanofibers Prepared by TEMPO-Mediated Oxidation of Native Cellulose. Biomacromolecules 2007, 8, 2485-2491. (6) Suk, J. W.; Piner, R. D.; An, J.; Ruoff, R. S. Mechanical Properties of Monolayer Graphene Oxide. ACS Nano 2010, 4, 6557-6564. (7) Sun, P.; Zhu, M.; Wang, K.; Zhong, M.; Wei, J.; Wu, D.; Xu, Z.; Zhu, H. Selective Ion Penetration of Graphene Oxide Membranes. ACS Nano 2013, 7, 428-437. (8) Yao, Q.; Fan, B.; Xiong, Y.; Jin, C.; Sun, Q.; Sheng, C. 3D Assembly Based on 2D Structure of Cellulose Nanofibril/Graphene Oxide Hybrid Aerogel for Adsorptive Removal of Antibiotics in Water. Sci. Rep. 2017, 7, 45914. (9) Fang, Q.; Zhou, X.; Deng, W.; Zheng, Z.; Liu, Z. Freestanding Bacterial Cellulose-Graphene Oxide Composite Membranes with High Mechanical Strength for Selective Ion Permeation. Sci. Rep. 2016, 6, 33185.

ACS Paragon Plus Environment

22

Page 23 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

(10) Chen, C.; Zhang, T.; Zhang, Q.; Chen, X.; Zhu, C.; Xu, Y.; Yang, J.; Liu, J.; Sun, D. Biointerface by Cell Growth on Graphene Oxide Doped Bacterial Cellulose/Poly(3,4-ethylenedioxythiophene) Nanofibers. ACS Appl. Mater. Interfaces 2016, 8, 10183-10192. (11) Sitko, R.; Musielak, M.; Zawisza, B.; Talik, E.; Gagor, A. Graphene Oxide/Cellulose Membranes in Adsorption of Divalent Metal Ions. RSC Adv. 2016, 6, 96595-96605. (12) Zhu, C.; Soldatov, A.; Mathew, A. P. Advanced Microscopy and Spectroscopy Reveal the Adsorption and Clustering of Cu(II) onto TEMPO-Oxidized Cellulose Nanofibers. Nanoscale 2017, 9, 7419-7428. (13) Melone, L.; Rossi, B.; Pastori, N.; Panzeri, W.; Mele, A.; Punta, C. TEMPO-Oxidized Cellulose Cross-Linked with Branched Polyethyleneimine: Nanostructured Adsorbent Sponges for Water Remediation. ChemPlusChem 2015, 80, 1408-1415. (14) Karim, Z.; Hakalahti, M.; Tammelin, T.; Mathew, A. P. In Situ TEMPO Surface Functionalization of Nanocellulose Membranes for Enhanced Adsorption of Metal Ions from Aqueous Medium. RSC Adv. 2017, 7, 5232-5241. (15) Hontoria-Lucas, C.; López-Peinado, A. J.; López-González, J. d. D.; Rojas-Cervantes, M. L.; MartínAranda, R. M. Study of Oxygen-Containing Groups in a Series of Graphite Oxides: Physical and Chemical Characterization. Carbon 1995, 33, 1585-1592. (16) He, H.; Klinowski, J.; Forster, M.; Lerf, A. A New Structural Model for Graphite Oxide. Chem. Phys. Lett. 1998, 287, 53-56. (17) Gao, W.; Majumder, M.; Alemany, L. B.; Narayanan, T. N.; Ibarra, M. A.; Pradhan, B. K.; Ajayan, P. M. Engineered Graphite Oxide Materials for Application in Water Purification. ACS Appl. Mater. Interfaces 2011, 3, 1821-1826. (18) Wang, B.; Lou, W.; Wang, X.; Hao, J. Relationship between Dispersion State and Reinforcement Effect of Graphene Oxide in Microcrystalline Cellulose-Graphene Oxide Composite Films. J. Mater. Chem. 2012, 22, 12859-12866. (19) Zhu, C.; Liu, P.; Mathew, A. P. Self-Assembled TEMPO Cellulose Nanofibers: Graphene OxideBased Biohybrids for Water Purification. ACS Appl. Mater. Interfaces 2017, 9, 21048-21058. (20) Liu, P.; Oksman, K.; Mathew, A. P. Surface Adsorption and Self-Assembly of Cu(II) Ions on TEMPO-Oxidized Cellulose Nanofibers in Aqueous Media. J. Colloid Interface Sci. 2016, 464, 175-182. (21) Zhu, C.; Dobryden, I.; Rydén, J.; Öberg, S.; Mathew, A.P.; Holmgren, A. Adsorption Behavior of Cellulose and Its Derivatives toward Ag(I) in Aqueous Medium: An AFM, Spectroscopic, and DFT Study. Langmuir 2015, 31, 12390-12400. (22) Paajanen, A.; Vaari, J. High-Temperature Decomposition of the Cellulose Molecule: a Stochastic Molecular Dynamics Study. Cellulose 2017, 24, 2713-2725. (23) Li, Y.; Zhu, H.; Zhu, S.; Wan, J.; Liu, Z.; Vaaland, O.; Lacey, S.; Fang, Z.; Dai, H.; Li, T.; Hu, L. Hybridizing Wood Cellulose and Graphene Oxide toward High-Performance Fibers. NPG Asia Mater. 2015, 7, e150. (24) Medhekar, N. V.; Ramasubramaniam, A.; Ruoff, R. S.; Shenoy, V. B. Hydrogen Bond Networks in Graphene Oxide Composite Paper: Structure and Mechanical Properties. ACS Nano 2010, 4, 2300-2306. (25) Achtyl, J. L.; Unocic, R. R.; Xu, L.; Cai, Y.; Raju, M.; Zhang, W.; Sacci, R. L.; Vlassiouk, I. V.; Fulvio, P. F.; Ganesh, P.; Wesolowski, D. J.; Dai, S.; Van Duin, A. C. T.; Neurock, M.; Geiger, F. M. Aqueous Proton Transfer Across Single-Layer Graphene. Nat. Commun. 2015, 6, 6539. (26) Singh, S. K.; Srinivasan, S. G.; Neek-Amal, M.; Costamagna, S.; Van Duin, A. C. T.; Peeters, F. M. Thermal Properties of Fluorinated Graphene. Phys. Rev. B Condens. Matter Mater. Phys. 2013, 87, 104114. (27) Neyts, E. C.; Van Duin, A. C. T.; Bogaerts, A. Formation of Single Layer Graphene on Nickel under Far-from-Equilibrium High Flux Conditions. Nanoscale 2013, 5, 7250-7255.

ACS Paragon Plus Environment

23

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 26

(28) Huang, L.; Seredych, M.; Bandosz, T. J.; Van Duin, A. C. T.; Lu, X.; Gubbins, K. E. Controllable Atomistic Graphene Oxide Model and Its Application in Hydrogen Sulfide Removal. J. Chem. Phys. 2013, 139, 194707. (29) Rahaman, O.; Van Duin, A. C. T.; Bryantsev, V. S.; Mueller, J. E.; Solares, S. D.; Goddard III, W. A.; Doren, D. J. Development of a ReaxFF Reactive Force Field for Aqueous Chloride and Copper Chloride. J Phys Chem A 2010, 114, 3556-3568. (30) Van Duin, A. C. T.; Bryantsev, V. S.; Diallo, M. S.; Goddard, W. A.; Rahaman, O.; Doren, D. J.; Raymand, D.; Hermansson, K. Development and Validation of a ReaxFF Reactive Force Field for Cu Cation/Water Interactions and Copper Metal/Metal Oxide/Metal Hydroxide Condensed Phases. J Phys Chem A 2010, 114, 9507-9514. (31) Senftle, T. P.; Hong, S.; Islam, M. M.; Kylasa, S. B.; Zheng, Y.; Shin, Y. K.; Junkermeier, C.; EngelHerbert, R.; Janik, M. J.; Aktulga, H. M.; Verstraelen, T.; Grama, A.; Van Duin, A. C. T. The ReaxFF Reactive Force-Field: Development, Applications and Future Directions. npj Computational Mater. 2016, 2, 15011. (32) Pan, S.; Aksay, I. A. Factors Controlling the Size of Graphene Oxide Sheets Produced via the Graphite Oxide Route. ACS Nano 2011, 5, 4073-4083. (33) Stankovich, S.; Piner, R. D.; Nguyen, S. T.; Ruoff, R. S. Synthesis and Exfoliation of IsocyanateTreated Graphene Oxide Nanoplatelets. Carbon 2006, 44, 3342-3347. (34) Liu,P. Zhu, C. Mathew, P. A. Graphene Oxide Coated Nanocellulosic Fibrous Membranes for Ultrafast Separation of Dyes from Water. J Hazar Mater 2018. Submitted. (35) Pasquarello, A.; Petri, I.; Salmon, P. S.; Parisel, O.; Car, R.; Tóth, É; Powell, D. H.; Fischer, H. E.; Helm, L.; Merbach, A. E. First Solvation Shell of the Cu(II) aquaion: Evidence for Fivefold Coordination. Science 2001, 291, 856-859. (36) Schwenk, C. F.; Rode, B. M. New Insights into the Jahn-Teller Effect through ab Initio QuantumMechanical/Molecular-Mechanical Molecular Dynamics Simulations of Cu(II) in water. ChemPhysChem 2003, 4, 931-943. (37) Li, C.; Wu, Z.; Liang, H.; Chen, J.; Yu, S. Ultralight Multifunctional Carbon-Based Aerogels by Combining Graphene Oxide and Bacterial Cellulose. Small 2017, 13, 1700453. (38) Kafy, A.; Sadasivuni, K. K.; Kim, H.; Akther, A.; Kim, J. Designing Flexible Energy and Memory Storage Materials Using Cellulose Modified Graphene Oxide Nanocomposites. Phys. Chem. Chem. Phys. 2015, 17, 5923-5931. (39) Si, H.; Luo, H.; Xiong, G.; Yang, Z.; Raman, S. R.; Guo, R.; Wan, Y. One-step In Situ Biosynthesis of Graphene Oxide-Bacterial Cellulose Nanocomposite Hydrogels. Macromol. Rapid Commun. 2014, 35, 1706-1711. (40) Liu, Y.; Zhou, J.; Zhu, E.; Tang, J.; Liu, X.; Tang, W. Facile Synthesis of Bacterial Cellulose Fibres Covalently Intercalated with Graphene Oxide by One-Step Cross-Linking for Robust Supercapacitors. J. Mater. Chem. C 2015, 3, 1011-1017. (41) Liu, Y.; Gordeyeva, K.; Bergström, L. Steady-Shear and Viscoelastic Properties of Cellulose Nanofibril–Nanoclay Dispersions. Cellulose 2017, 24, 1815-1824. (42) Charlier, L.; Mazeau, K. Molecular Modeling of the Structural and Dynamical Properties of Secondary Plant Cell Walls: Influence of Lignin Chemistry. J Phys Chem B 2012, 116, 4163-4174. (43) Alqus, R.; Eichhorn, S. J.; Bryce, R. A. Molecular Dynamics of Cellulose Amphiphilicity at the Graphene-Water Interface. Biomacromolecules 2015, 16, 1771-1783. (44) Nishiyama, Y.; Langan, P.; Chanzy, H. Crystal Structure and Hydrogen-Bonding System in Cellulose Iß from Synchrotron X-ray and Neutron Fiber Diffraction. J. Am. Chem. Soc. 2002, 124, 90749082. (45) Monti, S.; Bronco, S.; Cappelli, C. Toward the Supramolecular Structure of Collagen: A Molecular

ACS Paragon Plus Environment

24

Page 25 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Dynamics Approach. J Phys Chem B 2005, 109, 11389-11398. (46) Monti, S.; Bramanti, E.; Porta, V. D.; Onor, M.; D'Ulivo, A.; Barone, V. Interaction of Collagen with Chlorosulphonated Paraffin Tanning Agents: Fourier Transform Infrared Spectroscopic Analysis and Molecular Dynamics Simulations. Phys. Chem. Chem. Phys. 2013, 15, 14736-14747. (47) http://atb.uq.edu.au/molecule.py?molid=34774 (accessed on 9 June 2018). (48) Dri, F. L.; Wu, X.; Moon, R.J.; Martini, A.; Zavattieri, P.D. Evaluation of Reactive Force Fields for Prediction of the Thermo-Mechanical Properties of Cellulose iß. Comput Mater Sci. 2015, 109, 330-340. (49) Chenoweth, K.; van Duin, A. C. T.; Goddard, W. A. ReaxFF Reactive Force Field for Molecular Dynamics Simulations of Hydrocarbon Oxidation. J Phys Chem A. 2008, 112, 1040-1053. (50) Zhang, H.; Iype, E.; Nedea, S. V.; Rindt, C. C. M. Molecular Dynamics Study on Thermal Dehydration Process of Epsomite (MgSO4·7H2O). Mol Simul. 2014, 40, 1157-1166. (51) Okita, Y.; Saito, T.; Isogai, A. Entire Surface Oxidation of Various Cellulose Microfibrils by TEMPO-Mediated Oxidation. Biomacromolecules 2010, 11, 1696-1700. (52) Baerends, E.J; Ziegler, T.; Atkins, A. J.; Autschbach, J.; Baseggio, O.; Bashford, D.; Bérces, A.; Bickelhaupt, F.M.; Bo, C.; Boerrigter, P.M.; Cavallo, L.; Daul, C.; Chong, D.P.; Chulhai, D.V.; Deng, L.; Dickson, R.M.; Dieterich, J.M.; Ellis, D.E.; van Faassen, M.; Fan, L.; et al. ADF2017, SCM, Theoretical Chemistry, Vrije Universiteit, Amsterdam, The Netherlands, http://www.scm.com. (53) Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1-19. (54) Berendsen, H. J. C.; Postma, J. P. M.; Van Gunsteren, W. F.; Dinola, A.; Haak, J. R. Molecular Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81, 3684-3690.

ACS Paragon Plus Environment

25

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 26

■ TABLE OF CONTENTS

ACS Paragon Plus Environment

26