Nonlinear Langevin dynamics via holography

In this work, we consider non-linear corrections to the Langevin effective theory of a heavy quark moving through a strongly coupled CFT plasma. In AdS/CFT, this system can be identified with that of a string stretched between the boundary and the horizon of an asymptotically AdS black brane solution. We compute the Feynman-Vernon influence phase for the heavy quark by evaluating the Nambu-Goto action on a doubled string configuration. This configuration is the linearised solution of the string motion in the doubled black brane geometry which has been proposed as the holographic dual of a thermal Schwinger-Keldysh contour of the CFT. Our expression for the influence phase passes non-trivial consistency conditions arising from the underlying unitarity and thermality of the bath. The local effective theory obeys the recently proposed non-linear fluctuation dissipation theorem relating the non-Gaussianity of thermal noise to the thermal jitter in the damping constant. This furnishes a non-trivial check for the validity of these relations derived in the weak coupling regime. ar X iv :1 90 6. 07 76 2v 3 [ he pth ] 2 9 Ja n 20 20

by the equation serves as a textbook example of non-equilibrium statistical mechanics. Here f 2 is the variance of the force per unit mass which controls the Gaussian ensemble for N i : The classic result of fluctuation dissipation theorem then asserts that this variance of the fluctuation f 2 is directly proportional to the coefficient of linear drag/dissipation γ. From the point of view of open quantum systems, the emergence of this linear Langevin behaviour has been well-explored beginning with the classic works of Schwinger [1], Feynman-Vernon [2,3] and Caldeira-Leggett [4,5]. It is now a paradigmatic example in the study of open quantum systems [6][7][8].
Despite the great theoretical and experimental successes of this model, it is clear that the picture of a linear dissipation with Gaussian noise can only be an approximation to real systems. Any actual system exhibits nonlinearities in dissipation and non-Gaussianities in its noise. Thus, it is of great physical interest to ask whether one could extend the linear Langevin theory to a non-linear Langevin theory. Rather than adding ad hoc terms into the equation above, one may ask, in the spirit of effective theory : what is the most general and universal non-linear extension which describes the dynamics of the Brownian particle ?
As we review below, this question can be systematically addressed using the Schwinger-Keldysh formalism. The set of all terms that could be added to the theory of a Brownian oscillator, up to one derivative order and till cubic order in amplitudes were classified in [9]. This has recently been extended to quartic dissipative oscillators by the authors of [10]. The resultant non-linear Langevin theory with these additional terms describes the universal dynamics of a Brownian oscillator weakly coupled to a thermal bath. In fact, general sum rules could be written down which relate the parameters of the non-linear Langevin theory and the thermal spectral functions/correlators of the bath degrees of freedom [9]. Further, this model provides an arena where the consequences of microscopic time-reversal invariance and the non-linear generalisations of Onsager relations could be explored. New non-linear generalisations of the fluctuation dissipation theorem emerge within this context relating the non-Gaussianities of the fluctuation ensemble to thermal jitter in the damping constant. This result puts in a systematic footing some of the previous attempts at generalising FDT (see [11][12][13][14][15]).
While the results of these works generalise Langevin theory and its characteristic features in interesting new directions, the validity of their analysis rest crucially on treating the system-bath coupling perturbatively (Born approximation) and assuming that bath correlators decay exponentially fast resulting in a local effective theory for the system (Markov approximation). A typical example discussed in these works is that of a system oscillator non-linearly coupled to a set of bath oscillators in thermal equilibrium. The frequency spectrum/coupling of bath oscillators should then be chosen appropriately and the number of bath oscillators should be very large so as to guarantee fast decay of correlations. It is interesting to ask whether the structure of non-linear Langevin theory holds true beyond such weakly coupled systems in the 'Born-Markov' regime.
In this note, we will answer this question in affirmative by studying a strongly coupled model in detail. The model we study is that of a heavy quark moving in a bath of strongly coupled, thermal Yang-Mills plasma (or more generally thermal plasma of a strongly coupled conformal field theory (CFT)). The method we use to tackle this strongly coupled system is AdS/CFT duality, which maps the quark in the plasma problem to that of a string stretching from AdS boundary and probing a high temperature black brane [16][17][18].
The boundary endpoint of the string then represents the heavy quark. The string exhibits a random motion due to the Hawking radiation of the transverse string modes. These string modes also fall into the blackbrane thus resulting in dissipation of energy. The combination of Hawking and in-falling string modes then induce a Brownian motion for the heavy quark endpoint. We will use this system, widely studied within the holography literature [16][17][18][19][20][21], to compute the non-linear corrections to the quark dynamics.
As described by Feynman-Vernon [2,5], the corrections to the quark dynamics are encoded in the 'influence phase' that arises by integrating out the bath degrees of freedom in the Schwinger-Keldysh path integral for the combined system. In our case, the computation of the influence phase reduces in large N, strong coupling limit to constructing a doubled string configuration that probes the holographic dual of the SK contour. The configuration we seek is the linearised solution of the non-linear sigma model with a double black-brane spacetime as the target space. This doubled black-brane target space is the dual of the Schwinger-Keldysh contour in the CFT.
Given that we still lack a complete formalism which can address real time questions within string theory in general and holography in particular, some clarifying remarks are in order. The doubled black-brane spacetime that we

II. NON-LINEAR LANGEVIN THEORY
Our goal in this section is to introduce the non-linear Langevin theory for the Brownian particle and its description within Schwinger-Keldysh formalism. The discussion here is based on the results of [9,10] which we will refer the reader for a more detailed description of the models under consideration. The reader will find in those works a detailed justification of how the non-linear Langevin effective theory arises generically out of weakly coupled baths along with examples drawn from simple bath models.
Consider a Brownian particle moving in a spacetime of dimension d. We denote by q i the spatial position of this particle, where i = 1, 2, . . . (d − 1) denote spatial indices. In our context, we are interested in a heavy quark coupled to a CFT d plasma at a temperature T. This quark then feels a variety of forces, which at a macroscopic effective theory level result in energy/momentum dissipation and fluctuations. We will describe this by a stochastic ordinary differential equation for the position of the quark : Here N i is thermal Non-Gaussian noise drawn from the probability distribution and N 2 ≡ N i N i . We use f 2 to denote the variance of the thermal noise, Z I characterises its frequency dependence and γ is the damping coefficient resulting in the drag force. The coefficient ζ γ is a thermal jitter in the damping coefficient and ζ N is the non-Gaussianity parameter in the noise distribution. We immediately note that the system under consideration has a variety of symmetries including translational and rotational invariance. These symmetries drastically reduce the number of non-linear corrections to the linear theory that can appear in general. This is in contrast to [9,10], where a non-linear Brownian oscillator moving in a potential was studied. As we will argue below, the above equation then represents the most general non-linear Langevin equation to quartic order in q and upto first order in time derivatives.

A. Stochastic path integral
We would like to now write down an effective action/path integral which computes correlators in this stochastic theory. This can be done by using the technique of Martin-Siggia-Rose [38]-de Dominicis-Peliti [39]-Janssen [40]. We will do this in three steps. In the first step, we write down a path integral where E i [q j a , N] are the stochastic ODEs for the position q j a of the Brownian particle and P[N] is the noise probability distribution defined above. Here δ(E i [q a , N]) indicates that we choose to average over all solutions of our stochastic ODE, with an appropriate weight for the noise field. This path integral can then be used to compute noise-averaged correlators of the non-linear Langevin system.
The formal expression above, as is usual with path integrals, should be thought of as a limit of a UV regulated expression in which appropriate counter-terms are added to cancel UV divergences. A common way to regulate the expression is to discretise the time steps, add counter-terms and take the time-step δt → 0 limit. For example, to leading order in the non-linear couplings, the following counter-terms are sufficient to guarantee finite answers : The second step in the MSR-DPJ method is to introduce an auxiliary field for each direction to deal with the functional delta function. This gives Here m p denotes the mass of the Brownian particle which is introduced for dimensional reasons and q i d is the auxiliary field enforcing the stochastic equation of motion. In the third step, we perform the integral over the noise to get an effective action for the non-linear Langevin theory of the form where the effective Lagrangian takes the form We have denoted the effective Lagrangian by the symbol L SK , since this Lagrangian can alternately be thought of as coming from solving the underlying quantum non-equilibrium dynamics of the heavy quark moving in a CFT plasma using the Schwinger-Keldysh formalism, as we describe below.

B. Heavy quark in SK formalism
A dual perspective on Brownian motion and Langevin theory is obtained by thinking about their quantum microscopic origins. The microscopic description of the heavy quark is in terms of a reduced density matrix which evolves non-unitarily due to the interaction with the CFT plasma bath. To keep track of the density matrix and its evolution, one needs to move to the Schwinger-Keldysh(SK) formalism [1, 2, [41][42][43]. One introduces in the SK formalism, two fields per degree of freedom to keep track of the evolution on the ket and the bra side of the density matrix. We define q R and q L as the right-moving/Ket and left-moving/Bra degrees of freedom in the Schwinger-Keldysh contour shown in Fig.1. The effective SK action for the heavy quark is obtained by starting with the SK path integral of the whole system and then systematically integrating out the CFT bath degrees of freedom. The process of integrating out the bath introduces corrections to the particle SK action. These corrections encode the effect of the particle-bath interaction in the macroscopic effective theory of the heavy quark. Following Feynman-Vernon [2], these corrections are referred to as the "influence phase" of the heavy quark. In general, the influence phase of the quark is determined by the connected parts (cumulants) of contour ordered correlators of the CFT operators that couple to the quark. In a strongly coupled CFT, a direct computation of the influence phase on the CFT side is a near impossible task. We circumvent this difficulty by using AdS/CFT and working on the gravity side instead.
Even if we cannot compute the influence phase on the CFT side, one could deduce, based on general considerations, various structural properties of the effective action. The unitarity and thermality of the CFT strongly constrain the influence phase. Later, we will see how the answer from gravity manages to satisfy these constraints.
On a related note, the local approximation of the influence phase also satisfies a set of constraints. While the influence phase due to the CFT bath is in general non-local, if the bath cumulants decay sufficiently fast, one can viably expand the influence phase in a derivative expansion. Even if the influence phase from the CFT is hard to compute, its local approximation can be characterised using effective theory. We can write down the most general local expression consistent with microscopic constraints and parameterise our ignorance in terms of a set of phenomenological parameters.
The relevant microscopic constraints are three in number [8,42,44,45] : 1. The microscopic unitarity in the dynamics of the (particle+bath) combined system which in particular implies that the action should be zero if q i R = q i L . This condition is necessary for the heavy quark's correlators to satisfy the largest time equation. This condition is also responsible for the Lindblad structure of master equation obeyed by the heavy quark's density matrix [44,45].
2. The reality conditions which involve q i R ↔ q i L exchange and 3. Invariance of the microscopic dynamics under constant translations and rotations of the q i 's.
If we define [8,41,42] the most general quartic Lagrangian (up to terms with single time derivatives) consistent with these constraints is given by L SK written down in the previous subsection. Thus by completely different arguments than before, we arrive at the same result for the effective action of the Brownian particle. The mapping between the stochastic variables q a and q d vs q R and q L of the Schwinger-Keldysh description given by the expression above establishes a duality between two different ways of viewing the Langevin theory : one as a stochastic theory and another as a macroscopic SK effective theory for an open quantum system. To these two perspectives, gravity adds a third perspective in the form of an open string probing a black-brane.

C. Thermality, time reversal and fluctuation-dissipation relations
In the previous subsection, we mentioned that the effects of the particle-bath interaction on the quark's dynamics are encoded in its influence phase. Such influence phases have been studied in detail for several microscopic models where a Brownian particle is weakly coupled to a thermal bath [4,9,10,46,47]. In such scenarios, one can perturbatively expand the influence phase as a power series in the particle-bath coupling. The n th degree term in this expansion receives contributions from the n-point cumulants of the bath operator that couples to the particle [2,6]. This relationship between the particle's influence phase and the bath's cumulants allows one to study the constraints imposed on the particle's effective dynamics due to microscopic symmetries of the bath.
An interesting example of such a symmetry is the time reversal invariance of the bath's microscopic dynamics. Such a symmetry in the bath relates the bath's correlators to their time-reversed counterparts. These relations between the bath's correlators, in turn, lead to relations between the the coefficient functions of different terms in the non-local influence phase. In the derivative expansion of the influence phase, these relations reduce to certain reciprocal relations between the effective couplings. Such reciprocal relations were first explored by Onsager [48,49] and were later extended by Casimir [50] for the quadratic couplings in the effective dynamics of a system with multiple degrees of freedom. The discovery of these relations laid the foundations for understanding how the dynamics of an open system, which manifestly violates reversibility due to dissipation, is nevertheless constrained by microscopic time-reversal invariance in the environment.
Recently, the Onsager-Casimir reciprocal relations were generalised to the cubic and quartic couplings of a Brownian particle in [9] and [10]. These papers demonstrated that under time reversal certain Schwinger-Keldysh correlators of the bath get mapped to its Out of Time Order Correlators (OTOCs) [71]. As argued in these works, a convenient way to study the effects of such relations between the bath's correlators on the particle's dynamics is to extend its effective theory by including the contributions of the bath's OTOCs [10,51]. Such an extension of the effective theory requires the introduction of some new OTO couplings [10,51] which encode information about the bath's OTOCs. The relations between these OTOCs and Schwinger-Keldysh correlators due to microscopic reversibility in the bath constrains the particle's effective theory by connecting some of the additional OTO couplings to couplings appearing in the Schwinger-Keldysh effective theory of the particle [9,10].
Apart from the bath's microscopic reversibility, there is another source of relations between the particle's effective couplings, viz. the thermality of the bath. The fact that the bath is in a thermal state implies the existence of certain relations between its correlators which were first discussed by Kubo [52], Martin and Schwinger [53]. These relations (now commonly known as the KMS relations) connect thermal correlators which can be obtained from each other by cyclic permutation of insertions as shown below: where O 1 , O 2 , · · · O n are n operators of the bath, and β is the inverse of the bath's temperature (in units where the Boltzmann constant k B = 1). In the high temperature limit, such KMS relations between 2-point correlators of the bath lead to the following relation between two quadratic couplings in the particle's Schwinger-Keldysh effective action (see (9)): From the point of view of the dual Langevin dynamics given in (3), this implies the well-known fluctuation-dissipation relation (FDR) [54][55][56] which connects the strength of the thermal random force acting on the particle to its damping coefficient.
In [9] and [10], certain generalisations of the above-mentioned FDR were derived for the cubic and the quartic couplings in the OTO effective theory of the particle. These generalised FDRs follow from the high temperature limit of KMS relations which connect Schwinger-Keldysh correlators of the bath to its OTOCs by analytic continuation [57]. Combining these OTO generalisations of the FDR with the generalised Onsager-Casimir relations mentioned above, the authors of [9] and [10] demonstrated the existence of generalised fluctuation-dissipation relations in the dual non-linear Langevin dynamics which connect the non-Gaussianity in the thermal noise to the thermal jitter in the particle's damping coefficient. For the stochastic dynamics introduced in (3) and (4), this generalised fluctuation-dissipation relation [72] is given by

III. HOLOGRAPHIC BROWNIAN MOTION
We will now move on to our main system of interest : an open string hanging from the AdS boundary probing a black-brane geometry. In particular, we want to compute the Schwinger-Keldysh path integral, which means that the black-brane background should be doubled in some appropriate way. We will begin by describing the geometry of this doubled spacetime constructed according to the prescriptions of [22,24,25,28]. We will later see how this target space yields sensible answers for the influence phase of the heavy quark.

A. The black brane and the holographic SK contour
Following [28], we begin by considering the asymptotically AdS d+1 black-brane, written in ingoing Eddington-Finkelstein co-ordinates Here r h is the horizon radius. We have set the AdS radius to unity. The inverse Hawking temperature of this black-brane is given by We are interested in a string probing a doubled avatar of this spacetime, whose mathematical description is in terms of an appropriate non-linear sigma model living on the string worldsheet. For this reason, we will be interested in re-writing this spacetime with other more convenient co-ordinates. It is often convenient to scale out the horizon radius by introducing a scaled radial co-ordinate ρ ≡ r r h . The above metric then becomes The horizon in this co-ordinates is at ρ = 1 and the AdS boundary is at ρ → ∞. From now on, we will regulate this boundary to be at a cutoff surface given by ρ = ρ c . For the purposes of studying the string non-linear sigma model, it is useful to shift further to a negative imaginary radial co-ordinate ξ defined in the exterior of the black-brane via where η is a negative imaginary constant defined by  The co-ordinate ξ shares many features of the tortoise co-ordinate and plays a similar role in the worldsheet theory. Note that at this juncture, we have not changed in anyway the geometry of the black-brane spacetime. For now, we are merely using an imaginary co-ordinate to label points on a real Lorentzian spacetime.
In the next step, we would now like to double this black-brane spacetime so as to construct a solution with CFT SK contour as its boundary [22,24,25,28]. Such a spacetime will then serve as the holographic SK contour which computes SK correlators. The geometry of this spacetime contains various non-trivial features both along the time direction and the radial direction.
Let us begin by commenting on the geometry along time direction. The AdS boundaries of the holographic SK contour at (ρ = ρ c ± iε) asymptote to the SK contour seen by the CFT. This boundary has a doubled time contour with a thermal initial state as illustrated in the fig.2. This structure is replicated at every fixed ρ radial slice in the bulk.
The geometry along the radial direction is illustrated in the fig.3. On the radial direction, the right and the left copies of the black-brane (denoted by M R and M L ) are stitched together at their stretched horizon (i.e, at r < r h ) by a 'horizon cap'. This horizon cap region regulates the outgoing Hawking modes which blow up at the horizon in the ingoing EF co-ordinates.
Let us now be more explicit in describing the geometry along the radial direction. This can be done by extending the definition of the co-ordinate ξ to the entire radial contour. One of the primary advantages of working with the ξ co-ordinate is the fact that it can be used to describe the radial contour clearly. We take the convention that ξ(ρ c + iǫ) = 0.
Given the differential equation for ξ, an explicit expression for ξ is given by the contour integral Here the contour integral is performed over counter-clockwise holographic SK radial contour starting from a cutoff ρ c + iǫ and going towards ρ c − iǫ shown in the fig.3. By looking at the above contour integral representation, we note that the integrand has a simple pole at the horizon with unit residue. It follows that the function ξ has a branch cut with a discontinuity Thus, ξ is a good co-ordinate on the radial contour which takes purely imaginary values on M L . The imaginary part logarithmically diverges as we approach the horizon. Its real part then jumps by unity as we traverse the horizon cap region to go from the stretched horizon of M L to the stretched horizon of M R . Thus, as we move along M R , it is ξ − 1 which is purely imaginary. In what follows, we will use ξ as a useful parameter to specify where we are on the radial contour. We will conclude this discussion with the following parenthetical remark : we can alternately use the clockwise holographic SK contour as our defining contour to write We can then use 1 − ξ as an equally good co-ordinate on the doubled radial contour. Apart from the various branches of the holographic SK contour, there are also edges which join those branches together. However, as shown by authors of [24,25], these edge contributions cancel out in the final answer and hence will be neglected from hereon.

B. The Nambu-Goto action
We will next describe the non-linear sigma model with target space as the holographic SK spacetime defined in the previous subsection. For simplicity, we shall choose the worldsheet co-ordinates to be either one of {v, r} or {v, ρ} or {v, ξ}. We will choose {v, ρ} for definiteness, though we will also sometimes find it convenient to switch to {v, ξ} instead . In the holographic SK contour, we remind the reader that ρ varies from ρ c + iǫ to ρ c − iǫ.
The rest of the (d − 1) transverse x i co-ordinates then become worldsheet fields that describe the transverse motion of the string. We will denote these worldsheet fields by X i (v, ρ). The Nambu-Goto action for the string can then be written as Here h represents the determinant of the induced metric in the {v, ρ} co-ordinates . The radial integral is over the holographic SK spacetime using the counter-clockwise radial contour that goes from ρ c + iǫ (or equivalently ξ = 0) to ρ c − iǫ (or equivalently ξ = 1) as shown in the fig.3. Using the target space metric we can compute the Nambu-Goto action. It is convenient to write down the Nambu-Goto action in a hybrid notation involving both ρ and ξ. The conversion between them is given by the differential equation Using these expressions, the Nambu-Goto action takes the form Here we have used the notation introduced before in the context of the black-brane This action has to be varied and the resultant differential equations have to be solved with the boundary conditions In order to derive a tractable set of equations of motion and solve them, we proceed by substituting into this action, an amplitude expansion of the form where we take X i 2k+1 ∼ O(q 2k+1 ). Note that given the X i → −X i symmetry of the equations of motion, under which the boundary conditions flip to q i → −q i , the expansion of X i should necessarily be odd in q i . So, no even terms occur in the above expansion. If the boundary displacements of the endpoints are not too large, we can solve the relevant equations by an amplitude expansion.
As a first step, we get the quadratic action at the second order in amplitude expansion as Varying this action, we get the linearised equation of motion The appropriate boundary conditions for X i 1 are given by We note that similarly the boundary conditions for the higher order corrections are given by We will obtain the solution of this equation in the subsequent sections. Assuming the solution satisfying appropriate boundary conditions is given, we would be interested in computing the on-shell effective action on the solution. This calculation can be simplified using the following identity established by the use of linearised equations of motion : This shows that the quadratic Lagrangian is a total derivative and the whole action can be reduced to boundary terms : We will now move to the next order in the amplitude expansion. The quartic action has two terms : the first contribution is given by substituting the linearised solution to the leading non-linear correction : This is O(q 4 ) in the amplitude expansion. Another contribution at the same order is given by the substituting the sub-leading solution into the quadratic action : Using the linearised equation of motion, this term can be rewritten as a boundary term : We use this to write S (4) as a purely boundary contribution Since the solution X i 3 satisfies the boundary conditions we conclude that S (4) = 0.
To summarise, to the order we require, it suffices to obtain the solutions of linearised equations of string motion X i 1 on the holographic SK contour satisfying and substitute it into the following quadratic and quartic terms of the Nambu-Goto action : Since our focus is on linearised solutions, from hereon we will drop the subscript on X i 1 . In what follows, we will find it convenient to write the answers in terms of an effective 't Hooft coupling For example, we can use to write the action in the following form ∂ξ . (44)

C. Ingoing quasi-normal modes
Before we try to solve the linearised string equation in the holographic SK contour, we will get our feet wet by solving the familiar problem of ingoing quasi-normal modes. This does not involve any doubling of spacetime or ideas from Schwinger-Keldysh contour and is a straightforward application of familiar AdS/CFT techniques. Thus, we will be brief, relegating the technical details to the appendix A.
We are interested in the solution to the linearised string equation satisfying ingoing boundary conditions at the horizon. Equivalently, our present interest is in solving for the retarded bulk to boundary Green's function with in-falling boundary conditions at the horizon.
We begin with the linearised string equation of motion and seek a solution to it of the form where q i (v) is a given evolution of the string endpoint at the AdS boundary. Here g ω (ρ) is the retarded bulk to boundary Green's function in frequency domain which satisfies with the boundary conditions being regularity at the horizon and lim ρ→ρ c The regularity at the horizon implies Thus, g ω is the unique solution of the second order ODE quoted above with Dirichlet boundary condition at the AdS boundary (ρ = ρ c ) and Neumann boundary condition at the horizon (ρ = 1). As we show in appendix A, we can solve this problem at arbitrary orders in derivative expansion by writing with the sequence of radial functions ψ k (ρ) defined recursively via the integrals Note that dψ 1 dξ = ρ 2 − 1 with the consequence that lim ρ→1 dψ 1 dξ = 0. Further, if we assume that lim ρ→1 dψ k dξ = 0 for all k < k 0 , then the integrand for ψ k 0 has no horizon poles and is hence analytic in the ingoing Eddington Finkelstein co-ordinates. This implies that lim ρ→1 dψ k 0 dξ = 0. Thus, by induction, all ψ k (ρ) and in turn, g ω are analytic in the ingoing Eddington Finkelstein co-ordinates. The properties lim ρ→ρ c guarantee that g ω satisfies the correct boundary conditions. In the special case of d = 2, i.e., in AdS 3 , the sequence of radial functions ψ k (ρ) take an especially simple form These can be summed to get a simple closed form expression for the retarded Green's function These AdS 3 expressions nicely illustrate the general properties argued above.
Once the causal bulk to boundary Green's function with in-falling boundary conditions has been constructed, the ingoing solution for the string with a given boundary evolution q i (v) is given by the convolution integral In frequency domain, we simply have X i (ω) = g ω q i (ω). The expression above describes the linear retarded response of the string to the motion of its endpoint. It also gives a linearised description of how the excitations on the string fall into the black brane and get dissipated away.

D. Outgoing Hawking modes
We will now shift to a discussion of the other solution to the equation which is not ingoing. These are needed to describe Hawking radiation of the string modes and their effects. As we shall see, a correct IR divergence-free description of these modes should really be done within the holographic Schwinger-Keldysh contour with doubled spacetime rather than in the single copy spacetime. Nevertheless, it is useful to set the stage by describing these modes in the language of the single copy spacetime, which we will do in this subsection.
A naive attempt at getting the outgoing modes from the ingoing modes would be to just reverse the time argument. Given that the ingoing solution of frequency ω has the form g ω e −iωv , the naive time reversed solution would be g ω e iωv . But this does not quite work, since the string equations of motion in the ingoing Eddington-Finkelstein co-ordinates is The customary solution to this impasse is to pass to a co-ordinate system (say Kruskal co-ordinates) where the underlying time-reversal invariance becomes manifest. Here, we will resist the temptation to do so for two reasons: one, it is not clear such a solution generalises to the fully dynamical case where black-brane itself has formed out of a collapsing shell of matter in the past. The second reason is as follows : we take the viewpoint that the hiding/spontaneous breaking of time-reversal invariance is in fact a physical fact of Brownian motion as it emerges from the microscopic dynamics. We prefer to keep this physics manifest rather than camouflaging it.
How then do we deal with the hidden/spontaneously broken time reversal invariance in this system ? We will characterise it by an order parameter field for time-reversal symmetry breaking, which we will denote by the letter χ. We will demand that time reversal invariance be non-linearly realised in our system via the map The normalisation factor of iβ here is chosen for notational convenience.
In the context of the solutions to linearised string equations, this practically amounts to demanding that g ω e βωχ e iωv be the time-reversed counterpart of the ingoing solution. Substituting the ansatz X ∼ g ω e βωχ e iωv into the equation above and using the fact that g ω satisfies the differential equation we get a differential equation for χ of the form We already know the ingoing mode analytic at the future horizon to be X ∼ g ω e βωχ e iωv = g −ω e iωv , which gives one solution for χ in the above equation. The other solution can be simply obtained by setting Here we have used the fact that dψ 1 dξ = ρ 2 − 1, where ψ 1 is the function that appeared in the last subsection as the one that governs the retarded response of the string to leading order in small ω expansion. Thus, we conclude that if we take χ = ξ + ψ 1 , then g ω e βωχ e iωv is the required time-reversed solution.
Given that the retarded Green's function takes the form its time-reversed counterpart is given by Thus, time reversal at the level of linearised solutions is achieved by a replacement ψ 1 → ξ − ψ 1 keeping all other ψ k 's intact. Another perspective on χ is obtained by writing its differential equation in terms of ρ derivative: We recognise here the inverse of the black-brane redshift factor which then identifies χ as the tortoise co-ordinate up to a pre-factor. A cognizant reader might have realised that this is to be expected : from the point of view of Kruskal co-ordinates, the time reversed solution that we have constructed is just g ω e iωu with u being the outgoing Eddington Finkelstein co-ordinate. Converting u to v gives a shift proportional to tortoise co-ordinate which explains the origin of χ from the Kruskal viewpoint.

E. The doubled string on holographic SK contour
We will now put together the above two solutions to determine the configuration that satisfies the correct boundary conditions. We begin with the most general solution of the linearised string equations, i.e., where we have redefined ω in the Hawking mode to have the standard sign for Fourier transform. In frequency domain, we can simply write This solution has a causal Green's function part c ω which represents the in-falling quasi-normal mode, superposed with the Hawking mode h ω . While the ingoing quasi-normal modes are regular and give a smooth description of the string dynamics, this is not true about outgoing Hawking modes which are divergent at the horizon due to the χ function. Their physical import is unclear a priori, nor is it evident how to extract finite answers from them staying within outgoing EF co-ordinates. As we will now argue, the correct interpretation of the Hawking modes is that they appear as saddle points in the Schwinger-Keldysh path integral of the string world-sheet. We can embed the string solution into our holographic SK contour by interpreting the ξ's in the above expression as varying over the holographic SK contour, i.e., the ξ function begins from zero in the left contour boundary (at ρ = ρ c + iε) and ends in becoming unity in the right contour boundary (at ρ = ρ c − iε) after traversing the contour. This is a simple lifting of the one copy solution to a configuration in the entire holographic SK contour, whereby the problematic near horizon region gets excised. The ξ and hence the χ function are now completely regular on the new contour.
One can now interpret the complex spacetime, i.e., the horizon cap region between two stretched horizons as the geometry which captures the Hawking radiation process within Schwinger-Keldysh formalism. Since Hawking radiation can often be thought of as a tunnelling process (see [58] and references therein), one can interpret the string configuration within the horizon cap region as a 'Schwinger-Keldysh instanton' that computes the corresponding tunnelling probability. This perspective is familiar in the activation type problems dealt within Schwinger-Keldysh formalism [5,8]. Such saddle points with the difference fields turned on, are crucial in capturing the physics of dissipative tunnelling processes, It is thus interesting that this tunnelling perspective of Hawking modes emerges naturally out of holographic SK formallism.
Let us return back to our solution to impose the appropriate boundary conditions. We would like to impose We need g ω (ρ c ± iǫ) = 1 along with and This yields the equations in the frequency domain. These are readily solved to give the Son-Teaney doubled string configuration [17] Here, f ω is the Bose-Einstein factor given by The Bose-Einstein factor f ω satisfies 1 + f ω = e βω f ω and 1 + f ω + f −ω = 0. Using this, the solution above can also be written equivalently as The solution satisfies a reality condition imposed by demanding invariance under the exchanges Thus, we demand that [X i (−ω)] * takes the same form as X i (ω) with an exchange of the right/left boundary value data. We have (using 1 + f ω = e βω f ω and 1 + f ω + f −ω = 0) Since, q i R , q i L are real functions in time domain, they satisfy [q(−ω)] * = q(ω). For the same reason, we have [g −ω ] * = g ω for the retarded Green's function. Complex conjugating the holographic SK contour maps ξ to 1 − ξ * and hence χ to 1 − χ * . With these, it follows that [X i (−ω)] * takes the same form as X i (ω) after a complex conjugation and right/left exchange.

IV. NON-LINEAR LANGEVIN THEORY
Before proceeding to use this configuration to compute the influence phase, it is useful to define the following combinations that recur in our solutions : These combinations correspond to the retarded advanced (RA) basis within the SK formalism [14,[29][30][31][32][33][34][35][36][37] which is the most convenient in making KMS conditions manifest. In terms of these combinations, we can then write Thus, we see that RA basis combinations appear naturally in holography :q i F sources the ingoing quasi-normal mode regular at the future horizon whereasq i P excites the Hawking mode.

A. Influence phase of the heavy quark
We will now compute the influence phase of the heavy quark in the retarded-advanced basis. To this end, it is useful to define the derivative combinations where Here A ω , B ω are analytic retarded functions which are well-behaved at the horizon with B ω vanishing at the horizon.
With these definitions, we can compute the influence phase of the heavy quark as The functions that appear in the above influence phase evaluate to These expressions describe the variety of four mode interactions on the world sheet which contribute to the influence phase. Here, the modes with −ω in their subscripts are the Hawking modes and those with +ω subscript are ingoing.
We note the remarkable fact that the terms of the form PP or FF do not appear in the quadratic influence phase, nor do the terms of the form PPPP or FFFF appear in the quartic influence phase. To understand this, we begin by looking at the PP and FF part of the quadratic influence functional which takes the form where we have used the identity e (1−χ)β ω = 1. Because of the cancellation of χ, the entire integrand here has the same value in both left and right contours. Further, the simple pole of at the horizon ρ 2 − 1 is cancelled by the zero of the B ω function which vanishes at the horizon. As a result, the integrand is in fact entirely analytic and the contributions from two contours cancel each other. A similar cancellation of contours lies behind the vanishing of PPPP or FFFF terms. The FFFF and PPPP part of the quartic influence functional is given by In this case, we again have a single valued function whose double pole dξ dρ 2 at the horizon is cancelled by the zeroes of the two B ω functions at the horizon. As a result, the integrand is entirely analytic and this contribution again cancels between the two contours.
The vanishing of the all P and all F terms is the hallmark of the retarded advanced basis in the Schwinger-Keldysh formalism [33][34][35][36][37]. On the CFT side, this is a combined consequence of both the largest time equation of SK formalism and the KMS conditions arising from the thermality of the bath. It is satisfying to see that this fact appears in our analysis based on general grounds of analyticity, independent of the underlying details.
On the gravity side, we see that all the contributions to the influence phase arise from interactions between the ingoing modes and the Hawking modes. Both the self-interaction between the ingoing modes as well as the selfinteraction effects of Hawking modes cancel out in the influence phase. This, perhaps, is an indication that on the CFT side the dominant physics we are looking at is the interaction between the quark and the CFT bath rather than effects arising from self-interactions. Further the influence phase, in general, receives contributions from the whole of the exterior spacetime (due to the branch cut in χ function along the radial direction). However, as we shall see below, at leading order in low frequency/derivative expansion, we get localised contributions from the horizon.

B. Influence phase in Keldysh basis and its derivative expansion
While the influence functional stated in terms of the retarded advanced basisq i P andq i F is very efficient at exhibiting the essential thermal structure of the influence phase, it is not convenient for studying the low frequency response. To this end, we will begin here by making the Bose-Einstein factors manifest so as to study the derivative expansion. This is achieved by substitutingq where in the expressions for the influence phase. Here q i a ≡ 1 2 (q i R + q i L ) and q i d ≡ q i R − q i L constitute the Keldysh or averagedifference basis [8,41,42]. As we did in the retarded-advanced basis, we will find it convenient to define where The influence phase in this basis takes the form S (2) where we have defined the coefficient functions in the previous subsection. The absence of the corresponding terms in the influence phase implies that the master equation that evolves the quark's density matrix has a Lindblad form [8,45].
We will now like to extract the low frequency behaviour from the full influence phase. Using the fact that f ω , N ω ∼ 1 βω , one can see that the low frequency expansions of the functions A ω , B ω , A ω and B ω take the following forms Substituting the above expressions in the quartic coefficient functions, we conclude that their low frequency behaviour is given by We will be interested in keeping terms till O(ω) which mean we should retain leading/sub-leading terms in G[d 4 ] and leading term in G[d 3 a] and ignore the rest.
We will now argue that the sub-leading O(ω) terms in G[d 4 ] do not contribute to the influence phase. The argument is as follows : given that G[d 4 ] multiplies a term permutation symmetric in ω's, i.e., q i d (ω 1 )q i a (ω 2 )q j d (ω 3 )q j d (ω 4 ), its contribution to the influence phase should also be permutation symmetric in ω's. But, the only permutation symmetric linear term in ω is ω 1 + ω 2 + ω 3 + ω 4 which vanishes due to the frequency conserving delta function. Thus, it suffices to examine the leading O(ω 0 ) term in G[d 4 ].
The leading order form is given by Since the integrand is analytic, we can perform this integral by picking up the residue from the double pole of dξ dρ 2 . We use the Laurent's expansion so that we can write We can similarly compute the leading order form of G[d 3 a] to get Thus, the derivative expansion for the quartic influence phase is S (4) We remind the reader that β is the inverse Hawking temperature of the black brane. The parameter λ is an effective CFT 't Hooft coupling related to the string tension via We notice that these results are consistent with the fluctuation-dissipation relations [52, 54-56, 59, 60] which connect γ and f 2 on one hand and ζ γ and ζ N on the other hand : We tabulate below the value of the mass corrections and Z I in various dimensions : Notice that the thermal mass correction vanishes in d = 2. Our results for d = 4 match with [17,18].

V. DISCUSSION AND CONCLUSION
In this work, we have studied the leading non-linear corrections to the Brownian motion of a heavy quark probing a strongly coupled CFT plasma. This is done using the holographic avatar of the Schwinger-Keldysh path integral applied to the world-sheet theory. Our computation relies on a doubled black-brane target space and a doubled string which probes it. This string stretches from the two AdS boundaries and loops around a region obtained by radial Wick rotation that connects the two stretched horizons. Evaluating the Nambu-Goto action on this string configuration results in the influence phase of the heavy quark.
We recognise this influence phase to be that of a particle with position q i , obeying a non-linear Langevin equation of the form Here N 2 ≡ N i N i and N i denotes the thermal noise which is drawn from a non-Gaussian distribution The parameters can be derived from holography and are given by Eqn.(97). We notice that these results are consistent with the fluctuation-dissipation relations [52, 54-56, 59, 60] which connect γ and f 2 on one hand and ζ γ and ζ N on the other hand : These two relations were derived earlier by studying weakly coupled baths [9,10]. It is satisfying to see that these relations continue to hold in holography. We remind the reader that both ζ γ and ζ N contributions to the influence phase of the heavy quark arise from the correct treatment of the near horizon region using the prescription of [28]. The fact that we get a IR finite contribution which is consistent with the non-linear FDT, is a nontrivial check of our construction.
The system under study permits many generalisations : one could study fermionic strings and the emergence of Fermi-Dirac distributions within holographic Brownian motion. Or one could look at strings in rotating blackhole backgrounds. The computation we do in this work can be generalised in different directions : one might ask whether one can replace Brownian strings with Brownian branes and derive their influence functional. One can imagine a sequence of Brownian branes probing the CFT plasma and their effective theories, with the space-filling Brownian brane essentially yielding the effective theory of the full plasma. This would give a way to check various features of such effective theories proposed in the literature [61]. An especially interesting problem is to study fluid-gravity correspondence within holographic SK formalism and see what it can say about the non-Gaussianities in fluid fluctuations [62][63][64][65][66][67]. Another direction would be to study field theories in the black-brane backgrounds and derive dual open quantum field theories for the corresponding single trace primary, by integrating out the effect of the black-brane. All these require extending holographic SK formalism to novel systems and would teach us more about the structure of real time correlations that emerge out of AdS/CFT.
Another obvious generalisation of our work would be to construct a string configuration which will allow us to extract the out-of-time-order influence phase of the heavy quark. This would require us to generalise our non-linear sigma model to a target space with many copies of the black brane smoothly connected to each other. The studies on weakly coupled quartic damped oscillator make definite predictions on the relation between the coefficients in the OTO effective theory vs those of the non-linear Langevin theory [10]. These relations arise from both microscopic time reversal invariance (generalised Onsager relations) and the thermality of the bath (KMS relations). It would be interesting to check these relations against a holographic computation. Getting the full OTO influence phase would also be a way to re-derive the result that the quark inherits the maximal Lyapunov exponent from the black-brane [68].
The structure of the influence we derive in this work also raises many deep conceptual questions that we do not have immediate answers to. We observe that a time reversal symmetry breaking order-parameter χ with a branch cut in its radial dependence plays a crucial role in our story. The CFT interpretation of this object is still unclear to us. Should we think about this as a series of UV modes in the bath which extend till thermal scale that contribute to the influence phase of the heavy quark ? Perhaps the branch cut is a large N remnant of a line of poles in CFT/bath correlators ? Another question is to trace where the entropy of the particle is created from ? Is there an entropy inflow from the complex spacetime that replaces the horizon in holographic SK contour [62,63,69] ? We see in this system the glimmer of many ideas which have recently appeared in the context of fluid SK effective theories [61].