Defective ... - ACS Publications

May 27, 2015 - Advanced Materials Department, Instituto Potosino de Investigación Científica y Tecnológica (IPICYT), Camino a la Presa San José. 2...
1 downloads 21 Views 2MB Size
Subscriber access provided by NEW YORK UNIV

Article

Differential Response of Doped/Defective Graphene and Dopamine to Electric Fields: A DFT Study Josue Ortiz-Medina, Florentino Lopez-Urias, Humberto Terrones, Fernando Jaime Rodríguez-Macías, Morinobu Endo, and Mauricio Terrones J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b04165 • Publication Date (Web): 27 May 2015 Downloaded from http://pubs.acs.org on June 4, 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 C 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 27

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

Differential Response of Doped/Defective Graphene and Dopamine to Electric Fields: A DFT Study J. Ortiz-Medina1,2, F. López-Urías2*, H. Terrones3, F. J. Rodríguez-Macías2,4, M. Endo5, M. Terrones1,5,6 1

Research Center for Exotic Nanocarbons (JST), Wakasato 4-17-1, Nagano 380-8553, Japan. Advanced Materials Department, IPICYT Camino a la Presa San José 2055, Col. Lomas 4ª sección, San Luis Potosí S.L.P., 78216, México. 3 Department of Physics, Applied Physics and Astronomy, Rensselaer Polytechnic Institute, Troy 12180, USA 4 Department of Fundamental Chemistry, Universidade Federal de Pernambuco, Av.Prof. Luiz Freire, Recife, PE, 50740-540, Brazil 5 Institute of Carbon Science and Technology, Shinshu University, 4-17-1 Wakasato, Nagano 380-8553, 2

Japan 6

Department of Physics, Chemistry, Materials Science and Engineering, Center for 2-Dimensional and Layered Materials and Materials Research Institute. The Pennsylvania State University, 104 Davey Lab., University Park, PA 16802-6300.

Keywords: doped graphene, dopamine, DFT, electric-field.

Abstract First-principles Density Functional Theory (DFT) calculations are performed on dopamine (DA)-graphene systems in the presence of an external electric field. The graphene lattice is also modified via substitutional boron- and nitrogen-doping, and via the introduction of defects (monovacancy and Thrower-Stone-Wales). Geometry optimization, electronic density of states, cohesive energy, electronic charge density, and wave functions are analyzed. Our results revealed that dopamine is anchored on the surface of graphene via a physisorption mechanism, and the cohesive strength varies as B-doped > N-doped > Vacancy defect > Thrower-StoneWales defect. Boron-doped graphene exhibits valence states with dopamine molecules, furthermore this system showed the strongest cohesive energy. When an electric field is applied, 1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

we observed shifts in the valence states near the Fermi level producing a decrease in the molecule-layer interaction. We envisage that the present contribution could help developing novel biosensors based on doped/defective graphene FET devices.

1. Introduction Graphene-based devices possess important advantages for detection of specific molecules. Several reports on theoretical1–3 and experimental4–9 research have demonstrated that graphene has outstanding properties such as high electrical conductivity, quantum Hall effect, and a band gap tunable by chemical and physical doping, including the application of external electric fields. Based on electrostatic interactions between different chemical species and graphene it is possible to efficiently sense molecules within graphene Field Effect Transistors (FET).10–12 In this context, Schedin et al.13 reported graphene-based sensors capable of detecting individual gas molecules. Such extremely high sensitivity has been attributed to charge transfer phenomena occurring at the graphene surface, combined with the large changes in carrier density caused by the injection of electrons or holes. Alternate explanations of the sensitivity of graphene consider chemical doping processes driven by electrostatic interactions between the graphene surface and charged sections of adsorbed molecules. For example, Nistor et al.14 refer to the 2D surface-induced electronegativity equalization produced by the conductive nature of graphene and special alignments between the adsorbate-graphene highest-occupied and lowest-unoccupied molecular orbitals (HOMO-LUMO).

2

ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27

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

The doping of graphene and other carbon nanomaterials can be a practical route to tailor some of their properties,15 as it can alter the available charge density by the incorporation of noncarbon atoms with excess electrons (e.g. nitrogen, for n-type doping) or electron deficiency (e.g. boron, for p-type doping).6,16–18 Doping of carbon nanomaterials also increases chemical activity near the doped sites due to localization of electronic states and reduced charge mobility caused by electron scattering, which can be useful in FET devices for high sensitivity detectors.19,20 Moreover, the electronic properties of graphene can also be changed when introducing defects within its lattice.21,22 For example, graphene lattice vacancies induce strong localization of electronic states close to the Dirac point,23 thus enhancing chemical reactivity similar to doping. In addition, theoretical studies show that 5-7 defects boundaries (i.e. geometrical flat arrangements of carbon atoms, consisting in consecutive pentagons and heptagons) in between two regular graphene hexagonal lattices exhibit different electronic, magnetic and spin dependent conductance properties.24

Further advancements on graphene-based sensors development can greatly benefit from a thorough understanding of how analytes interact with doping or point defect on graphene. Such understanding could even enable tailoring the response of graphene sensors by controlling its electronic properties via defect/doping engineering. It could be argued that interactions with different substances produce different effects on the graphene electronic behavior depending on how doping or defects change its electronic properties. This in turn could be advantageous in detection of biomolecules, since such chemical species are more complex when compared to simple gaseous molecules. The feasibility of biosensors based on carbon nanomaterials has been

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

demonstrated by reports on interactions between graphene or carbon nanotubes with biomolecules such as aminoacids,25,26 neurotransmitters27,28 and drugs,29,30 but the use of doping and defect structures in such sensors has not yet been studied in detail.

In this work, we investigate the interactions of dopamine (DA) with pure (undoped), boron (B), and nitrogen (N) doped graphene, as well as defects including monovacancies and 5775 defects (or Thrower-Stone-Wales (TSW) defects) within the graphene sheet. We selected DA because it is an important biomolecule with strong regulation effects on a wide range of biological processes, including neural processes.31 In addition, DA is also representative of catecholamines, fundamental biomolecules with structures based on benzene rings with two hydroxyl groups and an amino-terminal side-chain. Dopamine is one of the most important neurotransmitters, implicated in several neurological disorders, including Parkinson’s disease, depression and addictions. Thus, reliable detection of DA is of great clinical importance. Several studies have used graphene based materials, including N-doped graphene, for electrochemical detection of dopamine, as recently reviewed.32 More recent reports33 show that nitrogen doped “graphene foam” 3D-nanostructures can result in highly efficient sensing of DA, down to nanomolar levels and with good selectivity. Studies for detection of dopamine via electrochemical techniques have shown that the HOPG surface may be active in oxidizing DA, even though it was previously believed that DA would only interact with edge sites during electrochemical reactions.34 This stresses the importance of understanding how DA interacts with graphene, and points to the need to assess whether electric fields enhance such interactions. In order to address this, we studied systematically the effect of

4

ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27

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

different doping and defect types on graphene-DA interactions, with and without an external electric field. To the best of our knowledge, this is the first study explore the combination of perpendicular E-fields with doped/defective graphene in order to follow biomolecule interactions, important when developing carbon-based biosensors, as well as to contribute towards the understanding of field effect transistor (FET) based biomolecule sensors using graphene.

2. Computational methods Electronic calculations of the interaction between dopamine and graphene were performed using density functional theory (DFT)35,36 in the framework of the general gradient approximation (GGA) within the Perdew-Burke-Ernzerhof (PBE) approach37 with a basis of linear combination of atomic orbitals (LCAO) as implemented in the SIESTA code.38,39 We used a double-ζ basis set plus polarization orbitals (DZP), the real-space grid used for charge and potential integration is equivalent to a plane wave cut-off energy of 250 Ry and the pseudo-potentials (pp) were built from 1, 3, 4, 5, and 6 valence electrons for the hydrogen (H: 1s1), boron (B: 2s22p1), carbon (C: 2s22p2), nitrogen (N: 2s22p3), and oxygen (O: 2s22p4) atoms respectively. The different systems were constructed using a 6x6 graphene supercell with minimum vertical cell dimension of 40 Å to avoid graphene layers interactions. All the calculations were performed using a MonkhorstPack grid of 25×25×1 k-points, and relaxed by conjugate gradient minimization until the maximum forces were lower than 0.06 eV Å-1, including control calculations using auxiliary Grimme potentials for accounting dispersion interactions (see Supporting Information). We also analyzed

the

electronic

structure

of

the

dopamine

5

ACS Paragon Plus Environment

molecule,

(2-(3,4-

The Journal of Physical Chemistry

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

dihydroxyphenyl)ethylamine, C8H11NO2 figure 1a), which has two hydroxyl groups and an ethylamine (-CH2-CH2-NH2) joined to a benzene aromatic ring. Figures 1b-f show the supercells of pristine-, doped-, and defective-graphene structures. The dopamine molecule is placed on the different graphene structures as illustrated in figure 1g. Finally, for analyzing the effects of perpendicular electric fields, a static and uniform external potential was defined for all the systems, with electric field oriented perpendicularly to the graphene surface (as indicated in figure 1g) with a magnitude of 0.1 V Å-1.

3. Results and discussion Figure 2 depicts the density of states (DOS) and wave functions around the Fermi level (EF) for DA, with the highest-occupied molecular orbital (HOMO, the molecule's valence orbitals) at 2.815 eV, whereas the lowest-unoccupied molecular orbital (LUMO) is located at 1.219 eV. Wave function isosurface plots corresponding to the HOMO-1, HOMO, LUMO, and LUMO+1 energy levels are also shown in figure 2, where it can be observed that the amino group shows electronic states mainly in the valence levels. It is noteworthy that the benzene ring possess states in both occupied and unoccupied levels that related to the π and π* states arising from its aromaticity, and are in fact the main source for the different interactions and charge transfer events with the graphene layer (see below).

In order to analyze the stability of the interactions between DA and the different doped/defective graphene layers, the systems were relaxed using a conjugated gradient algorithm. After geometry optimization, total energy and different electronic properties were calculated. The

6

ACS Paragon Plus Environment

Page 6 of 27

Page 7 of 27

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

cohesive energy (Ec) was calculated as: Ec = EG–DA– EDA– EG, where EG–DA is the energy for the complex system (graphene and dopamine), EDA for the isolated single dopamine molecule, and EG the energy of the undoped or doped graphene used in the calculations. The counterpoise correction for the basis set superposition error was applied for all the systems and considered in the energy calculations.40,41 Figure 3 shows the results for each particular system and when applying an E-field noticeable differences appear on the cohesion energy, which can be attributed to alterations in the electronic local distribution at the doping/defect sites resulting in different configurations of occupied states in which DA valence orbitals can interact. For systems analyzed without an E-field, B-doped graphene has the most stable interaction with DA (Ec = –436 meV), while N-doped graphene exhibits weaker DA physisorption (Ec = –192 meV). Less energetically favorable interactions occur between DA and vacancy defective graphene (Ec = –111 meV), and the most favorable defective graphene possess TSW defects (Ec = –55 meV). These defects disrupt π-π interactions between the aromatic ring of DA and the graphene sheet. An external E-field produces significant changes (ca. 140 meV) in Ec between DA with pristine or B-doped graphene. For N-doping the difference with or without E-field is negligible (Ec ≈ – 190 meV). Finally, for DA and graphenes with structural defects, monovacancy produces slight stability changes as a function of the E-field (∆Ec ≈ 67 meV) and much smaller ones for TSW defects (20 meV). Our results are comparable with published data on interactions between graphene and other biomolecules. For example, theoretical studies establish that for L-leucine, the adsorption energies (equivalent to our cohesion energies) can go from -170 to -310 meV, depending on the orientation of the aminoacid molecule over the graphene sheet,12 whereas for other organic molecules such as tetracyanoquinodimethane, these energies range from -151 to -

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

635 meV.10 Moreover, we show here how these cohesion energies can be tuned by applying an external E-field, especially for the undoped and B-doped graphene, an effect which has not been studied before.

In order to evaluate the relaxed structures of DA adsorbed on graphene, we calculated perpendicular distances from nitrogen (labeled as dN) and oxygen atoms (labeled as dO1 and dO2) to the graphene basal plane for all systems (see the model in figure 4). The distances are listed on table 1, where the effects of E-field can be compared. Distance changes are well related with cohesive energies: the B-doped graphene case, which was the most stable system, also shows the closest approach (0.7 Å) to graphene when the E-field is applied. However, N-doped and TSW-defected cases, which are insensitive to E-field, also show no DA-graphene distance variations. All these results could be explained by the fact that any effective interaction between DA and graphene involves a certain amount of charge transfer, which polarizes the system and produces local electrostatic repulsive-attractive forces.

Figure 5 shows the calculated DOS for undoped (a-b), B- (c-d) and N-doped (e-f) graphene with DA systems, with and without an applied E-field. The filled plots correspond to pure DA states, whereas the line plots represent the whole system (i.e. DA + graphene) states. Figure 6 depicts the corresponding molecular models with and without an applied E-field, including electron density isolines for each case: (a-b) undoped, (c-d) B-doped and (e-f) N-doped graphenes + DA. Only the electron density isolines for states ±0.25 eV around DA HOMO are shown for the sake of clarity. From DOS plots and electronic density distribution, it is clear that interactions

8

ACS Paragon Plus Environment

Page 8 of 27

Page 9 of 27

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

between DA and B-doped graphene are stronger when compared with N-doped or undoped graphene in terms of orbital mixing, even though both B-doped and undoped graphene upshift DA to energetic levels close to EF. The local electron deficiency produced by B-doping promotes charge transfer (see table 2) from DA to graphene, thus the electronic states spread as depicted in the DOS plots and in the overlapping electron density cloud in figs. 6c and d.

The DOS and electronic density isolines for defected graphene systems (figures 7 and 8) show that vacancy defects shift DA states closer to EF when compared to TSW defects (compare figures 7a-b against 7c-d), inducing some degree of orbital mixing, visible on figure 8a. Nevertheless, such interactions are smaller than with B-doping, since charge transfer is considerably smaller (0.13 vs 0.04 e-), and only when an E-field is applied (see table 2). For the TSW-defected graphene, as well as for undoped and N-doped cases, there is neither orbital mixing nor significant charge transfer, since DA states are too far from EF to promote strong electronic interactions.

Table 2 also shows how the E-field changes the intrinsic DA energy gap (∆Eg), i.e. the energy difference between DA HOMO and LUMO. While B-doped and vacancy-defected graphene reduce Eg of adsorbed DA molecules; pristine, N-doped and TSW-defected graphene increase it when the perpendicular E-field is applied. Table 2 also shows that pristine, N-doped and TSWdefected graphene produce the largest upshifts in DA HOMO energy when the electric field is applied. These differential effects on the electronic structure of adsorbed DA suggest an additional route for optical detection,42 since graphene doping appears to improve the optical

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

gap tuning of adsorbed molecules. In this case, for example, B-doping produces significant gap reduction when an E-field is applied, reducing the DA optical gap (found in the UV range) by 60 meV.

It is important to note that topological defects have been observed in graphene membranes,39 and we show here that such features can have an important effect in the interaction of graphene with organic molecules. Our calculations of DA adsorbed on graphene systems, with and without an applied electric field, display unique results not shown in previous studies where different types of organic molecules are placed in the surroundings of graphene sheets. For example, most of the reports that have addressed the changes occurring in the graphene’s electronic properties by electric fields focus on pure or doped graphene, with one or two layers typically2 and without studying the interaction with other molecules. The effects of applied electric fields over doped graphene systems towards adsorption/desorption have so far only been reported for simple molecules such as CO and H2.43,44 Different studies addressing biomolecules adsorption in graphene are commonly limited to the analysis of the interaction energies and changes in graphene properties.45–47 In here, we demonstrate that the use of electric fields, in combination with different type of doping or defect sites in graphene, has significant effects on the sensitivity for specific molecules, and the adsorption strength of molecules on the active sites, which can affect the performance of high sensitivity sensors.

4. Conclusions Our findings show that both doping and structural defects have important effects in the

10

ACS Paragon Plus Environment

Page 10 of 27

Page 11 of 27

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

interaction of graphene with biomolecules. This could be important in experimental studies where it is commonly assumed that the graphene sheets are free from defects, and also points the way to a potentially simple way of controlling the selectivity and sensitivity of graphene sensors via defect engineering. In summary, using DFT-calculations we investigated the dopamine molecule (DA) interacting with pristine-graphene, B-, N-doped graphene, vacancy- and TSW-defects within graphene. The electric field effects on the adsorption energy were analyzed. Our DFT results revealed a physisorption mechanism for DA over the different graphene sheets, attributable to relatively stable interactions between π electrons of the DA aromatic ring and the graphene's surface. Interestingly, the DA-graphene interaction can be controlled by an external E-field, especially for pristine and B-doped graphene, whereas N-doped graphene-DA interactions are insensitive to it. Defective graphene can be utilized as well for energetically stable adsorption of DA, with a smaller sensitivity to E-fields when compared to doped graphene. All these phenomena are related to the way DA valence levels mix with graphene's valence band, which is stronger when localized electronic states are available below the graphene EF. These results also demonstrate the feasibility of using doped and defective graphene based FET or optical devices for the development of biosensors, for DA and potentially other biomolecules.

Corresponding Author *To whom correspondence should be addressed: [email protected]. Tel/Fax: +52 4448342000 (ext. 7238)/+52 4448342040 Acknowledgments 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

This work was sponsored by the JST-Japan by funding the Research Center for Exotic Nanocarbons, under the Japanese regional Innovation Strategy Program by the Excellence. M.T. acknowledges the financial support from the MURI project award No. FA9550–12–1–0471. FLU acknowledges the partial support by CONACYT-México grant 60218-F1. FJRM acknowledges CONACYT-Mexico grant CB-2008-SEP-107082 for parts of this work made at IPICYT and the support of Rede NANOBIOTEC-Brasil (Edital 04/CII-2008 CAPES/MEC), and of FACEPE (BFP-0132-1.06/13) for a stay as visiting professor at UFPE, as well as IPICYT for the leave of absence granted for this stay.

References (1)

Castro Neto, A. H.; Peres, N. M. R.; Novoselov, K. S.; Geim, A. K. The Electronic Properties of Graphene. Rev. Mod. Phys. 2009, 81, 109-162.

(2)

Champlain, J. G. A First Principles Theoretical Examination of Graphene-Based Field Effect Transistors. J. Appl. Phys. 2011, 109, 84515.

(3)

Wang, L.; Zheng, J.; Zhou, J.; Qin, R.; Li, H.; Mei, W.-N.; Nagase, S.; Lu, J. Tuning Graphene Nanoribbon Field Effect Transistors Via Controlling Doping Level. Theor. Chem. Acc. 2011, 130, 483-489.

(4)

Novoselov, K. S.; Geim, A. K.; Morozov, S. V; Jiang, D.; Zhang, Y.; Dubonos, S. V; Grigorieva, I. V; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666-669.

(5)

Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Dubonos, S. V.; Zhang, Y.; Jiang, D. Room-Temperature Electric Field Effect and Carrier-Type Inversion in Graphene Films. Science 2004, 306, 660-662.

(6)

Late, D. J.; Ghosh, A.; Subrahmanyam, K. S.; Panchakarla, L. S.; Krupanidhi, S. B.; Rao, C. N. R. Characteristics of Field-Effect Transistors Based on Undoped and B- and NDoped Few-Layer Graphenes. Solid State Commun. 2010, 150, 734-738. 12

ACS Paragon Plus Environment

Page 12 of 27

Page 13 of 27

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

(7)

Lin, M.-W.; Ling, C.; Zhang, Y.; Yoon, H. J.; Cheng, M. M.-C.; Agapito, L. A.; Kioussis, N.; Widjaja, N.; Zhou, Z. Room-Temperature High On/Off Ratio in Suspended Graphene Nanoribbon Field-Effect Transistors. Nanotechnology 2011, 22, 265201.

(8)

Reddy, D.; Register, L. F.; Carpenter, G. D.; Banerjee, S. K. Graphene Field-Effect Transistors. J. Phys. D-Applied Phys. 2011, 44, 313001.

(9)

Schwierz, F. Graphene Transistors. Nat. Nano. 2010, 5, 487-496.

(10)

Chi, M.; Zhao, Y.-P. First Principle Study of The Interaction And Charge Transfer Between Graphene And Organic Molecules. Comp. Mater. Sci. 2012, 56, 79-84.

(11)

Jiang, H. Chemical Preparation of Graphene-Based Nanomaterials and Their Applications in Chemical and Biological Sensors. Small 2011, 7, 2413-2427.

(12)

Qin, W.; Li, X.; Bian, W.-W.; Fan, X.-J.; Qi, J.-Y. Density Functional Theory Calculations And Molecular Dynamics Simulations of The Adsorption of Biomolecules on Graphene Surfaces. Biomaterials 2010, 31, 1007-1016.

(13)

Schedin, F.; Geim, A. K.; Morozov, S. V; Hill, E. W.; Blake, P.; Katsnelson, M. I.; Novoselov, K. S. Detection of Individual Gas Molecules Adsorbed on Graphene. Nat. Mater. 2007, 6, 652-655.

(14)

Nistor, R. A.; Newns, D. M.; Martyna, G. J. The Role of Chemistry in Graphene Doping for Carbon-Based Electronics. ACS Nano 2011, 5, 3096-3103.

(15)

Lv, R.; Li, Q.; Botello-Méndez, A. R.; Hayashi, T.; Wang, B.; Berkdemir, A.; Hao, Q.; Elías, A. L.; Cruz-Silva, R.; Gutiérrez, H. R.; et al. Nitrogen-Doped Graphene: Beyond Single Substitution and Enhanced Molecular Sensing. Sci. Rep. 2012, 2, 586.

(16)

Martins, T.; Miwa, R.; da Silva, A.; Fazzio, A. Electronic and Transport Properties of Boron-Doped Graphene Nanoribbons. Phys. Rev. Lett. 2007, 98.

(17)

Wei, D.; Liu, Y.; Wang, Y.; Zhang, H.; Huang, L.; Yu, G. Synthesis of N-Doped Graphene by Chemical Vapor Deposition and Its Electrical Properties. Nano Lett. 2009, 9, 1752-1758.

(18)

Yu, S. S.; Zheng, W. T.; Jiang, Q. Electronic Properties of Nitrogen-/Boron-Doped Graphene Nanoribbons With Armchair Edges. IEEE Trans. Nanotechnol. 2010, 9, 78-81.

(19)

Wei, X.-L.; Chen, Y.-P.; Liu, W.-L.; Zhong, J.-X. Enhanced Gas Sensor Based on Nitrogen-Vacancy Graphene Nanoribbons. Phys. Lett. A 2012, 376, 559-562.

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(20)

Cruz-Silva, E.; Barnett, Z. M.; Sumpter, B. G.; Meunier, V. Structural, Magnetic, and Transport Properties of Substitutionally Doped Graphene Nanoribbons From First Principles. Phys. Rev. B 2011, 83, 155445.

(21)

Araujo, P. T.; Terrones, M.; Dresselhaus, M. S. Defects and Impurities in Graphene-Like Materials. Mater. Today 2012, 15, 98-109.

(22)

Terrones, H.; Lv, R.; Terrones, M.; Dresselhaus, M. S. The Role of Defects and Doping in 2D Graphene Sheets and 1D Nanoribbons. Rep. Prog. Phys. 2012, 75, 062501.

(23)

Palacios, J.; Fernández-Rossier, J.; Brey, L. Vacancy-Induced Magnetism in Graphene and Graphene Ribbons. Phys. Rev. B 2008, 77.

(24)

Botello-Méndez, A. R.; Cruz-Silva, E.; López-Urías, F.; Sumpter, B. G.; Meunier, V.; Terrones, M.; Terrones, H. Spin Polarized Conductance in Hybrid Graphene Nanoribbons Using 5-7 Defects. ACS Nano 2009, 3, 3606-3612.

(25)

Wu, S.; Lan, X. Q.; Huang, F. F.; Luo, Z. Z.; Ju, H. X.; Meng, C. G.; Duan, C. Y. Selective Electrochemical Detection of Cysteine in Complex Serum by Graphene Nanoribbon. Biosens. Bioelectron. 2012, 32, 293-296.

(26)

Sun, W.; Bu, Y.; Wang, Y. Interaction Between Glycine/Glycine Radicals and Intrinsic/Boron-Doped (8,0) Single-Walled Carbon Nanotubes: A Density Functional Theory Study. J. Phys. Chem. B 2008, 112, 15442-15449.

(27)

Zhang, Y.; Yuan, R.; Chai, Y.; Li, W.; Zhong, X.; Zhong, H. Simultaneous Voltammetric Determination for DA, AA and NO2- Based on Graphene/Poly-Cyclodextrin/MWCNTs Nanocomposite Platform. Biosens. Bioelectron. 2011, 26, 3977-3980.

(28)

Alwarappan, S.; Erdem, A.; Liu, C.; Li, C.-Z. Probing the Electrochemical Properties of Graphene Nanosheets for Biosensing Applications. J. Phys. Chem. C 2009, 113, 88538857.

(29)

Chaban, V. V; Prezhdo, O. V. Water Boiling Inside Carbon Nanotubes: Toward Efficient Drug Release. ACS Nano 2011, 5, 5647-5655.

(30)

Ghalandari, B.; Monajjemi, M.; Mollaamin, F. Theoretical Investigation of Carbon Nanotube Binding to DNA in View of Drug Delivery. J. Comput. Theor. Nanosci. 2011, 8, 1212-1219.

(31)

Roe, D. L. From DOPA to Parkinson’s Disease: The Early History of Dopamine Research. J. Hist. Neurosci. 1997, 6, 291-301.

14

ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27

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

(32)

Pandikumar, A.; Soon How, G. T.; See, T. P.; Omar, F. S.; Jayabal, S.; Kamali, K. Z.; Yusoff, N.; Jamil, A.; Ramaraj, R.; John, S. A.; et al. Graphene and Its Nanocomposite Material Based Electrochemical Sensor Platform for Dopamine. RSC Adv. 2014, 4, 63296-63323.

(33)

Feng, X.; Zhang, Y.; Zhou, J.; Li, Y.; Chen, S.; Zhang, L.; Ma, Y.; Wang, L.; Yan, X. Three-Dimensional Nitrogen-Doped Graphene As An Ultrasensitive Electrochemical Sensor for The Detection of Dopamine. Nanoscale 2015, 7, 2427-2432.

(34)

Patel, A. N.; Tan, S. Y.; Miller, T. S.; MacPherson, J. V.; Unwin, P. R. Comparison and Reappraisal of Carbon Electrodes for The Voltammetric Detection of Dopamine. Anal. Chem. 2013, 85, 11755-11764.

(35)

Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133-A1138.

(36)

Hohenberg, P. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864-B871.

(37)

Perdew, J. P. Generalized Gradient Approximation Made Simple Phys. Rev. Lett. 1997, 78, 1396.

(38)

Soler, J. M.; Artacho, E.; Gale, J. D.; García, A.; Junquera, J.; Ordejón, P.; SánchezPortal, D. The SIESTA Method for Ab Initio Order-N Materials Simulation. J. Phys. Condens. Matter 2002, 14, 2745.

(39)

Junquera, J.; Paz, O.; Sánchez-Portal, D.; Artacho, E. Numerical Atomic Orbitals for Linear-Scaling Calculations. Phys. Rev. B 2001, 64, 235111.

(40)

Van Duijneveldt, F. B.; van Duijneveldt-van de Rijdt, J. G. C. M.; van Lenthe, J. H. State of the Art in Counterpoise Theory. Chem. Rev. 1994, 94, 1873-1885.

(41)

Note: For BSSE correction, additional energy calculations are performed considering ghost orbitals for individual dopamine and graphene systems, which prevents the orbital superposition error. Specifically, an extra energy term is subtracted to the typical cohesion energy calculation: Ec = EG-DA – EDA – EG - Ecc, where Ecc equals to the sum of the differences between energies of individual dopamine or graphene with and without ghost orbitals.

(42)

Wang, H. Study on Fluorescence Property of Dopamine and Determination of Dopamine by Fluorimetry. Talanta 2002, 57, 899-907.

(43)

Ao, Z. M.; Li, S.; Jiang, Q. Correlation of The Applied Electrical Field and CO Adsorption/Desorption Behavior on Al-Doped Graphene. Solid State Commun. 2010, 150, 680-683. 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(44)

Ao, Z. M.; Peeters, F. M. Electric Field Activated Hydrogen Dissociative Adsorption to Nitrogen-Doped Graphene. J. Phys. Chem. C 2010, 114, 14503-14509.

(45)

Ma, F.; Zhang, Z.; Jia, H.; Liu, X.; Hao, Y.; Xu, B. Adsorption of Cysteine Molecule on Intrinsic and Pt-Doped Graphene: A First-Principle Study. J. Mol. Struct-Theochem 2010, 955, 134-139.

(46)

Cazorla, C. Ab Initio Study of The Binding of Collagen Amino Acids to Graphene and A-Doped (A=H, Ca) Graphene. Thin Solid Films 2010, 518, 6951-6961.

(47)

Ou, L.; Luo, Y.; Wei, G. Atomic-Level Study of Adsorption, Conformational Change, and Dimerization of An α-Helical Peptide at Graphene Surface. J. Phys. Chem. B 2011, 115, 9813-9822.

16

ACS Paragon Plus Environment

Page 16 of 27

Page 17 of 27

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

FIGURES AND CAPTIONS Figure 1: J. Ortiz-Medina et al.

Figure 1: (a) Dopamine molecule. Unit cell representations for: (b) pristine or undoped graphene, (c) boron-doped graphene, (d) nitrogen-doped graphene, (e) vacancy-defect into graphene, and (f) Thrower-Stone-Wales (TSW) defect into graphene. (g) Optimized structure of adsorbed dopamine molecule on pristine graphene under the presence of an external electric field

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 2: J. Ortiz-Medina et al.

Figure 2: Electronic density of states (DOS) showing the energies distribution for an isolated dopamine (DA) molecule. Right side depicts the isosurfaces plots of wave functions corresponding to HOMO-1, HOMO, LUMO and LUMO+1.

18

ACS Paragon Plus Environment

Page 18 of 27

Page 19 of 27

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

Figure 3: J. Ortiz-Medina et al.

Figure 3: Cohesion energy (in meV) of the different systems studied with an applied electric field, perpendicular to the graphene plane with a magnitude of 0.1 V Å-1, (red squared marks) and without applied E-field (blue rounded marks). Top panels depict the corresponding relaxed structures.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 4: J. Ortiz-Medina et al.

Figure 4: Lateral and front molecular model views of dopamine (DA) molecule anchored on graphene systems. The different distances (dO1, dO2, dN) represent the separation between DA molecule (from O and N) and graphene basal plane. The indicated distances are reported table 1.

20

ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27

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

Figure 5: J. Ortiz-Medina et al.

Figure 5: Density of states (DOS) for dopamine (DA) molecule anchored on (a) pristine (pure), (c) B-doped and (e) N-doped graphene systems. For each system, results with an external applied electric field are also shown (b, d, f, respectively). The states for isolated DA are highlighted by the shaded DOS plot, where differential shiftings towards EF as the E-field is applied can be observed.

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 6: J. Ortiz-Medina et al.

Figure 6: Contour lines for electron density of (a, b) pure graphene and dopamine (DA) molecule, (c, d) B-doped graphene and DA, and (e, f) N-doped graphene and DA. The models show differential electron density distribution without electric field (a, c, and e) and with electric field applied (b, d, and f). The white arrows indicate the sites of the boron doping (c and d) and N doping (e and f) in the graphene lattice.

22

ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27

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

Figure 7: J. Ortiz-Medina et al.

Figure 7: Density of states (DOS) for dopamine (DA) molecule and (a) vacancy defected graphene and (c) Thrower-Stone-Wales defected graphene systems. The corresponding DOS plot for each system under an applied electric field is presented as well (b, and d, respectively). The DA states are highlighted by the shaded DOS plot, where differential shiftings towards EF are observed with the application of the E-field.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 8: J. Ortiz-Medina et al.

Figure 8: Contour lines for electron density of (a, b) vacancy defected graphene and dopamine (DA), and (c, d) Thrower-Stone-Wales defected graphene and DA. The models show differential electron density distribution without electric field (a and c) and with electric field (b and d). The white arrows indicate the site of the vacancy (a and b) and Thrower-Stone-Wales defect (c and d) into the graphene lattice.

24

ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27

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

TABLES AND CAPTIONS

Table 1: J. Ortiz-Medina et al. dO1 (Å)

dO2 (Å)

dN (Å)

System

Without Field

With E-Field

Without E-Field

With E-Field

Without E-Field

With E-Field

Undoped

2.84

2.87

2.9

2.88

4.87

4.93

B-doped

2.83

2.76

2.94

2.87

4.46

4.41

N-doped

2.87

2.91

2.91

2.92

4.69

4.69

Vacancy-defected

2.95

3.01

2.71

2.77

4.75

4.61

TSW-defected

2.86

2.87

2.85

2.86

4.85

4.85

Table 1: Distances (in Å) for three heteroatoms from dopamine (DA) molecule to graphene after geometry relaxation, measured perpendicularly to graphene basal plane (see figure 4). The distances are presented for each system, and compared with and without applied E-field.

25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 27

Table 2: J. Ortiz-Medina et al. System

Charge transfer DA→graphene (e-)

∆Eg (meV)

HOMO Upshiftings (meV)

without field

with E field

Undoped

-0.01

(*)

63

291

B-doped

0.15

0.13

-60

30

N-doped

(*)

0.01

90

270

Vacancy-defected

(*)

0.04

-15

120

TSW-defected

0.01

0.03

24

203

Table 2: Calculated results for electronic interaction between DA and graphene films. For the charge transfer (in e-) from dopamine (DA) molecule to graphene (*) indicates a negligible (