Chains Of Tangent Hard Spheres Research Articles

Entropy-Driven Crystallization of Hard Colloidal Mixtures of Polymers and Monomers.

The most trivial example of self-assembly is the entropy-driven crystallization of hard spheres. Past works have established the similarities and differences in the phase behavior of monomers and chains made of hard spheres. Inspired by the difference in the melting points of the pure components, we study, through Monte Carlo simulations, the phase behavior of athermal mixtures composed of fully flexible polymers and individual monomers of uniform size. We analyze how the relative number fraction and the packing density affect crystallization and the established ordered morphologies. As a first result, a more precise determination of the melting point for freely jointed chains of tangent hard spheres is extracted. A synergetic effect is observed in the crystallization leading to synchronous crystallization of the two species. Structural analysis of the resulting ordered morphologies shows perfect mixing and thus no phase separation. Due to the constraints imposed by chain connectivity, the local environment of the individual spheres, as quantified by the Voronoi polyhedron, is systematically more spherical and more symmetric compared to that of spheres belonging to chains. In turn, the local environment of the ordered phase is more symmetric and more spherical compared to that of the initial random packing, demonstrating the entropic origins of the phase transition. In general, increasing the polymer content reduces the degree of crystallinity and increases the melting point to higher volume fractions. According to the present findings, relative concentration is another determining factor in controlling the phase behavior of hard colloidal mixtures based on polymers.

Densest packing of flexible polymers in 2D films.

How dense objects, particles, atoms, and molecules can be packed is intimately related to the properties of the corresponding hosts and macrosystems. We present results from extensive Monte Carlo simulations on maximally compressed packings of linear, freely jointed chains of tangent hard spheres of uniform size in films whose thickness is equal to the monomer diameter. We demonstrate that fully flexible chains of hard spheres can be packed as efficiently as monomeric analogs, within a statistical tolerance of less than 1%. The resulting ordered polymer morphology corresponds to an almost perfect hexagonal triangular (TRI) crystal of the p6m wallpaper group, whose sites are occupied by the chain monomers. The Flory scaling exponent, which corresponds to the maximally dense polymer packing in 2D, has a value of ν = 0.62, which lies between the limits of 0.50 (compact and collapsed state) and 0.75 (self-avoiding random walk).

Crystallization of Flexible Chains of Tangent Hard Spheres under Full Confinement.

We present results from extensive Monte Carlo simulations on the crystallization of athermal polymers under full confinement. Polymers are represented as freely jointed chains of tangent hard spheres of uniform size. Confinement is applied through the presence of flat, parallel, and impenetrable walls in all dimensions. We analyze crystallization as the summation of two contributions: one that occurs in the bulk volume of the system (bulk crystallization), and one on the wall surfaces (surface crystallization). Depending on volume fraction initially amorphous (disordered) hard-sphere chain packings transit to the stable crystal phase. The established ordered morphologies consist primarily of hexagonal close-packed (HCP) crystals in the bulk volume and of triangular (TRI) crystals on the surface. As in the case of athermal packings in the bulk (without confinement), a structural competition is observed between the 5-fold local symmetry and the formation of close-packed crystallites. Effectively, the full confinement inside a cube favors the growth of the HCP crystal, as the FCC one is quite incompatible with the imposed spatial constraints. Consequently, we observe the formation of noncompact ordered motifs which grow from the surface to the inner volume of the simulation cell. We further compare the 2D and 3D crystals formed by monomeric hard spheres under the same simulation conditions. Significant differences are observed at low densities that tend to diminish as concentration increases.

Entropy-Driven Heterogeneous Crystallization of Hard-Sphere Chains under Unidimensional Confinement.

We investigate, through Monte Carlo simulations, the heterogeneous crystallization of linear chains of tangent hard spheres under confinement in one dimension. Confinement is realized through flat, impenetrable, and parallel walls. A wide range of systems is studied with respect to their average chain lengths (N = 12 to 100) and packing densities (φ = 0.50 to 0.61). The local structure is quantified through the Characteristic Crystallographic Element (CCE) norm descriptor. Here, we split the phenomenon into the bulk crystallization, far from the walls, and the projected surface crystallization in layers adjacent to the confining surfaces. Once a critical volume fraction is met, the chains show a phase transition, starting from regions near the hard walls. The established crystal morphologies consist of alternating hexagonal close-packed or face-centered cubic layers with a stacking direction perpendicular to the confining walls. Crystal layer perfection is observed with an increasing concentration. As in the case of the unconstrained phase transition of athermal polymers at high densities, crystal nucleation and growth compete with the formation of sites of a fivefold local symmetry. While surface crystallites show perfection with a predominantly triangular character, the morphologies of square crystals or of a mixed type are also formed. The simulation results show that the rate of perfection of the surface crystallization is not significantly faster than that of the bulk crystallization.

Phase Behavior of Lyotropic Liquid-Crystalline Polymers under Varying Solvent Conditions: Effect of External Fields on the Phase Diagram.

In this work, we report a Density Functional Theory based study of phase behavior of lyotropic liquid-crystalline polymers under both good and varying solvent conditions in the presence of external electric or magnetic field. Our microscopic model for the good solvent case is based on the tangent hard-sphere chain with bond-bending potential to account for the chain stiffness; the variable solvent quality is modeled by adding attractive monomer-monomer interactions. The phase diagrams are constructed in three intensive variables (temperature, pressure, and field strength), and are characterized by the presence of critical and triple lines, which originate from the critical and triple points of the corresponding zero-field case. The merging of critical and triple lines results in the appearance of the "double critical" and "critical triple" points, already known from the earlier studies of the phase behavior of spin fluids in magnetic fields. The important difference of the present model from the spin fluids is due to the finite stiffness of the polymer chains (characterized by their persistence length), which adds an additional parameter controlling the morphology of the phase diagrams.

Crystal, Fivefold and Glass Formation in Clusters of Polymers Interacting with the Square Well Potential.

We present results, from Monte Carlo (MC) simulations, on polymer systems of freely jointed chains with spherical monomers interacting through the square well potential. Starting from athermal packings of chains of tangent hard spheres, we activate the square well potential under constant volume and temperature corresponding effectively to instantaneous quenching. We investigate how the intensity and range of pair-wise interactions affected the final morphologies by fixing polymer characteristics such as average chain length and tolerance in bond gaps. Due to attraction chains are brought closer together and they form clusters with distinct morphologies. A wide variety of structures is obtained as the model parameters are systematically varied: weak interactions lead to purely amorphous clusters followed by well-ordered ones. The latter include the whole spectrum of crystal morphologies: from virtually perfect hexagonal close packed (HCP) and face centered cubic (FCC) crystals, to random hexagonal close packed layers of single stacking direction of alternating HCP and FCC layers, to structures of mixed HCP/FCC character with multiple stacking directions and defects in the form of twins. Once critical values of interaction are met, fivefold-rich glassy clusters are formed. We discuss the similarities and differences between energy-driven crystal nucleation in thermal polymer systems as opposed to entropy-driven phase transition in athermal polymer packings. We further calculate the local density of each site, its dependence on the distance from the center of the cluster and its correlation with the crystallographic characteristics of the local environment. The short- and long-range conformations of chains are analyzed as a function of the established cluster morphologies.

Self-Avoiding Random Walks as a Model to Study Athermal Linear Polymers under Extreme Plate Confinement.

Monte Carlo (MC) simulations, built around chain-connectivity-altering moves and a wall-displacement algorithm, allow us to simulate freely-jointed chains of tangent hard spheres of uniform size under extreme confinement. The latter is realized through the presence of two impenetrable, flat, and parallel plates. Extreme conditions correspond to the case where the distance between the plates approaches the monomer size. An analysis of the local structure, based on the characteristic crystallographic element (CCE) norm, detects crystal nucleation and growth at packing densities well below the ones observed in bulk analogs. In a second step, we map the confined polymer chains into self-avoiding random walks (SAWs) on restricted lattices. We study all realizations of the cubic crystal system: simple, body centered, and face centered cubic crystals. For a given chain size (SAW length), lattice type, origin of SAW, and level of confinement, we enumerate all possible SAWs (equivalently all chain conformations) and calculate the size distribution. Results for intermediate SAW lengths are used to predict the behavior of long, fully entangled chains through growth formulas. The SAW analysis will allow us to determine the corresponding configurational entropy, as it is the driving force for the observed phase transition and the determining factor for the thermodynamic stability of the corresponding crystal morphologies.

The isotropic-nematic and nematic-nematic phase transition of binary mixtures of tangent hard-sphere chain fluids: an analytical equation of state.

An analytical equation of state (EoS) is derived to describe the isotropic (I) and nematic (N) phase of linear- and partially flexible tangent hard-sphere chain fluids and their mixtures. The EoS is based on an extension of Onsager's second virial theory that was developed in our previous work [T. van Westen, B. Oyarzún, T. J. H. Vlugt, and J. Gross, J. Chem. Phys. 139, 034505 (2013)]. Higher virial coefficients are calculated using a Vega-Lago rescaling procedure, which is hereby generalized to mixtures. The EoS is used to study (1) the effect of length bidispersity on the I-N and N-N phase behavior of binary linear tangent hard-sphere chain fluid mixtures, (2) the effect of partial molecular flexibility on the binary phase diagram, and (3) the solubility of hard-sphere solutes in I- and N tangent hard-sphere chain fluids. By changing the length bidispersity, two types of phase diagrams were found. The first type is characterized by an I-N region at low pressure and a N-N demixed region at higher pressure that starts from an I-N-N triphase equilibrium. The second type does not show the I-N-N equilibrium. Instead, the N-N region starts from a lower critical point at a pressure above the I-N region. The results for the I-N region are in excellent agreement with the results from molecular simulations. It is shown that the N-N demixing is driven both by orientational and configurational/excluded volume entropy. By making the chains partially flexible, it is shown that the driving force resulting from the configurational entropy is reduced (due to a less anisotropic pair-excluded volume), resulting in a shift of the N-N demixed region to higher pressure. Compared to linear chains, no topological differences in the phase diagram were found. We show that the solubility of hard-sphere solutes decreases across the I-N phase transition. Furthermore, it is shown that by using a liquid crystal mixture as the solvent, the solubility difference can by maximized by tuning the composition. Theoretical results for the Henry's law constant of the hard-sphere solute are in good agreement with the results from molecular simulation.

The isotropic-nematic phase transition of tangent hard-sphere chain fluids—Pure components

An extension of Onsager's second virial theory is developed to describe the isotropic-nematic phase transition of tangent hard-sphere chain fluids. Flexibility is introduced by the rod-coil model. The effect of chain-flexibility on the second virial coefficient is described using an accurate, analytical approximation for the orientation-dependent pair-excluded volume. The use of this approximation allows for an analytical treatment of intramolecular flexibility by using a single pure-component parameter. Two approaches to approximate the effect of the higher virial coefficients are considered, i.e., the Vega-Lago rescaling and Scaled Particle Theory (SPT). The Onsager trial function is employed to describe the orientational distribution function. Theoretical predictions for the equation of state and orientational order parameter are tested against the results from Monte Carlo (MC) simulations. For linear chains of length 9 and longer, theoretical results are in excellent agreement with MC data. For smaller chain lengths, small errors introduced by the approximation of the higher virial coefficients become apparent, leading to a small under- and overestimation of the pressure and density difference at the phase transition, respectively. For rod-coil fluids of reasonable rigidity, a quantitative comparison between theory and MC simulations is obtained. For more flexible chains, however, both the Vega-Lago rescaling and SPT lead to a small underestimation of the location of the phase transition.

Liquid-vapor equilibrium and surface properties of short rigid chains with one long range attractive potential

Liquid-vapor coexistence and interfacial properties of short lineal rigid vibrating chains with three tangent monomers in two and three dimensions are calculated. The effect of the range and position of a long ranged square well attractive potential is studied. Orthobaric densities, vapor pressures, surface tensions, and interfacial widths are reported. Two types of molecules are studied. Chains of three tangent hard sphere monomers and chains of three and five tangent hard sphere monomers interacting with a square well attractive potential with λ(∗) = λ∕σ = 1.5 in units of the hard core diameter σ. The results are reported in two and three dimensions. For both types of chains, a long ranged square well attractive potential is located at various positions in the chain to investigate its effect in the properties of the corresponding systems. Results for hard sphere chains are presented for a series of different sizes of λ(∗) between 2.5 and 5. For square well chains the position in the chain of the long ranged potential has no influence in the coexistence and interfacial properties. Critical temperatures increase monotonically with respect to λ(∗) and critical densities decrease systematically for both types of chains. When the long ranged potential is located in the middle monomer of the hard sphere chains no critical point is found for λ(∗) < 2.4. No critical point is found when the long ranged potential is located in one of the extremes of the hard sphere chains.

An equation of state for the isotropic phase of linear, partially flexible and fully flexible tangent hard-sphere chain fluids

A new equation of state is developed that accurately describes the isotropic phase behaviour of linear, partially flexible and fully flexible tangent hard-sphere chain fluids and their mixtures. The equation of state is based on the equation of state of Liu and Hu [H. Liu and Y. Hu, Fluid Phase Equilibr. 122, 75 (1996)] for fully flexible chain fluids. The effect of molecular flexibility is described by a pure-component parameter that is introduced in the theory at the level of the cavity correlation function of next-to-nearest neighbour segments in a chain molecule. The equation of state contains a total of three adjustable model constants. The extension to partially flexible- and linear chain fluids is based on a refitting of the first model constant to numerical data of the second virial coefficient of partially flexible and linear tangent hard-sphere chain fluids. The numerical data were obtained from an analytical approximation for the pair-excluded volume. The other two parameters were adjusted to molecular simulation data for the pressure of linear tangent hard-sphere chain fluids. For both, pure component systems and mixtures of chains of variable flexibility, the pressure and second virial coefficient obtained from the equation of state, are in excellent agreement with the results from Monte Carlo simulations. A significant improvement to TPT1, TPT2, generalised Flory-dimer theory and scaled particle theory is observed.

Jamming and crystallization in athermal polymer packings

Dense packings of chains of hard spheres possess characteristic features that do not have a counterpart in corresponding packings of monomeric spheres especially near the maximally random jammed (MRJ) state. From the modelling perspective the additional requirement that spheres keep their connectivity while maximizing the occupied volume fraction imposes severe constraints on generation algorithms of dense chain configurations. The extremely sluggish dynamics imposed by the uncrossability of chains precludes the use of deterministic or stochastic dynamics to generate all but dilute polymer packings. As a viable alternative, especially tailored chain-connectivity-altering Monte Carlo (MC) algorithms have been developed that bypass this kinetic hindrance and have actually been able to produce packings of hard-sphere chains in a volume fraction range spanning from infinite dilution up to the MRJ state. Such very dense athermal polymer packings share a number of structural features with packings of monomeric hard spheres, but also display unique characteristics due to the constraints imposed by connectivity. We give an overview of the most relevant results of our recent modeling work on packings of freely-jointed chains of tangent hard spheres about the MRJ state, local structure, chain dimensions and their scaling with density, topological constraints in the form of entanglements and knots, contact network at jamming, and entropically driven crystallization.

The equation of state of flexible chains of tangent hard spheres at high-density region from simulation and thermodynamic perturbation theory

Radial and triplet correlation functions of the reference hard sphere system are determined at several solid densities by canonical Monte Carlo (MC) simulations. These customized data are used to extend the second order thermodynamic perturbation theory (TPT) to the solid phase of flexible hard chain systems. In order to test the accuracy of the TPT equation of state (EOS) for hard chains, MC simulations are carried out for systems of chain length 4 to 15. Several simulations are performed in the isobaric-isothermal ensemble to obtain the high-density EOS of hard chains in the fluid and solid phases. To determine solid-fluid equilibrium (SFE), Helmholtz free energies of solid crystals at a reference density are determined in a series of canonical MC simulations. As the chain length increases, asymptotic behaviors are observed in the coexistence pressure and densities of fluid and solid phases. It is found that the accuracy of TPT for EOS and SFE in systems of hard chains greatly improves by extending it to second order.

Spontaneous crystallization in athermal polymer packings.

We review recent results from extensive simulations of the crystallization of athermal polymer packings. It is shown that above a certain packing density, and for sufficiently long simulations, all random assemblies of freely-jointed chains of tangent hard spheres of uniform size show a spontaneous transition into a crystalline phase. These polymer crystals adopt predominantly random hexagonal close packed morphologies. An analysis of the local environment around monomers based on the shape and size of the Voronoi polyhedra clearly shows that Voronoi cells become more spherical and more symmetric as the system transits to the ordered state. The change in the local environment leads to an increase in the monomer translational contribution to the entropy of the system, which acts as the driving force for the phase transition. A comparison of the crystallization of hard-sphere polymers and monomers highlights similarities and differences resulting from the constraints imposed by chain connectivity.

An analytical approximation for the orientation-dependent excluded volume of tangent hard sphere chains of arbitrary chain length and flexibility

Onsager-like theories are commonly used to describe the phase behavior of nematic (only orientationally ordered) liquid crystals. A key ingredient in such theories is the orientation-dependent excluded volume of two molecules. Although for hard convex molecular models this is generally known in analytical form, for more realistic molecular models that incorporate intramolecular flexibility, one has to rely on approximations or on computationally expensive Monte Carlo techniques. In this work, we provide a general correlation for the excluded volume of tangent hard-sphere chains of arbitrary chain length and flexibility. The flexibility is introduced by means of the rod-coil model. The resulting correlation is of simple analytical form and accurately covers a wide range of pure component excluded volume data obtained from Monte Carlo simulations of two-chain molecules. The extension to mixtures follows naturally by applying simple combining rules for the parameters involved. The results for mixtures are also in good agreement with data from Monte Carlo simulations. We have expressed the excluded volume as a second order power series in sin (γ), where γ is the angle between the molecular axes. Such a representation is appealing since the solution of the Onsager Helmholtz energy functional usually involves an expansion of the excluded volume in Legendre coefficients. Both for pure components and mixtures, the correlation reduces to an exact expression in the limit of completely linear chains. The expression for mixtures, as derived in this work, is thereby an exact extension of the pure component result of Williamson and Jackson [Mol. Phys. 86, 819-836 (1995)].

Thermodynamic properties of the freely jointed tangent homonuclear Yukawa chain fluid

Recently, a simple analytical expression for the first coordination shell (1 ≤ x ≤ 2) of the average site–site intermolecular radial distribution function (ASSIRDF) of the freely jointed tangent hard sphere chain (FJTHSC) fluid has been proposed. This expression was used in the framework of the second-order Barker–Henderson perturbation theory together with the site–site intermolecular hard-core Yukawa potential in order to obtain an equation of state (EOS) for the freely jointed tangent Yukawa chain (FJTYC) fluid. The derived EOS was used for calculating the thermodynamic properties of the FJTYC fluid, with the attraction of variable range, such as compressibility factor, vapor–liquid equilibria and critical constants. The calculated compressibility factors, above the critical temperature, were compared with the available simulation data and the results of the well-known theories for a wide range of densities, temperatures and different values of attractive ranges. The overall agreement was good and the accuracy of the results could compete with the accuracy of those obtained using very famous theories and equations of state in the literature such as the statistical associate fluid theory with the interaction of variable range (SAFT-VR).

A new simple analytical expression for the average site–site radial distribution function of hard sphere chain fluid

A new simple analytical expression for the first coordination shell (1 ≤ x ≤ 2) of the average site–site intermolecular radial distribution function (ASSIRDF) of freely jointed tangent hard sphere chain fluid (FJTHSC) is proposed. The proposed expression is a polynomial expansion of packing fraction. The dependence of each polynomial coefficient on the radial coordinate and segment number is described by a polynomial expansion of nonlinear base function and a simple relation proposed by Liu and Hu, respectively. The coefficients involved in the expression of ASSIRDF are optimized by fitting the numerical values of Percus–Yevick (PY) integral equation theory of Chiew. The proposed expression of ASSIRDF is applied to an analytical perturbative equation of state (EOS) in order to modify it and improve its accuracy in predicting the compressibility factor and vapor–liquid equilibria of freely jointed tangent square-well chain (FJTSWC) fluids of variable well-widths (1 ≤ λ ≤ 2). The predicted results are in good agreement with the simulation data and also compared with that of other well-known equations of state.

The influence of the polymer chain stiffness on tracer diffusion in polymeric matrices

AbstractThe diffusion of penetrants in polymers is of technological importance in many areas including chromatography and fuel cell membranes. In this work, the effect of chain conformations on tracer diffusion is studied using molecular simulations and a percolation theory. The polymeric matrix is composed of tangent hard sphere chains that are fixed in space; conformations are changed by tuning the stiffness of the chains. The tracer diffusion coefficient is relatively insensitive to the chain stiffness when polymer chains are frozen as in polymer glasses with the local chain dynamics switched off. An analysis of the matrix using percolation theory shows that the polymer volume fraction at the free volume percolation threshold is also relatively insensitive to the chain stiffness, consistent with the diffusion results. This is surprising because the site‐site intermolecular pair correlation functions in the matrix are quite sensitive to the chain stiffness. In contrast, the tracer diffusion coefficient in a melt of mobile chains decreases significantly as the chain stiffness is increased. We conclude that tracer diffusion is only weakly correlated with the chain conformations and local chain dynamics plays an important role. © 2011 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2011

Excluded-volume potential for rigid molecules endowed with C2v symmetry

The excluded volume of a pair of molecules is proportional to the second virial coefficient in hard-core models that represent, for instance, the reference model in perturbation approaches to statistical theories of fluids [see, e.g. Chap. 5 of Kalikmanov, V. (2001), Statistical Physics of Fluids. Texts and Monograph in Physics, Springer, Berlin]. In three space dimensions, there exist exact results for convex molecules and in fact lack of convexity has been a major obstacle in applying the mathematical techniques employed in the convex case. In this paper, we illustrate how a mixed—analytical and numerical—method can be used to obtain exact expressions of the excluded volume for a pair of non-convex molecules conceived as aggregates of hard spheres; these can model van der Waals regions associated to the atoms forming each molecule. To compute the excluded volume for molecules endowed with C2v symmetry, modelled as chains of tangent hard spheres, we adapt a numerical code available to the scientific community. Because the result is a rather cumbersome expression in term of the relative orientation of the interacting molecules, we expand it over a suitable set of symmetry adapted Wigner functions to build up approximate, but faithful expressions, and we also prove analytical results announced elsewhere [Bisi, F., Durand, G., Rosso, R. &amp; Virga, E. (2008), Polar steric interactions for v-shaped molecules. Phys. Rev. E, 78, 011705].

Excluded-Volume Short-Range Repulsive Potential for Tetrahedral Molecules

Tetrahedral symmetry in nematic liquid crystals has been studied for several years, since the seminal paper by Fel [1]. Amongst statistical theories of nematic liquid crystals, a molecular mean-field model has proven to be quite effective in predicting phase sequences for thermotropic biaxial nematic molecules endowed with D2h symmetry, as a function of the parameters entering their interaction potential [2,3]. This model has highlighted the role of a partially repulsive quadrupolar potential of mean torque in promoting condensed phases. It has been shown that the quadrupolar approximation to the excluded-volume interaction between hard spherocuboids can be written precisely as the superposition of two London interactions: one repulsive and one attractive [4]; furthermore, polar steric interactions have been shown to be capable of inducing orientationally ordered states possibly unexpected [5]. By adapting a numerical code available to the scientific community, we evaluate the analytical excluded volume function for non-convex tetrahedral molecules, modelled as chains of tangent hard spheres. Since this function is overly complicated, we expand it over a suitable set of Symmetry Adapted Wigner Functions (SAWFs).