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.

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.

Through extensive Monte Carlo simulations, we systematically study the effect of chain stiffness on the packing ability of linear polymers composed of hard spheres in extremely confined monolayers, corresponding effectively to 2D films. First, we explore the limit of random close packing as a function of the equilibrium bending angle and then quantify the local and global order by the degree of crystallinity and the nematic or tetratic orientational order parameter, respectively. A multi-scale wealth of structural behavior is observed, which is inherently absent in the case of athermal individual monomers and is surprisingly richer than its 3D counterpart under bulk conditions. As a general trend, an isotropic to nematic transition is observed at sufficiently high surface coverages, which is followed by the establishment of the tetratic state, which in turn marks the onset of the random close packing. For chains with right-angle bonds, the incompatibility of the imposed bending angle with the neighbor geometry of the triangular crystal leads to a singular intra- and inter-polymer tiling pattern made of squares and triangles with optimal local filling at high surface concentrations. The present study could serve as a first step toward the design of hard colloidal polymers with a tunable structural behavior for 2D applications.

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

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.