Electrostatic charge distribution in armchair and zigzag carbon nanotubes: a numerical comparison of CNT charge models

In engineering—especially in mechanics, lightweight engineering, aerospace engineering, electrical engineering as well as bioengineering—there is a desire of developing materials enabling an excellent performance with respect to mechanical, thermal, and electrical properties. One of the most promising materials are carbon nanotubes (CNTs), as they show excellent mechanical and electrical properties. To improve the understanding of the electrical behavior, i.e., the charge distribution in single-walled carbon nanotubes, both open ended armchair and zigzag types are investigated. In the present research, three different modeling approaches, a classic electrostatic model, the model by Li and Chou and the model by Mayer, are analyzed and compared with respect to their further applicability. In the numerical investigations, different test cases are performed: (i) the carbon nanotubes are charged with an overall charge, (ii) the CNTs are exposed to an external electric field, and (iii) the test cases (i) and (ii) are combined. Furthermore, the influence of different geometric parameters is investigated. It is shown that the charge applied to the CNTs distributes over the whole CNT having maxima at the ends of the tubes. These maxima can be influenced by both, the geometric parameters and the electric field strength.


Introduction
One of the highest driving factors of mankind is to be better, faster, and stronger. The same mindset applies to the materials humans are using and developing constantly. For that reason, a new structure formed out of carbon atoms was first discovered in 1952 by Radushkevich and Lukyanovich [1]. They discovered that carbon atoms not only can shape diamonds, graphene or fullerenes but also a tube, forming a carbon nanotube (CNT). Iijima [2] described later how this structure can be synthesized. This was the starting point for several investigations about this new and promising material. The different methods of fabrication were reported by Lu et al. [3], and the outstanding mechanical, thermal, and electrical properties are referred in [4][5][6].
As a result of different investigations about the mechanical and thermal properties, various applications have been proposed. For example, carbon nanotubes are applied in chemical sensors [7,8] or for the detection of gases [9,10].
Despite of the numerous investigations, a large research effort is still necessary. Especially the reaction of a carbon nanotube to an applied voltage or an applied electric field is still of interest. Therefore, different K. Lönnecke · O. Eberhardt · T. Wallmersperger (B) TU Dresden, Institut für Festkörpermechanik, Professur für Mechanik multifunktionaler Strukturen, George-Bähr-Str. 3c, 01069 Dresden, Germany E-mail: thomas.wallmersperger@tu-dresden.de models for simulating the charge distribution of CNTs have been proposed [21][22][23][24]. In the present research, the electrostatic behavior of carbon nanotubes is investigated by applying three different modeling approaches.
In order to get a better understanding of the electric properties of carbon nanotubes, the modeling approaches are used to obtain results of the behavior of the nanotube toward (i) increasing overall charges and (ii) different applied electric fields. The motivation of the present investigation is to understand the developed modeling approaches, to compare them, and to get a first impression on the deformation of the CNTs, due to electromechanical coupling. Furthermore, it is interesting to understand the differences and decide which model can be preferred for which specific application. Additionally, the comparison will show if the modeling approaches are equally suitable for armchair nanotubes and zigzag CNTs. In conclusion, the presented electrical models are the first step toward the electromechanical coupling-i.e., the deformation of CNTs [25] due to an applied charge or an applied electric field-of carbon nanotubes.

Fundamentals and modeling approaches for the charge distribution in carbon nanotubes
In order to understand (i) the mechanical structure of the carbon nanotubes and (ii) the computation of the charge distribution, a general introduction to CNTs is given. A carbon nanotube is formed by carbon atoms which are arranged in a very specific way. To visualize the structure, the easiest way is to start with a sheet of graphene. Graphene is a single layer of carbon atoms arranged in hexagons. In this structure, each carbon atom is sp 2 -hybridized. This means that three σ orbitals are arranged in a 2D structure having three 120 • angles. Therefore, each of these carbon atoms in the sheet of graphene has three covalent bonds to its neighboring atoms.
In order to obtain a carbon nanotube, the graphene layer has to be rolled up until a cylinder is formed. Therefore, it is called the rollup-model. This procedure describes a single-walled carbon nanotube, consisting of one layer of atoms. There also exist multi-walled CNTs which are made out of more than one layer.
Depending on the direction in which the graphene sheet is rolled up, different types of nanotubes can be categorized: (i) armchair, (ii) zigzag, and (iii) chiral nanotubes, see Fig. 1. To define the rollup direction of the graphene sheet, the chiral vector is used. For more details, the interested reader is, e.g., referred to Dresselhaus et al. [4]. To classify the type of the nanotube, an integer pair (n, m), resembling the chiral vector, is used. For the armchair configuration, n and m have the same value; for m = 0, a zigzag tube is formed; all other combinations are called chiral nanotubes. The higher the values of n and/or m are, the larger is the diameter of the nanotube. The height or length of the nanotube is defined by the quantity of unit cells stacked together. A unit cell (uc) is the basic construction element of the structure.
In the present research, the polyhedral model according to Cox and Hill [26], which is more sophisticated and accurate than the rollup-model, is applied.
After giving a short introduction of the geometry of a carbon nanotube, the electrical properties will be investigated more detailed. Therefore, the charge distribution of the tube has to be calculated. Various models investigating this topic can be found in the literature, e.g., electrostatic capacitance calculations by Rotkin et al. [27,28], ab initio calculations based on a density functional theory by Keblinski et al. [21] or energy-based methods by Ke et al. [29].
In the present research, three models are presented due to (i) their relative low computational effort, (ii) their good possibility to be calibrated with experimental results, and (iii) the ability to be easily compared to each other. In section 2.1, the fundamental electrostatic model [21] is presented. Based on this, an approach by Li and Chou [22,30] is discussed in Sect. 2.2. In Sect. 2.3, the model by Mayer [24] is outlined.

Classic electrostatic model
This Section shortly discusses the origin of the models applied in the present paper. The more advanced versions of the charge models are referred to the work by Keblinski et al. [21]. They used a density functional calculation and local density approximation (LDA) for their simulations. These results are compared with a classical electrostatic approach. More detailed information about the procedure can be found in their publication [21].
By the time this the work was published, the experimental distance between two carbon atoms with a sp 2hybridization was recorded approximately as d c−c = 1.2 Å; but Keblinski et al. [21] calculated the distance in their research as d c−c = 1.31 Å.
More recent investigations show that the C-C bond length is larger than previously calculated. In the works by Eberhardt and Wallmersperger [31,32], the bond distance is assumed to be d c−c = 1.421 Å. Therefore, in the present research the effect of the bond length on the charge distribution will be investigated.
In the publication by Keblinski et al. [21], the charge distribution of an (5,5) open-ended armchair nanotube consisting of 110 atoms is calculated. In order to have a reference configuration, the same type of CNT is used in the present research as in Ref. [21].
To calculate the electrostatic potential V * 1 , the following ansatz is used: Every carbon atom i of the nanotube is surrounded by a spherically distributed charge, see Fig. 2b. r i, j describes the distance between the atoms i and j, q j and q i are the charges of the atoms j and i, respectively. The first term in Eq. (1) represents the self potential of each atom. Keblinski [21] et al. fitted the parameter A to a good accordance with the quantum calculations of A = 0.78 Å −1 . In the following, this approach will be referred to as the classical model.

Electrostatic model by Li and Chou
Li and Chou [22,33] used a different approach than the classical one. They applied an atomistic moment method to calculate the distribution of the charges along the carbon nanotube. The main difference to the classical model is the area in which the electrostatic charge around the nucleus of the carbon atoms is spread. In Section 2.1, the charge of each atom is distributed in a sphere around the center, see Fig. 2b. In the model by Li and Chou, referred to as Li model, it is distributed within a triangle on the nanotube surface, see Fig. 2c. Similar to the previous approach, a description of the electric potential V (r i ) for each atomic position r i for each of the N atoms in the nanotube can be found. This can be formulated as where 0 is the permittivity of the free space, q j is the charge of the atom j, and r i, j again describes the distance between atom i and j in a tube consisting of N atoms. These N equations, one for each atom, can  [26], Here, {q} and {V 0 } denote the vector of the charge and the vector of the electric potential, respectively. The matrix [A] includes the interaction coefficients a i j . In the case of two different atoms i = j, Eq. (2) shows that a i j can be expressed as Please note that only the distance r i, j between two considered atoms i and j is relevant. The following expression considers the self-potential of each atom, i.e., if i equals j. Then, the triangle shape of the charge distributions becomes important. In this case, where i = j, the term for a i j results in The term s = 3 √ 3b 2 4 with b = 0.142 nm describes the triangular area in which the charge of the considered atom is distributed. A comparison of the model by Li and Chou and the model by Keblinski et al. [21] is performed in [22]. A good correlation between these two models is found. Therefore, in the present research, the newer and more sophisticated model by Li and Chou [22] is applied.

Charge-Dipole model by Mayer
In this Section, the model published by Mayer [24] is reviewed. In the following, it is referred as Mayer model. A short overview of this model is given in order to understand which parameters were used and how the solution procedure works.
This model describes every atom with a charge q i and also with a dipole p i . Subsequently, all calculations are performed with both, charges and dipoles. The consideration starts with the total energy E tot of all atoms in the nanotube at their given positions: Please note that Eq. (6) sums up the energy contributions of different interactions. The first term refers to the interaction between the charge of an atom and the charges of the other atoms. The sum is multiplied with 1/2 to avoid the double counting of each pair of atoms. The term is weighted with the interaction coefficient T i, j q−q . This coefficient can be calculated with Eq. (7) for i = j, Equation (7) looks similar to Eq. (4), which is considered for point charges. Mayer extended this approach and assumed that the charges are spread by a Gaussian distribution. The interaction coefficient T i, j q−q is affected by the vacuum permittivity 0 and the distance r i, j between atom i and j in a nanotube consisting of N atoms. The term "erf()" denotes the error function, where erf(x) = 2 √ π x 0 e −τ 2 dτ . The most important parameter in Eq. (7) is the width of the Gaussian distribution R. If R is small, the system reacts in a similar way as if the computation were performed with point charges, see Fig. 2a. If R is larger, the charge located on each nucleus is distributed with a Gaussian curve, see Fig. 2d. This results in a more evenly distribution of the charges over the whole nanotube. At this point, it has to be mentioned that the used indices at the interaction terms are not tensor notations. i and j in this case just denote the atoms for which the interaction coefficient is computed, () q−q denotes the interaction of two charges, () p−q describes the interaction of a charge and a dipole, and () p− p stands for the interaction of two dipoles.
For i = j, the self-energy coefficient is given by The different contributions in Eq. (6) are explained now in more detail: The second term in Eq. (6) describes the interaction of two dipoles, and the third term represents the interaction of a charge and a dipole. This can be calculated by for the second term and by for the third term.
Please note that r i, j is the vector from atom i to atom j, and I represents the identity matrix. These terms can also formulated for i = j by: The second to the last term in Eq. (6) for the total electrochemical energy describes the reaction of a charge to an applied external potential V i . This is the energy which is required to hold or press the charge in the position given by the external potential. In the fourth term, the interaction of the additional charge q i with the electrons of the atom i placed at the given position is accounted by χ i . In the publication by Mayer [24], χ i is defined as the electron affinity. One year later, in 2008, in the research by Mayer and Åstrand [34], it is called electronegativity. On the Pauling scale this would mean χ i = 2.55 for a carbon atom [35].
In the publication by Mulliken [36], it is shown that the electronegativity can be calculated by the ionization energy and the electron affinity. Unlike to the Pauling scale, this refers to an absolute electronegativity. Bultnick et al. [37] present a chart with different values for the electronegativity from χ = 2.47 to χ = 8.58, based on different possible calculation approaches.
Please note that the exact value of χ i is nonrelevant as it cancels out if only one kind of atoms is considered. The last term of Eq. (6) defines the interaction of the dipoles and the applied external electric field E.

Comparison of the different modeling approaches
In order to show the details of the given models, different test cases are investigated. In the first test case, a carbon nanotube which is charged by an overall charge, i.e., a total electrical load, is considered. There, it can be observed how this total electrical load is distributed between the atoms of the carbon nanotube.
In Sect. 3.2, the second test case is formulated. In this case, the tube is exposed to an external electric field. In the last Section, the third test case will combine both effects, i.e., a charged carbon nanotube is exposed to an electric field.

Test case I: overall charge
In this first test case, a carbon nanotube is charged with an overall charge. At first the behavior of the CNT having an overall charge of 5e is investigated, where e = 1.6022 · 10 −19 C is the elementary electric charge. Afterward, different overall charges are compared in order to understand the electrical behavior of CNTs. Thereafter, the different models, mentioned in Sect. 2, are compared to analyze the differences. Both armchair and zigzag nanotubes are considered. Each of these CNTs has a length of five unit cells. One unit cell defines the basic construction element of a structure. Please note that for the zigzag CNT not exactly 5 unit cells are shown. In order to represent the symmetric behavior, the first row of atoms on the left side is omitted. The considered CNTs are open ended, i.e., they do not have any cap. In Fig. 3b, it can be seen that after two rows of atoms which are very close together the distance to the next row is slightly larger. The displayed results show that the charge distribution depends on the arrangement of the atoms in the nanotube.
Two graphs are depicted in Fig. 3 (bottom) directly underneath the tubes. Each of these graphs shows the charge in elementary electric charge e with respect to the position of the atoms in the CNT in Å. The data points are aligned with the atoms above. Atoms placed on the same position in z-direction have the same charge, due The nanotubes shown in Fig. 3 have an overall charge of 5 elementary charges. This total charge is distributed over the CNT as shown in the graph. Please note that, due to the distributed charges, also fractions of the elementary charge occur at certain local positions in the tube. The applied test case is just an exemplary case in order to demonstrate the general behavior. Therefore, we decided to use this first example as a fundamental example which can be extended to higher overall charges. In Fig. 3, it can be clearly seen that the atoms at the end of the tube have a higher charge compared to the average value. In the middle region of the nanotube, the charges are distributed almost evenly. The atoms of the considered nanotube are sp 2 -hybridized, i.e., each atom strives to be bonded to three other atoms. In the middle part of the nanotube, this criterion is fulfilled, i.e., all bindings are bounded. At both ends of the nanotube, due to the open end configuration, not all bindings are complete, i.e., there are free electrons on the carbon atoms which do not have all binding partners. The free electrons attract the additionally applied charge; this means that the charge is more concentrated at the ends of the nanotube. Please note that the total charge has been kept constant for the entire numerical simulation. If all charges are multiplied with the respective number of atoms and summed up, the result is exactly 5e.
In the next investigation, an overall charge of 20e is applied, see Fig. 4. As in Fig. 3, it can be seen that the relative distribution of the charge is the same, i.e., it is independent of the total charge. If the overall charge increases while the total number of atoms remains constant, the charge on each atom has to be larger. For example: the value of the charge of the very last atoms on the tube (a) increases from 0.0857e in Fig. 3 to 4 · 0.0857e = 0.3429e in Fig. 4; for the ten atoms in the middle it increases from 0.0287e to 0.1146e.
It can be seen that the (5,0) zigzag CNT shows small oscillations of the charge between neighboring atoms and the charge distribution is not as homogeneous as for the armchair nanotube. This can be justified with the value of the parameter R. A more detailed discussion of this will follow later on.
The correlation between higher overall charge and higher charges on each atom is depicted in Fig. 5. As shown in Figs. 3 and 4, again the charge distribution with respect to the longitudinal z-axis position of the atoms in the carbon nanotube is displayed. Five different curves are plotted in the graph: the same tube geometry is loaded with an overall charge of 1e, 2e, 5e, 10e, and 20e. As expected, it can be seen that with an increasing overall charge the charge of each atom increases. This behavior is independent of the tube configuration, i.e., whether it is an armchair tube or a zigzag one.
It can also be observed that the relative charge distribution between the ends and the middle of the tube increases with a higher prescribed overall charge.
As shown in Sect. 2, different approaches are used to represent the charge distribution on a carbon nanotube. Figure 6 gives the numerical results for a (5,5) armchair and a (5,0) zigzag CNT according to the model by Li and Chou. In the graph, the normalized charge over the length is given. Normalized charge means that the actual charge of the atom is divided by the average charge per atom q = q tot N where N is the number of atoms in the nanotube. Figure 6a shows identical results compared to the work by Li and Chou [22]. This comparison was made to test the functionality of the programmed MATLAB application. In order to compare the modeling approaches by Li and Chou and by Mayer, both results are shown in Fig.  7. For the armchair configuration as well as for the zigzag type, the charge distribution in the middle of the CNT is nearly the same for both modeling approaches.
At both ends of the nanotubes, a distinct difference can be observed. For both tube types, the charges of the last two rows of atoms show significant different values. This can be explained with the distribution approach of the charges at the nuclei in the different modeling approaches. As described in Sect. 2, the model by Li and Chou [22] distributes the charge over a triangular area. Mayer [24] on the other hand distributes the charge by using a Gaussian curve with the parameter R.
The approach by Mayer [24] seems to provide more homogeneous results. The increase in the charge of the last two rows of atoms in the armchair nanotube is less high than in the numerical simulation based on the model by Li and Chou. Also for zigzag nanotubes, the course of the charge seems to be more realistic. At both ends of the (5,0) zigzag CNT, the last two rows of atoms are close together. The model by Li and Chou shows the highest charge value at the last atom row. For the middle part of the nanotube, the value decreases and increases again toward the end of the CNT. Nonetheless, the model by Mayer shows a more comprehensible trend. The last rows of atoms at the ends have the highest value, and the next row has a slightly lower value because of the small atom distance. Between the second row and third row, the decrease is higher due to the arrangement of the atoms in the nanotube.
Longer nanotubes show the same qualitative behavior as short nanotubes. Figure 8 gives the results of the same setup as Fig. 7, but with a different nanotube length of 20 unit cells. Likewise the values of the charge of the atoms in the middle of the nanotubes are the same for both models. The differences, shown in Fig. 8a, occur in the two rows of atoms facing both ends and in Fig. 8b in the last four rows of atoms. This is a similar amount of atoms as in the investigations performed for shorter nanotubes. The modeling approaches are therefore similar except of the charge accumulation at the ends of the tubes. Regardless of the length of the carbon nanotubes, only the last few rows of atoms are affected by the differences of the charge distribution.
In Fig. 9, the effects of different values of the parameter R, used in the Mayer model, are shown. Also the charge distribution with respect to the position of the atoms of the nanotube is given. The results for the (5,5) armchair nanotubes are given in Fig. 9a and for the (5,0) zigzag nanotubes in Fig. 9b. The value of R, which shows the range in which the charge is located around the nucleus, is varied from R = 0.2 Å to R = 0.9 Å.
The parameter R is the decisive factor of the model by Mayer [24] in order to demonstrate how the charge on each atom is spread over the area around it. If R therefore is low, the Gaussian curve is narrow, and the atoms behave more like a point charge, i.e., the curve in Fig. 9 flattens. If R is nearly zero, each atom carries the same amount of charge, which equals to the total charge over the tube divided by the amount of atoms of the CNT. All points of the curve then align to a simple horizontal line.
On the other hand, if R is quite large, the Gaussian curve flattens. That means that the areas in which the charges are spread overlap each other. This results in alternating charge values in each row of atoms as shown in Fig. 9a and b for R = 0.9 Å. This effect and therefore the maximal value of the numerically obtained charges increase with an increasing R.
The curve for R = 0.6863 Å is indicated by the symbol " * ". This is the usual value for R as used for the calculation of all test cases. Please note that this value is also used in the research by Keblinski [21], Mayer [24], and Li [22]. In their research, this was found to be a good agreement between ab initio calculations and the performed molecular mechanics simulations.
In Fig. 10, the influence of the bond length d c−c -as mentioned in Sect. 2.1-is depicted. For a smaller bond length, the whole tube is shorter if the same configuration is used. Therefore, the distance between the rows of atoms is shorter. In the curves given in Fig. 10, the model by Mayer is used. Therefore, the interaction between the width of the Gaussian curve with the parameter R and the bond length d c−c can be observed. As mentioned previously, R describes how the charge belonging to one nucleus is distributed over the nanotube. If R is fixed and the distances between rows of atoms is changed, the same effects as shown in Fig. 9 occur due to the ratio between R and d c−c . If R is fixed at R = 0.6862 Å and the bond length d c−c decreases, a similar behavior like in the case of a fixed distance and an increasing R is observed. Therefore, if the bond length decreases, the charge begins to fluctuate as the area in which the charges are spread, overlaps. Please note, that-due to a better visibility-in Fig. 10b some data points are not displayed. On the other hand, with an increasing atomic distance the interaction declines. If the distance is further increased-which is not shown in order to improve the clarity of the graph, the charge distribution flattens like it does if R is decreased. An overview of the numerical investigations (test case I) and the used parameters are given in Table 1.

Test case II: applied external field
In order to understand the electric behavior of the CNTs, not only the reaction due to an overall charge but also due to an applied electric field has to be investigated. This will be discussed in the present Section. With the knowledge of the electric properties, it is possible to create an electromechanical model of a carbon nanotube. This could lead to a more efficient construction for, e.g., actuators made out of CNTs. In the different investigations, electric fields of E = 0.1 V/Å, 0.5 V/Å, 1 V/Å, 5 V/Å and 10 V/Å are applied to an (5,5) armchair CNT and a (5,0) zigzag nanotube with five unit cells, respectively. In the following, it is shown how   Table 2.
In Fig. 11, the charge on each atom with respect to the position on the z-axis of a tube of 20 unit cells is shown. In Fig. 11a, an armchair tube and in Fig. 11b a zigzag CNT is shown. Both tubes have been charged with an ideal electric field of E = 0.1 V/Å in (longitudinal) z-direction. In difference to Sect. 3.1, there is no overall external charge on the tube.
In Fig. 11a, the general behavior can be observed. As already shown in the simulations using a total charge, the values of the charges in the middle of the tube are not as high as at the ends. It can be seen that-due to the applied electric field-there is a slight decrease in the charge versus z. The last two rows of atoms behave differently. They show an increase in the charge values. All charges summed up equal to zero, as prescribed.
In Fig. 11b, the absolute values of the charges in the middle of the tube show an alternating behavior. The explanation for this can be found in the parameter R describing the width of the Gaussian distribution which characterizes the interacting of the charges in the model by Mayer [24]. In Refs. [22,24], R is stated for the armchair nanotubes but not for the zigzag tubes. All the simulations were performed with R = 0.6862Å. It might be possible that a slightly lower value of R would fit better for a zigzag CNT.
Next, the effect of different external electric fields is investigated. The results of the numerical simulation are shown in Fig. 12. A shorter tube is used consisting of only five unit cells. It is expected that the general behavior remains the same for a longer nanotube. If the electric field is larger, the values of the charges increase. They maintain the point symmetry around the exact middle of the tubes. As mentioned before, in Fig. 11(b) the values of the charge around the middle oscillate slightly, due to the influence of the parameter R. As already explained in Sect. 3.1, the reaction of the charges to the external field or to the total charge can be seen best in the last two rows of atoms toward the ends of the tubes. Figures 9 and 11 confirm that the length of the CNT does not affect this correlation.

Test case III: overall charge and applied external field
In most of the real applications, the factor of an overall charge and an applied external field cannot be considered completely separated. Therefore, in the last test case, both effects will be combined, in order to accomplish a complete impression of the electric behavior of armchair and zigzag nanotubes. First, the CNT is charged with an overall charge of 5 e. Then, an external electric field of E = 0.1 V/Å in (longitudinal) z-direction is applied to the carbon nanotube. Again, like in Sect. 3.1 both, armchair and zigzag carbon nanotubes are investigated. Furthermore, the effect of increasing the overall charge combined with an applied external field will be discussed.
The parameters used for this test case (test case III) are given in Table 3.
The charge with respect to the position of the atoms of the CNT is depicted for armchair in Fig. 13a and for zigzag tubes in Fig. 13b. According to the total charge of 5 e, the atoms in the middle region are charged similarly as in Fig. 3. A larger difference can be observed at both ends of the nanotubes. It can be seen that the effects of the external field and of the total charge are combined. At the left end of the tubes, both influences are summed up. On the right side, the contributions work against each other. The combination of a decreased charge due to the electric field and an increase in the charge due to the applied overall charge tends to result in lower values of charges at the ends of the tube.
In the next step, the external electric field is kept constant E = 0.1 V/Å and the overall charge on the tubes varies, see Fig. 14. As expected, the same behavior like in Fig. 5 occurs, with the additional effect depicted in Fig. 13.

Preliminary assessment of the test cases with respect to electromechanical coupling
In test case I, the nanotubes were charged with an overall charge. Therefore, atoms at the same Z-position have the same value and polarity of charge. Since the same charges repel each other, it is expected that the nanotube will expand radially. Additionally, the tube will expand along the z-axis, due to the distributed charges. Due to the higher charges, and therefore the higher repulsive forces, the extensions at both ends of the CNT are expected to be larger than in the middle region. This extension of the nanotube will again lead to a redistribution of the charges, i.e., especially at the ends of the CNT, the charge accumulation is expected to be reduced, due to the larger resulting atomic distances.
In the last test case, an electric field was applied to an already charged nanotube. As mentioned before, the nanotube expands radially due to the same polarity of the charges at each atom of the same z-position. Also in each half of the tube the charges are polarized in the same way; therefore, they react repulsive, and the tube will get longer. For a high electric field strength, both ends of the tube are charged in the opposite way resulting in an attractive effect. This could compress or deform the tube.
In further studies, it will be investigated which of both effects-repelling or attracting-are stronger. Therefore, a coupled electromechanical formulation will be derived and studied intensively.

Conclusions
In the present research, a systematic analysis of the behavior of charged carbon nanotubes exposed to electric fields is performed. Two different modeling approaches [22,24] are compared with each other. The main focus of this work is the reaction of armchair and zigzag CNTs: (i) to an applied overall charge, (ii) to an applied external electric field, and (iii) to the combination of both effects. Additionally, the effect of the carbon bond length and the effect of the width R of the Gaussian distribution in which the charge on the nucleus is spread is investigated. Concluding, the present work investigates the charge distribution over the tubes.
An applied overall total charge of the nanotube shows that the charge is symmetrically distributed on the CNT, independent of the tube configuration. The highest values of charge can be found at the ends of the nanotube. Even if the CNT is constructed by a large number of unit cells, i.e., it is quite long, only the last few atoms are higher charged.
This research compares two different modeling approaches, the model by Li and Chou [22,33] and the model by Mayer [24]. Both models give similar results with respect to the overall behavior, but there are significant differences in the values of the charges at the nanotube ends. The model by Li and Chou shows a more extreme step to higher values of the charge at the last row of atoms of the tube. The model by Mayer is more homogeneous, i.e., the increase is not as high as in the model by Li and Chou. It is also observed that the choice of the width of the Gaussian distribution R in the model by Mayer [24] is important. If it is too low, the behavior is similar to point charges. If R is too high, the charge is distributed quite evenly all over the nanotube, i.e., the charges interact largely with each other, and also oscillations of the charge occur. Besides, the parameter R has to be adjusted to the length of the carbon bond d c−c and to the configuration of the CNT for both, armchair and zigzag tubes. In the present research, R = 0.6862 Å is used, according to other investigations available in the literature. This value fits very well for armchair nanotubes. For a zigzag configuration, a smaller value of R should be chosen in order to eliminate the small oscillations of the charge.
If an electric field is applied, the charges are distributed point symmetrically with respect to the middle cross-sectional area of the nanotube. The effect of an overall charge and an electric field superpose each other if both are applied.
In further studies, the electromechanical coupling, i.e., the deformation of the tube due to the charge distribution, will be investigated. By applying an electromechanical model, e.g., for nanotweezers, the superposition of attracting and repelling effects, leading to deformation (expansion, compression or bending) as a consequence of the applied charge and the applied electric field, will be investigated.