Showing papers by "Steven J. Plimpton published in 2010"

1D-to-3D transition of phonon heat conduction in polyethylene using molecular dynamics simulations

Abstract: The thermal conductivity of nanostructures generally decreases with decreasing size because of classical size effects. The axial thermal conductivity of polymer chain lattices, however, can exhibit the opposite trend, because of reduced chain-chain anharmonic scattering. This unique feature gives rise to an interesting one-dimensional-to-three-dimensional transition in phonon transport. We study this transition by calculating the thermal conductivity of polyethylene with molecular dynamics simulations. The results are important for designing inexpensive high thermal-conductivity polymers.

1D-to-3D transition of phonon heat conduction in polyethylene using molecular dynamics simulations

Abstract: The thermal conductivity of nanostructures generally decreases with decreasing size because of classical size effects. The axial thermal conductivity of polymer chain lattices, however, can exhibit the opposite trend, because of reduced chain-chain anharmonic scattering. This unique feature gives rise to an interesting one-dimensional-to-three-dimensional transition in phonon transport. We study this transition by calculating the thermal conductivity of polyethylene with molecular dynamics simulations. The results are important for designing inexpensive high thermal-conductivity polymers.

Software components for parallel multiscale simulation: an example with LAMMPS

Abstract: Multiscale simulation is a promising approach for addressing a variety of real-world engineering problems. Various mathematical approaches have been proposed to link single-scale models of physics into multiscale models. In order to be effective, new multiscale simulation algorithms must be implemented which use partial results provided by single-scale software. This work considers aspects of software design for interfacing to existing single-scale simulation code to perform multiscale simulations on a parallel machine. As a practical example, we extended the large-scale atomistic/molecular massively parallel simulator (LAMMPS) atomistic simulation software to facilitate its efficient use as a component of parallel multiscale-simulation software. This required new library interface functions to LAMMPS that side-stepped its dependence on files for input and output and provided efficient access to LAMMPS’s internal state. As a result, we were able to take advantage of LAMMPS’s single-scale performance without adding any multiscale-specific code to LAMMPS itself. We illustrate the use of LAMMPS as a component in three different modes: as a stand-alone application called by a multiscale code, as a parallel library invoked by a serial multiscale code, and as a parallel library invoked by a parallel multiscale code. We conclude that it is possible to efficiently re-use existing single-scale simulation software as a component in multiscale-simulation software.

Mesoscale hydrodynamics via stochastic rotation dynamics: comparison with Lennard-Jones fluid.

Abstract: Stochastic rotation dynamics (SRD) is a relatively recent technique, closely related to lattice Boltzmann, for capturing hydrodynamic fluid flow at the mesoscale. The SRD method is based on simple constituent fluid particle interactions and dynamics. Here we parametrize the SRD fluid to provide a one to one match in the shear viscosity of a Lennard-Jones fluid and present viscosity measurements for a range of such parameters. We demonstrate how to apply the Muller-Plathe reverse perturbation method for determining the shear viscosity of the SRD fluid and discuss how finite system size and momentum exchange rates effect the measured viscosity. The implementation and performance of SRD in a parallel molecular dynamics code is also described.

Novel statistical ensemble analysis for simulating extrinsic noise-driven response in NF-{\kappa}B signaling network

Abstract: Cellular responses in the single cells are known to be highly heterogeneous and individualistic due to the strong influence by extrinsic and intrinsic noise. Here, we are concerned about how to model the extrinsic noise-induced heterogeneous response in the single cells under the constraints of experimentally obtained population-averaged response, but without much detailed kinetic information. We propose a novel statistical ensemble scheme where extrinsic noise is regarded as fluctuations in the values of kinetic parameters and such fluctuations are modeled by randomly sampling the kinetic rate constants from a uniform distribution. We consider a large number of signaling system replicates, each of which has the same network topology, but a uniquely different set of kinetic rate constants. A protein dynamic response from each replicate should represent the dynamics in a single cell and the statistical ensemble average should be regarded as a population-level response averaged over a population of the cells. We devise an optimization algorithm to find the correct uniform distribution of the network parameters, which produces the correct statistical distribution of the response whose ensemble average and distribution agree well with the population-level experimental data and the experimentally observed heterogeneity. We apply this statistical ensemble analysis to a NF-{\kappa}B signaling system and (1) predict the distributions of the heterogeneous NF-{\kappa}B (either oscillatory or non-oscillatory) dynamic patterns and of the dynamic features (e.g., period), (2) predict that both the distribution and the statistical ensemble average of the NF-{\kappa}B dynamic response depends sensitively on the dosage of stimulant, and lastly (3) demonstrate the sigmoidally shaped dose-response from the statistical ensemble average and the individual replicates.

Porting LAMMPS to GPUs.

Abstract: LAMMPS is a classical molecular dynamics code, and an acronym for Large-scale Atomic/Molecular Massively Parallel Simulator. LAMMPS has potentials for soft materials (biomolecules, polymers) and solid-state materials (metals, semiconductors) and coarse-grained or mesoscopic systems. It can be used to model atoms or, more generically, as a parallel particle simulator at the atomic, meso, or continuum scale. LAMMPS runs on single processors or in parallel using message-passing techniques and a spatial-decomposition of the simulation domain. The code is designed to be easy to modify or extend with new functionality.

Noise-induced oscillatory shuttling of NF-{\kappa}B in a two compartment IKK-NF-{\kappa}B-I{\kappa}B-A20 signaling model

Abstract: NF-{\kappa}B is a pleiotropic protein whose nucleo-cytoplasmic trafficking is tightly regulated by multiple negative feedback loops embedded in the NF-{\kappa}B signaling network and contributes to diverse gene expression profiles important in immune cell differentiation, cell apoptosis, and innate immunity. The intracellular signaling processes and their control mechanisms, however, are susceptible to both extrinsic and intrinsic noise. In this article, we present numerical evidence for a universal dynamic behavior of NF-{\kappa}B, namely oscillatory nucleo-cytoplasmic shuttling, due to the fundamentally stochastic nature of the NF-{\kappa}B signaling network. We simulated the effect of extrinsic noise with a deterministic ODE model, using a statistical ensemble approach, generating many copies of the signaling network with different kinetic rates sampled from a biologically feasible parameter space. We modeled the effect of intrinsic noise by simulating the same networks stochastically using the Gillespie algorithm. The results demonstrate that extrinsic noise diversifies the shuttling patterns of NF-{\kappa}B response, whereas intrinsic noise induces oscillatory behavior in many of the otherwise non-oscillatory patterns. We identify two key model parameters which significantly affect the NF-{\kappa}B dynamic response and deduce a two-dimensional phase-diagram of the NF-{\kappa}B response as a function of these parameters. We conclude that if single-cell experiments are performed, a rich variety of NF-{\kappa}B response will be observed, even if population-level experiments, which average response over large numbers of cells, do not evidence oscillatory behavior.