The Spatiotemporal Oscillations of Order Parameter for Isothermal Model of the Surface-Directed Spinodal Decomposition in Bounded BinaryMixtures

The asymptotical behavior of order parameter in confined binary mixture is considered in one-dimensional geometry. The interaction between bulk and surface forces in the mixture is investigated. Its established conditions are when the bulk spinodal decomposition may be ignored and when the main role in the process of formation of the oscillating asymptotic periodic spatiotemporal structures plays the surface-directed spinodal decomposition which is modelled by nonlinear dynamical boundary conditions.


Introduction
The problem of self-organization for order parameter in binary mixtures, which are confined with two flat walls, is important in physics of condense matter [1,2].At first, for the binary mixtures, this problem was considered by Binder, Puri, and Frisch.The existence of oscillating distributions of the order parameter (on a space variable for each fixed moment time) has been proved [2,3].Binder and Frisch considered also the Cahn-Hilliard (CH) equation for conserved-order parameter with nonlinear dynamical boundary conditions [1] and described a phase diagram for the problem.However, it is unknown what is being an attractor of the considered dynamical binary system, that is, what is being an asymptotic behavior of the solutions of such imposed boundary-value problem as a time t → ∞.
The aim of the present paper is to establish how the attractor of trajectories for this problem is constructed.The achievements of the paper are as follows: (i) we show that the solutions tend to the limit functions as t → ∞ with a finite set of discontinuities which are determined by topological properties of the functions F k , k = 1, 2 in the boundary conditions (in this case, we have the solutions of relaxation non-Fickian diffusion type); (ii) there are limit solutions with infinite (countable) and infinite (uncountable) number of discontinuities (in such case, we have the solutions of preturbulent type and turbulent type; accordingly, the set of points of discontinuities for the solutions of turbulent type is familiar to a Cantor-like set); (iii) there exist the solutions whose limit functions are discontinuities at all points.This solutions take up a particular place among the solutions and they lead to the appearance of the so-called "deterministic randomness," namely, in order to characterize the long-time behavior of the solutions, it is necessary to proceed to a probabilistic description [4].
In such way, we consider the problem which was introduced in [1], but instead of the classical (CH) equation, we research the so-called modified Cahn-Hilliard (mCH) equation which is proposed by Galenko et al. [5,6].The modified equation is the model for diffusion and phase separation which takes into account hyperbolic relaxation of Research Letters in Physics Figure 1: Asymptotically periodic distributions of the order parameter of the relaxation type.The map Φ : g → g 2 + μ.Here, β 1 and β 3 are stable ordered phases, and β 2 is unstable disordered phase.If the amplitude of initial bubbles is greater than β 2 , then there are the oscillations.
If the amplitude of initial bubbles is less than β 2 , then asymptotic distributions of order parameter tend to the constant value β 1 [4].
the solute diffusion flux I.The flux can be determined by the formula (see [5,6]) where O(τ D ) → 0 as τ D → 0, u is the order parameter, τ D is the delay time, M is the atomic mobility, F u is the u-derivative of the free energy F .We denote that 2 , where F 0 is connection energy [6] which is normalized by k B T. Here, k B is the Boltzmann constant, and ζ b is the bulk correlation length.A few advancements were made for strongly nonequilibrium phase separation.In 1986, Galenko and Lebedev [6] generalized the linearized (CH) equation to the case of existence of slowly relaxing variable.Their calculations show that the instability of the system is not standard diffusive type but rather it is controlled by the relaxation of a slow structural variable such as the diffusion flux determined by the formula (1).
In this paper, we consider the case only when the order parameter varies slowly in the bulk, and we determine the condition when this supposition is true.In this case, asymptotically stable equilibrium states of the system may be stabilized by nonlinear boundary conditions which are described by surface-directed spinodal decomposition (SDSD) in the binary system.As a result, the SDSD determines the limit distributions of order parameter which have the form of traveling waves and its graphes are asymptotic periodic functions with finite, infinite (countable), or infinite (uncountable) number of discontinuities (Figure 1).

Formulation of the Problem
We consider binary mixtures with diffusion flux (1) and, additionally, we suppose that the equation u t + I x = 0 is true.Then, from (1), we obtain the (mCH) equation: In [6], it is shown that (2) takes place for bulk spinodal decompositions with locally nonequilibrium diffusion.In what follows, we consider the boundary conditions in the common form [1]: for the solutions of (2).Here, x 1 = 0, x 2 = l and F k are given functions, h k describes surface attractive potentials (for magnetic liquid, it may be surface magnetic field [1]), and g k = J (k) s /J (k) , where J (•) is modelled the bulk forces and J (•)   s is modelled surface forces (in detail see [1,2]).Next, ν is a constant (or a function) which is proportional to the bulk correlation length ζ b .In [1], for example, it was considered (3) for In the simplest case, if F 1 = h 1 + g 1 u(0, t), then u(0, t) = u(0, 0) exp(g 1 t)−h 1 /g 1 .Thus, if g 1 > 0, then the surface order parameter increases in time, consequently, and the wall x = 0 is attractive.If g 1 < 0, the wall x = 0 is repulsive.In [2], it was considered (3) for where τ s and τ v are the surface and bulk relaxation times accordingly.In typical cases τ s τ v [1].Further, the self-organized structures (as in Figure 1) appear if and only if the Pfaffian form F 2 du 1 − F 1 du 2 = 0 is integrable, that is, it can be reduced to the form dW(u 1 , u 2 ) = 0, where W is the known function, and u k = u(x k , t).From the last equality, it follows that where Φ μ is some function.So, if the Pfaffian form is integrable, then (3) may be reduced to (4).Next, we consider the initial conditions u(x, 0) = u (0) (x), u t (x, 0) = u (1) (x).
In common case, analytical solutions of this problem are unknown; but the boundary-value problem may be reduced to some simpler problem for spatial functions F k and constant mobility M. The problem has unique solution such as in Figure 1.We confined ourselves to the polymerpolymer systems with the components A, B and the power of polymerization N A , N B , accordingly.The free energy (see [5]) is where u A + u B = 1, and κ is Flory-Huggins parameter.We denote u A = u, u B = 1 − u, and we consider, for simplicity, the case N A = N B = N.Then, we find from ( 5) by the linear approximation that at a neighborhood of u = 1/2.Then, perturbations u → u − 1/2 of the equilibrium state u = 1/2 can be considered as the order parameter.As a result, due to (6), equation ( 2) can be rewritten in the form Here, l is the spatial length, τ D = Mα j , α J is the coefficient specifying non-Fickian diffusion (see [6]).So, τ D is the relaxation time of the diffusion flux to its steady state.Further, we suppose that D 2 1 and D 4 1, in particular, it is possible when M(u) ∼ u(1 − u).Hence, we get from (7) that Next, we consider the boundary conditions for the solutions of (8), where h k and g k are some constant (in common case, they may be given functions).Moreover, we suppose that

Phase Diagram
The phase diagram for the solutions of ( 20) is shown on Figure 2.This diagram is similar to corresponding diagram which was obtained by numeric simulation in [1,2], where it was considered similar boundary problem to (CH) equation [2].

2 Figure 2 :
Figure 2: Phase diagram is obtained for the problem with nonlinear boundary conditions, which cause the surface-directed spinodal decomposition, by the theoretical method.Here (a) is the region of unique stable phase; (b) is the region of the spinodal decomposition; (c) is the region of the lamellar phase; (d) is the region of stochastic distributions of the order parameter.