期刊
SIAM JOURNAL ON APPLIED MATHEMATICS
卷 69, 期 5, 页码 1277-1308出版社
SIAM PUBLICATIONS
DOI: 10.1137/070695186
关键词
generalized Langevin equation; maximum likelihood; Kalman filter; microrheology; anomalous diffusion; time series analysis
资金
- NIH [R01 HL077546-01A2, R01-GM078994]
- ARO [47089-MS-SR]
- NSF [DMS-0604891, DMS-0502266, DMS-0403040, DMS-0308687]
- NATIONAL CANCER INSTITUTE [R33CA155618] Funding Source: NIH RePORTER
- NATIONAL HEART, LUNG, AND BLOOD INSTITUTE [R01HL077546] Funding Source: NIH RePORTER
- NATIONAL INSTITUTE OF BIOMEDICAL IMAGING AND BIOENGINEERING [P41EB002025] Funding Source: NIH RePORTER
- NATIONAL INSTITUTE OF ENVIRONMENTAL HEALTH SCIENCES [RC1ES018686] Funding Source: NIH RePORTER
- NATIONAL INSTITUTE OF GENERAL MEDICAL SCIENCES [R01GM078994] Funding Source: NIH RePORTER
Passive microrheology [T. G. Mason and D. A. Weitz, Phys. Rev. Lett., 74 (1995), pp. 1250-1253] utilizes measurements of noisy, entropic fluctuations (i.e., diffusive properties) of micron-scale spheres in soft matter to infer bulk frequency-dependent loss and storage moduli. Here, we are concerned exclusively with diffusion of Brownian particles in viscoelastic media, for which the Mason-Weitz theoretical-experimental protocol is ideal and the more challenging inference of bulk viscoelastic moduli is decoupled. The diffusive theory begins with a generalized Langevin equation (GLE) with a memory drag law specified by a kernel. We start with a discrete formulation of the GLE as an autoregressive stochastic process governing microbead paths measured by particle tracking. For the inverse problem (recovery of the memory kernel from experimental data) we apply time series analysis (maximum likelihood estimators via the Kalman filter) directly to bead position data, an alternative to formulas based on mean-squared-displacement statistics in frequency space. For direct modeling, we present statistically exact GLE algorithms for individual particle paths as well as statistical correlations for displacement and velocity. Our time-domain methods rest upon a generalization of well-known results for a single-mode exponential kernel to an arbitrary M-mode exponential series, for which the GLE is transformed to a vector Ornstein-Uhlenbeck process.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据