Motion-Induced Radiation Due to an Atom in the Presence of a Graphene Plane

: We study the motion-induced radiation due to the non-relativistic motion of an atom, coupled to the vacuum electromagnetic ﬁeld by an electric dipole term, in the presence of a static graphene plate. After computing the probability of emission for an accelerated atom in empty space, we evaluate the corrections due to the presence of the plate. We show that the effect of the plate is to increase the probability of emission when the atom is near the plate and oscillates along a direction perpendicular to it. On the contrary, for parallel oscillations, there is a suppression. We also evaluate the quantum friction on an atom moving at constant velocity parallel to the plate. We show that there is a threshold for quantum friction: friction occurs only when the velocity of the atom is larger than the Fermi velocity of the electrons in graphene.


Introduction
Casimir and Casimir-Polder forces are physical manifestations of the vacuum fluctuations of the electromagnetic field which involve, respectively, the interaction between static bodies and between an atom and a body. The dependence of those effects on the geometry of the system, as well as on the electromagnetic properties of the material media and the atom, has been intensively investigated in the last decades, at both zero and nonzero temperatures [1][2][3].
Graphene, a single layer of carbon atoms, can be effectively described as a twodimensional material. It owes its remarkable physical properties to its planar hexagonal crystal structure and to the fact that the electronic degrees of freedom can be described, at low energies, as Dirac fermions: indeed, they satisfy a linear dispersion relation, i.e., they behave as massless fermions that propagate with the Fermi velocity v F 0.003c [4]. This yields an unusual behavior for the conductivity, as well as peculiar transport and optical properties [5,6]. These features have understandably raised the interest in the analysis of the interaction of graphene with the vacuum electromagnetic field fluctuations [7].
In natural units, the response function of graphene is determined by two dimensionless quantities, v F , and the fine structure constant. Therefore, in the simplest configuration of two planar graphene sheets separated by a distance d at zero temperature, simple dimensional analysis implies that the static Casimir force has the same distance dependence as for perfect conductors; namely, it is proportional to 1/d 4 (the force is weaker than for ideal conductors). However, new interesting effects appear at finite temperature [8]. Besides, the phenomenology becomes richer when considering more realistic descriptions for graphene, for example by including a gap in the dispersion relation or a non-vanishing chemical potential [9,10].
From a theoretical standpoint, the system is often amenable to a full ab-initio description in the context of a continuum quantum field theory, treating the microscopic degrees of freedom as Dirac fields in two dimensions [11]. For the alluded Casimir and Casimir-Polder forces, the results of this approach have been shown to be equivalent to the ones derived from Lifshitz theory [10].
When atoms or macroscopic bodies are set in motion, new phenomena may appear. Motion induced radiation, or dynamical Casimir effect (DCE), consists in the production of photons from the vacuum due to acceleration [12][13][14][15]. A related but somewhat different effect is quantum friction (QF), which manifests itself on objects with a constant relative sliding velocity. Its origin is in the excitation of microscopic degrees of freedom of the bodies, mediated by the electromagnetic field. This results in a contactless dissipative force [16][17][18][19]. The interplay between atomic motion and quantum vacuum fluctuations has been analyzed in different contexts [20][21][22][23].
The DCE has been investigated in different models and physical setups, starting from the simplest case of a one-dimensional field theory in a cavity with a moving boundary [24], then a single accelerated mirror [25], and more recent calculations which focus on resonant effects [12][13][14][15]. The microscopic origin of the DCE [26,27] can be traced back to photon emission processes in accelerated atoms: indeed, even when the oscillation frequency is not sufficient to excite the atom, the coupling to the quantum electromagnetic field allows for the emission of photon pairs. Production of photons out of the vacuum can also occur in situations where the electromagnetic properties of a media change with time [28,29]. A practical implementation of this, involving superconducting circuits, has led to the first experimental observation of the DCE [30].
QF has also been discussed in a variety of situations, most of them involving planar sheets or semi-spaces electromagnetically described in terms of their dielectric properties. Due to its short range and very small magnitude, QF has eluded detection so far. Despite this, there could be observable traces of this effect in the velocity dependence of corrections to the accumulated geometric phase of a neutral particle which moves with constant velocity in front of an imperfect mirror [31]. It has also been shown that QF could influence the coherences of a two-level atom [32]. Furthermore, an innovative experiment was designed to track traces of QF by measuring this corrections to the geometric phase [33]. This experimentally viable scheme can spark, we believe, hope for the detection of noncontact friction.
Within the previous context of quantum dissipative effects, in this paper, we study a particular system which consists of an (externally driven) atom which moves non-relativistically (with or without acceleration) in the presence of a planar, static graphene sheet. Both systems are coupled to the vacuum electromagnetic (EM) field, which mediates the correlation between the quantum fluctuations of the charged degrees of freedom, located in the atom and on the graphene sheet. The atom is coupled to the EM field through its electric dipole moment. Because of the special features of graphene, for example, the dependence of its response function on just two dimensionless parameters, we expect an interesting manifestation of microscopic (i.e., one of the bodies involved is an atom) DCE and QF.
In a previous work [34], we studied QF between two sliding graphene sheets and found that there was a velocity threshold for the occurrence of friction, given by the Fermi velocity v F . This fact could be eventually relevant, in future technological applications, in order to avoid dissipative effects. Note, however, that by the same token it makes QF almost impossible to detect in that kind of system. As we show below, a similar threshold appears for the case of an atom moving with constant velocity over a graphene sheet; we note in passing that the threshold may be less difficult to overcome for the motion of a single atom than for the collective motion of a whole plate.
We also previously studied the microscopic DCE due to an atom, with an internal structure modelled by a quantum harmonic oscillator, coupled to a scalar field, in the presence of an imperfect mirror [27]. The internal degrees of freedom of the mirror were described in terms of a set of harmonic oscillators. We found an interesting effect when the mechanical frequency equals to the sum of the internal frequency of the atom and that of the microscopic oscillators, with regions of enhancement and suppression of the vacuum persistence amplitude. This behavior is similar to that of spontaneous emission for an atom immersed into a dielectric [35,36]. We expect a rather different behavior in the presence of the graphene plate, due to the already mentioned absence of a parameter with dimensions in the response function.
Many tools can be used to study the existence and magnitude of dissipative phenomena in this context; among them is the imaginary part of the in-out effective action, related with the vacuum persistence amplitude, and obtained by integrating out the quantum degrees of freedom. This results in an imaginary part which is a functional of the trajectory of the atom, the object that we evaluate here. This paper is organized as follows. In Section 2, we define the model at a microscopic level. Then, in Section 3, we integrate out the charged degrees of freedom, to obtain the effective action for the gauge field, which also depends functionally on the trajectory of the atom. Results for the remaining functional integral (i.e., over the gauge field) are presented in Section 4, in a perturbative approach. To first order, the imaginary part of the effective action gives the probability of emission of the atom oscillating in vacuum that, up to this order, and to the quadratic order in the oscillation amplitude, is non-vanishing only when the frequency of oscillation is larger than the frequency of the harmonic oscillator that describes the atom. At higher orders in the amplitude, the threshold is reduced.
The second-order contribution contains information about the influence of the presence of the graphene sheet on the empty space vacuum persistence amplitude. In Section 5, we present results about the imaginary part of the in-out effective action for different kinds of motion. We first consider the QF on an atom moving with constant velocity, parallel to the plane, and show that there is a threshold for QF. Then, we analyze the case of an oscillatory motion, perpendicular or parallel to the plane. We find that the effect of the plane on the emission probability depends crucially on the direction of motion. In Section 6, we comment on the changes in the expression for the effective action when the exact propagator for the gauge field in the presence of the graphene plane is used, and the way to implement the perfect conductor limit. Section 7 contains the conclusions of our work.

The Microscopic Model
We begin by defining the model in terms of its action S, depending on the intervening degrees of freedom: A µ , a 4-potential corresponding to the EM field, a Dirac field ψ,ψ in 2 + 1 dimensions, for the electronic degrees of freedom on the graphene sheet, r(t), the center of mass of the atom, and x(t), the position of the electron (relative to that center of mass).
S has the structure: where S 0 (A) denotes the free EM action, while S g and S a are the graphene and atom actions, respectively, each including its coupling to the gauge field. The free EM field action S 0 (A), including a gauge fixing term, is given by Indices from the middle of the Greek alphabet (µ, ν, . . .) run from 0 to 3, with x 0 ≡ ct. In our conventions, c ≡ 1, and we use the metric signature (+ − −−). The graphene action, S g , on the other hand, is localized on the region occupied by the sheet which, in our choice of coordinates, corresponds to x 3 = 0. Fields restricted to such region will therefore depend on the reduced, 2 + 1-dimensional space-time coordinates (x 0 , x 1 , x 2 ), which we denote collectively as x . For an unstrained plate, the form of the action for a single fermionic flavor becomes: where D α = ∂ α + ieA α (x , 0). Indices of the beginning of the Greek alphabet (α, β, . . .) run over the values 0, 1, and 2, and (ρ α β ) = diag(1, v F , v F ), with v F the Fermi velocity. A single flavor corresponds to two 2-component spinors, and the representation chosen for the Dirac's γ-matrices is such that parity is preserved even when the mass m = 0.
The atom, on the other hand, is described by a simplified model in which a single electron of charge e and mass m moves in the presence of a central potential V with origin at the nucleus, where most of the mass M of the atom, and a charge −e, are located.
The trajectory of the center of mass of the atom is practically identical to the one of the nucleus, and we describe it by the vector function t → r(t). Denoting by x(t) the position of the electron with respect to the nucleus, the interaction action between the atom and the EM field is: In the dipole approximation, i.e., assuming the fields vary smoothly on the spatial region where the electron's wavefunction is spatially concentrated, we expand the interaction action as After making an integration by parts and also using that |ṙ| << |ẋ|, the last expression can be written as with d the electric dipole moment d(t) = ex(t) and m the magnetic dipole m = e/2 (x(t) × x(t)). Neglecting the magnetic dipole term, one gets S The total action for the atom, S a , consistent with the approximations above, is then: where we ignore kinetic and potential terms for r, since its dynamics is externally determined.

Effective Action
The effective action Γ[r(t)] is obtained by functional integration of all the degrees of freedom, with Feynman's i prescription implicitly assumed for all the integrals, since we are interested in the in-out effective action, namely [37], where r 0 is some 'reference' trajectory. In most applications, we use a time-independent r 0 . This is useful when considering a bounded motion, where it is natural to identify it with the time average of r(t).
We find it useful to decompose the integral into two successive steps, where S eff (A; r) is assumed to be defined modulo an irrelevant constant, which cancels out in Γ. It is convenient to set: S eff (A; r) = S 0 (A) + S I (A; r) where we separate from S eff (A; r) the free EM field contribution from the one due to the charged sources. Therefore, we now need to integrate the fermionic field and the electron trajectories. In general, neither of them can be performed exactly; thus, consistently with the assumptions mentioned above, we retain the terms up to order e 2 . Namely, both produce terms which are quadratic in A.
Let us find the explicit form of those terms. In an obvious notation: In order to find S (a) I , one needs to perform the functional integral which, as usual, is assumed to vanish in the noninteracting case, i.e. S where .
This correlation function, due to time-translation and rotational invariances, regardless of the specific form of V, may be written as follows: with the precise form of ∆ a (ν) to be determined by the potential. For example, for a harmonic potential of frequency Ω, V = m 2 Ω 2 x 2 , one has the exact result: For a general potential it is interesting to note that, since x i (t)x j (t ) is the (exact) Feynman propagator of the coordinate operator, we can always invoke the spectral decomposition theorem to write: Therefore, for a general potential V, S (a) where ∆ ω (t) is the inverse Fourier transform of ∆ ω (ν).
It is noteworthy that we obtain qualitatively similar results, to the quadratic order in the charge, if we assume that atom is described by a two-level system. Indeed, to that order, we only need to know the correlation function: where T denotes the time-ordered product and |0 is the ground state of the atom. If now we assume that we have a two-level system, we may produce a more explicit expression for the correlation function. Indeed, we may insert an identity operator I in the time-ordered product, written in terms of the two eigenstates of the atom's Hamiltonian. One of those eigenstates is of course |0 with energy E 0 , and there is just one an excited level, |1 . For a rotationally symmetric potential, the excited state must have a degeneracy consistent with the conservation of angular momentum. The simplest non-trivial one is 3. Therefore, we assume that the excited level has an energy E 1 , with the three degenerate eigenstates |1 (α) , with α = 1, 2, 3. We then insert the identity written as: it is important to note that Wigner-Eckart's theorem implies thatx i can only have non-vanishing matrix elements between the ground state and the excited ones. Thus, where We see that ξ ij behaves as a second-order tensor under rotations and, because of the rotational symmetry of the system, can only be proportional to δ ij : ξ ij = χ 2 δ ij , for some constant χ. Thus, We conclude that, for this quadratic approximation to the effective action due to the atom, a two-level system produces essentially the same effect as a harmonic potential of frequency Regarding the graphene contribution, we see that: with the tensor kernel Π αβ denoting the vacuum polarization tensor for the Dirac field on the plane. Using a tilde to denote Fourier transformation, Π αβ may be decomposed into two irreducible tensors (projectors) adapted to the symmetries of the system. Introducing first the ingredients:k α ≡ k α − k 0 n α ,g αβ ≡ g αβ − n α n β and n α = (1, 0, 0), a convenient pair of projectors is P l αβ and P t αβ , where The result, for the most relevant case of m = 0, is as follows: where we introduced α N ≡ e 2 N 16 , with N the number of 2-component Dirac fermion fields. We thus find the explicit form of the two terms into which S I may be decomposed.
It is worth noting that, assuming gauge, translation, and rotation invariances, one can generalize the previous expression to other planar media as follows: where g t and g l are functions of k 0 and |k |.

Perturbative Expansion
We now consider the evaluation of the effective action, by proceeding to integrate the vacuum field, in a perturbative expansion approach.
Γ[r(t)] is given by: with Γ may be expanded in powers of S int , producing a series of terms, namely, Γ = Γ 1 + Γ 2 + . . .. Let us consider the first two of them:

First-Order Contribution
We see that the first-order contribution, Γ 1 , contains two terms, with Γ It is clear that Γ (g) 1 corresponds to a vacuum energy contribution, which is a well-known result for graphene.
We now consider Γ (a) 1 . We see that: which involves the free electric field correlation function. This may be expressed as follows: The first line in the electric-field correlation, when introduced in the effective action, induces a divergent energy shift. If l is a length characterizing the size of the atom, the energy shift reads E div = 3e 2 4mΩl 3 .
Aside from this contribution, the second, finite one, may be conveniently studied in Fourier space, where it can be related to our previous results on the scalar model [27]. Introducing the result for the finite part is, where Π(ν, p) ≡ − 2ie 2 m p 2 dp 0 2π Therefore, where From (35), we note that there is a threshold in the frequency ν, below which there is no photon emission. That is a threshold for a Fourier component of frequency ν in f (p, ν), which is not necessarily the same as the frequency of motion, unless one expands the expressions in powers of the oscillation amplitudes. Indeed, recalling the definition of f , we see that for a small, bounded motion r(t) = r 0 + y(t), where r 0 is the average position and y(t) the departure, a series expansion yields, up to the lowest non trivial order: whereỹ(ν) is the Fourier transform of y(t). Since the time average of y(t) vanishes, y(0) = 0. Upon insertion of this expansion for f into the imaginary part, we can integrate the spatial components of the momentum, to get, at the second order in the departure, the spectral form: where Due to dimensional reasons, the spectrum has a different power law than that of the scalar counterpart of this model [27], as well as different coefficient and factor of 2 (due to the polarizations of the EM field).
In particular, for a small linear oscillatory motion with frequency ν 0 and amplitude s, where T is the total time. Thus far, this corresponds to the second order in the oscillation amplitude, always within the Γ (a) 1 contribution. When considering the same oscillatory motion without expanding in powers of the amplitude, we note that f contains not just the oscillation frequency ν 0 but also its harmonics. Indeed, using the Jacobi-Anger expansion, (where J n denotes a Bessel function of the first kind), we see that all the multiples of ν 0 are going to be present in f , what means that the threshold Ω may be surpassed with a lower oscillation amplitude (if its amplitude is increased): It is interesting to note that, if one considers oscillations involving more than one direction, for example when performing a circular motion, then the phenomenon above involves sums of the respective frequencies, since the expansion above produces a series for each factor.
For the single linear oscillation, we find: A few comments are in order. On the one hand, the term quadratic in the oscillation amplitude is recovered if one assumes that ν 0 is above the threshold, and one expands in powers of the amplitude each Bessel function. The n = 1 terms produces the term quadratic in the departure. On the other hand, if the oscillation frequency is lower than the threshold, one gets a contribution starting from n > 1. The relevance of the harmonics of the fundamental frequency in the DCE was previously pointed out for the case of semitransparent mirrors in [38].

Second Order
The second-order term is found to be given by: of which, in a self-explanatory notation, only the last one involves correlations between the atom and the graphene plate (the first one deals with the atom in free space, and the second one produces a contribution to the graphene self-energy): and keeping the latter in terms of a yet unspecified vacuum polarization tensor, after some algebra, we get: In terms of the functions g t and g l , For the particular case of graphene, those functions may be conveniently represented by means of an integral representation: which, when used in (47), leads to: where we wrote explicitly the i terms in the denominators. The next step is to integrate k 0 , which appears in a similar fashion as the momentum to integrate in a loop integral in a quantum field theory system, although in this case in 0 + 1 dimensions. This can be done in several ways, for example by the method of residues. After that integration, the imaginary part of B may be written as follows: where we introduced the shorthand notation: p ≡ k 2 + p 2 3 , k ≡ k 2 + k 2 3 , and q ≡ v 2 F k 2 + q 2 3 . Note that the third term inside the integral can be integrated, by taking advantage of the δ-function involving q 3 .
Regarding the first and second terms, we see that they differ just in the interchange k 3 ↔ p 3 , so we discuss just the first one: the sum of the two δ-functions may be replaced by the product: θ(|ν| − Ω) δ(|ν| − Ω − p). Therefore, using the previous condition into the integral for the (only) factor involving q 3 , where the last equality follows from the use of a principal value prescription for the integral. Note however that the condition p = |ν| − Ω forced by the δ-function implies that v F |k | < |ν| − Ω for v F < 1. Therefore, the first and second terms in Equation (50) do vanish. Using this result, we see that the full imaginary part may be written as follows: . (52) In what follows, we apply this result to the analysis of the dissipative effects on the motion of the atom induced by the presence of a graphene sheet.

Quantum Friction
The first example that we consider corresponds to an atom moving with constant velocity u, parallel the plane. A non-vanishing imaginary part of the effective action would indicate the presence of quantum friction produced by the graphene.
The trajectory of the atom is given bay r(t) = (0, ut, a). From Equation (31), we obtain Therefore, denoting by T the total extent of the time interval, we see that the imaginary part of the effective action per unit time becomes: We see that, due to the presence of the threshold in the imaginary part of B, there will not be friction for u < v F .
Note, however, that u > v F , there will be a non-vanishing imaginary part. The physical reason is that, for that velocity, the frequency involved (k 1 u) can excite physical Dirac fermions on the graphene sheet. Note also that, as the velocity of the atom is u < 1, the integrals over p 3 and k 3 are well-defined, since the denominator in Equation (52) never vanishes.
After integrating those variables, we get: It is easy to verify that, under the assumptions mentioned above (|k | > |k 1 u| − Ω > v F |k |), the integrand is positive.
In Figure 1, we plot the imaginary part of the effective action as a function of the velocity of the atom. We see that, after the threshold and an initial plateau, quantum friction increases with the velocity. The results are qualitatively similar to those obtained for the quantum friction between graphene sheets in [34]. In that reference, the authors also showed how to obtain the frictional force from the imaginary part of the in-out effective action (54). Besides, the argument for the existence of a threshold presented in [34], for the case of two sliding graphene sheets, can be translated here in a rather direct way: consider the momentum and energy balance during a small time interval δt, assuming that both the frictional force and the dissipated energy are driven by pair creation. The only relevant component of the total momentum P of the pair for the (momentum) balance is the one along the direction of the velocity v. Relating that component of P to the frictional force, we see that: On the other hand, the energy balance reads where E is the energy of the pair. However, since the fermions are both on-shell, we have (the equal sign corresponds to a pair with momentum along the direction of v). Dividing Equations (57) and (56), and taking into account Equation (58), we see that a necessary condition for friction to happen is v ≥ v F . The situation considered here is more interesting for an eventual experimental observation of the effect, because frictional effects only exist if the speed of the sliding motion is larger than the Fermi velocity of the charge carriers in graphene, a condition that is more easily achieved for an atom than for a whole sheet (for example, for the implementation of the experimental proposal in [33]). There is no dissipative effect until u > v F = 0.003, and QF increases as the parallel speed of the atom grows. We set Λ = 8a 3 mΩ/(πe 2 α N ) and Γ (a,g) 2 = Γ.

Small Departures Normal to the Graphene Plane
In Section 4.1, we computed the vacuum persistence amplitude for an atom oscillating in vacuum. The result to lowest order in the amplitude is given in Equation (39), and without expanding in the amplitude in Equation (42).
Here, we consider the corrections to the vacuum persistence amplitude induced by the presence of a graphene sheet. To this end, we evaluate the imaginary part of Γ (ag) 2 in a situation where the atom undergoes a bounded motion, always in a direction normal to the graphene. Although the expression for the imaginary part is not perturbative in the amplitude of the atom's motion, we find here the result to the lowest order in that amplitude.
We assume r(t) = r 0 + y(t), where the mean position of the atom is r 0 = (0, 0, a), with a > 0. The departure, on the other hand, is y(t) = (0, 0, y ⊥ (t)). We just need to insert the corresponding f for this situation in the general expression; to the order we want to work, it is sufficient to use: The resulting imaginary part may be arranged, in a similar fashion as in (37), as follows: where Unlike the case of quantum friction described in the previous section, in the calculation of m (A) ⊥ (ν), the evaluation of the k 3 and p 3 integrals should be handled with care, using a principal value prescription. We include some details of the calculation in Appendix A. The result of performing the integrals allows us to write where ϕ(x) is a dimensionless function (of a dimensionless variable), given by: This integral can be computed analytically in terms of elementary functions. To do this, it is useful to perform the change of variables y = |ρ − 1|. The explicit expression reads where and Some properties of the function ϕ(x) can be derived from its analytic expression. For instance, in the limit x → 0, it tends to the large finite value ϕ(0) 1/(3 v 6 F ), while, in the large x limit, ϕ(x) x −1 sin(2x).
In Figure 2, we plot the function ϕ as a function of x = a |ν| − Ω . As anticipated, it reaches large amplitudes at very small values of x, but then approaches zero, oscillating with a small amplitude, for large x. The behavior of m ⊥ (ν) at extremely small distances a between the atom and the graphene sheet is determined by that of x 5 ϕ(x), as can be seen in Figure 3.

Small Departures Parallel to the Graphene Plane
We now consider oscillations of the atom that are parallel to the graphene sheet. The trajectory of the particle is now given by r(t) = r 0 + y(t), where, as above, r 0 = (0, 0, a), with a > 0. The departure from the mean position is y(t) = (y (t), 0) and the corresponding f function reads: The imaginary part of the effective action reads where m (ν) = 1 2 The formal expression for m (ν) is rather similar to that of m ⊥ (ν), Equation (61), with the replacement k 3 p 3 → −1/2|k | 2 . Its evaluation proceeds along rather similar steps, and we therefore omit the details. The final result reads: where The integral that defines the function φ(x) has a singularity at ρ = 1. However, it is well defined with the principal value prescription.
We highlight the relation between this result and the function ϕ(x) in Equation (63). Note that the difference is the extra factor ρ/(1 − ρ) in the integrand of Equation (71). Therefore, we write The function χ(x) can be computed analytically and can be written in terms of cosine integral and exponential integral functions. We omit this long expression here, but quote its main properties. For small values of x, we have χ(0) −1/(2 v 4 F ). This can be explicitly checked by setting x = 0 in the integral of Equation (72). On the other hand, for large values of x, it oscillates with an amplitude O(x −1 ).
In Figures 4 and 5, we plot φ and x 5 φ as functions of x = a |ν| − Ω . There is an important difference between the results for parallel and perpendicular motions. At short distances from the plane, in the case of parallel motion, the effect of the graphene is to decrease the probability of emission, while for perpendicular motion the probability is enhanced. Note that, at short distances, |χ(x)| ϕ(x), and therefore φ(x) −ϕ(x), as clearly shown in Figures 3 and 5. On the other hand, at large distances, there is a partial cancellation between the O(x −1 ) oscillations of ϕ(x) and χ(x), and therefore φ(x) oscillates with a smaller amplitude, which is found to be of order O(x −2 ).  The fact that m (ν) and m ⊥ (ν) have different signs also occurs for a moving atom in front of a perfect conductor. In that case, the different signs come from the evaluation of the two-point functions using the image method.

The Exact Propagator for the Gauge Field in the Presence of the Sheet
In Section 4.2, we calculated the imaginary part of the effective action up to the first order in both S I . Indeed, one can for example sum over all the terms which appear when expanding the coupling to the sheet. Equivalently, one can calculate the exact gauge field propagator in the presence of the sheet and subtract the free propagator (which has already been considered).
Both approaches amount, at the level of the expression for B, to an identical expression as (47), but with different functions G t and G l replacing g t and g l , respectively, and defined as follows: (73) Therefore, it is possible to consider, within the same approach, some interesting phenomena. For example, we may take the limit α N → ∞, under which both G t and G l tend to the same limit: Recalling the functions g t and g l for graphene, we see that this limit can be treated by using the graphene expressions, but with v F = 1 and α N = 2. In particular, which is the kernel determining the dissipative effects in the case of an atom moving in front of a perfectly conducting plane. Of course, the velocity threshold means that now there will not be friction, since that would require the atom to move at superluminal speeds.
For intermediate values of α N , we have: This shows that, in the small α N limit, we recover the singularities (cuts) on k 0 = v F |k | of the perturbative calculation. Note, however, that that contribution is overcome by the cut k 0 = |k | for bigger values of α N .
The effects that would result from considering the exact gauge field propagator (a resummation of the all terms involving a coupling to the medium) will, for graphene and other systems, be explored elsewhere.

Conclusions
We studied a model with an atom which moves non-relativistically, at both constant and non-constant speeds, in the presence of a planar graphene sheet. The model used for the atom is based on a dipole coupling which, as we showed, yields essentially the same results when one uses a harmonic coupling of the electron to the nucleus or when one instead implements a two-level description. This result appears when one considers the coupling of the atom to the EM field and integrates out the electron's degrees of freedom, to the lowest non-trivial order, which is quadratic in the electric field. Note that the order is exact in the harmonic case but approximate for a two-level system.
The quantum dissipative effects were studied by first deriving a general expression for a kernel which determines them (the imaginary part of the effective action) in terms of the motion of the atom. This general expression does not rely on the approximation of small amplitude motion and therefore allows us to study both QF and the DCE on the same footing.
In the next step, we considered and evaluated the effects for particular states of motion: constant speed and bounded motion. For the former, we showed that QF exhibits the same threshold as for two graphene sheets moving at a constant relative speed [34]. It is interesting to remark that there is also a threshold in QF for non-dispersive dielectrics [39]. Indeed, when considering half-spaces described by a real constant dielectric function in relative motion, a frictional force arises between them when the velocity of moving halfspaces, in their center of mass frame, is larger than the phase speed of light in the medium (this is a quantum analog of the well-known classical Cherenkov radiation). In this respect, graphene behaves as a non-dispersive dielectric medium. For atoms moving near a metallic surface, although there is no threshold, QF is exponentially small at low velocities [40].
For bounded motion, we particularized to the case of oscillatory motions with small amplitudes along directions which are either normal or parallel to the sheet, at the second order in those amplitudes. Note that, at this order, the effects of those two motions superpose. We found that, in tune with the result for a moving atom near a perfect conductor, the effect of the graphene on the imaginary part of the effective action has different signs for normal and parallel motions.
We also pinpoint an effect mentioned above for when the atom moves in vacuum, but which is of more general validity: when the approximation of small amplitudes is not made for a simple harmonic motion, the system receives excitations of not only the fundamental frequency but also from its harmonics (of course with decreasing amplitudes); remember that the position of the atom appears in the exponent of a Fourier integral. Now, we implement the principal value prescription as I − (0) = Re 1 0 dx dk 3 2π e −iak 3 dp 3 2π The calculation is now standard and yields: The integrals for n = 1, 2 can be computing using the following trick: − dI n (0) da 2 = PV dk 3 2π dp 3 2π Writing one can easily see that only the last two terms in the integral are not vanishing for a = 0. Therefore, Using a similar argument that involves the fourth derivative of I − (0), one can show that From these results, it is rather straightforward to obtain Equation (62) for m ⊥ (ν).