Preprint
Article

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

Altmetrics

Downloads

155

Views

72

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

29 December 2022

Posted:

03 January 2023

You are already at the latest version

Alerts
Abstract
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.
Keywords: 
Subject: Chemistry and Materials Science  -   Polymers and Plastics

1. Introduction

Over the last decades, developments in the synthesis of novel polymers and the fabrication of polymer-based materials have turned them into key components of our daily lives. The research is ever-growing in the pursuit of polymer-based materials with enhanced properties. Chain connectivity endows macromolecules with unique properties compared to monoatomic systems [1,2]. Describing polymer conformations and configurations statistically and understanding how these are connected to macroscopic properties are of paramount importance in technology and industry [3,4]. In parallel, many aspects of the phase behavior and self-organization of general atomic and particulate systems remain unknown or poorly understood. The analysis becomes a great deal harder when macromolecular systems are tackled. Advances in experimental, theoretical, and simulation methods continually enrich our fundamental knowledge of the phenomenon in a wide range of physical systems [5,6,7,8,9,10,11]. For example, theoretical models have been developed to predict the effect of the processing conditions on the phase transition in order to analyze the final morphologies [12,13]. Advances in the synthesis and characterization of colloidal and granular systems have provided significant insights on the phase behavior of monomeric and polymeric systems, especially given their simplicity and large size compared to traditional polymers [14,15,16,17,18,19,20,21,22]. Simulations can further aid in the research studies of phase transition given the constant advances at the level of hardware and software, starting from the traditional approaches at the molecular level based on Molecular Dynamics (MD) or Monte Carlo (MC) algorithms [23,24,25,26,27].
The hard-sphere model is frequently used in the study of crystallization due to its simplicity and athermal nature, despite the obvious disadvantage of lacking chemical detail. Through molecular dynamics (MD) simulations, early and pioneering works have reported the entropy-driven crystallization of systems of hard spheres [12,28] once the melting point φ m o n o m e r s M = 0.545 is reached [29] and given sufficient time for the observation of the phenomenon [30,31]. For individual hard spheres, the crystal expected to be obtained in experiments or simulations should be the face centered cubic (FCC), as it is slightly more stable thermodynamically than the hexagonal close-packed (HCP) one [32,33,34]. However, random hexagonal close packing (RHCP) is almost always observed due to the small free energy of stacking faults [35]. Thus, crystal perfection in the form of pure FCC crystals is rarely encountered in computer simulations [36,37].
The mechanism of crystallization of hard colloidal polymers [37,38,39,40,41,42,43,44,45] differs substantially from that of traditional polymers, the latter as revealed in x-ray scattering studies of short alkane chains [46,47,48], and in MD [49,50,51,52] or MC [39,53] simulations. Just as for monomeric athermal hard spheres, off-lattice MC simulations have shown that random packings of fully flexible linear chains of tangent hard spheres are able to crystallize at high volume fractions through an entropy-driven mechanism, very similar to the one of monomeric analogs [37,40,41,42,54]. Extremely long simulations [37] have shown the FCC polymorph to be the stable one for fully flexible chains of tangent hard spheres. In shorter simulations, an RHCP phase with a unique stacking direction [40,41,42] appears most often. In spite of the similarities to monomeric counterparts, athermal polymer crystallization shows unique characteristics: The critical volume fraction for the phase transition (melting point) for the fully flexible (freely-jointed, FJ) chains of tangent hard spheres, φ c h a i n s M ( F J ) (> 0.56), is higher than the melting point of monomers. Monomeric systems present diverse crystal morphologies, ranging from single FCC or HCP structures to close-stacked packings with random directions where twin defects, associated with the formation of fivefold (FIV) symmetry structures, kinetically frustrate crystallization [55,56]. However, both the FCC or RHCP crystals of tangent hard-sphere chains are usually free of twin defects due to a conformational entropic barrier [57]. It has been established that the phase transition and the ordered morphologies of hard-sphere chains are affected insignificantly by chain length, but are sensitive to factors such as gaps in bond lengths [43,58], the presence of surfaces/confinement [44,59,60] and chain stiffness [61,62,63,64].
For semi-flexible polymers, the rigidity of the chains is generally implemented in molecular simulations through bending [65,66,67,68,69,70,71,72] and/or torsional [71] potentials. A common topic of study about the phase behavior of semi-flexible polymers is the orientational (nematic) ordering of the chains. According to Onsager, thin and infinitely long rigid rods undergo a transition from the isotropic to the nematic phase which is driven by entropy, resulting from the competition between translational and orientational ordering [73]. Building on this, Bolhuis and Frenkel mapped out the complete phase diagram of hard spherocylinders as a function of shape anisotropy [74]. The long-range, nematic order in solutions of semi-flexible polymers depends on the size ratio of the chain [75,76,77]. For tangent hard-sphere chains in the rod limit, the transition occurs at high concentrations, still lower than the solidification of the system [78,79]. Theories have been developed [80,81,82] and studies have been conducted [83,84,85,86,87] under various conditions on the isotropic-to-nematic transition of semi-flexible chains. Nguyen et al. studied the effect of chain stiffness and temperature on the competition between crystallization and glass formation for unentangled semi-flexible polymer melts [61,62]. Shakirov and Paul studied the crystallization in melts of short, semi-flexible hard-sphere chains employing Monte Carlo simulations [63,64]. In their works, rod-like polymer melts undergo a first-order transition at a density lower than a melt of hard-sphere monomers or flexible hard-sphere chains [64]. During crystallization, the orientational ordering (nematic) transition is accompanied by a 2-D translational ordering transition synchronously in the plane perpendicular to the chains. These orderings induce the formation of hexagonal crystal planes [63].
Although rod-like chains have been extensively researched in the literature, the corresponding body of work on semi-flexible chains with different bending angles is very limited. Recent works studying the jamming and solidification of semi-flexible polymers employing the freely-rotating model varying the fixed bending angle through MD simulations have revealed the different mechanisms in solidification depending on the equilibrium bending angle, θ 0 [88,89].
The present manuscript analyzes the phase behavior of semi-flexible linear chains of tangent hard spheres, as a function of equilibrium bending angle and packing density. This is achieved through extensive Monte Carlo (MC) simulations using the Simu-D simulator-descriptor [90], built around polymer-oriented algorithms and especially chain-connectivity-altering moves (CCAMs) [91,92]. Motivated by previous studies on the crystallization of linear, fully flexible chains of hard spheres [40,41,42,43,44,93,94], we extend the research by the inclusion of chain stiffness via a bending potential that depends on a bending constant, k θ , and an equilibrium bending angle, θ 0 . Global order, at the level of chains, is quantified through the orientational (nematic) order parameter, q. Local order, at the level of monomers, is gauged through the Characteristic Crystallographic Element (CCE) norm descriptor [95,96]. Compiling the results of all the simulated systems, we present a phase diagram that reflects the combined effect of packing density and chain stiffness on the local and global order of the semi-flexible systems.
This article is organized as follows: in Section 2 we describe the molecular model, and the simulated systems and we briefly explain the descriptors to gauge long-range order and local structure. Section 3 presents the results on the phase behavior of the simulated semi-flexible systems. Finally, we summarize the main conclusions and discuss current efforts in Section 4.

2. Model and Methods

2.1. Molecular Model and Systems studied

Polymers are modeled as linear chains of identical hard spherical monomers with a collision diameter σ, which is considered the characteristic length. The systems are composed of Nat monomers distributed in Nch chains of average chain length of Nav (in monomers). All systems consist of Nch = 100 chains of average chain length Nav = 12, resulting in Nat = 1200. Chains present dispersity in their lengths following a uniform distribution in the range N 6 ,   18 due to the presence of chain-connectivity-altering MC moves (see below).
The non-overlapping condition of the hard monomers is adopted by employing the Hard Sphere (HS) potential to describe all non-bonded interactions between monomers. According to the HS potential, the pair-wise energy, U H S r i j , is determined by
U H S r i j =   0 ,   r i j σ   ,   r i j < σ
where rij is the distance between the centers of the monomers i and j. Periodic boundary conditions are applied in all dimensions, corresponding to a bulk, unconstraint system. For non-overlapping objects packing density (volume fraction), φ, is defined as the total volume occupied by monomers over the volume of the simulation cell. This is given by:
φ = V m o n V c e l l = π 6 N a t V c e l l σ 3
where Vmon is the total volume occupied by monomers and Vcell is the volume of the cubic simulation cell.
Figure 1 shows a sketch of the bond geometry encountered in linear chains, consisting of bond lengths, l, bending angles, θ, and torsion angles, ϕ. Bond lengths can fluctuate uniformly in the interval l σ ,   σ + d l , where dl is the maximum gap between successive monomers of the chain. For all systems simulated here, d l = 6.5 × 10 4 (in units of σ), practically enforcing tangency between bonded monomers [43]. Chain stiffness is introduced through a potential controlling the bending angle, θ, which is the angle formed by two successive bond vectors (Figure 1). The bending potential is given by
U b e n d θ = k θ θ θ 0 2
where k θ is the bending constant and θ 0 is the equilibrium bending angle. Setting k θ = 0 with fixed bond lengths allows the simulation of freely-jointed chains while k θ corresponds to freely-rotating chains. The equilibrium bending angle can vary from fully extended ( θ 0 = 0 ° ) to fully compact ( θ 0 = 120 ° ) configurations of triplets. In this work, we study five specific equilibrium bending angles: 0°, 60°, 90°, 108°, and 120°, as they are represented schematically in Figure 1. Torsion angles, φ , are allowed to fluctuate freely.
The bending constant is set equal to k θ / k B T = 9   rad-2 (for simplicity in the continuation we will report it as k θ = 9 ) in all simulated systems, where k B and T are the Boltzmann constant and temperature, respectively. The specific choice of the harmonic constant was made considering computational expediency and the fact that this range of constraints leads to orientational ordering for rod-like chains as demonstrated in [63]. This degree of stiffness limits the deviations from the equilibrium bending angle to, approximately, ±20°, as shown in Figure 2, where for comparison, we also present the corresponding curves as obtained by increasing the bending constant to k θ = 50 .

2.2. Simulation Algorithm

Simulations are carried out by means of the Simu-D simulator-descriptor suite [90]. The simulations are performed under constant volume and temperature ( T = 1 / k B ), employing Periodic Boundary Conditions (PBCs). Equilibrium simulations are performed mainly at high packing densities in the range φ = 0.54 , 0.60 . Additionally, simulations are also carried out at significantly lower volume fractions to study the long-range orientational order, especially for rod-like chains ( θ 0 = 0 ° ). The choice of volume fractions is driven by the fact that, in past studies, we showed that crystallization for systems of freely-jointed chains of tangent hard spheres occurs once a packing density of φ 0.58 is reached while at φ = 0.56 polymer packings remain amorphous [40,41,42]. For comparison, monomeric analogs crystallize at φ m o n o m e r s M = 0.545 . In Figure S1 (Supporting Information), computer-generated configurations are visualized corresponding to dilute conditions ( φ = 0.10) for all semi-flexible systems studied here.
Systems are generated and equilibrated employing a mix of localized moves and chain-connectivity-altering moves (CCAMs), as implemented in the Simu-D simulator [90]. The MC mix is composed of the following moves: reptation (10%), rotation (10%), flip (34,8%), intermolecular reptation (25%), end-segment re-arrangement (or CCB as in Ref. [97,98]) (20%), simplified end-bridging (sEB) (0.1%), and simplified intramolecular end-bridging (sIEB) (0.1%). The attempt probabilities of each MC move are reported in parentheses and remain unaltered for all systems. Due to the high packing densities, all local moves are executed according to a configurational bias (CB) pattern (see more technical explanations in Refs. [90,91]), with the number of trials per attempt depending on the packing fraction according to Table 1 of Ref. [91].
Initial configurations are generated at very dilute conditions ( φ = 0.001 ) as fully flexible systems. Under these conditions the desired equilibrium bending angle, θ 0 , and bending constant, k θ are activated, followed by rapid constant-volume equilibration. The dilute systems of semi-flexible chains are then shrunk until the desired volume fraction is reached through isotropic volume compressions attempted at regular intervals [91]. Equilibrium simulations are conducted at constant volume for every packing fraction between a minimum of 5 × 10 11 and a maximum of 8 × 10 11 MC steps, depending on the system. Systems configurations (“frames”) are recorded every 10 7 MC steps.

2.3. Global (Long-Range) Orientational Order

As we simulate semi-flexible chains, it is important to be able to quantify the long-range orientational order due to the self-organization of chains [73]. In nematic order, chains exhibit a long-range orientation with their long axis aligned along a preferent direction, denominated nematic director n. The isotropic-to-nematic transition is influenced by parameters such as packing density [78,79] and ratio of chain size [75,76,77].
Chain orientational order is determined through averages of a second-order invariant [99]. The orientation of each molecule is defined by the unit vector u along the long axis of the molecule. The long axis of each chain is determined by the inertia tensor I, a second-order tensor calculated according to Equation (4), where N(j) is the number of monomers of the chain j, mi is the mass of each monomer (considered here as unity), xi is the coordinate vector of the monomer, xcm is the coordinate vector of the centre of mass of the chain and δ is the second order isotropic tensor. The long axis of each chain corresponds to the eigenvector, v 3 , of the lowest eigenvalue of the inertia tensor I so that the unit vector u is calculated normalizing the eigenvector v 3 .
I = i = 1 N ( j ) m i x i x c m σ 2 δ x i x i
u = v 3 v 3
As an average of all molecular orientations, we define the second-order tensor Q according to the next expression:
Q = 1 N c h i = 1 N c h u i u i 1 3 δ .
The Q tensor presents some characteristic properties: 1) is symmetric, and 2) is traceless (TrQ = 0). These two properties reduce the independent components of the second-order tensor from 9 to 5.
The tensor Q can be calculated for the ideal different phases (isotropic, prolate mesogen, nematic mesophase, and oblate mesogen, nematic mesophase; for their schemes see for example Fig. 4 in [99]). For the isotropic phase, denoted here as “ISO”, the molecules can present any random orientation with equal probability, so no orientation is preferred. Thus, all components of the resulting matrix QISO are zero. In contrast, the perfectly aligned prolate mesogen, nematic mesophase case represents a system where the long axis of every molecule is aligned along the preferred orientation. In a coordinate system in which the preferred orientation of the system is along the x-axis, the matrix representation of Q for the perfect prolate mesogen, nematic mesophase, denoted here as “PRO”, takes the form:
Q P R O = 2 3 0 0 0 1 3 0 0 0 1 3 .
In the oblate mesogen, nematic mesophase, it is the short axis of the molecules that tends to be aligned with the nematic director n. In the same coordinate system, the matrix representation of Q for the oblate mesogen, nematic mesophase, denoted here as “OBL”, takes the form:
Q O B L = 1 3 0 0 0 1 6 0 0 0 1 6 .
In order to identify similarities with these ideal cases, we diagonalize Q and use its normalized eigenvalues as unit vectors of a new Cartesian coordinate system where Q’ is diagonal and its diagonal elements are its eigenvalues (λ1, λ2, and λ3) ordered by absolute value. The eigenvector of the largest absolute eigenvalue λ1, in the case of a nematic mesophase (denoted here as “NEM”), would correspond to the preferred orientation of the system or nematic director n.
A scalar order parameter q is obtained by comparing the diagonal tensor Q’ with the corresponding tensor QPRO of a perfect prolate mesogen, nematic mesophase. This scalar q represents the degree of alignment. In an isotropic system, where Q is null, q → 0. For a perfectly aligned prolate mesogen, nematic mesophase system, q = 1 , while q = -1/2 for the ideal oblate mesogen, nematic mesophase. The orientational order parameter, thus, takes values in the range 1 / 2 q 1 .
Q = λ 1 0 0 0 λ 2 0 0 0 λ 3 = q 2 3 0 0 0 1 3 0 0 0 1 3   λ 1 > λ 2 λ 3
Then, the long-range orientational order is characterized by the scalar orientational order parameter q and the nematic director n.

2.4. Local Order: Characteristic Crystallographic Element Norm

The metric used in this work to identify the disorder-order transition at the local level and quantify the degree of crystallinity of the simulated systems is the Characteristic Crystallographic Element (CCE) norm, which is explained in detail in Refs. [95,96]. The CCE norm is integrated into the descriptor part of the Simu-D suite, which is also used for the MC simulations [90].
The CCE norm descriptor gauges the crystallinity of an atomic or particulate system in two (2D) or three (3D) dimensions by comparing the local environment around each site with a set of ideal reference crystals under the main concept that each ideal crystal is uniquely identified by a set of symmetry operations [100,101,102,103]. So, for every site i of a system, the CCE norm descriptor identifies the nearest neighbors and quantifies the orientational and radial deviations of the “real” local environment with respect to the “ideal” environment of each reference crystal. This comparison provides, for the given site i, a CCE norm value with respect to each X reference crystal, ε i X . The closer the X-CCE norm is to zero, the higher the similarity of the local environment to the respective reference crystal X. Site i is identified as an X-type crystal when the calculated CCE norm is lower than a critical threshold, ε i X ε t h r e s . In past studies, an empirical threshold of 0.245 was determined for packings of non-overlapping spheres [41,95,96] in the bulk [40,41,42,43,55,56] and under confinement [44,45]. Due to the strict concept behind the descriptor, the CCE norm is highly discriminatory, so the value of a site cannot be simultaneously very low for two different reference crystals. The current version of the CCE norm descriptor can identify similarity with respect to the following crystals: hexagonal close packed (HCP), face centered cubic (FCC), body centered cubic (BCC), and hexagonal (HEX) for 3D systems and triangular (TRI), square (SQU), and honeycomb (HON) for 2D systems. The fivefold (FIV) and pentagonal (PEN) local symmetries can be identified as well in three and two dimensions, respectively. Sites that cannot be assigned to any of the previous reference ideal environments are labeled as amorphous (AMO) or, more precisely, as “unidentified”.
The process explained before is repeated over all sites of the system for each reference crystal. Once the CCE norm has been evaluated for every site and reference crystal, an order parameter for each reference crystal X, S X ( S X 0,1 ), can be calculated for the snapshot as
S X = 0 ε t h r e s P ε X d ε X
where P ( ε X ) is the probability function of norms of the reference crystal X.
Additionally, a degree of crystallinity or total crystallinity, τ c , can be calculated as the sum of the order parameters of all reference crystals. For a bulk, 3D system, the degree of crystallinity is calculated as:
τ c = k = 1 4 0 ε t h r e s P ε X d ε X = k = 1 4 S k = S H C P + S F C C + S B C C + S H E X .
In the present work, no appreciable population of sites with HEX or BCC similarity was detected (SHEX, SBCC → 0). Thus, in the continuation, we consider only the HCP and FCC crystals as well as fivefold local symmetry (FIV).

3. Results

The phase behavior of polymer chains under the effect of the chain stiffness is the main focus of this work, emphasizing the long-range, orientational (nematic) ordering and the local structure.
Prior to simulation data being analyzed in a post-processing step, a preliminary visual inspection of the initial and resulting system configurations is performed. In the case of rod-like chains ( θ 0 = 0 ° ), the visual inspection of the initial configurations, as presented in Figure 3, suggests a transition from an isotropic (φ = 0.10 and 0.20) to a nematic phase (φ = 0.30 and 0.50) with increasing packing density, i.e. the rod-like chains align in the system along the common nematic director. This transition happens at a packing density that is significantly lower than for solidification, in perfect qualitative agreement with past independent works on similar systems [73,78,79]. Figure S2 (in Supporting Information) shows computer-generated snapshots of bulk system configurations at the end of the simulation for semi-flexible hard-sphere systems at φ = 0.60.

3.1. Global Orientational Order

Following the preliminary visual inspection of the initial configurations for the rod-like chains ( θ 0 = 0 ° ), the long-range, orientational (nematic) ordering is analyzed through the second-order tensor Q and the scalar order parameter q, as explained in the previous section. The dependence of the nematic order parameter q on packing density, φ, is presented in Figure 4 for all semi-flexible systems under study ( θ 0 = 0, 60, 90, 108, and 120°), averaged over all frames of the equilibrated part of the simulation trajectory. In line with Figure 3, rod-like chains of average length Nav = 12 show a well-defined isotropic-nematic transition as concentration increases. This transition takes place in the interval 0.15 ≤ φ ≤ 0.20 and very closely resembles the prediction of Onsager’s theory. At higher volume fractions, chains form a prolate mesogen, nematic mesophase (PRO), as seen by the monotonically increasing value of q. A perfect PRO phase (q → 1) is established at packing densities approximately equal to φ = 0.45. For the average lengths studied here (Nav = 12), the transition from the isotropic to the nematic phase, ISO → PRO (i.e. nematic mesophase with prolate mesogens) for rod-like athermal chains occurs at significantly lower packing densities (φnem ≈ 0.20) than the freezing ( φ m o n o m e r s F 0.495 ) and melting ( φ m o n o m e r s M 0.545 ) points of monomeric hard spheres.
Besides this expected isotropic-to-nematic transition exhibited by rod-like chains, a nematic mesophase with oblate mesogens, OBL, (q < 0) appears for semi-flexible chains with θ 0 = 90°. In contrast to the long-range orientational order of rod-like chains, the isotropic-to-oblate, ISO → OBL, transition takes place at φobl ≈ 0.57, which is higher than the melting point of HS.
For the remaining values of the equilibrium bending angles, each chain system remains in the isotropic phase in the whole concentration range. Still, a small peak is produced at very high packing densities that can be associated with chains primarily of θ 0 = 60°, taking some preferred directions, corresponding to values of the nematic order parameter of around q ≈ 0.15. This can be explained as the establishment of the FCC and HCP crystallites (see later discussion on local order) enforces specific bending angles which are compatible with these crystals. If such bending angles are not available through intra-chain arrangements, then the only other option is inter-chains ones, thus inducing alignment among chains at a local level.
In the case of semi-flexible chains with θ 0 = 90 ° , the left panel of Figure 5 hosts the evolution of the exponential running average of the orientational order parameter, q, as a function of MC steps for the packing densities where a certain degree of nematic order was observed in Figure 4. A tendency for the formation of an oblate mesogen, nematic mesophase (OBL) is observed after sufficient simulation time. The degree of ordering of the oblate mesogens increases with volume fraction as indicated by q approaching the ideal value 1 / 2 . Compared to crystallization at the level of monomers (see below), the ISO → OBL transition for θ 0 = 90° takes places with an equal or even slower rate. We should note here that “rate” used in the present manuscript has no physical meaning and it corresponds to the number of MC steps required to observe the phase transition. Thus, the nematic phase for the semi-flexible chains with θ 0 = 90° practically coincides with crystallization, i.e. self-organization at the local and global levels appears synchronized.
The right panels of Figure 5 show semi-flexible chains with θ 0 = 90° at φ = 0.59 at the end of the MC simulation, once the system has reached an OBL state. The top snapshot presents the unwrapped representation of the semi-flexible chains of the complete system. As in the previous snapshots, monomers are color-coded according to the parent chain. For clarity, the bottom snapshot contains only 10 randomly selected chains of the same system. The semi-flexible chains tend to form flat layers, interrupted occasionally by right-angle jumps between planes, consistent with θ 0 = 90°. These flat chain configurations have, in general, a common behavior that further explains the OBL phase in Figure 4 and the left panel of Figure 5.

3.2. Local Order: Crystal Nucleation and Growth

Local structure is gauged through the CCE norm descriptor [95,96]. For consistency with our past studies on flexible chains of hard spheres, we employ a threshold of ε t h r e s = 0.245 to label a site as X-type, where X is the reference crystal or local symmetry. In the present work given that no appreciable population of BCC or HEX sites is detected in any of the simulated systems, X corresponds to HCP, FCC, and FIV. In the continuation and throughout the manuscript, the corresponding colors to be used for the representation of the HCP, FCC, and FIV sites (in snapshots) and curves (in figures) are blue, red, and green, respectively. The results from the semi-flexible systems are also compared with the ones of fully flexible (freely-jointed) chains of tangent hard spheres simulated and analyzed through the Simu-D software under the same conditions of volume fraction and average chain length [40,41,42,93,94].
Figure 6 presents the CCE-based order parameters for HCP ( S H C P ), FCC ( S F C C ), and FIV ( S F I V ), along with the degree of crystallinity, τ c , as a function of MC steps at a packing density φ = 0.59 for all semi-flexible systems including the fully flexible one as reference. This volume fraction is higher than the melting transition of freely-jointed hard-sphere chains as established in [37,40,41,42,54]. An inspection of all panels shows clearly that the mechanism of phase transition, here in the form of crystallization, is very similar and rather independent of the equilibrium bending angle. Initial packings are amorphous (AMO) as the fraction of sites with non-ordered local structure is vastly dominant: the percentage of close packed sites does not exceed 5% and, in most cases, the FIV population is commensurate or even exceeds the combined HCP+FCC one. This initial state of the AMO athermal chain packings agrees with past works on dense monomeric and polymer assemblies [42,45,55,56,57]. Out of all systems, the rod-like one is characterized by the lowest population of FIV sites at the beginning as FIV appears to be incompatible with the perfect nematic ordering exhibited by the chains at this range of volume fractions.
Once crystallites of HCP and FCC characters start growing, at almost identical rates, the fraction of FIV-like sites decreases gradually until it practically disappears. The structural competition, observed here for all semi-flexible chain systems, between the FIV local symmetry and crystallization in the form of HCP and FCC sites is in perfect match with identical observations in simulations of fully flexible [37,40,41,42,43,44,45] and monomeric HS systems [55,56], including bead-spring chains under quenching [61,105]. Once the population of FIV sites is eliminated, a trend manifestly valid for all simulated packings, the relative fractions of the HCP and FCC crystals undergo a sharp variation which is followed by the establishment of the final, stable ordered structures of crystalline (CRY) character. The average degree of crystallinity in the final CRY phase ranges between 0.65 and 0.85 with the lowest and highest values corresponding to the fully flexible and rod-like chains, respectively. At φ = 0.59, as seen in Figure 6, all crystals contain appreciable fractions of FCC and HCP sites and no perfection is registered. Accordingly, one expects that the formed morphologies correspond to fivefold-free but defect-ridden RHCP crystals.
Figure 7 shows, for all systems whose phase transition is presented in Figure 6, the snapshots at the beginning (top panels) and the end (bottom panels) of the corresponding MC simulations. As already quantified by the data in Figure 6, all initial configurations show amorphous structures with a remarkable population of FIV-like sites. At the end of MC simulations, the crystalline phase of every system, both freely-jointed and semi-flexible ones, shows a stable configuration of mixed HCP/FCC structures, being defect-ridden and fivefold-free, exactly as expected by the fractions of sites in Figure 6. In their majority, the semi-flexible chain packings form RHCP structures of alternating layers of unique FCC and HCP character with a single stacking direction, just like the freely-jointed systems. Additionally, semi-flexible systems can also form crystal structures with multiple, random stacking directions, as it is observed for the system with θ 0 = 90 ° .
The formation of RHCP crystals of mixed HCP/FCC structures, with unique or multiple stacking directions, is also observed for all the equilibrium bending angles studied at all packing densities where crystallization takes place.
Although the majority of the semi-flexible systems crystallize in RHCP crystals of mixed HCP/FCC layers, two important exceptions exist for the rod-like chains ( θ 0 = 0 ° ) at φ = 0.58 and 0.60. Figure 8 shows the CCE-based local order parameter, S X , and the total crystallinity, τ c , as a function of MC steps, while the corresponding snapshots at various simulation instances can be found in Figure 9. Both systems show very similar trends: after a very short initial period characterized by the rapid increase of HCP and FCC sites and the parallel reduction of the FIV population, the growth rates stop being the same and one type grows in favor of the other. At φ = 0.58, the resulting morphology is an almost perfect FCC crystal ( τ c = 0.83 , S F C C = 0.83 ) while the opposite happens at φ = 0.60 where a (less) perfect HCP crystal emanates ( τ c = 0.77 , S H C P = 0.74 ). In the case of the HCP crystal at φ = 0.60, the resulting structure is not so stable as the one of the FCC crystal, alternating between the almost perfect HCP crystal ( S F C C 0 ) and an HCP crystal with a small population of FCC-like sites ( S F C C 0.1 ). For the latter case where FCC impurities appear in the HCP crystal, these are produced in the border of the crystal with a regime of defects in the form of amorphous (AMO) sites. Figure S4 (in Supporting Information) hosts the CCE-based snapshots of the final configuration for the corresponding systems.
While perfection in the form of FCC crystal has been reported very recently from extremely long simulations on freely-jointed chains of hard spheres [37], it is the first time that an almost perfect HCP crystal, made of hard-sphere chains, in extended rod-like conformations, is observed. The easiness with which crystal perfection is observed could be related to the nematic ordering exhibited by the rod-like chains. This trend, especially compared to the fully flexible model, will be explored in more detail in future studies.

3.3. Local Order: Total Crystallinity

From the local order parameters S X , the total degree of crystallinity, τ c , can be gauged for all the simulated semi-flexible systems. As in the previous analysis, present results are compared against the ones of fully flexible (freely-jointed) analogs. Figure S5 (Supporting Information) presents the evolution of the degree of crystallinity as a function of MC steps at different packing densities, φ.
Based on these results and utilizing the value of the total crystallinity as established in the final, stable part of the MC simulation, we can further extract the one-dimensional phase diagram of local order, quantified by the degree of crystallinity, as a function of packing density and equilibrium bending angle (Figure 10). We remind here that the melting point for monomers is φ m o n o m e r s M = 0.545 while freely-jointed chains of tangent hard spheres crystallize at higher volume fractions, φ c h a i n s M (> 0.56). Some important conclusions can be drawn from the diagram. First, all semi-flexible chain packings eventually crystallize, independently of the equilibrium angle. Second, the equilibrium angle affects profoundly the phase transition of the systems. For example, the rod-like chains present semi-crystalline and crystalline phases at volume fractions significantly lower than the melting point for freely-jointed chains, even lower than the one for monomeric hard spheres, due to the effect of the nematic (global) ordering that precedes local order (as will be demonstrated in the continuation). On the other hand, chains with θ 0 = 60 ° , crystallize later than all other systems as at φ = 0.57 no crystal nucleation and growth is observed, even after 7 × 10 11 MC steps. Finally, the melting transition shows the trend: φ c h a i n s M 0 ° < φ m o n o m e r s M = 0.545 < φ c h a i n s M 120 ° < φ c h a i n s M 108 ° < φ c h a i n s M 90 ° = φ c h a i n s M F J < φ c h a i n s M 60 ° .
In general, with the sole exception of the rod-like chains, the presence of constraints, related to bond geometry in the form of bonds or bending angles, increases the melting transition with respect to monomeric analogs. Still, it is quite surprising that the systems with obtuse bending angles ( θ 0 = 120   a n d   108 ° ) crystallize earlier and the right angle ( θ 0 = 90 ° ) at a very similar volume fraction as the fully flexible chains, which are free of bending constraints. Obtuse angles are favored at high volume fractions as they minimize the local volume compared to acute ones. However, one should further consider the fact that only specific geometric arrangements of polymer chains are compatible with the sites of ideal FCC and HCP crystals. For example, the FCC and HCP share the bending angles of 0, 60, 90, and 120 ° , while the angles at 33.5 and 70.5 ° exist exclusively on the HCP crystal [37]. The lack of bending angles compatible with these crystals would require the adjacent sites to be covered by monomers belonging to other chains, imposing specific arrangements at the intermolecular level.
The phase behavior at the local (CCE-based crystallinity, τ c ) and global (nematic orientational order parameter, q) levels as a function of density for rod-like chains and right-angle chains is presented in Figure 11. As was explained previously and based on the data of Figure 4, rod-like chains present an isotropic-nematic transition at φnem ≈ 0.20, a volume fraction significantly lower than the freezing point for monomeric hard spheres ( φ m o n o m e r s F = 0.494 ), in accordance with the trends presented in [78,79]. At volume fractions close to the freezing point, the rod-like systems reach a practically perfect PRO phase. According to Shakirov and Paul [63], the rod-like chains self-arrange into a 2D hexagonal crystal structure in the plane perpendicular to the nematic director at a volume fraction φ = 0.47, driven by 2D translational entropy. Thus, those 2D hexagonal crystal structures evolve into a 3D semi-crystalline phase of RHCP structures with the increase in concentration in the range between the freezing and melting point for monomeric hard spheres. After reaching the melting point for monomeric hard spheres, the nematic systems start to crystallize in HCP, FCC, or mixed HCP/FCC structures, keeping the nematic order. So, rod-like systems pack at high densities into Nematic Close Packed (NCP) structures, as was observed by independent researchers through MD simulations [61,62].
Based on the local and global phase behavior, the one-dimensional diagram as a function of packing density can be split into three distinct regions, marked by XX-YY, where the first index XX corresponds to the local structure (XX = AMO or CRY) and the second index YY to the global structure (YY = ISO, PRO or OBL). Accordingly, for the rod-like chains studied here (Nav = 12, k θ = 9 ), we have the following approximate domains of phase behavior (left panel of Figure 11): i) AMO-ISO (φ ≤ 0.15), where the system is amorphous at the local level and isotropic at the global, ii) AMO-PRO (0.20 ≤ φ ≤ 0.45), where the packing is amorphous locally and nematic globally, and iii) CRY-PRO (0.50 ≤ φ) where the system shows crystallinity in the form of HCP and FCC close packed morphologies of varied level of perfection, and perfect nematic ordering at the global orientational level. Comparing the local and global trends, a hysteresis is clearly evident as local order is established at volume fractions significantly higher than the ones for the ISO → PRO transition.
In the case of right-angle chains ( θ 0 = 90 ° ), the phase behavior is quite different compared to rod-like chains as seen in the right panel of Figure 11. The small decrease in the long-range orientational order can be attributed to higher statistical uncertainty as the higher the packing density, and the closer to jamming, the more difficult is the sampling even when such advanced MC algorithms are employed. The right-angle systems do not crystallize until after reaching volume fractions higher than the melting point for monomeric hard spheres ( φ m o n o m e r s M = 0.545 ). Then, the ISO → OBL transition occurs simultaneously to crystal nucleation and growth. The formation of RHCP structures requires a semi-nematic OBL phase for the right-angle systems as the right-angle chain exists but is not dominant in the FCC and HCP crystals. Accordingly, for their creation, specific chain alignments are required. As a result, the local and global orders are simultaneously established and no hysteresis is observed. The phase diagram consists of two regimes: i) AMO-ISO (φ ≤ 0.56), where the system is amorphous at the local level and isotropic at the global, and ii) CRY-OBL (φ > 0.56) where the system shows RHCP morphologies and nematic order of oblate mesogens.

4. Discussion and Conclusions

We present results from extensive Monte Carlo simulations on the phase behavior of semi-flexible chains of tangent hard spheres as a function of packing density and equilibrium bending angle keeping temperature, spring bending constant, and average chain length constant. Local structure is quantified through the Characteristic Crystallographic Element (CCE) norm, while global structure is gauged through the nematic order parameter. A rich one-dimensional phase diagram as a function of packing density is identified where chains crystallize in close-packed morphologies, including random hexagonal close (RHCP) ones of single or multiple stacking directions, or almost perfect HCP and FCC crystals in the case of rod-like chains. Analysis of the long-range orientational tensor reveals the formation of prolate mesogen, nematic mesophase (PRO) for rod-like chains at rather low volume fractions and of oblate mesogen, nematic mesophase (OBL) at high packing densities. Although all systems of semi-flexible chains crystallize, the equilibrium bending angle significantly affects the melting point. While equilibrium angles of 108° and 120° degrees favor crystallization compared to the freely-jointed model, chains with 90° show a behavior that almost coincides with the fully flexible chains, and the acute angle of 60° hinders crystallization enforcing nucleation and growth to take place at higher concentrations.
Particularly, for rod-like chains, three distinct regimes can be identified in the one-dimensional phase diagram where hysteresis is observed between local and global self-organization: i) AMO-ISO (φ ≤ 0.15), where the system is amorphous at the local level and isotropic at the global, ii) AMO-PRO (0.20 ≤ φ ≤ 0.45), where the packing is amorphous locally and nematic globally, and iii) CRY-PRO (0.50 ≤ φ), where the system shows crystallinity in the form of HCP and FCC close packed morphologies of varied level of perfection and perfect nematic ordering at the global orientational level. Right-angle systems show a synchronous establishment of long-range nematic orientational order and formation of RHCP crystallites, thus splitting the behavior into two distinct regimes: i) AMO-ISO (φ ≤ 0.56), where the system is amorphous at the local level and isotropic at the global, and ii) CRY-OBL (φ > 0.56).
Present simulations are currently expanded to treat semi-flexible chains of tangent hard spheres in composites with nanofillers, under confinement, and in mixtures with different species in the form of linear chains and monomeric counterparts.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Figure S1: Snapshots of systems of semi-flexible chains of tangent hard spheres (Nch = 100, N = 12) at a packing density φ = 0.10 with a bending constant kb = 9 and different equilibrium bending angles, θ 0 ; Figure S2: Snapshots of the final configurations from MC simulations on semi-flexible chains of hard spheres at φ = 0.60 and different equilibrium bending angles, θ 0 ; Figure S3: Exponential moving average of the orientational order parameter, q, as a function of MC steps for 100-chains Nav = 12 semi-flexible system with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0 ° at different packing densities, φ; Figure S4: Final configuration for the 100-chains Nav = 12 semi-flexible system chains with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0 ° at φ = 0.58 and 0.60; Figure S5: Degree of crystallinity, τ c , as a function of MC steps at different packing densities, φ.

Author Contributions

Conceptualization, N.C.K.; methodology, N.C.K., K.F., D.M.F. and M.L; software (MC), D.M.F., M.H., N.C.K. and M.L.; software (CCE), K.F., N.C.K., M.H.; data curation, D.M.F.; writing—original draft preparation, D.M.F. and N.C.K.; writing—review and editing, M.H., K.F. and M.L.; visualization, D.M.F.; funding acquisition, K.F., N.C.K. and M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by MICINN/FEDER (Ministerio de Ciencia, Innovación y Universidades, Fondo Europeo de Desarrollo Regional), grant numbers “RTI2018-097338-B-I00” and “PID2021-127533NB-I00” and by UPM and Santander Bank, “Programa Propio UPM Santander”.

Acknowledgments

D.M.F. deeply appreciates the kind hospitality of the research group of Prof. M. Tripathy (IIT, Bombay, India). The authors acknowledge support through projects “RTI 2018-097338-B-I00” and “PID2021-127533NB-I00” of MICINN/FEDER (Ministerio de Ciencia e Innovación, Fondo Europeo de Desarrollo Regional). M.H. and D.M.F. acknowledge financial support through “Programa Propio UPM Santander” of Universidad Politécnica de Madrid and Santander Bank. The authors gratefully acknowledge the Universidad Politécnica de Madrid for providing computing resources on Magerit Supercomputer through projects “p208”, “r553”, “s341”, “u242” and “t736”.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used throughout the manuscript:
2D Two dimensions
3D Three dimensions
AMO Amorphous
BCC Body Centered Cubic
CB Configurational Bias
CCAM Chain-Connectivity-Altering Move
CCB End-segment re-arrangement move
CCE Characteristic Crystallographic Element (norm)
FCC Face Centered Cubic
FIV Fivefold
HCP Hexagonal Close Packed
HEX Hexagonal
HON Honeycomb
HS Hard Sphere
ISO Isotropic phase
MC Monte Carlo
MD Molecular Dynamics
NCP Nematic Close Packed
NEM Nematic mesophase
OBL Oblate mesogen, Nematic mesophase
PBC Periodic Boundary Condition
PEN Pentagonal
PRO Prolate mesogen, Nematic mesophase
RHCP Random Hexagonal Close Packing
SQU Square
TRI Triangular

References

  1. Flory, P.J. Statistical Mechanics of Chain Molecules Hanser-Verlag: Munchen, 1989.
  2. deGennes, P.G. Scaling Concepts in Polymer Physics; Cornell University Press: Ithaca, 1980.
  3. Mullin, J.W. 8 - Industrial Techniques and Equipment. In Crystallization (Fourth Edition), Mullin, J.W., Ed.; Butterworth-Heinemann: Oxford, 2001; pp. 315-402.
  4. Reiter, G.; Sommer, J.U. Polymer Crystallization: Obervations, Concepts and Interpretations; Springer Berlin Heidelberg: 2008.
  5. Zhou, J.H.; Yang, Y.S.; Yang, Y.; Kim, D.S.; Yuan, A.; Tian, X.Z.; Ophus, C.; Sun, F.; Schmid, A.K.; Nathanson, M.; et al. Observing Crystal Nucleation in Four Dimensions Using Atomic Electron Tomography. Nature 2019, 570, 500-+. [CrossRef]
  6. Hu, Y.C.; Tanaka, H. Physical Origin of Glass Formation from Multicomponent Systems. Science Advances 2020, 6. [CrossRef]
  7. Wang, S.Z.; Park, S.S.; Buru, C.T.; Lin, H.X.; Chen, P.C.; Roth, E.W.; Farha, O.K.; Mirkin, C.A. Colloidal Crystal Engineering with Metal-Organic Framework Nanoparticles and DNA. Nat. Commun. 2020, 11. [CrossRef]
  8. Dijkstra, M.; Luijten, E. From Predictive Modelling to Machine Learning and Reverse Engineering of Colloidal Self-Assembly. Nat. Mater. 2021, 20, 762-773. [CrossRef]
  9. Zhang, W.; Mazzarello, R.; Wuttig, M.; Ma, E. Designing Crystallization in Phase-Change Materials for Universal Memory and Neuro-Inspired Computing. Nature Reviews Materials 2019, 4, 150-168. [CrossRef]
  10. Li, M.M.; Mangalore, D.K.; Zhao, J.B.; Carpenter, J.H.; Yan, H.P.; Ade, H.; Yan, H.; Mullen, K.; Blom, P.W.M.; Pisula, W.; et al. Integrated Circuits Based on Conjugated Polymer Monolayer. Nat. Commun. 2018, 9. [CrossRef]
  11. Xia, X.X.; Lau, T.K.; Guo, X.Y.; Li, Y.H.; Qin, M.C.; Liu, K.; Chen, Z.; Zhan, X.Z.; Xiao, Y.Q.; Chan, P.F.; et al. Uncovering the out-of-Plane Nanomorphology of Organic Photovoltaic Bulk Heterojunction by Gtsaxs. Nat. Commun. 2021, 12. [CrossRef]
  12. Alder, B.J.; Wainwright, T.E. Phase Transition for a Hard Sphere System. J. Chem. Phys. 1957, 27, 1208-1209. [CrossRef]
  13. Frenkel, D. Entropy-Driven Phase Transitions. Physica A 1999, 263, 26-38. [CrossRef]
  14. Pusey, P.N.; Vanmegen, W. Phase-Behavior of Concentrated Suspensions of Nearly Hard Colloidal Spheres. Nature 1986, 320, 340-342. [CrossRef]
  15. Miskin, M.Z.; Jaeger, H.M. Adapting Granular Materials through Artificial Evolution. Nat. Mater. 2013, 12, 326-331. [CrossRef]
  16. Athanassiadis, A.G.; Miskin, M.Z.; Kaplan, P.; Rodenberg, N.; Lee, S.H.; Merritt, J.; Brown, E.; Amend, J.; Lipson, H.; Jaeger, H.M. Particle Shape Effects on the Stress Response of Granular Packings. Soft Matter 2014, 10, 48-59. [CrossRef]
  17. Miskin, M.Z.; Jaeger, H.M. Evolving Design Rules for the Inverse Granular Packing Problem. Soft Matter 2014, 10, 3708-3715. [CrossRef]
  18. Zou, L.-N.; Cheng, X.; Rivers, M.L.; Jaeger, H.M.; Nagel, S.R. The Packing of Granular Polymer Chains. Science 2009, 326, 408-410. [CrossRef]
  19. Sacanna, S.; Irvine, W.T.M.; Chaikin, P.M.; Pine, D.J. Lock and Key Colloids. Nature 2010, 464, 575-578. [CrossRef]
  20. Vutukuri, H.R.; Demirors, A.F.; Peng, B.; van Oostrum, P.D.J.; Imhof, A.; van Blaaderen, A. Colloidal Analogues of Charged and Uncharged Polymer Chains with Tunable Stiffness. Angew. Chem. Int. Ed. 2012, 51, 11249-11253. [CrossRef]
  21. Brown, E.; Nasto, A.; Athanassiadis, A.G.; Jaeger, H.M. Strain Stiffening in Random Packings of Entangled Granular Chains. Phys. Rev. Lett. 2012, 108. [CrossRef]
  22. Feng, L.; Pontani, L.L.; Dreyfus, R.; Chaikin, P.; Brujic, J. Specificity, Flexibility and Valence of DNA Bonds Guide Emulsion Architecture. Soft Matter 2013, 9, 9816-9823. [CrossRef]
  23. Allen, M.P.; Tildesley, D.J. Computer Simulation of Liquids, 2nd ed.; Oxford University Press: Oxford (UK), 2017.
  24. Frenkel, D.; Smit, B. Understanding Molecular Simulation: From Algorithms to Applications, second ed.; Academic Press: San Diego, 2002.
  25. Landau, D.P.; Binder, K. A Guide to Monte Carlo Simulations in Statistical Physics, 4th Edition; Cambridge University Press: Cambridge, 2014.
  26. Rapaport, D.C. The Art of Molecular Dynamics Simulation, 2nd Edition; Cambridge University Press: Cambridge, 2004.
  27. Leach, A. Molecular Modelling: Principles and Applications, 2nd Edition; Pearson: 2001.
  28. Wood, W.W.; Jacobson, J.D. Preliminary Results from a Recalculation of the Monte Carlo Equation of State of Hard Spheres. J. Chem. Phys. 1957, 27, 1207-1208. [CrossRef]
  29. Hoover, W.G.; Ree, F.H. Melting Transition and Communal Entropy for Hard Spheres. J. Chem. Phys. 1968, 49, 3609-&. [CrossRef]
  30. Rintoul, M.D.; Torquato, S. Metastability and Crystallization in Hard-Sphere Systems. Phys. Rev. Lett. 1996, 77, 4198-4201. [CrossRef]
  31. Rintoul, M.D.; Torquato, S. Computer Simulations of Dense Hard-Sphere Systems. J. Chem. Phys. 1996, 105, 9258-9265. [CrossRef]
  32. Bolhuis, P.G.; Frenkel, D.; Mau, S.C.; Huse, D.A. Entropy Difference between Crystal Phases. Nature 1997, 388, 235-236. [CrossRef]
  33. Woodcock, L.V. Entropy Difference between the Face-Centred Cubic and Hexagonal Close-Packed Crystal Structures. Nature 1997, 385, 141-143. [CrossRef]
  34. Bruce, A.D.; Wilding, N.B.; Ackland, G.J. Free Energy of Crystalline Solids: A Lattice-Switch Monte Carlo Method. Phys. Rev. Lett. 1997, 79, 3002-3005. [CrossRef]
  35. Auer, S.; Frenkel, D. Prediction of Absolute Crystal-Nucleation Rate in Hard-Sphere Colloids. Nature 2001, 409, 1020-1023. [CrossRef]
  36. Luchnikov, V.; Gervois, A.; Richard, P.; Oger, L.; Troadec, J.P. Crystallization of Dense Hard Sphere Packings - Competition of Hcp and Fcc Close Order. J. Mol. Liq. 2002, 96-7, 185-194. [CrossRef]
  37. Herranz, M.; Foteinopoulou, K.; Karayiannis, N.C.; Laso, M. Polymorphism and Perfection in Crystallization of Hard Sphere Polymers. Polymers 2022, 14. [CrossRef]
  38. Espinosa, J.R.; Vega, C.; Valeriani, C.; Frenkel, D.; Sanz, E. Heterogeneous Versus Homogeneous Crystal Nucleation of Hard Spheres. Soft Matter 2019, 15, 9625-9631. [CrossRef]
  39. Vega, C.; McBride, C.; MacDowell, L.G. Liquid Crystal Phase Formation for the Linear Tangent Hard Sphere Model from Monte Carlo Simulations. J. Chem. Phys. 2001, 115, 4203-4211. [CrossRef]
  40. Karayiannis, N.C.; Foteinopoulou, K.; Laso, M. Entropy-Driven Crystallization in Dense Systems of Athermal Chain Molecules. Phys. Rev. Lett. 2009, 103. [CrossRef]
  41. Karayiannis, N.C.; Foteinopoulou, K.; Abrams, C.F.; Laso, M. Modeling of Crystal Nucleation and Growth in Athermal Polymers: Self-Assembly of Layered Nano-Morphologies. Soft Matter 2010, 6, 2160-2173. [CrossRef]
  42. Karayiannis, N.C.; Foteinopoulou, K.; Laso, M. Spontaneous Crystallization in Athermal Polymer Packings. Int. J. Mol. Sci. 2013, 14, 332-358. [CrossRef]
  43. Karayiannis, N.C.; Foteinopoulou, K.; Laso, M. The Role of Bond Tangency and Bond Gap in Hard Sphere Crystallization of Chains. Soft Matter 2015, 11, 1688-1700. [CrossRef]
  44. Ramos, P.M.; Herranz, M.; Foteinopoulou, K.; Karayiannis, N.C.; Laso, M. Entropy-Driven Heterogeneous Crystallization of Hard-Sphere Chains under Unidimensional Confinement. Polymers 2021, 13. [CrossRef]
  45. Ramos, P.M.; Herranz, M.; Martinez-Fernandez, D.; Foteinopoulou, K.; Laso, M.; Karayiannis, N.C. Crystallization of Flexible Chains of Tangent Hard Spheres under Full Confinement. J. Phys. Chem. B 2022, 126, 5931-5947. [CrossRef]
  46. Sirota, E.B. Polymer Crystallization: Metastable Mesophases and Morphology. Macromolecules 2007, 40, 1043-1048. [CrossRef]
  47. Sirota, E.B.; Singer, D.M. Phase-Transitions among the Rotator Phases of the Normal-Alkanes. J. Chem. Phys. 1994, 101, 10873-10882. [CrossRef]
  48. Sirota, E.B.; King, H.E.; Singer, D.M.; Shao, H.H. Rotator Phases of the Normal Alkanes - an X-Ray-Scattering Study. J. Chem. Phys. 1993, 98, 5809-5824. [CrossRef]
  49. Bourque, A.; Locker, C.R.; Rutledge, G.C. Molecular Dynamics Simulation of Surface Nucleation During Growth of an Alkane Crystal. Macromolecules 2016, 49, 3619-3629. [CrossRef]
  50. Morthomas, J.; Fusco, C.; Zhai, Z.Q.; Lame, O.; Perez, M. Crystallization of Finite-Extensible Nonlinear Elastic Lennard-Jones Coarse-Grained Polymers. Phys. Rev. E 2017, 96. [CrossRef]
  51. Luo, C.F.; Kroger, M.; Sommer, J.U. Entanglements and Crystallization of Concentrated Polymer Solutions: Molecular Dynamics Simulations. Macromolecules 2016, 49, 9017-9025. [CrossRef]
  52. Zhang, W.; Zou, L. Molecular Dynamics Simulations of Crystal Nucleation near Interfaces in Incompatible Polymer Blends. Polymers 2021, 13, 347.
  53. Isobe, M.; Krauth, W. Hard-Sphere Melting and Crystallization with Event-Chain Monte Carlo. J. Chem. Phys. 2015, 143. [CrossRef]
  54. Karayiannis, N.C.; Foteinopoulou, K.; Laso, M. The Structure of Random Packings of Freely Jointed Chains of Tangent Hard Spheres. J. Chem. Phys. 2009, 130. [CrossRef]
  55. Karayiannis, N.C.; Malshe, R.; Kroger, M.; de Pablo, J.J.; Laso, M. Evolution of Fivefold Local Symmetry During Crystal Nucleation and Growth in Dense Hard-Sphere Packings. Soft Matter 2012, 8, 844-858. [CrossRef]
  56. Karayiannis, N.C.; Malshe, R.; de Pablo, J.J.; Laso, M. Fivefold Symmetry as an Inhibitor to Hard-Sphere Crystallization. Phys. Rev. E 2011, 83. [CrossRef]
  57. Karayiannis, N.C.; Foteinopoulou, K.; Laso, M. Twinning of Polymer Crystals Suppressed by Entropy. Symmetry-Basel 2014, 6, 758-780. [CrossRef]
  58. Ni, R.; Dijkstra, M. Effect of Bond Length Fluctuations on Crystal Nucleation of Hard Bead Chains. Soft Matter 2013, 9, 365-369. [CrossRef]
  59. Ramos, P.M.; Herranz, M.; Martínez-Fernández, D.; Foteinopoulou, K.; Laso, M.; Karayiannis, N.C. Crystallization of Flexible Chains of Tangent Hard Spheres under Full Confinement. The Journal of Physical Chemistry B 2022, 126, 5931-5947. [CrossRef]
  60. Luo, C.F.; Kroger, M.; Sommer, J.U. Molecular Dynamics Simulations of Polymer Crystallization under Confinement: Entanglement Effect. Polymer 2017, 109, 71-84. [CrossRef]
  61. Nguyen, H.T.; Smith, T.B.; Hoy, R.S.; Karayiannis, N.C. Effect of Chain Stiffness on the Competition between Crystallization and Glass-Formation in Model Unentangled Polymers. J. Chem. Phys. 2015, 143. [CrossRef]
  62. Nguyen, H.T.; Hoy, R.S. Effect of Chain Stiffness and Temperature on the Dynamics and Microstructure of Crystallizable Bead-Spring Polymer Melts. Phys. Rev. E 2016, 94, 052502. [CrossRef]
  63. Shakirov, T.; Paul, W. Crystallization in Melts of Short, Semiflexible Hard Polymer Chains: An Interplay of Entropies and Dimensions. Phys. Rev. E 2018, 97. [CrossRef]
  64. Shakirov, T. Crystallisation in Melts of Short, Semi-Flexible Hard-Sphere Polymer Chains: The Role of the Non-Bonded Interaction Range. Entropy 2019, 21. [CrossRef]
  65. Li, Y.; Huang, Q.; Shi, T.; An, L. Effects of Chain Flexibility on Polymer Conformation in Dilute Solution Studied by Lattice Monte Carlo Simulation. J. Phys. Chem. B 2006, 110, 23502-23506. [CrossRef]
  66. Lamura, A.; Burkhardt, T.W.; Gompper, G. Semiflexible Polymer in a Uniform Force Field in Two Dimensions. Phys. Rev. E 2001, 64, 061801. [CrossRef]
  67. Hsu, H.-P.; Paul, W.; Binder, K. Polymer Chain Stiffness Vs. Excluded Volume: A Monte Carlo Study of the Crossover Towards the Worm-Like Chain Model. Epl 2010, 92, 28003. [CrossRef]
  68. Cifra, P. Differences and Limits in Estimates of Persistence Length for Semi-Flexible Macromolecules. Polymer 2004, 45, 5995-6002. [CrossRef]
  69. Zierenberg, J.; Marenz, M.; Janke, W. Dilute Semiflexible Polymers with Attraction: Collapse, Folding and Aggregation. Polymers 2016, 8. [CrossRef]
  70. Auhl, R.; Everaers, R.; Grest, G.S.; Kremer, K.; Plimpton, S.J. Equilibration of Long Chain Polymer Melts in Computer Simulations. J. Chem. Phys. 2003, 119, 12718-12728. [CrossRef]
  71. Bulacu, M.; van der Giessen, E. Molecular-Dynamics Simulation Study of the Glass Transition in Amorphous Polymers with Controlled Chain Stiffness. Phys. Rev. E 2007, 76, 011807. [CrossRef]
  72. Plaza-Rivera, C.O.; Nguyen, H.T.; Hoy, R.S. Isostaticity and the Solidification of Semiflexible Polymer Melts. Soft Matter 2017, 13, 7948-7952. [CrossRef]
  73. Onsager, L. The Effects of Shape on the Interaction of Colloidal Particles. Annals of the New York Academy of Sciences 1949, 51, 627-659. [CrossRef]
  74. Bolhuis, P.; Frenkel, D. Tracing the Phase Boundaries of Hard Spherocylinders. J. Chem. Phys. 1997, 106, 666-687. [CrossRef]
  75. Khokhlov, A.R.; Semenov, A.N. On the Theory of Liquid-Crystalline Ordering of Polymer-Chains with Limited Flexibility. Journal of Statistical Physics 1985, 38, 161-182. [CrossRef]
  76. Khokhlov, A.R.; Semenov, A.N. Liquid-Crystalline Ordering in the Solution of Partially Flexible Macromolecules. Physica A 1982, 112, 605-614. [CrossRef]
  77. Khokhlov, A.R.; Semenov, A.N. Liquid-Crystalline Ordering in the Solution of Long Persistent Chains. Physica A 1981, 108, 546-556. [CrossRef]
  78. Yethiraj, A.; Fynewever, H. Isotropic to Nematic Transition in Semiflexible Polymer Melts. Mol. Phys. 1998, 93, 693-701. [CrossRef]
  79. Fynewever, H.; Yethiraj, A. Phase Behavior of Semiflexible Tangent Hard Sphere Chains. J. Chem. Phys. 1998, 108, 1636-1644. [CrossRef]
  80. Jaffer, K.M.; Opps, S.B.; Sullivan, D.E.; Nickel, B.G.; Mederos, L. The Nematic-Isotropic Phase Transition in Semiflexible Fused Hard-Sphere Chain Fluids. J. Chem. Phys. 2001, 114, 3314-3324. [CrossRef]
  81. Ivanov, V.A.; Martemyanova, J.A.; Rodionova, A.S.; Stukan, M.R. Computer Simulation of Stiff-Chain Polymers. Polymer Science Series C 2013, 55, 4-22. [CrossRef]
  82. Tang, X.L.; Chen, W.; Li, L.B. The Tough Journey of Polymer Crystallization: Battling with Chain Flexibility and Connectivity. Macromolecules 2019, 52, 3575-3591. [CrossRef]
  83. Midya, J.; Egorov, S.A.; Binder, K.; Nikoubashman, A. Phase Behavior of Flexible and Semiflexible Polymers in Solvents of Varying Quality. J. Chem. Phys. 2019, 151. [CrossRef]
  84. Egorov, S.A.; Milchev, A.; Virnau, P.; Binder, K. Semiflexible Polymers under Good Solvent Conditions Interacting with Repulsive Walls. J. Chem. Phys. 2016, 144. [CrossRef]
  85. Egorov, S.A.; Milchev, A.; Binder, K. Semiflexible Polymers in the Bulk and Confined by Planar Walls. Polymers 2016, 8. [CrossRef]
  86. Kos, P.I.; Ivanov, V.A.; Chertovich, A.V. Crystallization of Semiflexible Polymers in Melts and Solutions. Soft Matter 2021, 17, 2392-2403. [CrossRef]
  87. Roy, S.; Chen, Y.L. Rich Phase Transitions in Strongly Confined Polymer-Nanoparticle Mixtures: Nematic Ordering, Crystallization, and Liquid-Liquid Phase Separation. J. Chem. Phys. 2021, 154. [CrossRef]
  88. Dietz, J.D.; Hoy, R.S. Two-Stage Athermal Solidification of Semiflexible Polymers and Fibers. Soft Matter 2020, 16, 6206-6217. [CrossRef]
  89. Hoy, R.S. Jamming of Semiflexible Polymers. Phys. Rev. Lett. 2017, 118. [CrossRef]
  90. Herranz, M.; Martínez-Fernández, D.; Ramos, P.M.; Foteinopoulou, K.; Karayiannis, N.C.; Laso, M. Simu-D: A Simulator-Descriptor Suite for Polymer-Based Systems under Extreme Conditions. Int. J. Mol. Sci. 2021, 22, 12464.
  91. Karayiannis, N.C.; Laso, M. Monte Carlo Scheme for Generation and Relaxation of Dense and Nearly Jammed Random Structures of Freely Jointed Hard-Sphere Chains. Macromolecules 2008, 41, 1537-1551. [CrossRef]
  92. Ramos, P.M.; Karayiannis, N.C.; Laso, M. Off-Lattice Simulation Algorithms for Athermal Chain Molecules under Extreme Confinement. J. Comput. Phys. 2018, 375, 918-934. [CrossRef]
  93. Karayiannis, N.C.; Foteinopoulou, K.; Laso, M. Jamming and Crystallization in Athermal Polymer Packings. Philosophical Magazine 2013, 93, 4108-4131. [CrossRef]
  94. Foteinopoulou, K.; Karayiannis, N.C.; Laso, M. Monte Carlo Simulations of Densely-Packed Athermal Polymers in the Bulk and under Confinement. Chem. Eng. Sci. 2015, 121, 118-132. [CrossRef]
  95. Karayiannis, N.C.; Foteinopoulou, K.; Laso, M. The Characteristic Crystallographic Element Norm: A Descriptor of Local Structure in Atomistic and Particulate Systems. J. Chem. Phys. 2009, 130. [CrossRef]
  96. Ramos, P.M.; Herranz, M.; Foteinopoulou, K.; Karayiannis, N.C.; Laso, M. Identification of Local Structure in 2-D and 3-D Atomic Systems through Crystallographic Analysis. Crystals 2020, 10. [CrossRef]
  97. Siepmann, J.I.; Frenkel, D. Configurational Bias Monte-Carlo - a New Sampling Scheme for Flexible Chains. Mol. Phys. 1992, 75, 59-70. [CrossRef]
  98. Laso, M.; de Pablo, J.J.; Suter, U.W. Simulation of Phase-Equilibria for Chain Molecules. J. Chem. Phys. 1992, 97, 2817-2819. [CrossRef]
  99. Andrienko, D. Introduction to Liquid Crystals. J. Mol. Liq. 2018, 267, 520-541. [CrossRef]
  100. Nye, J.F. Physical Properties of Crystals: Their Representation by Tensors and Matrices; Oxford Science Publications: Oxford 2010.
  101. Malgrange, C.; Ricolleau, C.; Schlenker, M. Symmetry and Physical Properties of Crystals; Springer: Dordrecht, 2014.
  102. Giacovazzo, C.; Monaco, H.L.; Artioli, G.; Viterbo, D.; Ferraris, G.; Gilli, G.; Zanotti, G.; Gatti, M. Fundamentals of Crystallography; Oxford Science: Oxford, 2005.
  103. Laso, M.; Jimeno, N. Representation Surfaces for Physical Properties of Materials: A Visual Approach to Understanding Anisotropic Materials; Springer: 2020.
  104. Humphrey, W.; Dalke, A.; Schulten, K. Vmd: Visual Molecular Dynamics. J. Mol. Graphics Modell. 1996, 14, 33-38. [CrossRef]
  105. Hoy, R.S.; Karayiannis, N.C. Simple Model for Chain Packing and Crystallization of Soft Colloidal Polymers. Phys. Rev. E 2013, 88. [CrossRef]
Figure 1. a) Scheme of the bond geometry of a linear sequence of four monomers, indicating bond lengths, l, bending angles, θ, and torsion angles, φ , b) Schematic representation of trimers having equilibrium bending angles ( θ 0 ) studied in the present work.
Figure 1. a) Scheme of the bond geometry of a linear sequence of four monomers, indicating bond lengths, l, bending angles, θ, and torsion angles, φ , b) Schematic representation of trimers having equilibrium bending angles ( θ 0 ) studied in the present work.
Preprints 66761 g001
Figure 2. Distribution of bending angle, θ, of a 100-chain system of Nav = 12 at φ = 0.30 in the bulk for k θ = 9 (solid lines) and 50 (dashed lines) and different equilibrium bending angles, θ 0 .
Figure 2. Distribution of bending angle, θ, of a 100-chain system of Nav = 12 at φ = 0.30 in the bulk for k θ = 9 (solid lines) and 50 (dashed lines) and different equilibrium bending angles, θ 0 .
Preprints 66761 g002
Figure 3. Snapshots of the initial configurations for the 100-chains, Nav = 12, k θ = 9 , θ 0 = 0° system showing at increasing packing densities. From left to right: φ = 0.10, 0.20, 0.30, and 0.50. Monomers are colored according to the parent chain. For clarity, only 50 chains (half of the total population) are shown. Chains are shown with the coordinates being fully unwrapped in space. Images created with VMD visualization software [104].
Figure 3. Snapshots of the initial configurations for the 100-chains, Nav = 12, k θ = 9 , θ 0 = 0° system showing at increasing packing densities. From left to right: φ = 0.10, 0.20, 0.30, and 0.50. Monomers are colored according to the parent chain. For clarity, only 50 chains (half of the total population) are shown. Chains are shown with the coordinates being fully unwrapped in space. Images created with VMD visualization software [104].
Preprints 66761 g003
Figure 4. Evolution of the orientational (nematic) order parameter, q, as a function of the packing density, φ, at different equilibrium bending angles, θ 0 . Also shown for comparison as dashed vertical lines are the freezing, φ m o n o m e r s F 0.495 , and melting, φ m o n o m e r s M 0.545 , points of monomeric hard spheres.
Figure 4. Evolution of the orientational (nematic) order parameter, q, as a function of the packing density, φ, at different equilibrium bending angles, θ 0 . Also shown for comparison as dashed vertical lines are the freezing, φ m o n o m e r s F 0.495 , and melting, φ m o n o m e r s M 0.545 , points of monomeric hard spheres.
Preprints 66761 g004
Figure 5. (Left) Exponential moving average of the orientational order parameter, q, as a function of MC steps 100-chains Nav = 12 semi-flexible system with bending constant k θ = 9 and equilibrium bending angle θ 0 = 90 ° at different packing densities, φ. (Right) Visual representations at φ = 0.59 showing all 100 chains (top) and 10 randomly selected chains (bottom) with coordinates of the monomers being fully unwrapped in space. Monomers are colored according to the parent chain. Images created with VMD visualization software [104].
Figure 5. (Left) Exponential moving average of the orientational order parameter, q, as a function of MC steps 100-chains Nav = 12 semi-flexible system with bending constant k θ = 9 and equilibrium bending angle θ 0 = 90 ° at different packing densities, φ. (Right) Visual representations at φ = 0.59 showing all 100 chains (top) and 10 randomly selected chains (bottom) with coordinates of the monomers being fully unwrapped in space. Monomers are colored according to the parent chain. Images created with VMD visualization software [104].
Preprints 66761 g005
Figure 6. Crystal nucleation and growth, S H C P and S F C C (left y-axis), and fraction of fivefold-like sites, S F I V (right y-axis), as a function of MC steps for 100-chains Nav = 12 systems of freely-jointed chains and semi-flexible chains with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0, 60, 90, 108, and 120° at a packing density φ = 0.59.
Figure 6. Crystal nucleation and growth, S H C P and S F C C (left y-axis), and fraction of fivefold-like sites, S F I V (right y-axis), as a function of MC steps for 100-chains Nav = 12 systems of freely-jointed chains and semi-flexible chains with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0, 60, 90, 108, and 120° at a packing density φ = 0.59.
Preprints 66761 g006
Figure 7. Snapshots of the initial (top) and final configuration (bottom) of the MC simulations for the 100-chains Nav = 12 systems at φ = 0.59 of (from left to right) freely-jointed chains and semiflexible chains ( k θ = 9 ) with equilibrium bending angle θ 0 = 0, 60, 90, 108, and 120°. Monomers are color-coded according to CCE norm. Blue, red, cyan, purple, and green correspond to HCP-, FCC-, BCC-, HEX-, and FIV-like sites, respectively. Amorphous (AMO) sites are colored in yellow and are shown with reduced dimensions for visual clarity.
Figure 7. Snapshots of the initial (top) and final configuration (bottom) of the MC simulations for the 100-chains Nav = 12 systems at φ = 0.59 of (from left to right) freely-jointed chains and semiflexible chains ( k θ = 9 ) with equilibrium bending angle θ 0 = 0, 60, 90, 108, and 120°. Monomers are color-coded according to CCE norm. Blue, red, cyan, purple, and green correspond to HCP-, FCC-, BCC-, HEX-, and FIV-like sites, respectively. Amorphous (AMO) sites are colored in yellow and are shown with reduced dimensions for visual clarity.
Preprints 66761 g007
Figure 8. CCE norm order parameter, S X , as a function of MC steps for the 100-chains Nav = 12 semi-flexible system chains with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0 ° at φ = 0.58 (left panel) and 0.60 (right panel). Inset corresponds to the early part of the MC trajectory (left: 10 × 109 and right: 35 × 109 MC steps). Dashed black line corresponds to the evolution of the degree of crystallinity, τ c , of the system.
Figure 8. CCE norm order parameter, S X , as a function of MC steps for the 100-chains Nav = 12 semi-flexible system chains with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0 ° at φ = 0.58 (left panel) and 0.60 (right panel). Inset corresponds to the early part of the MC trajectory (left: 10 × 109 and right: 35 × 109 MC steps). Dashed black line corresponds to the evolution of the degree of crystallinity, τ c , of the system.
Preprints 66761 g008
Figure 9. Computer-generated representations of the crystal growth of the 100-chains Nav = 12 semi-flexible system chains with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0 ° at (top) φ = 0.58 after (from left to right) 1, 25, 60 and 100   × 10 8 MC steps, and at (bottom) φ = 0.60 after (from left to right) 1, 240, 300 and 350 × 10 8 MC steps. Monomers are color-coded according to CCE norm. Blue, red, and green correspond to HCP-, FCC-, and FIV-like sites, respectively. Amorphous (AMO) sites are colored yellow. For visual clarity, monomers are shown with reduced dimensions.
Figure 9. Computer-generated representations of the crystal growth of the 100-chains Nav = 12 semi-flexible system chains with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0 ° at (top) φ = 0.58 after (from left to right) 1, 25, 60 and 100   × 10 8 MC steps, and at (bottom) φ = 0.60 after (from left to right) 1, 240, 300 and 350 × 10 8 MC steps. Monomers are color-coded according to CCE norm. Blue, red, and green correspond to HCP-, FCC-, and FIV-like sites, respectively. Amorphous (AMO) sites are colored yellow. For visual clarity, monomers are shown with reduced dimensions.
Preprints 66761 g009
Figure 10. Degree of crystallinity, τ c , as a function of packing density, φ, and equilibrium bending angle, θ 0 , for 100-chains Nav = 12 systems of semi-flexible chains with bending constant k θ = 9 (filled symbols) and freely-jointed chains (open square). Dashed lines connecting points are used as visual support. Vertical dashed line corresponds to the melting transition ( φ m o n o m e r s M 0.545 ) for monomeric hard spheres [29].
Figure 10. Degree of crystallinity, τ c , as a function of packing density, φ, and equilibrium bending angle, θ 0 , for 100-chains Nav = 12 systems of semi-flexible chains with bending constant k θ = 9 (filled symbols) and freely-jointed chains (open square). Dashed lines connecting points are used as visual support. Vertical dashed line corresponds to the melting transition ( φ m o n o m e r s M 0.545 ) for monomeric hard spheres [29].
Preprints 66761 g010
Figure 11. Orientational order parameter, q (left y-axis, blue color), and crystallinity, τ c (right y-axis, red color), as a function of packing density, φ, for 100-chains Nav = 12 systems of semi-flexible chains with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0 ° (left) and θ 0 = 90 ° (right). Dashed lines connecting points are used as visual support. Vertical dashed lines correspond to the freezing ( φ m o n o m e r s F 0.494 ) and melting ( φ m o n o m e r s M = 0.545 ) transition for monomeric hard spheres [29].
Figure 11. Orientational order parameter, q (left y-axis, blue color), and crystallinity, τ c (right y-axis, red color), as a function of packing density, φ, for 100-chains Nav = 12 systems of semi-flexible chains with bending constant k θ = 9 and equilibrium bending angle θ 0 = 0 ° (left) and θ 0 = 90 ° (right). Dashed lines connecting points are used as visual support. Vertical dashed lines correspond to the freezing ( φ m o n o m e r s F 0.494 ) and melting ( φ m o n o m e r s M = 0.545 ) transition for monomeric hard spheres [29].
Preprints 66761 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated