Fine-tuning of colloidal polymer crystals by molecular simulation.

Through extensive molecular simulations we determine a phase diagram of attractive, fully flexible polymer chains in two and three dimensions. A rich collection of distinct crystal morphologies appear, which can be finely tuned through the range of attraction. In three dimensions these include the face-centered cubic, hexagonal close packed, simple hexagonal, and body-centered cubic crystals and the Frank-Kasper phase. In two dimensions the dominant structures are the triangular and square crystals. A simple geometric model is proposed, based on the concept of cumulative neighbors of ideal crystals, which can accurately predict most of the observed structures and the corresponding transitions. The attraction range can thus be considered as an adjustable parameter for the design of colloidal polymer crystals with tailored morphologies.

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).

Polymorph Stability and Free Energy of Crystallization of Freely-Jointed Polymers of Hard Spheres.

The free energy of crystallization of monomeric hard spheres as well as their thermodynamically stable polymorph have been known for several decades. In this work, we present semianalytical calculations of the free energy of crystallization of freely-jointed polymers of hard spheres as well as of the free energy difference between the hexagonal closed packed (HCP) and face-centered cubic (FCC) polymorphs. The phase transition (crystallization) is driven by an increase in translational entropy that is larger than the loss of conformational entropy of chains in the crystal with respect to chains in the initial amorphous phase. The conformational entropic advantage of the HCP polymer crystal over the FCC one is found to be ΔschHCP-FCC≈0.331×10-5k per monomer (expressed in terms of Boltzmann's constant k). This slight conformational entropic advantage of the HCP crystal of chains is by far insufficient to compensate for the larger translational entropic advantage of the FCC crystal, which is predicted to be the stable one. The calculated overall thermodynamic advantage of the FCC over the HCP polymorph is supported by a recent Monte Carlo (MC) simulation on a very large system of 54 chains of 1000 hard sphere monomers. Semianalytical calculations using results from this MC simulation yield in addition a value of the total crystallization entropy for linear, fully flexible, athermal polymers of Δs≈0.93k per monomer.

Local and Global Order in Dense Packings of Semi-Flexible Polymers of Hard Spheres.

The local and global order in dense packings of linear, semi-flexible polymers of tangent hard spheres are studied by employing extensive Monte Carlo simulations at increasing volume fractions. The chain stiffness is controlled by a tunable harmonic potential for the bending angle, whose intensity dictates the rigidity of the polymer backbone as a function of the bending constant and equilibrium angle. The studied angles range between acute and obtuse ones, reaching the limit of rod-like polymers. We analyze how the packing density and chain stiffness affect the chains' ability to self-organize at the local and global levels. The former corresponds to crystallinity, as quantified by the Characteristic Crystallographic Element (CCE) norm descriptor, while the latter is computed through the scalar orientational order parameter. In all cases, we identify the critical volume fraction for the phase transition and gauge the established crystal morphologies, developing a complete phase diagram as a function of packing density and equilibrium bending angle. A plethora of structures are obtained, ranging between random hexagonal closed packed morphologies of mixed character and almost perfect face centered cubic (FCC) and hexagonal close-packed (HCP) crystals at the level of monomers, and nematic mesophases, with prolate and oblate mesogens at the level of chains. For rod-like chains, a delay is observed between the establishment of the long-range nematic order and crystallization as a function of the packing density, while for right-angle chains, both transitions are synchronized. A comparison is also provided against the analogous packings of monomeric and fully flexible chains of hard spheres.

Polymorphism and Perfection in Crystallization of Hard Sphere Polymers.

We present results on polymorphism and perfection, as observed in the spontaneous crystallization of freely jointed polymers of hard spheres, obtained in an unprecedentedly long Monte Carlo (MC) simulation on a system of 54 chains of 1000 monomers. Starting from a purely amorphous configuration, after an initial dominance of the hexagonal closed packed (HCP) polymorph and a transitory random hexagonal close packed (rHCP) morphology, the system crystallizes in a final, stable, face centered cubic (FCC) crystal of very high perfection. An analysis of chain conformational characteristics, of the spatial distribution of monomers and of the volume accessible to them shows that the phase transition is caused by an increase in translational entropy that is larger than the loss of conformational entropy of the chains in the crystal, compared to the amorphous state. In spite of the significant local re-arrangements, as reflected in the bending and torsion angle distributions, the average chain size remains unaltered during crystallization. Polymers in the crystal adopt ideal random walk statistics as their great length renders local conformational details, imposed by the geometry of the FCC crystal, irrelevant.

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.

Simu-D: A Simulator-Descriptor Suite for Polymer-Based Systems under Extreme Conditions.

We present Simu-D, a software suite for the simulation and successive identification of local structures of atomistic systems, based on polymers, under extreme conditions, in the bulk, on surfaces, and at interfaces. The protocol is built around various types of Monte Carlo algorithms, which include localized, chain-connectivity-altering, identity-exchange, and cluster-based moves. The approach focuses on alleviating one of the main disadvantages of Monte Carlo algorithms, which is the general applicability under a wide range of conditions. Present applications include polymer-based nanocomposites with nanofillers in the form of cylinders and spheres of varied concentration and size, extremely confined and maximally packed assemblies in two and three dimensions, and terminally grafted macromolecules. The main simulator is accompanied by a descriptor that identifies the similarity of computer-generated configurations with respect to reference crystals in two or three dimensions. The Simu-D simulator-descriptor can be an especially useful tool in the modeling studies of the entropy- and energy-driven phase transition, adsorption, and self-organization of polymer-based systems under a variety of conditions.

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.

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.

The role of bond tangency and bond gap in hard sphere crystallization of chains.

We report results from Monte Carlo simulations on dense packings of linear, freely-jointed chains of hard spheres of uniform size. In contrast to our past studies where bonded spheres along the chain backbone were tangent, in the present work a finite tolerance in the bond is allowed. Bond lengths are allowed to fluctuate in the interval [σ, σ + dl], where σ is the sphere diameter. We find that bond tolerance affects the phase behaviour of hard-sphere chains, especially in the close vicinity of the melting transition. First, a critical dl(crit) exists marking the threshold for crystallization, whose value decreases with increasing volume fraction. Second, bond gaps enhance the onset of phase transition by accelerating crystal nucleation and growth. Finally, bond tolerance has an effect on crystal morphologies: in the tangent limit the majority of structures correspond to stack-faulted random hexagonal close packing (rhcp). However, as bond tolerance increases a wealth of diverse structures can be observed: from single fcc (or hcp) crystallites to random hcp/fcc stackings with multiple directions. By extending the simulations over trillions of MC steps (10(12)) we are able to observe crystal-crystal transitions and perfection even for entangled polymer chains in accordance to the Ostwald's rule of stages in crystal polymorphism. Through simple geometric arguments we explain how the presence of rigid or flexible constraints affects crystallization in general atomic and particulate systems. Based on the present results, it can be concluded that proper tuning of bond gaps and of the connectivity network can be a controlling factor for the phase behaviour of model, polymer-based colloidal and granular systems.

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.

Numerical simulation of bubble dynamics in a Phan-Thien-Tanner liquid: non-linear shape and size oscillatory response under periodic pressure.

We study non-linear bubble oscillations driven by an acoustic pressure with the bubble being immersed in a viscoelastic, Phan-Thien-Tanner liquid. Solution is provided numerically through a method which is based on a finite element discretization of the Navier-Stokes flow equations. The proposed computational approach does not rely on the solution of the simplified Rayleigh-Plesset equation, is not limited in studying only spherically symmetric bubbles and provides coupled solutions for the velocity, stress fields and bubble interface. We present solutions for non-spherical bubbles, with asphericity being addressed by means of Legendre polynomials or associated Legendre functions. A parametric investigation of the bubble dynamical oscillatory response as a function of the fluid rheological properties shows that the amplitude of bubble oscillations drastically increases as liquid elasticity (quantified by the Deborah number) increases or as liquid viscosity decreases (quantified by the Reynolds number). Extensive numerical calculations demonstrate that increasing elasticity and/or viscosity of the surrounding liquid tend to stabilize the shape anisotropy of an initially non-spherical bubble. Results are shown for pressure amplitudes 0.2-2MPa and Deborah, Reynolds numbers in the intervals of 1-8 and 0.094-1.256, respectively.

Topological analysis of polymeric melts: chain-length effects and fast-converging estimators for entanglement length.

Primitive path analyses of entanglements are performed over a wide range of chain lengths for both bead spring and atomistic polyethylene polymer melts. Estimators for the entanglement length N_{e} which operate on results for a single chain length N are shown to produce systematic O(1/N) errors. The mathematical roots of these errors are identified as (a) treating chain ends as entanglements and (b) neglecting non-Gaussian corrections to chain and primitive path dimensions. The prefactors for the O(1/N) errors may be large; in general their magnitude depends both on the polymer model and the method used to obtain primitive paths. We propose, derive, and test new estimators which eliminate these systematic errors using information obtainable from the variation in entanglement characteristics with chain length. The new estimators produce accurate results for N_{e} from marginally entangled systems. Formulas based on direct enumeration of entanglements appear to converge faster and are simpler to apply.

Contact network in nearly jammed disordered packings of hard-sphere chains.

We present salient results of the analysis of the geometrical structure of a large fully equilibrated ensemble of nearly jammed packings of linear freely jointed chains of tangent hard spheres generated via extensive Monte Carlo simulations. In spite the expected differences due to chain connectivity, both the pair-correlation function and the contact network for chain packings are found to strongly resemble those in packings of monomeric hard spheres at the maximally random jammed (MRJ) state. A remarkable finding of the present work is the tendency of chains to form closed loops at the MRJ state as a consequence of chain collapse. Our simulations on disordered nearly jammed chain packings yield an average coordination number of 6, which fulfills the isostaticity condition and is in excellent agreement with the corresponding simulation [A. Donev, S. Torquato, and F. H. Stillinger, Phys. Rev. E 71, 011105 (2005)] and experimental [T. Aste, M. Saadatfar, and T. J. Senden, Phys. Rev. E 71, 061302 (2005)] findings for jammed packings of monatomic hard spheres. An exact correspondence between the statistical-mechanical ensembles of monomeric spheres and of hard-sphere chains offers insights regarding the structure and topology of the contact network of hard-sphere systems at the MRJ state.

Entropy-driven crystallization in dense systems of athermal chain molecules.

We describe the direct observation of entropy-driven crystallization in simulations of dense packings of linear hard-sphere chains. Crystal nuclei form spontaneously in the phase coexistence region independently of chain length. Incipient nuclei consistently develop well defined, stack-faulted layered morphologies with a single stacking direction. These morphologies deviate markedly from those of monomeric analogs. The ordering transition is driven by the increase in translational entropy: ordered sites exhibit enhanced mobility as their local free volume becomes more spherical and symmetric.

The structure of random packings of freely jointed chains of tangent hard spheres.

We analyze the structure of dense random packings of freely jointed chains of tangent hard spheres as a function of concentration (packing density) with particular emphasis placed on the behavior in the vicinity of their maximally random jammed (MRJ) state. Representative configurations over the whole density range are generated through extensive off-lattice Monte Carlo simulations on systems of average chain lengths ranging from N=12 to 1000 hard spheres. Several measures of order are used to quantitatively describe either local structure (sphere arrangements and bonded geometry) or global behavior (chain conformations and statistics). In addition, the employed measures are used to elucidate the effect of connectivity on structure, by comparing monatomic and chain assemblies of hard spheres at the MRJ state.

The characteristic crystallographic element norm: A descriptor of local structure in atomistic and particulate systems.

We introduce the characteristic crystallographic element (CCE) norm as a powerful descriptor of local structure in atomistic and particulate systems. The CCE-norm is sensitive both to radial and orientational deviations from perfect local order. Unlike other measures of local order, the CCE-norm decreases monotonically with increasing order, is zero for a perfectly ordered environment, and is strictly discriminating among different, competing crystal structures in imperfectly ordered systems. The CCE-norm descriptor can be used as a sensitive, quantitative measure to detect and track changes in local order in atomistic and general particulate systems. In a specific example we show the ability of the CCE-norm to monitor the onset and evolution of order in an initially amorphous, densely packed assembly of hard-sphere chains generated through extensive Monte Carlo simulations [Phys. Rev. Lett. 100, 050602 (2008)].

Structure, dimensions, and entanglement statistics of long linear polyethylene chains.

This work elucidates the effect of both temperature and molecular length on the conformational and structural properties as well as on the entanglement statistics of long amorphous, polydisperse, and molten linear polyethylene (PE). A large number of PE samples are modeled in atomistic detail, with average molecular lengths ranging from C24 up to C1,000 over a wide range of temperatures in the interval of 300

Modeling the effect of cell-associated polymeric fluid layers on force spectroscopy measurements. Part II: experimental results and comparison with model predictions.

In this paper, experimentally obtained force curves on Staphylococcus aureus are compared with a previously developed model that incorporates hydrodynamic effects of extracellular polysaccharides together with the elastic response of the bacterium and cantilever. Force-displacement curves were predicted without any adjustable parameters. It is demonstrated that experimental results can be accurately described by our model, especially if viscoelastic effects of the extracellular polysaccharide layer are taken into account. Polysaccharide layer viscoelasticity was treated by means of a multimode Phan-Thien/Tanner (PTT) constitutive equation. Typical maximum relaxation times range from 0.2 to 2 s, whereas the corresponding zero-shear-rate viscosities are 6-9 Pa.s, based on published, steady-state rheological measurements on Staphylococcus aureus polysaccharide extracted from its native environment. The bacterial elastic constant is found to be in the range 0.02-0.4 N/m, corresponding to bacterial wall Young's moduli in the range of a few hundred MPa. Repeatability of measurements performed on different bacteria is found to be only fair, due to large individuum variability, whereas repetitions of measurements on the same bacterium showed high reproducibility. Improved force-indentation curve predictions are expected if transient rheological characterization of extracellular polysaccharides is available. More desirable however is the direct, in vivo rheological characterization of the extracellular polysaccharide. A model-based analysis of experimental force-indentation curves shows that appreciable further experimental improvements are still necessary to achieve this goal.

Modeling the effect of cell-associated polymeric fluid layers on force spectroscopy measurements. Part I: model development.

The mechanical response, the force-indentation relationship, in normal force spectroscopy measurements carried out on individual polysaccharide encapsulated bacteria is modeled using three increasingly refined approaches that consider the elastic response of the bacterium and cantilever in combination with a fluid (hydrodynamic) model for the polysaccharide layer. For the hydrodynamic description of the polysaccharide layer, several increasingly realistic models are described in detail, together with numerical solution techniques. These models range from one-dimensional, Newtonian, to two-dimensional, axisymmetric, fully viscoelastic (Phan-Thien/Tanner). In all cases, the models rigorously consider the time-dependent rheological-mechanical coupling between the elastic and fluid viscoelastic physical components of the experimental setup. Effects of inherent variability in geometrical and material properties of the bacterium and polysaccharide layer on the measurable response are quantified. A parametric investigation of the force-indentation relationship highlights the importance of accurate knowledge of the rheology of the extracellular polysaccharides. We also draw conclusions about the design and evaluation of force spectroscopy experiments on single encapsulated bacteria. Supported by model calculations, we also point the way to methods of in vivo rheological characterization of the extracellular polysaccharide as a preferable alternative to characterization after its removal from the native environment.