---
title: Differentiable Lithography Techniques
url: https://www.emergentmind.com/topics/differentiable-lithography
type: topic
---

# Differentiable Lithography Techniques

Differentiable lithography is a family of methods that represents lithographic image formation, resist response, fabrication transfer, or related manufacturing distortions as differentiable operators inside an optimization loop, so gradients can be propagated to masks, sources, optical parameters, device geometries, or fabrication layouts. In semiconductor computational lithography, this usually takes the form of differentiable Abbe or Hopkins imaging, differentiable resist surrogates, and gradient-based OPC or ILT. In computational optics, it includes learned fabrication “digital twins” that map intended diffractive layouts to predicted as-fabricated \(2.5\)D or \(3\)D structures, allowing optical objectives to be evaluated on fabricated rather than ideal geometry. In photonics and EUV mask design, it extends to differentiable fabrication twins, physics-grounded inverse lithography, and rigorous electromagnetic solvers embedded in automatic differentiation loops.

## 1. Scope and taxonomy

Differentiable lithography is not a single algorithmic template. The literature spans at least four distinct but connected formulations: differentiable lithography imaging for semiconductor manufacturing, differentiable mask-geometry optimization, fabrication-aware differentiable optics for manufactured diffractive elements, and domain-specific extensions in photonics, EUV masks, and computed axial lithography. A useful distinction is whether the differentiable object is the imaging pipeline, the mask parameterization, the fabrication process, or the full electromagnetic forward model [2409.15306] [2505.22313] [2602.15762].

| Strand | Core differentiable object | Representative papers |
|---|---|---|
| Semiconductor imaging and resist | Illumination–projector–mask–resist graph | [2409.15306], [2502.06838] |
| Mask optimization | Pixels, edges, level sets, B-splines | [2408.08969], [2504.11962], [2308.12299], [2405.03574] |
| Fabrication-aware computational optics | Layout-to-fabricated-topography digital twins | [2309.17343], [2505.22313] |
| Photonics, EUV, CAL | Fabrication twins, rigorous EM solvers, blur-calibrated tomography | [2602.15762], [2606.25753], [2606.04174] |

A second distinction concerns scope. Some works differentiate through the full physical or learned process used in optimization; others are adjacent rather than core. The metalens study for zone-plate-array lithography is explicitly an instance of inverse electromagnetic design for lithography hardware, but it does not differentiate through resist chemistry, exposure thresholding, pattern fidelity, or process models [2212.06441]. Likewise, LithoDreamer is a learned world model with differentiable latent dynamics for the “Layout-Mask-Resist Image-After Development Image” pipeline, but it does not embed explicit optical imaging equations or resist PDEs in the classical sense [2606.26713].

## 2. Differentiable forward models and gradient propagation

A canonical semiconductor formulation decomposes lithography into parameterized differentiable modules for illumination \(f_i\), projector transfer \(f_h\), mask \(f_m\), and resist \(f_r\), and writes the overall forward model as
\[
f = \eta \circ f_i \circ f_h \circ f_m \circ f_r \Rightarrow Z = f(\cdot, \mathbf{\theta}) + \gamma,
\]
with \(Z\) the final resist image and \(\eta,\gamma\) denoting signal-dependent and signal-independent noise. Within this formulation, Abbe imaging, Hopkins imaging, and SOCS/TCC-SVD approximations can all be placed inside automatic differentiation. The same framework states source optimization as
\[
(\theta_{i\_opt}, \theta_{h\_opt}) := \underset{\mathbf{\theta}_i, \mathbf{\theta}_h}{\operatorname{argmin}}|f(\cdot;\mathbf{\theta}_i, \mathbf{\theta}_h) - Z_t|^2,
\]
mask optimization as
\[
{M}_{opt} = \underset{\mathbf{\theta}_m}{\operatorname{argmin}}|f(\cdot, \mathbf{\theta}_m) - Z_t|^2,
\]
and SMO as
\[
(M_{opt}, J_{opt}) = \underset{\mathbf{\theta}_i, \mathbf{\theta}_m}{\operatorname{argmin}}|f(\cdot;\mathbf{\theta}_i, \mathbf{\theta}_m) - Z_t|^2.
\]
In this setting, Hopkins/SOCS reduces complexity from \(\mathcal{O}(n^6)\) to \(\mathcal{O}(Q \times n^4)\), while Abbe remains attractive when gradients with respect to source variables are required [2409.15306].

A major unresolved component in many imaging-only formulations is the resist stage. TorchResist addresses this by making a compact analytical resist simulator differentiable. Its exposure model uses
\[
\frac{\partial I(h,t)}{\partial h} = - I(h,t)\,[A\,m(h,t)+B], \qquad
\frac{\partial m(h,t)}{\partial t} = - m(h,t)\, I(h,t)\, C,
\]
and its development model uses the Mack-style rate
\[
r = r_{\max} \frac{(a+1)(1-m)^n}{a+(1-m)^n} + r_{\min}.
\]
To calibrate the model against binary wafer labels, it replaces non-differentiable thresholding and \(L_0\) mismatch by the surrogate
\[
\mathrm{loss} = \mathrm{BCE}\bigl(\sigma(s f_\theta(X)-s\tau),\,Y\bigr).
\]
On LithoBench MetalSet, TorchResist reports **0.22** pixel difference, **0.73 nm** EPE-mean, and **2.87 nm** EPE-max, versus **0.49**, **1.21 nm**, and **3.95 nm** for the variable-threshold baseline; it also provides depth simulation and runs in **0.04 s** per \(2\,\mu m \times 2\,\mu m\) patch at **7 nm/pixel** on a 3090 GPU [2502.06838].

This family of formulations establishes the basic differentiable lithography principle: the objective is evaluated on the final printed or resist image, while gradients are propagated through FFTs, convolutions, source summations, threshold surrogates, and optionally resist kinetics, instead of stopping at the aerial image. A plausible implication is that wafer-level objectives become substantially more trustworthy when the resist stage is no longer reduced to fixed thresholding.

## 3. Mask parameterizations and inverse optimization

Early differentiable lithography work often optimized raster masks directly, but several later methods shift the optimization variable from pixels to geometry. DiffOPC is an explicit example of differentiable edge-based OPC. It represents a mask by boundary segments
\[
\mathbf{S} = \{\mathbf{s}_1, \mathbf{s}_2, \ldots, \mathbf{s}_{N_s}\}, \qquad
\mathbf{s}_i = \bigl[[x_1, y_1], [x_2, y_2]\bigr],
\]
constrains each segment to move along a legal normal-direction velocity \(\mathbf{v}_i\) satisfying \(\mathbf{v}_i \cdot \mathbf{d}_i = 0\), uses straight-through rounding for integer-grid coordinates, and backpropagates lithography gradients to segment motion via midpoint sampling. Its total objective combines nominal fidelity, process variation band, and a differentiable EPE surrogate,
\[
\mathcal{L}_{total} = w_1 \mathcal{L}_2 + w_2 \mathcal{L}_{pvb} + w_3 \mathcal{L}_{epe},
\]
with \(w_1=1\), \(w_2=0.9\), and \(w_3=100\). On ICCAD13 metal layers, DiffOPC reports \(L2=28280\), \(PVB=48217\), \(EPE=2.2\), and \(\#shots=106.1\), compared with MultiILT at \(28704.6\), \(42132.0\), \(2.7\), and \(276.1\); it reports zero MRC violations on ICCAD13, the larger dataset, and the via set [2408.08969].

A curvilinear alternative is the B-spline/Delaunay formulation, which makes the chain
\[
\text{control points} \to \text{sampled boundary points} \to \text{triangulation vertices} \to \text{Gaussian quadrature points} \to \text{image intensity} \to \text{loss}
\]
explicitly differentiable. The boundary is represented by periodic B-splines,
\[
\mathring{\mathbf P}(\xi) = \sum_{k=1}^{n+p}\mathbf P_k \,\mathring N_{k,p}(\xi),
\]
the enclosed domain is approximated by Delaunay triangles, and the coherent image amplitude is evaluated by Gaussian quadrature over those triangles. The loss is
\[
J := \left\| \operatorname{sig}(I)-I_{\text{target}}\right\|_2^2,
\]
with \(\operatorname{sig}(x)=\frac{1}{1+\exp(-a(x-tr))}\). The distinguishing contribution is the derivation of explicit formulas for \(\partial J/\partial P_{kx}\) and \(\partial J/\partial P_{ky}\), so control-point gradients can be computed analytically rather than by finite differences [2504.11962].

Level-set parameterizations form a third line. ILDLS represents the mask boundary by a signed-distance level set, uses a Hopkins model with \(K=24\) coherent kernels,
\[
I(x,y; h_\mu)=\sum_{k=1}^{K} w_k(h_\mu)\left|M(x,y)\otimes h_k(x,y; h_\mu)\right|^2,
\]
replaces hard thresholding by
\[
Z(x, y; h_\mu, t_q) = \frac{1}{1+\exp\left\{-\theta_z\left[\frac{I(x,y;h_\mu)}{1+t_q}-I_{th}\right]\right\}},
\]
and injects the level-set ILT correction gradient into UNet training. The final training objective is
\[
Loss = L_{fit} + \alpha L_{LS},
\]
with \(\alpha=0.008\). In test results, ILDLS reports **5.89 nm** AEDE versus **7.79 nm** for pure UNet, while ILDLS+ILT reports **4.02 nm** versus **5.67 nm** for ILT; for process window, ILDLS reports **133.9** DOF@5%EL versus **116.0** for ILT-PV [2308.12299].

A different position in the design space is ILILT, which frames inverse lithography as a learned recurrent optimizer with lithography-conditioned recurrence,
\[
\mathbf{M}_{t+1}=g(\mathbf{M}_t,\mathbf{Z}_t,\mathbf{Z}^\ast,\mathbf{w}), \qquad \mathbf{Z}_t=f_l(\mathbf{M}_t).
\]
It is relevant because the lithography simulator is queried at every step, but it is not primarily a differentiable forward-lithography paper. The method is trained mainly by supervision from solver-produced masks rather than by direct wafer-loss minimization, although it reports **EPE = 0.08** and **PVB = 4695** for ILILT-P2PHD-8, compared with **0.21** and **4656** for GPU-ILT [2405.03574].

## 4. Fabrication-aware differentiable optics and neural lithography

A major expansion of differentiable lithography occurred in computational optics, where the central problem is not mask-to-wafer imaging but layout-to-fabricated-device mismatch. “Neural Lithography” introduces a learned fabrication simulator \(\hat g_\theta\) so that inverse design solves
\[
M^\star(x,y) = \operatorname*{arg\,min}_{M(x,y)} \mathcal{L}_{\text{optics}}\!\left(f_{\text{optics}}\!\left(\hat g_{\theta}(M(x,y))\right)\right),
\]
rather than assuming \(g(M)=M\). The simulator is trained from **96** layout–print pairs fabricated on a **Nanoscribe Photonic Professional GT2** two-photon lithography system and measured by AFM, with a **72/24** train/validation split. Its PBL architecture decomposes photolithography into optical exposure, thresholding, diffusion, shrinkage, and mismatch correction; quantized \(13\)-level layouts are handled with the Gumbel–Softmax reparameterization. The reported PBL forward prediction error is \(\mu = 24.35\,\text{nm}\), while printed flat-surface roughness is \(\sigma = 15.96\,\text{nm}\) and clean substrate roughness is \(\sigma = 1.25\,\text{nm}\). The framework improves fabricated HOEs and MDLs relative to conventional design and to simpler lithography-aware baselines [2309.17343].

The large-area extension makes this fabrication-aware view explicitly mass-production compatible. In that pipeline, a master DOE is fabricated by direct-write grayscale lithography using **positive AZ® 4562 photoresist** on a **Heidelberg Instruments DWL 66+**, then replicated by nanoimprint lithography in **UV-curable OrmoComp resin** using an **Obducat Eitre 3**. The optical design is reformulated from
\[
h^* = \arg\min_h \sum_\lambda \mathcal{L}_p\!\left(\left\|E^{\text{dest}_\lambda(h)}\right\|^2\right)
\]
to
\[
l^* = \arg\min_l \sum_\lambda \mathcal{L}_p\!\left( \left\|E^{\text{dest}_\lambda\bigl(\mathcal{G}\{l\}\bigr)}\right\|^2 \right),
\]
where \(\mathcal{G}\) is a fully differentiable learned map from lithography layout \(l\) to fabricated height map \(h\). The model is super-resolved: calibration maps a \(1\,\mu\mathrm{m}\)-sampled design layout to AFM-measured topography at about **200 nm** sampling pitch, i.e. **\(5\times\)** spatial super-resolution. On held-out patterns, the neural lithography model reaches **35.20 dB PSNR** and **2.45% NRMSE**, compared with **27.36 dB PSNR** and **7.27% NRMSE** for an MTF-based physical model. The system is combined with distributed FFTs, tensor-parallel convolutions, **JAX**, and **GSPMD**, allowing simulation grids up to **128,640 \times 85,760** for a **32.16 mm \(\times\) 21.44 mm** DOE on **16 A100 GPUs**. Fabrication-aware optimization yields a nearly speckle-free coherent hologram, a beam splitter with **53% higher overall spot intensity** than the conventional design, and a single-DOE broadband imaging system whose measured PSF is close enough to the simulated PSF to support one-step Wiener restoration [2505.22313].

Taken together, these papers redefine differentiable lithography as fabrication-aware co-design: the optimized variable is no longer the nominal structure alone, but the input actually sent to the fabrication system, and the objective is evaluated on the predicted as-fabricated geometry.

## 5. Photonics, EUV masks, and volumetric lithography

In photonic integrated circuits, PRISM adapts inverse lithography to the fact that geometry fidelity is not a sufficient proxy for optical function. Its differentiable fabrication twin factorizes printed geometry into a continuous exposure field and a soft threshold:
\[
W=\Phi_{\theta = (\theta_I,\theta_T)}(M) =\psi\!\bigl(I_{\theta_I}(M)-T_{\theta_T}(M)\bigr), \qquad \psi(x)=\frac{\tanh(2\beta x)+1}{2}.
\]
For DUV, the exposure model is Hopkins/SOCS-based; for EBL, it is PSF-based. The ILT objective is not plain contour matching but a weighted \(\ell_4\) wafer loss,
\[
\mathcal{L}_{\mathrm{ILT}} = \|\Gamma \odot (W-W^{*})\|_4^4,
\]
where the weight map is derived from normalized photonic adjoint sensitivity,
\[
\Gamma = \left|\frac{\partial \mathrm{FoM}}{\partial M_0}\right|,
\]
normalized to \([0.5, 1.5]\). Under **193 nm DUV**, PRISM-SOCS improves yield\(_{80\%}\)/yield\(_{90\%}\) for crossing from **11%/0%** to **100%/100%**, for optical diode from **0%/0%** to **90%/85%**, and for \(1\times 2\) MDM mux from **0%/0%** to **98%/98%**; harder devices such as WDM mux remain challenging even after correction [2602.15762].

For EUV masks, a distinct line of work treats the full diffraction engine itself as the differentiable object. The waveguide-method framework solves inverse lithography by differentiating through a rigorous layered periodic Maxwell solver. The absorber permittivity is parameterized by a density field,
\[
\varepsilon(x)=\varepsilon_{\rm v} + \rho(x)\left(\varepsilon_{\rm ab}-\varepsilon_{\rm v}\right),
\]
with either pixel-wise or Fourier parameterization, and optimized against
\[
\mathcal{L}(\varepsilon)=\left.\left(\left|E^{(r)}\right|^2-\left|E^{(d)}\right|^2\right)^2\right|_S.
\]
The experiments use a realistic EUV mask at \(\lambda = 11.2\) nm with a **30 Ru/Be/Sr period** multilayer, **60 nm** absorber thickness, and **6^\circ** incidence angle. For absorber materials, the paper reports normalized central peak intensities \(\max I / I_0\) of **0.64** for TaBN, **1.99** for La, and **0.45** for U; La gives the strongest central maximum, while U gives the closest overall field match. A 3D extension is also demonstrated [2606.25753].

Computed Axial Lithography provides a volumetric analogue. There, the standard backprojected dose \(R^*g\) is augmented by a blur kernel,
\[
\tilde{R}^{*}g(x,y,z) = R^*g(x,y,z) * K(x,y,z),
\]
and the adjoint is modified consistently:
\[
\tilde{R} f(r,\theta,z) = \int_{L(r,\theta,z)} \left(f(x,y,z) * K^*(x,y,z)\right)\, ds.
\]
The kernel is fitted experimentally from micro-CT data and uncorrected prints, with a notional diffusivity near
\[
D = 1\times 10^{-10}\ \text{m}^2/\text{s},
\]
while corrected-print trends suggest a peak near
\[
D \approx 3.4\times 10^{-10}\ \text{m}^2/\text{s}.
\]
The paper shows that co-optimizing projections under this blurred forward model improves fidelity over Richardson–Lucy target deconvolution; for the RL baseline it reports **ASSD \(=42.5\,\mu\text{m}\)** and **MS-SSIM \(=0.696\)** [2606.04174].

These extensions show that differentiable lithography is not limited to planar semiconductor mask optimization. It also encompasses photonics-informed mask correction, rigorous EUV electromagnetic inverse design, and volumetric dose optimization under experimentally identified blur.

## 6. Limitations, misconceptions, and research directions

A persistent misconception is that any gradient-based design related to lithography is “differentiable lithography.” The literature is more specific. The high-NA metalens study for zone-plate-array lithography performs adjoint electromagnetic optimization of a lithography-critical optical component and reports **85.50%** transmission normalized focusing efficiency at **0.60 NA** and **405 nm**, but it does not differentiate through resist chemistry, pattern fidelity, or lithographic process transfer; it is best classified as inverse-designed lithography hardware rather than end-to-end differentiable lithography [2212.06441]. LithoDreamer likewise provides differentiable latent rollouts and inverse planning in a world-model framework, but its “physics-informed” character derives from stage decomposition, latent subspace priors, and process conditioning rather than from explicit Abbe, Hopkins, or resist equations [2606.26713].

Another recurring limitation is process specificity. The learned fabrication simulators in computational optics are calibrated on small, process-specific datasets: **96** layout–print pairs in the TPL-based neural lithography work, **20** fabricated random patterns in the large-area grayscale-lithography study, and compact but process-specific calibration reticles in PRISM [2309.17343] [2505.22313] [2602.15762]. Those papers explicitly note intrinsic stochastic variability, process dependence, and the possibility of out-of-distribution corrected masks. PRISM argues that physics-grounded twins are often more reliable than unconstrained neural predictors under limited data and OOD ILT masks, especially for DUV [2602.15762].

Computational cost remains another fault line. The large-area diffractive-optics pipeline requires custom distributed operators and multi-GPU infrastructure; the rigorous EUV mask solver differentiates through eigenproblems, scattering-matrix assembly, and field reconstruction; CAL co-optimization introduces additional FFT-based object-space convolution and experimentally fitted blur. These costs do not invalidate differentiability, but they move practical performance questions from pure algorithm design to systems engineering [2505.22313] [2606.25753] [2606.04174].

A further limitation is model fidelity at the chemistry and process-transfer level. Semiconductor imaging frameworks centered on scalar Abbe or Hopkins models remain strongest on optics and mask/source optimization, while resist modeling may be simplified, abstract, or delegated to compact surrogates [2409.15306] [2502.06838]. In computational optics, the fabrication-aware models capture effective fabrication well enough for demonstrated tasks, but they remain empirical surrogates rather than mechanistic simulators with explicit latent variables for beam shape, material batch variation, or stochastic development fluctuations [2309.17343] [2505.22313]. This suggests that future progress will likely depend on combining process-calibrated differentiable twins with stronger physical priors and, where available, hardware-in-the-loop correction.

The process chain itself may also change. Resistless EUV patterning on HF-treated Si(100) replaces the usual photoresist route with EUV-induced surface oxidation and selective TMAH etching, achieving SiO\(_2\)/Si gratings with **75 nm half-pitch** and **31 nm height**. The paper is experimental rather than differentiable, but it shows that the relevant forward model could shift from “aerial image \(\to\) resist chemistry \(\to\) development” to “aerial dose \(\to\) oxidation-state field \(\to\) etch-resistance field \(\to\) topography” [2310.01268]. This suggests a broader research direction: differentiable lithography need not be tied to chemically amplified resist models if the manufacturing physics itself changes.

Across these lines of work, the unifying idea is stable: lithography or fabrication is treated as an optimization-native physical process rather than as a post hoc verification step. What differs is the differentiable state variable—pixels, edges, level sets, spline control points, absorber permittivity, projection images, or fabrication layouts—and the fidelity of the forward model through which gradients are propagated.

Source: https://www.emergentmind.com/topics/differentiable-lithography