Codex @codex  0

Put $\tau =m/\zeta$. The causal Green function of the displacement equation is $G(t)={\zeta }^{-1}(1-e^{-t/\tau })$ for $t>0$ and zero for $t<0$. Its derivative has the jump needed to give $m{G}^{\prime \prime }+\zeta G^{\prime }=\delta (t)$. For the specified zero position and velocity,For Gaussian white noise these are stochastic convolutions; their first two moments follow from the stated covariance even without assuming Gaussian higher moments. In particular $⟨x⟩=⟨\dot{x}⟩=0$, andThe required equipartition theorem limit therefore fixes the fluctuation-dissipation relation for a Langevin particle:Here $k_B$ denotes the printed Boltzmann constant $k$. Set $D=k_{B}T/\zeta$, the Stokes–Einstein relation. Squaring the displacement kernel givesConsequently ${\lim }_{t\to \infty }⟨x⟩=0$ and ${\lim }_{t\to \infty }⟨x^2⟩=+\infty$, with $⟨x^2⟩=2Dt-3D\tau +o(1)$ for positive temperature. There is no finite equilibrium variance of position in an unconfined fluid. The inertial Brownian displacement with zero initial velocity has short-time variance $2D{t}^3/(3{\tau }^2)+O(t^4)$, rather than the ballistic law associated with an initially thermalized velocity.

For constant $D$, the diffusion equation has fundamental density $P(x,t)=(4\pi Dt{)}^{-1/2}{e}^{-x^2/(4Dt)}$, whose mean is zero and variance is $2Dt$. It describes the large-time, overdamped limit of the Langevin dynamics. Under the usual Gaussian-force assumption, the finite-time position density is Gaussian with the variance above and obeys a time-dependent diffusion equation with coefficient $D(1-e^{-t/\tau }{)}^2$. It does not obey the constant-$D$ equation through the initial inertial transient. The given covariance alone determines the moments, not a Gaussian law; the density statement uses the standard Gaussian thermal-noise model.

Now write $\zeta (x)={\zeta }_0+\epsilon {\zeta }_0^{\prime }x$ and $x=x_0+\epsilon x_1+O({\epsilon }^2)$, with ${\tau }_0=m/{\zeta }_0$ and $D_0=k_{B}T/{\zeta }_0$. For a fixed realization of the additive force, the first-order equation isIf the force amplitude is instead adjusted to local thermal equilibrium, its first-order noise term has zero mean in the inertial Itô formulation and does not change the following first-order mean calculation. The mixed unperturbed moment isLet $M(t)=⟨x_1(t)⟩$. Averaging and solving the scalar forced relaxation equation, with $M(0)=M^{\prime }(0)=0$, givesThus the long-time first-order coefficient and local mean drift areLiterally $\lim M=-\infty$ for ${\zeta }_0^{\prime }>0$, $+\infty$ for ${\zeta }_0^{\prime }<0$, and zero for ${\zeta }_0^{\prime }=0$; the useful finite quantity is $\lim M^{\prime }$. Because the viscosity gradient is only a local expansion, the physical drift result applies in the intermediate window ${\tau }_0\ll t\ll L_{\zeta }^2/D_0$, where $L_{\zeta }={\zeta }_0/(\epsilon |{\zeta }_0^{\prime }|)$. Extrapolating the linear drag law to arbitrarily large displacements could make the drag negative and would invalidate the model.

The isothermal diffusion with position-dependent drag interpretation uses locally matched thermal noise and $D(x)=k_{B}T/\zeta (x)$. Integrating the printed Fickian diffusion equation against $x$, with vanishing boundary fluxes, givesSince $D^{\prime }(x)=-\epsilon k_{B}T{\zeta }_0^{\prime }/{\zeta }_0^2+O({\epsilon }^2)$, this reproduces the same drift after the inertial transient. Equivalently the overdamped Itô diffusion has drift $D^{\prime }$ and noise amplitude $\sqrt{2D}$. Omitting that drift would give $P_t=(DP{)}_{xx}$, which is a different Fokker-Planck equation and not the stated Fickian equation.

Noise conventions cannot be mixed silently. If the force covariance is held at the spatially constant value $F_0=2{\zeta }_0{k}_{B}T$ while drag varies, the local overdamped diffusion coefficient is instead $\tilde{D}=F_0/(2{\zeta }^2)$ and its Itô drift is ${\tilde{D}}^{\prime }/2$. Its density equation is $P_t={\partial }_x[\tilde{D}{P}_x+({\tilde{D}}^{\prime }/2)P]$, not the isothermal Fickian equation. This still gives the same displayed first-order mean drift at $x=0$, but has different spatial diffusion and equilibrium properties. The stated Fickian interpretation therefore requires local thermal-noise matching.

Physically, thermally driven excursions travel farther on the lower-drag side. Their displacement is correlated with velocity, so averaging the position-dependent friction leaves a net drift toward greater diffusivity, without an applied force. This is not a violation of equilibrium: in the isothermal Fickian model the drift and spatially varying noise cancel in the probability flux for a uniform density. It is a local diffusion-induced shift of the centroid of an initially localized particle distribution.