Source author record

Subir K. Das

Subir K. Das appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2020arXiv

Active Particles in Explicit Solvent: Dynamics of clustering for alignment interaction

We study dynamics of clustering in systems containing active particles that are immersed in an explicit solvent. For this purpose we have adopted a hybrid simulation method, consisting of molecular dynamics and multi-particle collision dynamics. In our model, overlap-avoiding passive interaction of an active particle with another active particle or a solvent particle has been taken care of via variants of Lennard-Jones potential. Dynamic interaction among the active particles has been incorporated via the Vicsek-like self-propulsion that facilitates clustering. We quantify the effects of activity and importance of hydrodynamics on the dynamics of clustering via variations of relevant system parameters. We work with low overall density of active particles. For this the morphology consists of disconnected clusters, the mechanism of growth switching among particle diffusion, diffusive coalescence and ballistic aggregation, depending upon the presence or absence of active and hydrodynamic interactions. Corresponding growth laws have been quantified and discussed in the background of appropriate theoretical pictures. Our results suggest that multi-particle collision dynamics is an effective method for investigation of hydrodynamic phenomena even in active matter systems.

preprint2020arXiv

Aging Exponents for Nonequilibrium Dynamics following Quenches from Critical Point

Via Monte Carlo simulations we study nonequilibrium dynamics in the nearest-neighbor Ising model, following quenches to points inside the ordered region of the phase diagram. With the broad objective of quantifying the nonequilibrium universality classes corresponding to spatially correlated and uncorrelated initial configurations, in this paper we present results for the decay of the order-parameter autocorrelation function for quenches from the critical point. This autocorrelation is an important probe for the aging dynamics in far-from-equilibrium systems and typically exhibits power-law scaling. From the state-of-the-art analysis of the simulation results we quantify the corresponding exponents ($\mathbfλ$) for both conserved and nonconserved (order parameter) dynamics of the model, in space dimension $d=3$. Via structural analysis we demonstrate that the exponents satisfy a bound. We also revisit the $d=2$ case to obtain more accurate results. It appears that irrespective of the dimension, $λ$ is same for both conserved and nonconserved dynamics.

preprint2020arXiv

Dependence of Cluster Growth on Coefficient of Restitution in a Cooling Granular Fluid

Starting from configurations having homogeneous spatial density, we study kinetics in a two-dimensional system of inelastically colliding hard particles, a popular model for cooling granular matter. Following an initial time period, the system exhibits a crossover to an inhomogeneous regime that is characterized by the formation and growth of particle-rich clusters. We present results on the time dependence of average mass of the clusters and that of average kinetic energy, obtained via event driven molecular dynamics simulations, for a wide range of values for the coefficient of restitution ($e$), by fixing the overall density of particles in the system to a constant number. The time of onset of crossover from homogeneous to the inhomogeneous regime, as is well known, strongly increases as one moves towards the elastic limit. Nevertheless, our presented results suggest that the asymptotic growth is independent of $e$, for uniform definition of cluster, onset of which has a different $e$-dependence than the onset of above mentioned crossover. In other words, not only the exponent but also the amplitude of the power-law growth, which is widely believed to be the form of the evolution, is at the most very weakly sensitive to the choice of $e$. While it is tempting to attribute this fact to the similar feature in the decay of energy, we caution that our current understanding is not matured enough to draw such a connection between cluster growth and energy decay in a meaningful manner.

preprint2020arXiv

Kinetics of Domain Growth and Aging in a Two-Dimensional Off-lattice System

We have used molecular dynamics simulations for a comprehensive study of phase separation in a two-dimensional single component off-lattice model where particles interact through the Lennard-Jones potential. Via state-of-the-art methods we have analyzed simulation data on structure, growth and aging for nonequilibrium evolutions in the model. These data were obtained following quenches of well-equilibrated homogeneous configurations, with density close to the critical value, to various temperatures inside the miscibility gap, having vapor-"liquid" as well as vapor-"solid" coexistence. For the vapor-liquid phase separation we observe that $\ell$, the average domain length, grows with time ($t$) as $t^{1/2}$, a behavior that has connection with hydrodynamics. At low enough temperature, a sharp crossover of this time dependence to a much slower, temperature dependent, growth is identified within the time scale of our simulations, implying "solid"-like final state of the high density phase. This crossover is, interestingly, accompanied by strong differences in domain morphology and other structural aspects between the two situations. For aging, we have presented results for the order-parameter autocorrelation function. This quantity exhibits data-collapse with respect to $\ell/\ell_w$, $\ell$ and $\ell_w$ being the average domain lengths at times $t$ and $t_w$ ($\leq t$), respectively, the latter being the age of a system. Corresponding scaling function follows a power-law decay: $~\sim (\ell/\ell_w)^{-λ}$, for $t\gg t_w$. The decay exponent $λ$, for the vapor-liquid case, is accurately estimated via the application of an advanced finite-size scaling method. The obtained value is observed to satisfy a bound.

preprint2020arXiv

Relaxation in a Phase-separating Two-dimensional Active Matter System with Alignment Interaction

Via computer simulations we study kinetics of pattern formation in a two-dimensional active matter system. Self-propulsion in our model is incorporated via the Vicsek-like activity, i.e., particles have the tendency of aligning their velocities with the average directions of motion of their neighbors. In addition to this dynamic or active interaction, there exists passive inter-particle interaction in the model for which we have chosen the standard Lennard-Jones form. Following quenches of homogeneous configurations to a point deep inside the region of coexistence between high and low density phases, as the systems exhibit formation and evolution of particle-rich clusters, we investigate properties related to the morphology, growth and aging. A focus of our study is on the understanding of the effects of structure on growth and aging. To quantify the latter we use the two-time order-parameter autocorrelation function. This correlation, as well as the growth, is observed to follow power-law time dependence, qualitatively similar to the scaling behavior reported for passive systems. The values of the exponents have been estimated and discussed by comparing with the previously obtained numbers for other dimensions as well as with the new results for the passive limit of the considered model. We have also presented results on the effects of temperature on the activity mediated phase separation.

preprint2019arXiv

Initial Correlation Dependence of Aging in Phase Separating Solid Binary Mixtures and Ordering Ferromagnets

Following quenches of initial configurations having long range spatial correlations, prepared at the demixing critical point, to points inside the miscibility gap, we study aging phenomena in solid binary mixtures. Results on the decay of the two-time order-parameter autocorrelation functions, obtained from Monte Carlo simulations of the two-dimensional Ising model, with Kawasaki exchange kinetics, are analyzed via state-of-the art methods. The outcome is compared with that obtained for the ordering in uniaxial ferromagnets. For the latter, we have performed Monte Carlo simulations of the same model using the Glauber mechanism. For both types of systems we provide comparative discussion of our results with reference to those concerning quenches with configurations having no spatial correlation. We also discuss the role of structure on the decay of these correlations.

preprint2016arXiv

Fractality in Persistence Decay and Domain Growth during Ferromagnetic Ordering: Dependence upon initial correlation

Dynamics of ordering in Ising model, following quench to zero temperature, have been studied via Glauber spin-flip Monte Carlo simulations in space dimensions $d=2$ and $3$. One of the primary objectives has been to understand phenomena associated with the persistent spins, viz., time decay in the number of unaffected spins, growth of the corresponding pattern and its fractal dimensionality, for varying correlation length in the initial configurations, prepared at different temperatures, at and above the critical value. It is observed that the fractal dimensionality and the exponent describing the power-law decay of persistence probability are strongly dependent upon the relative values of nonequilibrium domain size and the initial equilibrium correlation length. Via appropriate scaling analyses, these quantities have been estimated for quenches from infinite and critical temperatures. The above mentioned dependence is observed to be less pronounced in higher dimension. In addition to these findings for the local persistence, we present results for the global persistence as well. Further, important observations on the standard domain growth problem are reported. For the latter, a controversy in $d=3$, related to the value of the exponent for the power-law growth of the average domain size with time, has been resolved.

preprint2016arXiv

Pattern, Growth and Aging in a Colony of Clustering Active Swimmers

Via molecular dynamics simulations, we study the kinetics in a phase separating active matter model. Quantitative results for the isotropic bicontinuous pattern formation, its growth and aging, studied, respectively, via the two-point equal-time density-density correlation function, the average domain length and the two-time density autocorrelation function, are presented. Both the correlation functions exhibit basic scaling properties, implying self-similarity in the pattern dynamics, for which the average domain size exhibits a power-law growth in time. The equal-time correlation has a short distance behavior that provides reasonable agreement of the corresponding structure factor tail with the Porod law. The autocorrelation decay is a power-law in the average domain size. Apart from these basic similarities, the quantitative behavior of the above mentioned observables are found to be vastly different from those of the corresponding passive limit of the model which also undergoes phase separation. The functional forms of these have been quantified. An exceptionally rapid growth in the active system occurs due to fast coherent motion of the particles, mean-squared-displacements of which exhibit multiple scaling regimes, including a long time ballistic one.

preprint2015arXiv

Dimensionality Dependence of Aging in Kinetics of Diffusive Phase Separation: Behavior of order-parameter autocorrelation

Behavior of two-time autocorrelation during the phase separation in solid binary mixtures are studied via numerical solutions of the Cahn-Hilliard equation as well as Monte Carlo simulations of the Ising model. Results are analyzed via state-of-the-art methods, including the finite-size scaling technique. Full forms of the autocorrelation in space dimensions $2$ and $3$ are obtained empirically. The long time behavior are found to be power-law type, with exponents unexpectedly higher than the ones for the ferromagnetic ordering. Both Chan-Hilliard and Ising models provide results consistent with each other.

preprint2015arXiv

Kinetics of Fluid Phase Separation

We review understanding of kinetics of fluid phase separation in various space dimensions. Morphological differences, percolating or disconnected, based on overall composition in a binary liquid or density in a vapor-liquid system, have been pointed out. Depending upon the morphology, various possible mechanisms and corresponding theoretical predictions for domain growth are discussed. On computational front, useful models and simulation methodologies have been presented. Theoretically predicted growth laws have been tested via molecular dynamics simulations of vapor-liquid transitions. In case of disconnected structure, the mechanism has been confirmed directly. This is a brief review on the topic for a special issue on coarsening dynamics, expected to appear in Comptes Rendus Physique.

preprint2015arXiv

Persistence in Ferromagnetic Ordering: Dependence upon initial configuration

We study the dynamics of ordering in ferromagnets via Monte Carlo simulations of the Ising model, employing the Glauber spin-flip mechanism, in space dimensions $d=2$ and $3$. Results for the persistence probability and the domain growth are discussed for quenches to various temperatures ($T_f$) below the critical one ($T_{c}$), from different initial temperatures $T_{i} \geq T_{c}$. In long time limit, for $T_{i} > T_{c}$, the persistence probability exhibits power-law decay with exponents $θ\simeq 0.22$ and $\simeq 0.18$ in $d=2$ and $3$, respectively. For finite $T_i$, the early time behavior is a different power-law whose life-time diverges and exponent decreases as $T_{i} \rightarrow T_{c}$. The crossover length between the two steps diverges as the equilibrium correlation length. $T_i=T_c$ is expected to provide a {\it{new universality class}} for which we obtain $θ\simeq 0.035$ in $d=2$ and $\simeq 0.10$ in $d=3$. The time dependence of the average domain size $\ell$, however, is observed to be rather insensitive to the choice of $T_i$.

preprint2014arXiv

Phase Behavior of Active Swimmers in Depletants: Molecular Dynamics and Integral Equation Theory

We study the structure and phase behavior of a binary mixture where one of the components is self-propelling in nature. The inter-particle interactions in the system were taken from the Asakura-Oosawa model, for colloid-polymer mixtures, for which the phase diagram is known. In the current model version the colloid particles were made active using the Vicsek model for self-propelling particles. The resultant active system was studied by molecular dynamics methods and integral equation theory. Both methods produce results consistent with each other and demonstrate that the Vicsek model based activity facilitates phase separation, thus broadening the coexistence region.

preprint2013arXiv

Effects of Density Conservation and Hydrodynamics on Aging in Nonequilibrium Processes

Aging in kinetics of three different phase transitions, viz., magnetic, binary solid and single component fluid, have been studied via Monte Carlo and molecular dynamics simulations in three space dimensions with the objective of identifying the effects of order-parameter conservation and hydrodynamics. We observe that the relevant autocorrelations exhibit power-law decay in ferromagnet and binary solid but with different exponents. At early time fluid autocorrelation function nicely follows that of binary solid, the order parameter being conserved for both of them, as opposed to a ferromagnet. At late time the fluid data crosses over to an exponential decay which we identify as a hydrodynamic effect and provide analytical justification for this behavior.

preprint2013arXiv

Temperature and Composition Dependence of Kinetics of Phase Separation in Solid Binary Mixtures

We present results for the kinetics of phase separation in solid binary mixtures ($A_1+A_2$) from Monte Carlo simulations of the Ising model in two dimensions. The simulation results are understood via appropriate application of the finite-size scaling theory. At moderately high temperatures, for symmetric $(50:50)$ compositions of $A_1$ and $A_2$ particles the average size of the domains exhibit power-law growth with the exponent having the Lifshitz-Slyozov value of 1/3 from very early time. However, our analysis shows that for low enough temperatures, the growth exponent at early time is smaller than the Lifshitz-Slyozov value. For composition dependence, we find that at moderate temperature, even for extreme off-critical composition, curvature dependent correction to the growth law is weak which is counterintuitive in this case that gives rise to droplet morphology. This is however consistent with the recent understanding on curvature dependence of surface tension. Results from rather general studies on the finite-size effects with the variations of temperature and composition have also been presented.

preprint2012arXiv

Aging and Crossovers in Phase-Separating Fluid Mixtures

We use state-of-the-art molecular dynamics simulations to study hydrodynamic effects on aging during kinetics of phase separation in a fluid mixture. The domain growth law shows a crossover from a diffusive regime to a viscous hydrodynamic regime. There is a corresponding crossover in the autocorrelation function from a power-law behavior to an exponential decay. While the former is consistent with theories for diffusive domain growth, the latter results as a consequence of faster advective transport in fluids for which an analytical justification has been provided.

preprint2012arXiv

Surface-Directed Spinodal Decomposition: A Molecular Dynamics Study

We use molecular dynamics (MD) simulations to study surface-directed spinodal decomposition (SDSD) in unstable binary ($AB$) fluid mixtures at wetting surfaces. The thickness of the wetting layer $R_1$ grows with time $t$ as a power-law ($R_1 \sim t^θ$). We find that hydrodynamic effects result in a crossover of the growth exponent from $θ\simeq 1/3$ to $θ\simeq1$. We also present results for the layer-wise correlation functions and domain length scales.

preprint2011arXiv

Diffusive Domain Coarsening : Early Time Dynamics and Finite-Size Effects

We study diffusive dynamics of phase separation in a binary mixture, following critical quench, both in spatial dimensions $d=2$ and $d=3$. Particular focus in this work is to obtain information about effects of system size and correction to the growth law via appropriate application of finite-size scaling method to the results obtained from Kawasaki exchange Monte Carlo simulation of Ising model. Observations of only weak size effects and very small correction to scaling in the growth law are significant. The methods used in this work and information thus gathered will be of paramount importance in the study of kinetics of phase separation in fluids and other problems of growing length scale. We also provide detailed discussion on standard methods of understanding simulation results which may lead to inappropriate conclusions.

preprint2011arXiv

Monte Carlo Methods for Estimating Interfacial Free Energies and Line Tensions

Excess contributions to the free energy due to interfaces occur for many problems encountered in the statistical physics of condensed matter when coexistence between different phases is possible (e.g. wetting phenomena, nucleation, crystal growth, etc.). This article reviews two methods to estimate both interfacial free energies and line tensions by Monte Carlo simulations of simple models, (e.g. the Ising model, a symmetrical binary Lennard-Jones fluid exhibiting a miscibility gap, and a simple Lennard-Jones fluid). One method is based on thermodynamic integration. This method is useful to study flat and inclined interfaces for Ising lattices, allowing also the estimation of line tensions of three-phase contact lines, when the interfaces meet walls (where "surface fields" may act). A generalization to off-lattice systems is described as well. The second method is based on the sampling of the order parameter distribution of the system throughout the two-phase coexistence region of the model. Both the interface free energies of flat interfaces and of (spherical or cylindrical) droplets (or bubbles) can be estimated, including also systems with walls, where sphere-cap shaped wall-attached droplets occur. The curvature-dependence of the interfacial free energy is discussed, and estimates for the line tensions are compared to results from the thermodynamic integration method. Basic limitations of all these methods are critically discussed, and an outlook on other approaches is given.

preprint2011arXiv

Universality in Fluid Domain Coarsening: The case of vapor-liquid transition

Domain growth during the kinetics of phase separation is studied following vapor-liquid transition in a single component Lennard-Jones fluid. Results are analyzed after appropriately mapping the continuum snapshots obtained from extensive molecular dynamics simulations to a simple cubic lattice. For near critical quench interconnected domain morphology is observed. A brief period of slow diffusive growth is followed by a linear viscous hydrodynamic growth that lasts for an extended period of time. This result is in contradiction with earlier inclusive reports of late time growth exponent 1/2 that questions the uniqueness of the non-equilibrium universality for liquid-liquid and vapor-liquid transitions.

preprint2010arXiv

Does Young's equation hold on the nanoscale? A Monte Carlo test for the binary Lennard-Jones fluid

When a phase-separated binary ($A+B$) mixture is exposed to a wall, that preferentially attracts one of the components, interfaces between A-rich and B-rich domains in general meet the wall making a contact angle $θ$. Young's equation describes this angle in terms of a balance between the $A-B$ interfacial tension $γ_{AB}$ and the surface tensions $γ_{wA}$, $γ_{wB}$ between, respectively, the $A$- and $B$-rich phases and the wall, $γ_{AB} \cos θ=γ_{wA}-γ_{wB}$. By Monte Carlo simulations of bridges, formed by one of the components in a binary Lennard-Jones liquid, connecting the two walls of a nanoscopic slit pore, $θ$ is estimated from the inclination of the interfaces, as a function of the wall-fluid interaction strength. The information on the surface tensions $γ_{wA}$, $γ_{wB}$ are obtained independently from a new thermodynamic integration method, while $γ_{AB}$ is found from the finite-size scaling analysis of the concentration distribution function. We show that Young's equation describes the contact angles of the actual nanoscale interfaces for this model rather accurately and location of the (first order) wetting transition is estimated.

preprint2010arXiv

Domain Coarsening in 2-d Ising Model: Finite-Size Scaling for Conserved Dynamics

We quantify the effect of system size in the kinetics of domain growth in Ising model with 50:50 composition in two spatial dimensions. Our estimate of the exponent, $α=0.334\pm0.004$, for the power law growth of linear domain size, from Monte Carlo simulation using small systems of linear dimensions L=16, 32, 64, and 128, is in excellent agreement with the prediction of Lifshitz-Slyozov (LS) theory, $α=1/3$. We find that the LS exponent sets in very early and continues to be true until average size of domains reaches three quarters of equilibrium limit.

preprint2010arXiv

Hydrodynamic Crossovers in Surface-Directed Spinodal Decomposition and Surface Enrichment

We present comprehensive molecular dynamics (MD) results for the kinetics of surface-directed spinodal decomposition (SDSD) and surface enrichment (SE) in binary mixtures at wetting surfaces. We study the surface morphology and the growth dynamics of the wetting and enrichment layers. The growth law for the thickness of these layers shows a crossover from a diffusive regime to a hydrodynamic regime. We provide phenomenological arguments to understand this crossover.

preprint2010arXiv

Kinetics of Phase Separation in Fluids: A Molecular Dynamics Study

We present results from extensive 3-d molecular dynamics (MD) simulations of phase separation kinetics in fluids. A coarse-graining procedure is used to obtain state-of-the-art MD results. We observe an extended period of temporally linear growth in the viscous hydrodynamic regime. The morphological similarity of coarsening in fluids and solids is also quantified. The velocity field is characterized by the presence of monopole-like defects, which yield a generalized Porod tail in the corresponding structure factor.

preprint2010arXiv

Kinetics of Surface Enrichment: A Molecular Dynamics Study

We use molecular dynamics (MD) to study the kinetics of surface enrichment (SE) in a stable homogeneous mixture (AB), placed in contact with a surface which preferentially attracts A. The SE profiles show a characteristic double-exponential behavior with two length scales: ξ_-, which rapidly saturates to its equilibrium value, and ξ_+, which diverges as a power-law with time (ξ_+ \sim t^θ). We find that hydrodynamic effects result in a crossover of the growth exponent from θ\simeq 0.5 to θ\simeq 1.0. There is also a corresponding crossover in the growth dynamics of the SE-layer thickness.