🔬 physics

Vibrational Dynamics of Fixed-Free Mass-Spring Chains: Insights into Longitudinal Oscillations of Carbon Nanotubes (CNTs)

This study demonstrates that while the longitudinal vibrational displacement of carbon nanotubes remains linear and amplitude-independent under fixed-free conditions, the emergence of harmonics in axial force variance serves as an early indicator of nonlinearity, validating a combined displacement-force analysis for characterizing nanoscale dynamics and developing NEMS.

Original authors: Ayouba BATOURE, Alio ISSOUFOU ARZIKA

Published 2026-08-04

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Ayouba BATOURE, Alio ISSOUFOU ARZIKA

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the tiniest, thinnest string you can possibly picture, made not of cotton or nylon, but of carbon atoms woven into a perfect, hollow tube. This is a carbon nanotube, a superstar of the nanoworld. Scientists are obsessed with these tubes because they might be the future of super-fast computers, ultra-sensitive sensors, and microscopic machines called NEMS (Nano-Electro-Mechanical Systems). But to make these machines work, we need to understand how they vibrate. Think of a guitar string: when you pluck it, it wiggles back and forth. If you hold one end tight and let the other end dangle free, the way it wiggles is very specific. Now, imagine that guitar string is so small that the "atoms" inside it are the individual beads of the string. How do those beads move? Do they all wiggle the same amount, or does the free end dance more wildly than the held end? This is the puzzle of vibrational dynamics. It's a mix of simple physics (like springs connecting balls) and complex reality (where atoms bump and push each other in tricky ways). Understanding this helps engineers build better, more reliable nano-machines that don't break or behave strangely when they start humming.

In this study, two researchers, Ayouba Batouré and Alio Issoufou Arzika, decided to solve this puzzle by playing a game of "compare and contrast." They wanted to see if a super-simple model could predict how a real, complex carbon nanotube behaves. First, they built a theoretical model: a line of identical balls connected by perfect springs, with one end glued to a wall and the other end free to swing. They did the math to see how these balls would wiggle. Then, they used a powerful computer program called LAMMPS to simulate a real carbon nanotube (specifically an armchair (10,10) type) made of 1,640 carbon atoms. They held one end fixed and shook the other end with tiny, controlled pushes, measuring how the atoms moved and how much force they exerted on each other.

Here is what they found, and it's a bit like discovering that a simple toy car can predict how a real race car handles a turn, but only up to a point.

The "Springy" Truth: The Free End Dances More
Both the simple spring model and the complex computer simulation agreed on one big thing: the location of the vibration matters. When you hold one end of the chain (or tube) still, the atoms near the fixed end barely move. But as you get closer to the free end, the wiggles get bigger and bigger. The researchers found that the "variance" (a fancy word for how much the atoms jitter around their spot) grows in a nearly straight line from the fixed wall to the free tip. It's like a line of people holding hands; if the person at the front is glued to a wall, the person at the back gets to swing their arms the widest. This happens whether the system is a simple line of springs or a real carbon nanotube. This suggests that the way the tube is held is the most important factor in how it vibrates.

The "Force" Surprise: The Hidden Nonlinearity
This is where things get interesting. The researchers shook the nanotube with different strengths, ranging from a tiny 0.01 Ångström to a slightly larger 0.10 Ångström (an Ångström is one ten-billionth of a meter, so these are incredibly tiny movements).

They looked at two things: how far the atoms moved (displacement) and how hard they pushed against each other (axial force).

The Movement: The atoms' movement stayed very "clean." Even when they shook it harder, the atoms just wiggled faster or a bit more, but they didn't change their rhythm. The main frequency stayed at 6.248 GHz, and the pattern of the wiggles didn't change shape. It was like a swing that goes higher but still takes the exact same amount of time to go back and forth.
The Force: However, the force inside the tube told a different story. As they shook it harder, the force signal started to get "messy." New frequencies, called harmonics, appeared. It's like if you plucked a guitar string gently, it sounds pure, but if you strum it hard, you start to hear a slightly distorted, buzzy sound. The computer simulation showed that the internal stresses (the forces between atoms) were starting to act non-linearly, meaning the atoms were pushing back in a way that wasn't perfectly proportional to how far they were pulled.

The "Softening" Effect
The researchers also calculated the "stiffness" of the tube, which they expressed as Young's Modulus. At the tiniest shake (0.01 Å), the tube was very stiff, with a value of 1.02 TPa (Terapascals). But as they shook it harder, the tube seemed to get softer. At 0.05 Å, the value dropped to 0.78 TPa, and at 0.10 Å, it plummeted to 0.40 TPa. This suggests that while the tube is incredibly strong, if you push it too hard, it starts to lose some of its rigidity. The paper notes that this "softening" is a sign that the atoms are starting to interact in more complex, non-linear ways that the simple spring model doesn't capture.

What This Means for the Future
The big takeaway is that the simple "balls and springs" model is actually a pretty good teacher. It correctly predicted that the free end of the tube would wiggle the most and that the vibrations would follow a specific pattern based on how the tube is held. This means scientists can use these simple, cheap-to-run models to get a good first guess about how nano-machines will behave.

However, the paper also warns that if you want to know exactly how the forces inside the tube behave when things get a bit more intense, you need the complex computer simulations. The simple model misses the "buzzy" harmonics and the softening effect that appear in the real atomic world. The researchers suggest that for building future nano-sensors and resonators, it's best to use the simple model to understand the basics, but rely on the detailed simulations to catch the subtle, non-linear tricks that happen when the vibrations get a little too energetic. They didn't prove that this solves all problems, but they showed a clear path: use the simple model for the big picture, and the complex one for the fine print.

Technical Summary: Vibrational Dynamics of Fixed-Free Mass-Spring Chains and Carbon Nanotubes

Problem Statement
The longitudinal vibrational dynamics of one-dimensional nanostructures, particularly carbon nanotubes (CNTs), remain less explored than their flexural modes, despite their potential applications in nano-resonators and nano-electromechanical systems (NEMS). While simplified mass-spring models have historically provided conceptual frameworks for wave propagation, there is a need to rigorously characterize the spatial and spectral organization of vibrations under fixed-free boundary conditions and to distinguish between boundary-induced effects and genuine atomistic nonlinearities in real nanostructures.

Methodology
The study employs a multi-scale approach combining analytical modeling with classical molecular dynamics (MD) simulations:

Analytical Model: A one-dimensional chain of $N$ identical point masses ( $m$ ) interconnected by linear springs (stiffness $K$ ) is analyzed under fixed-free boundary conditions (left end fixed, right end free). The system is solved via diagonalization of the dynamical matrix to identify normal modes, eigenfrequencies, and eigenvectors. The spatial distribution of displacement variance is derived analytically.
Molecular Dynamics Simulations: MD simulations were performed on a single-walled armchair (10,10) CNT (1640 atoms, radius $\approx 6.78$ Å) using the LAMMPS software. The interatomic interactions were modeled using the AIREBO potential, which accounts for covalent bonds, van der Waals forces, and torsion, allowing for the capture of anharmonicities.
Simulation Protocol: The CNT was configured with a closed, fixed end and an open, free end. After energy minimization and thermalization at 300 K, controlled mechanical excitations were applied to the free end. The system evolved in the microcanonical ensemble (NVE) to preserve intrinsic nonlinear effects without artificial dissipation.
Analysis: Trajectories were analyzed for spatial variance of displacements and axial forces. Spectral analysis (FFT) was conducted on displacement and force signals across excitation amplitudes ranging from 0.01 to 0.10 Å. Young's modulus was estimated via linear regression of force versus elongation in the small deformation regime.

Key Results

Spatial Variance in Harmonic Chains: Analytical and numerical results for the mass-spring chain demonstrate that the time-averaged squared displacement variance ( $\langle \delta u^2_n \rangle$ ) increases nearly linearly from the fixed end to the free end. This growth is a robust property of fixed-free boundary conditions, independent of excitation intensity in the linear regime. The spatial profile of fluctuations remains invariant, scaling only with the amplitude of excitation.
Energy Conservation: In the harmonic model, total mechanical energy is strictly conserved, with continuous exchange between kinetic and potential forms. In the CNT simulations, total energy conservation was excellent, validating the NVE integration.
CNT Vibrational Dynamics:
- Displacement: The spatial variance of axial displacement in the (10,10) CNT mirrors the linear growth observed in the harmonic chain. Crucially, this variance remains practically independent of the excitation amplitude (0.01–0.10 Å), and the displacement spectrum is dominated by the fundamental frequency (6.248 GHz) with negligible harmonics.
- Axial Force: In contrast to displacement, the variance of axial force exhibits a significant dependence on excitation amplitude. Spectral analysis reveals the appearance of harmonics (at 12.497 GHz and 18.745 GHz) in the force signal as amplitude increases, indicating that internal stresses capture nonlinear effects earlier and more distinctly than displacement.
- Young's Modulus: The effective Young's modulus decreases with increasing excitation amplitude (from 1.02 TPa at 0.01 Å to 0.40 TPa at 0.10 Å), reflecting a softening of the mechanical response due to anharmonic contributions in the AIREBO potential.
Modal Structure: Low-frequency modes dominate collective displacements, while high-frequency modes contribute more to energy distribution. The spatial growth of variance is directly linked to the modal decomposition of the system.

Significance and Claims
The paper claims that the joint analysis of displacements and forces provides an effective framework for characterizing early nonlinear effects in CNTs. The primary contribution is the demonstration that the fixed-free mass-spring model faithfully reproduces the main dynamic properties observed at the atomic scale, specifically the spatial distribution of fluctuations and the influence of boundary conditions.

The authors assert that the first signatures of nonlinearity in CNTs appear in internal stresses (manifested as harmonics in the force spectrum and amplitude-dependent softening of the modulus) before significantly altering the displacement dynamics. Consequently, the simple harmonic mass-spring model serves as a relevant, low-cost theoretical tool for interpreting longitudinal vibrations in one-dimensional nanostructures and offers a conceptual basis for understanding the role of boundary conditions in fluctuation propagation and energy distribution. This work supports the development of NEMS by clarifying the distinction between boundary-induced spatial asymmetry and intrinsic atomistic nonlinearities.

More like this