---
title: Anisotropic Spinodal Infills
url: https://www.emergentmind.com/topics/anisotropic-spinodal-infills
type: topic
---

# Anisotropic Spinodal Infills

Searching arXiv for recent and relevant papers on anisotropic spinodal structures, spinodoids, and related phase-field/mechanics work.
Anisotropic spinodal infills are spinodal-decomposition-derived internal structures in which phase-separated, shell-based, or compositionally modulated regions are directionally organized. In one formulation, they are “non-uniform, directionally structured networks” of amorphous polytetrahedral aggregates and bcc crystallites generated during liquid–solid spinodal decomposition [1012.4331]. In another, they are “self-assembled, compositionally modulated channels and interfaces” in spinodal superlattices of topological insulators [1803.06289]. In architected materials, they appear as stochastic, aperiodic, bicontinuous “spinodoid” microarchitectures whose directional stiffness, orientation, continuity, and anisotropy are explicitly controlled or optimized [2012.15744, 2506.23420]. Across these settings, the common theme is that a spinodal instability sets an internal length scale, while kinetics, elasticity, diffusion, geometry, or optimization determines whether the resulting infill is isotropic only in the statistical sense or strongly anisotropic in local structure and effective response.

## 1. Thermodynamic and kinetic basis

A central description of liquid–solid spinodal decomposition uses a Landau free energy functional for a scalar order parameter $\psi(\mathbf{r})$, the deviation of particle density from the uniform liquid density,
\[
\Phi[\psi] = \int dV \left[ W(\psi) + D (\nabla \psi)^2  + b \int dV' \frac{\psi(\mathbf{r}) \, \psi(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|} \right],
\]
with
\[
W(\psi) = -\tau \psi^2 + \nu \psi^3 + \lambda \psi^4.
\]
The dynamics is purely relaxational, conserving the order parameter, and is modeled by a Cahn–Hilliard equation,
\[
\frac{\partial \psi}{\partial t} = M \nabla^2 \left( \frac{\delta \Phi}{\delta \psi} \right).
\]
Linearization around the uniform liquid gives the amplification factor
\[
\lambda(k) = -D k^4 + \tau k^2 - b,
\]
so there is a finite band of unstable wavevectors and the most unstable mode occurs near
\[
k_0 \sim \sqrt{\frac{\tau}{D}}
\]
for shallow quenches near the spinodal line $\tau_S = 2\sqrt{bD}$ [1012.4331]. This finite-wavelength instability is the basic source of characteristic spinodal length scales.

In intercalation materials, anisotropy enters directly through elasticity. A phase-field model for Li concentration $c(\mathbf{x},t)$ uses
\[
F[c] = \int_\Omega\left\{f(c) + \frac{\rho}{2}\nabla c \cdot ({\bf K}\nabla c) + \frac{1}{2}{\bf E}:{\bf T}\right\}\mathrm{d}{\bf x},
\]
with eigenstrain
\[
\mathbf{E} = \frac{1}{2}\left(\nabla\mathbf{u} + \nabla\mathbf{u}^{\mathrm{T}}\right) - c\,\mathbf{M},
\qquad
\mathbf{M} = \mathrm{diag}(M_{11},M_{22},M_{33}).
\]
For FePO$_4$/LiFePO$_4$, the mismatch values are
\[
M_{11} = 0.052,\quad M_{22} = 0.036,\quad M_{33} = -0.019,
\]
so insertion produces expansion along two axes and contraction along the third. The resulting growth rate $\sigma(\mathbf{q})$ is direction dependent, and the neutral curve $\sigma(\mathbf{q})=0$ is anisotropic in $(q_x,q_y,q_z)$ [1202.1626].

Anisotropy can also arise from transport. In dynamical density functional theory with a constant diffusion tensor $\mathbf{D}$, the Fourier-space evolution of density fluctuations obeys
\[
\frac{\partial \hat\rho(\mathbf{k},t)}{\partial t}
= - (\mathbf{k}\cdot \mathbf{D}\cdot \mathbf{k})\left[1 - \rho_b \hat c^{(2)}(\mathbf{k})\right]\hat\rho(\mathbf{k},t)
+ \frac{1}{(2\pi)^3}\int d\mathbf{k}'\,(\mathbf{k}\cdot \mathbf{D}\cdot \mathbf{k}')\,\hat\rho(\mathbf{k}-\mathbf{k}',t)\,\hat c^{(2)}(\mathbf{k}')\,\hat\rho(\mathbf{k}',t).
\]
The linear growth rate is
\[
\omega(\mathbf{k}) = - (\mathbf{k} \cdot \mathbf{D}\cdot \mathbf{k})\left[1 - \rho_b \hat c^{(2)}(k)\right].
\]
In this setting, the instability band in $k$ is unchanged, but the rate at which each direction’s mode grows depends on orientation through $\mathbf{k}\cdot\mathbf{D}\cdot\mathbf{k}$, and anisotropic diffusion strongly reduces the coupling between different Fourier modes [1810.01714].

## 2. Morphological classes and the emergence of anisotropy

In liquid–solid spinodal decomposition with an isotropic scalar free energy, the early-stage density field is a superposition of random density waves with wavevectors near $k_0$. For near-critical quenches, the band is very narrow, so the pattern is a strongly correlated filamentary network of density fluctuations, and the two-point correlation function behaves like a Bessel function decaying as
\[
C(r)\sim \frac{1}{r}.
\]
The paper explicitly notes that there is “no built-in directional anisotropy in the free energy,” yet “anisotropic spatial organization” emerges through a filamentary random-wave morphology, surface-driven anisotropy of crystallization, and the network topology of amorphous infills [1012.4331]. This directly addresses a common misconception: explicit anisotropy in the free energy is not required for locally anisotropic spinodal infills.

The same model produces a sequence of structural motifs. Strong length scale selectivity induces the appearance of small precursors for crystallization with icosahedral order. These precursors grow in densely packed clusters of tetrahedra; about $90\%$ of particles participate in at least one regular tetrahedron, and the local tetrahedral order parameter distribution peaks at $q_3\sim 0.73$. As the average size of the amorphous nuclei becomes large enough to reduce geometric frustration, crystalline particles with a body center cubic symmetry heterogeneously nucleate on the growing clusters. Crystalline particles begin to appear at radii $r\sim r_c \approx (3–5)\,a$, with $a$ the mean interparticle spacing, and the crystalline volume fraction obeys
\[
fc_T \sim \tau_r^{-1/2}
\]
in the studied regime [1012.4331]. The resulting infill is composite: amorphous cores, bcc shells or patches, and dense polytetrahedral aggregates.

In coherent crystalline solids, anisotropy is more explicit. For LiFePO$_4$/FePO$_4$, simultaneous expansion along $a,b$ and contraction along $c$ leads to “tilted, striped phase boundaries in equilibrium.” Three-dimensional simulations reveal skewed, striped phases whose boundaries align along $\{101\}$ or $\{100\}$ planes, in qualitative agreement with experiments by Laffont et al., Ramana et al., and Chen et al. [1202.1626]. This is not merely a geometric detail: the preferred instability direction is selected by the tensorial mismatch $\mathbf{M}$ and the elastic kernel in Fourier space.

In fluids with anisotropic diffusion, the morphology becomes directionally focused. With sufficiently strong anisotropy, only modes along the fast direction show significant growth, the spectrum becomes almost one-dimensional, and real-space domains become almost filamentary or layered along the fast axis [1810.01714]. In topological-insulator alloys, the same broad logic appears in a crystalline setting: the TlBiX$_2$ systems are modeled with composition modulation along rhombohedral $(111)$, so interfaces lie in $(111)$ planes and form lamellar or stripe-like spinodal superlattices. These interfaces act as quasi-two-dimensional channels embedded in a three-dimensional matrix [1803.06289].

## 3. Spinodoids and architected anisotropic infills

Spinodoid microstructures are synthetic, spinodal-inspired architectures generated not by direct solution of the Cahn–Hilliard equations but by compact stochastic parametrizations. In one framework, a Gaussian random field is written as
\[
\varphi(x) = \sqrt{\frac{2}{N}} \sum_{i=1}^{N\gg 1} \cos\!\big(\beta\, n_i \cdot x + \gamma_i\big),
\]
and thresholded through
\[
\xi(x) =
\begin{cases}
1 & \text{if } \varphi(x) \le \varphi_0,\\
0 & \text{if } \varphi(x) > \varphi_0.
\end{cases}
\]
The design parameters are
\[
\Theta = (\rho,\theta_1,\theta_2,\theta_3),
\]
with $\rho$ the average relative density and $\theta_1,\theta_2,\theta_3$ cone angles controlling the anisotropic orientation distribution function. Adding a rotation angle $\alpha$ produces
\[
p = (\rho,\theta_1,\theta_2,\theta_3,\alpha)^T,
\]
and the rotated homogenized stiffness is
\[
\hat{\mathbb{C}}_v(p)=T(\alpha)\,\mathbb{C}_v(f(\chi))\,T(\alpha)^T.
\]
A deep neural network $F_\omega:\Theta\mapsto S$ is trained on $n_{\text{train}}=19{,}170$ homogenized samples, achieves $R^2\ge 0.996$ on test data for each stiffness component, and reduces single predictions to approximately $0.001\ \mathrm{s}$ instead of minutes for finite-element homogenization [2012.15744].

A second framework uses Spectral Density Functions. A target phase field $\phi_T(\mathbf{X})$ with SDF $S_T(\mathbf{k})$ is combined with white noise $\phi_W(\mathbf{X})$ to reconstruct
\[
S_R(\mathbf{k}) = S_T(\mathbf{k}) \cdot S_W(\mathbf{k}),
\qquad
\phi_R(\mathbf{X}) = \mathcal{F}^{-1}\!\big( \mathcal{F}[\phi_T(\mathbf{X})] \cdot \mathcal{F}[\phi_W(\mathbf{X})] \big).
\]
Thresholding then gives
\[
A(\mathbf{X}) =
\begin{cases}
1,& \phi_R(\mathbf{X}) \le \phi_{\text{cut}}(p_m),\\
0,& \text{otherwise}.
\end{cases}
\]
Within this representation, anisotropy is explicit. Isotropic spinodoids are described by two Lamé constants. Monoclinic spinodoids have one plane of symmetry and six independent stiffness components $[C_1,\dots,C_6]$. Orthotropic spinodoids have two orthogonal symmetry planes and four independent components $[C_1,\dots,C_4]$. The anisotropy controls are
\[
a_{\text{mon}} = \sin(\theta_{\text{mon}}),
\qquad
a_{\text{ort}} = \cos(\theta_{\text{ort}}),
\qquad
a_{\text{grt}} = \cos(\theta_{\text{grt}}),
\]
together with a frequency rotational angle $\gamma$. A Gaussian Process surrogate $Q(\mathbf{s}) = m(\mathbf{s})^T\beta + \xi(\mathbf{s})$ maps descriptors to effective moduli, with prediction errors $(\mathrm{RRMSE}) \sim 10^{-3}$ for each spinodoid type [2506.23420].

These formulations also embed anisotropy into topology optimization. In the two-scale framework, the objective is compliance minimization subject to equilibrium and a volume constraint,
\[
J(\theta)=\mathbf{F}^T\mathbf{U},
\qquad
K(\theta)\mathbf{U}=\mathbf{F},
\]
while in the data-driven multiscale framework the design variables per element are
\[
\theta_e = \{ \mathbf{t}_e,\; p_M^e,\; p_m^e,\; k^e,\; \gamma^e,\; a_{\text{mon}}^e,\; a_{\text{ort}}^e,\; a_{\text{grt}}^e\}.
\]
A shallow neural network with 2 input neurons, 20 hidden neurons, and 10 output neurons maps element centroid coordinates to these descriptors, and gradients are computed through JAX automatic differentiation [2506.23420]. A recurring result is that monoclinic spinodoids are selected where stress trajectories are predominantly unidirectional, orthotropic spinodoids where two dominant orthogonal load paths coexist, and isotropic spinodoids where stress directions are highly multi-directional [2506.23420].

## 4. Mechanical response, curvature, and robustness

The mechanical literature distinguishes solid spinodal models from shell spinodal models. In numerical and experimental studies of cellular materials with spinodal topologies, solid spinodal models in the density range of $30$–$70\%$ are relatively inefficient, with
\[
\frac{E}{E_s} \sim \rho^n,\quad n \approx 2.0 - 2.6,
\qquad
\frac{\sigma_y}{\sigma_{ys}} \sim \rho^m,\quad m \approx 1.7 - 2.3.
\]
Shell spinodal models in the density range of $0.01$–$1\%$ are exceptionally stiff and strong, with
\[
\frac{E}{E_s} \sim \rho^n,\quad n \approx 1.2 - 1.3,
\qquad
\frac{\sigma_y}{\sigma_{ys}} \sim \rho^m,\quad m \approx 1.2 - 1.3,
\]
and they are remarkably imperfection insensitive [1904.06733]. This distinction is fundamental for the notion of infill: shell-based spinodal infills are not a minor variation on solid bicontinuous networks but a different mechanical regime.

Interpenetrating phase composites with spinodal shell reinforcements show the same pattern. For spinodal shell IPCs,
\[
E \propto V_f^{1.4},
\qquad
\sigma_y \propto V_f^{1.6},
\qquad
U \propto V_f^{1.1},
\]
whereas spinodal solid IPCs follow
\[
E \propto V_f^{2.2},
\qquad
\sigma_y \propto V_f^{2.7},
\qquad
U \propto V_f^{1.8}.
\]
Spinodal shell IPCs have comparable compressive strength and stiffness to IPCs reinforced by the Schwarz P TPMS and the octet truss-lattice, but exhibit far less catastrophic failure and greater damage resistance, particularly at high volume fraction of reinforcing phase [2102.06707]. The stochastic topology produces stochastically distributed cracks rather than shear-band-like localization, and the large surface area of the shell promotes matrix constraint and redistribution of load.

A more recent shell-based theory makes the curvature dependence explicit. In shell-based spinodal architected materials, each local patch is described by principal curvatures $\kappa_1,\kappa_2$, mean curvature
\[
H = \frac{\kappa_1+\kappa_2}{2},
\]
and Gaussian curvature
\[
K = \kappa_1\kappa_2.
\]
The areal stretching energy density under an affine macroscopic deformation in direction $\mathbf{e_d}$ is
\[
\overline{W}_s = U_g^2 \frac{E h}{2(1-\nu^2)}\left(1 - (\mathbf{e_d} \cdot \hat{\mathbf{x}}_3)^2\right)^2,
\]
and the geometric proxy for stretch-to-bend ratio is
\[
\Gamma_p(\mathbf{e_d}) =
\frac{\sum_{k=1}^N \sin^2\alpha_k S_k}{\sum_{k=1}^N \beta_{p,k} S_k}.
\]
A corresponding strength proxy is
\[
\Gamma_{s,p}(\mathbf{e_d}) =
\frac{\left(\sum_{k=1}^N \sin^2\alpha_k S_k\right)^{3/2}}{\sum_{k=1}^N \beta_{p,k} S_k}.
\]
The paper shows that directional effective stiffness scales approximately as $E^*(\mathbf{e_d}) \propto \Gamma^2(\mathbf{e_d})$, and that normalized stiffness and yield strength correlate strongly with the geometric proxies [2505.21509]. This provides a direct route from curvature distributions to anisotropic infill design.

## 5. Functional material realizations

One realization of anisotropic spinodal infills is electronic rather than structural. In TlBiX$_2$ ($X=$ S, Se, Te), all three binary alloys considered have positive formation enthalpy over the entire composition range, are immiscible at $0\ \mathrm{K}$, and display miscibility gaps at finite temperature. Reported consolute temperatures and critical compositions are approximately $162\ \mathrm{K}$ and $x_c\approx0.4$ for TlBiS$_2$–TlBiSe$_2$, $400\ \mathrm{K}$ and $x_c\approx0.4$ for TlBiSe$_2$–TlBiTe$_2$, and $1040\ \mathrm{K}$ and $x_c\approx0.35$ for TlBiS$_2$–TlBiTe$_2$ [1803.06289]. In the TlBiS$_2$–TlBiTe$_2$ system, the eutectic regime around $x_e\sim64\%$ TlBiTe$_2$ at $\sim810\ \mathrm{K}$ allows spinodal superlattices and two-dimensional eutectic microstructures. The resulting interfaces host gapless or quasi-gapless states, and the work explicitly frames them as “quasi-2D metallic channels inside a 3D insulating or weakly conducting matrix” without vacuum interfaces [1803.06289].

A second realization is electro-chemo-mechanical. In LiFePO$_4$/FePO$_4$, anisotropic coherency strain narrows the spinodal gap, selects preferred wavevectors, and yields tilted, striped phase boundaries in equilibrium [1202.1626]. The material-specific eigenstrain values
\[
M_{11}=0.052,\quad M_{22}=0.036,\quad M_{33}=-0.019
\]
produce direction-dependent suppression of instability and favor fluctuations along the direction of greatest expansion. The resulting striped or lamellar infills are therefore elasticity selected rather than merely diffusion limited.

A third realization is magnetic. In alnico, spinodal decomposition of a bcc solid solution into Fe–Co–rich $\alpha_1$ and Ni–Al–rich $\alpha_2$ is biased by magnetic annealing, which aligns $\alpha_1$ rods along the $\langle100\rangle$ direction closest to the applied field. The microstructure evolves from an initially isotropic, interpenetrating network to highly anisotropic rods, followed by Cu-rich spheres, rods, and blades and additional $\alpha_1$ rod precipitation during draw. The intrinsic coercivity rises from $35\ \mathrm{Oe}$ in the solutionized condition to $476\ \mathrm{Oe}$ after $30\ \mathrm{s}$ magnetic annealing, to about $600\ \mathrm{Oe}$ after $10\ \mathrm{min}$ magnetic annealing, and to $1707\ \mathrm{Oe}$ after full draw [1810.12580]. In this case, “anisotropic spinodal infills” are the ensemble of rods, blades, and interfacial phases that control shape anisotropy and magnetic isolation.

## 6. Misconceptions, limitations, and open directions

A first misconception is that spinodal infills are necessarily periodic. The architected-material literature explicitly defines spinodoids as stochastic and aperiodic, and the mechanically efficient shell topologies generated from Cahn–Hilliard evolution are likewise non-periodic [2012.15744, 1904.06733]. A second misconception is that anisotropy requires a built-in anisotropic free energy. In the liquid–solid Landau model, the gradient and nonlocal terms are isotropic, yet local anisotropy emerges through random-wave morphology, curvature-dependent heterogeneous nucleation, and the topology of the amorphous network [1012.4331].

The present theories also have clear regime limits. The curvature-guided shell framework is developed for the thin-shell regime
\[
0.01 \lesssim \kappa h \lesssim 0.1,
\]
and outside this range transverse shear and more complex three-dimensional stress states become important [2505.21509]. The two-scale topology-optimization framework assumes separation of scales and is restricted to linear elastic stiffness; nonlinear behavior, yielding, large strains, and damage are not part of the surrogate [2012.15744]. The topological-insulator work does not solve the full time-dependent Cahn–Hilliard equation but assumes sinusoidal composition modulation along $(111)$ and periodic superlattice structures [1803.06289].

Several open directions are already identified inside the literature. The isotropic liquid–solid Landau model states that anisotropic gradient coefficients $D_{ij}$ or anisotropic nonlocal kernels would produce direction-dependent correlation lengths and genuinely anisotropic spinodal morphological length scales [1012.4331]. The architected-material literature points to nonlinear and multi-physics surrogate models, full 3D rotations, adaptive resolution, and manufacturing-aware optimization as natural extensions [2012.15744]. In filled polymer blends, the current Sanchez–Lacombe framework treats fillers as isotropic spherical inclusions, but the paper explicitly argues that an anisotropic generalization would require orientation distribution functions, tensorial stress responses, and a direction-dependent instability analysis; its approximate spinodal formula already gives deviations in the spinodal temperature of less than $4\ \mathrm{K}$ relative to the exact calculation in the low-compressibility regime [2601.10433].

Taken together, these works suggest that “anisotropic spinodal infills” is not a single microstructure class but a family of spinodal-derived internal organizations. They range from amorphous polytetrahedral skeletons with bcc patches, to tilted striped phases selected by coherency strain, to lamellar topological channels, to shell-based spinodoid architectures whose anisotropy is optimized against local stress trajectories. The unifying principle is that a spinodal instability supplies a bicontinuous or finite-wavelength scaffold, while anisotropy emerges from elastic mismatch, diffusion tensors, magnetic or crystallographic bias, surface-energy design, or data-driven multiscale optimization [1012.4331, 1202.1626, 1803.06289, 2506.23420].

Source: https://www.emergentmind.com/topics/anisotropic-spinodal-infills