Ray M. Sehgal

A Smoluchowski model of crystallization dynamics of small colloidal clusters

We investigate the dynamics of colloidal crystallization in a 32-particle system at a fixed value of interparticle depletion attraction that produces coexisting fluid and solid phases. Free energy landscapes (FELs) and diffusivity landscapes (DLs) are obtained as coefficients of 1D Smoluchowski equations using as order parameters either the radius of gyration or the average crystallinity. FELs and DLs are estimated by fitting the Smoluchowski equations to Brownian dynamics (BD) simulations using either linear fits to locally initiated trajectories or global fits to unbiased trajectories using Bayesian inference. The resulting FELs are compared to Monte Carlo Umbrella Sampling results. The accuracy of the FELs and DLs for modeling colloidal crystallization dynamics is evaluated by comparing mean first-passage times from BD simulations with analytical predictions using the FEL and DL models. While the 1D models accurately capture dynamics near the free energy minimum fluid and crystal configurations, predictions near the transition region are not quantitatively accurate. A preliminary investigation of ensemble averaged 2D order parameter trajectories suggests that 2D models are required to capture crystallization dynamics in the transition region.

Fokker-Planck analysis of separation dependent potentials and diffusion coefficients in simulated microscopy experiments

Total internal reflection microscopy (TIRM) and video microscopy (VM) are methods for nonintrusively measuring weak colloidal interactions important to many existing and emerging applications. Existing analyses of TIRM measured single particle trajectories can be used to extract particle-surface potentials and average particle diffusion coefficients. Here we develop a Fokker-Planck (FP) formalism to simultaneously extract both particle-surface interaction potentials and position dependent diffusion coefficients. The FP analysis offers several advantages including capabilities to measure separation dependent hydrodynamic interactions and nonequilibrium states that are not possible with existing analyses. The FP analysis is implemented to analyze Brownian dynamic simulations of single particle TIRM and VM experiments in several configurations. Relative effects of spatial and temporal sampling on the correct interpretation of both conservative and dissipative forces are explored and show a broad range of applicability for accessible experimental systems. Our results demonstrate the ability to extract both static and dynamic information from microscopy measurements of isolated particles near surfaces, which provides a foundation for further investigation of particle ensembles and nonequilibrium systems.

Colloidal cluster crystallization dynamics

The crystallization dynamics of a colloidal cluster is modeled using a low-dimensional Smoluchowski equation. Diffusion mapping shows that two order parameters are required to describe the dynamics. Using order parameters as metrics for condensation and crystallinity, free energy, and diffusivity landscapes are extracted from brownian dynamics simulations using bayesian inference. Free energy landscapes are validated against Monte Carlo simulations, and mean first-passage times are validated against dynamic simulations. The resulting model enables a low-dimensional description of colloidal crystallization dynamics.

Phase behavior of the 38-atom Lennard-Jones cluster.

We have developed a coarse-grained description of the phase behavior of the isolated 38-atom Lennard-Jones cluster (LJ38). The model captures both the solid-solid polymorphic transitions at low temperatures and the complex cluster breakup and melting transitions at higher temperatures. For this coarse model development, we employ the manifold learning technique of diffusion mapping. The outcome of the diffusion mapping analysis over a broad temperature range indicates that two order parameters are sufficient to describe the clusters phase behavior; we have chosen two such appropriate order parameters that are metrics of condensation and overall crystallinity. In this well-justified coarse-variable space, we calculate the clusters free energy landscape (FEL) as a function of temperature, employing Monte Carlo umbrella sampling. These FELs are used to quantify the phase behavior and onsets of phase transitions of the LJ38 cluster.

Onset of the crystalline phase in small assemblies of colloidal particles

We report the findings of a computational study designed to determine the onset of a stable crystalline phase in assemblies of small numbers (13–32) of colloidal particles that interact via a depletion-based short-ranged attractive potential. Using Monte Carlo umbrella sampling with coarse graining in two order parameters, we generate free-energy landscapes that can indicate coexistence between fluid-like and crystalline phases. The emergence of a stable crystalline phase is observed as the number of particles in the assembly increases beyond a critical value. We find that the critical cluster size for crystallization onset decreases with increasing strength of the interparticle attraction.