Fivefold symmetry as an inhibitor to hard-sphere crystallization

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.

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.

New phase transition pattern of fivefold twins transformed into lamellar structure in Ti3Al alloy

A new phase transition pattern of fivefold twins into a lamellar structure leading to a second phase transition was found in Ti3Al alloy.

Identification of Local Structure in 2-D and 3-D Atomic Systems through Crystallographic Analysis

In the present work, we revise and extend the Characteristic Crystallographic Element (CCE) norm, an algorithm used to simultaneously detect radial and orientational similarity of computer-generated structures with respect to specific reference crystals and local symmetries. Based on the identification of point group symmetry elements, the CCE descriptor is able to gauge local structure with high precision and finely distinguish between competing morphologies. As test cases we use computer-generated monomeric and polymer systems of spherical particles interacting with the hard-sphere and square-well attractive potentials. We demonstrate that the CCE norm is able to detect and differentiate, between others, among: hexagonal close packed (HCP), face centered cubic (FCC), hexagonal (HEX) and body centered cubic (BCC) crystals as well as non-crystallographic fivefold (FIV) local symmetry in bulk 3-D systems; triangular (TRI), square (SQU) and honeycomb (HON) crystals, as well as pentagonal (PEN) local symmetry in thin films of one-layer thickness (2-D systems). The descriptor is general and can be applied to identify the symmetry elements of any point group for arbitrary atomic or particulate system in two or three dimensions, in the bulk or under confinement.

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.

Slowing down supercooled liquids by manipulating their local structure

Glasses remain an elusive and poorly understood state of matter. It is not clear how we can control the macroscopic dynamics of glassy systems by tuning the properties of their microscopic building blocks. In this paper, we propose a simple directional colloidal model that reinforces the optimal icosahedral local structure of binary hard-sphere glasses. We show that this specific symmetry results in a dramatic slowing down of the dynamics. Our results open the door to controlling the dynamics of dense glassy systems by selectively promoting specific local structural environments.

Complex Crystals from Size-Disperse Spheres

Colloids are rarely perfectly uniform but follow a distribution of sizes, shapes, and charges. This dispersity can be inherent (static) or develop and change over time (dynamic). Despite a long history of research, the conditions under which nonuniform particles crystallize and which crystal forms is still not well understood. Here, we demonstrate that hard spheres with Gaussian radius distribution and dispersity up to 19% always crystallize if compressed slowly enough, and they do so in surprisingly complex ways. This result is obtained by accelerating event-driven simulations with particle swap moves for static dispersity and particle resize moves for dynamic dispersity. Above 6% dispersity, AB_{2} Laves, AB_{13}, and a region of Frank-Kasper phases are found. The Frank-Kasper region includes a quasicrystal approximant with Pearson symbol oS276. Our findings are relevant for ordering phenomena in soft matter and alloys.

Confined Polymers as Self-Avoiding Random Walks on Restricted Lattices

Polymers in highly confined geometries can display complex morphologies including ordered phases. A basic component of a theoretical analysis of their phase behavior in confined geometries is the knowledge of the number of possible single-chain conformations compatible with the geometrical restrictions and the established crystalline morphology. While the statistical properties of unrestricted self-avoiding random walks (SAWs) both on and off-lattice are very well known, the same is not true for SAWs in confined geometries. The purpose of this contribution is (a) to enumerate the number of SAWs on the simple cubic (SC) and face-centered cubic (FCC) lattices under confinement for moderate SAW lengths, and (b) to obtain an approximate expression for their behavior as a function of chain length, type of lattice, and degree of confinement. This information is an essential requirement for the understanding and prediction of entropy-driven phase transitions of model polymer chains under confinement. In addition, a simple geometric argument is presented that explains, to first order, the dependence of the number of restricted SAWs on the type of SAW origin.

Coupling between criticality and gelation in "sticky" spheres: a structural analysis

We combine experiments and simulations to study the link between criticality and gelation in sticky spheres. We employ confocal microscopy to image colloid-polymer mixtures and Monte Carlo simulations of the square-well (SW) potential as a reference model. To this end, we map our experimental samples onto the SW model. We find an excellent structural agreement between experiments and simulations, both for locally favored structures at the single particle level and large-scale fluctuations at criticality. We follow in detail the rapid structural change in the critical fluid when approaching the gas-liquid binodal and highlight the role of critical density fluctuations for this structural crossover. Our results link the arrested spinodal decomposition to long-lived energetically favored structures, which grow even away from the binodal due to the critical scaling of the bulk correlation length and static susceptibility.

disLocate: tools to rapidly quantify local intermolecular structure to assess two-dimensional order in self-assembled systems

Order classification is particularly important in photonics, optoelectronics, nanotechnology, biology, and biomedicine, as self-assembled and living systems tend to be ordered well but not perfectly. Engineering sets of experimental protocols that can accurately reproduce specific desired patterns can be a challenge when (dis)ordered outcomes look visually similar. Robust comparisons between similar samples, especially with limited data sets, need a finely tuned ensemble of accurate analysis tools. Here we introduce our numerical Mathematica package disLocate, a suite of tools to rapidly quantify the spatial structure of a two-dimensional dispersion of objects. The full range of tools available in disLocate give different insights into the quality and type of order present in a given dispersion, accessing the translational, orientational and entropic order. The utility of this package allows for researchers to extract the variation and confidence range within finite sets of data (single images) using different structure metrics to quantify local variation in disorder. Containing all metrics within one package allows for researchers to easily and rapidly extract many different parameters simultaneously, allowing robust conclusions to be drawn on the order of a given system. Quantifying the experimental trends which produce desired morphologies enables engineering of novel methods to direct self-assembly.

Anisotropic stress inhibits crystallization in Cu-Zr glass-forming liquids

Liquids attain a metastable state without crystallizing by cooling rapidly to a given temperature below the melting point. With increasing supercooling, the nucleation rate would show an increase based on the prediction of the classical nucleation theory. It is generally thought that the nucleation rate will reach the maximum upon approaching the glass transition temperature, T_g, for glass-forming liquids. We report that there exists a supercooled region above T_g in which the crystallization has actually been severely suppressed. Our molecular dynamics simulations show that the growth of embryos in the supercooled Cu₆₀Zr₄₀ melt is subjected to a strong anisotropic stress associated with the dynamic heterogeneity. Its long-range effect drives the embryo to grow into a ramified morphology so that the interface energy dominates over the embryo growth, leading to the suppression of nucleation.

Local rotational symmetry in the packing of uniform spheres

Local rotational symmetry (LRS) of a particulate system is important for understanding its structure and phase transition. However, how to properly characterize LRS for this system is still a challenge as the system normally includes both ordered and disordered local structures. Herein, based on the so-called common neighbour subcluster (CNS), we proposed a method to characterize the LRS of uniform spheres packings with the packing fraction ρ ranging within 0.20 and 0.74. It was found that different fold LRSs coexist in most packings, and their maximum degree increases at ρ < 0.64, except for the 2-fold LRS held by 6-sphere CNS that continuously increases to form the fcc crystal at ρ = 0.74. The overall LRS involving all the CNSs monotonically increases with two critical changes at ρ = (0.35-0.40) and 0.64; the evolution of individual LRSs held by specific CNS groups critically changes at ρ ≈ (0.35-0.40), 0.50, 0.55-0.60, and 0.64. The physics corresponding to these critical changes has also been discussed. The findings will significantly enrich the understanding of the structural symmetry of materials including atoms and particles.

The role of fivefold symmetry in suppressing crystallization

Although long assumed to have an important role in the suppression of crystallization and the development of glassformers, the effect of local fivefold symmetry has never been directly tested. Here we consider whether such suppression of crystallization has a kinetic or thermodynamic nature and investigate its mechanism. We introduce a model in which the degree of fivefold symmetry can be tuned by favouring arrangements of particles in pentagonal bipyramids. We thus show that fivefold symmetry has both kinetic and thermodynamic effects on the mechanism of crystallization to a face-centred cubic crystal. Our results suggest that the mechanism of crystallization suppression is related to the surface tension between fluid and crystal. Interestingly, the degree of fivefold symmetry has little effect on crystal growth rate, suggesting that growth may be only weakly coupled to fluid structure in hard sphere like systems. Upon increasing the fivefold symmetry, we find a first-order transition to an alternative icosahedra-rich phase. At intermediate bias strengths we find a one-component glassformer.

The suppression of crystallization due to the appearance of structures with fivefold symmetry is widely adopted, but its kinetic and thermodynamic origin remains elusive. Taffs et al. show that fivefold symmetry substantially slows down the nucleation rate but not the crystal growth rate as expected.

Effect of chain stiffness on the competition between crystallization and glass-formation in model unentangled polymers

We map out the solid-state morphologies formed by model soft-pearl-necklace polymers as a function of chain stiffness, spanning the range from fully flexible to rodlike chains. The ratio of Kuhn length to bead diameter (lK/r0) increases monotonically with increasing bending stiffness kb and yields a one-parameter model that relates chain shape to bulk morphology. In the flexible limit, monomers occupy the sites of close-packed crystallites while chains retain random-walk-like order. In the rodlike limit, nematic chain ordering typical of lamellar precursors coexists with close-packing. At intermediate values of bending stiffness, the competition between random-walk-like and nematic chain ordering produces glass-formation; the range of kb over which this occurs increases with the thermal cooling rate |Ṫ| implemented in our molecular dynamics simulations. Finally, values of kb between the glass-forming and rodlike ranges produce complex ordered phases such as close-packed spirals. Our results should provide a useful initial step in a coarse-grained modeling approach to systematically determining the effect of chain stiffness on the crystallization-vs-glass-formation competition in both synthetic and colloidal polymers.

Crystallization and arrest mechanisms of model colloids

We performed dynamic simulations of spheres with short-range attractive interactions for many values of interaction strength and range. Fast crystallization occurs in a localized region of this parameter space, but the character of crystallization pathways is not uniform within this region. Pathways range from one-step, in which a crystal nucleates directly from a gas, to two-step, in which substantial liquid-like clusters form and only subsequently become crystalline. Crystallization can fail because of slow nucleation from either gas or liquid, or because of dynamic arrest caused by strong interactions. Arrested states are characterized by the formation of networks of face-sharing tetrahedra that can be detected by a local common neighbor analysis.

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 structure transition in nanoparticle aggregates: from amorphous clusters to super-crystals

Spontaneous structure transition is studied in a real NP system, which reveals some important details of this transition.

Geometrical frustration in amorphous and partially crystallized packings of spheres

We study the persistence of a geometrically frustrated local order inside partially crystallized packings of equal-sized spheres. Measurements by x-ray tomography reveal previously unseen grain scale rearrangements occurring inside large three-dimensional packings as they crystallize. Three successive structural transitions are detected by a statistical description of the local volume fluctuations. These compaction regimes are related to the disappearance of densely packed tetrahedral patterns of beads. Amorphous packings of monodisperse spheres are saturated with these tetrahedral clusters at Bernal's limiting density (φ≈64%). But, no periodic lattice can be built upon these patterns; they are geometrically frustrated and are thus condemned to vanish while the crystallization occurs. Remarkably, crystallization-induced grain rearrangements can be interpreted in terms of the evolution of key topological features of these aggregates.

Simple model for chain packing and crystallization of soft colloidal polymers

We study a simple bead-spring polymer model exhibiting competing crystallization and glass transitions. Constant-pressure molecular dynamics simulations are employed to study phase behavior and morphological order. For adequately slow quench rates, chain systems exhibit a first-order phase transition (crystallization) below a critical temperature T=T(cryst). We observe the formation of close-packed crystallites of FCC and/or HCP order, separated by domain walls, twin defects, and amorphous regions. Such crystal structures closely resemble the corresponding ordered morphologies of athermal polymer packings: fully flexible chains retain random-walk-like configurations in the crystalline state and do not form lamellae, while semiflexible chains do form lamellae. The model presented here is well suited to the modeling of granular and colloidal polymers, in particular for elucidating the factors that dictate the formation of specific ordered morphologies.

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.

The microscopic pathway to crystallization in supercooled liquids

Despite its fundamental and technological importance, a microscopic understanding of the crystallization process is still elusive. By computer simulations of the hard-sphere model we reveal the mechanism by which thermal fluctuations drive the transition from the supercooled liquid state to the crystal state. In particular we show that fluctuations in bond orientational order trigger the nucleation process, contrary to the common belief that the transition is initiated by density fluctuations. Moreover, the analysis of bond orientational fluctuations shows that these not only act as seeds of the nucleation process, but also i) determine the particular polymorph which is to be nucleated from them and ii) at high density favour the formation of fivefold structures which can frustrate the formation of crystals. These results can shed new light on our understanding of the relationship between crystallization and vitrification.

Structure of finite sphere packings via exact enumeration: implications for colloidal crystal nucleation

We analyze the geometric structure and mechanical stability of a complete set of isostatic and hyperstatic sphere packings obtained via exact enumeration. The number of nonisomorphic isostatic packings grows exponentially with the number of spheres N, and their diversity of structure and symmetry increases with increasing N and decreases with increasing hyperstaticity H≡N_{c}-N_{ISO}, where N_{c} is the number of pair contacts and N_{ISO}=3N-6. Maximally contacting packings are in general neither the densest nor the most symmetric. Analyses of local structure show that the fraction f of nuclei with order compatible with the bulk (rhcp) crystal decreases sharply with increasing N due to a high propensity for stacking faults, five- and near-fivefold symmetric structures, and other motifs that preclude rhcp order. While f increases with increasing H, a significant fraction of hyperstatic nuclei for N as small as 11 retain non-rhcp structure. Classical theories of nucleation that consider only spherical nuclei, or only nuclei with the same ordering as the bulk crystal, cannot capture such effects. Our results provide an explanation for the failure of classical nucleation theory for hard-sphere systems of N≲10 particles; we argue that in this size regime, it is essential to consider nuclei of unconstrained geometry. Our results are also applicable to understanding kinetic arrest and jamming in systems that interact via hard-core-like repulsive and short-ranged attractive interactions.