Source author record

Mykyta V. Chubynsky

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

ResearcherUnclaimed source record

cond-mat.stat-mech physics.comp-ph cond-mat.soft

Catalog footprint

What is connected

3works

3topics

2close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2014arXiv

"Diffusing diffusivity": A model for anomalous and "anomalous yet Brownian" diffusion

Wang et al. [PNAS 106 (2009) 15160] have found that in several systems the linear time dependence of the mean-square displacement (MSD) of diffusing colloidal particles, typical of normal diffusion, is accompanied by a non-Gaussian displacement distribution (DisD), with roughly exponential tails at short times, a situation they termed "anomalous yet Brownian" diffusion. The diversity of systems in which this is observed calls for a generic model. We present such a model where there is "diffusivity memory" but no "direction memory" in the particle trajectory, and we show that it leads to both a linear MSD and a non-Gaussian DisD at short times. In our model, the diffusivity is undergoing a (perhaps biased) random walk, hence the expression "diffusing diffusivity". The DisD is predicted to be exactly exponential at short times if the distribution of diffusivities is itself exponential, but an exponential remains a good fit to the DisD for a variety of diffusivity distributions. Moreover, our generic model can be modified to produce subdiffusion.

preprint2012arXiv

Monte Carlo Approaches for Simulating a Particle at a Diffusivity Interface and the "Ito--Stratonovich Dilemma"

The possibility of different interpretations of the stochastic term (or calculi) in the overdamped Langevin equation for the motion of a particle in an inhomogeneous medium is often referred to as the "Ito--Stratonovich dilemma," although there is, in fact, a continuum of choices. We introduce two Monte Carlo (MC) simulation approaches for studying such systems, with both approaches giving the choice between different calculi (in particular, Ito, Stratonovich, and "isothermal"). To demonstrate these approaches, we study the diffusion on a 1D interval of a particle released at an interface (in the middle of the system) between two media where this particle has different diffusivities (e.g., two fluids with different viscosities). For reflecting boundary conditions at the ends of the interval, a discontinuity at the interface in the stationary-state particle distribution is found, except for the isothermal case, as expected. We also study the first-passage problem using absorbing boundary conditions. Good agreement is found when comparing the MC approaches to theory as well as Brownian and Langevin dynamics simulations. Also, the results themselves turn out to be interesting. For instance: 1) for some calculi, there are more particles on the low-viscosity side (LVS) at earlier times and then fewer particles at later times; 2) there is no preference to end up on a particular wall for the Ito variant, but a bias towards the wall on the LVS in all other cases; 3) the mean first-passage time to the wall on the LVS grows as the viscosity on the other side is increased, except for the isothermal case where it approaches a constant; 4) when the viscosity ratio is high, the first-passage-time distribution for the wall on the LVS is much broader than for the other wall, with a power-law dependence of the former in a certain time interval whose exponent depends on the calculus.

preprint2010arXiv

Optimizing the accuracy of Lattice Monte Carlo algorithms for simulating diffusion

The behavior of a Lattice Monte Carlo algorithm (if it is designed correctly) must approach that of the continuum system that it is designed to simulate as the time step and the mesh step tend to zero. However, we show for an algorithm for unbiased particle diffusion that if one of these two parameters remains fixed, the accuracy of the algorithm is optimal for a certain finite value of the other parameter. In one dimension, the optimal algorithm with moves to the two nearest neighbor sites reproduces the correct second and fourth moments (and minimizes the error for the higher moments at large times) of the particle distribution and preserves the first two moments of the first-passage time distributions. In two and three dimensions, the same level of accuracy requires simultaneous moves along two axes ("diagonal" moves). Such moves attempting to cross an impenetrable boundary should be projected along the boundary, rather than simply rejected. We also treat the case of absorbing boundaries.