---
title: Liquid Metal Dealloying (LMD)
url: https://www.emergentmind.com/topics/liquid-metal-dealloying-lmd
type: topic
---

# Liquid Metal Dealloying (LMD)

Liquid metal dealloying (LMD) is the corrosion of a two-component solid alloy by contact with a molten third component, leading to the selective dissolution of one species into the liquid and the self-organization of the retained species into a nano-porous or bicontinuous network with ultra-high interfacial area. In current theoretical and computational treatments, LMD is a diffusion-coupled moving-boundary problem in which local interfacial equilibrium, finite solubility of the retained element in the melt, interfacial instability, grain-boundary transport, and solvent chemistry jointly determine front kinetics, morphology, topology, and the resulting mechanical response [2202.01984, 2202.04025, 2206.01234, 2509.20770, 2604.22726].

## 1. Thermodynamic basis and governing fields

A standard LMD setting consists of a precursor solid alloy \(A_xB_{1-x}\), where \(A\) is the immiscible element to be retained and \(B\) is the dissolvable element, immersed in a liquid metal melt \(M\). At the solid–liquid interface, local equilibrium imposes chemical potentials \(\mu_i^s(\phi,c)=\mu_i^l(\phi,c)\). The small but finite solubility of \(A\) in \(M\), measured from ternary phase diagrams, enables leakage of the nominally retained species into the melt; values such as \(c_{\mathrm{Ta}}^l\approx 10^{-3}\!-\!10^{-2}\) are explicitly invoked in LMD analyses [2202.04025].

Phase-field descriptions introduce a non-conserved field \(\phi(x,t)\), with \(\phi\approx 1\) in the solid and \(\phi\approx 0\) in the liquid, together with conserved composition fields \(c_A(x,t)\) and \(c_B(x,t)\), with \(c_C=1-c_A-c_B\). One representative free-energy functional is
\[
F[\phi,c_A,c_B]=\int_\Omega \Bigl[\frac\eta2|\nabla\phi|^2 + f_\mathrm{int}(\phi)
+ g(\phi)\,f_\mathrm{chem}(c_A,c_B)\Bigr]\,dV,
\]
with \(f_\mathrm{int}(\phi)=W\phi^2(1-\phi)^2\) and
\[
f_\mathrm{chem}(c_A,c_B)=\sum_i c_i\ln c_i + \Omega_{AB}c_A c_B + \Omega_{AC}c_A(1-c_A-c_B)+\ldots
\]
The corresponding dynamics combine a modified Allen–Cahn evolution for \(\phi\) and coupled Cahn–Hilliard transport for the conserved species,
\[
\frac{\partial \phi}{\partial t}
= -\tilde M_\phi\,\frac{\pi^2}{8\eta}\,\frac{\delta F}{\delta\phi},
\qquad
\frac{\partial c_i}{\partial t}
= \nabla\!\cdot\!\sum_{j=\{A,B\}} M_{ij}(\phi)\,\nabla\!\Bigl(\frac{\delta F}{\delta c_j}\Bigr),
\]
or, in shorthand,
\[
\frac{\partial c}{\partial t}=\nabla\!\cdot\!\bigl[M(c)\,\nabla\mu\bigr].
\]
These formulations are central in both mechanistic studies and surrogate-model development for LMD [2509.20770].

Alternative phase-field formulations enrich the chemical and interfacial terms. In the treatment summarized by Lai–Geslin–Karma, the free energy contains regular-solution contributions, latent-heat terms, and composition-gradient penalties, while the mobility matrix takes the form \(M_{ij}(\phi)=D_{\mathrm{liq}}(1-\phi)(V_a/kT)c_i(\delta_{ij}-c_j)\). Across these formulations, the same structural point recurs: selective dissolution is thermodynamically driven by the difference in chemical free energies of the species in the solid and liquid, while the retained network is stabilized and reshaped by gradient penalties and interfacial energy [2202.01984, 2406.17119].

A persistent computational consequence is stiffness. Because the conserved-species equations contain fourth-order spatial derivatives through terms such as \(\nabla\cdot[M\nabla(\delta F/\delta c)]\), explicit solvers require \(\Delta t\approx 10^{-12}\,\mathrm{s}\) for stability. This stiffness is one of the defining features of LMD phase-field computation and motivates reduced-order or learned surrogates [2509.20770].

## 2. Front propagation and instability selection

The kinetics of LMD have been analyzed through one-dimensional dissolution models and higher-dimensional stability analyses. In diffusion-limited form, the dealloying front position obeys
\[
x_i(t)\simeq \sqrt{4 p D_l t},
\]
where \(D_l\) is the diffusivity of the dissolving species in the melt and \(p\) is a Péclet number depending on interface reaction rates and diffusivities; at long times, the front velocity scales as \(V\equiv dx_i/dt\sim \sqrt{D_l/t}\). This \(\sqrt{t}\) law is a central scaling for LMD [2202.04025].

A more detailed one-dimensional treatment distinguishes two regimes. In the early stage, the low solubility of the retained element in the melt causes accumulation of that element at the interface, generating a passivation layer. The peak retained-species fraction satisfies \(c_A^p(t)\simeq (c_A^0)x_i(t)/\xi\), and a simple phenomenology for the local front velocity is \(v_i=v_0\exp[-c_A^p/c^\*]\), with \(c^\*\sim 0.04\!-\!0.05\). This yields
\[
v_i(t)=\frac{v_0}{1+t/\tau},
\qquad
x_i(t)=v_0\tau\ln(1+t/\tau).
\]
Once the interfacial composition reaches its tiny equilibrium solubility in the melt, the system enters a stationary regime in which both elements dissolve at fixed interfacial compositions and the classical diffusion-controlled \(\sqrt{t}\) law emerges [2202.01984].

Pattern formation originates in instability of the dealloying front. Linear stability analysis of the interfacial alloy layer gives a Cahn–Hilliard-type dispersion relation,
\[
\omega(k)=-Mk^2\,[f''+\Gamma k^2],
\]
and, in the ternary-reduced binary formulation,
\[
\omega(k)=\frac12 k^2 \Bigl[-(A+D)+\sqrt{(A+D)^2-4(AD-BC)}\Bigr].
\]
Spinodal instability occurs if
\[
f_s \equiv \det M \cdot (f_{12}^2-f_{11}f_{22})>0,
\]
and the fastest-growing wavenumber \(k_{\max}\) determines the initial ligament spacing \(\lambda_0=2\pi/k_{\max}\), reported as \((1.5\!-\!3)w\) in phase-field units. The instability criterion also predicts a transition between planar dissolution and spinodal-driven ligamentation as the precursor and melt compositions vary [2202.01984].

An important correction to earlier intuition is the role of solid-state diffusion. Although the solid diffusivity is four to five orders of magnitude smaller than the liquid diffusivity, it affects both dissolution kinetics and ligament morphologies. Finite \(D_s\) broadens the retained-species interfacial peak, increases the effective Péclet number toward the phase-diagram prediction, and can promote parallel laminæ rather than highly tortuous bicontinuous networks. This directly contradicts the common simplification that solid-state transport may be neglected without qualitative consequence [2202.01984].

## 3. Connectivity, topology, and solvent chemistry

Topology control in LMD is commonly discussed in terms of the solid volume fraction \(\rho\) in the dealloyed layer. In mass-conservation form, \(\rho=c_A^0/c_A^s\) in the absence of leakage of the retained element. With retained-element leakage at a rate determined by \(c_A^l\) and \(p\),
\[
\rho=\frac{2p c_A^0-c_A^l}{2p c_A^s-c_A^l}.
\]
Genus number \(g\) correlates with \(\rho\) by topological invariants, and highly bicontinuous structures require \(\rho\) above a percolation threshold \(\rho_c\approx 0.3\). In this framework, bulk diffusive transport of the immiscible element in the liquid melt strongly influences the evolution of the solid fraction and topology during dealloying [2202.04025].

A compact topology criterion is the ratio
\[
\Lambda \equiv \frac{c_A^l}{2p c_A^0}.
\]
For \(\Lambda\gtrsim 1\), corresponding to high leak or low front velocity, \(\rho\) at the front falls below threshold and connectivity is lost. For \(\Lambda\lesssim 1\), connectivity is maintained. This parameterization clarifies why melt chemistry exerts such strong morphological control. Adding Ti to Cu increases \(B\) solubility in the melt, lowers front velocity, but increases \(A\) leak, which lowers \(\rho(t)\) below the percolation threshold and yields structures that detach, fragment, or dissolve. Adding Ag to Cu lowers both \(c_A^l\) and \(c_B^l\), reduces \(A\) leakage, keeps \(\rho\) above \(\rho_c\), and preserves high-genus bicontinuous topology [2202.04025].

LMD is often conflated with electrochemical dealloying (ECD), but the two processes differ in their controlling transport and topology-selection mechanisms.

| Aspect | ECD | LMD |
|---|---|---|
| Front rate | Constant interface-controlled \(V\) | Diffusion-limited \(V(t)\sim \sqrt{D_l/t}\) |
| Retained species in fluid | Negligible bulk diffusion, \(c_A^l\approx 0\) | Finite solubility and bulk diffusion of \(A\) |
| Pattern formation | Surface-spinodal decomposition | Unstable diffusion-coupled growth |
| Solid fraction \(\rho\) | Remains uniform | Spatially varying due to leakage |

This distinction matters for both interpretation and design. In ECD, topologically connected bicontinuous structures occur only at low \(B\) content and low voltage. In LMD, the front velocity rapidly drops because of diffusion-limited kinetics, and bulk diffusion of the retained species generates spatially varying \(\rho\) and topology [2202.04025].

Experimental comparisons support the topology criterion. For Ta–Ti precursors dealloyed for \(10\,\mathrm{s}\) at \(1240\,^\circ\mathrm{C}\), pure Cu produced rapid dealloying depth of \(\sim 270\,\mu\mathrm{m}\), an \(A\)-leak layer \(h_L\sim 20\,\mu\mathrm{m}\), a dealloyed layer \(h_D\sim 250\,\mu\mathrm{m}\), and a connected core \(h_C\sim 42\,\mu\mathrm{m}\), with a predominantly disconnected, fine porous layer near the interface. Under otherwise comparable conditions, \(\mathrm{Cu}_{70}\mathrm{Ag}_{30}\) produced a slower depth of \(\sim 220\,\mu\mathrm{m}\), no leak layer, \(h_D\sim h_C\sim 220\,\mu\mathrm{m}\), and a fully connected bicontinuous morphology with uniform \(\sim 200\,\mathrm{nm}\) ligament size. The same study also reports that the classical coarsening law \(\lambda^n=k t_c+\lambda_0^n\) with \(n=4\) fails to capture LMD data; instead, bulk diffusive transport of the retained species and concurrent dissolution/redeposition dominate ligament evolution [2202.04025].

## 4. Grain-boundary-mediated dealloying

Grain boundaries are not passive defects in LMD. A multi-phase-field model extending the original approach of Geslin et al. treats two solid grains and one liquid using three non-conserved phase fields \(\phi_1,\phi_2,\phi_3\) and three conserved concentrations \(c_A,c_B,c_C\). In that formulation, the local diffusivity is
\[
D(\{\phi\})=D_\ell\,\phi_3 +16\,D_{gb}\,\phi_1^2\phi_2^2,
\]
so diffusion is fast in the liquid, negligible in the bulk solid, and enhanced along the diffuse grain boundary. The model is implemented with explicit Euler for the phase fields and a semi-implicit spectral scheme, following Badalassi et al., for the concentrations, enabling large \(2\mathrm{D}\) and \(3\mathrm{D}\) simulations while maintaining a realistic interfacial thickness [2206.01234].

The governing picture is one of coupled grain-boundary migration and dealloying. Grain-boundary migration velocity obeys \(v_{gb}=M_{gb}f\), where the driving force \(f\) is the local excess grand-potential gradient. As the liquid–solid triple line at the grain boundary dissolves \(B\) more rapidly because of enhanced boundary diffusion, the local solid becomes enriched in \(A\). Uphill diffusion of \(A\) along the migrating boundary provides an extra evacuation path for the more noble element rejected by the advancing liquid channel. Two reinforcing mechanisms are reported: a solute-flux mechanism and a curvature mechanism. The triple junction maintains the Young’s-law force balance to within numerical tolerance [2206.01234].

The morphological consequences are pronounced. In \(2\mathrm{D}\), for initial \(c_{A,0}=40\%\), bicontinuous dealloying forms a uniform ligament network within each grain, but at the grain boundary a single deep channel penetrates \(\gtrsim 2\times\) further. In \(3\mathrm{D}\), for \(c_{A,0}=20\%\), regions away from the boundary develop \(3\mathrm{D}\) nodules that grow into ligaments, whereas the migrating boundary develops \(2\mathrm{D}\) ridges and linear corrosion channels. The boundary migrates in opposite \(Z\)-directions at different \(Y\)-locations, a result presented as a possible origin of experimentally observed waviness [2206.01234].

Grain-boundary migration also disrupts final connectivity. The \(A\)-rich block formed by the migrating boundary, termed the GBMD region, attaches selectively to one grain’s ligament network and leaves a widening liquid gap to the other grain, asymmetrically disrupting connectivity. This behavior is reported to be in qualitative agreement with SEM observations of Ta–Ti LMD by McCue et al. and Fe–Ni LMD by Joo et al. A further result is that grain-boundary migration can locally bypass the parting-limit constraint: even for \(c_{A,0}=50\%\), where no bulk bicontinuous front forms in single-crystal runs, the boundary still nucleates a singular deep liquid channel via GBMD [2206.01234].

## 5. Phase-field computation and surrogate extrapolation

Direct numerical simulation of LMD phase-field models is computationally expensive because the stiffness of the conserved-species equations forces extremely small time steps and because the evolving microstructure must be resolved on fine grids. Representative DNS setups use \(O(N_x\times N_y)\) degrees of freedom such as \(512\times 512\) to resolve sub-nanometer features, with \(\Delta t\approx 10^{-12}\,\mathrm{s}\). Reaching \(\sim 1\,\mu\mathrm{s}\) then requires \(\gtrsim 10^6\) time steps, and a single \(512\times 512\), \(1\,\mu\mathrm{s}\) run is reported to take \(\sim 7.5\) days on a large CPU cluster; extending to \(1024\times 512\) or to times \(\sim 3\times\) longer becomes intractable, with cost scaling roughly as \(O(N^2\cdot T/\Delta t)\) [2509.20770].

One response has been operator-learning surrogates. The U-Shaped Adaptive Fourier Neural Operator (U-AFNO) couples a U-Net encoder–decoder to an Adaptive Fourier Neural Operator bottleneck implemented as a frequency-space vision transformer. In the reported configuration, the latent AFNO contains \(L=12\) blocks and \(H=16\) heads, and the model is trained on 87 high-fidelity runs with snapshots every \(0.05\,\mu\mathrm{s}\). After 100 auto-regressive leaps, U-AFNO achieves \(e_{AC}(\phi)\lesssim 15\%\), \(e_{AC}(c_A)\lesssim 15\%\), and \(e_{AC}(c_B)\lesssim 20\%\), on par with the intrinsic high-fidelity-to-high-fidelity discrepancy across different initial seeds, while the relative \(\ell^2\) error in total solid mass is \(\sim 1\%\). The solver cost is reported as \(1500\,\mathrm{s}\) per \(5\times 10^4\) steps on 128 CPU cores, versus \(0.116\,\mathrm{s}\) per leap on one A100 GPU [2406.17119].

A later direction uses purely convolutional, conditionally parameterized U-Nets. In one version, the network has 5 encoder levels with \(3\times 3\) convolutions and stride 2, 5 decoder levels with transposed \(3\times 3\) convolutions, skip connections, convolutional self-attention in the bottleneck, and physics-aware padding: circular in the horizontal direction, zero at the top boundary, and value replication at the bottom. Parameter conditioning uses an MLP to inject the time-skipping interval \(\Delta\tau\) and a concentration parameter, enabling leaps of \(50\mathrm{k}\!-\!100\mathrm{k}\) solver steps and adaptation to diverse alloy systems. During inference, translational invariance allows the same weights to process larger domains such as \(1024\times 512\), and auto-regressive rollout applies the learned map repeatedly in time. Within the training regime, quantities-of-interest errors are reported as \(<5\%\); when extrapolating to domains twice as tall and time horizons \(\gtrsim 3\times\) longer, the error remains \(<10\%\). For a \(512\times 512 \rightarrow 1\,\mu\mathrm{s}\) reference simulation, the surrogate completes the same horizon in \(38.5\,\mathrm{s}\) on one A100 GPU, corresponding to a \(\sim 16{,}000\times\) speed-up [2509.20770].

A related 2026 framework extends the convolution-only strategy with a flood-fill field corrector and a conditional diffusion model for initial conditions. The U-Net conditions on \(\theta=(\Delta\tau,c_A^\mathrm{init})\), enforces fields in \([0,1]\) through a final sigmoid, and uses physically informed padding to match periodic side boundaries, a top Dirichlet condition, and a bottom Neumann condition. The flood-fill corrector resets predicted uncorroded regions to the exact initial concentration, while a DDPM generates synthetic deformed-interface initial conditions conditioned on \(c_A\). On the reported benchmarks, QoI relative errors are \(\lesssim 5\%\) in the training regime and \(\lesssim 15\%\) in large-scale, long-time extrapolation, with worst case \(\sim 20\%\) for high \(c_A\). One U-Net forward pass with \(\Delta\tau=50{,}000\) takes \(\approx 0.15\,\mathrm{s}\) on A100 versus 59 minutes for the direct solver, and end-to-end speed-ups reach \(\approx 35{,}900\times\). Synthetic initial conditions degrade QoI errors by only \(1\!-\!2\%\) in training and negligibly in extrapolation [2601.04510].

These surrogate studies also delimit current failure modes. Because LMD dynamics are chaotic, the surrogates do not reproduce pixel-wise microstructures exactly; instead they target quantities of interest such as interface curvature, penetration depth, ligament statistics, and spatial auto-correlations. The ligament height metric is repeatedly identified as the most challenging. Extreme auto-regressive extrapolation may lead to unphysical patterns or drift, and the models do not exactly satisfy PDE residuals away from the explicitly corrected regions [2509.20770, 2601.04510].

## 6. Mechanical scaling and materials design

LMD is used not only to generate morphology but also to tune mechanical response through topology and connectivity. Ramallo et al. study single-crystalline nanoporous tantalum produced by liquid metal dealloying of Ti65Ta35 in a Cu40Bi60 melt at \(\sim 800\,^\circ\mathrm{C}\) for 5 minutes under Ar. The resulting network has ligament diameter \(l\approx 200\pm 100\,\mathrm{nm}\), solid volume fraction \(\Phi\approx 0.35\), pure bcc Ta by XRD, and \(<2\) at.% oxygen by SEM/EDS [2604.22726].

The mechanical behavior is expressed relative to Gibson–Ashby scaling for open-cell foams,
\[
\frac{E}{E_{\rm lig}}=C_E\,\Phi^{\,n_E},
\qquad
\frac{H}{H_{\rm lig}}=C_H\,\Phi^{\,n_H}.
\]
For nanoporous Ta with \(\Phi=0.24\!-\!0.42\), nanoindentation yields \(E=10\!-\!30\,\mathrm{GPa}\) and \(H=0.3\!-\!1.1\,\mathrm{GPa}\), with the data best fit by
\[
E \approx E_{\rm bulk}\,\Phi^2,\qquad
H \approx H_{\rm bulk}\,\Phi^2.
\]
The paper explicitly contrasts this with nanoporous Au, where \(n_H\approx 1.5\), and attributes the difference to enhanced ligament connectivity in the Cu–Bi LMD system [2604.22726].

The mechanistic basis of that enhanced connectivity is solvent chemistry. In a pure Cu bath, Ta has finite solubility at the solid/fluid interface, leading to broken ligaments and lower connectivity. Adding Bi lowers the activity of Cu by forming a Cu–Bi eutectic with \(T_{\mathrm{eut}}\approx 271\,^\circ\mathrm{C}\), so that
\[
c_{\rm Ta}^{\ell}(\mathrm{in}\ \mathrm{Cu\!-\!Bi}) \ll c_{\rm Ta}^{\ell}(\mathrm{in}\ \mathrm{Cu}).
\]
The reported interpretation is that reduced Ta leak results in a more bicontinuous network, with fewer dead-ends, higher genus, and more uniform load paths. This directly parallels the Ag-in-Cu topology-control strategy reported for Ta–Ti dealloying, where reduced leak of the retained species preserves connectivity and suppresses fragmentation [2202.04025, 2604.22726].

Molecular-dynamics simulations indicate that the deformation mechanism is not anomalous. Using LAMMPS with an extended Finnis–Sinclair potential and an atomistic diamond Berkovich tip, the simulations show surface-nucleated glissile \( \tfrac12\langle 111\rangle \) dislocations, annihilation at free surfaces, occasional \(a\langle 100\rangle\) sessile lines, and localized twinning below the contact. Densification under the indenter is minimal, \(\Delta\Phi\lesssim 2\%\), and the plastic zone extends only \(\sim 3\!-\!4\times\) the penetration depth, in agreement with post-indent SEM. The stiffness-density response is therefore attributed to connectivity rather than to unusual deformation mechanisms [2604.22726].

The resulting design rules are consistent across the literature. Melt compositions that minimize solubility of the retained element, such as Ag additions to Cu or Bi-containing Cu baths, preserve bicontinuity. Temperature and time control coarsening and volume fraction: higher \(T\) or longer \(t\) coarsen ligaments and reduce \(\Phi\), whereas lower \(T\) or shorter \(t\) yield finer ligaments and higher \(\Phi\), although incomplete dissolution becomes a competing risk. Precursor composition tunes the amount of dissolvable species and thus the final solid fraction. Taken together, these results establish solvent chemistry, front kinetics, and retained-species leakage as coupled levers for topology and mechanical scaling in LMD-derived nanoporous metals [2202.04025, 2604.22726].

Source: https://www.emergentmind.com/topics/liquid-metal-dealloying-lmd