Abstract
We study the crystallisation processes occurring in a nonvibrating two-dimensional magnetic granular system at various fixed values of the effective temperature. In this system, the energy loss due to dissipative effects is compensated by the continuous energy input coming into the system from a sinusoidal magnetic field. When this balance leads to high values of the effective temperature, no aggregates are formed, because particles’ kinetic energy prevents them from aggregating. For lower effective temperatures, formation of small aggregates is observed. The smaller the values of the applied field’s amplitude, the larger the number of these disordered aggregates. One also observes that when clusters form at a given effective temperature, the average effective diffusion coefficient decreases as time increases. For medium values of the effective temperature, formation of small crystals is observed. We find that the sixth bond-orientational order parameter and the number of bonds, when considering more than two, are very sensitive for exhibiting the order in the system, even when crystals are still very small.
Similar content being viewed by others
Introduction
For a long time, crystallisation has attracted scientific attention, particularly through attempts to understand the nucleation mechanisms. Crystallisation under certain conditions can be seen as a phase transition from a disordered structure to a periodically ordered one. This process can be observed in some systems while cooling down, going from a fluid phase to a solid one. It is also observed in systems going from disordered solid states to ordered states when they are subjected to annealing. Crystallisation begins with the formation of a small core, a nucleus. The classical homogeneous nucleation theory (CNT)1,2 states that due to particle concentration fluctuations, small aggregates are formed. If these are smaller than a certain critical size, they disintegrate, but if they are bigger than the critical size, they become stable and grow. According to CNT, once an aggregate becomes stable, it grows by further aggregation, and its crystalline structure will determine the structure of the final aggregate. The gradual growth of the structure is due to the incorporation of new particles in minimal energy positions of the lattice. Crystallisation processes can also start with heterogeneous nucleation, whereby the nucleus is formed on impurities or on the borders of the system. Heterogeneous nucleation occurs more frequently than homogeneous nucleation, and it has been observed in confined systems3,4 and in cases where the walls that limit the system are modified to induce nucleation5.
It is well known that better understanding of crystallisation will facilitate the control of a wide variety of phenomena. At the atomic scale, crystallisation processes are studied by indirect methods due to their microscopic nature. Furthering our in-depth understanding of crystallisation processes requires direct studies capable of describing them at the level of single particles. These kinds of studies would allow us to validate nucleation theories; in particular, they could provide evidence to support or discard non-classical nucleation theories. No studies provide sufficient spatial and temporal resolution to describe crystallisation processes at the molecular level, although some notable progress in that direction has been achieved6,7,8,9. For many years, colloids have served as models to describe the crystallisation process in two10 and three dimensions11,12,13,14. For instance, it has been found that under certain conditions, non-classical two-step nucleation occurs15,16. This result is contrary to what is proposed by CNT, where the crystallisation process occurs in a single step. A two-step mechanism non-classical theory of nucleation establishes that a cluster of particles of sufficiently large size may arise with a disordered structure, which subsequently rearranges and gives rise to a crystalline nucleus17,18.
Even though colloidal systems provide some insight into the true nature of crystallisation, they present some limitations. For instance, particle concentration and effective temperature are coupled: the effective temperature is proportional to the inverse of the particle concentration. It is desirable to control particle concentration and effective temperature independently. To overcome some of the issues arising when modelling crystallisation with colloidal systems, and to improve our knowledge, we have been using granular matter systems, in which particles are clearly athermal, because they are not affected by thermal fluctuations19. To introduce something analogous to an effective temperature in these systems, it is necessary to inject energy through some mechanisms.
Most experiments involving granular systems use mechanical vibration to fluidise the set of particles, thus controlling the effective temperature. In some other cases, periodic shearing has been used to induce structural changes in the system20,21,22. Despite being dissipative systems, granular models present structural configurations equivalent to systems in equilibrium23,24. Colloidal systems and granular systems allow us to follow the individual trajectories of the constituent particles and therefore to obtain the structural and dynamic properties. In Ref.21, a combination of vibration and shearing was used to produce reversible phase transitions between crystalline phases. In Refs.23,25, a vibrating granular system was used to measure the structural changes during the crystallisation and glass transition processes in a two-dimensional system.
Previously, we have used a nonvibrating granular system under a sinusoidal magnetic field as a model to investigate, at a particle level, some pattern formation processes in fluids. According to the theory of Ornstein–Uhlenbeck processes, the resulting motion has all the Brownian motion characteristics, as we have shown in Refs.24,26. We have studied the behaviour of the system when it suffers a sudden quenching. In those experiments, particles were under a sinusoidal magnetic field with an offset which produces effective repulsive interaction between particles. This effective repulsion prevents particles from aggregating, and then the system is able to model a glass-transition-like transformation. We have also studied crystallisation by using a lens27, and a tilted planar cell28. In those experiments, we have produced a cooling down process by decreasing the magnetic field using a step-down-like signal.
A special characteristic of the granular nonvibrating system is that particle concentration and effective temperature can be independently controlled for low and medium particle concentrations, which is a very attractive advantage over other macroscopic systems. We can carry out experiments at a fixed particle concentration while varying the effective temperature. It has been shown that at low particle concentrations, the effective temperature is directly proportional to the amplitude of the applied magnetic field24. By adding a constant magnetic field, the intensity of the interactions between particles can be controlled. When the system is under a sinusoidal magnetic field with an offset, the effective repulsion prevents particles from aggregating. Without the offset and below a magnetic threshold, particles tend to aggregate. As the particle concentration increases, the size of the aggregates grows. Depending on the effective temperature, the aggregates are more or less ordered. If the temperature is relatively high, these aggregates disintegrate, but below a certain threshold, the aggregates are stable and tend to grow.
In the present work, we analyse aggregates formed under a long-term regime at a constant temperature; we considered several effective temperature values. Because the aggregates are formed by a few particles, we determine the structural characteristics by using the sixth bond-orientational order parameter and the number of bonds. To standardise, we start the experiments at a high temperature and then abruptly set the desired temperature in the system. The experiments are recorded during a relatively long period to observe whether aggregates are formed or not, and to determine the structural characteristics of the clusters if aggregation occurs.
Results
The main advantage of using this granular system to study crystallisation processes is that simple optical microscopy allows observation, in real time, of complex phenomena in which energetic and entropic interactions generate the structure patterns. For instance, one may observe at fixed particle concentrations the crystallisation process for several fixed effective temperature values. Under these conditions, we observe that between high temperatures, at which the kinetic energy of particles prevents stable aggregates forming, and low temperatures, at which numerous disordered aggregates occur, medium temperatures promote the formation of aggregates with a more ordered structure. It is observed that most of the aggregates are small crystals.
We carried out three series of experiments (series 1, series 2 and series 3) for each value of the effective temperature. The configurations of the granular system after 30 mins at three different temperatures (from series 3) are shown in Fig. 1. The circular walls of the cell prevent particles from escaping. It is observed that a ring of particles is formed at the wall container. Particles near the circular wall suffer more collisions from particles inside the sample than those near the wall. Therefore, the entropic interactions are the origin of the formation of the ring29,30,31,32. No aggregates are formed at a relatively high effective temperature, but below an effective threshold temperature, they start to form. The number of aggregates increases as the effective temperature decreases. It is observed that the lower the effective temperature, the bigger the disorder of the inner cluster structure. Some aggregates are anchored to the cell border, when they are crystalline this is due to heterogeneous nucleation33,34,35,36. This effect increases as the temperatures are increased. For lower temperatures, the formation of aggregates separated from the wall is observed, when the aggregates are ordered this is preceded by a two-dimensional homogeneous nucleation37,38,39.
We characterised the dynamical nature of the system by using mean square displacement, \(<\Delta r^{2}\hbox {(t)}>=<|r(t)-r(0)|^2>\), \(\Delta r\) is measured in particle diameter units \(\sigma \), for several values of the applied magnetic field. The results are shown in Fig. 2. The first six points correspond to the quasi-ballistic regime40, and these points are not used to calculate the diffusion coefficient, which was determined by using the following 30 points. In the range of the effective temperature here explored, the \(<\Delta r^{2}\hbox {(t)}>\) curves show that particles’ behaviour is diffusive, except for the first six points. Moreover, as the effective temperature increases, the slope of the curve also increases. This means that when the effective temperature increases, the particles move at a higher speed. To show this, we determined the diffusion coefficient. The effective diffusion coefficient (D) was obtained by means of a linear fit of the \(<\Delta r^{2}\hbox {(t)}>\) (ignoring the first points) and from the relation \(<\Delta r^{2}\hbox {(t)}> \sim 4D t\). The time interval in which the effective diffusion coefficient is calculated is shown in Fig. 2, delimited by vertical dotted lines. Figure 3 shows the time evolution of the diffusion coefficient over the range of magnetic fields (effective temperature) analysed for one of the sets of experiments (series 1). At high temperatures, the diffusion coefficient remains almost constant, except for some fluctuations. At lower temperatures, it decreases as a function of time. For low and medium temperatures, we observe that after 200 s, D decreases slightly and then remains constant until the final measurement time.
Figure 4 shows the behaviour of the average diffusion coefficient as a function of the amplitude of the applied field—namely, our effective temperature. At least three regions in the behaviour of the average D are observed, clearly distinguished by their slope at high, medium and low temperatures. At low temperatures (black symbols), the particles move more slowly and form aggregates all over the surface, mostly disordered. In the medium region (green symbols), the particles move at a higher speed, with the formed aggregates being more organised and tending to form compact hexagonal arrangements, which are mostly located in the central region of the cell. Finally, in the high-temperature region (red symbols), the particles move quickly, and their kinetic energy almost prevents cluster formation, which is observed only in the formation of very few small ordered aggregates. At even higher temperatures, the only observed aggregates are anchored to the wall of the cell, forming part of a ring-like structure. The inset in Fig. 4 shows the raw data, and its arithmetic averaged for the three values of the effective temperature.
We have found that for granular systems such as ours, to evaluate the degree of order of the aggregates as a function of the number of neighbours (n) around a given particle, the sixth bond-orientational order parameter \(\Psi _{6}\) works well27,41. Details are discussed in the “Methods” section.
To discuss the order and compactness of the aggregates’ structure, we choose an effective temperature in the medium range, where aggregates are more ordered, and most of them are crystalline. In Fig. 1b, it can be seen several ordered aggregates, both at the centre and at the cell’s border. In Fig. 5a, particles are coloured according to the number of neighbours; the dark blue colour represents those particles with six neighbours. Comparing the picture and Fig. 1b, we observed that many particles have six neighbours even though some of them are far away from the particle. This is a characteristic of determining neighbours using the Delaunay triangulation or its dual graph, the Voronoi diagram42. Based on these neighbours, the sixth bond-orientational order parameter \(\Psi _6\) is evaluated. Figure 5b shows particles coloured according to this quantity. It is observed that some particles are indicated as disordered even though they belong to an ordered cluster. This is a direct consequence of taking into account all the neighbours obtained through the Delaunay triangulation. We observed that these parameters are not sensitive to small ordered aggregates. For characterisation of those small ordered aggregates that form a compact hexagonal network, the bounded neighbours (\(n_B\)) are calculated—i.e., those particles in contact with a particle within a distance of one particle diameter \(r=\sigma \). Based on these bounded neighbours, a new parameter is introduced (\(\Psi ^{'}_{6}\)); this is similar to \(\Psi _{6}\) but measured considering only bounded neighbours. The perfect sixth bond-orientational order configuration is reached when angles between the nearest neighbours are multiples of \(\pi /3\). Considering the above configuration, we plot the number of the bounded neighbours, \(n_B\), and the parameter, \(\Psi ^{'}_{6}\), in Fig. 5c and d, respectively. Those particles having six neighbours in contact are coloured dark blue, while particles without neighbours in contact are red coloured. Also, we found that the parameter \(\Psi ^{'}_{6}\) adequately recognises the ordered aggregates. These clusters appear enclosed in orange circles in Fig. 1b.
When the number of bounds is just one, a particle could be unstable and is likely to abandon the aggregate. That is why we further modified the sixth bond-orientational order parameter by considering only those particles whose number of bounded neighbours is equal to or greater than two, represented by the parameter \(n^{''}_{B}\). The arithmetic average of the three values of \(n^{''}_{B}\) for the three series, at different magnetic field values, is shown in Fig. 6a. We compare the initial configuration (measured at short time intervals) to the final configuration (measured at long time intervals). This figure shows behaviour similar to the effective diffusion coefficient as a function of effective temperature. At low temperatures, the number of particles with two or more neighbours in contact is small; as the temperature increases, the number of the bounded neighbours in contact increases, which allows the formation of compact hexagonal networks that subsequently lead to bigger ordered aggregates. Finally, at high temperatures, the particles move with a higher velocity and form a few aggregates; however, these are made up of hexagonal networks. The corresponding sixth bond-orientational order parameter, considering two or more bounded neighbours, is denoted \(\Psi ^{''}_{6}\). This parameter allows the identification of the most stable and well-ordered formations. Figure 6b shows the average of \(\Psi ^{''}_{6}\) as a function of the magnetic field. From this figure, it is more evident that there are three regimes in the behaviour of the structure of the aggregates. Thus, it supports the observation that at medium temperatures, the aggregates are well ordered, forming hexagonal lattices.
As mentioned, the aggregates are stable and only increase in size as a function of time. To demonstrate this fact, an intermediate temperature is considered, and the time evolution of the average number of aggregates, \(N_a\), and the average size of aggregates, \(<N>\), were analysed. Figure 7a shows that the number of aggregates remains constant over the measurement interval, which means that the aggregates that overpass the critical size will never disaggregate. On the other hand, Figure 7b shows the temporal behaviour of \(<N>\). The largest observed clusters are about N = 30, where N is the particle number. We notice that the aggregates at the early stage grow quickly and then hold their average size, before finally growing again. Experimentally, we observe that the disordered aggregates diminish in size while the ordered ones increase, in such a way that the average size remains constant as it leads the plateau of the curve. The above is confirmed when we only analyse the time evolution of the ordered aggregates (\(\Psi ^{''}_{6}\) > 0.8). Figure 8a shows the behaviour of \(<N>\) with a high degree of order as a function of time. Two characteristic growth ratios are found. In the first one (0–50 s), the aggregates form rapidly and reach the critical size; in the second one, they exhibit a lower rate of growth. Figure 8b shows the time evolution of \(\Psi ^{''}_{6}\). The aggregates experience a fast ordering until 150 s, and then they rearrange, reaching a higher degree of order. After that (335 s), the aggregates only grow epitaxially.
This idea that ordered aggregates only grow in size by the adhesion of free particles is consistent with the effective diffusion coefficient decreasing when the aggregates’ size increases. To demonstrate this statement, Fig. 9a and b show the relationship between the effective diffusion coefficient D and the bond-orientational order parameters that characterise the degree of order of the aggregates, \(n^{''}_{B}\) and \(\Psi ^{''}_{6}\), respectively. In both cases, the diffusion coefficient decreases as the number of bounded neighbours and the degree of order of the aggregates increase. This is attributable to the fact that as the size of the aggregates increases, the number of free particles decreases. To represent the decrease, an exponential fit was done in both cases; in Fig. 9a and b, the solid line corresponds to this fit.
In summary, we have studied the particle kinetics and aggregation processes in a two-dimensional granular system at constant effective temperature. The essential founds of this experimental analysis are the following: first, we have observed that the diffusion coefficient decreases in time due to the formation and growth of clusters at a constant effective temperature. A second one is that the sixth bond-orientational order parameter and the number of bounded neighbours, when considering more than two, are appropriate descriptors for evolution of the degree order of the aggregates. Third, we have observed three regimes of temperature: low, medium, and high. At high temperature the formation of aggregates is hardly observed. At medium and low temperatures, we found the formation of several aggregates. Particularly, at medium temperature the order of the aggregates is higher, the aggregates are crystalline and numerous. In Supplementary Movies 1–3, one observes the temporal evolution of the aggregates at low ( first 30 s), medium (at 300 s), and high temperature (at 300 s), respectively.
The classical theory of nucleation assumes that once the nucleus is formed, the ulterior aggregation process occurs practically as an epitaxial growth, propagating the nucleus crystalline structure. From this analysis, we may conclude that, after the energetic balance that allows the formation of small disordered aggregates, these structures evolve, becoming more compact and ordered. Thus, our observations indicate that there is a threshold condition in which the volumetric potential becomes predominant over the surface potential, propitiating from that point and beyond the formation of disordered aggregates. Then, a second threshold is needed for the flyover to achieve order of these aggregates. In this sense, the growth of the aggregates observed in our experiments is consistent with a two-step nucleation theory43,44.
Methods
To characterise the different structural characteristics’ phases, the radial distribution function (RDF) is commonly used. This quantity is useful when the measured characteristic phase is sufficiently large to capture a robust average by the RDF. The domain of the RDF must be spanned by at least some tens of the aggregate size. When the sample contains small parts of a different phase, their characteristics are not clearly shown in the RDF. To deal with small aggregates, some other quantities have been developed. That is the case with the sixth bond-orientational order parameter, which is used when aggregates tend to be ordered in the hexagonal close-packed arrangement. Several authors have used this parameter to quantify the two-dimensional crystallisation process45,46,47. It also allows us to identify phase transitions—for instance, from a liquid to a crystalline lattice48,49.
To prepare our system, we used 1050 spherical particles that settled on the horizontal glass surface bounded by a circular wall of 70 mm in diameter. The particle concentration is the area fraction occupied by the particles \(\phi _{2D}\)=0.19. The particles are steel balls of 1 mm in diameter, ANSI 420 grade 1000 by Gimex S.A. The glass plane is allocated in the centre of a pair of Helmholtz coils, which produce a vertical magnetic field. The coils are fed by a Kepco BOP 36-6 M power amplifier. This system of magnetic spheres is subjected to a sinusoidal magnetic field \(B=B_{0}\sin {(2\pi ft)}\), whose amplitude (\(B_{0}\)) was kept constant during each experiment and takes values from 22 to 50 G. We carried out three sets of experiments: series 1, series 2 and series 3. The frequency (f) was fixed at 9.24 Hz. The applied magnetic field plays an effective temperature role. Henceforth, we refer to the amplitude of the magnetic field as the effective temperature24,40. Figure 10 shows the experimental setup.
The dynamics of the system are explained as follows. Each particle has a magnetic moment. As the magnetic field changes its direction (up and down), a particle’s magnetic moment tends to change its direction following the magnetic field’s direction, and minimising its potential energy. When the particle’s magnetic moment points opposite to the field direction, it is in a very unstable state. Thus, it rotates in a random direction to align with the magnetic field. In following the magnetic field direction, the particle rotates and rolls, experiencing repulsive contact forces due to the other particles that tend to separate them.
A HandyCam digital camera was used to record the dynamics of the system. We obtained 5 mins of the video at a standard resolution of 720,540 pixels at 30 fps in AVI interlaced format. At this capture rate, the particle positions were not always well defined; thus, we used a filter to deinterlace the video frames to obtain well-defined image sequences with a resolution of \(\tau =1/60\)s. Subsequently, the image sequences underwent a cleaning process to calculate the particles’ trajectories using the program ImageJ and its plugin Mosaic50,51. The mean square displacement was calculated from these trajectories. Figure 11 shows the superimposed trajectories of all the particles during a time interval for the case when the magnetic field’s amplitude is 44 G. Those particles on the circular boundary do not collide with other particles, only they slightly move about their positions. Similarly, particles that form stable and ordered aggregates, such as those enclosed in white circles in Fig. 11, do not move either. Only free particles can move over the entire surface. In the case of Fig. 11, the effective temperature is high, so the particles follow long trajectories. The diffusion coefficient, D, was obtained by the Einstein fluctuation-dissipation relation in a diffusive two-dimensional system,\(<\Delta r^{2}\hbox {(t)}>=4Dt\).
To determine structural parameters, particles are first detected and separated from each other by using ImageJ and some of its plugins—see Fig. 12a. The plugin Delaunay-Voronoi in ImageJ was used to generate the Delaunay triangulation, as shown in Fig. 12b. Lastly, the sixth bond-orientational order parameter is calculated by using the expression:
where \(\theta _{ij}\) is the angle between the line formed by the reference particle i and neighbour j and the x-axis, and \(n_i\) is the number of neighbours determined on the basis of the Delaunay triangulation27. The value of the average of the sixth bond-orientational order parameter, \(<\Psi _{6}>\), tells us whether an aggregate is stable and ordered or if it is in a disordered phase—i.e., if \(<\Psi _{6}>\)= 1, the aggregates form perfect hexagonal networks, and if \(<\Psi _{6}>\)= 0, it is a strongly disordered system.
Figure 12c shows the number of neighbours, n, for a particle array, and Fig. 12d shows the sixth bond-orientational order parameter \(\Psi _{6}\). These parameters did not fully represent the hexagonal compact arrangement. Therefore, the bounded neighbours \(n_B\) calculation and the sixfold bond-orientational order parameter from the bounded neighbours \(\Psi ^{'}_{6}\) were introduced. These quantities are represented in Fig. 12e and f, respectively. The bounded neighbours were determined considering the nearest neighbours within a distance of one particle diameter. The sixth bond-orientational order parameter \(\Psi ^{'}_{6}\) is calculated considering only the bounded neighbours.
References
Oxtoby, D. W. Nucleation of first-order phase transitions. Acc. Chem. Res. 31, 91. https://doi.org/10.1021/ar9702278 (1998).
Kelton, K. F. Crystal nucleation in liquids and glasses. Solid State Phys. 45, 75. https://doi.org/10.1016/S0081-1947(08)60144-7 (1991).
Sandomirski, K., Allahyarov, E., Löwen, H. & Egelhaaf, S. Heterogeneous crystallization of hard-sphere colloids near a wall. Soft Matter 7, 8050. https://doi.org/10.1039/C1SM05346A (2011).
Block, B. et al. Computer simulation of heterogeneous nucleation of colloidal crystals at planar walls. Eur. Phys. J. Spec. Top. 223, 347. https://doi.org/10.1140/epjst/e2014-02095-0 (2014).
Assoud, L., Messina, R. & Löwen, H. Heterogeneous crystallization near structured walls in quenched two-dimensional binary colloidal suspensions. Mol. Phys. 109, 7. https://doi.org/10.1080/00268976.2011.562871 (2011).
Sauter, A. et al. On the question of two-step nucleation in protein crystallization. Faraday Discuss. 179, 41. https://doi.org/10.1039/C4FD00225C (2015).
Schreiber, R. E. et al. Real-time molecular scale observation of crystal formation. Nat. Chem. 9, 369. https://doi.org/10.1038/nchem.2675 (2017).
Chen, P. Z., Niu, L., Zhang, H., Chen, Y. Z. & Yang, Q. Z. Exploration of the two-step crystallization of organic micro/nano crystalline materials by fluorescence spectroscopy. Mater. Chem. Front. 2, 1323. https://doi.org/10.1039/C8QM00118A (2018).
Nakamuro, T., Skakibara, M., Nada, H., Harano, K. & Nakamura, E. Capturing the moment of emergence of crystal nucleus from disorder. J. Am. Chem. Soc. 143, 1763. https://doi.org/10.1021/jacs.0c12100 (2021).
Gasser, U. Crystallization in three- and two-dimensional colloidal suspensions. J. Phys. Condens. Matter 21, 203101. https://doi.org/10.1088/0953-8984/21/20/203101 (2009).
van Megen, W. & Pusey, P. N. Dynamic light-scattering study of the glass transition in a colloidal suspension. Phys. Rev. A 43, 5429. https://doi.org/10.1103/physreva.43.5429 (1991).
Palberg, T. Crystallization kinetics of colloidal model suspensions: Recent achievements and new perspectives. J. Phys. Condens. Matter 26, 333101. https://doi.org/10.1088/0953-8984/26/33/333101 (2014).
Palberg, T., Stipp, A. & Bartsch, E. Unusual crystallization kinetics in a hard sphere colloid-polymer mixture. Phys. Rev. Lett. 102, 038302. https://doi.org/10.1103/PhysRevLett.102.038302 (2009).
Ojeda-Mendoza, G. J., Moncho-Jordá, A., González-Mozuelos, P., Haro-Pérez, C. & Rojas-Ochoa, L. F. Evidence of electrostatic-enhanced depletion attraction in the structural properties and phase behavior of binary charged colloidal suspensions. Soft Matter 14, 1355. https://doi.org/10.1039/C7SM02220D (2018).
Zhang, T. & Liu, X. T. How does a transient amorphous precursor template crystallization. J. Am. Chem. Soc. 129, 13520. https://doi.org/10.1021/ja073598k (2007).
Savage, J. R. & Dinsmore, A. D. Experimental evidence for two-step nucleation in colloidal crystallization. Phys. Rev. Lett. 102, 198302. https://doi.org/10.1103/PhysRevLett.102.198302 (2009).
Russo, J. & Tanaka, H. Nonclassical pathways of crystallization in colloidal systems. MRS Bull. 41, 369. https://doi.org/10.1557/mrs.2016.84 (2016).
Cölfen, H. Nonclassical nucleation and crystallization. Crystals 10, 61. https://doi.org/10.3390/cryst10020061 (2020).
Jaeger, H. M., Nagel, S. R. & Behringer, R. P. Granular solids, liquids, and gases. Rev. Mod. Phys. 68, 1259. https://doi.org/10.1103/RevModPhys.68.1259 (1996).
Huerta, D. A., Sosa, V., Vargas, M. C. & Ruiz-Suárez, J. C. Archimedes’ principle in fluidized granular systems. Phys. Rev. E 72, 031307. https://doi.org/10.1103/PhysRevE.72.031307 (2005).
Carvente, O. & Ruiz-Suárez, J. C. Crystallization of confined non-Brownian spheres by vibrational annealing. Phys. Rev. Lett. 95, 018001. https://doi.org/10.1103/PhysRevLett.95.018001 (2005).
Panaitescu, A., Reddy, K. A. & Kudrolli, A. Nucleation and crystal growth in sheared granular sphere packings. Phys. Rev. Lett. 108, 108001. https://doi.org/10.1103/PhysRevLett.108.108001 (2012).
Reis, P. M., Ingale, R. A. & Shattuck, M. D. Crystallization of a quasi-two-dimensional granular fluid. Phys. Rev. Lett. 96, 258001. https://doi.org/10.1103/PhysRevLett.96.258001 (2006).
Tapia-Ignacio, C., Garcia-Serrano, J. & Donado, F. Nonvibrating granular model for a glass-forming liquid: Equilibration and aging. Phys. Rev. E 94, 062902. https://doi.org/10.1103/PhysRevE.94.062902 (2016).
Reis, P. M., Ingale, R. A. & Shattuck, M. D. Caging dynamics in a granular fluid. Phys. Rev. Lett. 98, 188301. https://doi.org/10.1103/PhysRevLett.98.188301 (2007).
Donado, F., Moctezuma, R. E., López-Flores, L. L., Medina-Noyola, M. & Arauz-Lara, J. L. Brownian motion in non-equilibrium systems and the Ornstein–Uhlenbeck stochastic process. Sci. Rep. 7, 12614. https://doi.org/10.1038/s41598-017-12737-1 (2017).
Escobar, A., Tapia-Ignacio, C., Donado, F., Arauz-Lara, J. L. & Moctezuma, R. E. Glass-and crystal-forming model based on a granular two-dimensional system. Phys. Rev. E 101, 052907. https://doi.org/10.1103/PhysRevE.101.052907 (2020).
Sánchez-Miranda, M., Carrillo-Estrada, J. L. & Donado, F. Crystallization processes in a nonvibrating magnetic granular system with short range repulsive interaction. Sci. Rep. 9, 3531. https://doi.org/10.1038/s41598-019-40062-2 (2019).
Castañeda Priego, R., Rodríguez-López, A. & Méndez-Alcaraz, J. M. Entropic forces in dilute colloidal systems. Phys. Rev. E 73, 051404. https://doi.org/10.1103/PhysRevE.73.051404 (2006).
Yodh, A. G. et al. Entropically driven self-assembly and interaction in suspension. Philos. Trans. R. Soc. A. 359, 921. https://doi.org/10.1098/rsta.2000.0810 (2001).
Ziese, F., Maret, G. & Gasser, U. Heterogeneous nucleation and crystal growth on curved surfaces observed by real-space imaging. J. Phys. Condens. Matter 25, 375105. https://doi.org/10.1088/0953-8984/25/37/375105 (2013).
Zöttl, A. & Stark, H. Emergent behavior in active colloids. J. Phys. Condens. Matter 28, 253001. https://doi.org/10.1088/0953-8984/28/25/253001 (2016).
Liu, X. Y., Maiwa, K. & Tsukamoto, K. Heterogeneous two-dimensional nucleation and growth kinetics. J. Chem. Phys. 106, 1870. https://doi.org/10.1063/1.473325 (1997).
Abedi, S., Chen, C. & Vanapalli, S. Collective nucleation dynamics in two-dimensional emulsions with hexagonal packing. Phys. Rev. E 101, 030602. https://doi.org/10.1103/PhysRevE.101.030602 (2020).
Espinosa, J., Vega, C., Valeriani, C., Frenkeld, D. & Sanz, E. Heterogeneous versus homogeneous crystal nucleation of hard spheres. Soft Matter 15, 9625. https://doi.org/10.1039/C9SM01142K (2019).
Sandomirski, K. et al. Heterogeneous crystallization of hard and soft spheres near flat and curved walls. Eur. Phys. J. Spec. Top. 223, 439. https://doi.org/10.1140/epjst/e2014-02101-7 (2014).
Gasser, U., Weeks, E., Schofield, A., Pusey, P. & Weitz, D. A. Real-space imaging of nucleation and growth in colloidal crystallization. Science 292, 258. https://doi.org/10.1126/science.1058457 (2001).
Rietz, F., Radin, C., Swinney, H. L. & Schröter, M. Nucleation in sheared granular matter. Phys. Rev. Lett. 120, 055701. https://doi.org/10.1103/PhysRevLett.120.055701 (2018).
Radin, C. & Swinney, H. Phases of granular matter. J. Stat. Phys. 175, 542. https://doi.org/10.1007/s10955-018-2144-4 (2019).
Tapia-Ignacio, C., Moctezuma, R. E., Donado, F. & Weeks, E. R. Brownian motion of ellipsoidal particles on a granular magnetic bath. Phys. Rev. E 102, 022902. https://doi.org/10.1103/PhysRevE.102.022902 (2020).
Zangi, R. & Rice, S. A. Phase transitions in a quasi-two-dimensional system. Phys. Rev. E 58, 7529. https://doi.org/10.1103/PhysRevE.58.7529 (1998).
Moctezuma, R. E., Arauz-Lara, J. L. & Donado, F. A. Structural characterization of a magnetic granular system under a time-dependent magnetic field: Voronoi tessellation and multifractal analysis. Physica A 496, 27. https://doi.org/10.1016/j.physa.2017.12.123 (2018).
Zhang, T. & Liu, X. T. Nucleation: What happens at the initial stage?. Angew. Chem. Int. Ed 48, 1308. https://doi.org/10.1002/anie.200804743 (2009).
Dillmann, P., Maret, G. & Keim, P. Two-dimensional colloidal systems in time-dependent magnetic fields. Eur. Phys. J. Spec. Top. 222, 2941. https://doi.org/10.1140/epjst/e2013-02068-9 (2013).
Wang, Z., Alsayed, A., Yodh, G. & Han, Y. Two-dimensional freezing criteria for crystallizing colloidal monolayers. J. Chem. Phys. 132, 154501. https://doi.org/10.1063/1.3372618 (2010).
Horn, T., Deutschländer, S., Löwen, H., Maret, G. & Keim, P. Fluctuations of orientational order and clustering in a two-dimensional colloidal system under quenched disorder. Phys. Rev. E 88, 062305. https://doi.org/10.1103/PhysRevE.88.062305 (2013).
Li, Y.-W. & Ciamarra, M. Attraction tames two-dimensional melting: From continuous to discontinuous transitions. Phys. Rev. Lett. 124, 218002. https://doi.org/10.1103/PhysRevLett.124.218002 (2020).
Li, B., Xiao, X., Wang, S., Wen, W. & Wang, Z. Real-space mapping of the two-dimensional phase diagrams in attractive colloidal systems. Phys. Rev. X 9, 031032. https://doi.org/10.1103/PhysRevX.9.031032 (2019).
Digregorio, P. et al. Full phase diagram of active Brownian disks: From melting to motility-induced phase separation. Phys. Rev. Lett. 121, 098003. https://doi.org/10.1103/PhysRevLett.121.098003 (2018).
Sbalzarini, I. & Koumoutsakos, P. Feature point tracking and trajectory analysis for video imaging in cell biology. J. Struct. Biol. 151, 182. https://doi.org/10.1016/j.jsb.2005.06.002 (2005).
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. Nih image to imagej: 25 years of image analysis. Nat. Method 9, 671. https://doi.org/10.1038/nmeth.2089 (2012).
Acknowledgements
Partial financial support by CONACyT, México, through grants A1-S-39909, 80629, 256176 (Ciencia Básica) and 440 (Fronteras de la Ciencia) is acknowledged. M.L-M. acknowledges the CONACyT postdoctoral fellowship through the grant A1-S-39909.
Author information
Authors and Affiliations
Contributions
F.D. conceived and developed the experimental system, F.D. and M.L.-M. performed the experiments, F.D., M.L.-M. and J.L.C. analysed and discussed the results. The three authors participated in the preparation of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Supplementary Movie 1.
Supplementary Movie 2.
Supplementary Movie 3.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ledesma-Motolinía, M., Carrillo-Estrada, J.L. & Donado, F. Crystallisation in a two-dimensional granular system at constant temperature. Sci Rep 11, 16531 (2021). https://doi.org/10.1038/s41598-021-96099-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-96099-9
This article is cited by
Two-step crystallisation in a 2D active magnetic granular system confined by a parabolic potential
Scientific Reports (2023)
Network efficiency of spatial systems with fractal morphology: a geometric graphs approach
Scientific Reports (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.