Consider the "short time" relaxation of a free Brownian particle. The Langevin equation for the velocity is \(m \mathrm{~d} v / \mathrm{d} t=-\gamma v+\xi(t)\). (a) Find the Fokker-Planck equation for the probability, \(P(v, t) \mathrm{d} v\), to find the Brownian particle with velocity, \(v \rightarrow v+\mathrm{d} v\), at time \(t\). (b) Solve the FokkerPlanck equation, assuming that at time, \(t=0\), the velocity is \(v=v_{0}\). (Hint: Use the transform in Section Z.3.3.1 to write the Fokker-Planck in terms of a Hermitian operator. The eigenfunctions of that operator will be Hermite polynomials.)

Short Answer

Expert verified
The Fokker-Planck equation is approximately solved using Hermitian operators and Hermite polynomials to yield the distribution, \( P(v,t) \).

Step by step solution

01

Start from the Langevin Equation

The given Langevin equation describing the velocity of a Brownian particle is: \[ m \frac{\mathrm{d}v}{\mathrm{d}t} = -\gamma v + \xi(t). \]
02

Rewrite the Langevin Equation

Divide both sides by the mass \(m\) to simplify: \[ \frac{\mathrm{d}v}{\mathrm{d}t} = -\frac{\gamma}{m} v + \frac{1}{m} \xi(t). \]
03

Identify Coefficients for the Fokker-Planck Equation

Recall that the general form of the Langevin equation corresponds to a Fokker-Planck equation with drift coefficient \( A_v(v) = -\frac{\gamma}{m} v \) and diffusion coefficient \( B_v(v) = \frac{1}{m} \xi(t) \).
04

Write the Fokker-Planck Equation

Using the identified coefficients, the Fokker-Planck equation for the probability distribution \( P(v,t) \) is given by: \[ \frac{\partial P(v,t)}{\partial t} = \frac{\partial}{\partial v} \left( \frac{\gamma}{m} v P(v,t) \right) + \frac{1}{2} \frac{\partial^2}{\partial v^2} \left( \frac{1}{m^2} \xi(t)^2 P(v,t) \right). \]
05

Set Initial Condition

Assume that at time \( t=0 \), the velocity \( v = v_0 \). This sets up the initial condition for the probability distribution: \[ P(v,0) = \delta(v - v_0). \]
06

Simplify the Fokker-Planck Equation

For short time relaxation, the noise term \( \xi(t) \) can be approximated to simplify the equation. Rewrite the Fokker-Planck equation based on this assumption.
07

Transform Using Hermitian Operator

Apply the suggested transform to write the Fokker-Planck in terms of a Hermitian operator. Use the eigenfunctions of that operator, which are Hermite polynomials, to solve the equation.
08

Solve for the Probability Distribution

Solve the resulting simplified Fokker-Planck equation using the Hermite polynomials to find: \[ P(v,t) = ... \] (write the detailed solution here).

Unlock Step-by-Step Solutions & Ace Your Exams!

  • Full Textbook Solutions

    Get detailed explanations and key concepts

  • Unlimited Al creation

    Al flashcards, explanations, exams and more...

  • Ads-free access

    To over 500 millions flashcards

  • Money-back guarantee

    We refund you if you fail your exam.

Over 30 million students worldwide already upgrade their learning with Vaia!

Key Concepts

These are the key concepts you need to understand to accurately answer the question.

Fokker-Planck Equation
The Fokker-Planck equation is key for describing the time evolution of the probability distribution of a system's state variables. In this context, it details how the probability density function of the velocity of a Brownian particle changes over time.

Given the Langevin equation, the first step is to determine the drift and diffusion coefficients. These coefficients help translate the stochastic differential equation into the deterministic partial differential equation known as the Fokker-Planck equation.

The drift coefficient represents the systematic part of the motion (e.g., friction), while the diffusion coefficient encapsulates the random forces acting on the particle. For our exercise:
  • The drift coefficient is given by \(A_v(v) = -\frac{\gamma}{m} v\).
  • The diffusion coefficient derived from the noise term \(\xi(t)\) is \(B_v(v) = \frac{1}{m} \xi(t)\).
Using these coefficients, we form the Fokker-Planck equation:
\[ \frac{\partial P(v,t)}{\partial t} = \frac{\partial}{\partial v} \left( \frac{\gamma}{m} v P(v,t) \right) + \frac{1}{2} \frac{\partial^2}{\partial v^2} \left( \frac{1}{m^2} \xi(t)^2 P(v,t) \right). \] Finally, this equation describes how the probability of finding a particle with velocity v at time t evolves. The initial condition is usually a delta function if the initial velocity is known exactly.
Brownian Motion
Brownian motion, named after botanist Robert Brown, refers to the random motion of particles suspended in a fluid. This motion results from their collision with fast-moving molecules in the fluid. In physics, it is a fundamental example of a stochastic process, characterized by random and unpredictable changes in a particle's position and velocity over time.

Analyzing a Brownian particle involves considering both deterministic forces (like friction) and random forces (stemming from molecular impacts).
  • Deterministic forces often lead to a drift, modeled as a function of velocity.
  • Random forces introduce a diffusion term, accounting for the spread in velocity due to stochastic impacts.
The Langevin equation describes the velocity of a Brownian particle, incorporating both these forces together. Switching from the Langevin description to a statistical description using the Fokker-Planck equation allows us to study the overall behavior of many particles over time.
Understanding Brownian motion is crucial across various scientific fields. In finance, it's used to model stock prices; in biology, to track the diffusion of molecules in cells; and in physics, it lays the groundwork for the study of stochastic processes.
Hermite Polynomials
Hermite polynomials arise in probability theory, particularly when dealing with Gaussian distributions and solutions to differential equations with Gaussian noise. In the context of the Langevin and Fokker-Planck equations, Hermite polynomials facilitate solving these equations through spectral methods.

Hermite polynomials, \(H_n(x)\), are orthogonal polynomials over the weight function \(e^{-x^2}\), making them suitable for problems involving Gaussian-like processes:
  • The zeroth Hermite polynomial is constant: \(H_0(x) = 1 \).
  • The first few polynomials include \(H_1(x) = 2x\) and \(H_2(x) = 4x^2 - 2\).
  • Orthogonality condition: \[ \int_{-\infty}^{\infty} H_m(x) H_n(x) e^{-x^2} dx = 0 \quad \text{for} \ \ m e n \]
In our specific problem, the Hint suggests transforming the Fokker-Planck equation into a form where these Hermite polynomials are the eigenfunctions of a Hermitian operator. This transformation reduces the original problem into simpler parts. The final probability distribution, \(P(v, t)\), is then often expressed in terms of these polynomials.
Hermite polynomials thus form a powerful toolset for tackling complex differential equations by converting them into solvable algebraic equations.

One App. One Place for Learning.

All the tools & learning materials you need for study success - in one app.

Get started for free

Most popular questions from this chapter

Consider a casual function, \(g(t)\), such that \(g(t)=0\) for \(t<0\) and \(g(t) \neq 0\) for \(t>0\). Let \(\bar{g}_{\mathrm{L}}(z)=\int_{0}^{\infty} \mathrm{d} t \cdot \mathrm{e}^{-z t} g(t)\) be the Laplace transform. Let \(\bar{g}_{\mathrm{F}}(\omega)=\int_{-\infty}^{\infty} \mathrm{d} t \mathrm{e}^{\mathrm{i} \omega t} g(t)\) be the Fourier transform. Show that $$ \bar{g}_{\mathrm{L}}(z)=\int_{\infty}^{\infty} \frac{\mathrm{d} \omega}{2 \pi} \frac{\tilde{g}_{\mathrm{F}}(\omega)}{z+\mathrm{i} \omega}=\bar{g}_{\mathrm{F}}(\mathrm{i} z) $$ and show that the inverse transforms give the correct behavior for \(g(t)\).

A paramagnetic system sits in a constant external field \(H_{0} z\), and therefore has constant magnetization, \(M_{0} \hat{z}\). Assume that this is the equilibrium configuration of the system and that a small uniform time- dependent magnetic field, \(H_{1}(t)=H_{1} \cos \left(\omega_{0} t\right) \hat{x}\) is applied to the system. The magnetization then becomes time- dependent. Its equation of motion can be written $$ \frac{\mathrm{d}(\boldsymbol{M}\rangle}{\mathrm{d} t}=-\gamma\langle\boldsymbol{M}\rangle \times H(t)-D(\langle\boldsymbol{M}\rangle) $$ (these are called Bloch equations), where \((M\rangle=\operatorname{Tr}(\hat{\rho} M)=M_{x}(t) \hat{e}_{x}+M_{y}(t) \hat{e}_{y}+\) \(\left(M_{0}-M_{z}(t)\right) \hat{e}_{z}, H(t)=H_{0} \hat{e}_{z}+\vec{H}_{1}(t) \hat{e}_{x}\), and \(D(\quad M)\) is a damping term due to interactions between particles in the medium. The constant \(y\) equals \(g \mu / h\), where \(g\) is the Lande \(g\)-factor and \(\mu\) is the magnetic moment of particles in the system. (a) Write the equations of motion for \(M_{x}(t), M_{y}(t)\), and \(M_{z}(t)\), assuming \(D_{x}=M_{x} / T_{2}, D_{y}=M_{y} / T_{2}\), and \(D_{z}=M_{z} / T_{1}\). The equations you derive will be nonlinear functions of \(H_{1}, M_{x}, M_{y}\), and \(M_{z}\). (b) Assume that \(H_{1}, M_{x}, M_{y}\), and \(M_{z}\) are small, so the equations can be linearized. If we define the response matrix, \(\bar{K}(t)\), from the equation \(\langle M(t)\rangle=\int_{-\infty}^{\infty} \mathrm{d} t^{r} \bar{K}\left(t-t^{\prime}\right) \cdot H\left(t^{\prime}\right)\), find \(X_{x x}(\omega)\) and \(K_{x x}(t) .\) Write the expression for \(\chi_{x x}(\omega)\) in the limits \(1 / T_{1}\) \(\rightarrow 0\) and \(1 / T_{2} \rightarrow 0 .\) (c) Compute the average power absorbed during the period \(T=2 \pi / \omega_{0}\). (Note that \(\hat{e}_{x} \hat{e}_{y \text { 'and }} \hat{e}_{z \text { are unit vectors in the } x-,} y\)-, and \(z\)-directions, respectively.)

Consider a 1D ballistic waveguide which is coupled to low-temperature electron reservoirs with Fermi distributions \(f_{\mathrm{L}}\left(E_{n}\right)\) and \(f_{\mathrm{R}}\left(E_{n}\right)\) on the left and right ends, respectively. Assume that the temperatures of the reservoirs are the same but that the left reservoir has an electric potential energy bias \(-e V\), where \(e\) is the electron charge. Assume the waveguide has a length \(L\) that is very large so you can take the limit \(L \rightarrow \infty\), and keep only the lowest-order term in \(V\). Compute the average current in the waveguide and show that the conductivity of the waveguide is given by \(G=e^{2} / h\), where \(h\) is Planck's constant.

Prove that $$ \langle\alpha(t)\rangle=2 \mathrm{i} \int_{-\infty}^{t} \mathrm{~d} t^{\prime} \bar{K}^{\prime \prime}\left(t-t^{\prime}\right) \cdot \bar{F}\left(t^{\prime}\right) $$ where \(\bar{K}^{\prime \prime}(t)\) is the Fourier transform of \(\bar{X}^{\prime \prime}(\omega) .\) This form of the linear response is commonly seen in the literature. [Note that the identity $$ \lim _{\eta \rightarrow 0} \frac{1}{\omega^{\prime}-\omega \mp \mathrm{i} \eta}=P \frac{1}{\omega^{\prime}-\omega} \pm \mathrm{i} \pi \delta\left(\omega^{\prime}-\omega\right) $$ and the spectral representation of the Heaviside function $$ \theta\left(t-t^{\prime}\right)=-\lim _{\eta \rightarrow 0} \int_{-\infty}^{\infty} \frac{\mathrm{d} \omega}{(2 \pi \mathrm{i})} \frac{\mathrm{e}^{-\mathrm{i} \omega\left(t-t^{\prime}\right)}}{\omega+i \eta} $$ are useful.]

The magnetization \(M(t)\) of a spin system is driven out of equilibrium by a magnetic field \(H(t)\) and obeys the Langevin equation of motion $$ \frac{\mathrm{d} M(t)}{\mathrm{d} t}+\Gamma M(t)=\xi(t)-\Gamma_{X_{0}} H(t) $$ where \(\mathcal{\xi}(t)\) is a delta-correlated white noise with zero mean that is due to background magnetic field fluctuations and \(\Gamma\) is a positive constant. (a) Compute the linear response function, \(K(t)\), which is defined through the equation \(\langle M(t)\rangle_{\varepsilon}=\int_{-\infty}^{\infty} \mathrm{d} t^{\prime} K\left(t-t^{\prime}\right) H\left(t^{\prime}\right),(M(t))_{\xi}=\) (b) Compute the equilibrium dynamic correlation function \((M(t) M(0)\rangle_{e q}+(c)\) What is \(\chi_{0}\) ?

See all solutions

Recommended explanations on Physics Textbooks

View all explanations

What do you think about this solution?

We value your feedback to improve our textbook solutions.

Study anywhere. Anytime. Across all devices.

Sign-up for free