- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Journal of Nanomaterials
Volume 2011 (2011), Article ID 805313, 10 pages
Analysis of Carbon Nanotubes on the Mechanical Properties at Atomic Scale
1Interdisciplinary Graduate School of Science and Technology, Shinshu University, 3-15-1 Tokida, Ueda, Nagano 386-8567, Japan
2Department of Functional Machinery and Mechanics, Shinshu University, 3-15-1 Tokida, Ueda, Nagano 386-8567, Japan
3Institute of Road and Bridge Engineering, Dalian Maritime University, 1 Linghai Road, Ganjingzi, Dalian 116026, China
Received 14 June 2010; Accepted 23 August 2010
Academic Editor: Junping Wang
Copyright © 2011 Xiaowen Lei et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This paper aims at developing a mathematic model to characterize the mechanical properties of single-walled carbon nanotubes (SWCNTs). The carbon-carbon (C–C) bonds between two adjacent atoms are modeled as Euler beams. According to the relationship of Tersoff-Brenner force theory and potential energy acting on C–C bonds, material constants of beam element are determined at the atomic scale. Based on the elastic deformation energy and mechanical equilibrium of a unit in graphite sheet, simply form ED equations of calculating Young's modulus of armchair and zigzag graphite sheets are derived. Following with the geometrical relationship of SWCNTs in cylindrical coordinates and the structure mechanics approach, Young's modulus and Poisson's ratio of armchair and zigzag SWCNTs are also investigated. The results show that the approach to research mechanical properties of SWCNTs is a concise and valid method. We consider that it will be useful technique to progress on this type of investigation.
Since the discovery by Iijima in 1991 , carbon nanotubes (CNTs) have generated huge activities in most areas of science and engineering due to their unprecedented mechanical, electrical, and thermal properties [2–13]. Especially in mechanical field, experimental measurements have determined that CNTs possess excellent mechanical properties [14–21]. Therefore, an effective method to analyze the basic characteristics of nanosized CNTs is essential.
In the past, researchers used experimental method to measure mechanical properties of CNTs. Treacy et al.  firstly measured the amplitude of intrinsic thermal vibrations observed in transmission electron microscopy (TEM). The average value of Young’s modulus of CNTs derived from this experimental technique is 1.8 TPa by 11 tubes, in which the lowest value and the highest value are 0.40 TPa and 4.15 TPa, respectively. Later, Poncharal et al.  obtained Young’s modulus of CNTs which is between 0.7 and 1.3 TPa by electromechanical resonant vibrations. In addition, based on an atomic force microscope (AFM), Wong et al.  in 1997 firstly directly measured the stiffness constant of armchair MWCNTs pinned at one end, from which the value of Young’s modulus of CNTs is 1.28 TPa. Salvetat et al.  used the AFM for experiment of bending an armchair (multi-walled carbon nanotubes) MWCNT pinned at each end over a hole and obtained an average modulus value of CNTs of 0.81 TPa. These experiments all promote the research of mechanical properties of CNTs. However, in description of nanoscale structures, the results are with experimental errors.
Meanwhile, for researching mechanical properties of CNTs, a number of researchers solved the difficulties in nanosized experiments in terms of computer simulation. For the analysis of nanostructural materials, atomic simulation methods such as first-principle quantum-mechanical methods , molecular dynamics (MD) [23, 24], and Monte Carlo  simulations have been routinely adopted. As early as 1993, Overney et al.  calculated Young’s modulus of rigid short SWCNTs which is 1.5 TPa, approximately equal to that of graphite. This was followed by a range of papers predicting that Young’s modulus of CNTs is close to 1 TPa independent of type and diameter . Yakobson et al.  fitted these results by MD simulations of the continuum shell model. Unlike the previous work that assumed a thickness of 0.34 nm, both the thickness and Young’s modulus were taken as the fitting parameters, yielding a thickness of 0.066 nm and Young’s modulus of 5.5 TPa. The MD approach was also used by Lu [27, 29] who reported that Young’s modulus is 1 TPa and claimed that chirality and the number of walls have little effect on the value of Young’s modulus. A different potential model was used by Yao and Lordi  who obtained Young’s modulus of CNTs as 1 TPa. Although MD method has been widely used in simulating the properties of nanostructural materials, it is complex and time consuming, especially for large amount atomic systems. Therefore, the continuum mechanics seems to be a better way to investigate the properties of CNTs. However, for the case of nanoreinforced adhesives, these models cannot accurately describe the influence of the relationship between carbon atoms upon the mechanical properties and their interactions in the composite systems because they lack the appropriate constitutive relations that govern material behavior at this scale . Another modeling approach is the atomic-based continuum technique, which has the unique advantage of describing atomic structured properties in a continuum framework for reducing the computational demand while employing the appropriated atomic constitutive relations.
Therefore, there is a demand of developing a modeling technique that could analyze the mechanical properties of CNTs at the atomic scale. Considering CNTs as a rolled cylindrical graphite sheet, we step from Young’s modulus of the C–C bonds counted as Euler beam at atomic scale and extend the theory of classical structural mechanics into the modeling of carbon graphite sheet. The effects of tube curvature on the mechanical properties of SWCNTs are considered in closed-form solutions. The mechanical properties of SWCNTs, including Young’s modulus, Poisson’s ratio, the length of C–C bonds and the angle between the adjacent C–C bonds are discussed as functions of nanosized structure.
2. The Structure of CNTs
CNTs can be considered as graphite sheets rolled into cylindrical shape. The one-atom-thick graphite sheet looks like chicken wire which is made of a single-carbon-atom thickness. The structure of CNTs as shown in Figure 1 is conveniently explained in terms of the chiral vector integers (, ) and the chiral angle :
CNTs are classified into three categories named as zigzag (, 0): ; armchair (, ): ; and chiral (, ): . The relationship between radius and integers (, ) is expressed as where is the length of C–C bonds, 0.142 nm, for SWCNTs.
3. Mechanics Model of Graphite Sheet
From the viewpoint of molecular mechanics, CNTs are treated as a large array of molecules consisting of carbon atoms. According to the Tersoff-Brenner force field theory [3, 32], the total molecular potential energy of the force field for nanostructured material can be defined as the sum of bonding and nonbonding energies: where is the energy induced by bond stretch interaction, is the energy by the bond angle bending, is the energy by the dihedral angle torsion, is the energy by improper torsion, is the energy by the nonbonding van der Waals (vdW) interaction, and is the energy by the electrostatic interaction. As the axial loading, improper torsion and nonbonding interactions subjected to CNTs are very small , for covalent systems, the main contributions to the total molecular potential energy come from the first four terms of (3). Simplified system potential energy of CNTs with C–C bonds is given as where , , and are the bond stretching, bond bending, and torsional resistance constants, respectively, while , , and represent the bond stretching increment, the bond angle change, and the angle change of bond twisting, respectively.
According to classical structural mechanics, the strain energy of a uniform beam in graphite sheet is expressed as where is the tensile resistance of beam (C–C bonds), is the flexural rigidity, is the torsional stiffness, and , , and are the axial stretching deformation, the rotational angle at the end of the beam, and the relative rotation between the ends of the beam, respectively.
Based on energy conservation law, a linkage between the force constants in molecular mechanics and the sectional stiffness parameters in structural mechanics is established. Equations (4) and (5) are equal to each others in one-to-one corresponding directions. The direct relationship between the structural mechanics parameters and the molecular mechanics parameters is deduced as follows:
As long as the force constants , , and are known, the sectional stiffness parameters , , and can be determined, then the deformation and elastic behavior of CNTs at atomic scale can be modeled. By further analysis of (4) and (5), three more constants of C–C bonds (diameter of C–C bonds), , and depend on , , and :
3.1. Young’s Modulus for Armchair Graphite Sheet
Knowledge of Young’s modulus (E) is the first step towards the material using as a structural element for various applications. SWCNTs can be regarded as a two-dimensional continuum shell which is composed of discrete molecular structures linked by C–C bonds. The unrolled graphite sheet of armchair SWCNTs is shown in Figure 2. Figure 2(b) plots the smallest unit of armchair graphite sheet, in which the force and moment are displayed. The unit of armchair graphite sheet can be analyzed based on solid mechanics, and the unit along BD bond is symmetry of both structure and force. Therefore, in terms of the elastic deformation energy, the energy of an armchair unit is written as where is Young’s modulus of beam and is the C–C bonds of armchair graphite sheet. When the relationship between force and moment is determined, the strain of the unit can be resolved by using energy .
According to Castigliano’s Law, the rotation angel on point is zero due to symmetry of structure and force that is shown in
Substituting (8) into (9), the relationship between force and moment is obtained: Meanwhile, the elastic deformation energy is rewritten as follows: Based on Castigliano’s law, the displacement of the unit is defined by the elastic deformation energy differential and the strain of armchair graphite sheet is defined as follows:
For an armchair graphite sheet being subjected to the tension stress , the equation of static equilibrium of the unit is expressed as
Based on Hooke’s law, the relationship of the tension stress and the strain is When , Young’s modulus of armchair graphite sheet yields where is the thickness of the graphite sheet.
3.2. Young’s Modulus for the Zigzag Graphite Sheet
For a zigzag graphite sheet, the analytical approach is similar to that of the armchair graphite sheet. The unit of a zigzag graphite sheet is drawn in Figure 3, and the energy of the unit is drawn in where is the C–C bonds of the zigzag graphite sheet. The displacement of the zigzag unit is defined as
Then the strain of per unit length is obtained as follows:
For a zigzag graphite sheet being subjected to the tension stress , the force acting on the unit can be written as
When , we can obtain Young’s modulus of zigzag graphite sheet:
4. Molecular Mechanics Model of SWCNTs
SWCNTs can be ideally constructed starting from a graphite sheet. According to the chiral vector, there are three kinds of structure, as shown in Figure 4. After the two-dimensional sheets are rolled into a three-dimensional tube, intrinsic properties of the structure will be changed, such as mechanical and physical characteristics.
4.1. Angle Measurement for SWCNTs
Figure 4(a) shows an armchair SWCNT. According to (2) the radius of armchair SWCNTs is . Meanwhile, by geometric relationships in cylindrical coordinates shown in Figure 4(a), we obtain the coordinates of carbon atoms:
Because of the effect of curvature, the lengths of the bonds between carbon atoms in vector space are shorter than those in plane and characterized in different numerical values, which are displayed in Figure 4(a) and given as follows:
In cylindrical coordinates, on account of the curvature, the included angles between the adjacent bonds in a carbon cycle of SWCNTs are different from those of graphite sheets. In accordance with geometrical relationships, the included angle between the adjacent bonds for an arbitrary unit of armchair SWCNTs plotted in Figure 4(a) relates with coordinate figures of carbon atoms in the unit:
The distances of the bonds in vector space are as follows:
For zigzag SWCNTs, the analysis method resembles armchair SWCNTs; we calculate the included angle for zigzag SWCNTs as follows:
4.2. Elastic Properties for SWCNTs
Considering space curvature, Young’s modulus of SWCNTs in cylindrical coordinates depends on bond length and included angle between two bonds. The method for calculating Young’s modulus of armchair SWCNTs in three-dimensions resembles that in two dimensions. Affecting factors in three dimensional coordinates of Young’s modulus of armchair SWCNTs are mentioned in (23) and (25); Young’s modulus of armchair SWCNTs is expressed in where , , , and .
We obtain Young’s modulus of zigzag SWCNTs given in the following in the same way where , , and .
4.3. Poisson’s Ratio for SWCNTs
For armchair SWCNTs, the stretch deformations of the bonds caused by concentrated force and bending moment are schematically signed in Figure 5(a) to analyze Poisson’s ratio. The equilibrium of force and bending moment can be computed from Figure 5(a)
According to geometric properties and elastic theory, taking the included angle between two bonds in three dimensions into consideration, we obtain the equilibrium equations about extension variation of the C–C bonds where and are the tensile resistance and the flexural rigidity of beam, and , , and are the deformations of beam caused by force and bending moment , respectively. The strains in axial direction and circumferential direction are defined as follows, respectively:
For zigzag SWCNTs, the analysis step is similar to that for armchair SWCNTs. To analyze Poisson’s ratio of zigzag SWCNTs, Figure 5(b) schematically signs the stretch deformations of the bonds caused by force and bending moment ; the relationship of force and bending moment is given in
The equilibrium equations about extension variation of the bond are described as
The strains in axial direction and circumferential direction are defined as follows, respectively:
Poisson’s ratio of armchair SWCNTs can be defined as the ratio between circumferential strain and axial strain; then we obtain
5. Results and Discussions
The atomic-based continuum mechanic approach described in the previous section was implemented for studying the effective elastic properties of graphite sheets and SWCNTs. In this section, the mechanical characteristics of graphite sheet and SWCNTs are examined.
5.1. Mechanics Model of Graphite Sheet
In the present simulation, kcal/mole/, kcal/mole/ and kcal/mole/, are taken . and are 0.142 nm and 0.34 nm, respectively. According to (6) and (7), the constants of beam are obtained and summarized in Table 1. Comparing the evaluated elastic moduli of graphite sheets with the literature results, Tserpes  and Kalamkarov  reported TPa and TPa, TPa and TPa of C–C bonds, respectively.
Comparing (16) to (21), Young’s moduli of the armchair graphite sheet and the zigzag graphite sheet have the same expression form except for and . Because the carbon cycle is a regular hexagon, graphite sheets have . We obtain Young’s moduli of both armchair and zigzag graphite sheets as 1.0424 TPa, which is close to the value of 1.033 TPa calculated by Li and Chou  and the value of 1.04 TPa computed by Shokrieh and Rafiee . It can be seen that the predicted values obtained from the model at the atomic scale agree well with those reported in the literatures.
5.2. Molecular Mechanics Model of SWCNTs
Since SWCNTs are defined as rolled graphite sheets, the lengths of bonds in vector space are changed owing to the effect of curvature. For armchair SWCNTs rolled by armchair graphite sheets, all of the lengths of C–C bonds become shorter because of the connection with curvature in circumferential direction, while, for zigzag SWCNTs rolled by zigzag graphite sheet, the length is only in vertical direction equal to the C–C bonds length of graphite sheets. According to (23) and (27), the variation of C–C bonds lengths with radius are shown in Figure 6. With increasing radius of SWCNTs, the bond lengths , , and approach to the length nm which is the bond length of graphite itself, owing to the curvature approaching to zero in infinite radius.
When graphite sheets are rolled into SWCNTs, the lengths and spatial relations of C–C bonds change obviously. Figure 7 shows the variation of included angles of two adjacent C–C bonds with radius according to (25) and (28). It expresses that the included angles are quite sensitive to SWCNTs in small radius. The included angle of zigzag SWCNTs is a little larger than that of armchair SWCNTs in the same radius. As being larger than 1.1 nm, the radius makes less influence on the included angles of both armchair and zigzag SWCNTs, and the two lines show a tendency to be 120 degree which is close to that of graphite sheets, due to the SWCNTs being regarded as graphite sheet when radius is infinite.
Comparing (29) to (30), Young’s moduli of armchair and zigzag SWCNTs have the same expression form except for coefficients , , and , , . As the two Young’s moduli shown in Figure 8, Young’s moduli predicted by the present theory decrease monotonically with the increase of radius. Young’s modulus of armchair SWCNTs is slightly higher than that of zigzag SWCNTs with the same smaller radius. When the radius is larger, the two declining curves gradually develop into horizontal lines, which tend to be a constant of 1.0424 TPa. The value is Young’s modulus of graphite sheet.
Considering the variations of length and included angles of two adjacent C–C bonds, from graphite sheets into SWCNTs, Figure 9 shows the variations of Poisson’s ratio of the armchair and zigzag SWCNTs as a function of nanotube radius. It is found that, in general, Poisson’s ratio of armchair SWCNTs is slightly higher than that of zigzag SWCNTs. Poisson’s ratios of armchair and zigzag SWCNTs show a smoothly monotonic decreasing versus increasing radius. When the radius of SWCNTs is larger, the two curves tend to be a straight line which means that Poisson’s ratios approach to be a steady value of 0.26, which makes a great agreement with those reported in the literatures [6, 33].
The mechanical properties of both armchair and zigzag SWCNTs are characterized by using continuum mechanics in the atomic scale. In terms of the conjunction of Tersoff-Brenner force field method and energy conservation law, the graphite sheet is of isotropic property and Young’s modulus of graphite sheet is obtained to be 1.04 TPa. Furthermore, considering the variations of the length and the included angle of two adjacent C–C bonds rolled from graphite sheets into SWCNTs, Young’s modulus and Poisson’s ratio of SWNCTs with armchair and zigzag structures are investigated as a function of SWCNT radius. We predict that Young’s modulus and Poisson’s ratio of SWNCTs are influenced obviously by relatively smaller radius while being little affected by larger radius. We are confident that this model provides a useful method to analyze mechanical properties of CNTs and other nanosized structures at the atomic scale.
This work was supported by Grant-in-Aid for Global COE Program by the Ministry of Education, Culture, Sports, Science and Technology.
- S. Iijima, “Helical microtubules of graphitic carbon,” Nature, vol. 354, no. 6348, pp. 56–58, 1991.
- F. Li, H. M. Cheng, S. Bai, G. Su, and M. S. Dresselhaus, “Tensile strength of single-walled carbon nanotubes directly measured from their macroscopic ropes,” Applied Physics Letters, vol. 77, no. 20, pp. 3161–3163, 2000.
- W. H. Chen, H. C. Cheng, and Y. L. Liu, “Radial mechanical properties of single-walled carbon nanotubes using modified molecular structure mechanics,” Computational Materials Science, vol. 47, pp. 985–993, 2010.
- Z. Spitalskya, D. Tasisb, K. Papagelisb, and C. Galiotis, “Carbon nanotube—polymer composites: chemistry, processing, mechanical and electrical properties,” Progress in Polymer Science, vol. 35, pp. 357–401, 2010.
- M. A. L. Manchado, L. Valentini, J. Biagiotti, and J. M. Kenny, “Thermal and mechanical properties of single-walled carbon nanotubes-polypropylene composites prepared by melt processing,” Carbon, vol. 43, no. 7, pp. 1499–1505, 2005.
- T. Natsuki, K. Tantrakarn, and M. Endo, “Effects of carbon nanotube structures on mechanical properties,” Applied Physics A, vol. 79, no. 1, pp. 117–124, 2004.
- M. R. Falvo, G. J. Clary, R. M. Taylor II et al., “Bending and buckling of carbon nanotubes under large strain,” Nature, vol. 389, no. 6651, pp. 582–584, 1997.
- K. I. Tserpes and P. Papanikos, “Finite element modeling of single-walled carbon nanotubes,” Composites Part B, vol. 36, no. 5, pp. 468–477, 2005.
- A. L. Kalamkarov, A. V. Georgiades, S. K. Rokkam, V. P. Veedu, and M. N. Ghasemi-Nejhad, “Analytical and numerical techniques to predict carbon nanotubes properties,” International Journal of Solids and Structures, vol. 43, no. 22-23, pp. 6832–6854, 2006.
- O. L. Blakslee, D. G. Proctor, E. J. Seldin, G. B. Spence, and T. Weng, “Elastic constants of compression-annealed pyrolytic graphite,” Journal of Applied Physics, vol. 41, no. 8, pp. 3373–3382, 1970.
- M. M. Shokrieh and R. Rafiee, “Prediction of Young's modulus of graphene sheets and carbon nanotubes using nanoscale continuum mechanics approach,” Materials and Design, vol. 31, no. 2, pp. 790–795, 2010.
- K. N. Kudin, G. E. Scuseria, and B. I. Yakobson, “ , BN, and C nanoshell elasticity from ab initio computations,” Physical Review B, vol. 64, no. 23, Article ID 235406, 10 pages, 2001.
- O. Lourie and H. D. Wagner, “Evaluation of Young's modulus of carbon nanotubes by micro-Raman spectroscopy,” Journal of Materials Research, vol. 13, no. 9, pp. 2418–2422, 1998.
- M. M. J. Treacy, T. W. Ebbesen, and J. M. Gibson, “Exceptionally high Young's modulus observed for individual carbon nanotubes,” Nature, vol. 381, no. 6584, pp. 678–680, 1996.
- J.-P. Salvetat, J.-M. Bonard, N. B. Thomson et al., “Mechanical properties of carbon nanotubes,” Applied Physics A, vol. 69, no. 3, pp. 255–260, 1999.
- N. Yao and V. Lordi, “Young's modulus of single-walled carbon nanotubes,” Journal of Applied Physics, vol. 84, no. 4, pp. 1939–1943, 1998.
- P. M. Ajayan, O. Stephan, C. Colliex, and D. Trauth, “Aligned carbon nanotube arrays formed by cutting a polymer resin-nanotube composite,” Science, vol. 265, no. 5176, pp. 1212–1214, 1994.
- P. Poncharal, Z. L. Wang, D. Ugarte, and W. A. de Heer, “Electrostatic deflections and electromechanical resonances of carbon nanotubes,” Science, vol. 283, no. 5407, pp. 1513–1516, 1999.
- M. R. Falvo, G. J. Clary, R. M. Taylor II et al., “Bending and buckling of carbon nanotubes under large strain,” Nature, vol. 389, no. 6651, pp. 582–584, 1997.
- E. W. Wong, P. E. Sheehan, and C. M. Lieber, “Nanobeam mechanics: elasticity, strength, and toughness of nanorods and nanotubes,” Science, vol. 277, no. 5334, pp. 1971–1975, 1997.
- J.-P. Salvetat, A. J. Kulik, J.-M. Bonard et al., “Elastic modulus of ordered and disordered multiwalled carbon nanotubes,” Advanced Materials, vol. 11, no. 2, pp. 161–165, 1999.
- F. Ding, “Theoretical study of the stability of defects in single-walled carbon nanotubes as a function of their distance from the nanotube end,” Physical Review B, vol. 72, Article ID 245409, 7 pages, 2005.
- K. M. Liew, C. H. Wong, X. Q. He, M. J. Tan, and S. A. Meguid, “Nanomechanics of single and multiwalled carbon nanotubes,” Physical Review B, vol. 69, no. 11, Article ID 115429, 8 pages, 2004.
- V. U. Unnikrishnan, D. Banerjee, and J. N. Reddy, “Atomistic-mesoscale interfacial resistance based thermal analysis of carbon nanotube systems,” International Journal of Thermal Sciences, vol. 47, no. 12, pp. 1602–1609, 2008.
- Y. Zhou, M. A. Baseer, H. Mahfuz, and S. Jeelani, “Monte Carlo simulation on tensile failure process of unidirectional carbon fiber reinforced nano-phased epoxy,” Materials Science and Engineering A, vol. 420, no. 1-2, pp. 63–71, 2006.
- G. Overney, W. Zhong, and D. Tománek, “Structural rigidity and low frequency vibrational modes of long carbon tubules,” Zeitschrift für Physik D, vol. 27, no. 1, pp. 93–96, 1993.
- J. P. Lu, “Elastic properties of single and multilayered nanotubes,” Journal of Physics and Chemistry of Solids, vol. 58, no. 11, pp. 1649–1652, 1997.
- B. I. Yakobson, C. J. Brabec, and J. Bernholc, “Nanomechanics of carbon tubes: instabilities beyond linear response,” Physical Review Letters, vol. 76, pp. 2511–2514, 1996.
- J. P. Lu, “Elastic properties of carbon nanotubes and nanoropes,” Physical Review Letters, vol. 79, no. 7, pp. 1297–1300, 1997.
- S. A. Meguid, J. M. Wernik, and Z. Q. Cheng, “Atomistic-based continuum representation of the effective properties of nano-reinforced epoxies,” The International Journal of Solids and Structures, vol. 47, pp. 1723–1736, 2010.
- C. Li and T.-W. Chou, “A structural mechanics approach for the analysis of carbon nanotubes,” The International Journal of Solids and Structures, vol. 40, no. 10, pp. 2487–2499, 2003.
- J.-L. Tsai, S.-H. Tzeng, and Y.-T. Chiu, “Characterizing elastic properties of carbon nanotubes/polyimide nanocomposites using multi-scale simulation,” Composites Part B, vol. 41, pp. 106–115, 2009.
- T. Natsuki, K. Tantrakarn, and M. Endo, “Prediction of elastic properties for single-walled carbon nanotubes,” Carbon, vol. 42, no. 1, pp. 39–45, 2004.
- W. D. Cornell, P. Cieplak, C. I. Bayly et al., “A second generation force field for the simulation of proteins, nucleic acids, and organic molecules,” Journal of the American Chemical Society, vol. 117, no. 19, pp. 5179–5197, 1995.