Takashi Uneyama

Density functional simulation of spontaneous formation of vesicle in block copolymer solutions

The author carries out numerical simulations of vesicle formation based on the density functional theory for block copolymer solutions. It is shown by solving the time evolution equations for concentrations that a polymer vesicle is spontaneously formed from the homogeneous state. The vesicle formation mechanism obtained by this simulation agrees with the results of other simulations based on the particle models as well as experiments. By changing parameters such as the volume fraction of polymers or the Flory-Huggins interaction parameter between the hydrophobic subchains and solvents, the spherical micelles, cylindrical micelles, or bilayer structures can also be obtained. The author also shows that the morphological transition dynamics of the micellar structures can be reproduced by controlling the Flory-Huggins interaction parameter.

Multi-chain slip-spring model for entangled polymer dynamics

It has been established that entangled polymer dynamics can be reasonably described by single chain models such as tube and slip-link models. Although the entanglement effect is a result of hard-core interaction between chains, linkage between the single chain models and the real multi-chain system has not been established yet. In this study, we propose a multi-chain slip-spring model where bead-spring chains are dispersed in space and connected by slip-springs inspired by the single chain slip-spring model [A. E. Likhtman, Macromolecules 38, 6128 (2005)]. In this model the entanglement effect is replaced by the slip-springs, not by the hard-core interaction between beads so that this model is located in the niche between conventional multi-chain simulations and single chain models. The set of state variables are the position of beads and the connectivity (indices) of the slip-springs between beads. The dynamics of the system is described by the time evolution equation and stochastic transition dynamics for these variables. We propose a simple model which is based on the well-defined total free-energy and detailed balance condition. The free energy in our model contains a repulsive interaction between beads, which compensate the attractive interaction artificially generated by the slip-springs. The explicit expression of linear relaxation modulus is also derived by the linear response theory. We also propose a possible numerical scheme to perform simulations. Simulations reproduced expected bead number dependence in transitional regime between Rouse and entangled dynamics for the chain structure, the central bead diffusion, and the linear relaxation modulus.

Concentration dependence of rheological properties of telechelic associative polymer solutions.

We consider concentration dependence of rheological properties of associative telechelic polymer solutions. Experimental results for model telechelic polymer solutions show rather strong concentration dependence of rheological properties. For solutions with relatively high concentrations, linear viscoelasticity deviates from the single Maxwell behavior. The concentration dependence of characteristic relaxation time and moduli is different in high- and low-concentration cases. These results suggest that there are two different concentration regimes. We expect that densely connected (well percolated) networks are formed in high-concentration solutions, whereas sparsely connected (weakly percolated) networks are formed in low-concentration solutions. We propose single chain type transient network models to explain experimental results. Our models incorporate the spatial correlation effect of micellar cores and average number of elastically active chains per micellar core (the network functionality). Our models can reproduce nonsingle Maxwellian relaxation and nonlinear rheological behavior such as the shear thickening and thinning. They are qualitatively consistent with experimental results. In our models, the linear rheological behavior is mainly attributable to the difference of network structures (functionalities). The nonlinear rheological behavior is attributable to the nonlinear flow rate dependence of the spatial correlation of micellar core positions.

Structure of entangled polymer network from primitive chain network simulations

The primitive chain network (PCN) model successfully employed to simulate the rheology of entangled polymers is here tested versus less coarse-grained (lattice or atomistic) models for what concerns the structure of the network at equilibrium (i.e., in the absence of flow). By network structure, we mean the distributions of some relevant quantities such as subchain length in space or in monomer number. Indeed, lattice and atomistic simulations are obviously more accurate, but are also more difficult to use in nonequilibrium flow situations, especially for long entangled polymers. Conversely, the coarse-grained PCN model that deals more easily with rheology lacks, strictly speaking, a rigorous foundation. It is therefore important to verify whether or not the equilibrium structure of the network predicted by the PCN model is consistent with the results recently obtained by using lattice and atomistic simulations. In this work, we focus on single chain properties of the entangled network. Considering the significant differences in modeling the polymer molecules, the results here obtained appear encouraging, thus providing a more solid foundation to Brownian simulations based on the PCN model. Comparison with the existing theories also proves favorable.

Accelerating the singular value decomposition of rectangular matrices with the CSK600 and the integrable SVD

We propose an approach to speed up the singular value decomposition (SVD) of very large rectangular matrices using the CSX600 floating point coprocessor. The CSX600-based acceleration board we use offers 50GFLOPS of sustained performance, which is many times greater than that provided by standard microprocessors. However, this performance can be achieved only when a vendor-supplied matrix-matrix multiplication routine is used and the matrix size is sufficiently large. In this paper, we optimize two of the major components of rectangular SVD, namely, QR decomposition of the input matrix and back-transformation of the left singular vectors by matrix Q, so that large-size matrix multiplications can be used efficiently. In addition, we use the Integrable SVD algorithm to compute the SVD of an intermediate bidiagonal matrix. This helps to further speed up the computation and reduce the memory requirements. As a result, we achieved up to 3.5 times speedup over the Intel Math Kernel Library running on an 3.2GHz Xeon processor when computing the SVD of a 100,000 × 4000 matrix.

Fluctuation analysis of time-averaged mean-square displacement for the Langevin equation with time-dependent and fluctuating diffusivity.

The mean-square displacement (MSD) is widely utilized to study the dynamical properties of stochastic processes. The time-averaged MSD (TAMSD) provides some information on the dynamics which cannot be extracted from the ensemble-averaged MSD. In particular, the relative standard deviation (RSD) of the TAMSD can be utilized to study the long-time relaxation behavior. In this work, we consider a class of Langevin equations which are multiplicatively coupled to time-dependent and fluctuating diffusivities. Various interesting dynamics models such as entangled polymers and supercooled liquids can be interpreted as the Langevin equations with time-dependent and fluctuating diffusivities. We derive a general formula for the RSD of the TAMSD for the Langevin equation with the time-dependent and fluctuating diffusivity. We show that the RSD can be expressed in terms of the correlation function of the diffusivity. The RSD exhibits the crossover at the long time region. The crossover time is related to a weighted average relaxation time for the diffusivity. Thus the crossover time gives some information on the relaxation time of fluctuating diffusivity which cannot be extracted from the ensemble-averaged MSD. We discuss the universality and possible applications of the formula via some simple examples.

Anisotropic mobility model for polymers under shear and its linear response functions.

We propose a simple dynamic model of polymers under shear with an anisotropic mobility tensor. We calculate the shear viscosity, the rheo-dielectric response function, and the parallel relaxation modulus under shear flow deduced from our model. We utilize recently developed linear response theories for nonequilibrium systems to calculate linear response functions. Our results are qualitatively consistent with experimental results. We show that our anisotropic mobility model can reproduce essential dynamical nature of polymers under shear qualitatively. We compare our model with other models or theories such as the convective constraint release model or nonequilibrium linear response theories.

Detailed balance condition and effective free energy in the primitive chain network model

We consider statistical mechanical properties of the primitive chain network (PCN) model for entangled polymers from its dynamic equations. We show that the dynamic equation for the segment number of the PCN model does not reduce to the standard Langevin equation which satisfies the detailed balance condition. We propose heuristic modifications for the PCN dynamic equation for the segment number, to make it reduce to the standard Langevin equation. We analyse some equilibrium statistical properties of the modified PCN model, by using the effective free energy obtained from the modified PCN dynamic equations. The PCN effective free energy can be interpreted as the sum of the ideal Gaussian chain free energy and the repulsive interaction energy between slip-links. By using the single chain approximation, we calculate several distribution functions of the PCN model. The obtained distribution functions are qualitatively different from ones for the simple slip-link model without any direct interactions between slip-links.

Primitive chain network simulations for entangled DNA solutions.

Molecular theories for polymer rheology are based on conformational dynamics of the polymeric chain. Hence, measurements directly related to molecular conformations appear more appealing than indirect ones obtained from rheology. In this study, primitive chain network simulations are compared to experimental data of entangled DNA solutions [Teixeira et al., Macromolecules 40, 2461 (2007)]. In addition to rheological comparisons of both linear and nonlinear viscoelasticities, a molecular extension measure obtained by Teixeira et al. through fluorescent microscopy is compared to simulations, in terms of both averages and distributions. The influence of flow on conformational distributions has never been simulated for the case of entangled polymers, and how DNA molecular individualism extends to the entangled regime is not known. The linear viscoelastic response and the viscosity growth curve in the nonlinear regime are found in good agreement with data for various DNA concentrations. Conversely, the molecular extension measure shows significant departures, even under equilibrium conditions. The reason for such discrepancies remains unknown.

Equilibrium statistics of weakly slip‐linked Gaussian polymer chains

We calculate the free energy and the pressure of a weakly slip-linked Gaussian polymer chains. We show that the equilibrium statistics of a slip-linked system is different from one of the corresponding ideal chain system without any constraints by slip-links. It is shown that the pressure of a slip-linked system decreases compared with the ideal system, which implies that slip-linked chains spontaneously form aggregated cluster like compact structures. These are qualitatively consistent with previous theoretical analyses or multi chain simulations. We also show that repulsive potentials between chains, which have been phenomenologically utilized in simulations, can cancel the artificial pressure decrease.