Source author record

Sanjib Sabhapandit

Sanjib Sabhapandit appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2022arXiv

Direction reversing active Brownian particle in a harmonic potential

We study the two-dimensional motion of an active Brownian particle of speed $v_0$, with intermittent directional reversals in the presence of a harmonic trap of strength $μ$. The presence of the trap ensures that the position of the particle eventually reaches a steady state where it is bounded within a circular region of radius $v_0/μ$, centered at the minimum of the trap. Due to the interplay between the rotational diffusion constant $D_R$, reversal rate $γ$, and the trap strength $μ$, the steady state distribution shows four different types of shapes, which we refer to as active-I & II, and passive-I & II phases. In the active-I phase, the weight of the distribution is concentrated along an annular region close to the circular boundary, whereas in active-II, an additional central diverging peak appears giving rise to a Mexican hat-like shape of the distribution. The passive-I is marked by a single Boltzmann-like centrally peaked distribution in the large $D_R$ limit. On the other hand, while the passive-II phase also shows a single central peak, it is distinguished from passive-I by a non-Boltzmann like divergence near the origin. We characterize these phases by calculating the exact analytical forms of the distributions in various limiting cases. In particular, we show that for $D_R\llγ$, the shape transition of the two-dimensional position distribution from active-II to passive-II occurs at $μ=γ$. We compliment these analytical results with numerical simulations beyond the limiting cases and obtain a qualitative phase diagram in the $(D_R,γ,μ^{-1})$ space.

preprint2022arXiv

Integer partitions and exclusion statistics: Limit shapes and the largest part of Young diagrams

We compute the limit shapes of the Young diagrams of the minimal difference $p$ partitions and provide a simple physical interpretation for the limit shapes. We also calculate the asymptotic distribution of the largest part of the Young diagram and show that the scaled distribution has a Gumbel form for all $p$. This Gumbel statistics for the largest part remains unchanged even for general partitions of the form $E=\sum_i n_i i^{1/ν}$ with $ν>0$ where $n_i$ is the number of times the part $i$ appears.

preprint2022arXiv

Run-and-tumble particle in one-dimensional confining potential: Steady state, relaxation and first passage properties

We study the dynamics of a one-dimensional run and tumble particle subjected to confining potentials of the type $V(x) = α\, |x|^p$, with $p>0$. The noise that drives the particle dynamics is telegraphic and alternates between $\pm 1$ values. We show that the stationary probability density $P(x)$ has a rich behavior in the $(p, α)$-plane. For $p>1$, the distribution has a finite support in $[x_-,x_+]$ and there is a critical line $α_c(p)$ that separates an active-like phase for $α> α_c(p)$ where $P(x)$ diverges at $x_\pm$, from a passive-like phase for $α< α_c(p)$ where $P(x)$ vanishes at $x_\pm$. For $p<1$, the stationary density $P(x)$ collapses to a delta function at the origin, $P(x) = δ(x)$. In the marginal case $p=1$, we show that, for $α< α_c$, the stationary density $P(x)$ is a symmetric exponential, while for $α> α_c$, it again is a delta function $P(x) = δ(x)$. For the special cases $p=2$ and $p=1$, we obtain exactly the full time-dependent distribution $P(x,t)$, that allows us to study how the system relaxes to its stationary state. In addition, in these two cases, we also study analytically the full distribution of the first-passage time to the origin. Numerical simulations are in complete agreement with our analytical predictions.

preprint2022arXiv

Universal framework for the long-time position distribution of free active particles

Active particles self-propel themselves with a stochastically evolving velocity, generating a persistent motion leading to a non-diffusive behavior of the position distribution. Nevertheless, an effective diffusive behavior emerges at times much larger than the persistence time. Here we develop a general framework for studying the long-time behaviour for a class of active particle dynamics and illustrate it using the examples of run-and-tumble particle, active Ornstein-Uhlenbeck particle, active Brownian particle, and direction reversing active Brownian particle. Treating the ratio of the persistence-time to the observation time as the small parameter, we show that the position distribution generically satisfies the diffusion equation at the leading order. We further show that the sub-leading contributions, at each order, satisfies an inhomogeneous diffusion equation, where the source term depends on the previous order solutions. We explicitly obtain a few sub-leading contributions to the Gaussian position distribution. As a part of our framework, we also prescribe a way to find the position moments recursively and compute the first few explicitly for each model.

preprint2021arXiv

Record statistics for random walks and Lévy flights with resetting

We compute exactly the mean number of records $\langle R_N \rangle$ for a time-series of size $N$ whose entries represent the positions of a discrete time random walker on the line. At each time step, the walker jumps by a length $η$ drawn independently from a symmetric and continuous distribution $f(η)$ with probability $1-r$ (with $0\leq r < 1$) and with the complementary probability $r$ it resets to its starting point $x=0$. This is an exactly solvable example of a weakly correlated time-series that interpolates between a strongly correlated random walk series (for $r=0$) and an uncorrelated time-series (for $(1-r) \ll 1$). Remarkably, we found that for every fixed $r \in [0,1[$ and any $N$, the mean number of records $\langle R_N \rangle$ is completely universal, i.e., independent of the jump distribution $f(η)$. In particular, for large $N$, we show that $\langle R_N \rangle$ grows very slowly with increasing $N$ as $\langle R_N \rangle \approx (1/\sqrt{r})\, \ln N$ for $0<r <1$. We also computed the exact universal crossover scaling functions for $\langle R_N \rangle$ in the two limits $r \to 0$ and $r \to 1$. Our analytical predictions are in excellent agreement with numerical simulations.

preprint2020arXiv

Run-and-tumble particles in two-dimensions : Marginal position distributions

We study a set of Run-and-tumble particle (RTP) dynamics in two spatial dimensions. In the first case of the orientation θ of the particle can assume a set of n possible discrete values while in the second case θ is a continuous variable. We calculate exactly the marginal position distributions for n = 3,4 and the continuous case and show that in all the cases the RTP shows a cross-over from a ballistic to diffusive regime. The ballistic regime is a typical signature of the active nature of the systems and is characterized by non-trivial position distributions which depends on the specific model. We also show that, the signature of activity at long-times can be found in the atypical fluctuations which we also characterize by computing the large deviation functions explicitly.

preprint2019arXiv

Entropy production for partially observed harmonic systems

The probability distribution of the total entropy production in the non-equilibrium steady state follows a symmetry relation called the fluctuation theorem. When a certain part of the system is masked or hidden, it is difficult to infer the exact estimate of the total entropy production. Entropy produced from the observed part of the system shows a significant deviation from the steady state fluctuation theorem. This deviation occurs due to the interaction between the observed and the masked part of the system. A naive guess would be that the deviation from the steady state fluctuation theorem may disappear in the limit of small interaction between both parts of the system. In contrast, we investigate the entropy production of a particle in a harmonically coupled Brownian particle system (say, particle A and B) in a heat reservoir at a constant temperature. The system is maintained in the non-equilibrium steady state using stochastic driving. When the coupling between particle A and B is infinitesimally weak, the deviation from the steady state fluctuation theorem for the entropy production of a partial system of a coupled system is studied. Furthermore, we consider a harmonically confined system (i.e., a harmonically coupled system of particle A and B in harmonic confinement). In the weak coupling limit, the entropy produced by the partial system (e.g., particle A) of the coupled system in a harmonic trap satisfies the steady state fluctuation theorem. Numerical simulations are performed to support the analytical results. Part of these results was announced in a recent letter, Europhys. Lett. 115, 60003 (2016).

preprint2019arXiv

Exact stationary state of a run-and-tumble particle with three internal states in a harmonic trap

We study the motion of a one-dimensional run-and-tumble particle with three discrete internal states in the presence of a harmonic trap of stiffness $μ.$ The three internal states, corresponding to positive, negative and zero velocities respectively, evolve following a jump process with rate $γ$. We compute the stationary position distribution exactly for arbitrary values of $μ$ and $γ$ which turns out to have a finite support on the real line. We show that the distribution undergoes a shape-transition as $β=γ/μ$ is changed. For $β<1,$ the distribution has a double-concave shape and shows algebraic divergences with an exponent $(β-1)$ both at the origin and at the boundaries. For $β>1,$ the position distribution becomes convex, vanishing at the boundaries and with a single, finite, peak at the origin. We also show that for the special case $β=1,$ the distribution shows a logarithmic divergence near the origin while saturating to a constant value at the boundaries.

preprint2016arXiv

Exact distributions of cover times for $N$ independent random walkers in one dimension

We study the probability density function (PDF) of the cover time $t_c$ of a finite interval of size $L$, by $N$ independent one-dimensional Brownian motions, each with diffusion constant $D$. The cover time $t_c$ is the minimum time needed such that each point of the entire interval is visited by at least one of the $N$ walkers. We derive exact results for the full PDF of $t_c$ for arbitrary $N \geq 1$, for both reflecting and periodic boundary conditions. The PDFs depend explicitly on $N$ and on the boundary conditions. In the limit of large $N$, we show that $t_c$ approaches its average value $\langle t_c \rangle \approx L^2/(16\, D \, \ln N)$, with fluctuations vanishing as $1/(\ln N)^2$. We also compute the centered and scaled limiting distributions for large $N$ for both boundary conditions and show that they are given by nontrivial $N$-independent scaling functions.

preprint2016arXiv

Fluctuation theorem for entropy production of a partial system in the weak coupling limit

Small systems in contact with a heat bath evolve by stochastic dynamics. Here we show that, when one such small system is weakly coupled to another one, it is possible to infer the presence of such weak coupling by observing the violation of the steady state fluctuation theorem for the partial entropy production of the observed system. We give a general mechanism due to which the violation of the fluctuation theorem can be significant, even for weak coupling. We analytically demonstrate on a realistic model system that this mechanism can be realized by applying an external random force to the system. In other words, we find a new fluctuation theorem for the entropy production of a partial system, in the limit of weak coupling.

preprint2015arXiv

Driven inelastic Maxwell gases

We consider the inelastic Maxwell model, which consists of a collection of particles that are characterized by only their velocities, and evolving through binary collisions and external driving. At any instant, a particle is equally likely to collide with any of the remaining particles. The system evolves in continuous time with mutual collisions and driving taken to be point processes with rates $τ_c^-{1}$ and $τ_w^{-1}$ respectively. The mutual collisions conserve momentum and are inelastic, with a coefficient of restitution $r$. The velocity change of a particle with velocity $v$, due to driving, is taken to be $Δv=-(1+r_w) v+η$, mimicking the collision with a vibrating wall, where $r_w$ the coefficient of restitution of the particle with the "wall" and $η$ is Gaussian white noise. The Ornstein-Uhlenbeck driving mechanism given by $\frac{dv}{dt}=-Γv+η$ is found to be a special case of the driving modeled as a point process. Using both the continuum and discrete versions we show that while the equations for the one-particle and the two-particle velocity distribution functions do not close, the joint evolution equations of the variance and the two-particle velocity correlation functions close. With the exact formula for the variance we find that, for $r_w\ne-1$, the system goes to a steady state. On the other hand, for $r_w=-1$, the system does not have a steady state. Similarly, the system goes to a steady state for the Ornstein-Uhlenbeck driving with $Γ\not=0$, whereas for the purely diffusive driving ($Γ=0$), the system does not have a steady state.

preprint2015arXiv

Dynamical transition in the temporal relaxation of stochastic processes under resetting

A stochastic process, when subject to resetting to its initial condition at a constant rate, generically reaches a non-equilibrium steady state. We study analytically how the steady state is approached in time and find an unusual relaxation mechanism in these systems. We show that as time progresses, an inner core region around the resetting point reaches the steady state, while the region outside the core is still transient. The boundaries of the core region grow with time as power laws at late times. Alternatively, at a fixed spatial point, the system undergoes a dynamical transition from the transient to the steady state at a characteristic space dependent timescale $t^*(x)$. We calculate analytically in several examples the large deviation function associated with this spatio-temporal fluctuation and show that generically it has a second order discontinuity at a pair of critical points characterizing the edges of the inner core. Our results are verified in the numerical simulations of several models, such as simple diffusion and fluctuating one-dimensional interfaces.

preprint2015arXiv

Exact probability distribution for the two-tag displacement in single-file motion

We consider a gas of point particles moving on the one-dimensional line with a hard-core inter-particle interaction that prevents particle crossings --- this is usually referred to as single-file motion. The individual particle dynamics can be arbitrary and they only interact when they meet. Starting from initial conditions such that particles are uniformly distributed, we observe the displacement of a tagged particle at time $t$, with respect to the initial position of another tagged particle, such that their tags differ by $r$. For $r=0$, this is the usual well studied problem of the tagged particle motion. Using a mapping to a non-interacting particle system we compute the exact probability distribution function for the two-tagged particle displacement, for general single particle dynamics. As by-products, we compute the large deviation function, various cumulants and, for the case of Hamiltonian dynamics, the two-particle velocity auto-correlation function.

preprint2015arXiv

Large deviations for Markov processes with resetting

Markov processes restarted or reset at random times to a fixed state or region in space have been actively studied recently in connection with random searches, foraging, and population dynamics. Here we study the large deviations of time-additive functions or observables of Markov processes with resetting. By deriving a renewal formula linking generating functions with and without resetting we are able to obtain the rate function of such observables, characterizing the likelihood of their fluctuations in the long-time limit. We consider as an illustration the large deviations of the area of the Ornstein-Uhlenbeck process with resetting. Other applications involving diffusions, random walks, and jump processes with resetting or catastrophes are discussed.

preprint2015arXiv

Random walk with random resetting to the maximum

We study analytically a simple random walk model on a one-dimensional lattice, where at each time step the walker resets to the maximum of the already visited positions (to the rightmost visited site) with a probability $r$, and with probability $(1-r)$, it undergoes symmetric random walk, i.e., it hops to one of its neighboring sites, with equal probability $(1-r)/2$. For $r=0$, it reduces to a standard random walk whose typical distance grows as $\sqrt{n}$ for large $n$. In presence of a nonzero resetting rate $0<r\le 1$, we find that both the average maximum and the average position grow ballistically for large $n$, with a common speed $v(r)$. Moreover, the fluctuations around their respective averages grow diffusively, again with the same diffusion coefficient $D(r)$. We compute $v(r)$ and $D(r)$ explicitly. We also show that the probability distribution of the difference between the maximum and the location of the walker, becomes stationary as $n\to \infty$. However, the approach to this stationary distribution is accompanied by a dynamical phase transition, characterized by a weakly singular large deviation function. We also show that $r=0$ is a special `critical' point, for which the growth laws are different from the $r\to 0$ case and we calculate the exact crossover functions that interpolate between the critical $(r=0)$ and the off-critical $(r\to 0)$ behavior for finite but large $n$.

preprint2015arXiv

Tagged particle diffusion in one-dimensional systems with Hamiltonian dynamics - II

We study various temporal correlation functions of a tagged particle in one-dimensional systems of interacting point particles evolving with Hamiltonian dynamics. Initial conditions of the particles are chosen from the canonical thermal distribution. The correlation functions are studied in finite systems, and their forms examined at short and long times. Various one-dimensional systems are studied. Results of numerical simulations for the Fermi-Pasta-Ulam chain are qualitatively similar to results for the harmonic chain, and agree unexpectedly well with a simple description in terms of linearized equations for damped fluctuating sound waves. Simulation results for the alternate mass hard particle gas reveal that - in contradiction to our earlier results [1] with smaller system sizes - the diffusion constant slowly converges to a constant value, in a manner consistent with mode coupling theories. Our simulations also show that the behaviour of the Lennard-Jones gas depends on its density. At low densities, it behaves like a hard-particle gas, and at high densities like an anharmonic chain. In all the systems studied, the tagged particle was found to show normal diffusion asymptotically, with convergence times depending on the system under study. Finite size effects show up at time scales larger than sound traversal times, their nature being system-specific.

preprint2014arXiv

First order transition for the optimal search time of Lévy flights with resetting

We study analytically an intermittent search process in one dimension. There is an immobile target at the origin and a searcher undergoes a discrete time jump process starting at $x_0\geq0$, where successive jumps are drawn independently from an arbitrary jump distribution $f(η)$. In addition, with a probability $0\leq r \leq1$ the position of the searcher is reset to its initial position $x_0$. The efficiency of the search strategy is characterized by the mean time to find the target, i.e., the mean first passage time (MFPT) to the origin. For arbitrary jump distribution $f(η)$, initial position $x_0$ and resetting probability $r$, we compute analytically the MFPT. For the heavy-tailed Lévy stable jump distribution characterized by the Lévy index $0<μ< 2$, we show that, for any given $x_0$, the MFPT has a global minimum in the $(μ,r)$ plane at $(μ^*(x_0),r^*(x_0))$. We find a remarkable first-order phase transition as $x_0$ crosses a critical value $x_0^*$ at which the optimal parameters change discontinuously. Our analytical results are in good agreement with numerical simulations.

preprint2014arXiv

High-Energy Tail of the Velocity Distribution of Driven Inelastic Maxwell Gases

A model of homogeneously driven dissipative system, consisting of a collection of $N$ particles that are characterized by only their velocities, is considered. Adopting a discrete time dynamics, at each time step, a pair of velocities is randomly selected. They undergo inelastic collision with probability $p$. With probability $(1-p)$, energy of the system is changed by changing the velocities of both the particles independently according to $v\rightarrow -r_w v +η$, where $η$ is a Gaussian noise drawn independently for each particle as well as at each time steps. For the case $r_w=- 1$, although the energy of the system seems to saturate (indicating a steady state) after time steps of $O(N)$, it grows linearly with time after time steps of $O(N^2)$, indicating the absence of a eventual steady state. For $ -1 <r_w \leq 1$, the system reaches a steady state, where the average energy per particle and the correlation of velocities are obtained exactly. In the thermodynamic limit of large $N$, an exact equation is obtained for the moment generating function. In the limit of nearly elastic collisions and weak energy injection, the velocity distribution is shown to be a Gaussian. Otherwise, for $|r_w| < 1$, the high-energy tail of the velocity distribution is Gaussian, with a different variance, while for $r_w=+1$ the velocity distribution has an exponential tail.

preprint2014arXiv

Universal Large Deviations for the Tagged Particle in Single File Motion

We consider a gas of point particles moving in a one-dimensional channel with a hard-core inter-particle interaction that prevents particle crossings --- this is called single-file motion. Starting from equilibrium initial conditions we observe the motion of a tagged particle. It is well known that if the individual particle dynamics is diffusive, then the tagged particle motion is sub-diffusive, while for ballistic particle dynamics, the tagged particle motion is diffusive. Here we compute exactly the large deviation function for the tagged particle displacement and show that this is universal, independent of the individual dynamics.

preprint2014arXiv

Work fluctuations for a Brownian particle driven by a correlated external random force

We have considered the underdamped motion of a Brownian particle in the presence of a correlated external random force. The force is modeled by an Ornstein-Uhlenbeck process. We investigate the fluctuations of the work done by the external force on the Brownian particle in a given time interval in the steady state. We calculate the large deviation functions as well as the complete asymptotic form of the probability density function of the performed work. We also discuss the symmetry properties of the large deviation functions for this system. Finally we perform numerical simulations and they are in a very good agreement with the analytic results.

preprint2013arXiv

Work fluctuations for a Brownian particle in a harmonic trap with fluctuating locations

We consider a Brownian particle in a harmonic trap. The location of the trap is modulated according to an Ornstein-Uhlenbeck process. We investigate the fluctuation of the work done by the modulated trap on the Brownian particle in a given time interval in the steady state. We compute the large deviation as well as the complete asymptotic form of the probability density function of the work done. The theoretical asymptotic forms of the probability density function are in very good agreement with the numerics. We also discuss the validity of the fluctuation theorem for this system.

preprint2012arXiv

Heat and work fluctuations for a harmonic oscillator

The formalism of Kundu et al. [J. Stat. Mech. (2011) P03007], for computing the large deviations of heat flow in harmonic systems, is applied to the case of single Brownian particle in a harmonic trap and coupled to two heat baths at different temperatures. The large-t form of the moment generating function <exp[-sQ]> ~ g(s) exp[t m(s)], of the total heat flow Q from one of the baths to the particle in a given time interval t, is studied and exact explicit expressions are obtained for both m(s) and g(s). For a special case of the single particle problem that corresponds to the work done by an external stochastic force on a harmonic oscillator coupled to a thermal bath, the large-t form of the moment generating function is analyzed to obtain the exact large deviation function as well as the complete asymptotic forms of the probability density function of the work.

preprint2012arXiv

Scaling behavior in the convection-driven Brazil-nut effect

The Brazil-nut effect is the phenomenon in which a large intruder particle immersed in a vertically shaken bed of smaller particles rises to the top, even when it is much denser. The usual practice, while describing these experiments, has been to use the dimensionless acceleration Γ=a ω^2/g, where a and ωare respectively the amplitude and the angular frequency of vibration and g is the acceleration due to gravity. Considering a vibrated quasi-two-dimensional bed of mustard seeds, we show here that the peak-to-peak velocity of shaking v= aω, rather than Γ, is the relevant parameter in the regime where boundary-driven granular convection is the main driving mechanism. We find that the rise-time τof an intruder is described by the scaling law τ~ (v-v_c)^{-α}, where v_c is identified as the critical vibration velocity for the onset of convective motion of the mustard seeds. This scaling form holds over a wide range of (a,ω), diameter and density of the intruder.

preprint2012arXiv

Tagged particle diffusion in one-dimensional gas with Hamiltonian dynamics

We consider a one-dimensional gas of hard point particles in a finite box that are in thermal equilibrium and evolving under Hamiltonian dynamics. Tagged particle correlation functions of the middle particle are studied. For the special case where all particles have the same mass, we obtain analytic results for the velocity auto-correlation function in the short time diffusive regime and the long time approach to the saturation value when finite-size effects become relevant. In the case where the masses are unequal, numerical simulations indicate sub-diffusive behaviour with mean square displacement of the tagged particle growing as t/ln(t) with time t. Also various correlation functions, involving the velocity and position of the tagged particle, show damped oscillations at long times that are absent for the equal mass case.

preprint2011arXiv

Application of importance sampling to the computation of large deviations in non-equilibrium processes

We present an algorithm for finding the probabilities of rare events in nonequilibrium processes. The algorithm consists of evolving the system with a modified dynamics for which the required event occurs more frequently. By keeping track of the relative weight of phase-space trajectories generated by the modified and the original dynamics one can obtain the required probabilities. The algorithm is tested on two model systems of steady-state particle and heat transport where we find a huge improvement from direct simulation methods.

preprint2011arXiv

Dynamics of a flexible polymer in planar mixed flow

We present exact spatio-temporal correlation functions of a Rouse polymer chain submerged in a fluid having planar mixed flow, in the steady state. Using these correlators, determination of the time scale distribution functions associated with the first-passage tumbling events is difficult in general; it was done recently in Phys. Rev. Lett.,101, 188301 (2008), for the special case of "simple shear" flow. We show here that the method used in latter paper fails for the general mixed flow problem. We also give many new estimates of the exponent $θ$ associated with the exponential tail of the angular tumbling time distribution in the case of simple shear

preprint2011arXiv

Large deviations of heat flow in harmonic chains

We consider heat transport across a harmonic chain connected at its two ends to white-noise Langevin reservoirs at different temperatures. In the steady state of this system the heat $Q$ flowing from one reservoir into the system in a finite time $τ$ has a distribution $P(Q,τ)$. We study the large time form of the corresponding moment generating function $<e^{-λQ}>\sim g(λ) e^{τμ(λ)}$. Exact formal expressions, in terms of phonon Green's functions, are obtained for both $μ(λ)$ and also the lowest order correction $g(λ)$. We point out that, in general a knowledge of both $μ(λ)$ and $g(λ)$ is required for finding the large deviation function associated with $P(Q,τ)$. The function $μ(λ)$ is known to be the largest eigenvector of an appropriate Fokker-Planck type operator and our method also gives the corresponding eigenvector exactly.

preprint2011arXiv

Record Statistics of Continuous Time Random Walk

The statistics of records for a time series generated by a continuous time random walk is studied, and found to be independent of the details of the jump length distribution, as long as the latter is continuous and symmetric. However, the statistics depend crucially on the nature of the waiting time distribution. The probability of finding M records within a given time duration t, for large t, has a scaling form, and the exact scaling function is obtained in terms of the one-sided Levy stable law. The mean of the ages of the records, defined as <t/M>, differs from t/<M>. The asymptotic behaviour of the shortest and the longest ages of the records are also studied.

preprint2011arXiv

Work fluctuations for a harmonic oscillator driven by an external random force

The fluctuations of the work done by an external Gaussian random force on a harmonic oscillator that is also in contact with a thermal bath is studied. We have obtained the exact large deviation function as well as the complete asymptotic forms of the probability density function. The distribution of the work done are found to be non-Gaussian. The steady state fluctuation theorem holds only if the ratio of the variances, of the external random forcing and the thermal noise respectively, is less than 1/3. On the other hand, the transient fluctuation theorem holds (asymptotically) for all the values of that ratio. The theoretical asymptotic forms of the probability density function are in very good agreement with the numerics as well as with an experiment.