Source author record

David C. Morse

David C. Morse appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

6 published item(s)

preprint2014arXiv

Strong scaling of general-purpose molecular dynamics simulations on GPUs

We describe a highly optimized implementation of MPI domain decomposition in a GPU-enabled, general-purpose molecular dynamics code, HOOMD-blue (Anderson and Glotzer, arXiv:1308.5587). Our approach is inspired by a traditional CPU-based code, LAMMPS (Plimpton, J. Comp. Phys. 117, 1995), but is implemented within a code that was designed for execution on GPUs from the start (Anderson et al., J. Comp. Phys. 227, 2008). The software supports short-ranged pair force and bond force fields and achieves optimal GPU performance using an autotuning algorithm. We are able to demonstrate equivalent or superior scaling on up to 3,375 GPUs in Lennard-Jones and dissipative particle dynamics (DPD) simulations of up to 108 million particles. GPUDirect RDMA capabilities in recent GPU generations provide better performance in full double precision calculations. For a representative polymer physics application, HOOMD-blue 1.0 provides an effective GPU vs. CPU node speed-up of 12.5x.

preprint2014arXiv

Universality of Block Copolymer Melts

Simulations of five different coarse-grained models of symmetric diblock copolymer melts are compared to demonstrate a universal (i.e., model-independent) dependence of the free energy on the invariant degree of polymerization $\overline{N}$, and to study universal properties of the order-disorder transition (ODT). The ODT appears to exhibit two regimes: Systems of very long chains ($\overline{N} \gtrsim 10^{4}$) are well described by the Fredrickson-Helfand theory, which assumes weak segregation near the ODT. Systems of smaller but experimentally relevant values, $\overline{N} \lesssim 10^4$, undergo a transition between strongly segregated disordered and lamellar phases that, though universal, is not adequately described by any existing theory.

preprint2011arXiv

Fluctuations in Symmetric Diblock Copolymers: Testing A Recent Theory

Composition fluctuations in disordered melts of symmetric diblock copolymers are studied by Monte Carlo simulation over a range of chain lengths and interaction strengths. Results are used to test three theories: (1) the random phase approximation (RPA), (2) the Fredrickson-Helfand (FH) theory, which was designed to describe large fluctuations near an order-disorder transition (ODT), and (3) a more recent renormalized one-loop (ROL) theory, which reduces to FH theory near the ODT, but which is found to be accurate over a much wider range of parameters.

preprint2011arXiv

Renormalized One-loop Theory of Correlations in Disordered Diblock Copolymers

A renormalized one-loop theory (ROL) is used to calculate corrections to the random phase approximation (RPA) for the structure factor $\Sc(q)$ in disordered diblock copolymer melts. Predictions are given for the peak intensity $S(q^{\star})$, peak position $q^{\star}$, and single-chain statistics for symmetric and asymmetric copolymers as functions of $χN$, where $χ$ is the Flory-Huggins interaction parameter and $N$ is the degree of polymerization. The ROL and Fredrickson-Helfand (FH) theories are found to yield asymptotically equivalent results for the dependence of the peak intensity $S(q^{\star})$ upon $χN$ for symmetric diblock copolymers in the limit of strong scattering, or large $χN$, but yield qualitatively different predictions for symmetric copolymers far from the ODT and for asymmetric copolymers. The ROL theory predicts a suppression of $S(q^\star)$ and a decrease of $q^{\star}$ for large values of $χN$, relative to the RPA predictions, but an enhancement of $S(q^{\star})$ and an increase in $q^{\star}$ for small $χN$ ($χN < 5$). By separating intra- and inter-molecular contributions to $S^{-1}(q)$, we show that the decrease in $q^{\star}$ near the ODT is caused by the $q$ dependence of the intermolecular direct correlation function, and is unrelated to any change in single-chain statistics, but that the increase in $q^{\star}$ at small values of $χN$ is a result of non-Gaussian single-chain statistics.

preprint2008arXiv

On the chain length dependence of local correlations in polymer melts and a perturbation theory of symmetric polymer blends

The self-consistent field (SCF) theory of dense polymer liquids assumes that short-range correlations are almost independent of how monomers are connected into polymers. Some limits of this idea are explored in the context of a perturbation theory for mixtures of structurally identical polymer species, A and B, in which the AB pair interaction differs slightly from the AA and BB interaction, and the difference is controlled by a parameter alpha Expanding the free energy to O(α) yields an excess free energy of the form alpha $z(N)ϕ_{A}ϕ_{B}$, in both lattice and continuum models, where z(N) is a measure of the number of inter-molecular near neighbors of each monomer in a one-component liquid. This quantity decreases slightly with increasing N because the self-concentration of monomers from the same chain is slightly higher for longer chains, creating a deeper correlation hole for longer chains. We analyze the resulting $N$-dependence, and predict that $z(N) = z^{\infty}[1 + β\bar{N}^{-1/2}]$, where $\bar{N}$ is an invariant degree of polymerization, and $β=(6/π)^{3/2}$. This and other predictions are confirmed by comparison to simulations. We also propose a way to estimate the effective interaction parameter appropriate for comparisons of simulation data to SCF theory and to coarse-grained theories of corrections to SCF theory, which is based on an extrapolation of coefficients in this perturbation theory to the limit $N \to \infty$. We show that a renormalized one-loop theory contains a quantitatively correct description of the $N$-dependence of local structure studied here.

preprint2008arXiv

Renormalized one-loop theory of correlations in polymer blends

The renormalized one-loop theory is a coarse-grained theory of corrections to the self-consistent field theory (SCFT) of polymer liquids, and to the random phase approximation (RPA) theory of composition fluctuations. We present predictions of corrections to the RPA for the structure function $S(k)$ and to the random walk model of single-chain statics in binary homopolymer blends. We consider an apparent interaction parameter $χ_{a}$ that is defined by applying the RPA to the small $k$ limit of $S(k)$. The predicted deviation of $χ_{a}$ from its long chain limit is proportional to $N^{-1/2}$, where $N$ is chain length. This deviation is positive (i.e., destabilizing) for weakly non-ideal mixtures, with $χ_{a} N \alt 1$, but negative (stabilizing) near the critical point. The positive correction to $χ_{a}$ for low values of $χ_{a} N$ is a result of the fact that monomers in mixtures of shorter chains are slightly less strongly shielded from intermolecular contacts. The depression in $χ_{a}$ near the critical point is a result of long-wavelength composition fluctuations. The one-loop theory predicts a shift in the critical temperature of ${\cal O}(N^{-1/2})$, which is much greater than the predicted ${\cal O}(N^{-1})$ width of the Ginzburg region. Chain dimensions deviate slightly from those of a random walk even in a one-component melt, and contract slightly with increasing $χ_{e}$. Predictions for $S(k)$ and single-chain properties are compared to published lattice Monte Carlo simulations.