Researcher profile

Theodore W. Burkhardt

Theodore W. Burkhardt contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

cond-mat.stat-mech cond-mat.soft math-ph math.MP

Trust snapshot

Quick read

Trust 21 - Emerging

10works

0followers

4topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

preprint2016arXiv

First Passage of a Randomly Accelerated Particle

In the random acceleration process, a point particle is accelerated according to $\ddot{x}=η(t)$, where the right hand side represents Gaussian white noise with zero mean. We begin with the case of a particle with initial position $x_0$ and initial velocity $v_0$ and review the statistics of its first arrival at the origin and its first return to the origin. Multiple returns to the origin, motion with a constant force in addition to a random force, and persistence properties for several boundary conditions at the origin are also considered. Next we review first-exit properties of a randomly accelerated particle from the finite interval $0<x<1$. Then the close connection between the extreme value statistics of a randomly accelerated particle and its first-passage properties is discussed. Finally some applications where first-passage statistics of the random acceleration process play a role are considered.

preprint2016arXiv

Occupation time statistics of the random acceleration model

The random acceleration model is one of the simplest non-Markovian stochastic systems and has been widely studied in connection with applications in physics and mathematics. However, the occupation time and related properties are non-trivial and not yet completely understood. In this paper we consider the occupation time $T_+$ of the one-dimensional random acceleration model on the positive half-axis. We calculate the first two moments of $T_+$ analytically and also study the statistics of $T_+$ with Monte Carlo simulations. One goal of our work was to ascertain whether the occupation time $T_+$ and the time $T_m$ at which the maximum of the process is attained are statistically equivalent. For regular Brownian motion the distributions of $T_+$ and $T_m$ coincide and are given by Lévy's arcsine law. We show that for randomly accelerated motion the distributions of $T_+$ and $T_m$ are quite similar but not identical. This conclusion follows from the exact results for the moments of the distributions and is also consistent with our Monte Carlo simulations.

preprint2011arXiv

Harmonically confined, semiflexible polymer in a channel: response to a stretching force and spatial distribution of the endpoints

We consider an inextensible, semiflexible polymer or worm-like chain which is confined in the transverse direction by a parabolic potential and subject to a longitudinal force at the ends, so that the polymer is stretched out and backfolding is negligible. Simple analytic expressions for the partition function, valid in this regime, are obtained for chains of arbitrary length with a variety of boundary conditions at the ends. The spatial distribution of the end points or radial distribution function is also analyzed.

preprint2010arXiv

Fluctuations of a long, semiflexible polymer in a narrow channel

We consider an inextensible, semiflexible polymer or worm-like chain, with persistence length $P$ and contour length $L$, fluctuating in a cylindrical channel of diameter $D$. In the regime $D\ll P\ll L$, corresponding to a long, tightly confined polymer, the average length of the channel $<R_\parallel>$ occupied by the polymer and the mean square deviation from the average vary as $<R_\parallel>=[1-α_\circ(D/P)^{2/3}]L$ and $<ΔR_\parallel^{\thinspace 2}\thinspace>=β_\circ(D^2/P)L$, respectively, where $α_\circ$ and $β_\circ$ are dimensionless amplitudes. In earlier work we determined $α_\circ$ and the analogous amplitude $α_\Box$ for a channel with a rectangular cross section from simulations of very long chains. In this paper we estimate $β_\circ$ and $β_\Box$ from the simulations. The estimates are compared with exact analytical results for a semiflexible polymer confined in the transverse direction by a parabolic potential instead of a channel and with a recent experiment. For the parabolic confining potential we also obtain a simple analytic result for the distribution of $R_\parallel$ or radial distribution function, which is asymptotically exact for large $L$ and has the skewed shape seen experimentally.

preprint2008arXiv

First-passage and extreme-value statistics of a particle subject to a constant force plus a random force

We consider a particle which moves on the x axis and is subject to a constant force, such as gravity, plus a random force in the form of Gaussian white noise. We analyze the statistics of first arrival at point $x_1$ of a particle which starts at $x_0$ with velocity $v_0$. The probability that the particle has not yet arrived at $x_1$ after a time $t$, the mean time of first arrival, and the velocity distribution at first arrival are all considered. We also study the statistics of the first return of the particle to its starting point. Finally, we point out that the extreme-value statistics of the particle and the first-passage statistics are closely related, and we derive the distribution of the maximum displacement $m={\rm max}_t[x(t)]$.

preprint2007arXiv

Free energy and extension of a semiflexible polymer in cylindrical confining geometries

We consider a long, semiflexible polymer, with persistence length $P$ and contour length $L$, fluctuating in a narrow cylindrical channel of diameter $D$. In the regime $D\ll P\ll L$ the free energy of confinement $ΔF$ and the length of the channel $R_\parallel$ occupied by the polymer are given by Odijk's relations $ΔF/R_\parallel=A_\circ k_BTP^{-1/3}D^{-2/3}$ and $R_\parallel=L[1-α_\circ(D/P)^{2/3}]$, where $A_\circ$ and $α_\circ$ are dimensionless amplitudes. Using a simulation algorithm inspired by PERM (Pruned Enriched Rosenbluth Method), which yields results for very long polymers, we determine $A_\circ$ and $α_\circ$ and the analogous amplitudes for a channel with a rectangular cross section. For a semiflexible polymer confined to the surface of a cylinder, the corresponding amplitudes are derived with an exact analytic approach. The results are relevant for interpreting experiments on biopolymers in microchannels or microfluidic devices.

preprint2006arXiv

Equilibrium statistics of an inelastically bouncing ball, subject to gravity and a random force

We consider a particle moving on the half line $x>0$ and subject to a constant force in the $-x$ direction plus a delta-correlated random force. At $x=0$ the particle is reflected inelastically. The velocities just after and before reflection satisfy $v_f=-rv_i$, where $r$ is the coefficient of restitution. This simple model is of interest in connection with studies of driven granular matter in a gravitational field. With an exact analytical approach and simulations we study the steady state distribution function $P(x,v)$.

preprint2005arXiv

Randomly accelerated particle in a box: mean absorption time for partially absorbing and inelastic boundaries

Consider a particle which is randomly accelerated by Gaussian white noise on the line segment $0<x<1$ and is absorbed as soon as it reaches $x=0$ or $x=1$. The mean absorption time $T(x,v)$, where $x$ and $v$ denote the initial position and velocity, was calculated exactly by Masoliver and Porrà in 1995. We consider a more general boundary condition. On arriving at either boundary, the particle is absorbed with probability $1-p$ and reflected with probability $p$. The reflections are inelastic, with coefficient of restitution $r$. With exact analytical and numerical methods and simulations, we study the mean absorption time as a function of $p$ and $r$.

preprint2004arXiv

Equilibrium of a confined, randomly-accelerated, inelastic particle: Is there inelastic collapse?

We consider the one-dimensional motion of a particle randomly accelerated by Gaussian white noise on the line segment 0<x<1. The reflections of the particle from the boundaries at x=0 and 1 are inelastic, with coefficient of restitution r. We have solved the Fokker-Planck equation satisfied by the equilibrium distribution function P(x,v) with a combination of exact analytical and numerical methods. Throughout the interval 0<r<1, P(x,v) remains extended, as opposed to collapsed. The particle is not localized at the boundary. However, for r<0.163 the equilibrium boundary collision rate is infinite, as predicted by Cornell et al., and all moments of the velocity just after reflection from the boundary vanish.

preprint1997arXiv

Surface Critical Behavior of Binary Alloys and Antiferromagnets: Dependence of the Universality Class on Surface Orientation

The surface critical behavior of semi-infinite (a) binary alloys with a continuous order-disorder transition and (b) Ising antiferromagnets in the presence of a magnetic field is considered. In contrast to ferromagnets, the surface universality class of these systems depends on the orientation of the surface with respect to the crystal axes. There is ordinary and extraordinary surface critical behavior for orientations that preserve and break the two-sublattice symmetry, respectively. This is confirmed by transfer-matrix calculations for the two-dimensional antiferromagnet and other evidence.