<?xml version='1.0'?>
<!DOCTYPE art SYSTEM 'http://www.biomedcentral.com/xml/article.dtd'>
<art><ui>1556-276X-7-268</ui><ji>1556-276X</ji><fm><dochead>Nano Express</dochead><bibl><title><p>Thinning and functionalization of few-layer graphene sheets by CF<sub>4</sub> plasma treatment</p></title><aug><au id="A1"><snm>Shen</snm><fnm>Chao</fnm><insr iid="I1"/><email>ssccsc@gmail.com</email></au><au id="A2" ca="yes"><snm>Huang</snm><fnm>Gaoshan</fnm><insr iid="I1"/><email>gshuang@fudan.edu.cn</email></au><au id="A3" ca="yes"><snm>Cheng</snm><fnm>Yingchun</fnm><insr iid="I2"/><email>yingchun.cheng@kaust.edu.sa</email></au><au id="A4"><snm>Cao</snm><fnm>Ronggen</fnm><insr iid="I1"/><email>caoronggen@fudan.edu.cn</email></au><au id="A5"><snm>Ding</snm><fnm>Fei</fnm><insr iid="I3"/><email>FDI@zurich.ibm.com</email></au><au id="A6"><snm>Schwingenschl&#246;gl</snm><fnm>Udo</fnm><insr iid="I2"/><email>udo.schwingenschlogl@kaust.edu.sa</email></au><au id="A7" ca="yes"><snm>Mei</snm><fnm>Yongfeng</fnm><insr iid="I1"/><email>yfm@fudan.edu.cn</email></au></aug><insg><ins id="I1"><p>Department of Materials Science, Fudan University, Shanghai, 200433, People's Republic of China</p></ins><ins id="I2"><p>PSE Division, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Kingdom of Saudi Arabia</p></ins><ins id="I3"><p>IBM Research&#8201;&#8722;&#8201;Z&#252;rich, S&#228;umerstrasse 4, R&#252;schlikon, CH-8803, Switzerland</p></ins></insg><source>Nanoscale Research Letters</source><section><title><p>Regular submissions</p></title></section><issn>1556-276X</issn><pubdate>2012</pubdate><volume>7</volume><issue>1</issue><fpage>268</fpage><url>http://www.nanoscalereslett.com/content/7/1/268</url><xrefbib><pubidlist><pubid idtype="doi">10.1186/1556-276X-7-268</pubid><pubid idtype="pmpid">22625875</pubid></pubidlist></xrefbib></bibl><history><rec><date><day>27</day><month>2</month><year>2012</year></date></rec><acc><date><day>25</day><month>4</month><year>2012</year></date></acc><pub><date><day>24</day><month>5</month><year>2012</year></date></pub></history><cpyrt><year>2012</year><collab>Shen et al.; licensee Springer.</collab><note>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</note></cpyrt><kwdg><kwd>Graphene sheets</kwd><kwd>Plasma</kwd><kwd>Thinning</kwd><kwd>Functionalization</kwd><kwd>Raman spectroscopy</kwd><kwd>Density functional theory</kwd><kwd>Dirac cone splitting</kwd></kwdg><abs><sec><st><p>Abstract</p></st><p>Structural changes of few-layer graphene sheets induced by CF<sub>4</sub> plasma treatment are studied by optical microscopy and Raman spectroscopy, together with theoretical simulation. Experimental results suggest a thickness reduction of few-layer graphene sheets subjected to prolonged CF<sub>4</sub> plasma treatment while plasma treatment with short time only leads to fluorine functionalization on the surface layer by formation of covalent bonds. Raman spectra reveal an increase in disorder by physical disruption of the graphene lattice as well as functionalization during the plasma treatment. The F/CF<sub>3</sub> adsorption and the lattice distortion produced are proved by theoretical simulation using density functional theory, which also predicts p-type doping and Dirac cone splitting in CF<sub>4</sub> plasma-treated graphene sheets that may have potential in future graphene-based micro/nanodevices.</p></sec><sec><st><p>PACS</p></st><p>81.05.ue; 73.22.Pr; 52.40.Hf.</p></sec></abs></fm><bdy><sec><st><p>Background</p></st><p>Graphene is one layer of C atoms, arranged in a hexagonal lattice <abbrgrp><abbr bid="B1">1</abbr><abbr bid="B2">2</abbr><abbr bid="B3">3</abbr></abbrgrp>. Since it was first produced by mechanical exfoliation in 2004 <abbrgrp><abbr bid="B4">4</abbr></abbrgrp>, graphene has been studied both theoretically and experimentally, and demonstrates highly attractive properties <abbrgrp><abbr bid="B5">5</abbr><abbr bid="B6">6</abbr><abbr bid="B7">7</abbr><abbr bid="B8">8</abbr><abbr bid="B9">9</abbr><abbr bid="B10">10</abbr></abbrgrp>. Especially, the Dirac equation predicts that unique electronic properties should arise from the hexagonal honeycomb lattice structure, making the electrons behave as massless relativistic fermions <abbrgrp><abbr bid="B1">1</abbr><abbr bid="B5">5</abbr><abbr bid="B9">9</abbr></abbrgrp>. The corresponding high mobility and velocity have great potential in future electronics <abbrgrp><abbr bid="B11">11</abbr><abbr bid="B12">12</abbr></abbrgrp>. However, several problems need to be solved before it can be ultimately employed in practical applications. For instance, the zero bandgap as well as bad wettability of graphene might cause problems in device fabrication <abbrgrp><abbr bid="B13">13</abbr><abbr bid="B14">14</abbr><abbr bid="B15">15</abbr></abbrgrp>. Under such circumstances, plasma treatment is considered to be one of the tricks to overcome the difficulties. The corresponding surface functionalization changes not only the surface status, but also the structure of graphene sheets <abbrgrp><abbr bid="B16">16</abbr><abbr bid="B17">17</abbr></abbrgrp>. Previous researches have already demonstrated that oxygen plasma has the ability to tune the properties of graphene sheets <abbrgrp><abbr bid="B14">14</abbr><abbr bid="B17">17</abbr></abbrgrp>. On the other hand, fluorine plasma, which may provide additional advantages <abbrgrp><abbr bid="B16">16</abbr><abbr bid="B18">18</abbr></abbrgrp>, is seldom experimentally investigated in detail, although a strong p-doping behavior was predicted <abbrgrp><abbr bid="B19">19</abbr><abbr bid="B20">20</abbr></abbrgrp>. In this work, structural changes of few-layer graphene sheets induced by CF<sub>4</sub> plasma treatment are studied by optical microscopy and Raman spectroscopy. Our results suggest an obvious thickness reduction effect in few-layer graphene sheets treated with CF<sub>4</sub> plasma as well as surface fluorine functionalization by formation of covalent bonds between the top graphene layer and the ions. The produced disorder in graphene lattice is well reflected in the Raman spectra, and the corresponding mechanism is studied theoretically. The results presented in this work provide a possible direction to obtain giant single-layer graphene sheets with necessary surface functionalization to realize the p-type doping level and to open the Dirac cone for future graphene-based micro/nanodevices.</p></sec><sec><st><p>Methods</p></st><p>The few-layer graphene sheets used in the current study were produced by mechanical exfoliation from highly ordered pyrolytic graphite and transferred on a silicon substrate covered with 300-nm SiO<sub>2</sub>. The samples were then cleaned by ultrasonication in acetone for 30 s to remove the residual of the scotch tape and unattached graphite pieces. The CF<sub>4</sub> plasma treatment was carried out in the reaction chamber of a Jupiter III reactive-ion etching setup. The samples were exposed to 0.8-Torr CF<sub>4</sub> under radio-frequency plasma (20 W) with different times ranging from 1 to 5 s with a step of 1 s. In our experiment, the number of layers in the few-layer graphene sheets can be estimated by a combination of optical microscopy and Raman spectroscopy. Raman spectra were taken on a Renishaw inVia micro-Raman spectrometer with the 514-nm line of an Ar<sup>+</sup> laser as light source (Shanghai, China). All the measurements were carried out at room temperature.</p></sec><sec><st><p>Results and discussion</p></st><p>As we know, morphology obtained by optical microscopy is a direct and convenient way to identify the few-layer graphene sheets. Optical images in Figure <figr fid="F1">1</figr> demonstrate the surface morphologies of the samples before and after CF<sub>4</sub> plasma treatment of different times. For a short plasma treatment (e.g., 2 s, see Figure <figr fid="F1">1</figr>a,b), no obvious change can be noticed at least in the optical images. However, with prolonged plasma treatment of up to 5 s, a thinning effect is obvious as the color contrast of the sample surface changes significantly. One can see in Figure <figr fid="F1">1</figr>c,d that dark purple regions appear lighter after treatment and some thin sheets can hardly be recorded by the camera (although they can be observed through eyepiece lens). During the process of plasma treatment, ions produced by an electromagnetic field attacked the sample surface and reacted with it. Therefore, the top layer of the sheet can be removed in longer treatment, which suggests an accumulation effect of the ion etching for such a low plasma power of 20 W.</p><fig id="F1"><title><p>Figure 1</p></title><caption><p>Optical images of samples before and after CF<sub>4</sub> plasma treatment</p></caption><text>
   <p><b>Optical images of samples before and after CF</b><sub><b>4</b></sub><b>plasma treatment.</b> ( <b>a</b>, <b>c</b>) Optical images demonstrating the morphologies of as-prepared samples before CF<sub>4</sub> plasma treatment. (<b>b</b>) Optical image of the sample in (a) after 3 s of CF<sub>4</sub> plasma treatment. (<b>d</b>) Optical image of the sample in (c) after 5 s of CF<sub>4</sub> plasma treatment.</p>
</text><graphic file="1556-276X-7-268-1"/></fig><p>However, a structural modification or surface functionalization may exist during CF<sub>4</sub> plasma treatment <abbrgrp><abbr bid="B14">14</abbr><abbr bid="B16">16</abbr><abbr bid="B17">17</abbr></abbrgrp>, even though they cannot be detected by normal optical microscopy. Previous researches revealed that the structural modification may induce a shift of the phonon frequencies in graphene sheets which can be optically probed by Raman spectroscopy <abbrgrp><abbr bid="B10">10</abbr><abbr bid="B21">21</abbr><abbr bid="B22">22</abbr></abbrgrp>. In order to guarantee the reproducibility of our experiment, the Raman signal was collected from the same spot with the assistance of optical microscopy. Figure <figr fid="F2">2</figr>a displays the Raman spectra of a graphene sheet before and after 2 s of CF<sub>4</sub> plasma treatment, and, as previously discussed, the optical microscopy shows no change. The lower panel in Figure <figr fid="F2">2</figr>a is the Raman spectrum of the sample before plasma treatment, where two strong peaks are noticeable. The G band located at approximately 1,580 cm<sup>&#8722;1</sup> is corresponding to optical E<sub>2g</sub> phonons at the Brillouin zone center. The sharp 2D band, which is the overtone of the D peak and the sum of two phonons with opposite momentum, appears at approximately 2,680 cm<sup>&#8722;1</sup><abbrgrp><abbr bid="B23">23</abbr></abbrgrp>. The intensity ratio between the 2D and G bands (<it>I</it><sub>2D</sub>/<it>I</it><sub>G</sub>&#8201;=&#8201;3.2) and the small full width at half maximum (FWHM) of the 2D band (26 cm<sup>&#8722;1</sup>) indicate that the graphene sheet measured contains only one layer (i.e., single-layer graphene) <abbrgrp><abbr bid="B14">14</abbr><abbr bid="B23">23</abbr></abbrgrp>. Besides, a weak D band is detected at approximately 1,350 cm<sup>&#8722;1</sup> and is connected to transverse optical phonons near the K point which require defects or lattice disorder (e.g., non-<it>sp</it><sup>2</sup> composition) for their activation via an inter-valley double-resonance Raman process <abbrgrp><abbr bid="B24">24</abbr></abbrgrp>. Therefore, the weak D band in the sample before plasma treatment suggests that the single-layer graphene sheet in our experiment is of high quality. It is interesting to note that obvious changes exist in the Raman spectra of the sample after 2 s of plasma treatment (as shown in the upper panel of Figure <figr fid="F2">2</figr>a): (1) The D peak at approximately 1,350 cm<sup>&#8722;1</sup> is remarkably intensified, suggesting an introduction of lattice disorder in the graphene sheet <abbrgrp><abbr bid="B24">24</abbr></abbrgrp>; (2) the G peak is broadened and a shoulder D' band arises, which originates from intra-valley resonant Raman scattering <abbrgrp><abbr bid="B24">24</abbr><abbr bid="B25">25</abbr></abbrgrp>; (3) a new band occurs at approximately 2,941 cm<sup>&#8722;1</sup> which is believed to be a combination of D and D' bands <abbrgrp><abbr bid="B17">17</abbr></abbrgrp> or G and D bands <abbrgrp><abbr bid="B26">26</abbr><abbr bid="B27">27</abbr></abbrgrp>. These features in the Raman spectra are reported in the graphene sheet treated with oxygen plasma <abbrgrp><abbr bid="B14">14</abbr></abbrgrp> and are the symbols of surface functionalization <abbrgrp><abbr bid="B14">14</abbr><abbr bid="B16">16</abbr></abbrgrp>. However, if the plasma treatment is elongated, the few-layer graphene sheet can be thinned and the evolution of Raman spectra becomes complicated. In our experiment, we found that single-layer graphene cannot survive the 5 s of CF<sub>4</sub> plasma treatment. Thus, the Raman measurement was carried out on a few-layer graphene sheet which contains more than one layer before plasma treatment. The obtained spectrum was plotted in the lower panel of Figure <figr fid="F2">2</figr>b. The FWHM of the 2D peak is broadened to 53 cm<sup>&#8722;1</sup> and the ratio of <it>I</it><sub>2D</sub>/<it>I</it><sub>G</sub> is 1.1, elucidating that this sheet consists of a few layers. After 5 s of plasma treatment, the Raman spectrum in the upper panel of Figure <figr fid="F2">2</figr>b shows that the FWHM of the 2D band is reduced to 32 cm<sup>&#8722;1</sup> and the ratio of <it>I</it><sub>2D</sub>/<it>I</it><sub>G</sub> is increased to 2.5. The spectral evolution is different from that in 2 s of plasma treatment and indicates that the few-layer graphene layer is thinned to one or two layers <abbrgrp><abbr bid="B14">14</abbr><abbr bid="B23">23</abbr><abbr bid="B28">28</abbr></abbrgrp>, as can be observed by optical microscopy.</p><fig id="F2"><title><p>Figure 2</p></title><caption><p>Raman spectra of graphene sheets before and after CF<sub>4</sub> plasma treatment</p></caption><text>
   <p><b>Raman spectra of graphene sheets before and after CF</b><sub><b>4</b></sub><b>plasma treatment.</b> ( <b>a</b>) Raman spectra of a single-layer graphene sheet before (lower panel) and after (upper panel) 2 s of CF<sub>4</sub> plasma treatment. (<b>b</b>) Raman spectra of a few-layer graphene sheet before (lower panel) and after (upper panel) 5 s of CF<sub>4</sub> plasma treatment. A thinning effect is obvious.</p>
</text><graphic file="1556-276X-7-268-2"/></fig><p>The evolution of the D band is the most important concern of many researches. The increase of its intensity is generally considered to be evidence that lattice disorder exists in the graphene sheet <abbrgrp><abbr bid="B14">14</abbr></abbrgrp>. An interesting phenomenon arises when we look close to the intensity of the D band in a few-layer graphene sheet subjected to CF<sub>4</sub> plasma for different times. The relative intensity of the D band (<it>I</it><sub>D</sub>/<it>I</it><sub>G</sub>) in the sample treated for 5 s is smaller than that for 2 s, implying that the disorder is even remarkable for shorter treatment. To understand this peculiar behavior, we must first make clear the mechanism of disorder production during plasma treatment. The disorder and the corresponding emergence of D and D' bands in our samples may arise from two processes: (1) physically, the order of pristine C atoms in graphene is disrupted by the CF<sub>4</sub> plasma and some of the C atoms may be sputtered out, which is named as ion bombardment effect <abbrgrp><abbr bid="B29">29</abbr><abbr bid="B30">30</abbr></abbrgrp>; and (2) chemically, covalent bonds of fluorine-related species to the graphene lattice form during plasma treatment, leading to corresponding surface modification and functionalization. Both processes contribute to the increase of the disorder. As for the chemical process, the bond energies need to be considered. It was disclosed that the C-C bond in graphene owns higher bond energy (607 kJ/mol) than the C-F bond (485 kJ/mol). Thus, the C-C bond can be hardly broken, and hence CF<sub>n</sub> (<it>n</it>&#8201;=&#8201;1 to 3) or F species were only adsorbed on the top graphene layer by formation of covalent bonds, which will become saturated because the number of available active C atoms decreases with time. It is worth noting that the formation of covalent bonds can be evidenced by the upshift of 2D peak from 2,685 to 2,691 cm<sup>&#8722;1</sup><abbrgrp><abbr bid="B31">31</abbr></abbrgrp>. One may infer that the few-layer graphene sheet subjected to CF<sub>4</sub> plasma treatment therefore possesses disorder features from both the physical and chemical interactions. The chemical interaction takes place only on the topmost layer, while the physical interaction accumulated with time and a long treatment can remove the top layer (i.e., thinning effect, as is reflected in the optical microscopy and Raman spectroscopy), exposing the beneath layer. The emergence and intensification of the D and D' bands in the sample subjected to 5 s of treatment thus originates from the disorder created in this new top layer and should increase gradually. Consequently, the disorder probed by Raman spectroscopy is even smaller in the sample treated by CF<sub>4</sub> plasma for 5 s than in the sample treated for 2 s (see Figure <figr fid="F2">2</figr>a,b).</p><p>Since the intensity of the D band reflects the disorder in the top layer of the few-layer graphene sheet, the relative intensity of the D band may give us a clue about the thickness of the sample. To prove this, we measured Raman spectra of plasma-thinned few-layer graphene sheets with different thicknesses. The three spectra in Figure <figr fid="F3">3</figr> were respectively collected from three spots labeled in the inset optical image. The color contrast among the three spots indicates that their thicknesses are rather different and increase from A to C, while none of them consist a single layer. One can see that the <it>I</it><sub>D</sub>/<it>I</it><sub>G</sub> is obviously relevant to the thickness of few-layer graphene sheet. For a thicker sheet, the weight of the topmost layer with remarkable disorder is smaller, and thus <it>I</it><sub>D</sub>/<it>I</it><sub>G</sub> should be correspondingly smaller. Therefore, the ration of <it>I</it><sub>D</sub>/<it>I</it><sub>G</sub> proves to be an effective Raman factor to compare the thicknesses of plasma-treated few-layer graphene sheets. Therefore, the <it>I</it><sub>D</sub>/<it>I</it><sub>G</sub> proves to be another effective factor to estimate the thicknesses of plasma-treated few-layer graphene sheets other than <it>I</it><sub>2D</sub>/<it>I</it><sub>G</sub>.</p><fig id="F3"><title><p>Figure 3</p></title><caption><p>Raman spectra of plasma-thinned few-layer graphene sheets with different thicknesses</p></caption><text>
   <p><b>Raman spectra of plasma-thinned few-layer graphene sheets with different thicknesses.</b> The three spectra were collected from three different spots respectively. The positions of the three spots are labeled in the inset optical image of a sample subjected to 5 s of CF<sub>4</sub> plasma treatment. Scale bar, 5 &#956;m.</p>
</text><graphic file="1556-276X-7-268-3"/></fig><p>Although a detailed experimental investigation is so far difficult, we managed to simulate the ion adsorption during the plasma treatment by employing density functional theory and the generalized gradient approximation of the exchange correlation functional with ultrasoft pseudopotentials <abbrgrp><abbr bid="B32">32</abbr><abbr bid="B33">33</abbr></abbrgrp> to reveal the disorder produced by fluorine functionalization and the underneath mechanism. A high cutoff energy of 800 eV and a k-point sampling with an 8&#8201;&#215;&#8201;8&#8201;&#215;&#8201;1 mesh are employed to achieve high accuracy in the calculations. Structural optimization is carried out on all systems until the residual forces are converged to 0.003 eV/&#197;. In order to avoid any drawback of the periodic boundary conditions, an over 20-&#197;-thick vacuum layer is included. We estimate that in the plasma, there are F, CF, CF<sub>2</sub>, and CF<sub>3</sub> ions. The energy barrier is calculated by the climbing-image nudged elastic band method <abbrgrp><abbr bid="B34">34</abbr></abbrgrp>, which enables us to find the minimum energy path between the given initial and final states of a transition. By calculating the structures of ions on the graphene sheet (see Figure <figr fid="F4">4</figr>), we find that CF and CF<sub>2</sub> can hardly be adsorbed, forming covalent bonds, while F and CF<sub>3</sub> ions can be adsorbed. The corresponding adsorption energies are 2.2 and 0.4 eV, respectively. The different absorption behaviors can be understood as follows: Because of the <it>sp</it><sup>3</sup> hybridization nature for C atoms in CF, CF<sub>2</sub>, and CF<sub>3</sub>, there are two un-bonded electrons in CF<sub>2</sub> and only one in CF<sub>3</sub>, while the C atoms in graphene are <it>sp</it><sup>2</sup> hybridized and there is only one un-bonded electron for each C atom. Thus, it is easier for CF<sub>3</sub> to be adsorbed on top of the C atom via a covalent bond (see Figure <figr fid="F4">4</figr>d). However, it is difficult for CF<sub>2</sub> to form two covalent bonds with two nearby C atoms in the graphene sheet. We checked different absorption configurations by structure relaxation and find that for all cases, the CF<sub>2</sub> will leave away from the graphene, as shown in Figure <figr fid="F4">4</figr>c. For different artificial absorption configurations, there would be large strain/stress around the absorption site, indicating that the difficulty of forming two covalent bonds between CF<sub>2</sub> and graphene is mainly due to the large strain/stress. A similar situation is also noted in the CF case. In addition, the calculated energy barrier for the CF<sub>3</sub> adsorbed on the graphene surface is only 0.04 eV, which demonstrates that the CF<sub>3</sub> adsorption is energetically favorable. Due to the existence of covalent bonds between F/CF<sub>3</sub> and C on the graphene sheet, there is distortion around the adsorption position, which gives birth to the appearance of D and D' bands in the Raman spectra of the plasma-treated samples.</p><fig id="F4"><title><p>Figure 4</p></title><caption><p>Atomic structures of (a) F, (b) CF, (c) CF<sub>2</sub>, and (d) CF<sub>3</sub> ions on graphene sheets</p></caption><text>
   <p><b>Atomic structures of (a) F, (b) CF, (c) CF</b><sub><b>2</b></sub><b>, and (d) CF</b><sub><b>3</b></sub><b>ions on graphene sheets.</b> C and F atoms are indicated by yellow and gray balls, respectively. Theoretical simulation indicates that CF and CF<sub>2</sub> can hardly be adsorbed.</p>
</text><graphic file="1556-276X-7-268-4"/></fig><p>Here, we would like to briefly discuss the potential applications of this fluorine functionalization. We calculated the band structures and densities of states of graphene sheets fluorinated with F and CF<sub>3</sub> with coverage of 5.6 %, as shown in Figure <figr fid="F5">5</figr>. Due to the Brillouin zone folding <abbrgrp><abbr bid="B20">20</abbr></abbrgrp>, the Dirac point is located at the &#915; point in the electronic band structure. It is worth noting that due to the symmetry breaking, the Dirac cone splits with a small gap (approximately 0.1 eV). Moreover, the Dirac cone shifts to 0.5/0.1 eV above the Fermi level after F/CF<sub>3</sub> functionalization due to electron transfer from the C atoms in the graphene sheet to F/CF<sub>3</sub>, indicating that the graphene sheets can be p-type doped by fluorinating with both F and CF<sub>3</sub>. However, it can be found that F is much more efficient to introduce p-type doping in graphene sheets than CF<sub>3</sub>. In addition, the adsorption of the F/CF<sub>3</sub> on graphene sheets also introduces pinning states around the Fermi levels, as demonstrated in Figure <figr fid="F5">5</figr>. Detailed experimental verification is required in the future, while the presented calculation has already proven the application potentials of graphene sheets functionalized by CF<sub>4</sub> plasma in micro/nanoelectronics.</p><fig id="F5"><title><p>Figure 5</p></title><caption><p>Band structures and densities of states of F- and CF<sub>3</sub>-functionalized graphene sheets</p></caption><text>
   <p><b>Band structures and densities of states of F- and CF</b><sub><b>3</b></sub><b>-functionalized graphene sheets.</b> ( <b>a</b>) Band structure and ( <b>b</b>) density of states of a F-functionalized graphene sheet. ( <b>c</b>) Band structure and ( <b>d</b>) density of states of a CF<sub>3</sub>-functionalized graphene sheet. The Fermi level is set at 0 eV.</p>
</text><graphic file="1556-276X-7-268-5"/></fig></sec><sec><st><p>Conclusion</p></st><p>In conclusion, few-layer graphene sheets were prepared by mechanical exfoliation and the structural evolution during CF<sub>4</sub> plasma treatment was studied in detail by optical microscopy and Raman spectroscopy. The experimental results indicate a thickness reduction under prolonged plasma treatment while short treatment leads only to fluorine functionalization on the surface layer. The combination of both physical and chemical reactions in the plasma treatment leads to structural modifications which can be well probed in Raman spectra. Theoretical simulation suggests a F/CF<sub>3</sub> functionalization by formation of covalent bonds and also predicts a corresponding p-type doping and Dirac cone opening after the F/CF<sub>3</sub> adsorption. Although further characterizations are needed to evaluate the electronic properties of treated samples, the current work of thinning and functionalizing few-layer graphene sheets by CF<sub>4</sub> plasma under control represents an integrative pathway to industrial fabrication of graphene-based micro/nanodevices.</p></sec><sec><st><p>Competing interests</p></st><p>The authors declare that they have no competing interests.</p></sec><sec><st><p>Authors' contributions</p></st><p>CS, GH, FD, and YM designed the study. CS performed the experiments with help from RC. YC and US carried out the theoretical study. CS, GH, YC, and YM contributed in drafting the manuscript. All the authors took part in the discussion of the results, and edited and approved the manuscript.</p></sec></bdy><bm><ack><sec><st><p>Acknowledgments</p></st><p>This work is supported by the Natural Science Foundation of China (Nos. 61008029 and 51102049), Program for New Century Excellent Talents in University (No. NCET-10-0345), and Shanghai Pujiang Program (No. 11PJ1400900).</p></sec></ack><refgrp><bibl id="B1"><title><p>Experimental observation of the quantum Hall effect and Berry's phase in graphene</p></title><aug><au><snm>Zhang</snm><fnm>YB</fnm></au><au><snm>Tan</snm><fnm>YW</fnm></au><au><snm>Stormer</snm><fnm>HL</fnm></au><au><snm>Kim</snm><fnm>P</fnm></au></aug><source>Nature</source><pubdate>2005</pubdate><volume>438</volume><issue>7065</issue><fpage>201</fpage><lpage>204</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1038/nature04235</pubid><pubid idtype="pmpid" link="fulltext">16281031</pubid></pubidlist></xrefbib></bibl><bibl id="B2"><title><p>The transport properties of graphene: an introduction</p></title><aug><au><snm>Peres</snm><fnm>NM</fnm></au></aug><source>Rev Mod Phys</source><pubdate>2010</pubdate><volume>82</volume><issue>3</issue><fpage>2673</fpage><lpage>2700</lpage><xrefbib><pubid idtype="doi">10.1103/RevModPhys.82.2673</pubid></xrefbib></bibl><bibl id="B3"><title><p>Andreev reflection and Klein tunneling in graphene</p></title><aug><au><snm>Beenakker</snm><fnm>CWJ</fnm></au></aug><source>Rev Mod Phys</source><pubdate>2008</pubdate><volume>80</volume><issue>4</issue><fpage>1337</fpage><lpage>1354</lpage><xrefbib><pubid idtype="doi">10.1103/RevModPhys.80.1337</pubid></xrefbib></bibl><bibl id="B4"><title><p>Electric field effect in atomically thin carbon films</p></title><aug><au><snm>Novoselov</snm><fnm>KS</fnm></au><au><snm>Geim</snm><fnm>AK</fnm></au><au><snm>Morozov</snm><fnm>SV</fnm></au><au><snm>Jiang</snm><fnm>D</fnm></au><au><snm>Zhang</snm><fnm>Y</fnm></au><au><snm>Dubonos</snm><fnm>SV</fnm></au><au><snm>Grigorieva</snm><fnm>IV</fnm></au><au><snm>Firsov</snm><fnm>AA</fnm></au></aug><source>Science</source><pubdate>2004</pubdate><volume>306</volume><issue>5696</issue><fpage>666</fpage><lpage>669</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1126/science.1102896</pubid><pubid idtype="pmpid" link="fulltext">15499015</pubid></pubidlist></xrefbib></bibl><bibl id="B5"><title><p>Two-dimensional gas of massless Dirac fermions in graphene</p></title><aug><au><snm>Novoselov</snm><fnm>KS</fnm></au><au><snm>Geim</snm><fnm>AK</fnm></au><au><snm>Morozov</snm><fnm>SV</fnm></au><au><snm>Jiang</snm><fnm>D</fnm></au><au><snm>Katsnelson</snm><fnm>MI</fnm></au><au><snm>Grigorieva</snm><fnm>IV</fnm></au><au><snm>Dubonos</snm><fnm>SV</fnm></au><au><snm>Firsov</snm><fnm>AA</fnm></au></aug><source>Nature</source><pubdate>2005</pubdate><volume>438</volume><issue>7065</issue><fpage>197</fpage><lpage>200</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1038/nature04233</pubid><pubid idtype="pmpid" link="fulltext">16281030</pubid></pubidlist></xrefbib></bibl><bibl id="B6"><title><p>Graphene cracking bilayers</p></title><aug><au><snm>Novoselov</snm><fnm>KS</fnm></au></aug><source>Nat Phys</source><pubdate>2009</pubdate><volume>5</volume><issue>12</issue><fpage>862</fpage><lpage>863</lpage><xrefbib><pubid idtype="doi">10.1038/nphys1471</pubid></xrefbib></bibl><bibl id="B7"><title><p>Coulomb interaction, ripples, and the minimal conductivity of graphene</p></title><aug><au><snm>Herbut</snm><fnm>IF</fnm></au><au><snm>Juricic</snm><fnm>V</fnm></au><au><snm>Vafek</snm><fnm>O</fnm></au></aug><source>Phys Rev Lett</source><pubdate>2008</pubdate><volume>100</volume><issue>4</issue><fpage>046403</fpage><xrefbib><pubid idtype="pmpid" link="fulltext">18352311</pubid></xrefbib></bibl><bibl id="B8"><title><p>Coulomb interaction at the metal-insulator critical point in graphene</p></title><aug><au><snm>Juricic</snm><fnm>V</fnm></au><au><snm>Herbut</snm><fnm>IF</fnm></au><au><snm>Semenoff</snm><fnm>GW</fnm></au></aug><source>Phys Rev B</source><pubdate>2009</pubdate><volume>80</volume><issue>8</issue><fpage>081405</fpage></bibl><bibl id="B9"><title><p>Gruneisen parameter of the G mode of strained monolayer graphene</p></title><aug><au><snm>Cheng</snm><fnm>YC</fnm></au><au><snm>Zhu</snm><fnm>ZY</fnm></au><au><snm>Huang</snm><fnm>GS</fnm></au><au><snm>Schwingenschl&#246;gl</snm><fnm>U</fnm></au></aug><source>Phys Rev B</source><pubdate>2011</pubdate><volume>83</volume><issue>11</issue><fpage>115449</fpage></bibl><bibl id="B10"><title><p>Stretchable graphene: a close look at fundamental parameters through biaxial straining</p></title><aug><au><snm>Ding</snm><fnm>F</fnm></au><au><snm>Ji</snm><fnm>HX</fnm></au><au><snm>Chen</snm><fnm>YH</fnm></au><au><snm>Herklotz</snm><fnm>A</fnm></au><au><snm>D&#246;rr</snm><fnm>K</fnm></au><au><snm>Mei</snm><fnm>YF</fnm></au><au><snm>Rastelli</snm><fnm>A</fnm></au><au><snm>Schmidt</snm><fnm>OG</fnm></au></aug><source>Nano Lett</source><pubdate>2010</pubdate><volume>10</volume><issue>9</issue><fpage>3453</fpage><lpage>3458</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1021/nl101533x</pubid><pubid idtype="pmpid" link="fulltext">20695450</pubid></pubidlist></xrefbib></bibl><bibl id="B11"><title><p>Carbon-based electronics</p></title><aug><au><snm>Avouris</snm><fnm>P</fnm></au><au><snm>Chen</snm><fnm>Z</fnm></au><au><snm>Perebeinos</snm><fnm>V</fnm></au></aug><source>Nat Nanotechnol</source><pubdate>2007</pubdate><volume>2</volume><issue>10</issue><fpage>605</fpage><lpage>615</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1038/nnano.2007.300</pubid><pubid idtype="pmpid">18654384</pubid></pubidlist></xrefbib></bibl><bibl id="B12"><title><p>Spin field effect transistor with a graphene channel</p></title><aug><au><snm>Semenov</snm><fnm>YG</fnm></au><au><snm>Kim</snm><fnm>KW</fnm></au><au><snm>Zavada</snm><fnm>JM</fnm></au></aug><source>Appl Phys Lett</source><pubdate>2007</pubdate><volume>91</volume><issue>15</issue><fpage>153105</fpage><xrefbib><pubid idtype="doi">10.1063/1.2798596</pubid></xrefbib></bibl><bibl id="B13"><title><p>Electrical control of the chemical bonding of fluorine on graphene</p></title><aug><au><snm>Sofo</snm><fnm>JO</fnm></au><au><snm>Suarez</snm><fnm>AM</fnm></au><au><snm>Usaj</snm><fnm>G</fnm></au><au><snm>Cornaglia</snm><fnm>PS</fnm></au><au><snm>Hernandez-Nieves</snm><fnm>AD</fnm></au><au><snm>Balseiro</snm><fnm>CA</fnm></au></aug><source>Phys Rev B</source><pubdate>2011</pubdate><volume>83</volume><issue>8</issue><fpage>081411</fpage></bibl><bibl id="B14"><title><p>Bandgap opening in oxygen plasma-treated graphene</p></title><aug><au><snm>Nourbakhsh</snm><fnm>A</fnm></au><au><snm>Cantoro</snm><fnm>M</fnm></au><au><snm>Vosch</snm><fnm>T</fnm></au><au><snm>Pourtois</snm><fnm>G</fnm></au><au><snm>Clemente</snm><fnm>F</fnm></au><au><snm>van der Veen</snm><fnm>MH</fnm></au><au><snm>Hofkens</snm><fnm>J</fnm></au><au><snm>Heyns</snm><fnm>MM</fnm></au><au><snm>De Gendt</snm><fnm>S</fnm></au><au><snm>Sels</snm><fnm>BF</fnm></au></aug><source>Nanotechnology</source><pubdate>2010</pubdate><volume>21</volume><issue>43</issue><fpage>435203</fpage><xrefbib><pubidlist><pubid idtype="doi">10.1088/0957-4484/21/43/435203</pubid><pubid idtype="pmpid" link="fulltext">20890016</pubid></pubidlist></xrefbib></bibl><bibl id="B15"><title><p>Surface-energy engineering of graphene</p></title><aug><au><snm>Shin</snm><fnm>YJ</fnm></au><au><snm>Wang</snm><fnm>Y</fnm></au><au><snm>Huang</snm><fnm>H</fnm></au><au><snm>Kalon</snm><fnm>G</fnm></au><au><snm>Wee</snm><fnm>ATS</fnm></au><au><snm>Shen</snm><fnm>Z</fnm></au><au><snm>Bhatia</snm><fnm>CS</fnm></au><au><snm>Yang</snm><fnm>H</fnm></au></aug><source>Langmuir</source><pubdate>2010</pubdate><volume>26</volume><issue>6</issue><fpage>3798</fpage><lpage>3802</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1021/la100231u</pubid><pubid idtype="pmpid" link="fulltext">20158275</pubid></pubidlist></xrefbib></bibl><bibl id="B16"><title><p>The functionalization of graphene using electron-beam generated plasmas</p></title><aug><au><snm>Baraket</snm><fnm>M</fnm></au><au><snm>Walton</snm><fnm>SG</fnm></au><au><snm>Lock</snm><fnm>EH</fnm></au><au><snm>Robinson</snm><fnm>JT</fnm></au><au><snm>Perkins</snm><fnm>FK</fnm></au></aug><source>Appl Phys Lett</source><pubdate>2010</pubdate><volume>96</volume><issue>23</issue><fpage>231501</fpage><xrefbib><pubid idtype="doi">10.1063/1.3436556</pubid></xrefbib></bibl><bibl id="B17"><title><p>Making graphene luminescent by oxygen plasma treatment</p></title><aug><au><snm>Gokus</snm><fnm>T</fnm></au><au><snm>Nair</snm><fnm>RR</fnm></au><au><snm>Bonetti</snm><fnm>A</fnm></au><au><snm>Boehmler</snm><fnm>M</fnm></au><au><snm>Lombardo</snm><fnm>A</fnm></au><au><snm>Novoselov</snm><fnm>KS</fnm></au><au><snm>Geim</snm><fnm>AK</fnm></au><au><snm>Ferrari</snm><fnm>AC</fnm></au><au><snm>Hartschuh</snm><fnm>A</fnm></au></aug><source>ACS Nano</source><pubdate>2009</pubdate><volume>3</volume><issue>12</issue><fpage>3963</fpage><lpage>3968</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1021/nn9012753</pubid><pubid idtype="pmpid" link="fulltext">19925014</pubid></pubidlist></xrefbib></bibl><bibl id="B18"><title><p>Surface analysis and bioreactions of F and Si containing a-C:H</p></title><aug><au><snm>Hauert</snm><fnm>R</fnm></au><au><snm>Muller</snm><fnm>U</fnm></au><au><snm>Francz</snm><fnm>G</fnm></au><au><snm>Birchler</snm><fnm>F</fnm></au><au><snm>Schroeder</snm><fnm>A</fnm></au><au><snm>Mayer</snm><fnm>J</fnm></au><au><snm>Wintermantel</snm><fnm>E</fnm></au></aug><source>Thin Solid Films</source><pubdate>1997</pubdate><volume>308</volume><fpage>191</fpage><lpage>194</lpage></bibl><bibl id="B19"><title><p>Highly p-doped epitaxial graphene obtained by fluorine intercalation</p></title><aug><au><snm>Walter</snm><fnm>AL</fnm></au><au><snm>Jeon</snm><fnm>K-J</fnm></au><au><snm>Bostwick</snm><fnm>A</fnm></au><au><snm>Speck</snm><fnm>F</fnm></au><au><snm>Ostler</snm><fnm>M</fnm></au><au><snm>Seyller</snm><fnm>T</fnm></au><au><snm>Moreschini</snm><fnm>L</fnm></au><au><snm>Kim</snm><fnm>YS</fnm></au><au><snm>Chang</snm><fnm>YJ</fnm></au><au><snm>Horn</snm><fnm>K</fnm></au><au><snm>Rotenberg</snm><fnm>E</fnm></au></aug><source>Appl Phys Lett</source><pubdate>2011</pubdate><volume>98</volume><issue>18</issue><fpage>184102</fpage><xrefbib><pubid idtype="doi">10.1063/1.3586256</pubid></xrefbib></bibl><bibl id="B20"><title><p>Origin of the high p-doping in F intercalated graphene on SiC</p></title><aug><au><snm>Cheng</snm><fnm>YC</fnm></au><au><snm>Kaloni</snm><fnm>TP</fnm></au><au><snm>Huang</snm><fnm>GS</fnm></au><au><snm>Schwingenschl&#246;gl</snm><fnm>U</fnm></au></aug><source>Appl Phys Lett</source><pubdate>2011</pubdate><volume>99</volume><issue>5</issue><fpage>053117</fpage><xrefbib><pubid idtype="doi">10.1063/1.3623484</pubid></xrefbib></bibl><bibl id="B21"><title><p>Compression behavior of single-layer graphenes</p></title><aug><au><snm>Frank</snm><fnm>O</fnm></au><au><snm>Tsoukleri</snm><fnm>G</fnm></au><au><snm>Parthenios</snm><fnm>J</fnm></au><au><snm>Papagelis</snm><fnm>K</fnm></au><au><snm>Riaz</snm><fnm>I</fnm></au><au><snm>Jalil</snm><fnm>R</fnm></au><au><snm>Novoselov</snm><fnm>KS</fnm></au><au><snm>Galiotis</snm><fnm>C</fnm></au></aug><source>ACS Nano</source><pubdate>2010</pubdate><volume>4</volume><issue>6</issue><fpage>3131</fpage><lpage>3138</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1021/nn100454w</pubid><pubid idtype="pmpid" link="fulltext">20496881</pubid></pubidlist></xrefbib></bibl><bibl id="B22"><title><p>Phonon softening and crystallographic orientation of strained graphene studied by Raman spectroscopy</p></title><aug><au><snm>Huang</snm><fnm>M</fnm></au><au><snm>Yan</snm><fnm>H</fnm></au><au><snm>Chen</snm><fnm>C</fnm></au><au><snm>Song</snm><fnm>D</fnm></au><au><snm>Heinz</snm><fnm>TF</fnm></au><au><snm>Hone</snm><fnm>J</fnm></au></aug><source>Proc Natl Acad Sci USA</source><pubdate>2009</pubdate><volume>106</volume><issue>18</issue><fpage>7304</fpage><lpage>7308</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1073/pnas.0811754106</pubid><pubid idtype="pmcid">2678610</pubid><pubid idtype="pmpid" link="fulltext">19380746</pubid></pubidlist></xrefbib></bibl><bibl id="B23"><title><p>Raman spectrum of graphene and graphene layers</p></title><aug><au><snm>Ferrari</snm><fnm>AC</fnm></au><au><snm>Meyer</snm><fnm>JC</fnm></au><au><snm>Scardaci</snm><fnm>V</fnm></au><au><snm>Casiraghi</snm><fnm>C</fnm></au><au><snm>Lazzeri</snm><fnm>M</fnm></au><au><snm>Mauri</snm><fnm>F</fnm></au><au><snm>Piscanec</snm><fnm>S</fnm></au><au><snm>Jiang</snm><fnm>D</fnm></au><au><snm>Novoselov</snm><fnm>KS</fnm></au><au><snm>Roth</snm><fnm>S</fnm></au><au><snm>Geim</snm><fnm>AK</fnm></au></aug><source>Phys Rev Lett</source><pubdate>2006</pubdate><volume>97</volume><issue>18</issue><fpage>187401</fpage><xrefbib><pubid idtype="pmpid" link="fulltext">17155573</pubid></xrefbib></bibl><bibl id="B24"><title><p>Studying disorder in graphite-based systems by Raman spectroscopy</p></title><aug><au><snm>Pimenta</snm><fnm>MA</fnm></au><au><snm>Dresselhaus</snm><fnm>G</fnm></au><au><snm>Dresselhaus</snm><fnm>MS</fnm></au><au><snm>Cancado</snm><fnm>LG</fnm></au><au><snm>Jorio</snm><fnm>A</fnm></au><au><snm>Saito</snm><fnm>R</fnm></au></aug><source>Phys Chem Chem Phys</source><pubdate>2007</pubdate><volume>9</volume><issue>11</issue><fpage>1276</fpage><lpage>1291</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1039/b613962k</pubid><pubid idtype="pmpid" link="fulltext">17347700</pubid></pubidlist></xrefbib></bibl><bibl id="B25"><title><p>Influence of the atomic structure on the Raman spectra of graphite edges</p></title><aug><au><snm>Cancado</snm><fnm>LG</fnm></au><au><snm>Pimenta</snm><fnm>MA</fnm></au><au><snm>Neves</snm><fnm>BRA</fnm></au><au><snm>Dantas</snm><fnm>MSS</fnm></au><au><snm>Jorio</snm><fnm>A</fnm></au></aug><source>Phys Rev Lett</source><pubdate>2004</pubdate><volume>93</volume><issue>24</issue><fpage>247401</fpage><xrefbib><pubid idtype="pmpid" link="fulltext">15697858</pubid></xrefbib></bibl><bibl id="B26"><title><p>Optoelectronic properties of graphene thin films prepared by thermal reduction of graphene oxide</p></title><aug><au><snm>Cuong</snm><fnm>TV</fnm></au><au><snm>Pham</snm><fnm>VH</fnm></au><au><snm>Tran</snm><fnm>QT</fnm></au><au><snm>Chung</snm><fnm>JS</fnm></au><au><snm>Shin</snm><fnm>EW</fnm></au><au><snm>Kim</snm><fnm>JS</fnm></au><au><snm>Kim</snm><fnm>EJ</fnm></au></aug><source>Mater Lett</source><pubdate>2010</pubdate><volume>64</volume><issue>6</issue><fpage>765</fpage><lpage>767</lpage><xrefbib><pubid idtype="doi">10.1016/j.matlet.2010.01.009</pubid></xrefbib></bibl><bibl id="B27"><title><p>Photoluminescence and Raman studies of graphene thin films prepared by reduction of graphene oxide</p></title><aug><au><snm>Cuong</snm><fnm>TV</fnm></au><au><snm>Pham</snm><fnm>VH</fnm></au><au><snm>Tran</snm><fnm>QT</fnm></au><au><snm>Hahn</snm><fnm>SH</fnm></au><au><snm>Chung</snm><fnm>JS</fnm></au><au><snm>Shin</snm><fnm>EW</fnm></au><au><snm>Kim</snm><fnm>EJ</fnm></au></aug><source>Mater Lett</source><pubdate>2010</pubdate><volume>64</volume><issue>3</issue><fpage>399</fpage><lpage>401</lpage><xrefbib><pubid idtype="doi">10.1016/j.matlet.2009.11.029</pubid></xrefbib></bibl><bibl id="B28"><title><p>Ultraviolet Raman microscopy of single and multilayer graphene</p></title><aug><au><snm>Calizo</snm><fnm>I</fnm></au><au><snm>Bejenari</snm><fnm>I</fnm></au><au><snm>Rahman</snm><fnm>M</fnm></au><au><snm>Liu</snm><fnm>G</fnm></au><au><snm>Balandin</snm><fnm>AA</fnm></au></aug><source>J Appl Phys</source><pubdate>2009</pubdate><volume>106</volume><issue>4</issue><fpage>043509</fpage><xrefbib><pubid idtype="doi">10.1063/1.3197065</pubid></xrefbib></bibl><bibl id="B29"><title><p>Photoluminescence and band gap modulation in graphene oxide</p></title><aug><au><snm>Luo</snm><fnm>Z</fnm></au><au><snm>Vora</snm><fnm>PM</fnm></au><au><snm>Mele</snm><fnm>EJ</fnm></au><au><snm>Johnson</snm><fnm>ATC</fnm></au><au><snm>Kikkawa</snm><fnm>JM</fnm></au></aug><source>Appl Phys Lett</source><pubdate>2009</pubdate><volume>94</volume><issue>11</issue><fpage>111909</fpage><xrefbib><pubid idtype="doi">10.1063/1.3098358</pubid></xrefbib></bibl><bibl id="B30"><title><p>Photoelectron spectroscopy studies of plasma-fluorinated epitaxial graphene</p></title><aug><au><snm>Sherpa</snm><fnm>SD</fnm></au><au><snm>Paniagua</snm><fnm>SA</fnm></au><au><snm>Levitin</snm><fnm>G</fnm></au><au><snm>Marder</snm><fnm>SR</fnm></au><au><snm>Williams</snm><fnm>MD</fnm></au><au><snm>Hess</snm><fnm>DW</fnm></au></aug><source>J Vac Sci Technol B</source><pubdate>2012</pubdate><volume>30</volume><issue>3</issue><fpage>03D102</fpage><xrefbib><pubid idtype="doi">10.1116/1.3688760</pubid></xrefbib></bibl><bibl id="B31"><title><p>Raman spectroscopic characterization of graphene</p></title><aug><au><snm>Tang</snm><fnm>B</fnm></au><au><snm>Hu</snm><fnm>GX</fnm></au><au><snm>Gao</snm><fnm>HY</fnm></au></aug><source>Appl Spectrosc Rev</source><pubdate>2010</pubdate><volume>45</volume><issue>5</issue><fpage>369</fpage><lpage>407</lpage><xrefbib><pubid idtype="doi">10.1080/05704928.2010.483886</pubid></xrefbib></bibl><bibl id="B32"><title><p>Soft self-consistent pseduopotentials in a generalized eigenvalue formalism</p></title><aug><au><snm>Vanderbilt</snm><fnm>D</fnm></au></aug><source>Phys Rev B</source><pubdate>1990</pubdate><volume>41</volume><issue>11</issue><fpage>7892</fpage><lpage>7895</lpage><xrefbib><pubid idtype="doi">10.1103/PhysRevB.41.7892</pubid></xrefbib></bibl><bibl id="B33"><title><p>Quantum espresso: a modular and open-source software project for quantum simulations of materials</p></title><aug><au><snm>Giannozzi</snm><fnm>P</fnm></au><au><snm>Baroni</snm><fnm>S</fnm></au><au><snm>Bonini</snm><fnm>N</fnm></au><au><snm>Calandra</snm><fnm>M</fnm></au><au><snm>Car</snm><fnm>R</fnm></au><au><snm>Cavazzoni</snm><fnm>C</fnm></au><au><snm>Ceresoli</snm><fnm>D</fnm></au><au><snm>Chiarotti</snm><fnm>GL</fnm></au><au><snm>Cococcioni</snm><fnm>M</fnm></au><au><snm>Dabo</snm><fnm>I</fnm></au><au><snm>Dal Corso</snm><fnm>A</fnm></au><au><snm>Fabris</snm><fnm>S</fnm></au><au><snm>Fratesi</snm><fnm>G</fnm></au><au><snm>De Gironcoli</snm><fnm>S</fnm></au><au><snm>Gebauer</snm><fnm>R</fnm></au><au><snm>Gerstmann</snm><fnm>U</fnm></au><au><snm>Gougoussis</snm><fnm>C</fnm></au><au><snm>Kokalj</snm><fnm>A</fnm></au><au><snm>Lazzeri</snm><fnm>M</fnm></au><au><snm>Martin-Samos</snm><fnm>L</fnm></au><au><snm>Marzari</snm><fnm>N</fnm></au><au><snm>Mauri</snm><fnm>F</fnm></au><au><snm>Mazzarello</snm><fnm>R</fnm></au><au><snm>Paolini</snm><fnm>S</fnm></au><au><snm>Pasquarello</snm><fnm>A</fnm></au><au><snm>Paulatto</snm><fnm>L</fnm></au><au><snm>Sbraccia</snm><fnm>C</fnm></au><au><snm>Scandolo</snm><fnm>S</fnm></au><au><snm>Sclauzero</snm><fnm>G</fnm></au><au><snm>Seitsonen</snm><fnm>AP</fnm></au><au><snm>Smogunov</snm><fnm>A</fnm></au><au><snm>Umari</snm><fnm>P</fnm></au><au><snm>Wentzcovitch</snm><fnm>RM</fnm></au></aug><source>J Phys: Condens Mater</source><pubdate>2009</pubdate><volume>21</volume><issue>39</issue><fpage>395502</fpage><xrefbib><pubid idtype="doi">10.1088/0953-8984/21/39/395502</pubid></xrefbib></bibl><bibl id="B34"><title><p>Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points</p></title><aug><au><snm>Henkelman</snm><fnm>G</fnm></au><au><snm>Jonsson</snm><fnm>H</fnm></au></aug><source>J Chem Phys</source><pubdate>2000</pubdate><volume>113</volume><issue>22</issue><fpage>9978</fpage><lpage>9985</lpage><xrefbib><pubid idtype="doi">10.1063/1.1323224</pubid></xrefbib></bibl></refgrp></bm></art>