A Real-Time Atmospheric Dispersion Modeling System

This paper describes a new 3-D multi-scale atmospheric dispersion modeling system and its on-going evaluation. This system is being developed for both real-time operational applications and detailed assessments of events involving atmospheric releases of hazardous material. It is part of a new, modernized Department of Energy (DOE) National Atmospheric Release Advisory Center (NARAC) emergency response computer system at Lawrence Livermore National Laboratory. This system contains coupled meteorological data assimilation and dispersion models, initial versions of which were described by Sugiyama and Chan (1998) and Leone et al. (1997). Section 2 describes the current versions of these models, emphasizing new features. This modeling system supports cases involving both simple and complex terrain, and multiple space and time scales from the microscale to mesoscale. Therefore, several levels of verification and evaluation are required. The meteorological data assimilation and interpolation algorithms have been previously evaluated by comparison to observational data (Sugiyama and Chan, 1998). The non-divergence adjustment algorithm was tested against potential flow solutions and wind tunnel data (Chan and Sugiyama, 1997). Initial dispersion model results for a field experiment case study were shown by Leone et al. (1997). A study in which an early, prototype version of the new modeling system was evaluated and compared to the current NARAC operational models showed that the new system provides improved results (Foster et al., 1999). In Section 3, we show example results from the current versions of the models, including verification using analytic solutions to the advection-diffusion equation as well as on-going evaluation using microscale and mesoscale dispersion field experiments

State-dependent diffusion: thermodynamic consistency and its path integral formulation

The friction coefficient of a particle can depend on its position as it does when the particle is near a wall. We formulate the dynamics of particles with such state-dependent friction coefficients in terms of a general Langevin equation with multiplicative noise, whose evaluation requires the introduction of specific rules. Two common conventions, the Ito and the Stratonovich, provide alternative rules for evaluation of the noise, but other conventions are possible. We show the requirement that a particle's distribution function approach the Boltzmann distribution at long times dictates that a drift term must be added to the Langevin equation. This drift term is proportional to the derivative of the diffusion coefficient times a factor that depends on the convention used to define the multiplicative noise. We explore the consequences of this result in a number examples with spatially varying diffusion coefficients. We also derive path integral representations for arbitrary interpretation of the noise, and use it in a perturbative study of correlations in a simple system.Comment: 18 pages, 8 figures, Accepted to PR

Viscosity Dependence of the Folding Rates of Proteins

The viscosity dependence of the folding rates for four sequences (the native state of three sequences is a beta-sheet, while the fourth forms an alpha-helix) is calculated for off-lattice models of proteins. Assuming that the dynamics is given by the Langevin equation we show that the folding rates increase linearly at low viscosities \eta, decrease as 1/\eta at large \eta and have a maximum at intermediate values. The Kramers theory of barrier crossing provides a quantitative fit of the numerical results. By mapping the simulation results to real proteins we estimate that for optimized sequences the time scale for forming a four turn \alpha-helix topology is about 500 nanoseconds, whereas the time scale for forming a beta-sheet topology is about 10 microseconds.Comment: 14 pages, Latex, 3 figures. One figure is also available at http://www.glue.umd.edu/~klimov/seq_I_H.html, to be published in Physical Review Letter

Like-charge attraction through hydrodynamic interaction

We demonstrate that the attractive interaction measured between like-charged colloidal spheres near a wall can be accounted for by a nonequilibrium hydrodynamic effect. We present both analytical results and Brownian dynamics simulations which quantitatively capture the one-wall experiments of Larsen and Grier (Nature 385, p. 230, 1997).Comment: 10 pages, 4 figure

Hydrodynamic Coupling of Two Brownian Spheres to a Planar Surface

We describe direct imaging measurements of the collective and relative diffusion of two colloidal spheres near a flat plate. The bounding surface modifies the spheres' dynamics, even at separations of tens of radii. This behavior is captured by a stokeslet analysis of fluid flow driven by the spheres' and wall's no-slip boundary conditions. In particular, this analysis reveals surprising asymmetry in the normal modes for pair diffusion near a flat surface.Comment: 4 pages, 4 figure

Influence of Hydrodynamic Interactions on Mechanical Unfolding of Proteins

We incorporate hydrodynamic interactions in a structure-based model of ubiquitin and demonstrate that the hydrodynamic coupling may reduce the peak force when stretching the protein at constant speed, especially at larger speeds. Hydrodynamic interactions are also shown to facilitate unfolding at constant force and inhibit stretching by fluid flows.Comment: to be published in Journal of Physics: Condensed Matte

A unified approach for the solution of the Fokker-Planck equation

This paper explores the use of a discrete singular convolution algorithm as a unified approach for numerical integration of the Fokker-Planck equation. The unified features of the discrete singular convolution algorithm are discussed. It is demonstrated that different implementations of the present algorithm, such as global, local, Galerkin, collocation, and finite difference, can be deduced from a single starting point. Three benchmark stochastic systems, the repulsive Wong process, the Black-Scholes equation and a genuine nonlinear model, are employed to illustrate the robustness and to test accuracy of the present approach for the solution of the Fokker-Planck equation via a time-dependent method. An additional example, the incompressible Euler equation, is used to further validate the present approach for more difficult problems. Numerical results indicate that the present unified approach is robust and accurate for solving the Fokker-Planck equation.Comment: 19 page

Stochastic Eulerian Lagrangian Methods for Fluid-Structure Interactions with Thermal Fluctuations

We present approaches for the study of fluid-structure interactions subject to thermal fluctuations. A mixed mechanical description is utilized combining Eulerian and Lagrangian reference frames. We establish general conditions for operators coupling these descriptions. Stochastic driving fields for the formalism are derived using principles from statistical mechanics. The stochastic differential equations of the formalism are found to exhibit significant stiffness in some physical regimes. To cope with this issue, we derive reduced stochastic differential equations for several physical regimes. We also present stochastic numerical methods for each regime to approximate the fluid-structure dynamics and to generate efficiently the required stochastic driving fields. To validate the methodology in each regime, we perform analysis of the invariant probability distribution of the stochastic dynamics of the fluid-structure formalism. We compare this analysis with results from statistical mechanics. To further demonstrate the applicability of the methodology, we perform computational studies for spherical particles having translational and rotational degrees of freedom. We compare these studies with results from fluid mechanics. The presented approach provides for fluid-structure systems a set of rather general computational methods for treating consistently structure mechanics, hydrodynamic coupling, and thermal fluctuations.Comment: 24 pages, 3 figure

Chaperone-assisted translocation of a polymer through a nanopore

Using Langevin dynamics simulations, we investigate the dynamics of chaperone-assisted translocation of a flexible polymer through a nanopore. We find that increasing the binding energy $\epsilon$ between the chaperone and the chain and the chaperone concentration $N_c$ can greatly improve the translocation probability. Particularly, with increasing the chaperone concentration a maximum translocation probability is observed for weak binding. For a fixed chaperone concentration, the histogram of translocation time $\tau$ has a transition from long-tailed distribution to Gaussian distribution with increasing $\epsilon$. $\tau$ rapidly decreases and then almost saturates with increasing binding energy for short chain, however, it has a minimum for longer chains at lower chaperone concentration. We also show that $\tau$ has a minimum as a function of the chaperone concentration. For different $\epsilon$, a nonuniversal dependence of $\tau$ on the chain length $N$ is also observed. These results can be interpreted by characteristic entropic effects for flexible polymers induced by either crowding effect from high chaperone concentration or the intersegmental binding for the high binding energy.Comment: 10 pages, to appear in J. Am. Chem. So

Effect of hydrodynamic interactions on the distribution of adhering Brownian particles

Brownian dynamics simulations were used to study the adhesion of hard spheres on a solid surface by taking the hydrodynamic interactions into account. Special attention was paid to analyze the configuration of the assembly of adsorbed particles. These results were compared to configurations generated by the extensively studied random sequential adsorption (RSA) model. In our case the adsorption probability for a particle is almost uniform over the entire available surfae. This surprising result shows that RSA provides a good approximation to generate adsorbed particle configurations