\section{Introduction}

 We address the time and spatial evolution of the ionization produced
 along a charged particle's track in a fluid (gas or liquid) subjected
 to a non-uniform, time-independent electric field. The prototypical
 calculation is for a localized ion density produced at some point
 ${\bf R}_{0}$ at time $t = 0$, described by a delta function $
 \delta({\bf r} - {\bf R}_{0})$. Of interest is the diffusion in both
 the direction of the electric field and transverse to it, which
 governs the distributions in arrival time and transverse position at
 some detecting surface. We focus on a volume that is bounded by
 cylindrical electrodes that produce a radial electric field, but
 begin with more general considerations.  Spherical geometry is
 considered in an appendix.

 The motion of electrons and atomic or molecular ions in
 electromagnetic fields has a huge experimental and theoretical
 literature. See, for example, the texts and cited literature of
 Huxley and Crompton\cite{Huxley} and Mason and McDaniel\cite{Mason}.
 An extensive literature about diffusion of charged particles in
 plasmas also exists, but conditions there are far from our concern of
 relatively few ions moving in a neutral medium.

















 
 
 We make the simplifying assumptions that the ion mobility $\mu $ and
 the longitudinal and transverse diffusion constants, $D_{L}$ and
 $D_{T}$, are independent of electric field strength $ E $ or rather
 of $ E/N $, where $ N $ is the particle number density of the
 fluid. The data compilations of transport properties of gaseous ions
 \cite{Ellis_1, Ellis_2, Ellis_3} show that the mobilities of ions in
 noble gases are field-independent (to $ \pm 10 \ \% $) typically for
 $E/N < 60-100 \ {\rm Td}$.\cite{notation} The longitudinal
 diffusion constants $ D_{L} $ for ions in noble gases are similarly
 roughly field-independent typically for $ E/N < 10-2 0 \ {\rm Td}
 $. For this range of $ E/N $, the ratio of the diffusion constants
 is $ D_{L}/D_{T} \approx 1$; only at higher values of $ E/N $ is $
 D_{L}/D_{T} $ greater than unity (and both are field-dependent). See, for example,
 figure 5-2-4 of Ref.~\cite{Mason}.
  
 Within our assumptions, we hypothesize a regime in which the
microscopic (atomic) processes that establish the local drift
velocity and diffusion parameters occur on length and time scales
small compared to the changes in magnitude or direction of the
electric field or the average drift position.  Such a regime permits
the use of the simple diffusion equation rather than the full
Boltzmann equation.

 For electrons in gases, field-independent mobilities occur typically
only for $ E/ N < 10^{-2} \ {\rm Td}$,with various field
dependencies at larger $ E/N $.\cite{Huxley} Diffusion constants for
electrons often exhibit rapidly varying field-strength dependence
because of resonant behaviour in the momentum transfer cross sections. Our approximations seem to have limited applicability for
electrons.
 

\section{Preliminaries and initial Ansatz}

Ions in a fluid and subject to an electric field can be described by a
charge density $\psi({\bf r},t)$ satisfying the diffusion equation

\begin{equation}
\frac {d\psi}{dt}=\frac{\partial\psi}{\partial t}+{\bf
  v}\cdot{\bf\nabla}\psi=D\nabla^2\psi.\label{eq:one}
\end{equation}

We take the velocity ${\bf v}$ to be given by the mobility $\mu$ times
the electric field ${\bf E}$.  We imagine that from the prescribed
electric field and mobility we can calculate the classical trajectory
${\bf r}_0(t)$ of a point charge located at ${\bf r}={\bf R}_0$ at time
$t=0$.  We indicate the velocity along the path by ${\bf v}_0(t)$ and,
in particular,
\begin{equation}
{\bf v}_0(t)=\mu {\bf E} ({{\bf r}}_0(t)).
\end{equation}
The trajectory is given implicitly as
\begin{equation}
{\bf r}_0(t)={\bf R}_0 +\int_0^tdt' \mu {\bf E}(r_0(t')).\label{eq:three}
\end{equation}
We are seeking a solution proportional to $\delta({{\bf r}} -{\bf R}_0)$
at $t=0$.  We begin with an Ansatz for $\psi$:

\begin{equation}
\psi=\frac {[\det(I+g)]^{1/2}}{(4\pi Dt)^{3/2}}\exp\left\{-[{\bf r}-{\bf r}_0(t)]_i[\delta_{ij}+g(t)_{ij}][{\bf r}-{\bf r}_0(t)]_j/(4Dt)\right\}.\label{eq:ansatz}
\end{equation}
The factor $[\det(I+g)]^{1/2}$ ensures that $\psi$ is normalized properly:

\begin{equation}
\int d^3{\bf r} \, \psi =1,
\end{equation}
which must be the case if we solve the initial equation, which derives from the continuity equation.  In the absence of the electric field, the familiar solution is given by Eq.(\ref{eq:ansatz}) with $g_{ij}=0$ and ${\bf r}_0(t)={\bf R}_0$.
The basic assumption is that the form is Gaussian with respect to the difference between the point at which the density is evaluated and the location of the classical trajectory at that time, ${\bf r}_0(t)$.  The form is motivated by the observation that, not only in the absence of the electric field, but also for a constant electric field, the solution is given by Eq.(\ref{eq:ansatz}) with $g_{ij}=0$ . 
 We calculate the time and spatial derivatives of $\psi$,  writing $N(t)$ for $[\det(I+g)]^{1/2}$ and 
 $\rho_i$ for $({\bf r}-{\bf r}_0(t))_i$:
\begin{eqnarray}
4Dt\frac{\partial\psi}{\partial t}&=&{{\rho}_i(\delta_{ij}+g_{ij}){\rho}_j}\frac\psi t -6D\psi +{2{\rho}_i (\delta_{ij}+g_{ij})v_{0j}}\psi \nonumber\\
&&\qquad-{{\rho}_i\frac{dg_{ij}}{dt}{\rho}_j}\psi
+4Dt\frac {d\ln N}{dt}\psi;\\
4Dt\nabla_k\psi&=&-{2(\delta_{jk}+g_{jk}){\rho}_j}\psi;\\
4D^2t\nabla_k\nabla_k \psi&=&(\delta_{kj}+g_{kj}){\rho}_j(\delta_{kl}+g_{kl}){\rho}_l\frac\psi t-2(\delta_{kk}+g_{kk})D\psi\nonumber\\
&=&{{\rho}_j(\delta_{jl}+2g_{jl}+g_{kj}g_{kl}){\rho}_l}\frac\psi t
-2(3+g_{kk})D\psi;\\
4Dt\mu{\bf E}\cdot{\bf \nabla}\psi&=&-{2\mu E_k(\delta_{kj}+g_{jk}){\rho}_j}\psi.
\end{eqnarray}
We now expand ${\bf E}({\bf r})$ around ${\bf r}_0(t)$, keeping first order terms,
\begin{equation}
{\bf v}_0(t)-\mu {\bf E}({\bf r})=\mu({\bf E}({\bf r}_0(t))-{\bf E}(r))=-\mu{\rho}_i\partial_i{\bf E}({\bf r}_0(t))+\ldots.
\end{equation}
The omitted terms in the expansion are of order $ (\rho/L)^{n}, \ n \geq 1$, times the first order term, where $L$ is the length scale over which the electric field changes appreciably. If the spread of the diffusion, $ O(\rho_{max})$, is small compared to $ L $, our approximation should be reliable.
We combine these into Eq.(\ref{eq:one})  to find the equations for $g_{ij}$:
\begin{eqnarray}
&&{\rho}_i\frac{(\delta_{ij}+g_{ij})}t{\rho}_j
 -{6D}+4Dt\frac {d\ln N}{dt}\psi
-{\rho}_i\frac{dg_{ij}}{dt}{\rho}_j\nonumber\\
&&\qquad\qquad={\rho}_j\frac{(\delta_{jl}+2g_{jl}+g_{kj}g_{kl})}t{\rho}_l
-2(3+g_{kk})D\nonumber\\
&&\qquad\qquad\qquad+2\mu{\rho}_i\partial_i E_k(\delta_{kj}+g_{jk}){\rho}_j.
\end{eqnarray}
The terms with and without quadratic dependence on ${\rho}$ must cancel separately.  In fact, we find that the equations that derive from terms with quadratic dependence ensure that the terms without that dependence cancel, as well.
The terms without the quadratic dependence give
\begin{equation}
1-\frac{2t}{3}\frac{d\ln N}{dt}=1+\frac 13g_{kk}\label{eq:no-quad}
\end{equation}
The terms with quadratic dependence on ${\rho}$ give
\begin{eqnarray}
&&\frac{(\delta_{ij}+g_{ij})}t
-\frac{dg_{ij}}{dt}=\frac{(\delta_{ji}+2g_{ji}+g_{kj}g_{ki})}t
+2\mu\partial_i E_k(\delta_{kj}+g_{jk})
\end{eqnarray}
 or
 \begin{equation}
 \frac{dg_{ij}}{dt}+\frac{g_{ij}+g_{ik}g_{kj}}t + 2\mu\partial_iE_k(\delta_{kj}+g_{kj})=0\label{eq:b}
 \end{equation}
 \section{Restriction to rectilinear motion and symmetric geometry}
 Only the symmetric part of $g_{ij}$ enters $\psi$.  Consistency requires that the antisymmetric part of Eq.(\ref{eq:b}) vanish:
 
 \begin{equation}
 (\delta_{ij}+g_{ij})\partial_kE_j=(\delta_{kj}+g_{kj})\partial_i E_j
 \end{equation}
 The vanishing of the curl of ${\bf E}$ is enough to take care of the $\delta_{ij}$ portion but we are left with the requirement
 \begin{equation}
 g_{ij}\partial_k E_j=g_{kj}\partial_i E_j\label{requirement}
 \end{equation}
 
 In general, the product of two symmetric matrices is not symmetric.  We need to impose additional constraints.  In particular, we suppose that the motion is rectilinear along the $x$ direction, i.e.
 ${\bf r}(t)=r(t){\bf \hat x}$, ${\bf E}(t)=E_x(t){\bf \hat x}$. Write the electrostatic potential in the neighborhood of the line of motion as an expansion in $y$ and $z$:
 \begin{eqnarray}
 \Phi(x,y,z)&=& \Phi(x,0,0) + y\Phi_y(x) +z\Phi_z(x) \nonumber\\
&&+\frac 12 y^2 \Phi_{yy}(x) +\frac 12z^2\Phi_{zz}(x)\nonumber\\
&& +xy\Phi_{xy}(x) + xz\Phi_{xz}(x)+yz\Phi_{yz}(x)+\ldots\nonumber\\\label{eq:expansion}
 \end{eqnarray}
 
 By hypothesis the electric field along $y=0,\ z=0$ is entirely in the $x$ direction. Then
 $\Phi_y(x)=\Phi_z(x)=\Phi_{xy}(x)=\Phi_{xz}(x)=0$ and Eq.(\ref{eq:expansion}) becomes, to second order inclusive,
 \begin{eqnarray}
 \Phi(x,y,z)&=& \Phi(x,0,0) 
+\frac 12 y^2 \Phi_{yy}(x) +\frac 12z^2\Phi_{zz}(x)
+yz\Phi_{yz}(x)\nonumber\\
 \end{eqnarray}
 
 If the surface charge distributions determining the electric field
 are symmetric under $y\to -y$ and $z\to -z$, the term $yz\Phi_{yz}$
 vanishes along with $\Phi_y,\Phi_z,\Phi_{xy}$ and $\Phi_{xz}$, and the
 only non-zero values of $\partial_iE_j$ along the line $(x,0,0)$
 occur for $i=j$.  Two clear examples are cylindrical symmetry with
 the $z$ direction taken as the axis of the cylinder and spherical
 symmetry.  In the former case all derivatives with respect to $z$
 vanish so $\partial_z E_z=0$ and in the latter, the derivatives with
 respect to $y$ and $z$ are equal so $\partial_zE_z=\partial_y E_y$.
 See figure \ref{fig:drift1}.
 
\begin{figure}\begin{center}
 \includegraphics[width=4in]{figure_1.eps}
 
 \caption[99]{The Ansatz for $\psi$, Eq.(\ref{eq:ansatz}), is applicable
   for rectilinear trajectories if the matrix $\partial_i E_j$ is
   diagonal, for a suitable choice of coordinates.  Taking $x$ along
   the direction of motion, a charge distribution (or potential)
   symmetric under $y \to -y$ and $z\to -z$ will ensure this.  The
   figure shows such a trajectory and charge distribution.  Here only
   the trajectory shown has the required symmetry.  For the
   cylindrical and spherical geometries all radial trajectories have
   the symmetry required by our Ansatz.}\label{fig:drift1}
\end{center}
\end{figure}
   Under any of these circumstances we can write 
 \begin{equation}
 \partial_i E_j=C_i(x)\delta_{ij}
 \end{equation}
 and we can satisfy Eq.(\ref{requirement}) with the ansatz
 \begin{equation}
 g_{ij}=g_i\delta_{ij}
 \end{equation}
 i.e. both $g_{ij}$ and $\partial_i E_j$ are diagonal along the linear trajectory path of Eq.(\ref{eq:three}).  Now Eq.(\ref{eq:b}) simplifies to
 \begin{equation}
  \frac{dg_{i}}{dt}+\frac{g_{i}(1+g_{i})}t + 2\mu\partial_iE_i(1+g_{i})=0\label{eq:c}
 \end{equation}  
 where there is no summation on $i$.
\section {Example of cylindrical geometry}

For a cylindrical geometry, uniform in the $z$ direction and having azimuthal symmetry, with a radial electric field, we have 
\begin{equation}
{\bf E} = \frac Ar {\bf\hat r},\label{eq:field}
\end{equation}
where ${\bf r} = x{\bf\hat x}+y{\bf\hat y}$ is now two-dimensional.  We assume here that the signs of the charge and field are such that the force is radially outward.  

We can anticipate the effect of the inhomogeneous electric field,
Eq.(\ref{eq:field}).  The portion of the cloud that is further out
than the center of the cloud experiences a weaker field and thus moves
more slowly than the outwardly moving center, while the portion of the
cloud closer in to the central axis experiences a greater field and
moves more rapidly than moving center.  Thus the spread of the cloud
in the radial direction should be smaller than in the absence of the
electric field. On the other hand, the diverging electric field will
tend to spread the cloud transversely, leading to a larger dispersion
than in the absence of the electric field.

Take the initial position ${\bf R}_0=(x_0,0,0)$.  The motion is along the $x$ axis where

\begin{equation}
\partial_xE_x=-\frac A{x^2};\quad \partial_y E_y=\frac A{x^2};\quad \partial_x E_y=0=
\partial_yE_x,
\end{equation}
satisfying the symmetry requirements for our Ansatz.  Thus we have the diagonal form $\partial_i E_j=C_i(x)\delta_{ij}$ 
where
\begin{eqnarray}
C_1(x)&=&-\frac A{x^2},\\
C_2(x)&=& \frac A{x^2},\\
C_3(x)&=&0.
\end{eqnarray}


 The trajectory is given by ${\bf r}_0(t)=r_0(t){\hat{\bf x}}$, with
\begin{equation}
r_0(t)^2=x_0^2+2\mu A t,
\end{equation}
so that
\begin{equation}
2\mu \partial_xE_x=-\frac{2\mu A}{x_0^2+2\mu A t}=-2\mu \partial_yE_y.
\end{equation}
  Introduce
\begin{equation}
T=\frac{x_0^2}{2\mu A}
\end{equation}
so that
\begin{equation}
r_0(t)=x_0\sqrt{1+\frac tT}.\label{eq:r_of_t}
\end{equation}
In the characteristic time $T$ the diffusion center would go a
distance $x_0/2$ if the field had the constant value $E=A/x_0$.  When
$t=3T$, the classical particle is twice as far from the axis  of the
cylinder in the $1/r$ field as it was at $t=0$.  From Eq.(\ref{eq:c}):

\begin{eqnarray}
\frac{dg_1}{dt}&=&-\frac{g_1(1+g_1)}t+\frac 1{T+t}(1+g_1);\nonumber\\
\frac{dg_2}{dt}&=&-\frac{g_2(1+g_2)}t-\frac 1{T+t}(1+g_2);\nonumber\\
\frac{dg_3}{dt}&=&-\frac{g_3(1+g_3)}t.\nonumber\\
\end{eqnarray}
Finally, writing $t=Ts$ 
\begin{eqnarray}
\frac{dg_1}{ds}&=&-\frac{g_1(1+g_1)}s+\frac 1{1+s}(1+g_1);\nonumber\\
\frac{dg_2}{ds}&=&-\frac{g_2(1+g_2)}s-\frac 1{1+s}(1+g_2);\nonumber\\
\frac{dg_3}{ds}&=&-\frac{g_3(1+g_3)}s.\nonumber\\\label{eq:nonlinear}
\end{eqnarray}
The equation for $g_3$ has the solution 

\begin{equation}
g_3=\frac 1{Bs-1}
\end{equation}
where $B$ is a constant.  When $s=0$, all the $g_i$ are zero.  It follows that $B=\infty$ and $g_3$ remains zero: Diffusion in
the $z$ direction is unchanged by the presence of the radial electric field.

Eqs.(\ref{eq:nonlinear}) for $g_1$ and $g_2$ are Ricatti equations and  can be solved analytically.
The details of the solutions are given in Appendix A.  The results are
\begin{eqnarray}
g_1&=& \frac s{2+s}\\
g_2&=&-1 +\frac s{(1+s)\ln(1+s)}
\end{eqnarray}
The behavior of the functions $g_1$ and $g_2$ is shown in figure \ref{fig:g1g2}.

\begin{figure}\begin{center}
\includegraphics[width=4in,angle=90]{figure_2.eps}
\caption{The functions appearing in the Gaussian description of the charge density $\psi$ for outward drift, as a function of the dimensionless time variable $s=t/T$ for cylindrical geometry.}\label{fig:g1g2}\end{center}
\end{figure}

The effective parallel and transverse spreads relative to $\sigma_0=\sqrt{2Dt}$ are
\begin{equation}
\frac{\sigma_\parallel}{\sigma_0}=\sqrt{\frac {1}{1+g_1}};\quad 
\frac{\sigma_\perp}{\sigma_0}=\sqrt{\frac{1}{1+g_2}},\label{eq:widths}
\end{equation}
shown in figure \ref{fig:sigmas}.  In the $z$-direction the spread is $\sigma_3=\sigma_0$.
\begin{figure}\begin{center}
\includegraphics[width=4in,angle=90]{figure_3.eps}
\caption[99]{The effective dimensionless widths Eq.(\ref{eq:widths}) for outward drift, in the directions parallel and perpendicular to the radial direction as functions of the dimensionless time variable $s=t/T$ for cylindrical geometry.  If the diffusion constants differ in the longitudinal (parallel to the field) and transverse directions, multiply the ratio $\sigma_\parallel/\sigma_0$ by $\sqrt{D_L/D_T}$ and take $\sigma_0=\sqrt{2D_T t}$. (see Sec. 6).}\label{fig:sigmas}\end{center}
\end{figure}
\clearpage
\section {Inward drift}
We have imagined that the radial drift is outward, towards larger $r$.  If the sign of the electric field is reversed, the derivation goes through much the same way except that
the sign of $T$ is negative.  Since $t=Ts$, positive times then correspond to negative values of $s$.  The center point of the diffusion cloud arrives at the central axis at $s=-1$.  Because the field increases in the direction of motion, we expect the effect of the field to be opposite to what we found above:  the dispersion should be greater in the radial direction and smaller in the transverse direction.  

We show the results obtained for inward drift in Figs. \ref{fig:inward} and \ref{fig:inward-sigmas}.


\begin{figure}\begin{center}
\includegraphics[width=4in,angle=90]{figure_4.eps}
\caption{The functions appearing in the Gaussian description of $\psi$ for inward drift in the cylinder. Negative values of $s$ correspond to positive values of $t$ according to $t=Ts$.  For $s=-1$, ${\bf r}_0(t)=0$.}
\label{fig:inward}\end{center}
\end{figure}

\begin{figure}\begin{center}
\includegraphics[width=4in,angle=90]{figure_5.eps}
\caption[99]{The effective dimensionless widths Eq.(\ref{eq:widths}) in
  the directions parallel and perpendicular to the radial direction
  for inward drift in the cylinder. Negative values of $s$ correspond
  to positive values of $t$ according to $t=Ts$.  If the diffusion
  constants differ in the longitudinal (parallel to the field) and
  transverse directions, multiply the ratio $\sigma_\parallel/\sigma_0$ by $\sqrt{D_L/D_T}$ and take $\sigma_0=\sqrt{2D_T t}$ (see
  Sec. 6).}\label{fig:inward-sigmas}\end{center}
\end{figure}

\section{Different coefficients of diffusion for longitudinal and transverse directions}

It is straightforward to include the effect of having different diffusion coefficients $D_L$ and $D_T$ for the longitudinal and transverse directions.  For the Ansatz we take

\begin{equation}
\psi=\frac {[\Pi_j(1+g_j)]^{1/2}}{(4\pi t)^{3/2}D_TD_L^{1/2}}\exp\left\{-\sum_i\frac{\rho_i^2(1+g_i)}{4D_it}\right\}.\label{eq:ansatztwo}
\end{equation}
with $D_1=D_L, D_2=D_3=D_T$.  Exactly the same differential equations emerge, but setting $\sigma_0=\sqrt{2D_Tt}$ we have
\begin{eqnarray}
\frac{\sigma_\parallel}{\sigma_0}&=&\sqrt{\frac{D_L}{D_T}}\frac 1{\sqrt{1+g_1}}\\
\frac{\sigma_\perp}{\sigma_0}&=&\frac 1{\sqrt{1+g_2}}\\
\frac{\sigma_3}{\sigma_0}&=&1
\end{eqnarray}
\section{Diffusion of a track segment}

To illustrate the possible use of our results, instead of a single
point source arising at ${\bf r} = {\bf R}_{0} $ at $t = 0$, we
consider the ionization along the path of a relativistic charged
particle through a cylindrical drift chamber.  The ionization trail is
assumed linear (in both senses) and is, on the time scale relevant for
the subsequent diffusion, deposited instantaneously. The resulting
charge distribution is thus a superposition of distributions of the
form Eq.(\ref{eq:ansatztwo}) for contributions along the particle's
track. For definiteness we consider a track segment with end points
${\bf R}_{1}, \ {\bf R}_{2} $, as shown in figure
(\ref{fig:line-drift}a), with a typical point on the track labeled by
${\bf R}(\xi)$.  For convenience, define the distance between the end
points as $L = |{\bf R}_{2} - {\bf R}_{1}| $.  Then we have ${\bf
  R}(\xi) = {\bf R}_{1} + (\xi/L) ({\bf R}_{2} - {\bf R}_{1}) $.  In
figure (\ref{fig:line-drift}b), the typical starting point ${\bf
  R}(\xi)$ is singled out.  We define $ R_{\parallel}(\xi) = \hat{{\bf
    z}} \cdot {\bf R}(\xi) $ and ${\bf R}_{\perp}(\xi) = {\bf R}(\xi)
- \hat{{\bf z}} R_{\parallel}(\xi) $.  Then the drift center's
coordinate ${\bf r}_{0}(t) $ can be written according to
Eq.(\ref{eq:r_of_t}) as
\begin{equation}
{\bf r}_{0}(t) \ = \ {\bf R}_{\perp}(\xi) [1 + t/T(\xi) ]^{1/2} \ + \ \hat{{\bf z}} R_{\parallel}(\xi) 
\end{equation}
where from (30), since $x_{0}(\xi) = |{\bf R}_{\perp}(\xi)| $, we have
\begin{equation}
 T(\xi) = \frac{|{\bf R}_{\perp}(\xi)|^{2}}{2 \mu A}
\end{equation}
In the exponent of $\psi$ the Cartesian components of the coordinate difference, $({\bf r} - {\bf r}_{0}(t))$, must be taken with respect to ${\bf R}_{\perp}(\xi) $ as the $x$-direction:
\begin{eqnarray}
\rho_1\equiv({\bf r} - {\bf r}_{0}(t))_{1} &  = & ({\bf r} - {\bf r}_{0}(t)) \cdot \hat{{\bf R}}_{\perp}(\xi) \\
\rho_2\equiv({\bf r} - {\bf r}_{0}(t))_{2} & =  & ({\bf r} - {\bf r}_{0}(t)) \cdot (\hat{{\bf z}} \times \hat{{\bf R}}_{\perp}(\xi))\\
\rho_3\equiv({\bf r} - {\bf r}_{0}(t))_{3} & =  &z-R_\parallel(\xi)
\end{eqnarray} 
\begin{figure}\begin{center}
\includegraphics[width=2.5in]{figure_6a.eps}\includegraphics[width=2.7in]{figure_6b.eps}
\caption{\label{fig:line-drift} Left (a): Track segment $({\bf R}_1,{\bf R}_2)$ with typical starting point ${\bf R}(\xi)$.  Right (b): Details of a typical starting point showing its decomposition ${\bf R}(\xi)={\bf R}_\perp(\xi)+R_\parallel {{\bf\hat z}}$ into a piece along the drift direction and a piece parallel to the axis of the cylinder.}
\end{center}
\end{figure}

\indent The normalized diffused charge density in space and time for the track segment can thus be written as a numerical integral over $\psi$ given by Eq.(\ref{eq:ansatztwo}),
\begin{equation}
\Psi({\bf r},  t| {\bf R}_{1}, {\bf R}_{2}) \ = \ \frac{1}{L}\int_{0}^{L} d\xi \ \psi({\bf r}, {\bf R}(\xi), t)
\end{equation}
The functions $g_{1}(s)$ and $g_{2}(s)$ are to be evaluated at $s = t/T(\xi)$; $g_3(s)=0.$\\

Acknowledgments: This work was supported by the Director, Office
of Science, High Energy Physics, U.S. Department of Energy under
Contract No. DE-AC02-05CH11231. We thank David Nygren for drawing our
attention to this class of ionic diffusion problems.


