Entropy Production from Spin–Vibrational Coupling in Endohedral-Fullerene Qubits Encapsulated in Suspended Carbon Nanotubes

Paramagnetic endohedral fullerenes consist of atoms or small clusters trapped inside a fullerene cage [26, 24]. In the families N@C₆₀ and P@C₆₀, the enclosed atom retains an unpaired electronic spin that is well shielded by the carbon cage from environmental perturbations, giving rise to narrow spin resonances and long coherence times [22, 38, 9, 44, 21]. The use of fullerene cages as molecular spin carriers has motivated extensive proposals for quantum information processing, including arrays of spin-bearing fullerenes assembled within carbon nanotubes [13].

Suspended CNT resonators provide exceptionally low effective mass, large vibrational motion, and high mechanical frequencies spanning MHz–GHz, thus enabling ground-state operation and strong coupling to spins [34, 12, 14, 40, 39, 10, 42, 36]. When inserted into CNTs, endohedral fullerenes form ordered or quasi-ordered 1D chains in which intermolecular spacing, host–guest interactions, and nanotube confinement can all influence spin dynamics (Fig. 1). Experimental and theoretical studies of fullerene peapods have highlighted several relevant features: (i) the CNT can stabilize the fullerene chain; (ii) it can modify electronic and magnetic interactions through confinement and charge redistribution; and (iii) it can serve as a transport channel or an electromechanical element. For quantum-control applications, the fullerene spin can therefore be modeled as interacting with a structured environment that is intermediate between a local environment and a continuum reservoir.

In this work we focus on a suspended segment of CNT hosting several PEF molecules, and concentrate fullerenes whose low-energy spin manifold defines an effective qubit. The CNT’s lowest flexural mode acts as the principal bosonic mode, while other vibrational modes and external reservoirs are treated as dissipative environments. This reduced description is appropriate when one flexural mode is tuned near resonance with a controllable spin splitting and remains spectrally well separated from the rest of the mechanical spectrum.

Both N@C₆₀ and P@C₆₀ exhibit high-spin electronic ground states originating from the encapsulated atom, together with hyperfine coupling to a nuclear spin [13]. For entropy-production analysis it is sufficient to project this richer manifold onto two states, $\{|g\rangle,|e\rangle\}$, that are selected by external magnetic fields and microwave-frequency addressing. The resulting splitting

can be tuned by a static magnetic field, although hyperfine structure and anisotropies may generate additional nearby transitions. Our treatment therefore assumes that either those additional levels are far detuned from the mechanical mode and external drive, or they contribute only perturbatively through renormalized rates and couplings.

This effective qubit approximation is standard in hybrid spin–oscillator descriptions [23, 45] and is particularly useful for identifying entropy pathways. In particular, it permits one to distinguish clearly between entropy associated with spin populations and coherences, on the one hand, and entropy associated with the vibrational phase-space distribution, on the other [5].

A suspended CNT behaves as a nanomechanical resonator with flexural frequencies ranging from the MHz to GHz regimes, depending on its length, tension, and electrostatic environment [34, 12, 14]. In the harmonic approximation, the relevant mechanical mode is described by the Hamiltonian [7]

where $a$ and $a^{\dagger}$ are bosonic annihilation and creation operators. Coherent driving can be implemented electrically or magnetically (for e.g. by coupling to an ac electric field on a nearby gate, Fig. 1), giving rise to a force term

with complex drive amplitude $\varepsilon(t)$.

Following earlier work on suspended CNT spin qubits, spin–motion coupling can arise when the nanotube displacement modulates the local magnetic field seen by the spin or alters a spin–orbit mediated quantization axis [23, 45]. In the fullerene setting considered in Fig. 1, an externally imposed magnetic field gradient is especially natural: if the fullerene spin is displaced together with the nanotube, then its Zeeman splitting becomes position dependent. Quantizing the displacement yields a spin–phonon interaction linear in the oscillator coordinate. Near resonance and after a rotating-wave approximation, this reduces to a Jaynes–Cummings type coupling, as we show in the next section.

The total Hamiltonian is written as

with

and with the effective spin–phonon interaction

where $g_{1}$ is the displacement-mediated coupling strength. Here the Pauli operators are defined on the $\{|g\rangle,|e\rangle\}$ manifold as $\sigma_{z}=|e\rangle\!\langle e|-|g\rangle\!\langle g|$ and $\sigma_{x}=|g\rangle\!\langle e|+|e\rangle\!\langle g|$.

Using $\sigma_{x}=\sigma_{+}+\sigma_{-}$, one may decompose $H_{\rm int}$ into resonant and counter-rotating contributions. Close to resonance, $\omega_{s}\approx\omega_{m}$, and for $g_{1}\ll\omega_{s},\omega_{m}$, the rotating-wave approximation (RWA) neglects the rapidly oscillating counter-rotating terms $a\sigma_{-}$ and $a^{\dagger}\sigma_{+}$. Moving to a frame rotating at $\omega_{m}$ and applying RWA yields the Jaynes–Cummings Hamiltonian (see Appendix A)

where $\Delta=\omega_{s}-\omega_{m}$ is the detuning frequency, and $\sigma_{+}=|e\rangle\!\langle g|$ , $\sigma_{-}=|g\rangle\!\langle e|$ are the spin ladder operators, and $\tilde{H}_{d}(t)$ denotes the drive term in the same rotating frame.. Equation (7) describes coherent exchange between the fullerene qubit and the CNT phonon mode.

In practice, the coupling $g_{1}$ may contain contributions from intrinsic spin–orbit–deflection mechanisms and from externally generated magnetic gradients. A useful estimate for the gradient-induced term is [23]

where $g_{e}$ is the electronic $g$–factor, $\mu_{B}$ the Bohr magneton, $\partial B/\partial x$ the local gradient, and $x_{\rm zpf}$ the zero-point fluctuation amplitude of the CNT flexural mode. This expression makes explicit the central design principle: increasing the magnetic gradient or softening the mechanical mode enhances the spin–vibrational hybridization.

The hybrid system interacts with external environments through both spin and mechanical channels. We model the reduced density operator $\rho$ of the coupled spin–mechanical system with the Lindblad master equation [35, 5]

For the mechanical mode we take the dissipator term

where $\gamma_{m}$ is the mechanical damping rate and $\bar{n}_{m}=[\exp(\hbar\omega_{m}/k_{B}T_{m})-1]^{-1}$ is the thermal phonon occupation of the mechanical reservoir [35, 5].

For the spin qubit we include relaxation, thermal excitation, and pure dephasing, and the dissipator is given by

The rates $\Gamma_{\downarrow}$ and $\Gamma_{\uparrow}$ denote spin relaxation and thermal excitation, respectively, and are determined by the bath spectral density at the spin transition frequency $\omega_{s}$ [7]. The rate $\Gamma_{\phi}$ accounts for pure dephasing induced, for example, by magnetic noise or slow charge fluctuations that shift the spin splitting. Equations (9)–(11) define the dynamical generator used throughout this work. They provide a minimal but flexible description of entropy generation in the driven hybrid system.

The full density operator contains both spin and oscillator correlations. For analytical insight, it is useful to identify two complementary limits. In the weak-hybridization regime $g_{1}\ll\gamma_{m},\Gamma_{2}$ (with $\Gamma_{2}=\Gamma_{\phi}+(\Gamma_{\uparrow}+\Gamma_{\downarrow})/2$), the mechanical mode behaves as a damped driven oscillator [7] that perturbs the spin through a structured noise spectrum. In the near-resonant strong cooperativity regime $g_{1}$ competes with the dissipative rates, so that coherent excitation exchange appreciably reshapes the entropy budget. The Wigner-function formulation developed below is especially convenient in the first regime, but it also remains useful in the second as a way to visualize the oscillator’s contribution to the joint state.

In the absence of coupling to the spin, the CNT flexural mode reduces to a damped driven quantum oscillator [7]. The Wigner-function formalism provides a particularly useful phase-space representation for this mode because it encodes the full information contained in the oscillator density matrix while offering an intuitive picture of the dynamics in terms of quasi-probability flow in position–momentum space [48, 6]. Unlike a classical probability distribution, the Wigner function can take negative values, thereby encoding genuinely quantum features such as interference and nonclassicality, while still allowing damping, diffusion, and coherent driving to be visualized in a form closely connected to classical trajectories [2, 20, 15, 7]. For linear systems coupled to Gaussian environments, this formalism is especially powerful because the time evolution of the mechanical state can be expressed exactly through a Gaussian propagator, with the deterministic drift of the phase-space center and the noise-induced broadening of the distribution appearing separately and transparently [11, 27]. This makes the Wigner representation particularly well suited for the present problem, where we wish to distinguish reversible transport from irreversible diffusion and damping, and later relate these features directly to entropy flow and entropy production.

In phase space, the reduced mechanical state is represented by the Wigner function

where $q$ and $p$ denote the canonical position and momentum variables of the CNT flexural mode in phase space, and $\rho_{m}$ is the oscillator density matrix. For linear damping and Gaussian reservoirs, the Wigner function evolves under an exact Gaussian propagator of the form [11]

with $\mathbf{x}=(q,p)^{T}$. The propagator can be written schematically as:

where $\bar{\mathbf{x}}(t;\mathbf{x}^{\prime})$ is the classical damped driven trajectory launched from $\mathbf{x}^{\prime}$ and $\Sigma(t)$ is the covariance matrix generated by environmental noise. This structure expresses the solution as a classical flow together with diffusive broadening [11, 27].

For a harmonic oscillator of mass $m$ and frequency $\omega_{m}$ with a damping kernel reducible to a Markovian rate $\gamma_{m}$, the phase-space center obeys

where $F_{d}(t)$ is the applied drive force. In the Markovian limit of quantum Brownian motion, the first moments obey damped-oscillator equations, while the second moments form a closed linear system that can be written in covariance-matrix form,

with drift matrix $A$ and diffusion matrix $D$ determined by the damping and bath temperature [5, 11].

This exact solution is especially useful because it cleanly separates three ingredients: (i) reversible phase-space transport generated by the oscillator Hamiltonian and external drive; (ii) irreversible contraction associated with damping; and (iii) noise-induced spreading governed by $D$. These ingredients map naturally onto the entropy balance discussed later in this paper. The explicit expressions for the covariance matrix $\Sigma(t)$ and the mean phase-space trajectory $\bar{\mathbf{x}}(t)$, which determine the Gaussian propagator in Eqs. (13) and (14), are derived in Appendix B.

When the spin–phonon coupling is restored, the oscillator no longer evolves autonomously. Nevertheless, the exact oscillator propagator remains useful in two ways. First, in perturbative treatments the spin can enter as a source term or as a conditioned displacement, so that the oscillator Wigner function becomes a weighted superposition of Gaussian packets associated with spin-conditioned trajectories. Second, even in the fully coupled regime, the reduced mechanical state

can be monitored through its Wigner function to diagnose the onset of irreversibility: phase-space broadening, loss of negativity, and displacement diffusion provide useful indicators of increasing decoherence and entropy production.

Writing the joint density operator in the spin basis,

one can define a matrix-valued Wigner function,

whose diagonal components describe oscillator phase-space distributions conditioned on the spin state, and whose off-diagonal components encode spin–mechanical coherence. The Lindblad equation induces coupled partial differential equations for these entries. In the weak-coupling limit, the diagonal components can remain approximately Gaussian, while the off-diagonal terms typically decay on a scale set by $\Gamma_{2}$ together with coupling-induced and mechanical dephasing.

The von Neumann entropy of the reduced mechanical state, $S_{m}=-\mathrm{Tr}(\rho_{m}\ln\rho_{m})$, generally cannot be written exactly as a simple functional of the Wigner function because $W(q,p)$ need not be positive. However, for a single-mode Gaussian state the entropy is completely determined by the symplectic eigenvalue of the covariance matrix $\Sigma$ [5]. Defining

with $\nu\geq 1$, one may identify an effective thermal occupation number through $\nu=2\bar{n}+1$, or equivalently $\bar{n}=(\nu-1)/2$. Because any single-mode Gaussian state is unitarily equivalent to a thermal state with this occupation number, its entropy is

Equation (21) makes the relation between the covariance matrix and the mechanical entropy explicit: for a single-mode Gaussian state, $S_{m}$ is a monotonic function of the symplectic eigenvalue $\nu$, and hence of $\det\Sigma$. Accordingly, whenever diffusion and dissipation increase $\det\Sigma$, the mechanical entropy also increases.

The total nonequilibrium entropy of the coupled spin–mechanical system is the von Neumann entropy

The reduced density operator evolves according to the Lindblad master equation

where $\mathcal{L}$ is the full dynamical generator and the $L_{j}$ are the jump operators [5]. Differentiating $S(t)$ with respect to time gives

Since $\mathrm{Tr}\rho=1$ at all times, the second term vanishes, so that

The Hamiltonian part of $\mathcal{L}$ does not contribute to $dS/dt$, because the trace of the commutator term vanishes by cyclicity of the trace and the fact that $\rho$ commutes with $\ln\rho$. Thus the entropy change is governed entirely by the dissipative part of the evolution.

To identify the irreversible contribution, we introduce a stationary reference state $\rho_{\rm ss}$ satisfying

and define the quantum relative entropy

Differentiating with respect to time yields

which motivates the definition of the entropy-production rate

Under the standard assumptions of complete positivity and stationarity, Spohn’s inequality ensures that [5]

Thus $\Pi(t)$ measures the irreversible entropy generated by the open-system dynamics. Expanding Eq. (29) and using Eq. (25), one obtains

This suggests defining the entropy-flux rate as

so that the entropy balance takes the form

For thermal reservoirs, $\Phi$ reduces to the usual heat flux form when $\rho_{\rm ss}$ is Gibbsian. In the driven CNT–PEF system, however, coherent forcing and multiple baths generally imply that the relevant stationary state is not a simple global Gibbs state, making Eq. (32) the more general definition. In a nonequilibrium steady state, $dS/dt=0$, and Eq. (33) reduces to

Because the dissipator is the sum of spin and mechanical contributions, the entropy production may be decomposed as

where $\Pi_{m}$ and $\Pi_{s}$ are associated with $\mathcal{L}_{m}$ and $\mathcal{L}_{s}$, respectively, and $\Pi_{\rm corr}$ accounts for the fact that the logarithm of the full state does not generally decompose additively once spin–mechanical correlations develop. In the weak-correlation regime one may approximate

which yields

This decomposition is useful for identifying whether irreversibility is dominated by vibrational diffusion or by spin relaxation and dephasing [17].

For a thermal mechanical bath with reference state

one finds

where where $J_{m}$ denotes the rate of energy transfer from the mechanical mode to its thermal reservoir in the weak-correlation (factorized) limit. An analogous expression applies to the spin bath when $\Gamma_{\uparrow}/{\Gamma_{\downarrow}}=e^{-\beta_{s}\hbar\omega_{s}}$ [17, 5].

For the oscillator sector, the Wigner formulation permits a semiclassical entropy analysis when the state remains Gaussian or nearly Gaussian. In the linear Markovian regime, the phase space dynamics can be written in continuity-equation form [5],

where $\mathbf{J}_{\rm rev}$ is the reversible current generated by Hamiltonian flow and $\mathbf{J}_{\rm irr}$ is the irreversible current arising from damping and diffusion. For Gaussian states, whose Wigner functions are positive, one typically introduces the Shannon-like phase space entropy

which provides a useful coarse-grained measure of the oscillator’s entropy. Although $S_{W}$ is not, in general, identical to the von Neumann entropy, for Gaussian states it is closely related to it through the covariance matrix and therefore serves as an intuitive diagnostic of diffusion-driven irreversibility [17, 35, 5].

In this representation, the irreversible phase-space current can be used to express entropy production. In particular, for linear Fokker–Planck dynamics one obtains a quadratic form involving the drift and diffusion matrices, showing that entropy is generated whenever damping-induced contraction and thermal broadening fail to balance at the instantaneous state [35, 11]. The exact propagator of Section IV therefore provides not only the state itself but also the ingredients needed to characterize the oscillator’s entropy budget.

Consider first a near-resonant regime $\omega_{s}\approx\omega_{m}$ with weak coherent drive and low bath temperature, such that $\bar{n}_{m}\ll 1$ and $\Gamma_{\uparrow}\ll\Gamma_{\downarrow}$. If the system is initialized in $|e,0\rangle$, coherent spin–phonon exchange produces Rabi-like oscillations between $|e,0\rangle$ and $|g,1\rangle$ [7]. Here $\lvert 0\rangle$ and $\lvert 1\rangle$ denote the zero- and one-phonon Fock states of the CNT flexural mode, respectively, so that $\lvert e,0\rangle$ and $\lvert g,1\rangle$ represent product states of the spin and mechanical sectors within the single-excitation manifold. In the absence of dissipation this exchange is reversible and the entropy remains low if the state is pure. Once mechanical damping and spin dephasing are included, these oscillations are damped and entropy production becomes positive.

A simple concrete example is obtained by neglecting the weak drive to leading order over a single exchange cycle and truncating the Hilbert space to the single-excitation manifold $\{|e,0\rangle,|g,1\rangle\}$. Under the resonant Jaynes–Cummings Hamiltonian, the pure state evolves as

Tracing over the spin gives the reduced mechanical state

so that the mechanical entropy is

Equation (44) shows explicitly that the reduced mechanical entropy oscillates between $0$ and $\ln 2$ , vanishing at $t=n\pi/2g_{1}$ and reaching its maximum at $g_{1}t=\pi/4+\pi n/2$. This behavior reflects reversible spin–phonon entanglement generated by coherent Jaynes–Cummings dynamics, rather than irreversible entropy production.

The corresponding reduced mechanical Wigner function also follows directly from Eq. (43). Since $\rho_{m}^{(0)}(t)$ is a mixture of the vacuum and one-phonon states (see Appendix  C):

where $W_{0}$ and $W_{1}$ are the Wigner functions of the harmonic-oscillator ground state and first excited state, respectively. Thus the phase-space distribution oscillates between a purely Gaussian profile and a nonclassical one-phonon profile, providing a direct visualization of the coherent transfer of a single excitation from the spin sector to the CNT mode Fig. 2.

To connect this coherent dynamics to entropy production, we now include weak dissipation perturbatively. In the low-temperature limit, the leading mechanical heat current into the bath is proportional to the instantaneous phonon occupation

while the corresponding entropy flux is

Likewise, the spin-bath channel contributes

and therefore

where we have neglected the thermally activated spin-excitation channel $\Gamma_{\uparrow}$.

These expressions make the redistribution of irreversibility transparent. At very short times,

showing that the mechanical entropy flux turns on quadratically as the coherent exchange transfers excitation from the spin to the CNT mode. Averaging over one Rabi period $T_{R}=\pi/g_{1}$ gives

Thus, even in this simplest regime, the formalism (Fig. 3) shows explicitly how increasing $g_{1}$ accelerates the transfer of excitation into the mechanical channel and thereby shifts part of the entropy flow from direct spin relaxation to dissipation through the vibrational mode. In the presence of weak dephasing, the oscillatory behavior in Eqs. (44)–(49) is progressively damped, and the Wigner function broadens (Fig. 4). In addition the entropy balance acquires a strictly positive irreversible component $\Pi(t)$ associated with both mechanical diffusion and spin decoherence (Fig. 5).

Under resonant or near-resonant coherent drive, the oscillator acquires a large displacement and the mechanical Wigner function is dragged along a damped orbit in phase space. In the uncoupled case the state remains Gaussian and the entropy increase is governed primarily by the covariance matrix $\Sigma(t)$. With spin–phonon coupling, the oscillator trajectory becomes spin conditioned, producing partial splitting or distortion of the Wigner distribution. These distortions reflect entangling dynamics and are accompanied by mutual-information buildup between spin and oscillator.

A simple concrete example can be obtained in the rotating frame of a resonant mechanical drive, $\omega_{d}=\omega_{m}$, in the regime $|\varepsilon|\gg g_{1}$, where $|\varepsilon|$ sets the strength of the external coherent drive acting on the oscillator (Eqn. (3)), and for times short compared with the spin-relaxation time, so that the spin may be treated as approximately frozen in the $\sigma_{x}$ basis. In that basis, the interaction $\hbar g_{1}(a+a^{\dagger})\sigma_{x}$ shifts the effective drive seen by the oscillator. If the spin is initially prepared in the state $|e\rangle$, which can be written in the $\sigma_{x}$ basis as

and the mechanical mode starts in its ground state, then to leading order the joint state evolves into a spin-conditioned superposition of displaced oscillator states,

In the above equations $|+\rangle_{x}$ and $|-\rangle_{x}$ are the eigenstates of $\sigma_{x}$, related to the $\{|e\rangle,|g\rangle\}$ basis by

In this basis the interaction is diagonal, so the strong-drive dynamics are naturally expressed in terms of spin-conditioned oscillator displacements. Here $|\alpha_{+}(t)\rangle$ and $|\alpha_{-}(t)\rangle$ are coherent states of the CNT flexural mode with spin-conditioned complex amplitudes $\alpha_{\pm}(t)$.

The coherent amplitudes obey

For a constant real resonant drive amplitude $\varepsilon$, this gives

The mean oscillator displacement is therefore large when $|\varepsilon|/\gamma_{m}\gg 1$, while the spin-dependent splitting is controlled by

Tracing over the spin gives the reduced mechanical state

so that the mechanical Wigner function becomes a sum of two Gaussian packets with common covariance matrix $\Sigma(t)$ but different centers,

where

and $\bar{\mathbf{x}}_{\pm}(t)$ are the phase-space centers associated with $\alpha_{\pm}(t)$. Equation (59) makes explicit how strong driving produces a large overall displacement, while the spin–phonon coupling generates a conditional splitting of the phase-space distribution. When the separation $|\delta\alpha(t)|$ becomes comparable to or larger than the Gaussian width set by $\Sigma(t)$, the Wigner function develops a resolved double-peak structure. For smaller separations, it appears as a broadened, distorted single packet (Fig. 6).

In the absence of diffusion, the entangling part of the dynamics can be quantified analytically. The overlap of the two coherent states is

so the nonzero eigenvalues of $\rho_{m}(t)$ are

The reduced mechanical entropy is therefore

This entropy grows from zero as the two spin-conditioned trajectories separate in phase space. In the purely coherent limit, Eq. (63) measures reversible entanglement between the spin and the oscillator. For the pure bipartite state in Eq. (53), the mutual information is

This shows explicitly that phase-space splitting is accompanied by the buildup of spin–oscillator correlations. The entropy flux to the mechanical bath can also be estimated directly. In the factorized weak-diffusion limit, the mechanical heat current is approximately

where the mean oscillator occupation is

The corresponding entropy flux is

Equations  (66) and  (67) show explicitly that in the strong-drive regime the entropy flux into the mechanical bath scales as $\varepsilon^{2}/\gamma_{m}$ and can become large even when the spin subsystem remains close to its local stationary state. This is the oscillator-dominated entropy-production regime.

Once diffusion and spin dephasing are restored, the two Gaussian packets in Eq. (59) broaden through the covariance matrix $\Sigma(t)$, while the overlap in Eq. (61) is further suppressed by decoherence (Fig. 6). As a result, the oscillatory entanglement dynamics cross over into genuine irreversibility: the Wigner function becomes increasingly broadened and distorted, the mutual information ceases to be purely reversible, and the entropy balance acquires a strictly positive production term $\Pi(t)$ associated with both finite-temperature diffusion and the continuous generation and decay of spin–oscillator correlations.

At higher temperatures, thermal excitation of the mechanical mode increases $\bar{n}_{m}$ and broadens the phase space distribution, so that the fullerene qubit interacts not with a nearly pure bosonic mode but with a thermally occupied one. Two effects then follow. First, the coherent signature of spin–vibrational hybridization is progressively washed out as thermal fluctuations randomize the spin dependence. Second, the total entropy production can initially increase with temperature, since a larger set of dissipative transitions becomes thermally accessible and the irreversible phase space current correspondingly grows.

When the drive injects work continuously into the system, the resulting nonequilibrium steady state typically violates global detailed balance even if each bath separately is thermal [17, 5]. In that case the stationary state supports a strictly positive entropy-production rate,

because $dS/dt=0$ in the steady state. This quantity provides an experimentally relevant summary of irreversibility. In the present platform, $\Pi_{\rm ss}$ can be tuned by the magnetic gradient (which changes $g_{1}$, see Eqn.  (8)), the drive amplitude, the bath temperature, and the quality factor of the suspended CNT.

A simple scaling analysis clarifies the competition between coherent hybridization and dissipation. Let the cooperativity-like ratio be defined as

For $C\ll 1$, the spin and oscillator mainly relax through their own baths, and entropy production is approximately additive. For $C\gtrsim 1$, coherent exchange competes strongly with dissipation, spin–oscillator correlations become significant, and the additive approximation fails. In this crossover regime one expects the correlation contribution $\Pi_{\rm corr}$ to become important, especially under continuous coherent driving. Introducing the dimensionless drive parameter

which provides a convenient measure of drive strength relative to mechanical damping, one distinguishes weakly displaced states ($\eta\ll 1$) from strongly nonequilibrium mechanical states ($\eta\gg 1$). The entropy budget then crosses over from relaxation dominated to drive sustained, with the mechanical entropy flux scaling as $\eta^{2}$ in the linear-response regime and remaining the dominant contribution once the resonator is held far from its thermal reference state. We provide representative values for the parameters, which indicate that this thermal crossover should be experimentally accessible in suspended CNT devices. Using the length-dependent estimates of Ref. [7] for $Q=10^{4}$, one finds $\omega_{m}/2\pi\simeq 5.37~\mathrm{GHz}$, $221~\mathrm{MHz}$, and $54~\mathrm{MHz}$ for $L=100$, $500$, and $1000~\mathrm{nm}$, with corresponding relaxation times $T_{1}\simeq 0.29$, $7.4$, and $29.6~\mu\mathrm{s}$ and, in the relaxation-limited regime, $T_{2}\simeq 2T_{1}$. At $T=10~\mathrm{mK}$ these frequencies correspond to thermal occupations $\bar{n}_{m}\simeq 6\times 10^{-12}$, $0.5$, and $3.4$, respectively, showing explicitly the crossover from a near-ground-state mechanical mode in short CNTs to a thermally occupied mode in longer devices  [7].

A complementary benchmark is provided by the suspended-CNT spin–phonon estimates in Ref. [23] who obtained $g_{1}/2\pi\simeq 0.56~\mathrm{MHz}$ for the spin qubit and $g_{1}/2\pi\simeq 0.49~\mathrm{MHz}$ for the transverse Kramers ($K_{x}$) qubit at $\omega_{m}/2\pi\simeq 500~\mathrm{MHz}$, together with a mechanical damping rate $\gamma_{m}\simeq 5\times 10^{4}~\mathrm{s}^{-1}$ for $Q\simeq 6.3\times 10^{4}$. Combining these values with a representative transverse decoherence rate $\Gamma_{2}\sim 10^{5}~\mathrm{s}^{-1}$ yields $C\sim 10^{4}$, i.e. well within the correlation-dominated regime. In the same work, the strong-drive response was analyzed at $T=50~\mathrm{mK}$ with drive scale $\omega/2\pi\simeq 0.027~\mathrm{MHz}$. Identifying $|\varepsilon|$ with this order of magnitude gives $\eta\simeq 1-3$, consistent with a driven regime in which the mechanical entropy flux can dominate the steady-state entropy budget.