Buckling behavior of carbon nanotubes under bending: from ripple to kink

This is the author accepted manuscript. The final version is available from Elsevier via the DOI in this record.This paper elucidates the buckling behavior of carbon nanotubes (CNTs) under bending. CNTs are modeled as continuous thin-wall circular tubes, and their buckling is governed by equations that take into account of the sectional Brazier effect and non-uniform structural deformation. The CNT governing equations (fourth-order ordinary differential nonlinear equations with integral conditions) are solved by introducing a continuation algorithm. In addition, the buckling behavior of CNTs under bending is simulated with objective molecular dynamics (OMD), and is used to verify the continuum results. We show that there exist low- and high-strain phases during the bending process of CNTs, and the transition in between may divide the whole bending process into three stages: low-curvature stage, mixed-curvature stage and high-curvature stage. Ripples are generated on the CNT surfaces before the formation of kinks. Compared to single-walled CNTs (SWCNTs), hydrogen-filled CNTs have a longer mixed-strain stage owing to the presence of internal pressure, and are therefore more inclined to exhibit a ripple morphology. Our results offer better understanding of the buckling behavior of CNTs, and may open up new opportunities for the design and applications of novel CNT-based nanoelectronics.National Natural Science Foundation of ChinaProgram for New Century 24 Excellent Talents in UniversityNatural Science Foundation of Heilongjiang Province of ChinaFundamental Research Funds for the Central Universitie

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives

Consequences of local gauge symmetry in empirical tight-binding theory

A method for incorporating electromagnetic fields into empirical tight-binding theory is derived from the principle of local gauge symmetry. Gauge invariance is shown to be incompatible with empirical tight-binding theory unless a representation exists in which the coordinate operator is diagonal. The present approach takes this basis as fundamental and uses group theory to construct symmetrized linear combinations of discrete coordinate eigenkets. This produces orthogonal atomic-like "orbitals" that may be used as a tight-binding basis. The coordinate matrix in the latter basis includes intra-atomic matrix elements between different orbitals on the same atom. Lattice gauge theory is then used to define discrete electromagnetic fields and their interaction with electrons. Local gauge symmetry is shown to impose strong restrictions limiting the range of the Hamiltonian in the coordinate basis. The theory is applied to the semiconductors Ge and Si, for which it is shown that a basis of 15 orbitals per atom provides a satisfactory description of the valence bands and the lowest conduction bands. Calculations of the dielectric function demonstrate that this model yields an accurate joint density of states, but underestimates the oscillator strength by about 20% in comparison to a nonlocal empirical pseudopotential calculation.Comment: 23 pages, 7 figures, RevTeX4; submitted to Phys. Rev.

Density functional theory calculations of the carbon ELNES of small diameter armchair and zigzag nanotubes: core-hole, curvature and momentum transfer orientation effects

We perform density functional theory calculations on a series of armchair and zigzag nanotubes of diameters less than 1nm using the all-electron Full-Potential(-Linearised)-Augmented-Plane-Wave (FPLAPW) method. Emphasis is laid on the effects of curvature, the electron beam orientation and the inclusion of the core-hole on the carbon electron energy loss K-edge. The electron energy loss near-edge spectra of all the studied tubes show strong curvature effects compared to that of flat graphene. The curvature induced $\pi-\sigma$ hybridisation is shown to have a more drastic effect on the electronic properties of zigzag tubes than on those of armchair tubes. We show that the core-hole effect must be accounted for in order to correctly reproduce electron energy loss measurements. We also find that, the energy loss near edge spectra of these carbon systems are dominantly dipole selected and that they can be expressed simply as a proportionality with the local momentum projected density of states, thus portraying the weak energy dependence of the transition matrix elements. Compared to graphite, the ELNES of carbon nanotubes show a reduced anisotropy.Comment: 25 pages, 15 figures, revtex4 submitted for publication to Phys. Rev.

Origin of Axial and Radial Expansions in Carbon Nanotubes Revealed by Ultrafast Diffraction and Spectroscopy

The coupling between electronic and nuclear degrees of freedom in low-dimensional, nanoscale systems plays a fundamental role in shaping many of their properties. Here, we report the disentanglement of axial and radial expansions of carbon nanotubes, and the direct role of electronic and vibrational excitations in determining such expansions. With subpicosecond and subpicometer resolutions, structural dynamics were explored by monitoring changes of the electron diffraction following an ultrafast optical excitation, whereas the transient behavior of the charge distribution was probed by time-resolved, electron-energy-loss spectroscopy. Our experimental results, and supporting density functional theory calculations, indicate that a population of the excited carriers in the antibonding orbitals of the nanotube walls drives a transient axial deformation in ∼1 ps; this deformation relaxes on a much longer time scale, 17 ps, by nonradiative decay. The electron-driven expansion is distinct from the phonon-driven dynamics observed along the radial direction, using the characteristic Bragg reflections; it occurs in 5 ps. These findings reveal the nonequilibrium distortion of the unit cell at early times and the role of the electron(phonon)-induced stress in the lattice dynamics of one-dimensional nanostructures

Formation of helices in graphene nanoribbons under torsion

We use objective boundary conditions and self-consistent charge density-functional-based tight-binding to simulate at the atomistic scale the formation of helices in narrow graphene nanoribbons with armchair edges terminated with fluorine and hydrogen. We interpret the microscopic data using an inextensible, unshearable elastic rod model, which considers both bending and torsional strains. When fitted to the atomistic data, the simple rod model uses closed-form solutions for a cubic equation to predict the strain energy and morphology at a given twist angle and the crossover point between pure torsion and a helix. Our modeling and simulation bring key insights into the origin of the helical graphene morphologies stored inside of carbon nanotubes. They can be useful for designing chiral nanoribbons with tailored properties