Synergistic Amination of Graphene: Molecular Dynamics and

Oct 21, 2015 - In this letter, we employ ab initio gas-phase thermochemistry and long-scale finite-temperature electronic-structure-based molecular dy...
2 downloads 12 Views 1000KB Size
Subscriber access provided by UNIV OF NEBRASKA - LINCOLN

Letter

Synergistic Amination of Graphene: Molecular Dynamics and Thermodynamics Vitaly V. Chaban, and Oleg V. Prezhdo J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.5b02206 • Publication Date (Web): 21 Oct 2015 Downloaded from http://pubs.acs.org on October 22, 2015

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry Letters 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 18

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

The Journal of Physical Chemistry Letters

Synergistic Amination of Graphene: Molecular Dynamics and Thermodynamics Vitaly V. Chaban1 and Oleg V. Prezhdo2 1

Instituto de Ciência e Tecnologia, Universidade Federal de São Paulo, 12231-280, São José dos

Campos, SP, Brazil 2

Department of Chemistry, University of Southern California, Los Angeles, CA 90089, United

States

Abstract. Functionalization of graphene using organic moieties constitutes an affordable way to modulate its physical and chemical properties. Finding an exact structural formula of functionalized graphene using experimental approaches is challenging. We studied in detail thermal stability and thermodynamics of amino- and ethylamino-graphene, and found a surprising synergistic effect: More amino-groups stabilize functionalized graphene favoring further amination, whereas a small concentration of amino-groups is unstable in many cases. The functional groups can be attached either on the same or simultaneously on different sides of the graphene sheet. Deformation of functionalized graphene is proportional to the number of amino groups. Complete amination leading to formation of the ultimate product, Cx(NH2)x, is hindered sterically. Our study assists in determination of structure of chemically modified graphene, and makes specific predictions that can be tested and validated experimentally.

Key words: graphene; amine; alkylamine; functionalization; simulation; thermodynamics.

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

TOC Graphic

ACS Paragon Plus Environment

Page 2 of 18

Page 3 of 18

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

The Journal of Physical Chemistry Letters

Graphene (GRA) exhibits many unusual and valuable properties, stimulating active investigation by chemists and physicists.1-33 GRA functionalization is of crucial importance for efficient solubilization and for tuning numerous chemical and physical characteristics.1, 9, 11-20, 2231, 33-37

Certain physical properties of GRA can be improved further if the grafted atoms or

groups weaken the impact of existing structural defects. Functionalization facilitates binding of GRA and organic molecules.11, 12, 18, 22, 27, 33 It also makes GRA an ideal platform for anchoring metallic nanoparticles.13, 17 Functionalized graphene (FG) can be obtained by both covalent and non-covalent modification. The first class of techniques usually delivers stronger changes in GRA properties. For instance, electrical conductivity of FG is much smaller than that of pristine GRA; surface area decreases due to oxidation, chemical reduction and sonication; chemical reactivity of graphene oxide (GO) appears much different than that of pure GRA or graphite.1 FG can be further processed by solvent-assisted techniques, such as drop-casting, layer-by-layer assembly, spraying, filtration and spin-coating. Sheets of FG do not agglomerate in solution, thus maintaining the properties of graphene as an isolated two-dimensional layer of carbon with a honeycomb structure.1, 2 The most widely-spread functionalization of GRA involves some kind of oxidation with a formation of multiple carbonyl, carboxyl, hydroxyl, epoxide and diol groups at the surface.3, 18, 30 Oxidation is actively employed for exfoliation directly from graphite. Alternative recipes to obtain FG rely on nucleophilic substitution, electrophilic addition, condensation, and other addition chemical reactions.1 Nucleophilic substitution mainly takes place via the epoxy groups of GO. The amine group, -NH2, bearing a lone pair of electrons attacks the epoxy group readily at room conditions in the aqueous medium. Aliphatic and aromatic amines, amino acids, roomtemperature ionic liquids and amine terminated biomolecules were proven to efficiently functionalize GRA and GO.1 This method exhibits a potential to be used for industrial scale applications. Spontaneous grafting of the aryl diazonium salt to GRA is a straightforward example of an electrophilic substitution.14 Condensation of GRA occurs with isocyanate,

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 18

diisocyanate and amine compounds via the formation of amides and carbamate ester linkages.26 Two or more molecules combine to form a larger molecule in organic addition reactions. An important illustration of the addition reactions involving GRA is the 1,3-dipolar cycloaddition of azomethine.12 Different substituted amines were embedded into GO for energy storage applications, electrocatalytic oxygen reduction reactions, and electrochemical sensors.29,

38

An improved

performance was observed in certain cases. Similarly, the amine functionalization of carbon nanotubes was used to enhance the controlled covalent bonding to polymers and biological molecules.27 These wide applications inspire research efforts towards synthesis of amino-FG with a fine control over its structure and properties. Detailed knowledge about amino-FG is required for this purpose. Navaee and Salimi proposed a simple method based on the Bucherer reaction to attach -NH2 to GO.31 In that method, the reaction of ammonia and GO is catalyzed by sodium bisulfate. According to the surface analysis and material characterization, the authors establish presence of multiple -NH2 groups on the GRA surface. An improved electrocatalytic activity of the product was recorded. Alkylamino-FG with oxygen containing groups was reported to exhibit a good dispersion in organic solvents and a relatively strong interfacial interaction with organic polymers.39 Despite the experimental successes, the structural formula of amino-FG remains largely unknown. A number of theoretical papers focused on FG and its properties. Tetsuka et al. used amino groups to tune optically properties of GRA quantum dots,38 while Chakraborti and Pal explored amino-GRA flakes for application as anode material in Li-ion batteries.40 Multiple papers investigated hydrogenated41-43 and fluorinated35,

37, 44

GRA. More complex aspects of gRA

chemistry, such as oxygen reduction45 and Diels-Alder reaction,46 were studied.

ACS Paragon Plus Environment

Page 5 of 18

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

The Journal of Physical Chemistry Letters

In this letter, we employ ab initio gas-phase thermochemistry and long-scale finitetemperature electronic-structure-based molecular dynamics (MD) in order to characterize in detail the composition and structure of amino-FG. We address the following and related questions. What is the maximum theoretically stable concentration of amino groups? Can amino groups be grafted in ortho and meta positions in view of steric conditions? Can -NH2 exist both above and below the GRA plane simultaneously? How does the formation free energy depend on the number and location of -NH2 groups? We discover, quite surprisingly, a strong synergistic effect. The presence of multiple amino-groups stabilize FG and favor further amination. At the same time, a small concentration of amino-groups leads, in many cases, to unstable products. The groups can be attached on one or both sides of the graphene sheet. Steric hindrances play an important role in preventing complete GRA amination. The reported results can be verified experimentally, complementing the experimental efforts aimed at characterization of FG. Figure 1 depicts simulated models of amino- and ethylamino-FG, whereas Table 1 provides a full list of the considered systems. We investigated a variety of structures involving different attachment and substitution positions – such as ortho (1,2); meta (1,3); para (1,4); neighboring rings (2,7); different arrangements of three, four and six (ethyl)amino groups per single 6-membered carbon ring – using a small model sheet (graphene quantum dot). According to our previous (unpublished) benchmark studies, dependence of functionalization on the hybridization state and chemical environment is much stronger than on the size of the model. Thermodynamic

data

are

provided

for

the

last

stage

of

amination,

i.e. GRA + °NH2 ↔ GRA-NH2 or GRA-H + °NH2 ↔ GRA-NH2 + °H at 298.15 K and 1.0 bar. We do not consider °NH2 liberation reactions, which depend on a particular method of synthesis, employed catalysts and available nitrogen source. Those reactions were well studied upon derivation of each particular method. If the last stage of amination exhibits a sufficiently positive

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 18

Gibbs free energy, such a reaction is, most likely, impossible. According to Table 1, all amination reactions bring modest changes in the free energy. In comparison, formation of most covalent bonds brings ca. 1000 kJ mol-1. Many of the considered amination positions appear thermodynamically forbidden, e.g. the ortho and meta positions, the edge position, etc. However, the para position is allowed, irrespective of the substituent. Grafting the ethylamino groups is generally more energetically favorable. For instance, attachment of the lone amino group to a central atom of GRA costs 13 kJ mol-1, whereas attachment of the lone ethylamino group at the same position brings 3 kJ mol-1. All presented values (Tables 1-2) were normalized per mole of -NH2. Substitution of °H by °NH2 or °C2H5-NH2 at the rim carbon atom without a catalyst is forbidden in line with the known amination of arenes.

Figure 1. Amino- and alkylamino-FG models: (a) numbering of the considered attachment sites; (b) 1,4-diaminographene; (c) 1,4-diethylaminographene. The functionalization at the “edge” position involves substitution of the terminating hydrogen atom of the functional group.

A clear synergistic effect is seen upon amination. While lone amino groups at the surface of GRA are unstable, addition of more -NH2 stabilizes the product. Two, three and four groups lead to systematically more favorable FG. Six amino groups per ring – four above the GRA plane and two below the GRA plane – constitute an unexpectedly favorable FG sample, -85 kJ mol-1. Note that the free energy is given per number of -NH2. A natural

ACS Paragon Plus Environment

Page 7 of 18

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

The Journal of Physical Chemistry Letters

supposition is that the synergistic stabilization occurs due to Coulombic attraction between -NH2 and electronically polarized GRA surface. This hypothesis is supported by several facts. First, amino-FG is more favorable than ethylamino-FG, -85 vs. -64 kJ mol-1. Second, formation of the ‘2,7’ amino-FG and ethylamino-FG systems is unfavorable, because the distance between these positions is significant to weaken the stabilizing Coulombic attraction. Many other interesting correlations can be derived from Tables 1-2.

Table 1. Considered configurations of amino- and alkylamino-FG, and their key characteristics Positions○ Stability□ ∆GM11, kJ mol-1 rmax(C-C),◊ nm 1, 2 stable +23 0.160 1, 3 stable +67 0.152 1, 4 stable -41 0.152 2, 7 stable +56 0.152 1, 2, 3 unstable ― ― 1, 2, 4 stable -33 0.160 1, 2, 4, 5 stable -27 0.160 1,2,4,5 (above) 8 -NH2 stable -85 0.161 3,10 (below) 9 -NH2 edge stable +30 0.143 10 -NH2 1 stable +13 0.151 11 -C2H5NH2 1, 2 unstable ― ― 12 -C2H5NH2 1, 3 unstable ― ― 13 -C2H5NH2 1, 4 stable -34 0.152 14 -C2H5NH2 2, 7 stable +36 0.151 15 -C2H5NH2 1, 3, 5 unstable ― ― 1,4 (above) 16 -C2H5NH2 stable -64 0.154 2,5 (below) 17 -C2H5NH2 edge stable +43 0.145 18 -C2H5NH2 1 stable -3 0.151 ○ See Figure 1 for the site designation. □ Stability indicates integrity of the FG during the 200 ps long PM7-MD simulations at 500 K, as described in the methodology section. ◊ The length of the longest carbon-carbon (C-C) bond neighboring the grafted (ethyl)amino group. The longest bond is provided, because it directly reflects FG stability of FG. The C-C bond is 0.141 nm long in pristine GRA. # 1 2 3 4 5 6 7

Group -NH2 -NH2 -NH2 -NH2 -NH2 -NH2 -NH2

Table 2 decomposes Gibbs free energies of the corresponding reactions into the enthalpic (∆H) and entropic (-T∆S) components. Entropy of the GRA amination is unfavorable in all cases. The entropic component, 298.15 K×∆SM11, of the reaction at the rim carbon atom amounts ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 18

to 16 (amino-FG) and 18 kJ mol-1 (ethylamino-FG). Compare this value to that for amination of benzene, 18 kJ mol-1, obtained using the same level of theory. Similarity of the entropic components confirms our earlier supposition that enthalpy (including Coulombic attraction) ultimately determines the direction of the amination process in the considered cases. Note, however, that the absolute value of the entropy change plays an extremely important role for this reaction at room conditions, constituting more than half of the reaction free energy in many cases.

Table 2. Standard enthalpic contribution (∆H), entropic contribution (-T∆S) and entropy (∆S) for the investigated amination reactions Positions ∆HM11, kJ mol-1 -T∆SM11, kJ mol-1 ∆SM11, J mol-1 K-1 ―NH2 1,2 -26 +50 -167 1,3 +18 +49 -165 1,4 -91 +50 -168 2,7 +7.0 +49 -163 1,2,4 -83 +50 -168 1,2,4,5 -78 +51 -173 1,2,4,5,3,6 -138 +52 -176 edge +14 +16 -54 ‘1’ -34 +46 -156 ―C2H5-NH2 ‘1,4’ -88 +54 -180 ‘2,7’ -19 +55 -185 ‘1,4,2,5’ -120 +56 -189 edge +26 +18 -60 ‘1’ -56 +53 -178

Thermal stability was investigated using PM7-MD simulations. The geometries of the constructed systems were optimized using PM7 and gradually heated from 0 to 500 K. 200 ps long trajectories were used to derive distance distributions (Figures 2-3) and dihedral distributions (Figure 4). Note that we study thermal stabilities of the FG samples in vacuum. If any aggressive solvent is present, the reported stabilities may alter drastically. The reported trends are expected to hold in inert solvents.

ACS Paragon Plus Environment

Page 9 of 18

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

The Journal of Physical Chemistry Letters

It was found that the ‘1,2,3’ amino-FG is unstable. This occurs for steric reasons: amino groups do not have enough space on the same side of the GRA plane. An average C-C-N angle, with both carbon atoms belonging to GRA, is 110 degrees in the ortho-FG, 107 degrees in the meta-FG, 106 degrees in the para-FG. The C-C-N angle is 106 degrees also in the ‘2,7’ aminoFG system. Geometry analysis also explains the thermodynamic observations (Tables 1-2). Note that the ‘1,2’ amino-FG system is perfectly stable over the entire MD run, even though its formation out of pristine GRA is not thermodynamically allowed. A more sophisticated synthetic approach should be used. Steric problems undermine thermal stabilities of the ‘1,2’ and ‘1,3’ and ‘1,3,5’ ethylaminoFG systems as well. We anticipate that longer alkylamino chains deliver new spatial restrictions, since a natural shape of the hydrocarbon chain has to accommodate to the available space. Calculations involving other groups can be readily performed applying our currently established methodology.

Figure 2. Distributions of the carbon-nitrogen distances (covalent bond lengths) in the amino-FG in the course of the PM7-MD simulations at 500 K in vacuum. See Table 1 for system designation.

The length of the carbon-nitrogen bond in aminobenzene is 0.144 nm. The corresponding distances in FG are different and slightly dependent on the position (Figure 2). When –NH2 gets linked at the ‘edge’ position, the bond length is 0.139 nm. If the attachment occurs at any other position (Table 1), the bond length is ca. 0.150 nm. The ‘1,4’ (para) position allows for ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 18

somewhat shorter C-N bonds than the meta and ortho positions. Compare these lengths to the lengths of the C-C bonds in the 6-membered rings (Table 1) for different positions of –NH2.

Figure 3. Distributions of the carbon-carbon (covalent bond) and carbon-nitrogen distances in the ethylamino-FG during the course of the PM7-MD simulations at 500 K in vacuum. See Table 1 for system designation.

Figure 3 characterizes the C-C and C-N distances in the case of ethylamino-FG. Note that the C-N distances correspond to the nitrogen atoms of –NH2 and the closest carbon atom of GRA. That is, these distances are not covalent bond lengths, unlike the C-C distances, which were calculated between the GRA plane and the closest ethylene group of the linked -C2H5NH2 moiety. The ethylamino groups are very mobile at 500 K, sampling 0.30-0.45 nm C-N distances to the GRA plane. A larger number of the -C2H5NH2 moieties makes each of them stiffer (Figure 3, right). The grafted –NH2 perturbs planarity of the GRA sheet considerably. The effect is stronger in the case of amino-FG than in the case of ethylamino-FG (Figure 4), because the nitrogencarbon groups are polar. In turn, the carbon-carbon bonds in ethylamino-FG are nearly nonpolar. Importantly, a clear dependence is observed between the number of –NH2 groups and nonplanarity of the corresponding 6-membered ring. An improper dihedral of the pristine GRA was also plotted for comparison, Figure 4. On the average, pristine GRA and GRA with a single – NH2 at the edge are planar. Thermal fluctuations are significant; deviations of ±9 degrees are ACS Paragon Plus Environment

Page 11 of 18

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

The Journal of Physical Chemistry Letters

seen in the pristine system at 500 K. The analogous deviations are somewhat larger in the edgeamino-FG. Although the ring deforms very significantly due to functionalization, the structure is maintained. The distribution of improper dihedrals in the ‘1,2,4,5&3,6’ system is quite narrow, suggesting nearly harmonic fluctuations from planarity during MD. One can expect that long range interactions present between polar groups in large graphene sheets will lead average out geometry fluctuations, reducing somewhat the non-planarity observed with the current systems.

Figure 4. Deviation from planarity of the 6-membered ring due to (left) amino- and (right) ethylamino-functionalization. The depicted improper dihedrals were computed using positions 2, 10, 4 and 1, see Figure 1.

Figure 5 depicts an interesting product of amination, which was obtained during the PM7MD simulations during the first five picoseconds. The starting configuration was 1,2ethylaminegraphene. This molecule readily cleaves two hydrogen molecules, formed out of two amino groups, whereas unsaturated nitrogen atoms create a double bond, thereby closing the ring. The 1,2-diethylamino-FG is, therefore, unstable, but gives an immediate rise to another stable compound.

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 18

Figure 5. The isomerization product obtained from 1,2-diethylaminographene (see Figure 1 for site designation).

Following the practice of Adeli and coworkers working with carbon nanotubes,47 we investigated hydrogen bonds between the nitrogen atom of the amino group and the hydrogen atom(s) of the neighboring amino groups. It was found that a single hydrogen bond, 2.0-2.2 Å long, was formed when the amino groups were attached to the same carbon ring. The hydrogen bond is longer in the ortho and para positions, ca. 2.2 Å, than in the meta position, ca. 2.0 Å. The ortho position is less favorable for hydrogen bonding due to a nitrogen-nitrogen electron-core repulsion, with the equilibrium distance of 2.6 Å. In turn, the meta position allows the nitrogennitrogen distance to be 2.9 Å, i.e. this distance is closer to the van der Waals diameter of the nitrogen atom. Hydrogen bonding between the grafted groups makes the free energy of the respective configuration more favorable. Figure 6 compares the accuracy of PM7-MD to the accuracy of the M11 hybrid density functional theory48 (HDFT) Born-Oppenheimer molecular dynamics (BOMD) for the ‘1,2,4,5’ system, i.e. one of the most complicated amino-FG. The derived pair correlation functions for the carbon-carbon, carbon-nitrogen and nitrogen-nitrogen atom pairs suggest minimal discrepancies between these methods, while all relevant peaks are reproduced perfectly.

ACS Paragon Plus Environment

Page 13 of 18

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

The Journal of Physical Chemistry Letters

Figure 6. C-C, C-N and N-N pair correlation functions, computed in the 1,2,4,5tetraaminographene using PM7-MD and M11-HDFT BOMD.

To recapitulate, we used the M11 HDFT electronic-structure calculations and PM7–MD simulations to unravel regularities and peculiarities of structure and thermodynamic properties of amino- and ethylamino-FG. We discovered an unusual, synergistic amination of GRA and characterized the products in detail. A large number of amino-groups stabilize functionalized graphene favoring further amination. In contrast, a small number of amino-groups typically leads to unstable products. The functional groups can be attached either on the same side or on different sides of the GRA plane. Steric hindrances play an important role in preventing complete amination and formation of the ultimate FG, Cx(NH2)x. The amination reactions bring modest changes in the Gibbs free energy. Many amination positions are thermodynamically forbidden, including ortho, meta, and edge positions. The para position is allowed, irrespective of the substituent. Grafting a single ethylamino group is more energetically favorable compared to a single amino group. Deformation of functionalized GRA is proportional to the number of attached groups. Testable experimentally, the presented results are particularly useful for structure determination experiments, and assist in new organic syntheses.

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 18

Methodology Thermodynamic potentials for the selected reactions were derived using molecular partition functions. The wave functions were calculated at the M11/6-311++G** level (hybrid density functional theory),48 followed by calculation of force constants and resulting vibrational frequencies. The self-consistent field convergence criterion was set to 10-8 Hartree. Prior to frequency analysis, all molecular geometries were thoroughly optimized following evolution of energy gradients. The reported enthalpy, entropy and Gibbs free energy correspond to the gasphase approximation. According to the developers, M11 provides an improved structure, electronic and thermochemical properties, as compared to earlier functionals.48 The thermal stabilities of amino- and alkylamino-functionalized graphene were investigated using the PM7-MD simulations49-53 at 500 K. After an initial geometry optimization, every system was gradually heated from 0 to 500 K with a pace of 0.025 K/fs. 200 ps production runs were used to investigate the dependence of FG structure and stability on –NH2 positions. The equations-of-motion were propagated with a 0.5 fs time-step using the velocity Verlet algorithm. The system temperature was maintained by means of the Andersen thermostat.54 PM7-MD was successfully applied to a set of problems in chemistry including ion-molecular systems,51 competitive solvation,52 nanoparticles,50 gas capture,53 etc. The present computations were parallelized using 4 processor cores with a speed up factor of 2.4. The M11 HDFT-BOMD simulation48 was used to record a 3.000 ps long molecular trajectory. The 6-31G* atom-centered split-valence double-zeta polarized Pople basis set was used. The constant temperature, 500 K, was maintained via a simple velocity rescaling every 50 integration steps to obtain the target temperature. The trajectory was propagated with a 0.5 fs time-step. The PM7-MD simulations were performed using the in-home code (developed by V.V.C). Thermochemistry and HDFT-BOMD were calculated by the GAMESS quantum chemistry ACS Paragon Plus Environment

Page 15 of 18

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

The Journal of Physical Chemistry Letters

program suite. Visualization of the PDB files was performed by the VMD 1.9 software55 and the Gabedit 2.4 software.56

Acknowledgments V.V.C. is CAPES fellow. OVP acknowledges support of the US Department of Energy grant DE-SC0014429.

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 18

REFERENCES (1) Kuila, T.; Bose, S.; Mishra, A. K.; Khanra, P.; Kim, N. H.; Lee, J. H., Chemical Functionalization of Graphene and Its Applications. Prog Mater Sci 2012, 57, 1061-1105. (2) Prezhdo, O. V., Graphene - the Ultimate Surface Material. Surf Sci 2011, 605, 1607-1610. (3) Bourgeteau, T.; Le Vot, S.; Bertucchi, M.; Derycke, V.; Jousselme, B.; Campidelli, S.; Cornut, R., New Insights into the Electronic Transport of Reduced Graphene Oxide Using Scanning Electrochemical Microscopy. J Phys Chem Lett 2014, 5, 4162-4166. (4) Long, R., Understanding the Electronic Structures of Graphene Quantum Dot Physisorption and Chemisorption onto the Tio2 (110) Surface: A First-Principles Calculation. Chemphyschem 2013, 14, 579582. (5) Duque, J. G.; Chen, H.; Swan, A. K.; Shreve, A. P.; Kilina, S.; Tretiak, S.; Tu, X. M.; Zheng, M.; Doorn, S. K., Violation of the Condon Approximation in Semiconducting Carbon Nanotubes. Acs Nano 2011, 5, 5233-5241. (6) Basiuk, E. V.; Martinez-Herrera, M.; Alvarez-Zauco, E.; Henao-Holguin, L. V.; Puente-Lee, I.; Basiuk, V. A., Noncovalent Functionalization of Graphene with a Ni(Ii) Tetraaza[14] Annulene Complex. Dalton T 2014, 43, 7413-7428. (7) Chen, H.; Golder, M. R.; Wang, F.; Doorn, S. K.; Jasti, R.; Tretiak, S.; Swan, A. K., Raman-Active Modes of Even-Numbered Cycloparaphenylenes: Comparisons between Experiments and Density Functional Theory (Dft) Calculations with Group Theory Arguments. J Phys Chem C 2015, 119, 28792887. (8) Wang, H. B.; Xie, M. S.; Thia, L.; Fisher, A.; Wang, X., Strategies on the Design of Nitrogen-Doped Graphene. J Phys Chem Lett 2014, 5, 119-125. (9) Long, R.; English, N. J.; Prezhdo, O. V., Photo-Induced Charge Separation across the GrapheneTio2 Interface Is Faster Than Energy Losses: A Time-Domain Ab Lnitio Analysis. J Am Chem Soc 2012, 134, 14238-14248. (10) Habenicht, B. F.; Kalugin, O. N.; Prezhdo, O. V., Ab Initio Study of Phonon-Induced Denhasina of Electronic Excitations in Narrow Graphene Nanoribbons. Nano Lett 2008, 8, 2510-2516. (11) Hammouri, M.; Jha, S. K.; Vasiliev, I., First-Principles Study of Graphene and Carbon Nanotubes Functionalized with Benzyne. J Phys Chem C 2015, 119, 18719-18728. (12) Georgakilas, V.; Bourlinos, A. B.; Zboril, R.; Steriotis, T. A.; Dallas, P.; Stubos, A. K.; Trapalis, C., Organic Functionalisation of Graphenes. Chem Commun 2010, 46, 1766-1768. (13) Mishra, A. K.; Ramaprabhu, S., Functionalized Graphene-Based Nanocomposites for Supercapacitor Application. J Phys Chem C 2011, 115, 14006-14013. (14) Bekyarova, E.; Itkis, M. E.; Ramesh, P.; Berger, C.; Sprinkle, M.; de Heer, W. A.; Haddon, R. C., Chemical Modification of Epitaxial Graphene: Spontaneous Grafting of Aryl Groups. J Am Chem Soc 2009, 131, 1336-1337. (15) Antipina, L. Y.; Sorokin, P. B., Converting Chemically Functionalized Few-Layer Graphene to Diamond Films: A Computational Study. J Phys Chem C 2015, 119, 2828-2836. (16) Pan, X. Y.; Yang, M. Q.; Tang, Z. R.; Xu, Y. J., Noncovalently Functionalized Graphene-Directed Synthesis of Ultralarge Graphene-Based Tio2 Nanosheet Composites: Tunable Morphology and Photocatalytic Applications. J Phys Chem C 2014, 118, 27325-27335. (17) Wu, Q. Y.; Lan, J. H.; Wang, C. Z.; Xiao, C. L.; Zhao, Y. L.; Wei, Y. Z.; Chai, Z. F.; Shi, W. Q., Understanding the Bonding Nature of Uranyl Ion and Functionalized Graphene: A Theoretical Study. J Phys Chem A 2014, 118, 2149-2158. (18) Vovusha, H.; Sanyal, S.; Sanyal, B., Interaction of Nucleobases and Aromatic Amino Acids with Graphene Oxide and Graphene Flakes. J Phys Chem Lett 2013, 4, 3710-3718. (19) Malic, E.; Appel, H.; Hofmann, O. T.; Rubio, A., Forster-Induced Energy Transfer in Functionalized Graphene. J Phys Chem C 2014, 118, 9283-9289. (20) Dai, J.; Zhao, Y.; Wu, X. J.; Zeng, X. C.; Yang, J. L., Organometallic Hexahapto-Functionalized Graphene: Band Gap Engineering with Minute Distortion to the Planar Structure. J Phys Chem C 2013, 117, 22156-22161.

ACS Paragon Plus Environment

Page 17 of 18

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

The Journal of Physical Chemistry Letters

(21) Colherinhas, G.; Fileti, E. E.; Chaban, V. V., The Band Gap of Graphene Is Efficiently Tuned by Monovalent Ions. J Phys Chem Lett 2015, 6, 302-307. (22) Kandagal, V. S.; Pathak, A.; Ayappa, K. G.; Punnathanam, S. N., Adsorption on EdgeFunctionalized Bilayer Graphene Nanoribbons: Assessing the Role of Functional Groups in Methane Uptake. J Phys Chem C 2012, 116, 23394-23403. (23) Al-Aqtash, N.; Vasiliev, I., Ab Initio Study of Boron- and Nitrogen-Doped Graphene and Carbon Nanotubes Functionalized with Carboxyl Groups. J Phys Chem C 2011, 115, 18500-18510. (24) Lin, Z. Y.; Liu, Y.; Yao, Y. G.; Hildreth, O. J.; Li, Z.; Moon, K.; Wong, C. P., Superior Capacitance of Functionalized Graphene. J Phys Chem C 2011, 115, 7120-7125. (25) Bissett, M. A.; Tsuji, M.; Ago, H., Mechanical Strain of Chemically Functionalized Chemical Vapor Deposition Grown Graphene. J Phys Chem C 2013, 117, 3152-3159. (26) Stankovich, S.; Piner, R. D.; Nguyen, S. T.; Ruoff, R. S., Synthesis and Exfoliation of IsocyanateTreated Graphene Oxide Nanoplatelets. Carbon 2006, 44, 3342-3347. (27) Ramanathan, T.; Fisher, F. T.; Ruoff, R. S.; Brinson, L. C., Amino-Functionalized Carbon Nanotubes for Binding to Polymers and Biological Systems. Chem Mater 2005, 17, 1290-1295. (28) Punckt, C.; Pope, M. A.; Aksay, I. A., On the Electrochemical Response of Porous Functionalized Graphene Electrodes. J Phys Chem C 2013, 117, 16076-16086. (29) Zhang, F.; Jiang, H. Y.; Li, X. Y.; Wu, X. T.; Li, H. X., Amine-Functionalized Go as an Active and Reusable Acid-Base Bifunctional Catalyst for One-Pot Cascade Reactions. Acs Catal 2014, 4, 394-401. (30) Mungse, H. P.; Khatri, O. P., Chemically Functionalized Reduced Graphene Oxide as a Novel Material for Reduction of Friction and Wear. J Phys Chem C 2014, 118, 14394-14402. (31) Navaee, A.; Salimi, A., Efficient Amine Functionalization of Graphene Oxide through the Bucherer Reaction: An Extraordinary Metal-Free Electrocatalyst for the Oxygen Reduction Reaction. Rsc Adv 2015, 5, 59874-59880. (32) Colherinhas, G.; Fileti, E. E.; Chaban, V. V., Can Inorganic Salts Tune Electronic Properties of Graphene Quantum Dots? Phys Chem Chem Phys 2015, 17, 17413-17420. (33) Ayissi, S.; Charpentier, P. A.; Farhangi, N.; Wood, J. A.; Palotas, K.; Hofer, W. A., Interaction of Titanium Oxide Nanostructures with Graphene and Functionalized Graphene Nanoribbons: A Dft Study. J Phys Chem C 2013, 117, 25424-25432. (34) Nelson, T. R.; Prezhdo, O. V., Extremely Long Nonradiative Relaxation of Photoexcited Graphane Is Greatly Accelerated by Oxidation: Time-Domain Ab Initio Study. J Am Chem Soc 2013, 135, 3702-3710. (35) Ueta, A.; Tanimura, Y.; Prezhdo, O. V., Infrared Spectral Signatures of Multilayered SurfaceFluorinated Graphene: A Molecular Dynamics Study. J Phys Chem C 2012, 116, 8343-8347. (36) Nelson, T.; Zhang, B.; Prezhdo, O. V., Detection of Nucleic Acids with Graphene Nanopores: Ab Initio Characterization of a Novel Sequencing Device. Nano Lett 2010, 10, 3237-3242. (37) Ueta, A.; Tanimura, Y.; Prezhdo, O. V., Infrared Spectral Signatures of Surface-Fluorinated Graphene: A Molecular Dynamics Study. J Phys Chem Lett 2012, 3, 246-250. (38) Tetsuka, H.; Asahi, R.; Nagoya, A.; Okamoto, K.; Tajima, I.; Ohta, R.; Okamoto, A., Optically Tunable Amino-Functionalized Graphene Quantum Dots. Adv Mater 2012, 24, 5333-5338. (39) Wu, D. Q.; Zhang, F.; Liu, P.; Feng, X. L., Two-Dimensional Nanocomposites Based on Chemically Modified Graphene. Chem-Eur J 2011, 17, 10804-10812. (40) Chakraborti, H.; Pal, S. K., Assessment of Amine Functionalized Graphene Nanoflakes for Anode Materials in Li-Ion Batteries: An Ab Initio Study. Chemical Physics Letters 2014, 600, 118-122. (41) Nguyen, M. T.; Phong, P. N.; Tuyen, N. D., Hydrogenated Graphene and Hydrogenated Silicene: Computational Insights. Chemphyschem 2015, 16, 1733-1738. (42) Valencia, H.; Gil, A.; Frapper, G., Trends in the Hydrogen Activation and Storage by Adsorbed 3d Transition Metal Atoms onto Graphene and Nanotube Surfaces: A Dft Study and Molecular Orbital Analysis. Journal of Physical Chemistry C 2015, 119, 5506-5522. (43) Ulman, K.; Bhaumik, D.; Wood, B. C.; Narasimhan, S., Physical Origins of Weak H-2 Binding on Carbon Nanostructures: Insight from Ab Initio Studies of Chemically Functionalized Graphene Nanoribbons. Journal of Chemical Physics 2014, 140. (44) Xiao, J.; Meduri, P.; Chen, H. H.; Wang, Z. G.; Gao, F.; Hu, J. Z.; Feng, J.; Hu, M.; Dai, S.; Brown, S., et al., Energetics of Defects on Graphene through Fluorination. Chemsuschem 2014, 7, 1295-1300.

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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 18

(45) Favaro, M.; Ferrighi, L.; Fazio, G.; Colazzo, L.; Di Vaentin, C.; Durante, C.; Sedona, F.; Gennaro, A.; Agnoli, S.; Granozzi, G., Single and Multiple Doping in Graphene Quantum Dots: Unraveling the Origin of Selectivity in the Oxygen Reduction Reaction. Acs Catalysis 2015, 5, 129-144. (46) Denis, P. A., Organic Chemistry of Graphene: The Diels-Alder Reaction. Chemistry-a European Journal 2013, 19, 15719-15725. (47) Adeli, M.; Soleyman, R.; Beiranvand, Z.; Madani, F., Carbon Nanotubes in Cancer Therapy: A More Precise Look at the Role of Carbon Nanotube-Polymer Interactions. Chem Soc Rev 2013, 42, 52315256. (48) Peverati, R.; Truhlar, D. G., Improving the Accuracy of Hybrid Meta-Gga Density Functionals by Range Separation. J Phys Chem Lett 2011, 2, 2810-2817. (49) Stewart, J. J. P., Optimization of Parameters for Semiempirical Methods Vi: More Modifications to the Nddo Approximations and Re-Optimization of Parameters. J Mol Model 2013, 19, 1-32. (50) Andreeva, N. A.; Chaban, V. V., Global Minimum Search Via Annealing: Nanoscale Gold Clusters. Chem Phys Lett 2015, 622, 75-79. (51) Chaban, V., Competitive Solvation of (Bis)(Trifluoromethanesulfonyl)Imide Anion by Acetonitrile and Water. Chem Phys Lett 2014, 613, 90-94. (52) Chaban, V., Competitive Solvation of the Imidazolium Cation by Water and Methanol. Chem Phys Lett 2015, 623, 76-81. (53) Chaban, V., The Thiocyanate Anion Is a Primary Driver of Carbon Dioxide Capture by Ionic Liquids. Chem Phys Lett 2015, 618, 89-93. (54) Andersen, H. C., Molecular-Dynamics Simulations at Constant Pressure and-or Temperature. J Chem Phys 1980, 72, 2384-2393. (55) Humphrey, W.; Dalke, A.; Schulten, K., Vmd: Visual Molecular Dynamics. J Mol Graph Model 1996, 14, 33-38. (56) Allouche, A. R., Gabedit-a Graphical User Interface for Computational Chemistry Softwares. J Comput Chem 2011, 32, 174-182.

ACS Paragon Plus Environment