Spinodal decomposition

Cahn–Hilliard model for spinodal decomposition

Free energies in the presence of small amplitude fluctuations, e.g. in concentration, can be evaluated using an approximation introduced by Ginzburg and Landau to describe magnetic field gradients in superconductors. This approach allows one to approximate the free energy as an expansion in terms of the concentration gradient ${\displaystyle \mathrm{\nabla }c}$, a vector. Since free energy is a scalar and we are probing near its minima, the term proportional to ${\displaystyle \mathrm{\nabla }c}$ is negligible. The lowest order term is the quadratic expression ${\displaystyle \kappa (\mathrm{\nabla }c{)}^2}$, a scalar. Here ${\displaystyle \kappa }$ is a parameter that controls the free energy cost of variations in concentration ${\displaystyle c}$.

The Cahn–Hilliard free energy is then

${\displaystyle F={\int }_v[f_b+\kappa (\mathrm{\nabla }c{)}^2]dV}$

where ${\displaystyle f_b}$ is the bulk free energy per unit volume of the homogeneous solution, and the integral is over the volume of the system.

We now want to study the stability of the system with respect to small fluctuations in the concentration ${\displaystyle c}$, for example a sine wave of amplitude ${\displaystyle a}$ and wavevector ${\displaystyle q=2\pi /\lambda }$, for ${\displaystyle \lambda }$ the wavelength of the concentration wave. To be thermodynamically stable, the free energy change ${\displaystyle \delta F}$ due to any small amplitude concentration fluctuation ${\displaystyle \delta c=a\sin (\overrightarrow{q}.\overrightarrow{r})}$, must be positive.

We may expand ${\displaystyle f_b}$ about the average composition c_o as follows:

${\displaystyle f_b(c)=f_b(c_0)+\left(c-c_0\right){\left(\frac{\mathrm{\partial }f}{\mathrm{\partial }c}\right)}_{c=c_0}+\frac{1}{2}{\left(c-c_0\right)}^2{\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}\right)}_{c=c_0}+\cdots }$

and for the perturbation ${\displaystyle \delta c=a\sin (\overrightarrow{q}.\overrightarrow{r})}$ the free energy change is

${\displaystyle f_b+\kappa (\mathrm{\nabla }c{)}^2=f_b(c_0)+a\sin (\overrightarrow{q}.\overrightarrow{r}){\left(\frac{\mathrm{\partial }f}{\mathrm{\partial }c}\right)}_{c=c_0}+\frac{1}{2}{a}^2{\sin }^2(\overrightarrow{q}.\overrightarrow{r}){\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}\right)}_{c=c_0}+a^2\kappa q^2{\cos }^2(\overrightarrow{q}.\overrightarrow{r})}$

When this is integrated over the volume ${\displaystyle V}$, the ${\displaystyle \sin (\overrightarrow{q}.\overrightarrow{r})}$ gives zero, while ${\displaystyle {\sin }^2(\overrightarrow{q}.\overrightarrow{r})}$ and ${\displaystyle {\cos }^2(\overrightarrow{q}.\overrightarrow{r})}$ integrate to give ${\displaystyle V/2}$. So, then

${\displaystyle \frac{\delta F}{V}=\frac{a^2}{4}\left[{\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}\right)}_{c=c_0}+2\kappa q^2\right]}$

As ${\displaystyle a^2>0}$, thermodynamic stability requires that the term in brackets be positive. The ${\displaystyle 2\kappa q^2}$ is always positive but tends to zero at small wavevectors, large wavelengths. Since we are interested in macroscopic fluctuations, ${\displaystyle q\to 0}$, stability requires that the second derivative of the free energy be positive. When it is, there is no spinodal decomposition, but when it is negative, spinodal decomposition will occur. Then fluctuations with wavevectors ${\displaystyle q<q_c}$ become spontaneously unstable, where the critical wave number ${\displaystyle q_c}$ is given by:

${\displaystyle q_c=\sqrt{\frac{-1}{2\kappa }{\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}\right)}_{c=c_0}}}$

which corresponds to a fluctuations above a critical wavelength

${\displaystyle {\lambda }_c=\sqrt{-8{\pi }^2\kappa /{\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}\right)}_{c=c_0}}}$

Dynamics of spinodal decomposition when molecules move via diffusion

Spinodal decomposition can be modeled using a generalized diffusion equation:

${\displaystyle \frac{\mathrm{\partial }c}{\mathrm{\partial }t}=M{\mathrm{\nabla }}^2\mu }$

for ${\displaystyle \mu }$ the chemical potential and ${\displaystyle M}$ the mobility. As pointed out by Cahn, this equation can be considered as a phenomenological definition of the mobility M, which must by definition be positive. It consists of the ratio of the flux to the local gradient in chemical potential. The chemical potential is a variation of the free energy and when this is the Cahn–Hilliard free energy this is

${\displaystyle \mu =\frac{\delta F}{\delta c}={\left(\frac{\mathrm{\partial }f}{\mathrm{\partial }c}\right)}_{c=c_0}-2\kappa {\mathrm{\nabla }}^{2}c}$

and so

${\displaystyle \frac{\mathrm{\partial }c}{\mathrm{\partial }t}=M{\mathrm{\nabla }}^2\mu =M\left[{\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}\right)}_{c=c_0}{\mathrm{\nabla }}^{2}c-2\kappa {\mathrm{\nabla }}^{4}c\right]}$

and now we want to see what happens to a small concentration fluctuation ${\displaystyle \delta c=a\exp (\omega t)\sin (\overrightarrow{q}.\overrightarrow{r})}$ - note that now it has time dependence as a wavevector dependence. Here ${\displaystyle \omega }$ is a growth rate. If ${\displaystyle \omega <0}$ then the perturbation shrinks to nothing, the system is stable with respect to small perturbations or fluctuations, and there is no spinodal decomposition. However, if ${\displaystyle \omega >0}$ then the perturbation grows and the system is unstable with respect to small perturbations or fluctuations: There is spinodal decomposition.

Substituting in this concentration fluctuation, we get

${\displaystyle \omega \delta c=M\left[-{\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}\right)}_{c=c_0}{q}^2-2\kappa q^4\right]\delta c}$

This gives the same expressions for the stability as above, but it also gives an expression for the growth rate of concentration perturbations

${\displaystyle \omega =M{q}^2\left[-{\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}\right)}_{c=c_0}-2\kappa q^2\right]}$

which has a maximum at a wavevector

${\displaystyle q_{\mathrm{m}\mathrm{a}\mathrm{x}}=\sqrt{-{\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}\right)}_{c=c_0}/(4\kappa )}}$

So, at least at the beginning of spinodal decomposition, we expect the growing concentrations to mostly have this wavevector.

Phase diagram

This type of phase transformation is known as spinodal decomposition, and can be illustrated on a phase diagram exhibiting a miscibility gap. Thus, phase separation occurs whenever a material transition into the unstable region of the phase diagram. The boundary of the unstable region sometimes referred to as the binodal or coexistence curve, is found by performing a common tangent construction of the free-energy diagram. Inside the binodal is a region called the spinodal, which is found by determining where the curvature of the free-energy curve is negative. The binodal and spinodal meet at the critical point. It is when a material is moved into the spinodal region of the phase diagram that spinodal decomposition can occur.

The free energy curve is plotted as a function of composition for a temperature below the convolute temperature, T. Equilibrium phase compositions are those corresponding to the free energy minima. Regions of negative curvature (∂²f/∂c² < 0 ) lie within the inflection points of the curve (∂²f/∂c² = 0 ) which are called the spinodes. Their locus as a function of temperature defines the spinodal curve. For compositions within the spinodal, a homogeneous solution is unstable against infinitesimal fluctuations in density or composition, and there is no thermodynamic barrier to the growth of a new phase. Thus, the spinodal represents the limit of physical and chemical stability.

To reach the spinodal region of the phase diagram, a transition must take the material through the binodal region or the critical point. Often phase separation will occur via nucleation during this transition, and spinodal decomposition will not be observed. To observe spinodal decomposition, a very fast transition, often called a quench, is required to move from the stable to the spinodal unstable region of the phase diagram.

In some systems, ordering of the material leads to a compositional instability and this is known as a conditional spinodal, e.g. in the feldspars.

Coherency strains

For most crystalline solid solutions, there is a variation of lattice parameters with the composition. If the lattice of such a solution is to remain coherent in the presence of a composition modulation, mechanical work has to be done to strain the rigid lattice structure. The maintenance of coherency thus affects the driving force for diffusion.

Consider a crystalline solid containing a one-dimensional composition modulation along the x-direction. We calculate the elastic strain energy for a cubic crystal by estimating the work required to deform a slice of material so that it can be added coherently to an existing slab of cross-sectional area. We will assume that the composition modulation is along the x' direction and, as indicated, a prime will be used to distinguish the reference axes from the standard axes of a cubic system (that is, along the <100>).

Let the lattice spacing in the plane of the slab be a_o and that of the undeformed slice a. If the slice is to be coherent after the addition of the slab, it must be subjected to a strain ε in the z' and y' directions which is given by:

${\displaystyle ϵ=\frac{a-a_0}{a_0}}$

In the first step, the slice is deformed hydrostatically in order to produce the required strains to the z' and y' directions. We use the linear compressibility of a cubic system 1 / ( c₁₁ + 2 c₁₂ ) where the c's are the elastic constants. The stresses required to produce a hydrostatic strain of δ are therefore given by:

${\displaystyle {\sigma }_{x^{\prime }}={\sigma }_{y^{\prime }}={\sigma }_{z^{\prime }}}$

The elastic work per unit volume is given by:

${\displaystyle W_E=\frac{1}{2}{\displaystyle \sum _i{\sigma }_i{ϵ}_i}}$

where the ε's are the strains. The work performed per unit volume of the slice during the first step is therefore given by:

${\displaystyle W_E(1)=\frac{3}{2}(c_{11}+2c_{12}){ϵ}^2}$

In the second step, the sides of the slice parallel to the x' direction are clamped and the stress in this direction is relaxed reversibly. Thus, ε_z' = ε_y' = 0. The result is that:

${\displaystyle W_E(2)=\frac{{ϵ}^2(c_{11}+2c_{22})}{2c_{11}}}$

The net work performed on the slice in order to achieve coherency is given by:

${\displaystyle W_E=W_E(1)-W_E(2)}$

or

${\displaystyle W_E=\frac{{ϵ}^2}{2}(c_{11}+2c_{12})\left(3-\left[\frac{c_{11}-2c_{12}}{c_{1^{\prime }{1}^{\prime }}}\right]\right)}$

The final step is to express c_1'1' in terms of the constants referred to the standard axes. From the rotation of axes, we obtain the following:

${\displaystyle c_{1^{\prime }{1}^{\prime }}=c_{11}+2(2c_{44}-c_{11}+c_{12})(l^2{m}^2+m^2{n}^2+l^2{n}^2)}$

where l, m, n are the direction cosines of the x' axis and, therefore the direction cosines of the composition modulation. Combining these, we obtain the following:

${\displaystyle W_E=Y{ϵ}^2}$

${\displaystyle Y=\frac{1}{2}(c_{11}+2c_{12})\left[3-\frac{c_{11}+2c_{12}}{c_{11}+2(2c_{44}-c_{11}+c_{12})(l^2{m}^2+m^2{n}^2+l^2{n}^2)}\right]}$

The existence of any shear strain has not been accounted for. Cahn considered this problem, and concluded that shear would be absent for modulations along <100>, <110>, <111> and that for other directions the effect of shear strains would be small. It then follows that the total elastic strain energy of a slab of cross-sectional area A is given by:

${\displaystyle W_E=4\int Y{ϵ}^{2}dx}$

We next have to relate the strain δ to the composition variation. Let a_o be the lattice parameter of the unstrained solid of the average composition c_o. Using a Taylor series expansion about c_o yields the following:

${\displaystyle a=a_0[1+\eta [c-c_0]+\cdots ]}$

in which

${\displaystyle \eta =\frac{1}{a_0}\frac{da}{dc}+\frac{d\ln a}{dc}}$

where the derivatives are evaluated at c_o. Thus, neglecting higher-order terms, we have:

${\displaystyle ϵ=\frac{a-a_0}{a_0}=\eta (c-c_0)}$

Substituting, we obtain:

${\displaystyle W_E=A\int {\eta }^{2}Y(c-c_0{)}^{2}dx}$

This simple result indicates that the strain energy of a composition modulation depends only on the amplitude and is independent of the wavelength. For a given amplitude, the strain energy W_E is proportional to Y. Consider a few special cases.

For an isotropic material:

${\displaystyle 2c_{44}-c_{11}+c_{12}=0}$

so that:

${\displaystyle Y[\mathrm{i}\mathrm{s}\mathrm{o}]=c_{11}+c_{12}-2\frac{c_{12}^2}{c_{11}}}$

This equation can also be written in terms of Young's modulus E and Poisson's ratio υ using the standard relationships:

${\displaystyle c_{11}=\frac{E(1-\nu )}{(1-2\nu )(1+\nu )}}$

${\displaystyle c_{12}=\frac{E\nu }{(1-2\nu )(1+\nu )}}$

Substituting, we obtain the following:

${\displaystyle Y[\mathrm{i}\mathrm{s}\mathrm{o}]=\frac{E}{1-\nu }}$

For most metals, the left-hand side of this equation

${\displaystyle 2c_{44}-c_{11}+c_{12}}$

is positive, so that the elastic energy will be a minimum for those directions that minimize the term: l²m² + m²n² + l²n². By inspection, those are seen to be <100>. For this case:

${\displaystyle Y[100]=c_{11}+c_{12}-2\frac{c_{12}^2}{c_{11}}}$

the same as for an isotropic material. At least one metal (molybdenum) has an anisotropy of the opposite sign. In this case, the directions for minimum W_E will be those that maximize the directional cosine function. These directions are <111>, and

${\displaystyle Y[111]=\frac{6c_{44}(c_{11}+2c_{12})}{c_{11}+2c_{12}+4c_{44}}}$

As we will see, the growth rate of the modulations will be a maximum in the directions that minimize Y. These directions, therefore, determine the morphology and structural characteristics of the decomposition in cubic solid solutions.

Rewriting the diffusion equation and including the term derived for the elastic energy yields the following:

${\displaystyle F_t=A\int f(c)+\eta Y(c-c_0{)}^2+K{\left(\frac{dc}{dx}\right)}^{2}dx}$

or

${\displaystyle \frac{\mathrm{\partial }c}{\mathrm{\partial }t}=\frac{M}{N_{\nu }}\left([f^{″}+2\eta Y]\frac{d^{2}c}{d{x}^2}-2K\frac{d^{4}c}{d{x}^4}\right)}$

which can alternatively be written in terms of the diffusion coefficient D as:

${\displaystyle \frac{\mathrm{\partial }c}{\mathrm{\partial }t}=\left[1+\frac{2\eta Y}{f^{″}}\right]\frac{d^{2}c}{d{x}^2}-\frac{2KF}{f^{″}}\frac{d^{4}c}{d{x}^4}}$

The simplest way of solving this equation is by using the method of Fourier transforms.

Fourier transform

The motivation for the Fourier transformation comes from the study of a Fourier series. In the study of a Fourier series, complicated periodic functions are written as the sum of simple waves mathematically represented by sines and cosines. Due to the properties of sine and cosine, it is possible to recover the amount of each wave in the sum by an integral. In many cases it is desirable to use Euler's formula, which states that e^2πiθ = cos 2πθ + i sin 2πθ, to write Fourier series in terms of the basic waves e^2πiθ, with the distinct advantage of simplifying many unwieldy formulas.

The passage from sines and cosines to complex exponentials makes it necessary for the Fourier coefficients to be complex-valued. The usual interpretation of this complex number is that it gives both the amplitude (or size) of the wave present in the function and the phase (or the initial angle) of the wave. This passage also introduces the need for negative "frequencies". (E.G. If θ were measured in seconds then the waves e^2πiθ and e^−2πiθ would both complete one cycle per second—but they represent different frequencies in the Fourier transform. Hence, frequency no longer measures the number of cycles per unit time, but is closely related.)

If A(β) is the amplitude of a Fourier component of wavelength λ and wavenumber β = 2π/λ the spatial variation in composition can be expressed by the Fourier integral:

${\displaystyle c-c_0=\int A(\beta )\exp (i\beta x)d\beta }$

in which the coefficients are defined by the inverse relationship:

${\displaystyle A(\beta )=\frac{1}{2\pi }\int (c-c_0)\exp (-i\beta x)dx}$

Substituting, we obtain on equating coefficients:

${\displaystyle \frac{dA(\beta )}{dt}=-\frac{M}{N_{\nu }}[f^{″}+2{\eta }^{2}Y+2Y{\beta }^2]{\beta }^{2}A(\beta )}$

This is an ordinary differential equation that has the solution:

${\displaystyle A(\beta ,t)=A(\beta ,0)\exp [R(\beta )t]}$

in which A(β) is the initial amplitude of the Fourier component of wave wavenumber β and R(β) defined by:

${\displaystyle R(\beta )=-\frac{M}{N_{\nu }}(f^{″}+2\eta Y+2k{\beta }^2){\beta }^2}$

or, expressed in terms of the diffusion coefficient D:

${\displaystyle R(\beta )=-\tilde{D}\left(1+\frac{2{\eta }^{2}Y}{f^{″}}+\frac{2K}{f^{″}}{\beta }^2\right){\beta }^2}$

In a similar manner, the new diffusion equation:

${\displaystyle \frac{\mathrm{\partial }c}{\mathrm{\partial }t}=M\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}{\mathrm{\nabla }}^{2}c-2MK{\mathrm{\nabla }}^{4}c)}$

has a simple sine wave solution given by:

${\displaystyle c-c_0=exp[R\overline{\beta }t]cos\beta \cdot r}$

where ${\displaystyle R(\beta )}$ is obtained by substituting this solution back into the diffusion equation as follows:

${\displaystyle R(\overline{\beta })-M{\beta }^2\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}+2K{\beta }^2\right)}$

For solids, the elastic strains resulting from coherency add terms to the amplification factor ${\displaystyle R(\beta )}$ as follows:

${\displaystyle R(\overline{\beta })=-M{\beta }^2\left(\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{c}^2}+2{\eta }^{2}Y+2K{\beta }^2\right)}$

where, for isotropic solids:

${\displaystyle Y=\frac{E}{1-\nu }}$,

where E is Young's modulus of elasticity, ν is Poisson's ratio, and η is the linear strain per unit composition difference. For anisotropic solids, the elastic term depends on the direction in a manner that can be predicted by elastic constants and how the lattice parameters vary with composition. For the cubic case, Y is a minimum for either (100) or (111) directions, depending only on the sign of the elastic anisotropy.

Thus, by describing any composition fluctuation in terms of its Fourier components, Cahn showed that a solution would be unstable concerning to the sinusoidal fluctuations of a critical wavelength. By relating the elastic strain energy to the amplitudes of such fluctuations, he formalized the wavelength or frequency dependence of the growth of such fluctuations, and thus introduced the principle of selective amplification of Fourier components of certain wavelengths. The treatment yields the expected mean particle size or wavelength of the most rapidly growing fluctuation.

Thus, the amplitude of composition fluctuations should grow continuously until a metastable equilibrium is reached with preferential amplification of components of particular wavelengths. The kinetic amplification factor R is negative when the solution is stable to the fluctuation, zero at the critical wavelength, and positive for longer wavelengths—exhibiting a maximum at exactly ${\displaystyle \sqrt{2}}$ times the critical wavelength.

Consider a homogeneous solution within the spinodal. It will initially have a certain amount of fluctuation from the average composition which may be written as a Fourier integral. Each Fourier component of that fluctuation will grow or diminish according to its wavelength.

Because of the maximum in R as a function of wavelength, those components of the fluctuation with ${\displaystyle \sqrt{2}}$ times the critical wavelength will grow fastest and will dominate. This "principle of selective amplification" depends on the initial presence of these wavelengths but does not critically depend on their exact amplitude relative to other wavelengths (if the time is large compared with (1/R). It does not depend on any additional assumptions, since different wavelengths can coexist and do not interfere with one another.

Limitations of this theory would appear to arise from this assumption and the absence of an expression formulated to account for irreversible processes during phase separation which may be associated with internal friction and entropy production. In practice, frictional damping is generally present and some of the energy is transformed into thermal energy. Thus, the amplitude and intensity of a one-dimensional wave decrease with distance from the source, and for a three-dimensional wave, the decrease will be greater.

Dynamics in k-space

In the spinodal region of the phase diagram, the free energy can be lowered by allowing the components to separate, thus increasing the relative concentration of a component material in a particular region of the material. The concentration will continue to increase until the material reaches the stable part of the phase diagram. Very large regions of material will change their concentration slowly due to the amount of material that must be moved. Very small regions will shrink away due to the energy cost of maintaining an interface between two dissimilar component materials.

To initiate a homogeneous quench a control parameter, such as temperature, is abruptly and globally changed. For a binary mixture of ${\displaystyle A}$-type and ${\displaystyle B}$-type materials, the Landau free-energy

${\displaystyle F=\int \left(\frac{A}{2}{\varphi }^2+\frac{B}{4}{\varphi }^4+\frac{\kappa }{2}{\left(\mathrm{\nabla }\varphi \right)}^2\right)dx.}$

is a good approximation of the free energy near the critical point and is often used to study homogeneous quenches. The mixture concentration ${\displaystyle \varphi ={\rho }_A-{\rho }_B}$ is the density difference of the mixture components, the control parameters which determine the stability of the mixture are ${\displaystyle A}$ and ${\displaystyle B}$, and the interfacial energy cost is determined by ${\displaystyle \kappa }$.

Diffusive motion often dominates at the length-scale of spinodal decomposition. The equation of motion for a diffusive system is

${\displaystyle {\mathrm{\partial }}_t\varphi =\mathrm{\nabla }(m\mathrm{\nabla }\mu +\xi (x)),}$

where ${\displaystyle m}$ is the diffusive mobility, ${\displaystyle \xi (x)}$ is some random noise such that ${\displaystyle ⟨\xi (x)⟩=0}$, and the chemical potential ${\displaystyle \mu }$ is derived from the Landau free-energy:

${\displaystyle \mu =\frac{\delta F}{\delta \varphi }=A\varphi +B{\varphi }^3-\kappa {\mathrm{\nabla }}^2\varphi .}$

We see that if ${\displaystyle A<0}$, small fluctuations around ${\displaystyle \varphi =0}$ have a negative effective diffusive mobility and will grow rather than shrink. To understand the growth dynamics, we disregard the fluctuating currents due to ${\displaystyle \xi }$, linearize the equation of motion around ${\displaystyle \varphi ={\varphi }_{in}}$ and perform a Fourier transform into ${\displaystyle k}$-space. This leads to

${\displaystyle {\mathrm{\partial }}_t\tilde{\varphi }(k,t)=-m((A+3B{\varphi }_{in}^2)k^2+\kappa k^4)\tilde{\varphi }(k,t)=R(k)\tilde{\varphi }(k,t),}$

which has an exponential growth solution:

${\displaystyle \tilde{\varphi }(k,t)=\exp (R(k)t).}$

Since the growth rate ${\displaystyle R(k)}$ is exponential, the fastest growing angular wavenumber

${\displaystyle k_{sp}=\sqrt{\frac{-(A+3B{\varphi }_{in}^2)}{2\kappa }},}$

will quickly dominate the morphology. We now see that spinodal decomposition results in domains of the characteristic length scale called the spinodal length:

${\displaystyle {\lambda }_{sp}=\frac{2\pi }{k_{sp}}=2\pi \sqrt{\frac{2\kappa }{-(A+3B{\varphi }_{in}^2)}}.}$

The growth rate of the fastest-growing angular wave number is

${\displaystyle R(k_{sp})=-m((A+3B{\varphi }_{in}^2)k_{sp}^2+\kappa k_{sp}^4)=\frac{m(A+3B{\varphi }_{in}^2{)}^2}{4\kappa }=\frac{1}{t_{sp}}}$

where ${\displaystyle t_{sp}}$ is known as the spinodal time.

The spinodal length and spinodal time can be used to nondimensionalize the equation of motion, resulting in universal scaling for spinodal decomposition.

Spinodal Architected Materials

Spinodal phase decomposition has been used to generate architected materials by interpreting one phase as solid, and the other phase as void. These spinodal architected materials present interesting mechanical properties, such as high energy absorption, insensitivity to imperfections, superior mechanical resilience, and high stiffness-to-weight ratio. Furthermore, by controlling the phase separation, i.e., controlling the proportion of materials, and/or imposing preferential directions in the decompositions, one can control the density, and preferential directions effectively tuning the strength, weight, and anisotropy of the resulting architected material. Another interesting property of spinodal materials is the capability to seamlessly transition between different classes, orientations, and densities, thereby enabling the manufacturing of effectively multi-material structures.