---
title: 'Graphene Helicoids: Mechanics, Thermal & Electronic Effects'
url: https://www.emergentmind.com/topics/graphene-helicoids
type: topic
---

# Graphene Helicoids: Mechanics, Thermal & Electronic Effects

Searching arXiv for relevant papers on graphene helicoids and helical graphene systems.
Graphene helicoids are graphene-based structures in which the material is organized with helicoidal or helical screw symmetry rather than as a flat sheet or a simple multilayer stack. In the literature surveyed here, the term covers at least two closely related regimes: a graphene nanoribbon wrapped into a screw-dislocation-like helicoid, often denoted GH, and a helical multilayer moiré system exemplified by helical trilayer graphene (HTG) with twist configuration $(\theta_1,\theta_2,\theta_3)=(\theta,0,-\theta)$ [1709.08329], [1808.02978], [2305.03031]. Across these regimes, the helicoidal geometry is not merely morphological. It governs tensile response, heat transport, and Dirac-carrier dynamics, and in the multilayer moiré setting it can self-organize into locally periodic $C_{2z}$-broken domains with topological flat bands and near-ideal quantum geometry [1610.03742], [1610.03742], [2305.03031].

## 1. Structural definitions and geometric realizations

In nanoscale mechanical and thermal studies, a graphene helicoid is constructed from a graphene nanoribbon arranged into a helicoid topology via a screw-dislocation geometry [1709.08329]. One formulation uses a single screw dislocation with Burgers vector magnitude $|b| = 3.4\ \text{Å}$, together with two graphene monolayers at the two ends to reduce edge effects [1709.08329]. The geometry is described by the outer radius $R$, inner radius $r$, and turn number $N$, with width
\[
w = R - r
\]
and total height
\[
h_{\text{tot} = (N+2)|b|
\]
in the molecular-dynamics nanospring model; the deformable region has
\[
N_{\text{eff} = N - 2
\]
turns [1709.08329].

A closely related thermal-transport construction also treats GH as a graphene nanoribbon wrapped into a continuous helical, screw-dislocation-like shape, with
\[
w = R - r, \qquad L_{\text{tot} = N \times |b|
\]
and a Burgers vector spacing of $|b| = 3.35\ \text{Å}$ [1808.02978]. In that setting, GH is contrasted with multilayer graphene (MLG): instead of discrete sheets held together only by van der Waals interactions, the helicoid is a continuous spiral with strong geometric continuity along its axis [1808.02978].

In continuum electronic treatments, the helicoid is a minimal surface. One helicoidal nanoribbon parameterization is
\[
\vec r = x\,\vec e_x + \xi\,[\cos(\omega x)\,\vec e_y + \sin(\omega x)\,\vec e_z],
\]
with twist rate
\[
\omega = \frac{2\pi n}{L},
\]
and induced metric
\[
ds^2 = (1+\omega^2 \xi^2)\,dx^2 + d\xi^2
\]
[1610.03742]. A related helicoid parameterization used for massless Dirac particles is
\[
\begin{pmatrix} x\\y\\z \end{pmatrix} = a\begin{pmatrix} \sinh u \cos v\\ \sinh u \sin v\\ v \end{pmatrix},
\]
with induced metric
\[
g_{ij}= \begin{pmatrix} a^2\cosh^2 u & 0\\ 0 & a^2\cosh^2 u \end{pmatrix}
\]
[1505.05159]. In both formulations the helicoid is emphasized as a minimal surface, so the electronic effects are attributed to geometry itself rather than to ordinary strain [1610.03742], [1505.05159].

In moiré graphene, the helical concept is extended to a three-layer stack with distributed twist. HTG is defined by
\[
(\theta_1,\theta_2,\theta_3)=(\theta,0,-\theta),
\]
with a special angle
\[
\theta \approx 1.5^\circ
\]
[2305.03031]. Unlike ordinary bilayer moiré systems, this unrelaxed structure contains two coupled moiré patterns and is therefore a supermoiré, or moiré-of-moiré, structure [2305.03031].

## 2. Mechanical response and nanospring behavior

The GH nanospring literature identifies a tensile response that is governed by delamination and subsequent ribbon stretching rather than by the mechanics of a conventional helical wire spring [1709.08329]. The deformation comprises four stages overall: initial delamination, stable delamination, elastic deformation, and failure [1709.08329]. The three elastic stages are central to its characterization as a nanospring.

During initial delamination, the force and strain energy increase sharply as adjacent turns begin to separate, overcoming interlayer van der Waals attraction. The force reaches a local maximum
\[
F_a \approx 8\ \text{nN},
\]
which is the critical force for initiating stable delamination [1709.08329]. In the stable-delamination stage, the force drops to
\[
F_c \approx 0.56\ \text{nN}
\]
and remains nearly constant with fluctuations, while the strain energy increases approximately linearly with strain [1709.08329]. After full delamination, the system behaves like a stretched graphene nanoribbon. The force rises again nonlinearly, the strain energy follows a parabolic trend, and the elastic-limit force is
\[
F_e \approx 6.04\ \text{nN},
\]
marking the onset of failure [1709.08329].

The reported tensile deformation capability is unusually large. The paper states that the yield strain can exceed $1000\%$, with a representative case approaching $1500\%$ [1709.08329]. The geometry dependence is explicit. With initial turn height
\[
h_0 = 3.4\ \text{Å}
\]
and fully delaminated height
\[
h = \sqrt{h_0^2 + 36 s_0^2},
\]
the delamination strain is estimated as
\[
\varepsilon_c = \frac{h - h_0}{h_0}.
\]
For a GH with eight effective turns and $s_0 \approx 7.1\ \text{Å}$, the predicted value is
\[
\varepsilon_c \approx 1157\%,
\]
in agreement with simulation [1709.08329]. The delamination strain is independent of turn number $N$ and outer radius $R$, and increases with increasing inner radius $r$ [1709.08329].

A central finding is that the elastic-limit force remains nearly constant across geometries, around $6\ \text{nN}$ [1709.08329]. The paper attributes this to highly non-uniform strain localization: inner edge atoms experience the largest tensile strain, atoms in the ribbon middle may experience slight compressive strain, and failure is governed by fracture of the graphene nanoribbon rather than by failure of the helicoid topology itself [1709.08329]. Crack initiation occurs at the armchair edge location, propagation follows the zigzag direction, and monoatomic chains and pentagon carbon rings form near fracture [1709.08329].

Hydrogen saturation changes these behaviors only modestly in force but more strongly in strain. The reported changes are a decrease of about $10\%$ in $F_a$, a decrease of about $7\%$ in $F_e$, and a slight increase of $F_c$ from $0.56$ to $0.66\ \text{nN}$; the strains decrease from about $13\%$ to $8\%$ for $\varepsilon_a$, from about $1137\%$ to $900\%$ for $\varepsilon_c$, and from about $1517\%$ to $1193\%$ for $\varepsilon_e$ [1709.08329]. This suggests that edge chemistry modifies stretchability more strongly than the characteristic force scales.

## 3. Thermal transport and strain-dependent heat conduction

Graphene helicoids have also been studied as thermal-transport systems whose axial conduction differs qualitatively from the through-plane response of multilayer graphene [1808.02978]. The key structural distinction is continuity: GH is a single continuous helical graphene ribbon, so axial heat flow contains both interlayer van der Waals transport and projected in-plane graphene conduction along the helical axis [1808.02978]. By contrast, MLG is composed of separate graphene sheets whose through-plane transport is governed mainly by weak van der Waals coupling, and whose cross-plane thermal conductivity saturates with increasing thickness [1808.02978].

The simulation framework adopts the Fourier-law-style definition
\[
k = -\frac{J}{\nabla T}, \qquad J = \frac{Q}{S},
\]
with
\[
\nabla T = \frac{\Delta T}{\Delta L}
\]
and $\Delta T = 20\ \text{K}$ [1808.02978]. Over the simulated range, GH exhibits thickness scaling
\[
k \sim L^{0.62},
\]
whereas MLG follows
\[
\frac{1}{k} = \frac{1}{K_0}\left(1+\frac{2}{L}\right)
\]
and converges to a thickness-independent limit [1808.02978]. The GH power law is attributed to the divergent in-plane thermal conductivity of 2D graphene, projected onto the helical axis [1808.02978].

Interlayer overlap and alignment are quantitatively important. The paper reports larger interlayer centroid misfit in MLG, around $\sim 1.5$ Å, and smaller misfit in GH, around $\sim 0.8$ Å [1808.02978]. The smaller misfit indicates better alignment between adjacent turns, which increases effective contact area for van der Waals-mediated transfer and reduces sliding or misalignment losses [1808.02978].

The strain response is nonstandard. In the small-strain regime, below roughly $10\%$, compressive strain enhances thermal conductivity in GH because adjacent turns are pushed closer and van der Waals coupling becomes stronger [1808.02978]. Under ultra-large tensile strain, on the order of
\[
100\% \text{ to } 500\%,
\]
and even approaching the broader elastic regime near $1000\%$, the heat current does not collapse as in generic solid-state materials [1808.02978]. After initial delamination, the structure contains a delaminated region behaving like a curved nanoribbon and an undelaminated region where van der Waals interactions remain. The reported heat current drops substantially only in the early stretching stage and then converges to a nearly constant value around
\[
0.5 \times 10^{-8}\ \text{W}
\]
in the large-strain regime [1808.02978].

The phonon interpretation is expressed through a phonon population variation ratio,
\[
A_n(\omega)= \frac{1+\int \mathrm{VDOS}_{\mathrm{GH},\epsilon}(\omega')\, d\omega'}{1+\int \mathrm{VDOS}_{\mathrm{MLG},\epsilon}(\omega')\, d\omega'} -1,
\]
with the conclusion that GH consistently has a larger out-of-plane phonon population than MLG [1808.02978]. This supports the view that GH thermal transport is a mixed in-plane/out-of-plane phonon transport problem rather than simple interlayer conduction.

## 4. Curved-surface Dirac physics and helicoidal nanoribbons

The electronic structure of graphene helicoids has been analyzed by treating graphene carriers as Dirac fermions constrained to a curved surface with screw symmetry [1505.05159]. In the general isothermal-coordinate construction, the curved-space Dirac equation is written using zweibeins and covariant derivatives,
\[
i\sigma^a e^\mu_a D_\mu \psi = E\psi,
\]
with
\[
g_{\mu\nu}=e_\mu^a e_\nu^b \eta_{ab}, \qquad D_\mu=\partial_\mu+\Gamma_\mu
\]
[1505.05159]. For the helicoid, the metric factor is $g(u)=\cosh u$, and the Dirac equation becomes
\[
\frac{i}{a\cosh u} \begin{pmatrix} 0 & \partial_u-i\partial_v+\frac{\tanh u}{2}\\ \partial_u+i\partial_v+\frac{\tanh u}{2} & 0 \end{pmatrix} \begin{pmatrix} \psi_+\\ \psi_- \end{pmatrix} = E \begin{pmatrix} \psi_+\\ \psi_- \end{pmatrix}
\]
[1505.05159].

Because the coordinate $v$ is cyclic, one separates variables as
\[
\psi_{\pm}(u,v)=e^{i\ell v}\phi_\pm(u),
\]
where $\ell$ is the conserved momentum along the helical direction [1505.05159]. After the transformation
\[
s=\sinh u,\qquad \psi_\pm(u)=\phi_\pm(u)\sqrt{\cosh u},
\]
the problem reduces to a one-dimensional Schrödinger-type equation,
\[
-\frac{d^2}{ds^2}\psi(s) + \left[ \frac{\ell^2}{1+s^2} -\frac{\ell s}{(1+s^2)^{3/2}} \right]\psi(s) = E^2\psi(s),
\]
with effective potential
\[
V_\ell(s)= \frac{\ell^2}{1+s^2} -\frac{\ell s}{(1+s^2)^{3/2}}
\]
[1505.05159].

A central conclusion is the absence of bound states. Since
\[
V_\ell(s)\to 0 \qquad \text{as } s\to \pm \infty,
\]
the asymptotic region is free and the spectrum consists of scattering states rather than localized bound states [1505.05159]. The paper concludes that massless Dirac electrons do not form bound states on the helicoid, and further shows that bound states remain absent even after adding a mass term [1505.05159].

Scattering therefore becomes the dominant phenomenon. Reflection and transmission amplitudes are written as
\[
\mathcal{R}=Re^{i\delta_R},\qquad \mathcal{T}=Te^{i\delta_T},
\]
and the potential obeys the symmetry
\[
V_\ell(s)=V_{-\ell}(-s)
\]
[1505.05159]. The local density of states is defined by
\[
\rho(\mathbf{r},E)=\frac{1}{\pi}\operatorname{Im}\operatorname{Tr}G(\mathbf{r},\mathbf{r},E),
\]
and decomposed as
\[
\rho(s,E)=\sum_\ell \rho_\ell(s,E).
\]
The reported numerical behavior is that the LDoS has large oscillations near the helicoid axis and becomes more uniform far from the axis [1505.05159]. The helicoid thus acts as a smooth curved scattering region rather than a trap.

## 5. Geometry-induced pseudo-electric fields and chiraltronics

A distinct continuum treatment of helicoidal graphene nanoribbons emphasizes an iso-spin-dependent separation mechanism generated by the twist itself [1610.03742]. In that formulation the stationary curved-surface Dirac equation is
\[
\left( \frac{\hbar v_F}{\sqrt{1+\omega^2 \xi^2} }\, \sigma_1\partial_x + \hbar v_F\, \sigma_2 \partial_\xi \right)\psi = E\sigma_3 \psi,
\]
with separated spinor components
\[
\psi_{A,B}(x,\xi)= e^{ik_x x}\chi_{A,B}(\xi)
\]
and quantized longitudinal momentum
\[
k_x = m\omega,\qquad m\in\mathbb Z
\]
[1610.03742].

The two iso-spin components satisfy effective Schrödinger-like equations,
\[
-\partial_{\xi}^2 \chi_A + U_A(\xi)\chi_A = -k_\xi^2 \chi_A,
\]
\[
-\partial_{\xi}^2 \chi_B + U_B(\xi)\chi_B = -k_\xi^2 \chi_B,
\]
with
\[
U_A=\frac{k_x^2}{1+\omega^2\xi^2} +\frac{k_x\omega^2\xi}{(1+\omega^2\xi^2)^{3/2}},
\qquad
U_B=\frac{k_x^2}{1+\omega^2\xi^2} -\frac{k_x\omega^2\xi}{(1+\omega^2\xi^2)^{3/2}}
\]
[1610.03742]. The second term is the twist-induced contribution whose sign depends on iso-spin and chirality.

In the thin-strip regime,
\[
W < \frac{1}{\omega\sqrt{8}},
\]
the potentials reduce to linear form,
\[
U_A \approx k_x^2 + k_x\omega^2\xi, \qquad U_B \approx k_x^2 - k_x\omega^2\xi,
\]
which the authors interpret as opposite effective potentials from a transverse electric field,
\[
\mathcal E = \frac{k_x\omega^2}{e}
\]
[1610.03742]. This is described as reminiscent of the Hall effect, except that the field is geometry-induced rather than externally applied [1610.03742].

The sign difference is the basis for iso-spin separation onto opposing rims of the ribbon. For one chirality convention and $m\ge 0$, iso-spin $A$ and $B$ are driven toward opposite edges, with the assignment reversed by changing the sign of $m$ or the helicoid chirality $\omega$ [1610.03742]. The effective potential also hinders rim-to-rim transport; in the Born approximation the backward-scattering probability is
\[
w(\theta)\propto \sin^2\!\left(\frac{\theta}{2}\right)
\]
[1610.03742].

The same framework predicts iso-spin transitions. The transition frequency is estimated as
\[
\nu \approx v_F \sqrt{|m||n|^3\,\frac{2\pi W}{L^3},
\]
and, using the thin-strip condition,
\[
\nu \approx |n| \frac{v_F}{L}\sqrt{\frac{|m|}{2\sqrt{2}.
\]
For a micron-sized ribbon with $L\sim 10^{-6}\,\text{m}$ and $v_F\sim 10^6\,\text{m/s}$, this gives
\[
\nu \sim 10^{12}\,\text{Hz},
\]
that is, the THz range [1610.03742]. The authors describe this program as “chiraltronics,” in which geometry replaces an external field and the relevant degree of freedom is graphene iso-spin or chirality [1610.03742].

## 6. Helical trilayer graphene, local reconstruction, and topological flat bands

Helical trilayer graphene extends the helicoid concept from a single twisted ribbon to a multilayer moiré structure with distributed twist [2305.03031]. The proposed configuration,
\[
(\theta_1,\theta_2,\theta_3)=(\theta,0,-\theta),
\]
contains two coupled moiré patterns and nominally forms a supermoiré structure [2305.03031]. The central structural result is that lattice relaxation reconstructs this supermoiré into large periodic single-moiré domains.

The local aperiodicity measure is defined by
\[
A(\mathbf r)\equiv \sum_{l=1,3}\left|K_{lx}(\mathbf r)/K_{2x}(\mathbf r)-1\right|,
\]
with
\[
A(\mathbf r)\approx 0
\]
in the bulk of a relaxed domain [2305.03031]. The relaxed state forms a triangular tiling of large domains separated by domain walls. The locally periodic domains are denoted h-HTG and $\overline{\mathrm{h}\text{-HTG}}$, related by $C_{2z}\mathcal{T}$ [2305.03031].

A defining feature of h-HTG is a finite lateral offset
\[
d=\pm \delta = \pm \frac{1}{3}(a_2-a_1),
\]
between the two moiré superlattices [2305.03031]. This shift breaks $C_{2z}$ symmetry locally, and that symmetry breaking is identified as crucial for nontrivial valley topology [2305.03031]. At $\theta\approx 1.5^\circ$, the relaxed h-HTG domain hosts a pair of nearly flat isolated bands with valley-contrasting Chern numbers
\[
C=\pm(1,-2)
\]
[2305.03031]. In the chiral limit $\kappa=0$, the model becomes analytically tractable and yields exactly flat bands with ideal quantum geometry; at the first magic angle the exact zero-mode bands have
\[
C_A=1,\qquad C_B=-2
\]
[2305.03031].

The bands are isolated from remote bands by a large gap. The abstract states a scale
\[
E_{\mathrm{gap}\sim 100~\mathrm{meV},
\]
and the detailed estimate in the continuum model is around $85$ meV [2305.03031]. This is large compared with the flat-band bandwidth
\[
W\approx 15\ \text{meV}
\]
[2305.03031]. The geometric diagnostic for ideality is the trace-condition violation
\[
\overline{T}=\int d^2k\,\big(\mathrm{tr}\,g_{\mathrm{FS}(\mathbf k)-|F(\mathbf k)|\big)\ge 0,
\]
with ideal quantum geometry corresponding to
\[
\overline{T}=0
\]
[2305.03031]. The paper reports that HTG is close to this ideal limit, with remarkably uniform Berry curvature and charge density [2305.03031].

The correlated and topological implications follow directly from this band structure. The paper highlights the prospect of integer and fractional quantum anomalous Hall states in $C=1$ and $C=2$ bands, generalized quantum Hall ferromagnets at integer filling, and fractional Chern insulators at fractional filling [2305.03031]. It also notes that the large domains, on the scale of hundreds of nanometers, are amenable to local probes such as scanning SETs or scanning SQUIDs, and that a very small uniform heterostrain of order $0.03\%$ could relax the entire device into a single h-HTG domain [2305.03031].

When multiple domains remain, adjacent h-HTG and $\overline{\mathrm{h}\text{-HTG}}$ regions carry opposite valley-Chern responses, and the domain walls host gapless counterpropagating edge modes, producing a Chalker–Coddington-like network on the supermoiré scale [2305.03031]. This suggests a real-space network realization of helical edge transport tied directly to the reconstructed helicoidal moiré texture.

## 7. Conceptual scope, common distinctions, and research significance

The literature does not use “graphene helicoid” in a single narrow sense. One well-established usage refers to a graphene nanoribbon-based helicoid produced by a screw dislocation and studied for its mechanical and thermal properties [1709.08329], [1808.02978]. Another usage concerns helicoidal or screw-symmetric curved graphene surfaces analyzed within Dirac theory, where geometry-induced potentials modify scattering and iso-spin transport [1505.05159], [1610.03742]. A broader helical-graphene extension appears in HTG, where the twist is distributed across three consecutive layers and relaxation produces locally periodic topological domains [2305.03031].

A common misconception is that the helicoid effects are reducible to ordinary strain engineering. The continuum electronic papers explicitly stress that the helicoid is a minimal surface and that the dominant effects arise from geometry itself, not from the usual strain-induced pseudo-magnetic-field mechanism [1610.03742], [1505.05159]. Another possible misconception is that helicoidal curvature should confine Dirac carriers. The curved-surface Dirac analysis reaches the opposite conclusion: neither massless nor massive Dirac electrons form bound states on the helicoid, and the principal signatures are scattering, phase shifts, and local density variations near the axis [1505.05159].

The combined body of work shows that helicoidal graphene architectures couple geometry, van der Waals interactions, and low-dimensional transport in ways unavailable to flat graphene or ordinary multilayers. In GH nanosprings, that coupling yields very large tensile deformation capability, stage-resolved delamination mechanics, and a nearly constant elastic-limit force set by inner-edge bond loading [1709.08329]. In GH thermal transport, the same continuity and overlap produce axial conductivity that scales as a power law with thickness and remains functional under extreme strain [1808.02978]. In helicoidal nanoribbon electronics, twist generates iso-spin-dependent effective potentials and THz-scale transition estimates [1610.03742]. In HTG, a helical multilayer geometry reconstructs into $C_{2z}$-broken periodic moiré domains with flat topological bands, valley Chern numbers $C=\pm(1,-2)$, near-ideal quantum geometry, and a large remote-band gap [2305.03031].

Taken together, these results establish graphene helicoids as a structural and theoretical class in which screw symmetry or helical stacking is the organizing principle. Depending on the realization, that principle yields nanospring mechanics, anomalous axial heat transport, geometry-controlled Dirac scattering, iso-spin separation, or correlated topological moiré bands.

Source: https://www.emergentmind.com/topics/graphene-helicoids