Generation and Propagation of Ultrasonic Waves in Piezoelectric Graphene Nanoribbon

Generation and propagation of ultrasonic waves in single layer Graphene Nanoribbon is studied using semi-classical approach. When piezoelectric Graphene Nanoribbon (GNR) is exposed to time varying light beam, ultrasonic waves are produced which propagate in the medium. At low frequencies, we observed oscillations of the ultrasonic observables, velocity change and attenuation which are characteristics of massless Dirac fermions in graphene. Exploiting this oscillatory behavior, we estimate graphene’s electronic mobility to be around 5 2 10 cm V s  . Propagating ultrasonic waves can be amplified, depending on the electric field amplitude. Specifically, amplification occurs when drift velocity exceeds sound velocity. This scheme can be employed for efficient ultrasonic amplifier device operation.


Introduction
Ultrasonic waves are elastic waves consisting of frequencies greater than 20 kHz and exist in excess of Tera Hertz in graphene.Graphene utilizes ultrasonic waves in wide range of applications including acoustics, superconductivity and sensing.Ultrasonic velocity and attenuation are the important parameters required for ultrasonic technique of material characterization.The velocity is related to the elastic constants and density of material.Hence, information about mechanical, anisotropic and elastic properties of a medium can be determined from knowledge of the velocity change.
The discovery of graphene in 2004 sparked a wide range of research endeavors in this thinnest material.Despite its rich physics, there has been less exploration of acoustic characteristics of graphene.In sharp contrast, acoustoelectric properties of conventional two dimensional electron gas systems (2DEGs) has been greatly studied.Charge carriers in the 2DEG systems couples effectively with sound waves.Propagation and amplification of acoustic weaves in piezoelectric semiconductors was earlier studied theoretically in 1962 [1] and experimentally under crossed electric and magnetic fields [2].We expect that, both the gapless and linear energy dispersion nature might allow graphene carriers to couple uniquely with sound waves.
At present, using semi-classical method to explore ultrasonic properties of graphene has not been studied greatly.However, a classical approach using nonlocal elasticity theory and incorporating small scale effects was employed to analyze ultrasonic wave propagation in graphene sheet [3].We assume our graphene system is piezoelectric.This is quite good assumption, because acoustoelectric effect occurs in piezoelectric semiconductors with more than one carrier type, or in a multivalley band having well defined maxima and minima or a coaxed graphene consisting of holes with right geometry [4].
In this paper, we employ inverse piezoelectric effect to generate ultrasonic waves which interacts with massless Dirac fermions in graphene.We will consider, instead of just d.c field, a superimposed alternating field with frequency comparable to (or stronger than) the sound frequency.Considering a zero magnetic field, we will investigate frequency dependence of ultrasonic properties.Specifically, amplification characteristics of ultrasonic waves through attenuation and velocity change in graphene will be studied employing classical equations.It is often a good approximation to treat the problem of sound propagation in graphene as purely a mechanical problem once the influence of the relativistic 2D electrons is ignored.

Model and Propagation Equations
The model we considered here, which can easily be realized experimentally, is an engineered homogeneous piezoelectric Graphene Nanoribbon (GNR) or a graphene nanoribbon resting on piezoelectric substrate.A uniform applied time varying light beam along the ribbon axis generates high frequency traveling acoustic waves within the medium.This effect is inverse-acoustoelectric effect.The direct effect will only be viewed as perturbations in the medium, as we will see shortly.In Equation (1), is an a.c frequency and 0 is static electric field amplitude.GNR is purely one-dimensional along propagation direction and has completely frozen out dynamics on the in-plane transverse and out-plain flexural modes.Actually, this makes our system slightly different from similar treatment using nonlocal elasticity theory [3,5,6].

Ω E
A system of coupled equations governing propagation of the generated acoustic waves, consists of Gauss, Poisson, and diffusion equations respectively; where   electronic density,   piezoelectric coup- ling constant, dielectric permittivity,   D  electric diffusion coefficient defined as n e   , C  elastic modulation constant and n   electronic density fluctuation.An auxiliary equation for the velocity is also defined.i.e.,   .
The differential Equations ( 2)-( 5) are often solve to yield velocity change and attenuation of sound waves in a medium, if one considers two cases for the piezo-electric coupling; zero and finite coupling.In the following, we will show how these limits result in the solution.

Limit of Zero Coupling
For very weak or zero coupling of acoustic waves to the carriers, one can set 0   .In addition, the absence of any source (such as gate voltage, g v often used in de- vice applications) to control the electronic fluctuation means that 0 n   .It turns out that, as in Equations ( 2) and ( 3), only plain wave solutions for are admitted and the electric field becomes constant, 0 .However, a trial solution u E     0 plugged in Equation (5) will produce carrier position with time

Finite Coupling
Now, for 0   sound waves couple effectively to dy- namics of the carriers.The traveling ultrasonic waves caused by the external field produced perturbations on the system parameters.Specifically, a scalar potential induced is associated with an electric field     1 0 , displacement vector,   1 and number fluctuation, In view of these, we search for solutions with complex harmonic forms where   0 is the equilibrium number of electrons which produces electrical neutrality.  1 is the deviation from equilibrium which has been defined as is the sinusoidal field amplitude due to the ultrasonic wave motion in the piezoelectric medium.  1 is the induced displacement amplitude.
is the wave vector along E u q x -direction.Substituting the harmonic solutions into Equations ( 2)-( 5) immediately yields We have defined Here, we do not make any distinction between longitudinal velocity,  


Obtaining exact solution to these complex coupled Equations ( 10) and (11) is often difficult.Usually, the problem is well handled by adopting simplifying small signal approximations.For instance, the sound wavelength can be taken to be far longer than the electronic mean free path.Also, we assume a simple time dependence, together with time averaging and retaining only order.In the argument of the Bessel function, zeta can be written as where we have used the relations In the following, we defined system parameters such that  is the drift frequency and  is the normalized electric field (frequency).In Equation ( 12), we demand that so that a non-trivial frequency dispersion solution for the ultrasonic waves in terms of frequencies

Iterative Solutions
In this section, we present approximate solution to Equation (13).The equation is complicated, since both sides contain  .We seek an iterative solution by considering a situation where 0 This is a zero order solution.It is very trivial and does not present any relevant physics.For second order iteration, the zero-th order approximation is substituted into the left hand side of the general solution to find Following this iterative procedure, and taking into account the fact that 1   , the  approximate solution takes the form It is immediately obvious from the preceding equation, that the presence of the complex term is an indication of wave attenuation in graphene.

Discussions and Conclusions
The velocity change and absorption of ultrasonic waves by Dirac fermions are defined as and respectively.There are quite a number of system constants appearing in the ultrasonic physical observables of Equations ( 18) and ( 19).These must be fixed in order to proceed with any meaningful discussions.In the next section, we will attempt to estimate some of these parameters.Specifically, the electron mobility which fixes .
In this section, we estimate electron mobility and the minimum electric field required to achieve the maximum reported electron mobility of 2 230, 000 cm V s  .In order to do this, it is appropriate that we fix the parameter that is valid for the conditions The Fermi energy F E  , so that thermal processes are minimal.In this limit, interactions of Dirac electrons with ultrasonic waves become an intra-band affair.Finally, we used m s as sound speed in graphene.For time scales longer than the relaxation time,  (and for Ω 1   ), the argument of the Bessel function in Equation (13) modifies to We plot attenuation (absorption) with q and identi- fied the roots as points of minimal or zero absorption of sound waves.A zero of Equation ( 19) is obtained at 2.389 q  . This yields mobility of about 2 125, 000 cm V s  for   2 0 . The numerical value has the same order as the one measured in [7].The value can go higher for strong a.c frequencies, or for 10 V m E  Ω .
  Since higher frequencies are possible in graphene as demonstrated in [8,9], where typical acoustic wave frequencies range from .The sensitivity of mobility to the electric field parameters indicates that 0.1 -10 GHz  can actually be tuned by those external factors, and

Velocity and Attenuation Change with Frequency
Here, we observe the behavior of physical observables as frequency and electric field amplitude are varied.Using Equations ( 18) and ( 19), we first plot velocity and attenuation with ultrasonic frequency in Figure 1.It is clear that at high frequencies velocity and attenuation disappears.Also, at small values of electric field amplitude, where 1   , both observables are critically dam- ped.However, as the drift frequency exceeds the sound frequency, i.e. 1   , the curves quickly falls before beginning to oscillate at very low velocities where maximum absorption at low frequencies occur.The reason for these behaviors may be explain as follows; the high fields create electronic bunching which initially absorbed more of the sound and slows down the electron motion.However, over time there is debauching and the ultrasonic waves can now propagate in the system.The fact that graphene has unique energy spectrum may be responsible for the oscillations of the physical observables.This unusual feature is missing in conventional semiconductors [10] and two dimensional electron gas (2DEG) systems.

Amplification Scheme
Apart from the presence of complex term in the dispersion (19) which indicates a damped wave, there is also a possibility to get amplified waves depending on the value of  .To get a feeling of the amplification, we plot ab- sorption versus  in .The equality holds at the beginning of the process.

Conclusion
We have studied ultrasonic properties of graphene sub-  ject to time varying external electric field.We observed unique oscillations of ultrasonic dispersion and absorption in contrast to conventional semiconductors that behave monotonically.This is due to the high fields creating electronic bunching which initially absorbed more of the sound energy produced and slows down the motion of Dirac electrons.However, debauching sets in after some time and the ultrasonic waves now propagate freely in the system.The unique electronic spectrum of gra-phene may account for such behaviors of the physical observables.We utilized the oscillatory nature of the absorption and velocity variation to compute electronic mobility which agrees quite well with reported values.The mobility can be tuned by the applied electric field amplitude.It is possible to have amplification of ultrasonic waves when drift velocity is larger than the sound velocity.An ultrasonic amplifier based on graphene devices and operating at higher frequencies can be very attractive.

Figure 2 .
The question now is how does one see amplification?In practical applications, amplification corresponds to negative attenuation.However, positive absorption will yield damping.From Fig- ure 2, the region 1  is where  is negative.The ultrasonic waves are amplified in this region.Thus, the condition for amplification is

Figure 1 .
Figure 1.Behavior of ultrasonic physical observables with frequency.

Figure 2 .
Figure 2. Amplification scheme for ultrasonic waves in grapheme.