Structure of an associating polymer melt in a narrow slit by molecular dynamics simulation

Molecular dynamics simulation has been used to study the equilibrium properties of a generic coarse-grained polymer melt with associating terminal groups, confined in a narrow slit by two atomically smooth walls. Simulations were carried out as a function of wall separation and attracting strength as well as polymer end-end interaction strength. We find that confinement has an important effect on the melt properties. In particular, strongly attracting walls can produce radical changes in chain conformation, the nature of the transient network, and the structure of the aggregates formed by the associating terminals.     

Analysis of shear-induced and extensional-induced associating polymer assemblies: Brownian dynamics simulation

In order to examine the difference between shear-induced and extensional-induced associating polymer assemblies at the molecular level, Brownian dynamics simulations with the bead-spring model were carried out for model DNA molecules with sticky spots. The radial distribution of molecules overestimates from that in the absence of flow and increases with increasing Weissenberg number in extensional flow, but slightly underestimates without regard to shear rate in shear flow. The fractional extension progresses more rapidly in extensional flow than in shear flow and the distribution of fractional extension at the formation time has a relatively sharper peak and narrower spectrum in extensional flow than in shear flow. In shear flow, the inducement of the assembly mainly results from the progress of the probability distribution of fractional extension. However, in extensional flow, the assembly is induced by both the progress of the probability distribution and increasing the values of the radial distribution.     

Characterization of polyethylene crystallization from an oriented melt by molecular dynamics simulation

Molecular dynamics is used to characterize the process of crystallization for a united atom model of polyethylene. An oriented melt is produced by uniaxial deformation under constant load, followed by quenching below the melting temperature at zero load. The development of crystallinity is monitored simultaneously using molecular-based order parameters for density, energy, and orientation. For crystallization temperatures ranging from 325 to 375 K, these simulations clearly show the hallmarks of crystal nucleation and growth. We can identify multiple nucleation events, lamellar growth up to the limit imposed by periodic boundaries of the simulation cell, and lamellar thickening. We observe a competition between the rate of nucleation, which results in multiple crystallites, the rate of chain extension, which results in thicker lamellae, and the rate of chain conformational relaxation, which is manifested in lower degrees of residual order in the noncrystalline portion of the simulation. The temperature dependence of lamellar thickness is in accord with experimental data. At the higher temperatures, tilted chain lamellae are observed to form with lamellar interfaces corresponding approximately to the [201] facet, indicative of the influence of interfacial energy.     

Equilibrium interaction of solid surfaces across a polymer melt

Forces across polymer melts are poorly understood despite their importance for adhesion and fabricating composite materials. Using an atomic force microscope (AFM), this interaction was measured for poly(dimethyl siloxane) (PDMS). The structure of the polymer at the surface changed during the first approximately 10 h. Afterward, short-range attractive forces were observed with short-chain PDMS (M(w) = 4200 g/mol). Using PDMS with a molecular weight (M(w) = 18 000 g/mol) above the entanglement limit, we measured a monotonically decaying repulsive force, which indicates that a quasi-immobilized layer had formed at the solid surface. Due to the small radius of curvature of the tip, forces could be measured in equilibrium.     

Rapid solidification of cobalt melt by molecular dynamics simulation

Molecular dynamics simulation was applied to investigating the evolvement rule of cobalt melt microstructure during solidification at different cooling rates. The cooling rate for the formation of amorphous phase is determined by analyzing the radial distribution function, the H–A bond-type index and the mean square displacement. The simulation results showed that the nucleation undercooling increases with the initial temperature, and in the undercooling versus temperature curve, there are two inflection points. Besides, when the initial temperature reaches 2450 K, the undercooling will be stabilized at 1061 K. As the cooling rate is less than 1.0?×?10^11.0 K s^?1, the FCC and HCP crystal structures will be obtained. Amorphous structure will be obtained if the cooling rate is more than 1.0?×?10^13.0 K s^?1. If the cooling rate of the Co melt is between 1.0?×?10^11.0 and 1.0?×?10^13.0 K s^?1, the crystal and amorphous structures will be coexistent, which indicates that the critical cooling rate of crystal–amorphous transition is 1.0?×?10^11.0 K s^?1.

     

A molecular dynamics computer simulation study of room-temperature ionic liquids. II. Equilibrium and nonequilibrium solvation dynamics

The molecular dynamics (MD) simulation study of solvation structure and free energetics in 1-ethyl-3-methylimidazolium chloride and 1-ethyl-3-methylimidazolium hexafluorophosphate using a probe solute in the preceding article [Y. Shim, M. Y. Choi and H. J. Kim, J. Chem. Phys. 122, 044510 (2005)] is extended to investigate dynamic properties of these liquids. Solvent fluctuation dynamics near equilibrium are studied via MD and associated time-dependent friction is analyzed via the generalized Langevin equation. Nonequilibrium solvent relaxation following an instantaneous change in the solute charge distribution and accompanying solvent structure reorganization are also investigated. Both equilibrium and nonequilibrium solvation dynamics are characterized by at least two vastly different time scales--a subpicosecond inertial regime followed by a slow diffusive regime. Solvent regions contributing to the subpicosecond nonequilibrium relaxation are found to vary significantly with initial solvation configurations, especially near the solute. If the solvent density near the solute is sufficiently high at the outset of the relaxation, subpicosecond dynamics are mainly governed by the motions of a few ions close to the solute. By contrast, in the case of a low local density, solvent ions located not only close to but also relatively far from the solute participate in the subpicosecond relaxation. Despite this difference, linear response holds reasonably well in both ionic liquids.     

Equilibrium conformational dynamics of a polymer in a solvent

Molecular dynamics simulations were used to study the conformational dynamics of a bead-spring model polymer in an explicit solvent under good solvent conditions. The dynamics of the polymer chain were investigated using an analysis of the time autocorrelation functions of the Rouse coordinates of the polymer chain. We have investigated the variation of the correlation functions with polymer chain length N, solvent density rho, and system size. The measured initial decay rates gamma(p) of the correlation functions were compared with the predictions from a theory of polymer dynamics which uses the Oseen tensor to describe hydrodynamic interactions between monomers. Over the range of chain lengths considered (N = 30-60 monomers), the predicted scaling of gamma(p) proportional to N(-3nu) was observed at high rho, where nu is the polymer scaling exponent. The predicted gamma(p) are generally higher than the measured values. This discrepancy increases with decreasing rho, as a result in the breakdown in the conditions required for the Oseen approximation. The agreement between theory and simulation at high rho improves considerably if the theoretical expression for gamma(p) is modified to avoid sum-to-integral approximations, and if the values of (R(p)2), which are used in the theory, are taken directly from the simulation rather than being calculated using approximate scaling relations. The observed finite-size scaling of gamma(p) is not quantitatively consistent with the theoretical predictions.     

一种超支化疏水缔合聚合物的制备与性能评价

首先制备了疏水单体2-丙烯酰胺基十四烷磺酸,在此基础上又以改性Si O2功能单体为反应核制备了超支化疏水缔合聚合物(HBPAM),结构经红外光谱(FTIR)和核磁共振氢谱(1H NMR)表征证实。HBPAM在低浓度时主要是分子内缔合,表观黏度低,随着浓度的增加,分子内缔合逐渐变为分子间缔合,又因其独特的三维立体网状结构,溶液黏度显著增加。与梳形KYPAM相比,HBPAM在耐温抗盐及抗剪切方面有较高的优势:升温到85℃时黏度保留率为62.2%;100000mg·L-1Na Cl、150000 mg·L-1Na Cl、500 mg·L-1Mg Cl2、1000 mg·L-1Mg Cl2、500 mg·L-1Ca Cl2、1000 mg·L-1Ca Cl2时的黏度保留率分别为233.0%、132.9%、64.4%、26.1%、66.2%、15.7%;3400rpm/min剪切30s后HBPAM的黏度保留率为69.8%,比KYPAM高10.9%。在60℃烘箱中的30d老化实验证明,HBPAM比KYPAM有明显抗老化能力,尤其是在高矿化度条件下优势更明显。     

The dynamics of single chains within a model polymer melt

Discontinuous molecular dynamics simulations are performed on a system containing 32 hard chains of length 192 at a volume fraction of phi = 0.45 to explore the idea that localized entanglements have a significant effect on the dynamics of the individual chains within an entangled polymer melt. Anomalous behavior can still be observed when studying the dynamics of the individual chains, although increased time averaging causes the anomalous relaxation-memory-release behavior that was observed previously in the system to smooth out. First, the individual chain mean squared displacements and apparent diffusion coefficients are calculated, and a wide distribution of diffusive behavior is found. Although the apparent diffusion coefficient curve averaged over all chains displays the predicted long-time diffusive behavior, the curves for the individual chains differ both qualitatively and quantitatively. They display superdiffusive, diffusive, and subdiffusive behavior, with the largest percentage of chains exhibiting superdiffusive behavior and the smallest percentage exhibiting the predicted diffusive behavior. Next, the individual chain end-to-end vector autocorrelation functions and relaxation times are determined, and a wide distribution of stress relaxation behavior is found. The times when the end-to-end vector autocorrelation functions relax completely span almost an order of magnitude in reduced time. For some chains, the end-to-end vector autocorrelation function relaxes smoothly toward zero similar to the system average; however, for other chains the relaxation is slowed greatly, indicating the presence of additional entanglements. Almost half of the chains exhibit the anomalous behavior in the end-to-end vector autocorrelation function. Finally, the dynamic properties are displayed for a single chain exhibiting anomalous relaxation-memory-release behavior, supporting the idea that the relaxation-memory-release behavior is a single-chain property.     

Molecular dynamics simulations of monodisperse/bidisperse polymer melt crystallization

We use large scale coarse‐grained molecular dynamics simulations to study the kinetics of polymer melt crystallization. For monodisperse polymer melts of several chain lengths under various cooling protocols, we show that short chains have a higher terminal crystallinity value compared to longer ones. They align at the early stages and then cease evolving. Long chains, however, align, fold into lamella structures and then slowly optimize their dangling ends for the remaining simulation time. We then identify the mechanism behind bidisperse blend crystallization. To this end, we introduce a new algorithm (called Individual Chain Crystallinity) that allows the calculation of the crystallinity separately for short and long chains in the blend. We find that, in general, bidispersity hinders crystallization significantly. At first the crystallinity of the long chain components exceeds that of the monodisperse melt, but subsequently falls below the corresponding monodisperse melt curve after a certain “crossover time.” The time of the crossover can be attributed to the time required for the full crystallization of the short chains. This indicates that at the early stages the short chains are helping long chains to crystallize. However, after all short chains have crystallized they start to hinder the crystallization of the long chains by obstructing their motion. Lastly, polymer crystallization upon various thermodynamic protocols is studied. Slower cooling is found to increase the crystallinity value. Upon an instantaneous deep quench and subsequent isothermal relaxation, the crystallinity grows rapidly with time at early stages and subsequently saturates. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 2318–2326     

Mixed solutions of an associating polymer with a cleavable surfactant

Mixtures of hydrophobically modified hydroxyethyl cellulose (HMHEC) and alkali-sensitive cleavable betaine ester surfactants have been studied by viscometry, 1H NMR, absorbance measurements, and birefringence determinations. Before the hydrolysis, the surfactants behaved as conventional nondegradable surfactants in terms of the effect on the viscosity of increasing surfactant concentration. As the surfactants were hydrolyzed, systems with time-dependent viscosity were obtained. The viscosity either decreased monotonically or went through a maximum as a function of time, depending on the initial surfactant concentration. Different surfactant chain lengths gave rise to different viscosity profiles. The rate of hydrolysis, and thus the time-dependency of the surfactant concentration, could be controlled by changing the pH of the solution.     

Molecular arrangements and conformations of liquid unbranched alkanes in narrow slits

Realistic, atomistic models of liquid tridecane in broad slits (>3 nm) and in narrow slits of thickness 1,2 nm and 1,0 nm have been obtained using the Monte Carlo technique. The setup of the models is such that the molecules in the slits are in equilibrium with the bulk liquid. The surfaces of the plates are modelled as two-dimensional arrays of hexagonally packed units having the same size and interaction parameters of a methylene group. The regions adjacent to the plates in slits with thickness > 3 nm are characterized by a well defined tendency to form partially ordered layer structures, while molecules at a distance from the plates larger than 1,5 nm are unperturbed. The simultaneous presence of two plates increases the tendency to form layer structures when their distance is 1,2 nm, while this tendency is almost totally destroyed when the slit is squeezed down to a thickness of 1,0 nm. This is also associated with a 10% decrease of the density in the latter slit.     

Brownian dynamics simulation of grafted polymer brushes

We present results of computer simulations by the method of Brownian dynamics of polymeric brushes attached to impenetrable planes. For testing both model and method we have used one polymer brush attached to a repulsive plane and compare some results with Monte Carlo results of Lai and Binder on the bond fluctuation model. We have also studied two polymeric brushes attached to two parallel planes at different distances between planes, and investigate the interplay between the interpenetration of the brushes and the configurational properties of the grafted chains.     

Decoupling effects in computer simulation of liquid-state molecular dynamics

Decoupling effects in liquid-state molecular dynamics are unambiguously proven by computer simulations. These can be discribed as the weakening of the dissipative interaction between the system if interest (e.g. molecular angular valocity) and its thermal bath caused by an intense external field of force. The new phenomenon may be used to confirm or challenge the validity of simple models of the liquids state of molecular matter and provide information on the microscopic time scale.     

A molecular dynamics computer simulation study of room-temperature ionic liquids. I. Equilibrium solvation structure and free energetics

Solvation in 1-ethyl-3-methylmidazolium chloride and in 1-ethyl-3-methylimidazolium hexafluorophosphate near equilibrium is investigated via molecular dynamics computer simulations with diatomic and benzenelike molecules employed as probe solutes. It is found that electrostriction plays an important role in both solvation structure and free energetics. The angular and radial distributions of cations and anions become more structured and their densities near the solute become enhanced as the solute charge separation grows. Due to the enhancement in structural rigidity induced by electrostriction, the force constant associated with solvent configuration fluctuations relevant to charge shift and transfer processes is also found to increase. The effective polarity and reorganization free energies of these ionic liquids are analyzed and compared with those of highly polar acetonitrile. Their screening behavior of electric charges is also investigated.