Katerina Foteinopoulou

We conduct extensive Monte Carlo simulations to investigate the factors that control the isotropic-to-nematic transition of hard colloidal polymers in bulk and under various conditions of confinement. Utilizing a highly idealized model, polymers are represented as linear chains of tangent hard spheres of uniform length, whose stiffness is controlled by a bending potential leading to rod-like configurations. Confinement is realized through the presence of flat, parallel, and impenetrable walls in one, two, or three dimensions while periodic boundary conditions are applied on the unconstrained dimensions. All simulations are performed through the Simu-D software, composed of conventional and advanced, chain-connectivity-altering Monte Carlo algorithms. We explore in detail how distinct factors, including chain length, stiffness, confinement, and packing density affect the isotropic-to-nematic transition exhibited by the polymer chains and identify with high precision the concentration range where this phase change takes place as a function of the applied 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.

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.

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

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.

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. 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 from acute to obtuse ones, reaching the limit of rod-like polymers. We analyze how packing density and chain stiffness affect the ability of chains 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 is obtained, ranging from 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, to nematic mesophases, with prolate and oblate mesogens at the level of chains. For rod-like chains, hysteresis is observed between the establishment of long-range nematic order and crystallization, 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 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.

We present results on 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 and after a transitory dominance of the hexagonal closed packed (HCP) polymorph, the system crystallizes in a final, stable, face centered cubic (FCC) crystal of very high perfection. Through refined metrics we gauge the degree of ordering and identify the regions of the phase transition and the corresponding morphologies. 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. Polymer chains 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.

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.