Ultrathin freestanding carbon nanotube (CNT) films Wu et al. (2004) possess a number of unique physical and mechanical properties, making them attractive for a number of interesting applications, in particular, in flexible electronics: transparent displays Park et al. (2013), mirrors for soft X-ray radiation Gubarev et al. (2019), wearable medical devices Wang and Loh (2017) etc. The major advantage of these materials, as compared to possible competitors, is their ability to withstand significant plastic deformations without irreversible structural changes. This property is conditioned by discontinuous nature of CNT films - they consist of individual CNTs that are chemically reactive and bound only by van der Waals forces. Due to molecular structure of CNTs, the interfaces between neighboring parallel CNTs resolve relative sliding of one CNT along the other, with very small static friction. Such sliding provides mechanism for plastic flow of CNT material. This mechanism is non-destructive in a sense that small strain plastic deformations do not lead to significant changes in microstructure of a CNT material, as well as in its mechanical, thermal, optical and electric properties. Another interesting feature is that such plastic deformations do not immediately localize in a small damage zone, but for a considerable strain diapason remain distributed over the length of a tested specimen.
The modern capabilities for experimental research on the properties of CNT materials are limited both in terms of technology and cost efficiency. Thus, it is highly desirable to have a numerical model that could reliably reproduce the mechanics of thin CNT films. Although molecular-level modeling techniques are capable to provide useful insights on the properties of individual nanotubes Yakobson et al. (1996); Dumitrica et al. (2006); Zhang and Dumitrică (2008); Nikiforov et al. (2010) they’re inefficient for modeling assemblies comprising large numbers of CNTs - the most extreme molecular dynamics calculations so far included only hundreds of CNTs Cornwell and Welch (2011).
In order to bypass the computational inefficiency of full molecular modeling techniques, a number of coarse-grained mesoscale models were developed Buehler (2006); Cranford and Buehler (2010); Mirzaeifar et al. (2015); Volkov and Zhigilei (2010a, b, c); Wittmaack et al. (2018)
. The general idea of these models is to reduce the number of model degrees of freedom (DOF) in order to reach the computationally tractable model size. However, such coarse-graining often leads to significant artifacts in model’s behavior. For example, pair vdW potential used inBuehler (2006); Cranford and Buehler (2010) leads to absence of CNT relative sliding, which makes the model inapplicable for modeling irreversible deformations of CNT assemblies. Similarly, giving up on torsional degrees of freedom in models Volkov and Zhigilei (2010a, b) significantly alters the mechanical properties of CNT assemblies under certain types of loading. Moreover, neither model so far provides a realistic description of energy dissipation during relative sliding of CNTs, which makes it impossible to describe rate-dependent deformation properties of CNT materials.
In our works Ostanin et al. (2013, 2014, 2015) we have developed the mesoscopic modeling technique based on distinct element method Cundall and Strack (1979). Within our approach, each CNT is represented as a chain of rigid bodies possessing tensile, bending, shear and torsional stiffnesses. These stiffnesses are calibrated based on molecular dynamics simulations. CNTs interact via vdW potential tailored to provide adequate description of relative CNT sliding. In our athermal model, energy dissipation was represented with viscous damping, appropriately calibrated in a top-down fashion Ostanin et al. (2014). The model showed to be efficient for the description of CNT self-folded configurations Wang et al. (2013, 2014), weakly bonded and cross-linked CNT ropes and bundles Ostanin et al. (2014, 2015) and qualitative description of large-strain deformations of CNT films Wang et al. (2018). Until recently, the model was lacking a scalable parallel implementation, which severely limited its capabilities in terms of achievable length and time scales of the simulation. We addressed this reservation with the development of a scalable parallel implementation of our model based on the rigid particle dynamics module of the waLBerla package Ostanin et al. (2018).
In our present work, utilizing our new parallel computation capabilities, we investigate the mechanics of representatively large assemblies of thin films. We present the study of the representative volume element (RVE) size, dependence of the mechanics of CNT films on dissipation rate, film thickness and length of individual CNTs. We restrict our consideration with the case of small strain deformation of CNT materials, with strains less than few percent. We could see that even at such strains CNT films exhibit complex nonlinear behavior, conditioned by CNT sliding within the tested specimen. Our observations are in good agreement with recent experimental studies on the mechanics of CNT films.
. For clarity, we present a brief description of our modeling technique here. The model is based on the Mesoscopic Distinct Element Method (MDEM), that computes the damped dynamics of a collection of interacting rigid bodies of given mass, shape and tensor of inertia – distinct elements. The state variables for each element include translational positions and velocities, as well as rotational positions and angular velocities. The bodies change their velocities and angular velocities due to contact forces and moments arising in pair interactions, as well as external forces and moments, acting at each body. The system is evolved in time with an explicit velocity Verlet time integration scheme.
Within our approach, undeformed CNTs are homogenized into cylindrical shells with finite thickness and then partitioned into identical segments of finite lengths (Fig. 1). Each distinct element represents the inertial properties of a CNT segment. Table 1 provides segment parameters for CNTs with diameter and length . Each segment contains approximately carbon atoms.
The inertial properties of a segment are represented with mass and the moment of inertia (a spherical tensor of inertia of a distinct element is assumed here for simplicity, however, a number of non-spherical shapes are immediately available in the framework). Parameters and are matched with the mass and moment of inertia of a cylindrical segment taken with respect to the CNT axis by proper choice of the radius of a distinct element:
The elements representing CNT segments are equispaced at a distance .
Such a coarse-graining leads to the elimination of atomistic degrees of freedom. Therefore, the dissipative microscopic processes associated with CNT sliding are not explicitly captured and should be included in a phenomenological manner. In this respect, two independent channels of energy dissipation are introduced in our model – local damping and viscous damping, controlled by the parameters and respectively.
PFC-style local damping group Inc. (2008) acts at each body. It is introduced to damp stiff interactions and stabilize time integration in dense particle assemblies. The components of damping force (moment ) are proportional to the corresponding components of unbalanced force (moment ) according to:
Here and are components of the translational and rotational velocity of an element, and is the sign function. In our simulations the local damping is kept relatively low () and its influence on CNT dynamics is insignificant.
Recent theoretical findings Brilliantov et al. (2017) predict a linear dependence between force and velocity in nanoscale thermal-induced friction. In our coarse-grained model we introduce the viscous damping that provides a simplified model for nanotribological energy losses during CNT sliding. The viscous damping forces, proportional to relative segment velocities, act in parallel with vdW contact forces. These forces are controlled by parameter . Normal and tangential viscous forces are proportional to normal and tangential relative velocities of elements in vdW contact:
Viscosity coefficients and are related to as follows:
where and are linearized stiffnesses of the contact model, taken , .
is kept as a free parameter in our simulations. Initial estimate ofis performed in a top-down way, as discussed in Ostanin et al. (2014). Unless otherwise noted is assumed in this work.
Elasticity of CNTs in our model is represented with the formalism of enhanced vector model (EVM) Kuzkin and Asonov (2012). Unlike incremental parallel bonds, widely used in DEM community Potyondy and Cundall (2004), EVM bonds explicitly conserve energy in symplectic time integration. EVM is based on a binding potential, describing the behavior of an elastic bond linking two rigid bodies. The formulation provides straightforward generalization for the case of large strains and accounts for a bending-twisting coupling.
Consider two equal-sized particles and with equilibrium separation and orientation described in terms of orthogonal vectors , as depicted in fig. 2(A). Then the EVM bond potential is given as follows:
Here is the radius vector connecting centers of bonded particles. Parameters are directly related to longitudinal, shear, bending, and torsional rigidities of a bond, according to Euler-Bernoulli beam theory (see Kuzkin and Asonov (2012) for more details):
where are area, moment of inertia and polar moment of inertia of a cylinder shell beam with radius and thickness :
Table 1 summarizes the elastic properties of CNT segments.
Van der Waals interactions between CNT segments are represented with the coarse-grained anisotropic vdW potential, as described in Ostanin et al. (2013, 2014). The potential is defined in terms of three variables , and , related to mutual position and orientation of two cylindrical segments (fig. 2(C)):
The potential describes interactions between CNT segments, taking into account their relative orientation, and providing symmetry for the interaction of two parallel infinite straight CNTs with respect to translation of one CNT along its axis. This property is crucially important for correct representation of shear interactions in CNT bundles and load transfer in CNT materials. It is important to note that such vdW potential does not create any barriers to CNT sliding, which corresponds to incommensurate mode of CNT slip. Such model does not create any static friction between CNTs and therefore will underestimate the strength of CNT materials in comparison with full atomistic model, where commensurate contacts between neighboring CNTs with nonzero static friction are possible.
The parameters of the anisotropic vdW potential for CNTs are provided in table 2.
Large-scale simulations of CNT assemblies are possible using the waLBerla framework. The development of the framework is focused on providing a highly parallel and highly optimized basis for multi-physics applications Godenschwager et al. (2013). It offers a rigid particle dynamics module which is capable of simulating up to particles with up to parallel processes Preclik and Rüde (2015). A detailed description of the algorithms and their implementation used in the framework can be found in Preclik and Rüde (2015); Schornbaum and Rüde (2016); Eibl and Rüde (2018). The framework is freely available under GPLv3 license at www.walberla.net. Here we only outline the basic features used for our simulations.
The simplified course of our simulation is presented in fig. 3. The simulation begins with the generation of the initial geometry of CNTs, and imposition of boundary conditions. Next, the simulation domain is partitioned using the distributed forest of octrees approach implemented in the waLBerla framework. Since no refinement is used all subdomains have exactly the same rectangular shape and are arranged on a regular three dimensional grid. Subsequently, the subdomains are distributed among the available processes in a balanced manner. Every process is now responsible for one or more subdomains. In the following, we will denote the numbers of subdomains in , and directions with a set of three integers . During the simulation each process stores the information of all rigid bodies that are within its associated subdomains. At the next stage, time integration cycles are performed on all MPI processes. The integration cycle consists of four parts. First, all particle contacts are detected. To reduce the number of expensive contact checks advanced algorithms based on hierarchical hash grids and bounding volumes are used. The contact detection scheme used for rigid body dynamics is adapted for potential-based interactions. This is done by decoupling the particle’s inertial radius, used to compute its moment of inertia, and its interaction radius, used in contact detection schemes. The latter in our case is equal to potential interaction cutoff radius. After all contacts are determined the contacts are used to calculate the forces and moments acting at each contact. For the calculation of the forces and moments the interaction model described above is used. In the next step the forces and moments are used to calculate accelerations and angular accelerations and the particle position and rotation gets updated using a velocity Verlet time integration scheme. Since the whole simulation is conducted in parallel information that is relevant for more than one processes has to be updated. Therefore the final step of the integration cycle is the synchronization. It not only takes care of synchronizing all information but also migrates particles between processes if the subdomain they belong to has changed. Throughout the whole simulation information about the current particle configuration is gathered on a regular interval and saved to log files for later analysis.
3 Numerical results
3.1 Simulations setup
Below we present the results of our numerical study on the mechanics of thin CNT films. Every test consists of two stages - i) self-assembly of the CNT film specimen, and ii) mechanical test on the specimen. Fig. 4 illustrates the geometry and boundary conditions used in our simulations. At self-assembly stage, initially separate, straight CNTs are deposited into a cuboid of the sizes with triple periodic boundary conditions (fig.4(A)), and then evolved to an equilibrated state. The computational domain is decomposed onto MPI blocks in manner. After equilibration, the obtained configuration is gathered to a master process, saved and then used in mechanical test simulations.
The mechanical test is performed in the following few steps.
First, we resize the computational domain to the size sufficient to fully contain the stretched specimen. For example, in uniaxial tension test, we resize domain borders to the size , where is tensile strain and is small margin.
Second, we replace the periodic boundary conditions with destructive boundary conditions in directions of stretching (-direction in uniaxial test and and directions in the biaxial test).
Then, we specify the new static subdivision into MPI blocks. For uniaxial tests we subdivide the domain into long blocks , oriented along the direction of stretching. For biaxial tests, we save the decomposition .
Next, we load the equilibrated configuration into a new MPI block configuration and scatter it over available MPI processes.
After that we specify grip regions and prescribe the velocities of grips. In a uniaxial tension test the grips are thin cuboid slices of the specimen (fig.4(B)), that are moving in opposite directions. They move with the constant velocity, except small constant acceleration period at the beginning of the simulation, introduced to avoid the initial inertial peak in the response:
In the case of biaxial tension test, the grips region is a margin surrounding the gage region, consisting of four grips (fig. 4(C)). The prescribed velocity field of the grip elements in this case is given by:
In the simulations described below the time is set to be 500 timesteps (10 ps).
In uniaxial test, the stress is measured as an arithmetic average of -projections of forces acting on two grips, related to a cross section area . For the case of biaxial tests four independently measured force projections are averaged to compute stress.
3.2 Self assembly of a CNT film specimen
Fig. 5 and fig. 6 summarize the initial relaxation of a CNT film specimen. Fig. 5 demonstrates the part of specimen at the beginning of the relaxation process (A) and after cycles ( ns) of relaxation (B,C). One can see that the initially straight CNTs are bent by vdW adhesive forces and most of the CNTs are joined into small bundles.
Fig. 6 demonstrates the evolution of vdW adhesion energy and elastic strain energy during the relaxation. One can see that the energy balance depends significantly on CNT length, whereas the size of the specimen almost does not affect it. The difference between random realizations appears to be negligible for large enough specimens.
The specimens generated for the mechanical tests described below had periodic out-of-plane boundary conditions, in-plane sizes from to and CNT lengths of , , and CNT segments. The same volumetric density ( was prescribed. Unless otherwise noted, the specimens have thickness of and surface density of .
3.3 Mechanical tests on CNT film specimen
Here we summarize the results of small-strain mechanical tests, performed as described above. The mechanics of a CNT film is characterized by the elastic deformation of individual CNTs and relative slip of the CNTs in the specimen. The resistance to slip is conditioned by change in surface energy and dissipation due to damping, introduced in our model to represent the effects of energy transfer to implicit (microscopic) degrees of freedom. Both local and viscous damping cause energy dissipation that is directly proportional to strain rate Ostanin et al. (2013, 2014). In our simulations local damping is kept low and has a negligible effect on the mechanics of CNT assembly. Overall mechanical response of the film is characterized by the initial viscoelastic deformation at small strains, conditioned by deformation of individual CNTs, followed by visco-plastic region, characterized by relative sliding of CNTs in the film. Due to viscous part of the response, the observed stress-strain curves start from nonzero stress.
Fig. 7 summarizes the mechanics of CNT films observed in uniaxial tests. The Young’s modulus of CNT films observed in our simulations varies in a large diapason of GPa, and conditioned by many factors, including strain rate (fig. 7(A)), CNT length (fig. 7(B)) and specimen size (fig. 7(C)). The dependence of the Young’s modulus on thickness is rather negligible due to out-of-plane periodic boundary conditions (fig. 7(D)).
The analysis of the velocity fields in CNT films during the mechanical test indicates that, given sufficient load transfer between CNTs, a CNT assembly deforms nearly uniformly, however, with some heterogeneity associated with random microstructure of the film. Relative slip of CNTs can be observed at large magnifications (fig. 8).
It appears that the dependence of the response on the size of the specimen is always associated with localization of strain close to grip regions. At large strain rates and low intertube sliding dissipation the specimen deformation appears to be nonuniform (fig. 9(A) gives the velocity field in the case of nonuniform, size-dependent deformation). Conversely, at slow enough strain rates and large enough dissipation the specimen deformation is uniform (fig. 9(B)).
Fig. 10(A-C) demonstrates the convergence to size-independent response at slow strain rates. It gives the comparison of stress-strain curves, obtained for two specimens of sizes and . At snapshots (A-C) the simulation is carried out at , and strain rate, respectively. The viscous dissipation coefficient is increased proportionally: it is , and , respectively. One can see that at slow enough strain rates the response becomes nearly size-independent which is associated with negligible role of inertia-dependent CNT rearrangement processes. It is therefore clear that strain localization and size-dependence of stress-strain curves is associated with inertial effects, that are negligible in slow and highly damped simulations.
It is interesting to compare the behavior of CNT films in uniaxial and biaxial tests. Figure 10(D) gives an example of such a comparison. It appears that the slope of stress-strain curve is slightly higher in a biaxial mechanical test. Assume validity of linear isotropic elastic relations for these initial slopes:
Fig. 10(D) indicates the values of and Therefore, it follows that the film behaves as linear in-plane isotropic elastic material with the Young’s modulus of GPa and Poisson’s ratio of . However, in simulations with long CNTs and low damping, where CNT sliding dominates the deformation, the observed Poisson’s ratio drops down to . The transitions between negative and positive Poisson’s ratios in thin CNT films depending on the specimen parameters has been observed experimentally Ma et al. (2010).
4 Related work and discussion
The mechanics of CNT films was studied with the toolkit of mesoscale modeling by a large community of authors. The pioneering works in the field are due to M. Buehler and co-workers Buehler (2006); Cranford and Buehler (2010). The use of bead-spring model allowed them to model large assemblies of CNTs and study stress-strain curves of a CNT material. Similar approaches were later used in a number of publications of other authors Xie et al. (2011); Li and Kröger (2012). However, due to the artifact of vdW potential corrugation, that is described in detail in Ostanin et al. (2013), such models are not capable to adequately represent shear interaction between CNTs, which led these authors to a conclusion that the only available mechanisms of deformation involve elastic deformations of individual CNTs, along with ”zipping and unzipping” of individual CNTs and bundles, resulting in histeretic, temperature-independent homogenized viscoelastic response Li and Kröger (2012). The other group of authors Volkov and Zhigilei (2010a, b, c); Wittmaack et al. (2018) combined symplectic finite temperature integration with ”slippery” vdW potential which resulted in nearly zero load transfer between separate CNTs and absent quasi-static mechanical response of a CNT material. This model assumptions were largely based on two observations: i) The CNT sliding observed in MD simulations does not feature significant viscous friction, and ii) given the absent static and dynamic friction between CNTs, the model predicts high mobility of CNTs, leading to the formation of large bundles, similar to ones observed in real CNT films and papers. The first argument is challenged by recent work Sinclair et al. (2018), that indicates that realistic MD simulation does indicate significant static and dynamic friction in graphitic systems. Regarding the second argument, it should be observed that the bundling of CNTs in real CNT film production occurs mostly on a gas/liquid stage and is driven by inelastic scattering and charge exchange processes, rather then by short-ranged vdW interactions. It should be also observed that the remaining rearrangements in CNT films, driven by vdW adhesion usually take hours and days, whereas the bundles formation in mesoscale simulations Volkov and Zhigilei (2010a, b, c); Wittmaack et al. (2018) take nanoseconds. Finally, it is clear that the real CNT films behave as a material, with the well-defined stress-strain curves, while frictionless simulations give zero quasi-static mechanical response.
In our work Ostanin et al. (2014) we have tried to model dissipation-based load transfer between CNTs, which led us to a conclusion, that, in spite of the absence of noticeable static shear interactions between CNTs, such systems demonstrate material-like response. In this work we conclude that the energy dissipation during CNT sliding does have the same effect in two-dimensional CNT systems.
Molecular dynamics simulations indicate that the incommensurate CNT tribology is well described by negligibly small static friction barrier and viscous energy dissipation during sliding. It is also important to note that the modern theoretical works Brilliantov et al. (2017) on nanofriction confirm linear dependence between friction force and relative sliding velocity in ”thermal” regime of friction between nanosurfaces. However, precise calibration of such dissipation is problematic, since the observed viscosity coefficient depends on multiple factors and require accurate statistic averaging. In this work we leave intertube viscosity as a free parameter – the calibration is performed in a top-down way, based on the macroscopic stress-strain curves.
Accurate bottom-up calibration of energy dissipation in CNT sliding would require the design of the mesoscale model for such dissipation calibrated with accurate ab initio and molecular dynamics simulations. Such model should also mimic static friction at commensurate CNT surface contacts.
5 Concluding remarks and future work
In this work we presented the computational study on the properties of thin CNT films under small strain deformations. Our study reveals the presence of a representative volume element in such structures and gives an idea on its size. Similarly to earlier studied CNT rope structures Ostanin et al. (2014), we observe an important role of energy dissipation, that is responsible for even load distribution along the specimen and its uniform deformation, leading to material-like mechanical response. The stress-strain curves observed in our simulations are qualitatively similar to the ones observed in real mechanical tests performed on CNT films Gubarev et al. (2019).
A number of questions remained beyond the scope of the present study. We provided the description of the mechanics for a narrow class of such films, produced by deposition of perfectly separate, length and type purified CNTs on a flat substrate, assembled and ordered only by short-range vdW interactions. In is clear, that realistic CNT assemblies have a lot more complicated morphology, conditioned not only by vdW adhesion on the deposition stage, but also kinetic interactions and aggregation on the aerosol/suspension stage.
The top-down calibration of energy dissipation in the model allowed us to qualitatively reproduce the behavior of realistic CNT film under mechanical loading. However, the most important deficiency of our model in its current shape is the absence of accurate model for the energy dissipation and its bottom-up calibration, that would allow us to reproduce the mechanics for CNT materials from the first principles. The development of such model, taking into account static friction for different CNT chiralities and commensurate/incommensurate contacts is the matter of future work.
I.O. acknowledges Russian Science Foundation (grant 17-73-10442) for the support of the development of the parallel system of modeling carbon nanotubes based on waLBerla framework. The financial support from the Russian Foundation for Basic Research under grant 18-29-19198 is deeply appreciated. The simulations presented in the paper were performed on computational clusters Pardus and Zhores (Skolkovo Institute of Science and Technology).
- Wu et al.  Zhuangchun Wu, Zhihong Chen, Xu Du, Jonathan M. Logan, Jennifer Sippel, Maria Nikolou, Katalin Kamaras, John R. Reynolds, David B. Tanner, Arthur F. Hebard, and Andrew G. Rinzler. Transparent, conductive carbon nanotube films. Science, 305(5688):1273–1276, 2004. ISSN 0036-8075. doi: 10.1126/science.1101243. URL http://science.sciencemag.org/content/305/5688/1273.
- Park et al.  Steve Park, Michael Vosguerichian, and Zhenan Bao. A review of fabrication and applications of carbon nanotube film-based flexible electronics. Nanoscale, 5:1727–1752, 2013. doi: 10.1039/C3NR33560G. URL http://dx.doi.org/10.1039/C3NR33560G.
- Gubarev et al.  V. M. Gubarev, V. Y. Yakovlev, M. G. Sertsu, O. F. Yakushev, V. M. Krivtsun, V. V. Medvedev, Yu. G. Gladush, I. A. Ostanin, A. Sokolov, F. Schafers, and A. G. Nasibulin. Single-walled carbon nanotube membranes for optical applications in the extreme ultraviolet range. Nanoscale, X(X):XX–XX, 2019. ISSN 0036-8075.
- Wang and Loh  Long Wang and Kenneth J Loh. Wearable carbon nanotube-based fabric sensors for monitoring human physiological performance. Smart Materials and Structures, 26(5):055018, 2017. URL http://stacks.iop.org/0964-1726/26/i=5/a=055018.
- Yakobson et al.  B. I. Yakobson, C. J. Brabec, and J. Bernholc. Nanomechanics of carbon tubes: Instabilities beyond linear response. Phys. Rev. Lett., 76:2511–2514, Apr 1996. doi: 10.1103/PhysRevLett.76.2511. URL https://link.aps.org/doi/10.1103/PhysRevLett.76.2511.
- Dumitrica et al.  Traian Dumitrica, Ming Hua, and Boris I. Yakobson. Symmetry-, time-, and temperature-dependent strength of carbon nanotubes. Proceedings of the National Academy of Sciences, 103(16):6105–6109, 2006. doi: 10.1073/pnas.0600945103. URL https://www.pnas.org/content/103/16/6105.
- Zhang and Dumitrică  D.-B. Zhang and T. Dumitrică. Elasticity of ideal single-walled carbon nanotubes via symmetry-adapted tight-binding objective modeling. Applied Physics Letters, 93(3):031919, 2008. doi: 10.1063/1.2965465. URL https://doi.org/10.1063/1.2965465.
- Nikiforov et al.  I. Nikiforov, D.-B. Zhang, R. D. James, and T. Dumitrică. Wavelike rippling in multiwalled carbon nanotubes under pure bending. Applied Physics Letters, 96(12):123107, 2010. doi: 10.1063/1.3368703. URL https://doi.org/10.1063/1.3368703.
- Cornwell and Welch  Charles F. Cornwell and Charles R. Welch. Very-high-strength (60-gpa) carbon nanotube fiber design based on molecular dynamics simulations. The Journal of Chemical Physics, 134(20):204708, 2011. doi: 10.1063/1.3594197. URL https://doi.org/10.1063/1.3594197.
- Buehler  Markus J. Buehler. Mesoscale modeling of mechanics of carbon nanotubes: Self-assembly, self-folding, and fracture. Journal of Materials Research, 21(11):2855–2869, 2006. doi: 10.1557/jmr.2006.0347.
- Cranford and Buehler  Steven W Cranford and Markus J Buehler. In silico assembly and nanomechanical characterization of carbon nanotube buckypaper. Nanotechnology, 21(26):265706, 2010. URL http://stacks.iop.org/0957-4484/21/i=26/a=265706.
- Mirzaeifar et al.  Reza Mirzaeifar, Zhao Qin, and Markus J. Buehler. Mesoscale mechanics of twisting carbon nanotube yarns. Nanoscale, 7:5435–5445, 2015. doi: 10.1039/C4NR06669C. URL http://dx.doi.org/10.1039/C4NR06669C.
- Volkov and Zhigilei [2010a] Alexey N. Volkov and Leonid V. Zhigilei. Structural stability of carbon nanotube films: The role of bending buckling. ACS Nano, 4(10):6187–6195, 2010a. doi: 10.1021/nn1015902. URL https://doi.org/10.1021/nn1015902. PMID: 20931973.
- Volkov and Zhigilei [2010b] Alexey N. Volkov and Leonid V. Zhigilei. Mesoscopic interaction potential for carbon nanotubes of arbitrary length and orientation. The Journal of Physical Chemistry C, 114(12):5513–5531, 2010b. doi: 10.1021/jp906142h. URL https://doi.org/10.1021/jp906142h.
- Volkov and Zhigilei [2010c] Alexey N. Volkov and Leonid V. Zhigilei. Scaling laws and mesoscopic modeling of thermal conductivity in carbon nanotube materials. Phys. Rev. Lett., 104:215902, May 2010c. doi: 10.1103/PhysRevLett.104.215902. URL https://link.aps.org/doi/10.1103/PhysRevLett.104.215902.
- Wittmaack et al.  Bernard K. Wittmaack, Alexey N. Volkov, and Leonid V. Zhigilei. Mesoscopic modeling of the uniaxial compression and recovery of vertically aligned carbon nanotube forests. Composites Science and Technology, 166:66 – 85, 2018. ISSN 0266-3538. doi: https://doi.org/10.1016/j.compscitech.2018.03.014. URL http://www.sciencedirect.com/science/article/pii/S0266353817330683. Carbon nanotube composites for structural applications.
- Ostanin et al.  Igor Ostanin, Roberto Ballarini, David Potyondy, and Traian Dumitrica. A distinct element method for large scale simulations of carbon nanotube assemblies. Journal of the Mechanics and Physics of Solids, 61(3):762 – 782, 2013. ISSN 0022-5096. doi: http://dx.doi.org/10.1016/j.jmps.2012.10.016. URL http://www.sciencedirect.com/science/article/pii/S0022509612002372.
- Ostanin et al.  Igor Ostanin, Roberto Ballarini, and Traian Dumitrica. Distinct element method modeling of carbon nanotube bundles with intertube sliding and dissipation. Journal of Applied Mechanics, 81(6):06004, 2014. ISSN 0021-8936. doi: 10.1115/1.4026484. URL http://dx.doi.org/10.1115/1.4026484.
- Ostanin et al.  Igor Ostanin, Roberto Ballarini, and Traian Dumitrica. Distinct element method for multiscale modeling of cross-linked carbon nanotube bundles: From soft to strong nanomaterials. Journal of Materials Research, xxx(xx):JMR-2014-0365, 2015. Accepted manuscript. To appear in 2015 Focus Issue on Soft Nanomaterials.
- Cundall and Strack  P. A. Cundall and O.D. Strack. A discrete numerical model for granular assemblies. Geotechnique, 29:47–65, 1979.
- Wang et al.  Yuezhou Wang, Cristian Gaidau, Igor Ostanin, and Traian Dumitrica. Ring windings from single-wall carbon nanotubes: A distinct element method study. Applied Physics Letters, 103(18):183902, 2013. doi: http://dx.doi.org/10.1063/1.4827337. URL http://scitation.aip.org/content/aip/journal/apl/103/18/10.1063/1.4827337.
- Wang et al.  Yuezhou Wang, Matthew R. Semler, Igor Ostanin, Erik K. Hobbie, and Traian Dumitrica. Rings and rackets from single-wall carbon nanotubes: manifestations of mesoscale mechanics. Soft Matter, 10:8635–8640, 2014. doi: 10.1039/C4SM00865K. URL http://dx.doi.org/10.1039/C4SM00865K.
- Wang et al.  Yuezhou Wang, Hao Xu, Grigorii Drozdov, and Traian Dumitrică. Mesoscopic friction and network morphology control the mechanics and processing of carbon nanotube yarns. Carbon, 139:94 – 104, 2018. ISSN 0008-6223. doi: https://doi.org/10.1016/j.carbon.2018.06.043. URL http://www.sciencedirect.com/science/article/pii/S0008622318306055.
- Ostanin et al.  I. Ostanin, P. Zhilyaev, V. Petrov, T. Dumitrica, S. Eibl, U. Ruede, and V. Kuzkin. Toward large scale modeling of carbon nanotube systems with the mesoscopic distinct element method. Letters on Materials, 8(3):240–245, 2018. doi: 10.22226/2410-3535-2018-3-240-245. URL https://lettersonmaterials.com/en/Readers/Article.aspx?aid=8700.
- group Inc.  Itasca Consulting group Inc. PFC3D (Particle Flow Code in Three Dimensions). Version 4.0. Itasca Consulting group Inc., 2008.
- Brilliantov et al.  N. V. Brilliantov, Yu. A. Budkov, and C. Seidel. Theoretical and numerical analysis of nano-actuators based on grafted polyelectrolytes in an electric field. Faraday Discuss., 199:487–510, 2017. doi: 10.1039/C6FD00240D. URL http://dx.doi.org/10.1039/C6FD00240D.
- Kuzkin and Asonov  Vitaly A. Kuzkin and Igor E. Asonov. Vector-based model of elastic bonds for simulation of granular solids. Phys. Rev. E, 86:051301, Nov 2012. doi: 10.1103/PhysRevE.86.051301. URL https://link.aps.org/doi/10.1103/PhysRevE.86.051301.
- Potyondy and Cundall  D.O. Potyondy and P.A. Cundall. A bonded-particle model for rock. International Journal of Rock Mechanics and Mining Sciences, 41(8):1329 – 1364, 2004. ISSN 1365-1609. doi: https://doi.org/10.1016/j.ijrmms.2004.09.011. URL http://www.sciencedirect.com/science/article/pii/S1365160904002874. Rock Mechanics Results from the Underground Research Laboratory, Canada.
- Godenschwager et al.  Christian Godenschwager, Florian Schornbaum, Martin Bauer, Harald Köstler, and Ulrich Rüde. A framework for hybrid parallel flow simulations with a trillion cells in complex geometries. In Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, page 35. ACM, ACM Press, 2013. doi: 10.1145/2503210.2503273.
- Preclik and Rüde  Tobias Preclik and Ulrich Rüde. Ultrascale simulations of non-smooth granular dynamics. Computational Particle Mechanics, 2(2):173–196, May 2015. doi: 10.1007/s40571-015-0047-6.
- Schornbaum and Rüde  Florian Schornbaum and Ulrich Rüde. Massively parallel algorithms for the lattice boltzmann method on NonUniform grids. SIAM Journal on Scientific Computing, 38(2):C96–C126, January 2016. doi: 10.1137/15M1035240.
- Eibl and Rüde  Sebastian Eibl and Ulrich Rüde. A local parallel communication algorithm for polydisperse rigid body dynamics. Parallel Computing, 80:36–48, December 2018. doi: 10.1016/j.parco.2018.10.002.
- Ma et al.  Yin Ji Ma, Xue Feng Yao, Quan Shui Zheng, Ya Jun Yin, Dong Jie Jiang, Guang Hui Xu, Fei Wei, and Qiang Zhang. Carbon nanotube films change poisson’s ratios from negative to positive. Applied Physics Letters, 97(6):061909, 2010. doi: 10.1063/1.3479393. URL https://doi.org/10.1063/1.3479393.
- Xie et al.  Bo Xie, Yilun Liu, Yiting Ding, Quanshui Zheng, and Zhiping Xu. Mechanics of carbon nanotube networks: microstructural evolution and optimal design. Soft Matter, 7:10039–10047, 2011. doi: 10.1039/C1SM06034A. URL http://dx.doi.org/10.1039/C1SM06034A.
- Li and Kröger  Ying Li and Martin Kröger. Viscoelasticity of carbon nanotube buckypaper: zipping–unzipping mechanism and entanglement effects. Soft Matter, 8:7822–7830, 2012. doi: 10.1039/C2SM25561H. URL http://dx.doi.org/10.1039/C2SM25561H.
- Sinclair et al.  Robert C. Sinclair, James L. Suter, and Peter V. Coveney. Graphene–graphene interactions: Friction, superlubricity, and exfoliation. Advanced Materials, 30(13):1705791, 2018. doi: 10.1002/adma.201705791. URL https://onlinelibrary.wiley.com/doi/abs/10.1002/adma.201705791.