Numerical Study of Duffing Nonlinearity in the Quantum Dot Embedded Nanomechanical Resonator

Geometry, electrostatics, and single-electron tunneling contribute to the nonlinearity in the quantum dot embedded nanomechanical resonator, while “Duffing term” is a kind of mathematics describing the third-order nonlinearity of the system as a whole. We study theoretically the influence of a variation of a mathematical parameter Fuffing term on the actual physical effect. -e position probability distribution, the average current, and the displacement fluctuation spectrum with the different Duffing parameter and electromechanical coupling are obtained through numerically calculating the Fokker Planck equation. -e mechanical bistability has been described by these quantities under different electromechanical coupling and Duffing parameters. We conclude that the nonlinearities of the nanotube resonator contribute to the mechanical bistability, which induces the asymmetry of the position probability distribution, compresses the current, and softens or stiffens the mechanical resonance frequency as the same as the electromechanical coupling to use it in mechanical engineering.


Introduction
e linear dynamic description of nanoelectromechanical systems (NEMS) is well understood when the electromechanical coupling is rather weak [1][2][3][4][5][6][7][8][9]; meanwhile, the existence of the nonlinearities induces too many interesting consequences such as the nonlinear dynamic response in the two-dimensional material membranes, the mechanical bistability in the carbon nanotube (CNT)-based resonator, [10][11][12][13], and unusual mechanical response [8,[14][15][16]. e CNT quantum dot embedded resonators have been widely investigated as ultrasensitive detectors and sensors [17][18][19]; the experimental research group reported higher record sensitivity in mass and force sensing [20,21]. Recently, people have found that strong coupling induces strong nonlinear signals in the mechanical resonance frequency and the mechanical noise in the single-or two-level quantum nanomechanical oscillator [22][23][24]. e strong coupling dominates the quadratic term in the restore force than the geometry, electrostatics, and single-electron tunneling. e nonlinear response of the system contains a mechanical resonator at nanoscale coupled to a transistor of single electron under the external drive which behaves like the Duffing oscillator [25]. e intrinsic nonlinearities of the carbon nanotube have been ignored by most of these research studies. For the given carbon nanotube-based resonators, the nonlinearities come not only from the electromechanical coupling but also from the parent materials themselves. It is significant to consider the intrinsic nonlinearities of the carbon nanotube. In this work, we pay attention to the term of the Duffing under the force of the restoring with the cubic structure of the considered displacement. We are also interested in the nonlinear CNT-based resonators which will be described by the position probability distribution, the current, and the displacement fluctuations with the different coupling and Duffing parameter [26].
In particular, nonlinearities examine the quantum coherence pure substance which is a pure state of matter that can be quite sensitive. When two identical particles do not overlap, they each have their perfect quantum system described by the quantum oscillator. However, the measured nonlinearities may have been in a miniature form or even stay purest after the quantum mechanical overlapping. It is related to the strength, density, and temperature of the considered system (quantum entanglement). During transition times, when the density concentration is high, bouncing may potentially play a role in quantum waveparticle duality. Hydrodynamics is frequently used to characterize the extremely dense medium formed, which implies local thermodynamic equilibrium and quantum decoherence well above thermal length scale [27].
In this curious and unique research work, we probe the nonlinearities, consequences of average current flow, and intrinsic fluctuations of the probability distributions, as well as the suppression of the displacement under the properties of the considered pioneer quantum system with nanotubes. e creation of a useable source which inflates relativistic interior of the considered field potential about the special oscillator that evaluates over time in the phase of quantum temperature regimes, and we also compute the various parameters for numerous order quantum characteristics by applying the one and two nanotubes for the mentioned phases at distinct temperature and momenta regimes. e nonlinearity of production systems influences the critical parameter significantly under the considered measurements at miscellaneous temperatures' dependence noise which are explored to use in thermal and mechanical applications of engineering. e study is organized as follows. In Section 2, we introduce our model with the description of the system which consists of the carbon nanotube quantum dot embedded resonator; then, the numerical method is introduced, the Langevin equation are given, and the Fokker Planck equation is discretized to obtain the stationary solution. In Section 3, we study three physical quantities such as the position probability distribution, the average current, and the displacement fluctuations spectrum with the presence of nonlinearity of the CNT and different electromechanical couplings. e mechanical bistability has been described by them. Finally, in Section 4, we present our conclusions.

Methods
e Hamiltonian is introduced by William Rowan Hamilton, who created a ground breaking reconceptualization of Newtonian mechanics called as Hamiltonian dynamics, which was substantially vital to the future of quantum physics. erefore, the Hamiltonian of the quantum dot embedded CNT mechanical resonator (see Figure 1).

We obtain
where the first term denotes the single-mode mechanical resonator. Specifically, in this study, we consider the influences of the force of restoring under the Duffing term which is cubic within displacement x. e restoring force is illustrated as where α is the Duffing parameter and then the first term reads: where the considered mass m, momentum p, and especially the resonating frequency which possess the inverse relation with the mass can be illustrated as ω 2 � (k/m). In the second term, the Hamiltonian is composed of three parts; two of them are with the specific energy of the quantum dot and electrons, respectively. We consider the CNTas a single-level quantum dot mentioned in [30] with the width of the single level with the dependence of density of the considered states. Such mathematical representations can be manipulated as Γ � Γ l + Γ r , and the last one is the energy of tunneling [28]. e detail expression of the second term is written as which possess the factor α � L(r), and ε αk determines the electrons energy summation in the exaggeration and the energy of chemical potential is μ α . However, the terms d as well as c represent the dot and destruction operator, respectively. When the tunnel of electron appears at the CNTs, the corresponding rate of tunneling can be mentioned as where ρ α represents the density of the considered states. We consider the symmetric system to simplify our findings, and in such techniques, we have to write Γ/2 � Γ l � Γ r . e last term that is more crucial and interesting in the nanomechanical systems gives the coupling strength of electromechanical for the aforementioned system, and one can read it as [30] One could find that the coupling strength is the function of the vibration displacement within the CNTs. In such quantum calculations, the factor n � d † d denotes the numbering of electrons at the average under the tunneling in CNTs, and F 0 represents the active force for the oscillator at that moment when an electron is tunneling just at the top of the CNTs.
In particular, we consider the electron tunnel in the coherent regime under the condition k B T ≪ Γ. In such circumstances, the factor Γ represents the rate of tunneling, and T measures the strength of the considered temperature. Especially, Z as well as k B denote the Planck and the Boltzmann important constants, respectively. In particular, the mechanical mode dynamics with the displacement can be illustrated by a famous equation of Langevin [14]: where A(x) represents the dissipation and η(t) is a force of stochastic that obeys the Because the effect of the time interval of the correlations functions can denote Γ − 1 , in such circumstances, the acquired Γ ≫ ω 0 . So, we approximate the considered calculations by using the function of Dirac delta and the force at an average manipulates as F e (x) � F 0 n d (x) and C(x) � F 2 0 S nn (x, ω ∼ 0), and similarly, the expression of A(x) � − F 2 0 (zS nn /zω)(x, ω)| ω�0 , within the limit of n d � 〈d † d〉 which examines the mean of the population of our considered dot. m also measures the effectiveness of the mass for the oscillator mode considered. e obvious expressions for n d , S nn , and A have been derived in detail in the Appendix of [15]. With the equation of Langevin, one can obtain the probability distributions Q(x, p, t) within the framework of the Fokker Planck form: where F(x) � F r + F e (x).

Results and Discussion
e main object of this research is to explore the nature of the nanomechanical system which is originated by quantum mechanical systems at an extraordinary momentum, and it can be feasible by using the system of carbon nanotubes' resonator with numerous cavities and their corresponding mechanical restoring force consistently. Here, we elucidate our model results about nonlinearities of the system, position probability distribution, average current, and displacement fluctuations. We are also interested to examine the characteristics of the position probability, mechanical noise, and average current which behaved meaningful significantly during the considered quantum systems at various energies and momenta, and we presume that the produced sources consist of the debris of frequencies that behave as the nanomechanical system to use in the thermal as well as in mechanical engineering.

Nonlinearities of the System.
In this section, we introduce nonlinearity in our study; the nonlinearity can have various sources. Here, we are interested to explore in the famous and unique sources which are the two within the main unit. One is directly from the electron tunneling, which contributes the F e in the total force and exhibits from the electromechanical coupling α SET . Another source is from the geometryα G and electrostatic α ES indirectly. In recent CNT resonator-based experiments [27], one can estimate the Duffing parameter α � − 1/(6m)d 3 F/dx 3 , when V g � 1V, α SET � 8 * 10 31 N /(kg * m 3 ), α G � 9.6 * 10 34 N/(kg * m 3 ), and α ES � − 3.2 * 10 28 N /(kg * m 3 ). One can write the Duffing term in the restoring force as the nonlinearities. Now, the total resultant force F(x) exerts at the considered oscillator, and such force can be seen in [30]. We also here considered that the vanishing temperature factor which implies k B T ≪ Γ, as well as the lowest voltages which are bias and demonstrate eV ≪ Γ, could be acquired with the summation of the force restoring factor within the mechanical system and the electronic contribution as mentioned in equation (2). e force due to the electronic contribution refers to [15] and can be illustrated as One could see that the Duffing parameter α, the variation of the gate energy ε 0 , and coupling ε P could contribute the nonlinearities for the system.

Position Probability Distribution.
In this section, we discuss the mechanical motion by the position probability distribution of the CNT with the Duffing parameter and the x C L C R + + + Γ L Γ R Figure 1: Schematics of the quantum dot embedded nanomechanical resonator. e mentioned CNT is clamped doubly between a drain electrode and source which has been suspended just above the given gate electrode; then, the CNT acts as a quantum dot, and single-electron tunneling is the dominant transport mechanism. electromechanical coupling. In order to solve the considered equation of Fokker Planck, we examine and investigate the given equation to solve sophistically and numerically within the linear problem. Defining the operator of the famous Fokker Planck as mentioned L in the acquired equation, For the stationary solution, z t Q � 0, we consider the oscillator of quantum mechanical systems in the presence of the symmetry, where the gate energy ε 0 � ε P /2 shows the special characteristics, and it also corresponds to adopting the most new point of electron hole symmetry for the aforementioned system [15]. Moreover, the parameters ω 0 /Γ � 10 − 3 and k B T/Γ � 10 − 2 also can be useful in the measured calculation. e units of the t, x, p, and α are ω − 1 0 , (mω 0 ) − 1/2 , (mω 0 ) 1/2 , and F 0 (mω 0 ) 3/2 in this paper. From Figure 2, the mechanical bistability of this system has been described when Duffing parameter, α � 0; these results work well with the conclusions of [15]. We discuss our results from the mechanical view when the Duffing parameter α � 0; the total force is where x � x − x 0 and x 0 � F 0 /2k. One can verify that the acquiring force is not symmetric corresponding to the acquiring points of x 0 � F 0 /2k and linear with x. e position probability distribution shows that the coupling increases over the critical coupling value ε c � πΓ. e position probability distribution displays a strong nonlinear signal which is that merely one curious sharp-peak turns to double peaks which show that the coupling induces strong nonlinear effect and the considered system tilts into the regimes of the bistable ε p > 1 which originate from basic single stable regime ε p ≤ ε c . For the critical value, our system distribution possess flat and wider sharp peaks. All these interesting consequences could be explained by the potential of the system. One could see that the nonlinear effect from the coupling is dominant in the tunneling process than the Duffing nonlinearity. e position probability distribution with different Duffing parameters shows that the nonlinearity of the nanotube resonator also contributes to the mechanical state's transition; it induces the asymmetric of the position probability distribution with x 0 . It could push or slow down the process of electron tunneling.

e Average Current.
Because of its mobility, a quantum particle like an electron generates an electric current. e flow of its likelihood is linked to that current. ese currents are determined by the wave function that explains a particle's state. At a higher level, the quantum operator that can operate between states or interact with a density matrix to determine currents even in scenarios such mixed levels of thermodynamic equilibrium which can be found. Here are the concepts and equations which were used to utilize these themes.
e current in our model under the coherent tunneling regime shows k B T ≪ Γ and under the certain condition implies V, T ≪ Γ [29]. In such circumstances, the limit of the lower voltage obeys the condition of eV ≪ k B T, and it can manipulate as [30] I � e 2 V 2π τ(z), (11) where τ � 1/(1 + z 2 ) and z � (μ − ϵ 0 + F 0 x)/Γ. It is easier to commute the average of the current ω ≪ Γ with the distribution of the probability P st by integrating x: From Figure 3, one could see that the strong electromechanical coupling strongly suppresses the current, which is called the current blockade. e presence of the Duffing nonlinearities could accelerate or slow down the current's declining, especially in the region of the bistable state.

Displacement Fluctuations.
In this section, we investigate the effect of the electromechanical coupling and the Duffing nonlinearity on the system from the displacement fluctuations spectrum. e displacement fluctuations' spectrum is as [13] Here, it is crucial to indicate that the whole terms which acquired a hat in (13) represent the super quantum operators those acting in the regime of the probability space and can be illustrated as x(t) � x(t) − 〈x〉.
From Figure 4, one could see that, in the mechanical frequency ω 0 region when the system contains only monostable state, there is appearing unit resonance frequency with the sharpness peak. Once the considering coupling increases over the certain value, the mentioned sharp-peak transference to the extremely low frequency and the corresponding width of mentioned peak becomes broaden and broaden with the softening of an interesting frequency which represents the mechanical phenomena. e Duffing nonlinearity can contribute to the resonance frequency shift under the specific condition that the considered system goes to that regime which is bistable. However, the appearance of the doublet peak originates as the electrons can transform from one identical state to distinct another state. erefore, the identical fluctuation of displacement spectrum features has also been obtained in [15] and [22]. Both situations appear for a system of nanomechanical within the regime of the semiclassical and the Frank--Condon principle [11]. When the Duffing term α is negative, the resonance frequency shifts to a higher frequency, meaning the stiffening spring behavior. Conversely, within the limit of α which took positive and the corresponding resonance certain peak shifts towards the critical lower frequency, meaning the softening-spring behavior.

Conclusions
We afford the research of the nonlinearities in the nanomechanical system of a carbon nanotube resonator strongly coupled to a quantum dot. e two main sources of the nonlinearities are given from the mechanical restoring force; one comes from directly the single-electron tunneling, especially with the presence of the strong electromechanical coupling; another nonlinearity is from geometry, electrostatics, and electron tunneling.
e Duffing term is considered that is cubic in displacement included in the restoring force in our paper. e position probability, the average current, and the mechanical noise with different Duffing parameters have been studied. It is found that the Duffing nonlinearities of the CNT-based resonator could induce the asymmetry of the probability distribution of the current blockade, soften, or stiffen the mechanical resonance frequency that is important as the electromechanical coupling in the nanomechanical systems.
We have explored the nonlinearities of the system, position probability distribution, average current, and displacement fluctuations. We have also investigated the characteristics of the position probability, mechanical noise, and average current with various factors which behaved meaningful during the considered quantum systems at distinct energies, momenta in the short, and wide temperature regimes. One persuasive explication of this ramification is that our considered quantum system possessed the asymmetry of the probability distribution.

Data Availability
e data used to support the findings of the study are available within the article. Disclosure is paper is new and not submitted to any other journal. All the results included are new and not reported before.

Conflicts of Interest
e authors declare that they have no conflicts of interest or personal relationships that could have appeared to influence the work reported in this paper.