Monte Carlo simulations of stress relaxation of entanglement-free Fraenkel chains. II. Nonlinear polymer viscoelasticity

The nonlinear viscoelastic behavior of the Fraenkel-chain model is studied with respect to the constitutive equation of the Rouse model. Distinctly different from the results of the Rouse model, the Fraenkel-chain model gives the following characteristic nonlinear behavior: (a) The two distinct dynamic modes in the relaxation modulus GS(t,lambda)--as observed in the linear region reported in Paper I [Y.-H. Lin and A. K. Das, J. Chem. Phys. 126, 074902 (2007), preceding paper]--or in the first normal-stress difference function GPsi1(t,lambda) are shown to have different strain dependences: strain hardening for the fast mode and strain softening for the slow mode. (b) The Lodge-Meissner relation GS(t,lambda)=GPsi1(t,lambda) holds over the whole time range, which has been shown both analytically and by simulation. (c) The second normal-stress difference is nonzero, being positive in the fast-mode region and negative in the slow-mode region. The comparisons between orientation and stress for all tensor components consistently confirm the strong correlation of the slow mode as well as its entropic nature with the segmental-orientation anisotropy as shown in the linear region studied in Paper I. A consequence of this correlation is the applicability of the stress-optical rule in the slow-mode region. This also leads to the expectation that the damping function h(lambda)=G(S)(t,lambda)/G(S)(t,lambda-->0) and the ratio between the first and second normal-stress differences, N2(t,lambda)/N1(t,lambda), are described by the orientation tensor which has the same form as that given by Doi and Edwards [J. Chem. Soc. Faraday Trans. 2 74, 1789 (1978); 74, 1802 (1978)] with independent-alignment approximation for an entangled system. The similarity between the slow mode of an entanglement-free Fraenkel-chain system and the terminal mode of an entangled polymer system as observed in the comparison of theory, simulation, and experiment suggests that the close correlation of the entropic nature of the mode with the orientation anisotropy--as of the Fraenkel segment or the primitive step in the Doi-Edwards theory--is a generally valid physical concept in polymer viscoelasticity.     

Lattice Monte Carlo simulations of three-dimensional charged polymer chains

The configurational properties of strongly charged polyelectrolytes accompanied by neutralizing counterions in dilute solutions are simulated using the cooperative motion algorithm on the face-centered-cubic lattice. The full Coulomb potential and the excluded volume condition between different ions/beads are taken into account and the reduced temperature T* is considered the main, variable parameter. The calculations that have been carried out for solutions of both single and several chains indicate a few regions of their behavior: (1) for T*--> infinity, it corresponds to that of neutral, self-avoiding polymers under good solvent conditions; (2) for T* approximately 1, due to the electrostatic interactions being effectively stronger, the chains are more outstretched compared to their size at other temperatures; (3) for T* well below one, the counterion condensation becomes more and more dominant, which gradually leads to strongly collapsed chains; and (4) at the lowest temperatures the chains and counterions assume low-energy configurations in the form of neutral, compact aggregates.     

Lattice Monte Carlo simulations of three-dimensional charged polymer chains. II. Added salt

The configurational properties of strongly charged polyelectrolytes accompanied by counterions and salt ions in dilute solutions are simulated using the cooperative motion algorithm on the face-centered-cubic lattice. The calculations indicate that both positive and negative ions condense on the chains at sufficiently low temperatures and their amount depends on the concentration of added salt. As the temperature decreases from high to low, the chains undergo conformational changes from neutral-like, self-avoiding polymers by more outstretched forms to compact globules. The observed extension of the chains at intermediate temperatures is also affected by the amount of salt. Furthermore, configurations with the lowest energies recorded at the lowest temperatures are aggregates of single or more entangled chains and ions of both types.     

Topological coarse graining of polymer chains using wavelet-accelerated Monte Carlo. I. Freely jointed chains

We introduce a new, topologically-based method for coarse-graining polymer chains. Based on the wavelet transform, a multiresolution data analysis technique, this method assigns a cluster of particles to a coarse-grained bead located at the center of mass of the cluster, thereby reducing the complexity of the problem by dividing the simulation into several stages, each with a fraction of the number of beads as the overall chain. At each stage, we compute the distributions of coarse-grained internal coordinates as well as potential functions required for subsequent simulation stages. In this paper, we present the basic algorithm, and apply it to freely jointed chains; the companion paper describes its applications to self-avoiding chains.     

Wavelet-accelerated Monte Carlo sampling of polymer chains

We introduce a new method for coarse-graining polymer chains, based on the wavelet transform, a multiresolution data analysis technique. This method, which assigns a cluster of particles to a coarse-grained bead located at the center of mass of the cluster, reduces the complexity of the problem significantly by dividing the simulation into several stages, each with a small fraction of the number of beads in the overall chain. At each stage, we compute the distributions of coarse-grained internal coordinates as well as potential functions required for subsequent simulation stages. We show that, with this wavelet-accelerated Monte Carlo method, coarse-grained Gaussian and self-avoiding random walks can reproduce results obtained from atomistic simulations to a high degree of accuracy in orders of magnitude less time. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 43: 897–910, 2005     

Monte Carlo simulations of chains confined inside a cube

Self-avoiding walks (SAWs) and random-flight walks (RFWs) of various lengths embedded on a simple cubic lattice have been computer generated inside cubes of varying side. If B is the side of the confining cube, we define the reduced cube side size B₀ as B₀ = (B − 1)/<r²>^1/2, where <r²>^1/2 is the root-mean-square end-to-end distance of the non-confined chains. Dimensionless diagrams are then given of the Monte Carlo estimates for the dimensions, the entropy, and the compressibility parameter PV/(kT) of the confined chains as a function of B₀. The comparative behaviour of the confined SAWs and RFWs is established, scaling properties are examined, and the Monte Carlo estimates compared with theory when such theory is available.     

Topological coarse graining of polymer chains using wavelet-accelerated Monte Carlo. II. Self-avoiding chains

In the preceding paper [A. E. Ismail, G. C. Rutledge, and G. Stephanopoulos J. Chem. Phys. (in press)] we introduced wavelet-accelerated Monte Carlo (WAMC), a coarse-graining methodology based on the wavelet transform, as a method for sampling polymer chains. In the present paper, we extend our analysis to consider excluded-volume effects by studying self-avoiding chains. We provide evidence that the coarse-grained potentials developed using the WAMC method obey phenomenological scaling laws, and use simple physical arguments for freely jointed chains to motivate these laws. We show that coarse-grained self-avoiding random walks can reproduce results obtained from simulations of the original, more-detailed chains to a high degree of accuracy, in orders of magnitude less time.     

Monte Carlo simulations of polymer dynamics: Recent advances

A brief review is given of applications of Monte Carlo simulations to study the dynamical properties of coarse-grained models of polymer melts, emphasizing the crossover from the Rouse model toward reptation, and the glass transition. The extent to which Monte Carlo algorithms can mimic the actual chain dynamics is critically examined, and the need for the use of coarse-grained rather than fully atomistic models for such simulations is explained. It is shown that various lattice and continuum models yield qualitatively similar results, and the behavior agrees with the findings of corresponding molecular dynamics simulations and experiments, where available. It is argued that these simulations significantly enhance our understanding of the theoretical concepts on the dynamics of dense macromolecular systems. © 1997 John Wiley & Sons, Inc.     

Monte Carlo simulations of single crystals from polymer solutions

A novel "anisotropic aggregation" model is proposed to simulate nucleation and growth of polymer single crystals as functions of temperature and polymer concentration in dilute solutions. Prefolded chains in a dilute solution are assumed to aggregate at a seed nucleus with an anisotropic interaction by a reversible adsorption/desorption mechanism, with temperature, concentration, and seed size being the control variables. The Monte Carlo results of this model resolve the long-standing dilemma regarding the kinetic and thermal roughenings, by producing a rough-flat-rough transition in the crystal morphology with increasing temperature. It is found that the crystal growth rate varies nonlinearly with temperature and concentration without any marked transitions among any regimes of polymer crystallization kinetics. The induction time increases with decreasing the seed nucleus size, increasing temperature, or decreasing concentration. The apparent critical nucleus size is found to increase exponentially with increasing temperature or decreasing concentration, leading to a critical nucleus diagram composed in the temperature-concentration plane with three regions of different nucleation barriers: no growth, nucleation and growth, and spontaneous growth. Melting temperatures as functions of the crystal size, heating rate, and concentration are also reported. The present model, falling in the same category of small molecular crystallization with anisotropic interactions, captures most of the phenomenology of polymer crystallization in dilute solutions.     

Monte Carlo simulations of radioactive waste embedded into polymer

Radioactive waste is generated from the nuclear applications and it should properly be managed according to the regulations set by the regulatory authority. Poly(carbonate urethane) and poly(bisphenol a-co-epichlorohydrin) are radiation-resistant polymers and they are possible candidate materials that can be used in the radioactive waste management. In this study, maximum allowable waste activity that can be embedded into these polymers and dose rate distribution of the waste drum (containing waste and the polymer matrix) were found via Monte Carlo simulations. The change of mechanical properties of above-mentioned polymers was simulated and their variations within the waste drum were determined for 15, 30 and 300 years after embedding.     

Water in nanopores. I. Coexistence curves from Gibbs ensemble Monte Carlo simulations

Coexistence curves of water in cylindrical and slitlike nanopores of different size and water-substrate interaction strength were simulated in the Gibbs ensemble. The two-phase coexistence regions cover a wide range of pore filling level and temperature, including ambient temperature. Five different kinds of two-phase coexistence are observed. A single liquid-vapor coexistence is observed in hydrophobic and moderately hydrophilic pores. Surface transitions split from the main liquid-vapor coexistence region, when the water-substrate interaction becomes comparable or stronger than the water-water pair interaction. In this case prewetting, one and two layering transitions were observed. The critical temperature of the first layering transition decreases with strengthening water-substrate interaction towards the critical temperature expected for two-dimensional systems and is not sensitive to the variation of pore size and shape. Liquid-vapor phase transition in a pore with a wall which is already covered with two water layers is most typical for hydrophilic pores. The critical temperature of this transition is very sensitive to the pore size, in contrast to the liquid-vapor critical temperature in hydrophobic pores. The observed rich phase behavior of water in pores evidences that the knowledge of coexistence curves is of crucial importance for the analysis of experimental results and a prerequiste of meaningful simulations.     

Monte Carlo simulations on short single-stranded oligonucleotides. I. Application to RNA trimers

Monte Carlo (MC) structural simulation of short RNA sequences has been carried out by random variations of the nucleotide conformational angles (i.e., phosphodiester chain torsional angles and sugar pucker pseudorotational angles). All of the chemical bond lengths and valence angles remained fixed during the structural simulation, except those of the sugar pucker ring. In this article we present the simulated structures of RNA trimers—r(AAA) and r(AAG)—obtained at 11°C and 70°C. The influence of various initial conformations (selected as starting points in the MC simulations) on the equilibrium conformations has been discussed. The simulated conformational angles have been compared with those estimated by nuclear magnetic resonance (NMR) spectroscopy. For both of the oligonucleotides studied here, the most stable structures are helical conformations with stacked bases, at 11°C and 70°C. However, when the starting point is a stretched chain, it is found that r(AAA) adopts a reverse-stacked structure at low temperature (11°C), in which the A3 base is located between the A1 and A2 bases. Although the energies of these conformations (helical and reverse stacked) are very close to each other, the potential barrier between them is extremely high (close to 30 kcal/mol). This hinders the conformational transition from one structure to the other at a given temperature (and in the course of a same MC simulation). However, it is possible to simulate this structural transition by heating the reverse-stacked structure up to 500°C and cooling down progressively to 70°C and 11°C: Canonical helical structures have been obtained by this procedure. © 1994 by john Wiley & Sons, Inc.     

Monte Carlo studies of tethered chains

Monte Carlo computer simulations of end-tethered chains grafted onto a hard wall have been performed. The chains were modeled as self-avoiding chains on a cubic lattice at athermal solvent conditions. The simulations spanned a wide range of chain lengths, N (100–1000, i.e., up to molecular weights of a few hundred thousands), and anchoring densities, σ (2 × 10⁻⁴ to 0.4), to properly chart the relevant parameter space. It is shown that the reduced surface coverage σ* = σπR is the most appropriate variable that quantitatively determines the mushroom, overlapping mushroom and brush regimes, where R_g is the radius of gyration of a free chain in solution. The simulation data are analyzed to determine the conformational characteristics and shape of the anchored chains and to compare them with the predictions of the analytical self consistent field theory. The strong stretching limit of the theoretical predictions is obtained only for σ* > 8. © 2009 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 47:2449–2461, 2009     

Nonlinear viscoelasticity of entangled polymer chains. I. Response of the slip-link model to elongation strains

The theoretical stress relaxation of polymer melts is derived from a new solution of the response of Doi's slip-link model in the plateau zone of the rubberlike part of the spectrum. The results of this treatment, based on the application of the statistical theory of elasticity to a virtual transient network of monodisperse entangled chains, are applied to the calculation of stress-strain relations for uniaxial extension performed at constant strain rate and to that of stress relaxation after such an elongation.     

Parallel Markov chain Monte Carlo simulations

With strict detailed balance, parallel Monte Carlo simulation through domain decomposition cannot be validated with conventional Markov chain theory, which describes an intrinsically serial stochastic process. In this work, the parallel version of Markov chain theory and its role in accelerating Monte Carlo simulations via cluster computing is explored. It is shown that sequential updating is the key to improving efficiency in parallel simulations through domain decomposition. A parallel scheme is proposed to reduce interprocessor communication or synchronization, which slows down parallel simulation with increasing number of processors. Parallel simulation results for the two-dimensional lattice gas model show substantial reduction of simulation time for systems of moderate and large size.     

Supramolecular polymer formation by metal-ligand complexation: Monte Carlo simulations and analytical modeling

Equilibrium metal-ligand complexation leading to formation of linear or ringlike supramolecular polymers is studied by means of Monte Carlo (MC) simulations and theoretical analysis. We found that in most of the cases high-molecular-weight polymers are formed over a rather narrow composition range (near the 2:1 ligand-metal ratio). Besides the imbalance in the number of metals and ligands, the molecular weight decrease in the metal-rich area is caused by an increase in 1:1 ligand-metal complex formation. The results of simulations and theoretical modeling show that the fraction of 1:1 complexes considerably decreases for metal-ligand pairs with a high cooperativity of complexation. On the basis of our analytical model, we suggest a simple criterion for choosing the metal/ligand pair to achieve high molecular weight complexes in a broad range of metal-rich compositions. Dilution of a solution of metallosupramolecular polymers is found to decrease the average molecular weight and to enhance ring formation, which otherwise is very limited.     

Monte Carlo simulation of tetrahedral chains, 6. Linear and star-branched polymers near to theta conditions

By use of the pivot algorithm, star-branched chains with F = 4, 8 and 12 arms of length n and linear chains (F = 2) are generated on a tetrahedral lattice (120 ⩽ nF ⩽ 3 840). By taking into account non-bonded nearest-neighbour interactions (each contact contributes an energy ϕ · kT to the total energy of the configuration) a variation of the thermodynamic quality of the solvent is simulated by a variation of the energy parameter ϕ in the range −0,425 to −0,525. The energy parameter ϕ_⊙ = −0,475, characteristic of theta conditions, was evaluated by use of an intramolecular criterion (proportionality between mean-square dimensions and total chain-length) and confirmed by a crossover scaling analysis. Theta dimensions are found to be larger than those of nonreversal random walks, the deviation growing with increasing number of arms.     

Computational methods in coupled electron-ion Monte Carlo simulations.

In the last few years, we have been developing a Monte Carlo simulation method to cope with systems of many electrons and ions in the Born-Oppenheimer approximation: the coupled electron-ion Monte Carlo method (CEIMC). Electronic properties in CEIMC are computed by quantum Monte Carlo rather than by density functional theory (DFT) based techniques. CEIMC can, in principle, overcome some of the limitations of the present DFT-based ab initio dynamical methods. The new method has recently been applied to high-pressure metallic hydrogen. Herein, we present a new sampling algorithm that we have developed in the framework of the reptation quantum Monte Carlo method chosen to sample the electronic degrees of freedom, thereby improving its efficiency. Moreover, we show herein that, at least for the case of metallic hydrogen, variational estimates of the electronic energies lead to an accurate sampling of the proton degrees of freedom.