Next Article in Journal
Polysarcosine-Based Lipids: From Lipopolypeptoid Micelles to Stealth-Like Lipids in Langmuir Blodgett Monolayers
Next Article in Special Issue
A Review of Multiscale Computational Methods in Polymeric Materials
Previous Article in Journal
Polyetheretherketone Hybrid Composites with Bioactive Nanohydroxyapatite and Multiwalled Carbon Nanotube Fillers
Previous Article in Special Issue
Atomistic Modelling of Confined Polypropylene Chains between Ferric Oxide Substrates at Melt Temperature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrodynamic Interactions and Entanglements of Polymer Solutions in Many-Body Dissipative Particle Dynamics

Department of Mechanical Engineering, State University of New York at Binghamton, Binghamton, NY 13902, USA
Polymers 2016, 8(12), 426; https://doi.org/10.3390/polym8120426
Submission received: 16 November 2016 / Revised: 5 December 2016 / Accepted: 6 December 2016 / Published: 9 December 2016
(This article belongs to the Special Issue Computational Modeling and Simulation in Polymer)

Abstract

:
Using many-body dissipative particle dynamics (MDPD), polymer solutions with concentrations spanning dilute and semidilute regimes are modeled. The parameterization of MDPD interactions for systems with liquid–vapor coexistence is established by mapping to the mean-field Flory–Huggins theory. The characterization of static and dynamic properties of polymer chains is focused on the effects of hydrodynamic interactions and entanglements. The coil–globule transition of polymer chains in dilute solutions is probed by varying solvent quality and measuring the radius of gyration and end-to-end distance. Both static and dynamic scaling relations for polymer chains in poor, theta, and good solvents are in good agreement with the Zimm theory with hydrodynamic interactions considered. Semidilute solutions with polymer volume fractions up to 0.7 exhibit the screening of excluded volume interactions and subsequent shrinking of polymer coils. Furthermore, entanglements become dominant in the semidilute solutions, which inhibit diffusion and relaxation of chains. Quantitative analysis of topology violation confirms that entanglements are correctly captured in the MDPD simulations.

1. Introduction

Because of the fast evolution of polymer science and nanotechnologies, polymeric materials with precisely tailored functionalities, such as block copolymer thin films and polymer nanocomposites, become central to a wide range of emerging applications. In order to achieve versatile properties and functionalities of the macroscopic materials, polymeric building units with increasingly smaller sizes and different types of chemistry are integrated into more and more complex structural hierarchy in controlled and synergistic manners. This demand challenges manufacturing processes, in which one must close the gap between adjusting macroscopic characteristics and controlling micro-nanostructures of building blocks.
Among many existing methods for polymer synthesis and processing, solution deposition techniques (e.g., drop casting, dip coating, spin coating, electrospray, and printing) [1,2,3,4,5,6] are frequently employed due to the ease of implementation and low cost. Namely, an initial liquid precursor with solutes and/or dispersed moieties transforms into a solid material upon evaporation. The thermodynamics and kinetics of the precursor solutions during the drying process govern the final structure of the system [7,8,9,10,11]. However, the complex interplay among polymer, solvent, and other components in the evaporating mixtures is not fully understood; in-situ experimental characterization of non-equilibrium multi-component systems across multiple length and time scales is still difficult.
Computational simulation is appealing for elucidating the intricate polymer–polymer and polymer–solvent interactions and revealing detailed dynamics of the multi-component systems on the molecular scale. Dissipative particle dynamics (DPD) emerges as a popular mesoscopic technique for modeling polymer solutions/melts and probing their structural and rheological properties [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]. A recently developed variation, many-body dissipative particle dynamics (MDPD), further extends DPD’s capabilities to producing the liquid–vapor coexistence and simulating evaporation process [30,31,32,33,34,35], which is a promising step toward modeling solution-based processing of polymeric materials. Despite many successful MDPD simulations [35,36,37,38,39,40], there is no report on the parameterization of MDPD interaction parameters [37,41] for the liquid–vapor coexistent systems. More importantly, a comprehensive study of the conformation and dynamics of polymer chains in solutions is lacking.
This work seeks to establish a general mapping of MDPD parameters onto the Flory–Huggins χ parameter for the liquid state with coexisting vapor, which enables physical representations of multi-component, multi-phase polymeric systems, e.g., polymer solutions, copolymers, and polymer blends. Based on the obtained relation, polymer solutions with a full range of concentration and solvent quality are modeled using MDPD, followed by systematic characterization of the structural and dynamic properties of solvated polymer chains. The analysis is focused on revealing hydrodynamic interactions and quantifying entanglements of chains.
The present article is organized as follows: Section 2 briefly reviews the MDPD method and describes the details of the polymer simulations. Section 3.1 presents the parameterization of the MDPD interactions based on the mean-field Flory–Huggins theory. In Section 3.2, the conformation and dynamics of the polymer chains in dilute solutions with varying solvent quality are probed and the simulation results are compared to the Zimm model. Section 3.3 focuses on semidilute solutions and the effect of polymer volume fraction on the chain morphology and dynamics, through which the roles of hydrodynamic interactions and chain entanglements are elucidated. Finally, a summary of the important results and conclusions is presented in Section 4.

2. Methods

Many-body dissipative particle dynamics [30,32,33,34,41,42] is used to model polymer solutions with volume fractions ranging from 0.008 to 0.7, representing dilute solutions to semidilute solutions. MDPD is a mesoscopic particle-based method that can effectively model multi-phase, multi-component systems and captures correct hydrodynamic behavior [35,39,40,43,44,45,46,47,48,49]. In MDPD, a volume of fluid is modeled by coarse-grained beads, and each bead represents a cluster of molecules. The evolution of the entire system over time is dictated by the motion of beads, which is governed by Newton’s equation of motion, m d v i / d t = f i . The thermodynamics and transport properties of the system are determined statistically via the ensemble of MDPD beads.
The force acting on each bead i from neighboring beads consists of three parts f i ( t ) = ( F i j C + F i j D + F i j R ) , each of which is pairwise additive. The three terms, respectively, describe the conservative, drag, and random forces. The sum runs over all beads j within a cutoff radius r c from bead i. The main advantage of MDPD compared to standard DPD is its ability to produce the coexistence of liquid and vapor phases in coarse-grained fluids, which is achieved by introducing a long-range attractive force that is responsible for surface tension. In particular, the conservative force F i j C is given by F i j C = A i j ( 1 r i j / r c ) r ^ i j + B i j ( ρ ¯ i + ρ ¯ j ) ( 1 r i j / r d ) r ^ i j with r i j = | r i r j | / r c and r ^ i j = ( r i r j ) / | r i r j | [32]. Here, The attraction strength A i j < 0 , the repulsion strength B i j > 0 , and the repulsion range r d < r c , making repulsion short-range and attraction long-range. The repulsive term depends on local densities ρ ¯ i and ρ ¯ j as well as inter-bead distance r i j ; the attractive term depends only on the distance. The interaction parameters A i j and B i j are given in terms of k B T / r c , where k B T is the Boltzmann constant and T is the temperature of the system. The local density for each particle is defined as ρ ¯ i = j i 15 / 2 π r d 3 ( 1 r i j / r d ) 2 [32]. This many-body conservative force was shown to produce cubic pressure-density equations of state with van der Waals loop, thereby permitting the liquid–vapor coexistence with a sharp interface [30,32,34,35]. Notably, the soft-core repulsive term allows a degree of overlap between beads, in contrast to hard-core potentials (e.g., the Lennard–Jones potential) in which the repulsion diverges at zero inter-bead distance. Thus, larger time steps than those typically used in molecular dynamics (MD) simulations, which commonly involve hard-core potentials, can be applied in the MDPD simulations. Together with the coarse-grained representation of multi-component, multi-phase systems, MDPD can capture physical phenomena occurring on relatively larger length and time scales than those normally captured by MD. These features make MDPD an ideal computational tool to resolve the polymer dynamics in liquid–vapor multi-phase systems.
The drag force is F i j D = λ ω D ( r i j ) ( r ^ i j v i j ) r ^ i j , where λ is a simulation parameter related to the viscosity arising from the interactions between the constituent beads of fluid. ω D is a weight function satisfying ω D = 0 at r i j = r c , and the relative velocity is v i j = v i v j . The random force is F i j R = σ ω R ( r i j ) ξ i j r ^ i j , where ξ i j is a zero-mean Gaussian random variable of unit variance and σ is the amplitude of the noise. The fluctuation-dissipation theorem relates σ to λ as σ 2 = 2 k B T λ [50]. Finally, the weight functions take the following form: ω D ( r i j ) = ω R ( r i j ) 2 = ( 1 r i j / r c ) 2 for r i j < r c [15]. The combination of drag force and random force serves as a thermostat applied on the MDPD system for producing the canonical (NVT) ensemble [14,15,50]. All three forces act on pairs of neighboring beads such that momentum is conserved locally, and hydrodynamic behavior emerges in relatively small systems [12,14].
The equation of motion is integrated in time using the velocity-Verlet algorithm. The simulation takes r c as the characteristic length scale with a dimensionless value as r c = 1 . The corresponding energy scale k B T at room temperature is chosen as the characteristic energy. Thus, k B T = 1 at T = 25   ° C is considered in this work. The characteristic mass is defined as the mass of a MDPD bead and takes the dimensionless value of 1. The characteristic time scale is then defined as τ = m r c 2 / k B T = 1 . A time step Δ t = 0.01 is used for all simulations characterizing static and dynamic properties of polymer [15,51]. Unless otherwise stated, the value of λ is chosen as 4.5 to obtain a relatively rapid equilibration of the system temperature and to ensure the numerical stability of the simulations for the specified time step. All simulations are performed by modifying and extending the particle dynamics software code LAMMPS [52].
In this work, the MDPD non-bonded parameters of A i j = 40 and B i j = 40 are used for any two beads of the same type (i.e., solvent–solvent or polymer–polymer interactions). Notably, MDPD applied to multi-component systems requires a constant B i j to ensure a conservative many-body force law [53]. Thus, the attraction parameter A i j is adjusted to control the polymer–solvent interactions, which can be mapped to the Flory–Huggins χ parameter as detailed below. The repulsion cutoff radius is set as r d = 0.8 for all non-bonded interactions. As used in our previous study, this parameter set yields a liquid–vapor system with a liquid number density of 3.926 and viscosity of 3.510 [35]. Therefore, the total density of polymer solutions is set to ρ sys = 3.926 for all simulations conducted in this study.
The widely used bead-spring model is applied to represent flexible polymer chains. Groups of repeating units in a linear polymer are modeled as consecutive MDPD beads connected by Hookean springs, whose potential is given by E bond = 1 2 K b ( r i j r 0 ) 2 . Here, K b = 128 is the elastic constant, r i j is the distance between bonded beads i and j, and r 0 is the equilibrium bond distance. An equilibrium bond distance of 0.685 is chosen such that the mean distance between connected beads coincides with the position of the first neighbor peak in the radial distribution function of the MDPD liquid with density 3.926. Initial configuration of each polymer chain is created by a self-avoiding random walk on a face-centered-cubic (fcc) lattice superimposed on the simulation box [54,55]. Polymer beads in a chain are consecutively inserted onto the lattice in a stepwise manner, mimicking chain-growth polymerization. The probability of successful insertion and chain growth depends on the generated bond angle and the density of unoccupied sites. Once the polymers are constructed, the corresponding number of solvent beads are added to the system at randomly positions to solvate the polymer and obtain the desired polymer volume fraction.
The dimensions of the simulation box vary from 10 × 10 × 10 to 40 × 40 × 40, with periodic boundary conditions imposed in all three directions. The total number of beads in our simulations ranges from 3000 to 251,264. The number of beads per chain is varied from N = 5 to 140. In order to minimize the finite-size effect on the equilibrium properties of the system [21,22,23], the length L of the cubic simulation box satisfies the condition L / R g > 5 , where R g is the radius of gyration of the polymer chains. Each simulation system is typically equilibrated for 5 × 105 time steps prior to a production run of at least 5 × 105 time steps.

3. Results and Discussion

3.1. Parameterization of MDPD Parameters

In order to represent physical systems of specific chemistry, one must establish a rigorous approach to the parameterization of the interaction parameters. Unlike the commonly applied DPD method [15,56], to date only a few MDPD studies have reported general methods for obtaining a corresponding parameter set capable of reproducing the thermodynamic properties of experimental fluids [37,41]. The central idea of the MDPD parameter derivation is mapping the configurational part of the free energy density of MDPD binary mixtures to the mixing free energy density from the mean-field Flory–Huggins lattice theory [57], similar to the approach applied to DPD method [15]. Due to the no-go theorem of MDPD ( B a a = B b b = B a b for a binary mixture of a and b fluids), the cubic repulsive term in the MDPD equation of state of a single component fluid [32] p = ρ k B T + α A ρ 2 + 2 α B r d 4 ( ρ 3 c ρ 2 + d ) does not contribute to the configurational part, where α, c, and d are fitting parameters [32]. Thus, the mapping yields the same expression of the Flory–Huggins χ parameter as the one for DPD, χ = 2 α ( ρ a + ρ b ) ( A ab A a a ) / k B T = 2 α ρ sys Δ A . A set of simulations of binary mixtures of monomers is carried out to establish the χA relation, where only the cross attraction parameter A ab is adjusted from −14 to −7. The χ parameter is calculated from the density profiles of the phase-separated systems. Only systems showing strong phase separation are sampled because the Flory–Huggins mean-field expression breaks down for systems with small χ [15]. A good linear relation between χ and excess repulsion ΔA is confirmed in Figure 1a; the fitting results in χ = ( 0.479 ± 0.007 ) Δ A for ρ sys = 3.926 . Notably, the MDPD method can be readily extended to simulate ternary systems given that the parameterization of the self and cross interaction parameters is carried out in a rigorous way. In addition to the Flory–Huggins approach employed in this work, Hildebrand solution theory approach [58,59] and more advanced ab-initio methodologies [60] may be adapted in MDPD to simulate ternary mixtures of beads dissimilar in both size and chemistry.
The surface tension of the binary mixtures can be also determined in the simulations. The Irving–Kirkwood expression is used to measure surface tension, given as σ = L z 2 p z z 1 2 ( p x x + p y y ) , where p x x , p y y , and p z z are the diagonal components of the microscopic pressure tensor [61] and L z is the dimension of the simulation box in the z direction. Here, the z direction is normal to the time-averaged interfaces. The factor of 1/2 arises from the existence of two interfaces in the periodic simulation domain. The angular bracket represents ensemble average. The dependence of surface tension on the χ parameter is subsequently compared to the classical van der Waals theory σ χ α ( 1 χ crit / χ ) 1.5 . Figure 1b shows good agreement between the simulation results and the theory. The variation of surface tension is given by the fitting σ ( 0.82 ± 0.04 ) χ 0.29 ± 0.02 [ 1 ( 2.34 ± 0.03 ) / χ ] 1.5 . Consistent with the DPD results [15], the extrapolation of the surface tension of the MDPD fluids leads to a non-classical critical point of a value of 2.34 higher than the theoretical value χ crit = 2 . Nevertheless, the critical point found from the MDPD simulations agrees well with the one obtained in the DPD simulations.

3.2. Dilute Polymer Solutions for Different Solvent Quality

Dilute polymer solutions with a volume fraction ϕ 0.008 are investigated. Each system has M chains with N beads per chain. The inter-molecular interactions between different chains are present in the system, closely representing physical systems. The number of chains is adjusted accordingly to keep the polymer volume fraction constant when the chain length is varied. The solvent quality is controlled by tuning the polymer–solvent attraction parameter A ps , which determines the polymer–solvent Flory–Huggins parameter χ ps according to the linear relation obtained above. The mean-square radius of gyration of the polymer chains and their end-to-end distance are calculated to characterize the equilibrium structure of polymer chains, defined by R g 2 = ( 1 / N ) i = 1 N ( r i r com ) 2 and R 1 N 2 = ( r 1 r N ) 2 , respectively. Here, r i is the position of the ith bead in a polymer chain and r com is the center-of-mass position of the chain. The coil–globule transition of the flexible chains having different lengths is monitored as the solvent quality varies. Figure 2 clearly shows the characteristic collapse of polymer chains from an expanded coil state to a globule state as A ps decreases, modeling the solvent quality change from good to poor. The observed continuous coil–globule transitions are consistent with theory predictions [62,63], where longer chains exhibit much more pronounced change in size. The spatial distribution of chains also reveals that the conformation transition is accompanied by the phase separation due to inter-molecular interactions. In particular, multiple collapsed chains aggregate into a bigger globule as shown in Figure 2b. Figure 2a also indicates that the Flory theta point of the transition corresponding to the ideal chain behavior does not coincide with the athermal condition of A ps = A ss = A pp = 40 .
To explicitly analyze the correspondence of solvent quality to A ps , Figure 3 plots the variation of mean-square radius of gyration as the chain length increases on the log–log scale for three different values of A ps . The linear behavior of R g 2 as a function of bond number confirms that the MDPD results can be well described by the classical scaling law [64,65], R g 2 ( N 1 ) 2 ν , where ν is the Flory scaling exponent. The good agreement between the simulation data and the theory indicates that, at low polymer volume fractions, the presence of polymer–polymer inter-molecular interactions have negligible effect on the polymer conformation. The power-law fitting of the simulation data for the athermal solvent ( A ps = 40 ) yields R g 2 = ( 0.066 ± 0.008 ) ( N 1 ) 1.22 ± 0.03 . The end-to-end distance also obeys the scaling relation, and the fitting results in R 1 N 2 = ( 0.49 ± 0.09 ) ( N 1 ) 1.17 ± 0.04 . The scaling exponents obtained from both measurements are consistent: v = 0.61 ± 0.02 from the radius of gyration and v = 0.59 ± 0.02 from the end-to-end distance. The exponents imply that the chains exhibit good solvent behavior in the athermal solvent, where effective excluded volume interactions are imparted through the MDPD conservative force. This result is consistent with the DPD simulations of dilute solutions of flexible “bead-spring” chains [21,23], but is in contrast to recent characterization by Jamali et al. [41] for MDPD fluids, in which the athermal solvent leads to the ideal-chain exponent of 0.5.
A prominent feature of MDPD is the many-body nature of the repulsive part of the conservative force, which depends on the local density. Due to the athermal condition of the attractive force, the excluded volume effect must be induced by the difference in repulsions, leading to energetically favorable polymer–solvent interactions. Moreover, constant repulsion strength suggests that the preferential interaction between polymer and solvent beads is caused by the local density difference. In order to provide quantitative evidence, the structural property of the polymer solution in the athermal condition is further examined by the coordination numbers and radial distribution functions (the calculations include the bonded, nearest neighbor beads). Figure 4 shows the time-averaged coordination numbers of polymer and solvent beads up to the repulsion cutoff r d = 0.8 . Although the polymer chains are formed by matching the equilibrium bond distance to the nearest neighbor distance of the solvent, the coordination number of the polymer beads is still slightly higher than that of the solvent beads, which is caused by the existence of the bonded neighbors for polymer beads. The local densities of a bead should scale with its coordination number. Thus, the repulsion between two polymer beads will be larger than between a polymer bead and a solvent bead, leading to unfavorable polymer–polymer interactions and the effective excluded volume effect. The radial distribution functions in Figure 4 also confirm that the polymer beads interact preferentially with the solvent beads.
The reduction of polymer–solvent attraction strength increases the χ parameter and cancels the excluded volume effect, leading to the collapse of the chains. The transition from good to poor solvent occurs near A ps = 38.5 , where the scaling exponent for radius of gyration becomes v = 0.51 ± 0.01 . The polymer coils behave like ideal chains in this theta solvent except for small N, as the ratio of the mean-square end-to-end distance to the mean-square radius of gyration R 1 N 2 / R g 2 6 (see Table 1). For even larger value of A ps , the scaling exponent decreases to v = 0.29 ± 0.03 , corresponding to the poor solvent condition and the associated globule state. Notably, the behavior of short polymers in theta and poor solvents deviates from the power-law predictions. In Figure 3, the slope of the R g 2 curve is significantly higher than 0.3 for short chains in the poor solvent ( A ps = 36 ). This implies that short chains are more extended than the ideal chains described in the theory, which is attributed to the intrinsic chain rigidity induced by the many-body repulsion. The results show that the scaling regime in MDPD generally requires chains having more than 10 beads. Thus, the effective persistence length of the bead-spring chains in MDPD is longer than that in DPD, where even chains with as few as five beads were found to follow the scaling relation well [20,22].
The excluded volume effect and associated entropic cost for solvation of polymer can be interpreted in the light of the scaled particle theory [66,67]. The theory considers a chain with N beads of diameter d in a solvent of beads with the same diameter, where the volume of a single bead is v = π d 3 / 6 . The polymer chains do not interact in dilute solutions—each independently occupies a volume characterized by R g . Therefore, the mean volume fraction of polymer is ϕ p N d 3 / R g 3 and the solvent volume fraction is ϕ s = 1 ϕ p . Solvation of polymer requires work to form a cavity in the solvent to accommodate the solute molecules. Considering consecutive insertions of chain beads, the total work of cavity formation for a chain is given by W p ( N , d ) = j = 0 N 1 W j ( d ) , where j is the number of previously inserted beads. Upon the insertion of the jth bead, the chain occupies volume fraction ϕ p , j = ϕ p j / N d j / ( N 1 ) 3 ν , and the total occupied volume faction is ϕ j = ϕ p , j + ϕ s . The work for individual bead W j ( d ) is given by the scaled particle theory [67] as W j ( d ) / k B T = ln ( 1 + y j ) + 9 y j + 15 y j 2 / 2 + 3 y j 3 , where the ratio of occupied to unoccupied volume is y j ϕ j ( 1 ϕ j ) . According this expression, the entropic cost for cavity formation increases as the solvent quality varies from good to poor, quantified by the scaling exponent ν. The increase of chain length N leads to smaller entropic cost regardless of the solvent quality.
Another measurement of the conformation and stiffness of polymer chains is the dimensionless characteristic ratio, defined as the ratio of the measured mean-square end-to-end distance to the value of an ideal chain C N 1 = R 1 N 2 / [ ( N 1 ) l 2 ] . Here, l 2 is the mean-square bond length. As the chain length increases, the characteristic ratio approaches a limiting value C as shown in Figure 5. Namely, the characteristic ratio decreases in a poor solvent and increases in theta and good solvents. Table 1 tabulates the characteristic ratio of chains with different lengths in a theta solvent. The asymptotic value C is approximately 1.4, which is slightly larger than the value of 1.27 in the DPD simulations [51]. This provides additional evidence of the enhanced excluded volume interactions caused by the MDPD force law compared to the DPD model, which effectively increases chain rigidity.
One major advantage of DPD and MDPD simulations over other coarse-grained techniques such as Monte-Carlo and self-consistent field theory is the ability to simulate dynamical systems. The dynamical properties of polymer chains in the theta solution are probed through diffusion coefficient and relaxation time. The diffusion coefficient can be measured by calculating the mean-square displacement of the center-of-mass coordinate for each chain, [ r com ( t ) r com ( 0 ) ] 2 , and utilizing the Einstein relation, [ r com ( t ) r com ( 0 ) ] 2 = 6 D t . The value is extracted from the linear fitting of the mean-square displacement. Notably, the diffusion behavior of polymer can still suffer from strong finite-size effect even when L / R g > 5 [23], but the influence on the scaling law of the diffusion coefficient is negligible. Figure 6a shows the diffusion of polymer chains in MDPD simulations with a fitting D N 0.52 ± 0.05 for the theta solvent. The scaling exponent is consistent with the Flory exponent v = 0.51 ± 0.01 obtained from the mean-square radius of gyration. The diffusion of polymer chains in dilute solutions obeys the Zimm model, D N υ [64]. Therefore, the MDPD simulations correctly capture hydrodynamic interactions considered in the Zimm model. In addition, the relaxation time is another important probe of the polymer dynamics, which can be extracted from the conformational autocorrelation functions. Herein, the autocorrelation function of the end-to-end vector of the chains R 1 N = r 1 r N is calculated, defined as C ( t ) = R 1 N ( t + t 0 ) R 1 N ( t 0 ) / R 1 N 2 . The angle bracket means average over all chains and many time origins, t 0 . The longest relaxation time τ is obtained by fitting the autocorrelation function to an exponential form, C ( t ) = C 0 exp ( t / τ ) . The relaxation times of polymer chains having different lengths in the theta solvent are shown in Figure 6b. Apparently, the linearity of the data point breaks down for long chain lengths, possibly due to the finite-size effect and interactions between image chains across the periodic boundaries. The best fitting of the relaxation time excluding chains of length N 80 yields τ = ( 0.73 ± 0.06 ) N 1.52 ± 0.02 . The Zimm model predicts that τ N 3 ν ~ N 1.5 for theta solvents, while the Rouse model overestimates the relaxation time as τ N 2 . Again, the simulation results confirm the Zimm dynamics for the modeled polymer chains.

3.3. Semidilute Solutions and Entanglements

In this section, the results of semidilute polymer solutions with volume fractions ranging from 0.02 to 0.7 are presented. The number of chains in this series of simulations is 100 and the chain length is 30. The chains of this length exhibit reptation behavior in previous simulations [51,68]. Due to considerably stronger density-dependent repulsion between MDPD beads than the standard DPD repulsion, an MDPD system exhibits thermodynamic behavior quite different from that of a standard DPD simulation [51]. Notably, the thermodynamics of polymer solutions with high volume fractions or even polymer melts is significant influenced by the MDPD force law. Consequently, the MDPD simulations typically require a stronger thermostat coupling to achieve stable temperature control. This demands elevated thermostat parameters σ and λ (obeying the fluctuation-dissipation relation). Using the standard parameter λ = 4.5 ( σ = 3 ), the maximum temperature increase of 33% is observed in a polymer melt system with ϕ = 1.0 , which completely shifts the thermodynamic state of the system. For semidilute polymer solutions, a value of λ = 50 ( σ = 10 ) is applied to maintain the temperature within 4% of the set value [51].
Figure 7 shows the dependence of the polymer conformational and dynamic properties on the volume fraction in semidilute solutions. The excluded volume interactions are screened when the polymer chains overlap with each other at high volume fractions. Therefore, the swelling of the coils in a good solvent gradually diminishes as φ increases and can finally vanish in the melt. The reduction of the degree of swelling, which is characterized by the mean-square radius of gyration, is clearly depicted in Figure 7a. The simulation observation agrees well with the theoretical prediction. In contrast, the chain morphology in a theta solvent does not exhibit φ dependence, since the coils always behave as ideal chains in the theta condition.
As the concentration increases, the polymer dynamics undergo a transition from the Zimm regime to the reptation regime, in which the onset of topological constraints (entanglements) between chains severely hinders thermal motion of polymers. Consequently, the diffusion and relaxation of chains are slowed down significantly, as shown in Figure 7b. The MDPD results show the relaxation time τ increases by a factor of 2 at the highest volume fraction investigated here. This increase is more pronounced than the factor of 1.4 observed in comparable DPD simulations [23]. However, the value is consistent with previous MD simulations [69]. Future analysis on entanglements is performed to provide additional insights into the dynamics behavior of polymer chains in the MDPD model. For coarse-grained simulations, there have been serious concerns about unphysical bond crossings of the polymer chains (referred to as topology violations below) associated with the use of soft-repulsive potentials [51,70,71,72,73,74,75]. An exemplary topology violation is illustrated in Figure 8. In order to accurately predict behavior of the polymers in semidilute and concentrated solutions, chain entanglements must be properly captured in the MDPD model.
Quantitative measurement of topology violations is conducted by calculating the minimum distance vector between bonds, d i j [23,51,70,76]. The method has been detailed elsewhere [70]. The angle α between the two minimum distance vectors defined at time t and t + Δt is given by the dot product, cos α = d ^ i j ( t ) d ^ i j ( t + Δ t ) , where d ^ i j = d i j / | d i j | is the normalized unit vector. A topology violation is detected when the direction of d ^ i j = d i j / | d i j | changes by 90° or more over a single time step, i.e., α > 90 ° [23,51,76]. As pointed out in previous studies [23,76], this cross product rule may declare a topology violation based on a “false positive” when the minimum distance vector point passes through the vicinity of one of the bonded beads. It is possible in this case that the angle α exceeds 90° but no bond crossing occurs. These false cases are carefully excluded in the present analysis [23]. Figure 9 plots the relative number of topology violations per 100τ for semidilute solutions with different polymer fractions. As evident in Figure 9, the number of topology violations increases as φ increases. The concentrated solution at ϕ = 0.7 generates topology violations 2.5 times as much as those in the dilute solutions. The inset of Figure 9 also demonstrates the number of topology violations occurring with a specific value of angle α. The results indicate that a majority of the bond crossings occur with a change in orientation by angles close to 180°. The distribution is consistent with the result of DPD simulations [75]. Importantly, the number of topology violations per 100τ in the MDPD simulations is on the order of 102 (see Table 2), which is two orders of magnitude smaller than the values measured in the standard DPD simulations (~104) [23,51]. This implies that the MDPD force law is more effective than the standard DPD force in preventing topology violations and capturing entanglements, which is attributed to the stronger repulsion between beads. Additional bond–bond repulsive interactions may be introduced to further reduce the number of topology violations [51,70,75,76].

4. Conclusions

This work focuses on modeling polymer solutions on the mesoscale using many-body dissipative particle dynamics. The polymer volume fraction of the solution is varied in a wide range from 0.008 to 0.7, corresponding to a gradual transition from dilute to semidilute regimes. A systematic parameterization of MDPD interactions is conducted for MDPD systems with liquid–vapor coexistence. By mapping onto the Flory–Huggins theory, a linear relation between the difference in MDPD attraction parameter Δ a and the Flory–Huggins χ parameter is obtained, similar to previous DPD and MDPD studies. The surface tension of MDPD binary mixtures as a function of the χ parameter agrees well with the classical van der Waals theory, despite a higher critical point found from the fitting than the classical value.
The static and dynamic scaling relations are characterized for dilute polymer solutions, where the solvent quality is controlled by the polymer–solvent attraction parameter. Through the measurements of the radius of gyration and end-to-end distance, the scaling exponents indicate that polymer chains exhibit good solvent behavior in an athermal solution. The effective excluded volume interactions are induced by the density-dependent repulsion of the MDPD force law. The dynamic scalings of the diffusion coefficient and relaxation time of chains are consistent with the Zimm theory, which verifies the existence of hydrodynamic interactions in MDPD.
As the volume fraction increases and the solution enters the semidilute regime, the spatial overlapping of polymer coils results in the screening of excluded volume effect and the decrease of the radius of gyration. Meanwhile, the increasingly dominant inter-molecular interactions and chain entanglements significantly inhibit polymer dynamics. Due to the soft repulsion in MDPD, quantitative analysis is carried out to measure unphysical topology violations. The results confirm that entanglements are effective captured in the MDPD simulations.
Finally, it is noteworthy that the MDPD model of polymer solutions provides a powerful tool for efficiently simulating polymeric systems with liquid–vapor coexistence. This is particularly appealing for exploring the evaporation-driven dynamics of polymer chains and other inclusions on the mesoscale, which has significant implications for exerting better control on the processing of polymer thin films and polymer nanocomposites.

Acknowledgments

The author acknowledges partial supports from the National Science Foundation under grant No. CMMI-1538090 and the American Chemical Society Petroleum Research Fund under grant No. 56884-DNI9.

Author Contributions

Xin Yong conceived the study, developed computer simulations, performed analysis, and wrote the manuscript.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Hench, L.L.; West, J.K. The sol-gel process. Chem. Rev. 1990, 90, 33–72. [Google Scholar] [CrossRef]
  2. Brinker, C.J.; Lu, Y.; Sellinger, A.; Fan, H. Evaporation-induced self-assembly: Nanostructures made easy. Adv. Mater. 1999, 11, 579–585. [Google Scholar] [CrossRef]
  3. Eslamian, M. Spray-on thin film PV solar cells: Advances, potentials and challenges. Coatings 2014, 4, 60–84. [Google Scholar] [CrossRef]
  4. Derby, B. Inkjet printing of functional and structural materials: Fluid property requirements, feature stability, and resolution. Annu. Rev. Mater. Res. 2010, 40, 395–414. [Google Scholar] [CrossRef]
  5. Faustini, M.; Boissière, C.; Nicole, L.; Grosso, D. From chemical solutions to inorganic nanostructured materials: A journey into evaporation-driven processes. Chem. Mater. 2014, 26, 709–723. [Google Scholar] [CrossRef]
  6. Han, W.; Lin, Z. Learning from “coffee rings”: Ordered structures enabled by controlled evaporative self-assembly. Angew. Chem. Int. Ed. 2012, 51, 1534–1546. [Google Scholar] [CrossRef] [PubMed]
  7. Rabani, E.; Reichman, D.R.; Geissler, P.L.; Brus, L.E. Drying-mediated self-assembly of nanoparticles. Nature 2003, 426, 271–274. [Google Scholar] [CrossRef] [PubMed]
  8. Zhang, X.; Douglas, J.F.; Jones, R.L. Influence of film casting method on block copolymer ordering in thin films. Soft Matter 2012, 8, 4980–4987. [Google Scholar] [CrossRef]
  9. Strawhecker, K.E.; Kumar, S.K.; Douglas, J.F.; Karim, A. The critical role of solvent evaporation on the roughness of spin-cast polymer films. Macromolecules 2001, 34, 4669–4672. [Google Scholar] [CrossRef]
  10. De Gennes, P.-G. Solvent evaporation of spin cast films: Crust effects. Eur. Phys. J. E 2001, 7, 31–34. [Google Scholar] [CrossRef]
  11. Jouault, N.; Zhao, D.; Kumar, S.K. Role of casting solvent on nanoparticle dispersion in polymer nanocomposites. Macromolecules 2014, 47, 5246–5255. [Google Scholar] [CrossRef]
  12. Español, P. Hydrodynamics from dissipative particle dynamics. Phys. Rev. E 1995, 52, 1734–1742. [Google Scholar] [CrossRef]
  13. Pivkin, I.V.; Caswell, B.; Karniadakisa, G.E. Dissipative particle dynamics. In Reviews in Computational Chemistry; John Wiley & Sons, Inc.: New York, NY, USA, 2010. [Google Scholar]
  14. Hoogerbrugge, P.J.; Koelman, J.M.V.A. Simulating microscopic hydrodynamic phenomena with dissipative particle dynamics. Europhys. Lett. 1992, 19, 155–160. [Google Scholar] [CrossRef]
  15. Groot, R.D.; Warren, P.B. Dissipative particle dynamics: Bridging the gap between atomistic and mesoscopic simulation. J. Chem. Phys. 1997, 107, 4423–4435. [Google Scholar] [CrossRef]
  16. Groot, R.D.; Madden, T.J. Dynamic simulation of diblock copolymer microphase separation. J. Chem. Phys. 1998, 108, 8713–8724. [Google Scholar] [CrossRef]
  17. Groot, R.D.; Madden, T.J.; Tildesley, D.J. On the role of hydrodynamic interactions in block copolymer microphase separation. J. Chem. Phys. 1999, 110, 9739. [Google Scholar] [CrossRef]
  18. Symeonidis, V.; Karniadakis, G.E.; Caswell, B. Dissipative particle dynamics simulations of polymer chains: Scaling laws and shearing response compared to DNA experiments. Phys. Rev. Lett. 2005, 95. [Google Scholar] [CrossRef]
  19. Li, Z.; Drazer, G. Hydrodynamic interactions in dissipative particle dynamics. Phys. Fluids 2008, 20, 103601. [Google Scholar] [CrossRef]
  20. Schlijper, A.G.; Hoogerbrugge, P.J.; Manke, C.W. Computer simulation of dilute polymer solutions with the dissipative particle dynamics method. J. Rheol. 1995, 39, 567–579. [Google Scholar] [CrossRef]
  21. Spenley, N.A. Scaling laws for polymers in dissipative particle dynamics. Europhys. Lett. 2000, 49, 534–540. [Google Scholar] [CrossRef]
  22. Kong, Y.; Manke, C.W.; Madden, W.G.; Schlijper, A.G. Effect of solvent quality on the conformation and relaxation of polymers via dissipative particle dynamics. J. Chem. Phys. 1997, 107, 592–602. [Google Scholar] [CrossRef]
  23. Jiang, W.; Huang, J.; Wang, Y.; Laradji, M. Hydrodynamic interaction in polymer solutions simulated with dissipative particle dynamics. J. Chem. Phys. 2007, 126, 044901. [Google Scholar] [CrossRef] [PubMed]
  24. Zhao, T.; Wang, X.; Jiang, L.; Larson, R.G. Dissipative particle dynamics simulation of dilute polymer solutions—Inertial effects and hydrodynamic interactions. J. Rheol. 2014, 58, 1039–1058. [Google Scholar] [CrossRef]
  25. Malfreyt, P.; Tildesley, D.J. Dissipative particle dynamics simulations of grafted polymer chains between two walls. Langmuir 2000, 16, 4732–4740. [Google Scholar] [CrossRef]
  26. Raos, G.; Moreno, M.; Elli, S. Computational experiments on filled rubber viscoelasticity: What is the role of particle-particle interactions? Macromolecules 2006, 39, 6744–6751. [Google Scholar] [CrossRef]
  27. Yong, X.; Kuksenok, O.; Matyjaszewski, K.; Balazs, A.C. Harnessing interfacially-active nanorods to regenerate severed polymer gels. Nano Lett. 2013, 13, 6269–6274. [Google Scholar] [CrossRef] [PubMed]
  28. Yong, X. Modeling the assembly of polymer-grafted nanoparticles at oil-water interfaces. Langmuir 2015, 31, 11458–11469. [Google Scholar] [CrossRef] [PubMed]
  29. Yong, X.; Simakova, A.; Averick, S.; Gutierrez, J.; Kuksenok, O.; Balazs, A.C.; Matyjaszewski, K. Stackable, covalently fused gels: Repair and composite formation. Macromolecules 2015, 48, 1169–1178. [Google Scholar] [CrossRef]
  30. Pagonabarraga, I.; Frenkel, D. Dissipative particle dynamics for interacting systems. J. Chem. Phys. 2001, 115, 5015–5026. [Google Scholar] [CrossRef]
  31. Trofimov, S.Y.; Nies, E.L.F.; Michels, M.A.J. Constant-pressure simulations with dissipative particle dynamics. J. Chem. Phys. 2005, 123, 144102. [Google Scholar] [CrossRef] [PubMed]
  32. Warren, P.B. Vapor-liquid coexistence in many-body dissipative particle dynamics. Phys. Rev. E 2003, 68, 066702. [Google Scholar] [CrossRef] [PubMed]
  33. Tiwari, A.; Abraham, J. Dissipative-particle-dynamics model for two-phase flows. Phys. Rev. E 2006, 74, 056701. [Google Scholar] [CrossRef] [PubMed]
  34. Arienti, M.; Pan, W.; Li, X.; Karniadakis, G. Many-body dissipative particle dynamics simulation of liquid/vapor and liquid/solid interactions. J. Chem. Phys. 2011, 134, 204114. [Google Scholar] [CrossRef] [PubMed]
  35. Yong, X.; Qin, S.; Singler, T.J. Nanoparticle-mediated evaporation at liquid–vapor interfaces. Extreme Mech. Lett. 2016, 7, 90–103. [Google Scholar] [CrossRef]
  36. Ghoufi, A.; Malfreyt, P. Mesoscale modeling of the water liquid-vapor interface: A surface tension calculation. Phys. Rev. E 2011, 83, 051601. [Google Scholar] [CrossRef] [PubMed]
  37. Ghoufi, A.; Malfreyt, P. Coarse grained simulations of the electrolytes at the water-air interface from many body dissipative particle dynamics. J. Chem. Theory Comput. 2012, 8, 787–791. [Google Scholar] [CrossRef] [PubMed]
  38. Ghoufi, A.; Emile, J.; Malfreyt, P. Recent advances in many body dissipative particles dynamics simulations of liquid-vapor interfaces. Eur. Phys. J. E 2013, 36, 1–12. [Google Scholar] [CrossRef] [PubMed]
  39. Chen, C.; Zhuang, L.; Li, X.; Dong, J.; Lu, J. A Many-body dissipative particle dynamics study of forced water-oil displacement in capillary. Langmuir 2012, 28, 1330–1336. [Google Scholar] [CrossRef] [PubMed]
  40. Wang, Y.; Chen, S. Numerical study on droplet sliding across micropillars. Langmuir 2015, 31, 4673–4677. [Google Scholar] [CrossRef] [PubMed]
  41. Jamali, S.; Boromand, A.; Khani, S.; Wagner, J.; Yamanoi, M.; Maia, J. Generalized mapping of multi-body dissipative particle dynamics onto fluid compressibility and the Flory-Huggins Theory. J. Chem. Phys. 2015, 142, 164902. [Google Scholar] [CrossRef] [PubMed]
  42. Trofimov, S.Y.; Nies, E.L.F.; Michels, M.A.J. Thermodynamic consistency in dissipative particle dynamics simulations of strongly nonideal liquids and liquid mixtures. J. Chem. Phys. 2002, 117, 9383–9394. [Google Scholar] [CrossRef]
  43. Zhu, Y.-L.; Liu, H.; Lu, Z.-Y. A highly coarse-grained model to simulate entangled polymer melts. J. Chem. Phys. 2012, 136, 144903. [Google Scholar] [CrossRef] [PubMed]
  44. Li, Y.C.; Liu, H.; Huang, X.R.; Sun, C.C. Evaporation- and surface-induced morphology of symmetric diblock copolymer thin films: A multibody dissipative particle dynamics study. Mol. Simul. 2011, 37, 875–883. [Google Scholar] [CrossRef]
  45. Tiwari, A.; Reddy, H.; Mukhopadhyay, S.; Abraham, J. Simulations of liquid nanocylinder breakup with dissipative particle dynamics. Phys. Rev. E 2008, 78, 1–11. [Google Scholar] [CrossRef]
  46. Stroberg, W.; Keten, S.; Liu, W.K. Hydrodynamics of capillary imbibition under nanoconfinement. Langmuir 2012, 28, 14488–14495. [Google Scholar] [CrossRef] [PubMed]
  47. Cupelli, C.; Henrich, B.; Glatzel, T.; Zengerle, R.; Moseler, M.; Santer, M. Dynamic capillary wetting studied with dissipative particle dynamics. New J. Phys. 2008, 10, 1–4. [Google Scholar] [CrossRef]
  48. Li, Z.; Hu, G.-H.; Wang, Z.-L.; Ma, Y.-B.; Zhou, Z.-W. Three dimensional flow structures in a moving droplet on substrate: A dissipative particle dynamics study. Phys. Fluids 2013, 25, 072103. [Google Scholar] [CrossRef]
  49. Chen, C.; Gao, C.; Zhuang, L.; Li, X.; Wu, P.; Dong, J.; Lu, J. A many-body dissipative particle dynamics study of spontaneous capillary imbibition and drainage. Langmuir 2010, 26, 9533–9538. [Google Scholar] [CrossRef] [PubMed]
  50. Español, P.; Warren, P. Statistical mechanics of dissipative particle dynamics. Europhys. Lett. 1995, 30, 191–196. [Google Scholar] [CrossRef]
  51. Sirk, T.W.; Slizoberg, Y.R.; Brennan, J.K.; Lisal, M.; Andzelm, J.W. An enhanced entangled polymer model for dissipative particle dynamics. J. Chem. Phys. 2012, 136, 134903. [Google Scholar] [CrossRef] [PubMed]
  52. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef]
  53. Warren, P.B. No-go theorem in many-body dissipative particle dynamics. Phys. Rev. E 2013, 87, 045303. [Google Scholar] [CrossRef] [PubMed]
  54. Binder, K. Monte Carlo and Molecular Dynamics Simulations in Polymer Science; Oxford University Press: New York, NY, USA, 1995. [Google Scholar]
  55. Hossain, D.; Tschopp, M.A.; Ward, D.K.; Bouvard, J.L.; Wang, P.; Horstemeyer, M.F. Molecular dynamics simulations of deformation mechanisms of amorphous polyethylene. Polymer 2010, 51, 6071–6083. [Google Scholar] [CrossRef]
  56. Groot, R.D.; Rabone, K.L. Mesoscopic simulation of cell membrane damage, morphology change and rupture by nonionic surfactants. Biophys. J. 2001, 81, 725–736. [Google Scholar] [CrossRef]
  57. Flory, P.J. Principles of Polymer Chemistry, 1st ed.; Cornell University Press: Ithaca, NY, USA, 1953. [Google Scholar]
  58. Maiti, A.; McGrother, S. Bead–bead interaction parameters in dissipative particle dynamics: Relation to bead-size, solubility parameter, and surface tension. J. Chem. Phys. 2004, 120, 1594. [Google Scholar] [CrossRef] [PubMed]
  59. Travis, K.P.; Bankhead, M.; Good, K.; Owens, S.L. New parametrization method for dissipative particle dynamics. J. Chem. Phys. 2007, 127, 014109. [Google Scholar] [CrossRef] [PubMed]
  60. Liyana-Arachchi, T.P.; Jamadagni, S.N.; Eike, D.; Koenig, P.H.; Siepmann, J.I. Liquid-liquid equilibria for soft-repulsive particles: Improved equation of state and methodology for representing molecules of different sizes and chemistry in dissipative particle dynamics. J. Chem. Phys. 2015, 142, 044902. [Google Scholar] [CrossRef] [PubMed]
  61. Irving, J.H.; Kirkwood, J.G. The statistical mechanical theory of transport processes. IV. The equations of hydrodynamics. J. Chem. Phys. 1950, 18, 817. [Google Scholar] [CrossRef]
  62. De Gennes, P.G. Collapse of a polymer chain in poor solvents. J. Phys. Lett. 1975, 36, 55–57. [Google Scholar] [CrossRef]
  63. De Gennes, P.G. Collapse of a flexible polymer chain II. J. Phys. Lett. 1978, 39, 299–301. [Google Scholar] [CrossRef]
  64. Zimm, B.H. Dynamics of polymer molecules in dilute solution: Viscoelasticity, flow birefringence and dielectric loss. J. Chem. Phys. 1956, 24, 269–278. [Google Scholar] [CrossRef]
  65. Rouse, P.E. A theory of the linear viscoelastic properties of dilute solutions of coiling polymers. J. Chem. Phys. 1953, 21, 1272–1280. [Google Scholar] [CrossRef]
  66. Pierotti, R.A. A scaled particle theory of aqueous and nonaqueous solutions. Chem. Rev. 1976, 76, 717–726. [Google Scholar] [CrossRef]
  67. Simmons, D.S.; Sanchez, I.C. Scaled particle theory for the coil–globule transition of an isolated polymer chain. Macromolecules 2013, 46, 4691–4697. [Google Scholar] [CrossRef]
  68. Pan, G.; Manke, C.W. Developments toward simulation of entangled polymer melts by dissipative particle dynamics (DPD). Int. J. Mod. Phys. B 2003, 17, 231–235. [Google Scholar] [CrossRef]
  69. Kaznessis, Y.N.; Hill, D.A.; Maginn, E.J. A molecular dynamics study of macromolecules in good solvents: Comparison with dielectric spectroscopy experiments. J. Chem. Phys. 1998, 109, 5078–5088. [Google Scholar] [CrossRef]
  70. Kumar, S.; Larson, R.G. Brownian dynamics simulations of flexible polymers with spring–spring repulsions. J. Chem. Phys. 2001, 114, 6937–6941. [Google Scholar] [CrossRef]
  71. Padding, J.T.; Briels, W.J. Uncrossability constraints in mesoscopic polymer melt simulations: Non-rouse behavior of C120H242. J. Chem. Phys. 2001, 115, 2846–2859. [Google Scholar] [CrossRef]
  72. Hoda, N.; Larson, R.G. Brownian dynamics simulations of single polymer chains with and without self-entanglements in theta and good solvents under imposed flow fields. J. Rheol. 2010, 54, 1061–1081. [Google Scholar] [CrossRef]
  73. Nikunen, P.; Vattulainen, I.; Karttunen, M. Reptational dynamics in dissipative particle dynamics simulations of polymer melts. Phys. Rev. E 2007, 75, 036713. [Google Scholar] [CrossRef] [PubMed]
  74. Sliozberg, Y.R.; Sirk, T.W.; Brennan, J.K.; Andzelm, J.W. Bead-spring models of entangled polymer melts: Comparison of hard-core and soft-core potentials. J. Polym. Sci. Part B 2012, 50, 1694–1698. [Google Scholar] [CrossRef]
  75. Goujon, F.; Malfreyt, P.; Tildesley, D.J. Mesoscopic simulation of entanglements using dissipative particle dynamics: Application to polymer brushes. J. Chem. Phys. 2008, 129, 034902. [Google Scholar] [CrossRef] [PubMed]
  76. Holleran, S.P.; Larson, R.G. Using spring repulsions to model entanglement interactions in brownian dynamics simulations of bead-spring chains. Rheol. Acta 2008, 47, 3–17. [Google Scholar] [CrossRef]
Figure 1. (a) Relation between effective χ parameter and excess repulsion between unlike species in a binary mixture of monomers defined as Δ A = A ab A a a , where a and b represent the types of interacting beads. The dash line is the least-square linear fitting. (b) Simulated surface tension of the binary mixture as a function of χ parameter. The smooth curve is a fit to a 3/2 power law with a non-classical critical point.
Figure 1. (a) Relation between effective χ parameter and excess repulsion between unlike species in a binary mixture of monomers defined as Δ A = A ab A a a , where a and b represent the types of interacting beads. The dash line is the least-square linear fitting. (b) Simulated surface tension of the binary mixture as a function of χ parameter. The smooth curve is a fit to a 3/2 power law with a non-classical critical point.
Polymers 08 00426 g001
Figure 2. (a) Mean-square radius of gyration as a function of solvent quality quantified by A ps for dilute polymer solutions with different chain lengths. Snapshots of dilute polymer solutions with chain length N = 100 in the: (b) globule; and (c) expanded coil states. Colors represent different chains. The lines are only guides to the eyes. The error bars represent the standard deviations of the ensemble averages.
Figure 2. (a) Mean-square radius of gyration as a function of solvent quality quantified by A ps for dilute polymer solutions with different chain lengths. Snapshots of dilute polymer solutions with chain length N = 100 in the: (b) globule; and (c) expanded coil states. Colors represent different chains. The lines are only guides to the eyes. The error bars represent the standard deviations of the ensemble averages.
Polymers 08 00426 g002
Figure 3. Log–log plot of mean-square radius of gyration as a function of the number of Hookean springs in each chain for different solvent quality. The dash-dot lines correspond to the power-law fittings of simulation data. The error bars are smaller than the symbols.
Figure 3. Log–log plot of mean-square radius of gyration as a function of the number of Hookean springs in each chain for different solvent quality. The dash-dot lines correspond to the power-law fittings of simulation data. The error bars are smaller than the symbols.
Polymers 08 00426 g003
Figure 4. Time-averaged coordination numbers for polymer and solvent beads in the athermal solution for A ps = 40 . The inset plots the time-averaged polymer–polymer and polymer–solvent radial distribution functions. The lines are only guide to the eyes.
Figure 4. Time-averaged coordination numbers for polymer and solvent beads in the athermal solution for A ps = 40 . The inset plots the time-averaged polymer–polymer and polymer–solvent radial distribution functions. The lines are only guide to the eyes.
Polymers 08 00426 g004
Figure 5. Characteristic ratio of the modeled chains as a function of chain length for different solvent quality. Poor, theta, and good solvents refer to the systems with polymer–solvent attraction parameter A ps = −36, −38.5, and −40, respectively. The lines are only guides to the eyes. The error bars represent the standard deviations of the ensemble averages.
Figure 5. Characteristic ratio of the modeled chains as a function of chain length for different solvent quality. Poor, theta, and good solvents refer to the systems with polymer–solvent attraction parameter A ps = −36, −38.5, and −40, respectively. The lines are only guides to the eyes. The error bars represent the standard deviations of the ensemble averages.
Polymers 08 00426 g005
Figure 6. Log–log plots of: (a) the diffusion coefficient; and (b) the longest relaxation time as functions of chain length in the theta solvent. The black dash-dot line in (a) indicates the diffusion coefficient of the solvent beads. The dashed lines in (a,b) represent the power-law fittings of the data. The error bars are smaller than the symbols.
Figure 6. Log–log plots of: (a) the diffusion coefficient; and (b) the longest relaxation time as functions of chain length in the theta solvent. The black dash-dot line in (a) indicates the diffusion coefficient of the solvent beads. The dashed lines in (a,b) represent the power-law fittings of the data. The error bars are smaller than the symbols.
Polymers 08 00426 g006
Figure 7. (a) Dependence of the normalized mean-square radius gyration (with respect to the values in the dilute limit R g 0 2 ) of the N = 30 polymer chains on their volume fraction for good ( A ps = 40 ) and theta ( A ps = 38.5 ) solvent conditions; (b) The diffusion coefficient and relaxation time of the chains in the theta solvent as functions of the volume fraction. The error bars represent the standard deviations of the ensemble averages, and their sizes are smaller than the symbols in (b).
Figure 7. (a) Dependence of the normalized mean-square radius gyration (with respect to the values in the dilute limit R g 0 2 ) of the N = 30 polymer chains on their volume fraction for good ( A ps = 40 ) and theta ( A ps = 38.5 ) solvent conditions; (b) The diffusion coefficient and relaxation time of the chains in the theta solvent as functions of the volume fraction. The error bars represent the standard deviations of the ensemble averages, and their sizes are smaller than the symbols in (b).
Polymers 08 00426 g007
Figure 8. Schematic of a topology violation, where the minimum distance vector between bonds changes orientation by an angle larger than 90° in one time step. Red and green colors represent different bonds.
Figure 8. Schematic of a topology violation, where the minimum distance vector between bonds changes orientation by an angle larger than 90° in one time step. Red and green colors represent different bonds.
Polymers 08 00426 g008
Figure 9. Relative number of topology violations as a function of polymer volume fraction in the good solvent with A ps = 40 . The reference number of topology violation is for the dilute limit. The inset shows the distribution of the α angle between 170° and 180° for the topology violation events in the concentrated solution ( ϕ = 0.7 ) during a total simulation time of t = 5000 τ .
Figure 9. Relative number of topology violations as a function of polymer volume fraction in the good solvent with A ps = 40 . The reference number of topology violation is for the dilute limit. The inset shows the distribution of the α angle between 170° and 180° for the topology violation events in the concentrated solution ( ϕ = 0.7 ) during a total simulation time of t = 5000 τ .
Polymers 08 00426 g009
Table 1. Structural properties of dilute polymer solutions in the theta solvent condition. The mean-square radius of gyration R g 2 and mean-square end-to-end distance R 1 N 2 were obtained by averaging over all chains in the system for each frame, and then averaged over 5000 frames with a separation of 1τ between consecutive frames.
Table 1. Structural properties of dilute polymer solutions in the theta solvent condition. The mean-square radius of gyration R g 2 and mean-square end-to-end distance R 1 N 2 were obtained by averaging over all chains in the system for each frame, and then averaged over 5000 frames with a separation of 1τ between consecutive frames.
N R g 2 R 1 N 2 R 1 N 2 / R g 2 C n
50.432.225.231.14
80.744.125.571.21
100.955.415.681.23
202.0612.065.851.30
303.2018.895.901.34
404.3625.875.941.36
505.3931.715.891.32
808.7452.646.031.37
10011.2166.915.971.38
14015.8292.875.81.37
Table 2. Average number of topology violations per 100τ for in semidilute solutions with different polymer volume fractions. The results are obtained in simulations having a total time of t = 5000 τ .
Table 2. Average number of topology violations per 100τ for in semidilute solutions with different polymer volume fractions. The results are obtained in simulations having a total time of t = 5000 τ .
φ0.020.050.10.20.30.40.50.60.7
N TV 169.7168.7191.3203.6237.4274.3310.9376.0432.8

Share and Cite

MDPI and ACS Style

Yong, X. Hydrodynamic Interactions and Entanglements of Polymer Solutions in Many-Body Dissipative Particle Dynamics. Polymers 2016, 8, 426. https://doi.org/10.3390/polym8120426

AMA Style

Yong X. Hydrodynamic Interactions and Entanglements of Polymer Solutions in Many-Body Dissipative Particle Dynamics. Polymers. 2016; 8(12):426. https://doi.org/10.3390/polym8120426

Chicago/Turabian Style

Yong, Xin. 2016. "Hydrodynamic Interactions and Entanglements of Polymer Solutions in Many-Body Dissipative Particle Dynamics" Polymers 8, no. 12: 426. https://doi.org/10.3390/polym8120426

APA Style

Yong, X. (2016). Hydrodynamic Interactions and Entanglements of Polymer Solutions in Many-Body Dissipative Particle Dynamics. Polymers, 8(12), 426. https://doi.org/10.3390/polym8120426

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop