Surface-directed spinodal decomposition in binary fluid mixtures

Surface-directed spinodal decomposition of fluids confined in a cylindrical pore

The surface-directed spinodal decomposition of a binary liquid confined inside a cylindrical pore is investigated using molecular dynamics simulations. One component of the liquid wets the pore surface while the other remains neutral. A variety of wetting conditions are studied. For the partial wetting case, after an initial period of phase separation, the domains organize themselves into pluglike structures and the system enters into a metastable state. Therefore, a complete phase separation is never achieved. Analysis of domain growth and the structure factor suggests a one-dimensional growth dynamics for the partial wetting case. As the wetting interaction is increased beyond a critical value, a transition from the pluglike to tubelike domain formation is observed, which corresponds to the full wetting morphology. Thus, a complete phase separation is achieved as the wetting species moves towards the pore surface and forms layers enclosing the nonwetting species residing around the axis of the cylinder. The coarsening dynamics of both the species are studied separately. The wetting species is found to follow a two-dimensional domain growth dynamics with a growth exponent 1/2 in the viscous hydrodynamic regime. This was substantiated by the Porod tail of the structure factor. On the other hand, the domain grows linearly with time for the nonwetting species. This suggests that the nonwetting species behaves akin to a three-dimensional bulk system. An appropriate reasoning is presented to justify the given observations.

Bicontinuous oxide heteroepitaxy with enhanced photoconductivity

Self-assembled systems have recently attracted extensive attention because they can display a wide range of phase morphologies in nanocomposites, providing a new arena to explore novel phenomena. Among these morphologies, a bicontinuous structure is highly desirable based on its high interface-to-volume ratio and 3D interconnectivity. A bicontinuous nickel oxide (NiO) and tin dioxide (SnO₂) heteroepitaxial nanocomposite is revealed here. By controlling their concentration, we fabricated tuneable self-assembled nanostructures from pillars to bicontinuous structures, as evidenced by TEM-energy-dispersive X-ray spectroscopy with a tortuous compositional distribution. The experimentally observed growth modes are consistent with predictions by first-principles calculations. Phase-field simulations are performed to understand 3D microstructure formation and extract key thermodynamic parameters for predicting microstructure morphologies in SnO₂:NiO nanocomposites of other concentrations. Furthermore, we demonstrate significantly enhanced photovoltaic properties in a bicontinuous SnO₂:NiO nanocomposite macroscopically and microscopically. This research shows a pathway to developing innovative solar cell and photodetector devices based on self-assembled oxides.

Universal fast mode regime in wetting kinetics

We present simulation results from a comprehensive molecular dynamics (MD) study of surface-directed spinodal decomposition (SDSD) in unstable symmetric binary mixtures at wetting surfaces. We consider long-ranged and short-ranged surface fields to investigate the early stage wetting kinetics. The attractive part of the long-ranged potential is of the form V(z)∼z^{-n}, where z is the distance from the surface and n is the power-law exponent. We find that the wetting-layer thickness R_{1}(t) at very early times exhibits a power-law growth with an exponent α=1/(n+2). It then crosses over to a universal fast-mode regime with α=3/2. In contrast, for the short-ranged surface potential, a logarithmic behavior in R_{1}(t) is observed at initial times. Remarkably, similar rapid growth is seen in this case too. We provide phenomenological arguments to understand these growth laws. Our MD results firmly establish the existence of universal fast-mode kinetics and settle the related controversy.

Microenvironment-mediated cancer dormancy: Insights from metastability theory

Dormancy is an evolutionarily conserved protective mechanism widely observed in nature. A pathological example is found during cancer metastasis, where cancer cells disseminate from the primary tumor, home to secondary organs, and enter a growth-arrested state, which could last for decades. Recent studies have pointed toward the microenvironment being heavily involved in inducing, preserving, or ceasing this dormant state, with a strong focus on identifying specific molecular mechanisms and signaling pathways. Increasing evidence now suggests the existence of an interplay between intracellular as well as extracellular biochemical and mechanical cues in guiding such processes. Despite the inherent complexities associated with dormancy, proliferation, and growth of cancer cells and tumor tissues, viewing these phenomena from a physical perspective allows for a more global description, independent from many details of the systems. Building on the analogies between tissues and fluids and thermodynamic phase separation concepts, we classify a number of proposed mechanisms in terms of a thermodynamic metastability of the tumor with respect to growth. This can be governed by interaction with the microenvironment in the form of adherence (wetting) to a substrate or by mechanical confinement of the surrounding extracellular matrix. By drawing parallels with clinical and experimental data, we advance the notion that the local energy minima, or metastable states, emerging in the tissue droplet growth kinetics can be associated with a dormant state. Despite its simplicity, the provided framework captures several aspects associated with cancer dormancy and tumor growth.

Impact of the prequench state of binary fluid mixtures on surface-directed spinodal decomposition

Using lattice Boltzmann simulations we investigate the impact of the amplitude of concentration fluctuations in binary fluid mixtures prior to demixing when in contact with a surface that is preferentially wet by one of the components. We find a bicontinuous structure near the surface for an initial, prequench state of the mixture close to the critical point where the amplitude of concentration fluctuations is large. In contrast, if the initial state of the mixture is not near the critical point and concentration fluctuations are relatively weak, then the morphology is not bicontinuous but remains layered until the very late stages of coarsening. In both cases, it is the morphology of a depletion layer rich in the nonpreferred component that dictates the growth exponent of the thickness of the fluid layer that is in direct contact with the substrate. In the early stages of demixing, we find a growth exponent consistent with a value of 1/4 for a prequench state away from the critical point, which is different from the usual diffusive scaling exponent of 1/3 that we recover for a prequench state close to the critical point. We attribute this to the structure of a depletion layer that is penetrated by tubes of the preferred fluid, connecting the wetting layer to the bulk fluid even in the early stages if the initial state is characterized by concentration fluctuations that are large in amplitude. Furthermore, we find that in the late stages of demixing the flow through these tubes results in significant in-plane concentration variations near the substrate, leading to dropletlike structures with a concentration lower than the average concentration in the wetting layer. This causes a deceleration in the growth of the wetting layer in the very late stages of the demixing. Irrespective of the prequench state of the mixture, the late stages of the demixing process produce the same growth law for the layer thickness, with a scaling exponent of unity usually associated with the impact of hydrodynamic flow fields.

Surface-directed spinodal decomposition on morphologically patterned substrates

This paper is the second in a two-part exposition on surface-directed spinodal decomposition (SDSD), i.e., the interplay of kinetics of wetting and phase separation at a surface which is wetted by one of the components of a binary mixture. In our first paper [P. Das, P. K. Jaiswal, and S. Puri, Phys. Rev. E 102, 012803 (2020)2470-004510.1103/PhysRevE.102.012803], we studied SDSD on chemically heterogeneous and physically flat substrates. In this paper, we study SDSD on a chemically homogeneous but morphologically patterned substrate. Such substrates arise in a vast variety of technological applications. Our goal is to provide a theoretical understanding of SDSD in this context. We present detailed numerical results for domain growth both inside and above the grooves in the substrate. The morphological evolution can be understood in terms of the interference of SDSD waves originating from the different surfaces comprising the substrate.

Two-dimensional bicontinuous structures from symmetric surface-directed spinodal decomposition in thin films

We present numerical simulations of symmetric surface-directed spinodal decomposition in thin films with varying film thickness and film composition. The simulations utilize a Cahn-Hilliard model to describe phase separation kinetics in confined film geometries. The systems consist of two phases: a wetting phase that completely wets the top and bottom surfaces, and a nonwetting phase. Three distinct morphologies emerge including a discrete nonwetting morphology, a discrete wetting morphology, as well as a unique two-dimensional bicontinuous morphology that forms for specific values of film thickness and composition. The morphologies are analyzed with a Hoshen-Kopelman algorithm to quantify the degree of continuity of the nonwetting phase, and a morphology map is presented to guide future work.

Surface-directed spinodal decomposition: a molecular dynamics study

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

Kinetics of phase separation in thin films: lattice versus continuum models for solid binary mixtures

A description of phase separation kinetics for solid binary (A,B) mixtures in thin film geometry based on the Kawasaki spin-exchange kinetic Ising model is presented in a discrete lattice molecular field formulation. It is shown that the model describes the interplay of wetting layer formation and lateral phase separation, which leads to a characteristic domain size l(t) in the directions parallel to the confining walls that grows according to the Lifshitz-Slyozov t;{13} law with time t after the quench. Near the critical point of the model, the description is shown to be equivalent to the standard treatments based on Ginzburg-Landau models. Unlike the latter, the present treatment is reliable also at temperatures far below criticality, where the correlation length in the bulk is only of the order of a lattice spacing, and steep concentration variations may occur near the walls, invalidating the gradient square approximation. A further merit is that the relation to the interaction parameters in the bulk and at the walls is always transparent, and the correct free energy at low temperatures is consistent with the time evolution by construction.

Transport theory at the nanoscale. II. Interface dynamics for film growth

Conditions for surface nucleation and growth of a film are determined in a diffuse interface model. A method is given, derived from a Fokker-Planck equation for the nonequilibrium particle distribution, which links atomic and mesoscopic events in a rheological description similar to the classical continuum theory of fluid flow. Film nucleation and growth are modeled by the spatially inhomogeneous continuous evolution of the instantaneous density profile which measures the average number of particles or molecules at given time and position. It is shown how an alteration in the distribution of particles in the vicinity of the boundary between parent and product phases induces transient film growth and damped vibrations at the surface. The method is general but as an illustration, the condensation of a simple classical fluid on cooling is considered in detail.

Molecular dynamics study of phase separation kinetics in thin films

We use molecular dynamics to simulate experiments where a symmetric binary fluid mixture (AB), confined between walls that preferentially attract one component (A), is quenched from the one-phase region into the miscibility gap. Surface enrichment occurs during the early stages, yielding a B-rich mixture in the film center with well-defined A-rich droplets. The droplet size grows with time as l(t) proportional t(2/3) after a transient regime. The present atomistic model is also compared to mesoscopic coarse-grained models for this problem.

Kinetics of phase separation in thin films: simulations for the diffusive case

We study the diffusion-driven kinetics of phase separation of a symmetric binary mixture (AB), confined in a thin-film geometry between two parallel walls. We consider cases where (i) both walls preferentially attract the same component (A), and (ii) one wall attracts and the other wall attracts (with the same strength). We focus on the interplay of phase separation and wetting at the walls, which is referred to as surface-directed spinodal decomposition (SDSD). The formation of SDSD waves at the two surfaces, with wave vectors oriented perpendicular to them, often results in a metastable layered state (also referred to as "stratified morphology"). This state is reminiscent of the situation where the thin film is still in the one-phase region but the surfaces are completely wet, and hence coated with thick wetting layers. This metastable state decays by spinodal fluctuations and crosses over to an asymptotic growth regime characterized by the lateral coarsening of pancakelike domains. These pancakes may or may not be coated by precursors of wetting layers. We use Langevin simulations to study this crossover and the growth kinetics in the asymptotic coarsening regime.

Early-stage compositional segregation in polymer-blend films

The existence of a transient period during the surface enrichment of a binary polymer blend by one of its components has been suggested by previous theoretical and experimental studies as well as computer simulations. Taking advantage of the high depth resolution of neutron reflectivity and the slow dynamics of polymers near their glass transition, we investigate this early-stage surface compositional enrichment in a phase separating polymer blend for the first time. Two stages of surface enrichment layer growth are observed. A rapid local surface enrichment at the chain segmental level occurs first, followed by a slower growth of a diffuse layer having a scale on the order of the bulk correlation length and the radius of gyration of the surface enriching polymer chains.

Hydrodynamics of binary fluid phase segregation

Starting with the Vlasov-Boltzmann equation for a binary fluid mixture, we derive an equation for the velocity field u when the system is segregated into two phases (at low temperatures) with a sharp interface between them. u satisfies the incompressible Navier-Stokes equations together with a jump boundary condition for the pressure across the interface which, in turn, moves with a velocity given by the normal component of u. Numerical simulations of the Vlasov-Boltzmann equations for shear flows parallel and perpendicular to the interface in a phase segregated mixture support this analysis. We expect similar behavior in real fluid mixtures.

Surface-directed phase separation with off-critical composition: analytical and numerical results

We study the interplay of wetting and phase separation in an unstable binary mixture (AB) with off-critical composition, placed in contact with a surface which prefers the component A. We consider surface potentials V(z) approximately z(-n), where z is the distance from the surface, and present analytical arguments and detailed numerical results to elucidate wetting-layer kinetics for arbitrary mixture compositions. If the preferred component is the minority phase, the wetting-layer thickness exhibits a potential-specific behavior at early times tau, R1 approximately tau(1/(n+2)), before crossing over to the universal growth law, R1 approximately tau(1/3). On the other hand, if the preferred component is the majority phase, there is a crossover from potential-specific growth (as before) to a slower growth regime.

Using patterned substrates to promote mixing in microchannels

Using a lattice Boltzmann model for fluid dynamics, we investigate the flow and phase behavior of a binary fluid moving over a patterned substrate within a microchannel. The binary fluid consists of two immiscible components, A and B, and this liquid is subjected to a Poiseuille flow. The substrate is decorated with a checkerboard pattern of A- and B-like patches. Through a coupling of hydrodynamics and thermodynamics, each component is driven to flow from the nonwettable domains to wettable regions. As a consequence, the A and B fluids undergo extensive mixing within the microchannels. We investigate how the degree of mixing depends on the size of the patches, the velocity of the imposed flow field, and the characteristics of the fluid. The results provide guidelines for creating localized "mixing stations" within microfluidic devices. The findings also reveal how a combination of imposed flow fields and surface patterning can be exploited to control the phase behavior of complex fluids.