WorksPhysics and Astronomy
Writing with Near-Field: towards Atomic Resolution in Material Processing
Abstract
Direct laser writing with ultra-short pulses is the most flexible route to two- and three-dimensional (2D/3D) material modification, yet its resolution is judged against the optical diffraction limit. This review argues the opposite: the finest, most reproducible modifications come not from shortening the wavelength but from localising the absorber. Through the threshold effect and a transient epsilon-near-zero (ENZ) state - the dielectric-to-metal (Die-Met) state of the excited solid - the deposited energy is confined to a dipole-like nano-seed whose near-field intensity falls as r^-6 (ideal point dipole). The localisation is thus set by the seed size, not the focal wavelength, and the final feature is shaped by the material response: post-pulse heat flow, cavitation and quenching. We show that polarisation controls the anisotropy of deposition: the quasi-static near-field and interface boundary conditions select the writing direction and the modification orientation, down to a single defect. We survey energy-deposition mechanisms from 100 nm to the single-nanometer scale and link them to working engineering: stitch-free large-area patterning and record aspect-ratio (10^3-10^4) dicing and wafering of transparent crystals. Finally, we conjecture guided volumetric deposition, in which laser-inscribed lossy nano-seeds act as antennas that let cavity microwaves anneal the modified nano-volumes from within, as a 2.45 GHz oven heats food through water-dipole absorption. Throughout, atomic resolution is used in a quantum-confined sense and denotes a direction of travel: 10-20 nm features are demonstrated, while deterministic single-defect and single-bond modification remain goals. The thesis is that ultrafast direct write beats diffraction-limited intuition because the material itself forms a transient nanoscale absorber whose near-field and boundary conditions govern energy deposition.
Introduction: The Two Moore’s Laws
The electronics (e-)Moore’s law of doubling the transistor count roughly every two years[1] has driven computing, telecommunications and, most recently, artificial intelligence, and is built on parallel, planar 2D nanolithography. Its optical (o-)Moore’s counterpart describes the doubling of the average power (photons per unit time) of ultra-short pulsed lasers (USPL), a trend that has held for the past years[2, 3]. The fundamental distinction between the two laws is that o-Moore’s scaling powers sequential direct laser writing (DLW) and material processing, whereas e-Moore’s success rests on parallelism[4]. The question this review addresses is when, and by what physics, sequential direct write can challenge parallel-mode lithography and how the energy of an ultra-short pulse can be localised to a nanoscale volume that is essentially independent of the illuminating wavelength.
The state of the art in commercial lithography is achieved with 193 nm (ArF) water-immersion and 13.5 nm extreme-ultraviolet (EUV) sources. EUV requires vacuum enclosures, reflective masks, and chemically amplified resists whose stochastic line-edge roughness limits resolution[5], and entails a cost of order $700 million per node change. Crucially, the photomasks themselves are written by electron-beam lithography (EBL), that is, in a direct-write, sequential mode. The energy efficiency of projection tools is poor: MW of wall-plug power, used to ablate Sn micro-droplets with a CO laser, yields only W of useful exposure on the wafer—an end-to-end efficiency of , while a planar CMOS flow requires processing steps per wafer. These realities motivate a paradigm shift toward maskless, energy-efficient, 3D-capable fabrication—a “3D-printing-like” democratisation of high-end micro- and nano-electronics, photonics and large-area patterning for photovoltaics[4, 6].
At the leading edge of semiconductor manufacturing, lithographic resolution is no longer the only limiting variable; the rate at which that resolution can be deployed has become equally consequential. EUV is required for selected critical layers in approximately 5-nm-class and smaller technologies, yet its availability is concentrated in a single commercial scanner supplier and a highly specialised subsystem network, including sole-source dependencies in critical projection optics.[7, 8, 9] The transition from 0.33-NA EUV to 0.55-NA High-NA EUV is intended to extend single-exposure scaling and reduce reliance on increasingly complex multipatterning, but industrial insertion requires the scanner, source, mask, resist, metrology, overlay and process-integration ecosystems to mature together.[10, 11] This qualification burden creates a substantial interval between demonstrated optical resolution and broadly replicated manufacturing capability. The resulting bottleneck is therefore not simply the minimum printable pitch, but the limited rate at which a complete, capital-intensive lithographic ecosystem can be qualified and distributed across the manufacturing base.
Two requirements define the frontier. First, feature size: a single modification as small as nm approaches the “atomic” or quantum scale, since a single dopant atom in a crystal, e.g. a P atom in Si, behaves as a quantum dot of nm diameter[12, 13]. We use “atomic/quantum scale” in this quantum-confined sense, distinct from atom-by-atom placement, and throughout separate feature size, placement/pitch, defect determinism and atomic selectivity, reserving “atomic” for cases with single-defect or single-bond evidence. By this usage, true atomic resolution in the crystallographic sense ( nm) remains beyond optical direct write; what is within reach, and demonstrated, is the quantum-confinement scale of –10 nm together with single-defect photophysics[13]. Second, resolution between two adjacent features, which sets the packing density. Both are required for basic quantum devices (single-photon emitters, single-electron transistors) and both, this review argues, are governed by the spatial localisation of energy deposition rather than by the wavelength alone.
The remainder of the paper is organised as a sequence of self-contained topics: far-field focusing (Sec. II), the threshold effect (Sec. III), a critical survey of deposition mechanisms (Sec. IV), the dipole/ENZ near-field that unifies 2D and 3D writing (Sec. V), the role of polarisation relative to the scan direction (Sec. VI), engineering that already works (Sec. VII), and a new laser-seeded microwave paradigm (Sec. VIII), the latter explicitly flagged as a forward-looking, hypothesis-level element of this review-perspective.
Far-Field Focusing and the Diffraction Barrier
For propagating (far-field) light the Heisenberg uncertainty relation (with the reduced Planck constant and the position and transverse-momentum spreads) translates into the waist-divergence relation ( the wavelength), fixing the trade-off between the lateral waist and the depth of focus (DoF)[14] [Fig. 1(a)]. A diffraction-limited spot is for near-infrared pulses at moderate numerical aperture (NA), and the aspect ratio is therefore bounded: lateral resolution and aspect ratio cannot be improved independently. A Bessel beam, synthesised with an axicon, replaces the short DoF with a long, non-diffracting core, raising the achievable aspect ratio while leaving the same uncertainty trade-off in force in the transverse plane[15].
The essential point for this review is that the focal spot is only the envelope within which energy is delivered. The size of the modified region is set not by but by where the deposited energy density crosses a material threshold [Fig. 1(b)], a distinction that the diffraction limit, framed for linear imaging, does not capture. Two further features of a tight focus matter below. First, at high-NA the focus is non-paraxial and acquires a longitudinal field component along the propagation axis, opening a resonant-absorption channel that is absent in the paraxial picture. Second, the same nonlinearity that sharpens the deposited dose also makes the process deterministic: only the peak of the Gaussian, where the intensity exceeds the modification-threshold intensity , is written.
The threshold effect: A tens-of-nm seed from a micron focus
The USPL modification is governed by a strongly nonlinear energy-deposition rate , with effective order –. The microscopic channels behind this phenomenological exponent were established in the classical analyses of photoionisation by Keldysh[16] and of laser-induced breakdown by Bloembergen[17]: the relative weight of multiphoton and tunnelling ionisation is set by the Keldysh parameter (bandgap , field amplitude ), which is of order unity at the – intensities relevant here, so the two channels coexist rather than compete cleanly. The carriers they seed are then multiplied by avalanche (impact) ionisation on the trailing edge of the pulse. The exponent is therefore an effective, intensity-range-dependent slope of the deposition rate, not a photon count. For a Gaussian profile (peak intensity , radius , radial coordinate ) and on-axis modification threshold , a threshold on the deposited dose is equivalent to a threshold on the local intensity, because is a monotonic function of . The modified radius is therefore
Equivalently, , where is the peak-to-threshold dose ratio. Expressed through the experimentally controlled ratio (equivalently ), the radius therefore does not depend on the nonlinear order . The nonlinearity acts elsewhere: it narrows the dose profile itself, so the dose gradient at the modification contour is times steeper (a sharper, more deterministic edge) and no dose is deposited outside the focal volume. This is consistent with the single-shot ablation study of Garcia-Lechuga et al.[18], who showed experimentally, across bandgaps of 3.8 to 8.8 eV and wavelengths of 515 and 1550 nm, that ablation craters map the beam contour at a strict threshold intensity irrespective of the nonlinear order. For , Eq. (1) gives at and at . A radius of 30 nm requires to be within 0.7% of , and 10 nm within 0.08%. The same logarithmic dependence simultaneously reduces the depth of the modified volume [Fig. 1(b)]. Deterministic ablation can therefore reach the tens-of-nanometre regime in both lateral and axial directions[19, 20, 21], provided that the pulse energy and focus are stable at the level or better. The energy required is modest: any phase transition: melting, evaporation, ionisation is induced with 0.01 eV to 10 eV per atom[22]. This is the first key idea of the review: the relevant length scale is set by where the deposited energy exceeds a material threshold, not by .
The threshold and near-field analyses of this review localise the deposited energy; the final modification is additionally broadened by post-pulse heat flow during electron-lattice relaxation and resolidification. A schematic scaling is
with the thermal diffusivity and the time over which the seed stays above the modification temperature (bounded from below by the electron-phonon relaxation time of –). In fused silica () heat spreads only in 10 ps and in 1 ns: a 10-nm deposited profile survives as a 10–20 nm feature provided quenching is sub-nanosecond, which confined microexplosion provides (). In crystalline silicon () the same 10 ps already gives nm of spreading, which is why the finest demonstrated features in Si rely on thin films, interfaces or stress/cavitation-mediated (rather than purely thermal) transformation[25]. Localising the deposition is therefore necessary but not sufficient: the quench rate decides how much of the localisation survives, and Eq. (2) gives the practical guidance: low- hosts (glasses) hold nanoscale features robustly, high- hosts (Si, metals) require confinement or ultrafast mechanical transformation.
The threshold effect alone, however, only carries the modification to before the logarithm in Eq. (1) saturates and shot-to-shot fluctuations dominate the edge definition. Reaching the single-nm scale requires a second mechanism that concentrates the field after the seed is formed, the dipole/ENZ near-field of § 5.
Energy deposition from 100 nm to the single-nanometer scale: A critical survey
We now scout the principal mechanisms by which energy is deposited into matter on scales from 100 nm to the single nanometer. Table I summarises the approaches, their demonstrated localisation, the underlying mechanism, and, most importantly, their dependence on wavelength. The trend we emphasise is that the finest, most wavelength-insensitive localisation is obtained when energy is deposited through a near-field or absorber-defined mechanism rather than by linear far-field absorption. This is not a claim of universal superiority. Competing high-resolution routes include electron- and ion-beam lithography and milling (EBL, FIB, and He/Ne-ion methods), electron-beam-induced deposition, plasmonic-antenna and near-field scanning optical lithography, microsphere-assisted structuring, AFM thermochemical or local-oxidation writing, scanning-tunnelling atom manipulation, and directed block-copolymer self-assembly. Each route wins on particular metrics. Recent surveys map this competitive landscape[26], as does the dedicated perspective on ultrafast-laser volume nanostructuring by Stoian[27], while microsphere-assisted schemes have reached feature sizes in ambient air[28]. Several far-field demonstrations now sit firmly in the 10-nm class and must be part of any balanced survey: dual-beam overlapping with orthogonal polarisations has produced –12 nm () features on Si surfaces in ambient air[25]; cavity-seeded evanescent-field structuring with Bessel beams has reached 7 nm () trenches of micrometre depth in glass[29]; features were written in van der Waals NbOI films[30]; and nanohole-clad waveguides with depth-to-diameter ratios up to have been written in single crystals[31]. Notably, the mechanisms invoked in these works - self-generated near-fields seeded by voids or cavities, thermal-stress localisation, and interference of incident with back-scattered light - overlap with, but are not identical to, the ENZ/dipole picture developed below. § 3 B weighs these alternatives explicitly. Our argument concerns the regime in which a transient, self-formed absorber governs the resolution.
Far-field nonlinear deposition
The pioneering demonstrations of sub-diffraction writing exploited nonlinear absorption. Two-photon polymerisation reached 120 nm features below the diffraction limit already in 2001[32], building on the discovery that tightly focused femtosecond pulses write permanent refractive-index changes and store data inside transparent solids[34, 35]. The mechanisms and parameter space of femtosecond micromachining in transparent materials were comprehensively reviewed[33]. Borrowing the stimulated-emission-depletion (STED) concept from super-resolution microscopy[37], STED-inspired lithography pushed polymer feature sizes to nm[38], the resolution improving as . These methods remain nonlinear but wavelength-referenced and, in the STED case, demanding in dual-beam overlap and resist chemistry; the damage/polymerisation threshold itself depends on wavelength and pulse duration, as quantified for 3D nanopolymerisation[36].
Alternative routes: cavitation, scattering centres and the material response
A balanced reading of the recent literature shows that the transient-absorber picture of this review is one member of a family of mechanisms, and we state the alternatives explicitly.
Threshold without nonlinearity. As noted in § 2, Garcia-Lechuga et al.[18] demonstrated that single-shot ablation contours map the beam profile at a strict threshold intensity independent of the absorption order. The nonlinearity therefore sharpens the dose, but the edge definition of an ablation crater can be a purely threshold-contour property. Eq. (1) shows that this holds for any monotonic dose threshold: at a given the contour radius does not depend on . The nonlinearity sharpens the dose gradient at the contour and confines absorption to the focal volume, but it does not by itself reduce the modified radius below the intensity-threshold contour.
Nanocavitation and rarefaction. A void or rarefied (low-density) region formed by pressure-gradient cavitation under strong confinement and fast quenching is itself a low-permittivity scattering centre, with no requirement that the carrier density ever approaches the critical value. The “slow microexplosion” in bulk silica[24] and the cavity-seeded evanescent-field structuring of Zhang et al.[29] show that, once such a centre exists, it sustains scattering and near-field confinement during subsequent writing. The same electromagnetic role that the ENZ/Die-Met seed plays in § 5, reached by a mechanical rather than electronic route. Because quenching can be arbitrarily fast under confinement, and mechanical processes such as interfacial sliding have interatomic action range, this pathway has no obvious lower length-scale bound. The two routes are complementary: the internal-field enhancement of Eq. (15) depends only on the permittivity contrast, however generated. The ENZ transit is thus a sufficient, self-sharpening route to a low- seed within a single pulse, but it is not a sine qua non for near-field-localised writing.
Surface instabilities and metals. Sub-surface cavitation and hydrodynamic instabilities, observed for instance in ripple formation, provide local field and energy-deposition enhancement even in metals, where free carriers pre-exist and no dielectric-to-metal conversion is needed. Polarisation-coupled scattering at such corrugations is the accepted origin of LIPSS. Our framework addresses the complementary case of initially transparent media, where the transient conversion is the self-sharpening step.
Transformation mechanisms. Field localisation fixes where energy arrives, i.e., the conversion into a permanent feature proceeds by melting and resolidification, densification/rarefaction, cavitation, defect accumulation or, in the extreme normal-field case, electrostatic (Coulomb-explosion-type) ion removal.[42] The latter requires the surface charge to survive longer than the ion response time before screening by conduction electrons neutralises the surface charge and suppresses Coulomb explosion, leaving thermal/hydrodynamic routes to dominate.[22] Experimental verification of the layer-by-layer electrostatic regime remains outstanding (Sec. VI).
The atomic and sub-atomic limits
Beyond nm, deterministic single-defect inscription becomes the goal. Polarisation-sensitive strong-field (tunnelling) ionisation[13, 41, 43, 44] can in principle direct bond breaking and create individual optically/magnetically active defects. The role of the tunnelling channel is that, for Keldysh parameter , the ionisation rate follows the instantaneous field direction, so bonds aligned with are ionised preferentially. Directionality of multiphoton absorption, which samples the cycle-averaged intensity, largely lacks in atomically thin h-BN the uncertainty was found to be set by the thermal motion of atoms rather than by the gradient of the optical envelope, and the emission turned from stochastic/polychromatic to deterministic/monochromatic.[13] At relativistic intensities () two distinct frontiers open. The
References
first is structural: confined micro-explosion synthesizes high-pressure phases and superdense materials that are quenched to ambient conditions: tetragonal Si polymorphs and superdense aluminium[45, 46, 47], a route to new phases. The second is genuinely nuclear and defines the ultimate single-atom limit of material processing: transmutation, the change of one element into another. Few-cycle relativistic pulses accelerate deuterons that drive the fusion reaction, and compact Extreme Light Infrastructure (ELI)-class sources have demonstrated laser-driven D-D fusion with fast-neutron yields already at the n/s level and a credible roadmap to n/s at kHz repetition rate using few-cycle pulses[48, 49]. Whereas every method above localises energy to a small volume, transmutation localises identity to a single nucleus, the true atomic endpoint of the resolution hierarchy developed in this review. The near-field and absorber-defined mechanisms, to which we now turn, occupy the crucial middle ground between these limits and the diffraction-bound far-field methods.
The unifying mechanism: dipole / ENZ near-field localisation
The concept that unifies the wavelength-independent entries of Table I is that a sub-wavelength ionised region of reduced permittivity behaves as an induced dipole whose near-field is strongly localised. Two terms recur below and are not synonyms: Die-Met denotes the entire transient conversion of the excited solid from dielectric () to metallic () response as the carrier density grows[22]. The ENZ denotes the single waypoint on that trajectory, reached near (not at) the critical density once the background permittivity, collision frequency and effective mass are accounted for. Each carrier-density milestone on the trajectory has distinct physical properties, catalogued below and, as a time sequence, at the end of this section.
| Approach | Localisation | Deposition mechanism | -dependence | Refs. |
| 193 nm immersion photolithography | 38 nm half-pitch | linear absorption in resist; phase + immersion | strong (/NA) | [5] |
| EUV (13.5 nm) projection | 13 nm | linear absorption; reflective mask, vacuum | strong | [5] |
| Electron-beam lithography | sub-10 nm | electron–resist energy loss | weak (de Broglie) | [4] |
| Multiphoton polymerisation (2PP/TPL) | 100–120 nm | -photon absorption above polymerisation threshold | nonlinear, sub- | [32, 33, 34, 35, 36] |
| STED-inspired DLW | 50–55 nm | depletion of excited state, | nonlinear | [37, 38] |
| Microsphere-assisted fs structuring | 38 nm (/27) | photonic-nanojet concentration; thermal-stress nanogrooves | moderate | [26, 28] |
| Dual-beam far-field patterning (Si) | 10–12 nm (/66) | orthogonal-polarisation double-beam overlap; thermal-stress localisation | moderate | [25] |
| Cavity-seeded evanescent structuring (Bessel, glass) | 7 nm (/100) | far-field-induced cavities seed evanescent near-fields; local thermal ablation | none (near-field) | [29] |
| Deterministic threshold ablation | tens of nm | above-threshold MPI/avalanche on Gaussian tip; contour mapping at strict threshold | weak | [19, 20] |
| Quantized interference structuring | threshold + standing-wave interference | moderate | [21] | |
| Controlled-avalanche ripplesingle-groove ablation (sapphire) | 50–70 nm (/6) | near-threshold avalanche; sub-surface plasma sphere-to-plane growth; ripple pattern collapses to one groove | moderate ( seed) | [39] |
| Far-field-induced near-field breakdown (O-FIB, “nano-pen”) | 20 nm | boundary-condition -field enhancement (normal continuity) | none (near-field) | [40] |
| Bessel nano-seed merging / super-stealth dicing (“nano-knife”) | 10 nm, AR | near-field merging of ionised seeds along filament | none | [14, 15] |
| ENZ / Die-Met dipole localisation | tens single nm | induced dipole, , | none | [4, 22] |
| Nanocavitation / void-seeded scattering | 10 nm | pressure-gradient cavitation, rarefaction and fast quenching; void sustains field confinement | none | [23, 24, 29] |
| Tunnelling-ionisation defect inscription | atomic / quantum | polarisation-directed bond breaking | none | [13, 41] |
| Laser-seeded microwave volumetric deposition (this work) | tens single nm | laser defines lossy nanovolume; cavity microwave Joule-heats it; thermal runaway | none (absorber-defined) | — |
Table I. Critical survey of energy-deposition methods from 100 nm to the single-nanometer/atomic scale. “-dependence” indicates how strongly the achievable localisation is tied to the illuminating wavelength.
The induced dipole and the law
Consider the ionised seed as a sphere of radius and (Drude) permittivity
embedded in a host of permittivity . Here is the optical angular frequency, the bound-electron (background) permittivity, the vacuum permittivity, the plasma frequency, the elementary charge, the laser-generated carrier density, the collision frequency and the carrier effective mass. In the quasi-static limit the sphere acquires a dipole moment with the Clausius-Mossotti factor[50, 51]
It is essential to follow the real part of the permittivity as the laser-driven carrier density grows during the pulse. From Eq. (3), falls monotonically from its bulk value. It first crosses below unity into the epsilon-near-zero (ENZ) window . The consequences for energy deposition are real but geometry-dependent and should not be reduced to a single velocity argument: near ENZ the wave impedance, the internal field distribution, the absorption length and the (dispersive, loss-corrected) energy density are all strongly modified, and ENZ media are known to enhance local fields and light-matter interaction under appropriate loss, angle and boundary conditions[22, 52, 53]. We do not invoke a simple slow-light relation (valid only for an idealised lossless, non-dispersive medium), nor a universal interface reflection minimum at a bulk dielectric/ENZ interface. The reflectance can be high unless loss, finite thickness or impedance matching intervene. The endpoint is doubly significant: it is the perfect-absorber, dielectric-to-metal (“Die-Met”) condition at which the deposited-energy density peaks[22], and it is reached near the critical density. For the Drude form, requires , i.e. with the bare vacuum-plasma value . At () the bare value is (for , ); in a real dielectric (–) the ENZ density is correspondingly higher and depends on , and band filling. Either way it is comparable to the carrier density of the laser-formed seed. The useful operating window is therefore the controlled descent , terminating at the ENZ/breakdown point. Past it (), the seed turns metallic and reflective, and the Fröhlich resonance lying in that regime makes the Clausius-Mossotti factor Eq. (4) diverge (limited only by losses) as sweeps transiently through it. Because the internal field of a sphere is , approaching ENZ concentrates the field inside the seed (by a factor in the limit) and suppresses the usual reflection phase so that incoming energy couples efficiently. The dipole near-field outside the seed[50],
yields an intensity
the quasi-static dipole field decays as , so intensity falls as . The prefactor depends on the seed polarisability, size and loss, so this is a localisation scaling (the dose width is a fixed fraction of the seed size , a length tied to , not to ), not a free-standing, resolution law[4]. We stress that the ENZ condition (, internal field only ), optical breakdown (the critical density), the Fröhlich dipole resonance (, a lossy/metallic state with large but loss-limited enhancement) and the perfect-absorber condition are distinct states at different carrier densities, not one continuous favourable regime. Moreover, as developed in § 3 B, the electromagnetic role of the low- seed can equally be played by a cavitation void or rarefied region generated mechanically: the ENZ transit is a self-sharpening route to the seed, not a prerequisite for the near-field localisation itself. Combining Eq. (6) with the threshold relation Eq. (1), a focus operated near threshold produces a -30 nm, in agreement with the demonstrated nano-pen and nano-knife features[14, 40].
Temporal sequence during the pulse and how to probe it
For a near-threshold pulse of duration fs the seed traverses these states in a definite order: (i) during the leading edge, multiphoton/tunnelling ionisation raises and descends from through unity via weak perturbation, little deposition; (ii) around the pulse peak the ENZ window is crossed, impedance mismatch collapses, the internal field rises to [Eq. (15)] and the volumetric deposition rate peaks. We expect the dominant share of the absorbed energy to be deposited in this stage. Then, (iii) if ionisation continues, sweeps through the Fröhlich value , giving a transient, loss-limited spike of the external near-field. Followed by (iv) deeper into the metallic branch the seed reflects, deposition migrates to a skin layer and to the external near-field lobes, and self-limits. The sequence is experimentally accessible: time-resolved reflectivity/transmission pump-probe traces the descent, THz emission tracks the transient currents in the focal volume[54], and four-polarisation near-field imaging resolves the anisotropy of the final deposition[55, 56]. A measured crossing time of the reflectivity minimum (ENZ) relative to the pulse peak would directly test stage (ii), the central claim of this section.
2D writing on surface: the nano-pen
The boundary conditions make this concrete in 2D. The continuity of the normal component of the displacement field, [50], enhances the -field on the low-permittivity side of an interface. This applies to lossless dielectrics; at a conducting interface the relevant continuity is that of the normal current density. Charge conservation across the boundary gives, for harmonic fields, continuity of the normal total (conduction + displacement) current,
i.e. continuity of the normal component of the complex permittivity [50]. In the quasi-static conduction limit Eq. (7) reduces to : the normal field is enhanced on the low-conductivity side, the conduction-current twin of the dielectric law, with the loss tangent deciding which term dominates.
Silicon and copper bracket the two regimes. At a 1030 nm optical frequency (), copper () has (and at 2.45 GHz): it is metallic at every relevant frequency, the field is expelled into a nanometric skin depth and reflected; the enhancement lives entirely in the adjacent dielectric. Intrinsic silicon (, ) has optically and behaves as a near-ideal dielectric, so the field penetrates and the dielectric law applies. This contrast is precisely why silicon is an especially favorable direct-write medium: it starts on the dielectric branch, and the pulse converts the focal nanovolume into a conducting ENZ/Die-Met seed. rising by as , so the seed flips locally onto the metallic branch where it was written. The key requirement is a large transient permittivity/conductivity change during the pulse, which yields a self-sharpening seed. Metals (skin-depth absorption, plasmonic field concentration, electron-phonon nonequilibrium) and passive dielectrics (high-NA nonlinear absorption) also localise energy, but by different, less self-sharpening routes. The self-induced near-field localisation emphasised here is strongest in media that traverse within the pulse.
In far-field-induced near-field breakdown (O-FIB), a diffraction-limited pulse first ablates a nm seed hole; at a reduced fluence (, i.e. of threshold) only the nanoscale vicinity of that hole is field-enhanced, so a nm groove is written by steering the far-field beam[40]. This “nano-pen” writes 2D etch-mask patterns in a dielectric nanofilm in ambient air, without vacuum or high-precision stages, and, decisively, its resolution is set by the interface enhancement, not by the wavelength. Crucially, the direction of writing is set by the laser polarisation, the subject of Sec. VI.
3D writing inside volume: the nano-knife
Three-dimensional recording, direct laser writing of structures within the volume of a transparent solid, is where the near-field, wavelength-independent localisation pays off most. It first reached the true nanoscale with Bessel-beam “super-stealth dicing”[Fig. 1(c)], and the same dipole/near-field mechanism carries over to an ordinary Gaussian (cylindrical) focus[57]. The same physics as for 2D surface patterning extends into the bulk. Along a Bessel-like filament, nanoscale seed modifications generated at the intensity maximum reduce the local permittivity; incident light back-scattered by the seeds interferes with the incoming beam in a “crawling” positive-feedback loop that homogenises the longitudinal deposition and confines the transverse field, merging the seeds into a continuous nano-plane of nm width[14][Fig. 1(c)]. A Gauss-law identity for the Bessel focus places the top 25% of the intensity within a nm-wide channel, the seed of the cut. This “super-stealth dicing” achieves a record aspect ratio length/width of order (with the longest filaments approaching )[14], complementary to axicon-doublet ultra-high-aspect machining in silicon[15], to the cavity-seeded -class trenching of Zhang et al.[29], and to nanohole-clad waveguides with depth-to-diameter ratios up to . Here the 3D localisation is again governed by the self-generated low-permittivity seed, a dipole-like absorber, rather than by . Why the energy concentrates at the reduced-permittivity centre and how the polarisation steers it for an ordinary Gaussian focus as well as a Bessel filament is developed in Sec. VI.F.
Polarisation control of dipole energy deposition
The seed oscillates as a dipole , so the deposited energy is polarisation-dependent. At a stationary focus the quasi-static near-field is strongest along , orienting the inscribed defect even without scanning (Figure 2). During scanning, the writing of an extended groove (on surface) is governed by the quasi-static near-field and the electrostatic boundary conditions at the groove wall, which place the strongest deposition transverse to (Figure 3). Both regimes are near-field; the far-field dipole radiation pattern is not the deposition mechanism. We treat the cases in turn.

Figure 2. Local field at a stationary focus. (a) The induced dipole’s quasi-static near-field Eq. (5) is twice as strong along (poles) as across it (equator), a 4:1 intensity ratio Eq. (9). A defect or modification inscribed without scanning is expected to be elongated along the polarisation. (b) FDTD simulated field maps in the XY-plane containing the polarisation axis () for three sphere permittivities (SiO host, , nm, nm). As decreases toward zero (ENZ), the field concentrates progressively inside the sphere, consistent with Eq. (15).

Figure 3. Polarisation and the written groove. (a) Profile of the line written by the nano-pen: a groove of width nm in a 20 nm titania film on silica, its width set by near-field boundary-condition enhancement and the seed size, not by [40]. (b) FDTD simulated field map ( groove, , groove nm, nm) in the plane: the normal field is strongly enhanced inside the low-permittivity air groove, with hotspots at the groove corners. (c) Horizontal profile of across the groove wall at mid-height ( nm), normalised to the TiO background. The measured peak intensity ratio across the groove wall is 36.22 (air side , TiO side ), consistent with the electrostatic prediction [Eq. (11)]; the shortfall reflects the finite groove geometry and off-corner measurement position.
Stationary focus: local-field anisotropy and oriented defects
With the beam held fixed, the relevant field is the quasi-static near-field of the induced dipole, Eq. (5). Its magnitude along the polarisation axis (the poles, ) and across it (the equator, ), at the same radius, is
so the field is twice as large along as across it and the deposited dose () obeys
This four-to-one anisotropy is exact for a point dipole[50] and is active during the inscription of a single modification: a colour centre, dangling bond, or amorphised voxel even with no scanning (Figure 2). The energy deposition, and hence the envelope within which the defect is generated, is elongated along . Two distinct effects should not be conflated here: (i) the deposition (damage-probability) envelope around the seed is elongated along by the 4:1 near-field anisotropy of Eq. (9), a statement about where a defect is likely to form, not about the shape of a point defect and (ii) the microscopic orientation of the created defect can follow the field direction through bond-aligned strong-field ionisation (§ 3). The h-BN experiments provide evidence for deterministic, polarisation-correlated defect creation at the quantum-confinement scale[13], while the deterministic polarisation control of ablation morphology was demonstrated on metal surfaces[58]. Once the modification extends into a groove, the complementary radiation/boundary-condition picture below takes over.
Near-field, not radiation, governs the local deposition
A clarification is essential. The far-field radiation of an oscillating dipole follows the toroidal pattern
with measured from , maximal perpendicular to and vanishing along it. At the nanometric distances that matter for writing (), however, the radiative term is negligible: the field is dominated by the quasi-static near-field [, Eq. (5)] and by the electrostatic boundary conditions at the seed/groove interface. The local deposition is therefore set by the near-field and boundary conditions, consistent with the seed being a sub-wavelength, non-radiating absorber, not by the radiation torus. Equation (10) is retained only to note that the far-field anisotropy and the near-field boundary enhancement share the same symmetry: both place the strongest transverse response across .
Boundary-condition enhancement inside the groove
The electrostatic boundary conditions point the same way and supply the quantitative enhancement. At the wall of an ablated groove, the normal component of is continuous, , so a field normal to the wall is enhanced inside the low-permittivity groove by[40]:
For an air groove in a titania film (, ) the internal field is enhanced -fold and the intensity -fold for the normal component, an order-of-magnitude electrostatic estimate for an idealised planar gap. For the nano-pen groove of Figure 3(a), FDTD gives a peak intensity ratio of 36 across the groove wall (Figure 3(b,Figure 3c)). The actual enhancement at a finite-curvature, lossy, evolving groove wall will differ. A field tangential to the wall is continuous and not enhanced. This normal-field boundary enhancement is the near-field mechanism we associate with the nano-pen[40], placing the strongest deposition transverse to .
The writing rule: polarisation perpendicular to the scan
A groove runs along the scan direction, so its long walls are parallel to the scan. To enhance the field inside the groove the polarisation must be normal to those walls, that is, perpendicular to the scan (Figure 3). This is exactly the demonstrated nano-pen: a nm groove ablated in a 20 nm titania film on silica, its width fixed by the normal- enhancement of Eq. (11) and the seed size rather than by the wavelength[40] [Fig. 3]. The enhanced normal field inside the groove deposits energy at the advancing wall and propagates the groove along the scan. With the polarisation instead parallel to the scan, the field at the long walls is tangential, continuous and not enhanced, and no narrow groove forms. The polarisation thus selects the writing direction: the nano-pen draws its finest line with perpendicular to the scan, and the same rule governs the seed-merging of the nano-knife[14, 40]. Rotating the polarisation by 90° switches the narrow-groove writing off, the same polarisation control was used to suppress periodic ripples on ablated Ni[58]. For completeness, we note that polarisation-independent nanostructuring is also possible when the deposition is dominated by thermal stress rather than by the near-field, as recently shown for microsphere-assisted femtosecond structuring[28]; the contrast between the two regimes underlines that the polarisation signature is a fingerprint of the near-field mechanism.
From ripples to a single groove: tuning the energy deposition. That the deposition can be continuously tuned between the two canonical surface morphologies, the periodic ripple (LIPSS) pattern and a single sub-diffraction groove, was demonstrated early on sapphire[39]. Ripples with period form via plasma nano-planes nucleating in a sub-surface region (a sphere-to-plane growth of the ionised volume, seeded by accumulated defects) and are imprinted on the surface by ablation. As the exposure is brought to within of the ablation threshold by reducing pulse energy or pulse-to-pulse overlap, the multi-groove pattern collapses to a single 50–70 nm groove (, one tenth of the nm diffraction-limited spot at the NA = 0.7, 400 nm conditions used), written by a controlled avalanche breakdown. This is the same near-threshold logic as Eq. (1), with the decisive practical advantage that the ripple-mediated route needs only pulse-energy stability, whereas confining a bare Gaussian to a 50 nm contour would demand . Polarisation enters here too: ripples parallel and perpendicular to the scan exhibit different formation thresholds and slightly different periods, attributed to the anisotropic heat transport and defect accumulation along the written texture[39], an early quantitative instance of the scan/polarisation coupling formalised by the writing rule above, and a direct precedent for the seed-and-merge writing of the nano-pen and nano-knife[14, 40].
The normal-field (out-of-plane) case and the single-bond limit
The third distinct geometry is polarisation normal to the surface. At high-NA the non-paraxial focus already carries a longitudinal , and a seed at an interface enhances the normal field by normal- continuity. Driving the field out of the plane opens the resonant-absorption channel and, in the extreme, a charged surface layer can undergo Coulomb explosion, removing material a monolayer at a time according to the theory of ablation in a field normal to the surface via the layer-by-layer ion explosion in the sub-ion-period regime[42]. This is the route to atomic-precision surface removal that has no diffraction-limited analogue. The experimental verification of this mechanism is still expected. The polarisation anisotropy persists all the way to the single bond. The quiver amplitude of a free electron in the optical field is . At () and (),
which exceeds the 0.16 nm Si-O bond length in silica[22]. This vacuum quiver amplitude is only a field-scale indicator, not a proof of bond removal: the displacement of a bound electron is reduced by local-field and screening factors and depends on the Keldysh parameter, collision time and band structure, while the ions barely move on the pulse timescale ( makes the 1030 nm field-scale × that at 800 nm). What it does show is that the optical field becomes comparable to the intra-bond field, so polarisation-sensitive strong-field (tunnelling) ionisation[16] can preferentially break the bond aligned with and inscribe an oriented defect, a plausible and experimentally indicated path to polarisation-controlled colour centres and single dopants[13, 41, 43, 44]. Demonstrating true single-bond determinism requires tunneling-rate/local-field modeling or atomistic simulation, not the quiver estimate alone.
3D field concentration at a reduced-permittivity core
The same dipolar localisation operates in the bulk, and a generic argument shows why the field pins itself to a region of depressed permittivity, whether that region is the core of a Bessel filament or the voxel of an ordinary Gaussian focus (Figure 4). In a charge-free medium , so a field crossing a permittivity gradient acquires the divergence

Figure 4. Field concentration at a Gauss-Bessel (GB) focus inside the bulk (as in ref.14); the same arguments are valid for a Gaussian beam. (a) The diffraction-limited GB focus () ionises a reduced-permittivity core; by Eqs. (13)–(15) the field is pulled into that core and, because only the component along is amplified, the near-field is anisotropic (lobes along the polarisation , transverse to the propagation ). (b) Transverse intensity cut: the deposited dose (red, “ENZ core”, FWHM ) is far narrower than the optical focal envelope (blue dashed, ) and is steered by for GB or Gaussian beam. This is a near-field effect localised to the core.
and the wave equation correspondingly carries the gradient term [14]; the component of along (normal to the seed boundary) is amplified where is small. For a field obliquely crossing a graded layer this is exact, the in-plane component obeys
with ( the host refractive index, the angle of incidence) and a slowly varying envelope that stays finite at the minimum. So, the field is pinned to the permittivity minimum and diverges at the ENZ point [14]. Equivalently, a sub-wavelength inclusion of reduced permittivity carries a uniform, enhanced internal field[50, 51],
rising to and respectively at ENZ, so the absorbed power self-focuses at the core. This is the bulk counterpart of the surface boundary-condition enhancement of Eq. (11); identical physics whether the low- region is a filament (cylinder) or an ordinary focal voxel (sphere). The concentration is intrinsically polarisation-steered, exactly as for the stationary near-field [Eq. (9)]: only the component along feels the amplification of Eq. (14), so even for an ordinary Gaussian focus (Figure 4) the core’s near-field and the next increment of deposition is strongest at the poles along the polarisation, the inscribed voxel is correspondingly anisotropic, and rotating rotates its axis. The in-volume permittivity contrast is gentler than that at an open surface (the host surrounds the seed on all sides rather than air on one), so the bulk anisotropy is weaker but of identical dipolar origin, the mechanism by which the cut plane and width track the rewriting polarisation in super-stealth dicing[14].
Two cautions complete the picture, and both follow from this being a near-field ENZ effect. First, the amplification [Eq. (14)] and the localisation [Eq. (6)] are quasi-static: they live within of the seed and do not radiate, they fix where energy is deposited, not how light propagates in the far field. Second, the responsible region has reduced permittivity, i.e., a negative index change, which is anti-guiding and cannot support a bound propagating mode. This mechanism, therefore, does not by itself create a propagating waveguide. A guiding (type-I) waveguide requires the opposite, a permanent positive index change from post-pulse densification or stress, a distinct material outcome from the transient ENZ energy deposition analysed here.
Real-life engineering: direct write that already works
The localization physics above translates into engineering (Figure 5) whose common denominator is that direct write removes constraints that parallel lithography cannot.

Figure 5. Practical implementation of polarisation controlled cutting, (stealth-)dicing, “wafering” (tomo-slicing) by Gauss-Bessel (GB) beam14 and cylindrical focusing57; see supplement in ref.57 for tomo-slicing conditions and mechanism of near-field energy deposition. Microwave energy deposition can be additional enhancement of this method by pre-heating during slicing or heating after laser inscription as hypothesized in this perspective review (§ 11).
Stitching. A central obstacle in large-area lithography is the stitching error between adjacent write fields, compounded by the cost of electron beam lithography (EBL)-made masks. Direct laser writing dissolves the stitching problem in two distinct ways[4]. (i) In a moving-beam moving-stage architecture the pattern is written continuously, avoiding field boundaries; combined with real-time surface tracking it scales to masks. (ii) In femtosecond direct write the written voxel is not a binary exposure but a multi-property material state (refractive index, birefringence, stress, crystallinity, conductivity); stitch-free fabrication is therefore a state-continuity problem, addressed by multi-dimensional holography and polarisation-resolved (four-polarisation) diagnostics that close the fabrication loop[4, 55, 56]. Parallel holographic writing with spatial light modulators or metalens focal-spot arrays eliminates stitching by writing the whole pattern at once[59, 60].
Solar cells and large-area masks. Bessel-beam direct write of – alumina etch masks enables photonic-crystal light-trapping textures on Si, targeting above-Lambertian absorption[6, 61, 62] and world-record-class cell efficiencies[63], with intentional non-ablated support pillars allowing undercut-tolerant masks over panel-scale () areas written by polygon scanners at surface speeds approaching [64, 65].
Functional 3D devices. The same toolbox writes domain structures in LiNbO for nonlinear optics[66], conductive graphitic tracks and diffractive optics in diamond and wide-bandgap semiconductors[43, 67], hyper-doped Si regions () for room-temperature single-electron transistors[12], and high-pressure phases via confined microexplosion[23, 46, 47]. Polarisation control of ablation via the normal -field component[58] and in-situ THz emission diagnostics of the focal currents[54] make these processes both manufacturable and observable (in situ monitored and controlled).
Outlook Hypothesis: Laser-Seeded Microwave Volumetric Deposition
This closing section is deliberately forward-looking: it is the hypothesis-generating element expected of an review–perspective, it builds directly on the absorber-defined localisation that is the review’s central theme (the laser defines where while a second, longer-wavelength source supplies energy there). We subject it to a quantitative feasibility analysis below, which identifies heat conduction as a severe constraint rather than claiming a worked scheme. We now develop the central hypothesis of this perspective. The threshold effect and the dipole/ENZ mechanism localise where energy is deposited to a tens-of-nm seed. Independently, microwave heating is the canonical volumetric and selective energy source in materials processing: just as a kitchen microwave oven deposits energy by shaking water molecules throughout the entire volume of the food rather than at its surface, a cavity field dissipates power only in lossy regions while leaving low-loss hosts cold[68, 69, 70]. The deposition is therefore intrinsically volumetric with no skin-depth or surface limitation at these scales, which is exactly what lets it fill a buried nanoseed uniformly.
We propose to combine the two: a femtosecond pulse, focused to a wavelength-scale spot ( m) inside a transparent sample, ionises a nanoscale seed at the tip of the Gaussian intensity profile by the threshold effect (Figure 6). The sample sits inside a microwave cavity whose field is non-resonant and hence transparent to the pristine host, but couples strongly to the laser-defined lossy seed (or inscribed pattern). The microwave then supplies volumetric energy preferentially to the laser-defined structure and modified nano-volume (Figure 5, Figure 6). Whether this can sustain a localised hotspot against heat conduction is the central question, examined quantitatively and found restrictive as discussed next.

Figure 6. Concept of fs-laser inscription of microwave absorbing antennas and nano-micro-volumes in the close proximity, which can be annealed in a microwave cavity. (a) Examples of THz Bessel beam generator (concentric grating) and polariser grating inscribed in a Nitrogen-doped (100 ppm) diamond. The width of the graphitised plane is m. (b) Sairem microwave cavity setup (2.45 GHz, 200 W) for sample annealing. The inset shows the concept of graphitised antenna heaters with inscribed single-pulse structured nano-volume modifications for controlled annealing. Laser writing inside a microwave-heated sample is also feasible.
Microwave Coupling to the Laser-Defined Nanovolume
Once the seed exists, it carries free or defect-band carriers and hence a finite conductivity , equivalently a loss at the microwave frequency , while the host has a negligible loss tangent. Three length-scale separations make the coupling clean:
Quasi-static. – cm , so the cavity field is uniform across the seed.
Field penetration. The skin depth is macroscopic: for S/m and GHz, m , so the field fully penetrates and heats the seed uniformly.
Per-volume selectivity. The deposited power density is larger in the seed than in the host by –. This does not mean the seed absorbs most of the total power: integrated host and cavity-wall absorption scale with their (vastly larger) volumes, so selective total heating of one seed requires many seeds, a microwave near-field concentrator of sub-m mode volume, or strong suppression of host/cavity loss.
The internal field of the quasi-static sphere is , and the time-averaged Joule (volumetric) power density is
For a moderately conducting (impedance-matched) seed, and , giving the working expression
deposited uniformly throughout the seed volume . The localisation is set by (tens-of-nm) whereas the microwave wavelength is cm: the energy is deposited into a region times smaller than the wavelength—sub-wavelength deposition by more than six orders of magnitude.
Thermal runaway: sharpening toward single nanometers
The conductivity of the seed rises with temperature in semiconductor-like fashion, , with a prefactor, the activation energy, Boltzmann’s constant and the temperature. Substituting into Eq. (17) gives : hotter sub-regions absorb more, a positive feedback that is the well-known thermal-runaway/selective-absorption instability of microwave processing[68, 69]. Microwave high-temperature annealing was demonstrated in formation of nanoscale gradient P-doping of Si beyond the thermodynamical solubility[71].
Within the dipole dose profile, Eq. (6), the central maximum runs away first, contracting the actively heated zone, in principle toward smaller scales; but lattice heat diffusion (length , thermal diffusivity ) imposes a stringent floor that, as § 17 shows, dominates at these parameters. Whether higher frequency helps is condition-dependent: at fixed conductivity , is frequency-independent; at fixed , rises with , but defect-relaxation and hopping losses are dispersive and roll off at high frequency. A 24 GHz to 30 GHz gyrotron source[72] thus helps only under a fixed-, fixed-field assumption that must be checked against the loss model. Deliberately driven, this runaway is also an opportunity rather than a nuisance: a self-localising, microwave-sustained hotspot seeded by the laser could assist the dicing and wafering of hard dielectrics and crystals: sapphire, SiC, diamond and optical glasses where purely optical cutting is slow and the kerf hard to control, with the microwave supplying the bulk of the separation energy volumetrically along the laser-defined plane (Figure 5).
Seed loss from point defects
How large is the seed loss in practice, and what sets the conductivity used below? The microwave conductivity is , so the range – S/m corresponds to – at 2.45 GHz (where S/m). Pristine silica is essentially loss-free here (, , S/m), whereas the seed/host contrast can exceed . The laser-written core, by contrast, is saturated with point defects: oxygen vacancies, silicon dangling bonds, non-bridging-oxygen and oxygen-interstitial (peroxy) centres[73] at densities – cm. The upper decade corresponding to heavily exposed, near-breakdown volumes. The number densities of individual defect species in irradiated and as-drawn silica are catalogued in refs. 73–75 and the local density in a laser-written core depends strongly on exposure dose and annealing history. Two channels make these lossy at GHz. (i) Debye relaxation of defect dipoles contributes ; with a defect dipole nm this gives at cm, peaking in near . (ii) Thermally activated hopping/small-polaron conduction between defect sites is the mechanism shown to raise the GHz dielectric loss of oxygen-vacancy- and Ti-rich titania-silica glass by orders of magnitude[76], and to make oxygen-deficient oxides strong microwave absorbers. Together these place – (hence – S/m) within reach for a defect-saturated seed, rising further if residual free carriers survive the microwave duration. A direct measurement of for laser-written volumes (Figure 6) is the key input needed for precise estimates. It also provides a real-time feedback of material modifications during annealing.
A practical measurement strategy for . The measurement is feasible with standard microwave metrology, provided the small modified-volume fraction is respected. A high- cavity-perturbation or split-post dielectric-resonator measurement at 1–10 GHz compares a reference wafer with an identical wafer carrying a laser-written grating of known volume fraction (e.g. – written planes, –): the frequency shift and degradation yield the composite and , from which the intrinsic follows by effective-medium inversion. Because , resolving at requires , demanding but routine in sapphire-loaded or superconducting cavities. The defect density in the same samples is quantified independently by electron spin resonance (EPR) and optical absorption of the known silica centres[74, 75, 77]. The written volumes can be re-measured after stepwise furnace or microwave annealing (Figure 6), giving versus on one sample set, and simultaneously testing the real-time-feedback idea above. The main experimental challenges are separating surface/adsorbate losses from bulk ones (addressed by measuring before/after writing on the same wafer) and the temperature dependence of the host loss during annealing runs. Feasibility of defect heating by microwaves is evaluated in § 17.
Conclusions and Outlook
Three messages run through this review. First, the highest-resolution energy deposition is achieved not by shrinking the wavelength but by localising the absorber. The threshold effect [Figure 1(b)], the interface boundary conditions, and the dipole/ENZ near-field with its tail (Figure 2 and Figure 4) confine the deposited energy to a tens-of-nm seed whose size is set by the laser nonlinearity and the seed itself, not by . Near-threshold exposure reaches the tens-of-nanometre scale and the nm modification scale when combined with the seed near-field mechanisms of § 5. The final feature size further depends on the post-pulse material response: heat diffusion and quench rate, Eq. (2), which favours low-diffusivity hosts and confinement geometries. This near-field route is intrinsically wavelength-insensitive, a conclusion that the latest surveys of ultrafast-laser nanofabrication[26, 27] and of resonant ENZ light-matter interaction[53] independently reinforce, and that holds whether the low-permittivity seed is created electronically (ENZ/Die-Met transit) or mechanically (cavitation/rarefaction, § 3 B).

Figure 1. Breaking the limit in aspect ratio of 3D modification (adapted from ref.14). (a) Equivalence of the Heisenberg’s uncertainty principle and the waist-divergence relation for the Gaussian pulse (standard notations: , , ). The aspect ratio DoF/waist of the laser focus is limited; the same is valid for a Bessel-beam axially long pulse, limited at a higher ratio. (b) The threshold effect: the damage/modification region defined by a threshold fluence reduces the transverse diameter and the depth of the laser-modified volume simultaneously. Therefore the final aspect ratio of laser-modified volume is limited by the shape of laser focus . (c) Side-view SEM image of laser cut (long focus was scanned with overlap in out-of-page direction). Nano-seed structural damage sites of silica glass are merging into a nano-line via near-field enhancement, which is still localised at nanoscale. (d) Breaking the limit: merged nano-seeds into nanoline cuts crystal at the record aspect ratio .
Second, polarisation controls the directionality and anisotropy of the deposition. The quasi-static near-field and the normal-/ boundary enhancement fix the writing direction, the finest lines are drawn with perpendicular to the scan (Figure 3) and the orientation of the inscribed modification down to the single-defect scale (Figure 2). Together, the threshold effect and polarisation-steered near-field deposition constitute a maskless 2D and 3D direct-write nanolithography operating at the nm scale, already practised as the nano-pen on surfaces and the nano-knife in the bulk [Figure 1(c,Figure 1d) and Figure 5]. To keep claims and goals distinct: 10–20 nm features, the polarisation-perpendicular writing rule, record-aspect-ratio dicing and polarisation-correlated single-emitter creation in h-BN are experimental results[13, 14, 40]. Deterministic single-defect placement on demand, single-bond selectivity and layer-by-layer electrostatic surface removal remain research goals, and the mechanism attribution of the writing rule awaits dedicated polarisation-angle-resolved studies (Sec. VI).
These principles are not merely theoretical. A clean, benign read-out of the ENZ energy-deposition physics is the multiphoton-polymerisation voxel, whose abrupt lateral-growth acceleration at the diffraction limit we analyse in Appendix § 16: below the focal-spot size the Gaussian threshold sets the width, while above it the ionised core, plausibly, though not yet provably, an ENZ-class volumetric absorber takes over and the two voxel boundaries widen independently. The same maskless, 3D-capable physics carries an economic message that is increasingly strategic: as set out in Appendix § 19, a M femtosecond writer can jump-start sub-10 nm and 3D-photonic prototyping in any economy without access to the M-class EUV steppers held by a handful of fabs, a route to sovereign nanofabrication capability that matters acutely amid today’s fragile, interruption-prone supply chains.
Third, and finally, the laser-seeded microwave hypothesis extends the central logic to its conclusion. The microwave is a volumetric, molecular-scale energy-delivery source: it couples to the very dipoles, defects and free carriers that the laser creates, and so realises guided energy deposition at the nano-/molecular scale. The laser defines where (a tens-of-nm, ENZ-defined volume) and the cavity supplies energy preferentially there (per unit volume), at a scale more than six orders of magnitude below its own wavelength (Figure 5). We present this as a hypothesis, not a worked scheme: the feasibility analysis of § 17 shows that a single isolated seed cannot be held hot ( mK at the assumed parameters, worsening as ), so any realisation must rely on field concentration, dense seed ensembles, sub-nanosecond pulsing, or a threshold-bias role rather than primary continuous heating. Future work must quantify the seed conductivity and its lifetime, the cavity-sample coupling, and the balance of microwave power against conductive loss; a concrete measurement path for is given in § 11 C. For fs-laser-inscribed macroscopic patterns of defects and seed-absorber antennas, the usual antenna arguments of coupling and microwave-field alignment would guide the localised heating inside the volume of the sample at the closest micro-nano proximity to the antenna. Such laser-inscribed microwave antennas would provide controlled thermal annealing of the neighbouring nano-/micro-scale structures and defects, from within the volume, beyond what currently available flash-annealing[78] approaches deliver from the outside.
Acknowledgments
The authors acknowledge discussions with Eugene G. Gamaly, Vladimir T. Tikhonchuk, Hong-Bo Sun, Zhenze Lee, Lei Wang, Sajeev John, Leonid V. Zhigilei. Support from the ARC and the Vilnius University Laser Research Center is gratefully acknowledged. S.J. is grateful to Gintas Šlekys and Workshop of Photonics Ltd., Lithuania for technology transfer project and first high-power femtosecond fab station at Swinburne Uni. This research was funded by ARC DP240103231 and Swinburne’s Proof-of-Concept “Silicon Carbide Single Photo Emitter Fabrication” grants.
AUTHOR DECLARATIONS
The authors have no conflicts to disclose.
DATA AVAILABILITY
Data sharing is not applicable to this review article as no new data were created or analyzed.
ENZ onset at the focus during multiphoton polymerisation
Multiphoton polymerisation (MPP) is a clean model system in which the transient ENZ/Die-Met onset can be read directly from the written voxel: a polymerisable resin is transparent and low-loss until the focal ionisation drives its instantaneous permittivity down. Using a broadly tunable source (700 nm to 1300 nm) with pulse-duration and dispersion control, the lateral and longitudinal voxel growth and the dynamic fabrication window (DFW, the damage-to-polymerisation threshold ratio) across the two- and three-photon bands was measured[36]. Two observations bear directly on the present picture: (i) the voxel growth shows an abrupt change in its power (fluence) dependence; and (ii) the DFW shrinks strongly toward longer wavelength, from ≈ 29 at 700 nm to ≈1 at 1250 nm. Both are interpreted as a consequence of a change in the instantaneous refractive index (plasma/ENZ formation) that lowers the polymerisation, but not the damage threshold[36].
This can be sharpened with the dipole/ENZ mechanism of the main text. Below the diffraction-limited focal-spot size, the modified lateral radius follows the threshold-on-Gaussian law of Eq. (1), : the voxel samples only the Gaussian tail above threshold, so it widens logarithmically, i.e. slowly and traces the beam profile at the true nanoscale. Once approaches , however, the ionised core fills the high-intensity focal volume and the Drude free carriers turn it into a volumetric absorber: the instantaneous permittivity falls toward ENZ, the internal field is enhanced [Eq. (15)]. The absorbed power now grows with the carrier density the pulse itself creates. The lateral size therefore accelerates with pulse energy rather than saturating as a passive Gaussian-threshold profile would. The abrupt change in the power dependence is seen in MPP. In short, below the beam sets the size; above the absorber does.
Status of this interpretation. We emphasise, as a matter of scientific hygiene, that a change in the slope of feature size versus pulse energy does not by itself establish a transient ENZ state. It establishes a change of growth regime. The ENZ/volumetric-absorber reading is a possible explanation, favoured here because it follows from the same equations [Eqs. (1), (15)] that account for the nano-pen and nano-knife, but at least three alternatives must be excluded by data: (i) radical/heat diffusion beyond the exposed volume (proximity effect), which also accelerates widening at high dose; (ii) cumulative thermal effects at high repetition rate; and (iii) focal aberrations that reshape the intensity distribution at high energy.
The wavelength trend follows from the same physics. Free-carrier (Drude) absorption scales as , while the critical density at which scales as . both make the ionised volume absorb more strongly, and reach ENZ at a lower carrier density, at longer wavelength. A focus at 1100 nm therefore couples energy into its ionised core more efficiently than one at 800 nm ( larger by ), giving a steeper lateral growth and a correspondingly narrower fabrication window; the DFW collapse toward longer wavelength is reported[36]. The polymerisation voxel is thus a benign, direct read-out of the same ENZ energy-deposition physics that, in harder materials, drives modifications discussed as the nano-pen (nanolithography)[40] and nano-knife (dicing and tomo-slicing or wafering)[14, 57].
Independent growth of the voxel boundaries at the diffraction limit
A complementary and, we argue, likely contribution to the accelerated lateral widening of Appendix A1 is purely geometric. While the modified region is smaller than the focal spot, its two transverse sides lie within one near-field range of each other (the dipole field decays over ), so they act as a single, coupled low- core with one effective deposition centre. Once the voxel grows to the diffraction-limited focal-spot size, however, the two lateral boundaries become separated by more than that near-field length and decouple: each is now an independent permittivity step that carries its own normal- enhancement [Eqs. (11) and (14)], exactly as the two walls of a 2D groove do in the nano-pen. The two fronts then deposit energy and advance independently of the voxel centre, so the polymerised line widens from both edges at once. This doubling and decoupling of the lateral growth channels sets in precisely at the focal-spot size and provides a natural, geometric reason why the side widening accelerates there rather than saturating, consistent with the abrupt change in the lateral voxel dynamics[36]. A falsifiable corollary is that, deep in this regime, the line cross-section should evolve from the single Gaussian-like core found below the diffraction limit toward a flat-bottomed or twin-lobed profile bounded by the two outward-moving edges.
Feasibility estimates
Table II collects representative power densities. Taking a cavity field V/m (10 kV/cm, attainable from a modest W source in a high- (cavity quality factor ), field-concentrating (re-entrant) cavity), a seed of radius nm ( m) and –100 S/m, Eq. (17) gives – W m and a deposited power per seed – W.
| (S/m) | (W m) | (W) | |
| 1 | 20 ms | ||
| 10 | 2 ms | ||
| 100 | 0.2 ms |
Table II. Representative power densities for laser-seeded microwave deposition. Seed radius nm, V/m, J/m. The times are the adiabatic (no-conduction) upper bound; the conduction loss W at K greatly exceeds (see text).
Heat transport is the binding constraint in its simplest form argues against single-seed heating. Indeed, a 25 nm seed is far from adiabatic: with thermal diffusivity m/s, the diffusion time is ns, so heat leaves the seed in nanoseconds, and over a 0.2 ms to 20 ms microwave dwell it spreads –140 μm, thousands of seed radii. In steady state the temperature rise at the surface of a continuously heated sphere is (the centre is 1.5 times hotter for equal conductivities): for S/m and nm this is only mK (silica thermal conductivity W/m/K), and the size scaling is unfavourable. The absorbed power falls as while conductive loss falls only as , so and smaller seeds heat less. Equivalently, the steady-state power needed to hold one seed at K is W, five to seven orders of magnitude above the – W the cavity delivers per seed. Reaching 1000 K in a single 25 nm seed with S/m would require V/m, a field enhancement of over the assumed V/m. For a dense ensemble that fills a region of radius with volume fraction , the same estimate gives µm for and µm for at K (for S/m), which is the length scale of the laser-inscribed planes and antennas of Figure 5 and Figure 6. A single isolated seed therefore cannot be held hot by the microwave at these parameters. The adiabatic times in Table II (with J/m) are upper bounds that ignore conduction.
The proposal can close only if one or more of the following holds: (i) the microwave acts as a transient bias that lowers a threshold the laser already nearly meets, rather than as the primary heat source; (ii) field concentration raises (hence ) by several orders; (iii) a dense ensemble of seeds is heated collectively, so that the warm region has size and the per-seed conductive loss falls by . For the laser-inscribed planes of Figure 5 with seed spacing this is the realistic operating mode, and macroscopic graphitised antennas (Figure 6) are its deliberately engineered limit; or (iv) heating is pulsed on sub- timescales. We therefore present the microwave route as a hypothesis whose single-seed energetics are unfavourable, whose ensemble/antenna energetics are plausible, and whose quantitative inputs (, seed lifetime, cavity coupling) remain to be measured.
Experimental implementation and diagnostics
The scheme requires (i) a high-NA femtosecond writer focusing inside a sample placed at an antinode of a microwave cavity (rectangular or cylindrical mode, GHz magnetron or 24 GHz to 30 GHz gyrotron), and (ii) synchronisation so that the cavity field is present during and after each pulse. Diagnostics already developed for USPL–matter interaction transfer directly: in-situ THz emission, which tracks currents and electron temperature in the focal volume[54], four-polarisation imaging of near-field anisotropy below the diffraction limit[55, 56], which directly visualises the polarisation-selected deposition of Sec. VI, and plasma/optical breakdown spectroscopy. Polymerisable resins without photoinitiators provide a clean model system in which the ENZ onset is visible as a step in the polymerised cross-section[36], allowing the laser-only and laser-plus-microwave regimes to be compared quantitatively.
Sovereign capability and the economics of maskless 3D direct write
The scientific case above has an economic counterpart that may be decisive for countries without access to extreme-ultraviolet (EUV) infrastructure. A High-NA EUV stepper now costs million per tool (low-NA EUV million)[79], requires vacuum, reflective masks and a -step planar flow, and is available to a mere handful of fabs. A femtosecond direct-write station, by contrast, is a million instrument, two to three orders of magnitude cheaper and, being maskless and 3D-capable, can serve four distinct roles from one platform: (i) writing the projection-lithography masks themselves, replacing slow electron-beam mask-making (the nano-pen)[4, 40]; (ii) inscribing 3D photonic chips - buried waveguides and quantum emitters/colour centres in the bulk of glasses and crystals[34, 41, 43]; (iii) on-chip photonic wiring by two-photon polymerisation[32, 59]; and (iv) sculpting 3D micro-optical elements: lenses, Fresnel and diffractive optics, phase plates - out of glasses and crystals[66, 67]. This is the “3D-printing-like” democratisation that lets a laboratory in Australia, Lithuania or any non-EUV economy prototype sub-10 nm and 3D devices at risk capital, in line with sovereign-capability and supply-chain-resilience goals[80].
A scaling argument shows why the maskless route wins precisely where modern devices are heading: 3D and photonic. A 2D projection stepper images one layer per exposure, so an intrinsically three-dimensional structure of functional layers must be emulated by maskexposeetchdeposit cycles: cost and time grow linearly in on top of a per-mask charge of order . Direct write instead lays the structure down volumetrically in a single pass, its throughput set not by layer count but by the optical Moore’s law: pulse repetition rate times the number of parallel beams. At 1 MHz to 1000 MHz rates with beams and polygon scanning at surface speeds approaching km/s[64, 65], – voxels per second are written, so a -voxel 3D photonic chip is completed in seconds to minutes with no masks at all. Three caveats keep this comparison honest. First, the quantities compared are of different kinds: tool prices are vendor specifications, the km/s scan speeds and multi-beam counts are demonstrated engineering performance on specific tasks (large-area ablation of thin etch-masks[61, 64, 65], parallel polymerisation[59, 60]), while the voxel-per-second figures are idealised projections that assume every pulse writes a useful voxel. Second, a voxel rate is an upper bound on exposure throughput, not an estimate of device completion time: stage overheads, trajectory planning, in-situ metrology, development/etch steps where applicable, and, above all, yield qualification are excluded, and these dominate real fabrication schedules exactly as mask qualification dominates EUV insertion. Third, the comparison is task-specific: the materials and tasks underlying the direct-write figures are glasses, polymers and dielectric films patterned as photonic/mask structures, not transistor-grade silicon. There is therefore a complexity cross-over. Below a few binary 2D layers, projection lithography is unbeatable for volume CMOS. Once becomes large and the geometry is genuinely three-dimensional, however, the sequential 3D writer can become faster and substantially cheaper. This regime includes photonic integrated circuits, whose market is projected to grow from approximately billion in 2024 to billion by 2030[81], and which industry roadmaps increasingly pursue through 3D heterogeneous integration and co-packaged optics[82].
The division of labour between the two technologies is then clear and should be stated plainly. High-volume planar logic and memory remain the domain of EUV projection lithography: at wafers per hour and transistors per die, no sequential tool competes on throughput, and the mature yield-learning ecosystem is as decisive as the optics. Direct write offers the greatest value where EUV’s strengths are irrelevant or its economics prohibitive: intrinsically 3D photonic devices (waveguide circuits, 3D photonic crystals), quantum-device prototyping (colour centres, single-emitter arrays, one-off chip variants where a mask set is unjustifiable), photomask writing and repair itself, low-volume/high-mix component manufacturing (micro-optics, diffractive elements, dicing/wafering of hard crystals), and large-area unique patterns (solar-cell light-trapping masks) where stitching, not pitch, is the binding constraint. Throughput, yield and device geometry decide the split: parallel projection wins when millions of identical planar layers amortise mask and qualification costs. Sequential direct write wins when the geometry is 3D, the volume is low, or the pattern is unique.
Why, finally, do the two Moore’s laws hold at all? Neither is a law of physics; both are learning curves. Empirically, cost falls as a power of cumulative production (Wright’s law) while production rises exponentially, and the two combine into the exponential-in-time improvement we call Moore’s law, a regularity demonstrated across dozens of technologies and shown to be genuinely forecastable[83]. The e-Moore engine has ridden this self-reinforcing investment-capability loop through parallel 2D scaling; the o-Moore engine rides the same loop in ultra-short-pulse laser power[2]. The thesis of this review is that, as value migrates from planar transistor count to 3D photonic and quantum function, the o-Moore, sequential direct-write engine is poised to convert its exponential power growth into exponential fabrication throughput at a fraction of the capital cost.
References
- [1]G. E. Moore, Electronics 38, 114 (1965).
- [2]M. Han, D. Smith, S.-H. Ng, V. Anand, T. Katkus, and S. Juodk azis, in Proceedings of the 2nd International Electronic Conference on Applied Sciences (MDPI, Basel, Switzerland, 2021).
- [3]A. Dubietis and A. Matijosius, Opto-Electronic Advances 6, 220046 (2023).
- [4]H.-H. Huang, H. Mu, E. Puig Vilardell, V. Anand, D. Gailevicius, and S. Juodkazis, Nanomaterials 16, 692 (2026).
- [5]T. Itani and T. Kozawa, Japanese Journal of Applied Physics 52, 010002 (2013).
- [6]J. Maksimovic, J.-W. Hu, S.-H. Ng, T. Katkus, G. Seniutinas, T. Pinedo Rivera, M. Stauber, Y. Nishijima, S. John, and S. Juodk azis, Opto-Electronic Advances 5, 210086 (2022).
- [7]J.-P. Kleinhans and N. Baisakova, “The global semiconductor value chain: A technology primer for policy makers,” Tech. Rep. (Stiftung Neue Verantwortung / interface, Berlin, 2025).
- [8]European Court of Auditors, “The EU’s strategy for microchips: Special report on semiconductors,” Tech. Rep. (European Court of Auditors, Luxembourg, 2023).
- [9]ASML Holding N.V., “Annual report 2024,” https://www.asml.com/en/investors/annual-report (2025).
- [10]EUV Lithography Working Group, “EUV lithography ecosystem readiness report,” Tech. Rep. (SEMI / SPIE Advanced Lithography Working Group, 2023).
- [11]H. J. Levinson, Japanese Journal of Applied Physics 61, SD0803 (2022).
- [12]M. Tabe, D. Moraru, M. Ligowski, M. Anwar, R. Jablonski, Y. Ono, and T. Mizuno, Phys. Rev. Lett. 105, 016803 (2010).
- [13]X.-J. Wang, H.-H. Fang, Z.-Z. Li, D. Wang, and H.-B. Sun, Light: Science & Applications 13, 6 (2024).
- [14]Z.-Z. Li, H. Fan, L. Wang, X. Zhang, X.-J. Zhao, Y.-H. Yu, Y.-S. Xu, Y. Wang, X.-J. Wang, S. Juodkazis, Q.-D. Chen, and H.-B. Sun, Nature Photonics 18, 799 (2024).
- [15]N. Ganguly, P. Sopenã, and D. Grojo, Light: Advanced Manufacturing 5, 22 (2024).
- [16]L. V. Keldysh, Sov. Phys. JETP 20, 1307 (1965).
- [17]N. Bloembergen, IEEE Journal of Quantum Electronics 10, 375 (1974).
- [18]M. Garcia-Lechuga, O. Utéza, N. Sanner, and D. Grojo, Optics Letters 45, 952 (2020).
- [19]X. Liu, D. Du, and G. Mourou, IEEE Journal of Quantum Electronics 33, 1706 (1997).
- [20]A. P. Joglekar, H. Liu, E. Meyhöfer, G. Mourou, and A. J. Hunt, Proceedings of the National Academy of Sciences 101, 5856 (2004).
- [21]K. Kumar, K. C. Lee, J. Li, J. Nogami, N. P. Kherani, and P. R. Herman, Light: Science & Applications 3, e157 (2014).
- [22]E. G. Gamaly and A. V. Rode, Applied Physics A 124, 278 (2018).
- [23]S. Juodkazis, K. Nishimura, S. Tanaka, H. Misawa, E. G. Gamaly, B. Luther-Davies, L. Hallo, P. Nicolai, and V. T. Tikhonchuk, Phys. Rev. Lett. 96, 166101 (2006).
- [24]M. K. Bhuyan, M. Somayaji, A. Mermillod-Blondin, F. Bourquard, J.-P. Colombier, and R. Stoian, Optica 4, 951 (2017).
- [25]Z. Lin, H. Liu, L. Ji, W. Lin, and M. Hong, Nano Letters 20, 4947 (2020).
- [26]Z. Lin, L. Ji, and M. Hong, Engineering (2025), 10.1016/j.eng.2025.03.017.
- [27]R. Stoian, Advanced Optical Technologies 12, 1237524 (2023).
- [28]J. Yin, H. Luo, T. Cao, and M. Hong, Light: Science & Applications 14, 2091 (2025).
- [29]G. Zhang, A. Rudenkо, R. Stoian, and G. Cheng, Ultrafast Science 5, 0103 (2025).
- [30]Y. Guan, Y. Ding, Y. Fang, J. Li, Y. Liu, R. Wang, J. Hao, H. Xie, C. Xu, L. Zhen, et al., Nature Communications 16, 4149 (2025).
- [31]B. Fan, Y. Wang, J. Hu, X. Han, L. Zhong, and J. Qiu, Light: Advanced Manufacturing 7, 336 (2026).
- [32]K. Kawata, H.-B. Sun, T. Tanaka, and K. Takada, Nature 412, 697 (2001).
- [33]R. R. Gattass and E. Mazur, Nature Photonics 2, 219 (2008).
- [34]K. M. Davis, K. Miura, N. Sugimoto, and K. Hirao, Optics Letters 21, 1729 (1996).
- [35]E. N. Glezer, M. Milosavljevic, L. Huang, R. J. Finlay, T.-H. Her, J. P. Callan, and E. Mazur, Optics Letters 21, 2023 (1996).
- [36]D. Samsonas, E. Sklutas, A. Čiburys, L. Kontenis, D. Gailevičius, J. Berzinš, D. Narbutis, V. Jukna, M. Vengris, S. Juodkazis, and M. Malinauskas, Nanophotonics 12, 1537 (2023).
- [37]S. W. Hell and J. Wichmann, Optics Letters 19, 780 (1994).
- [38]J. Fischer and M. Wegener, Laser & Photonics Reviews 7, 22 (2013).
- [39]R. Buividas, S. Rekštytė, M. Malinauskas, and S. Juodkazis, Optical Materials Express 3, 1674 (2013).
- [40]Z.-Z. Li, L. Wang, H. Fan, Y.-H. Yu, Q.-D. Chen, S. Juodkazis, and H.-B. Sun, Light: Science & Applications 9, 41 (2020).
- [41]S. Castelletto, A. F. M. Almutairi, K. Kumagai, T. Katkus, Y. Hayasaki, B. C. Johnson, and S. Juodkazis, Opt. Lett. 43, 6077 (2018).
- [42]E. G. Gamaly and S. Juodkazis, Nanomaterials 13, 1796 (2023).
- [43]S. Castelletto, J. Maksimovic, T. Katkus, T. Ohshima, B. C. Johnson, and S. Juodkazis, Nanomaterials 11, 72 (2021).
- [44]S. Choi, T.-T. Tran, C. Elbadawi, L. Lobo, X. Wang, S. Juod- kazis, G. Senjutin̄as, M. Toth, and I. Aharonovich, ACS Applied Materials & Interfaces 8, 29642 (2016).
- [45]L. Rapp, T. Matsuoka, K. L. Firestein, D. Sagae, H. Habara, K. Mukai, K. A. Tanaka, E. G. Gamaly, R. Kodama, Y. Seto, T. Shobu, A. Tominaga, L. Smillie, T. Pikuz, B. Haberl, T. Yabuuchi, T. Togashi, Y. Inubushi, M. Yabashi, S. Juodkazis, D. V. Golberg, A. V. Rode, and N. Ozaki, Phys. Rev. Research 6, 023101 (2024).
- [46]L. Rapp, B. Haberl, C. J. Pickard, J. E. Bradby, E. G. Gamaly, J. S. Williams, and A. V. Rode, Nature Communications 6, 7555 (2015).
- [47]A. Vailionis, E. G. Gamaly, V. Mizeikis, W. Yang, A. V. Rode, and S. Juodkazis, Nature Communications 2, 445 (2011).
- [48]K. Osvay, L. Stuhl, P. Varmazyar, T. Gilinḡer, Z. Elekes, A. Fenyvesi, P. K. Singh, S. Hussain, A. Börzsönyi, Z. Fülöp, T. Tajima, and G. Mourou, Scientific Reports 14, 25302 (2024).
- [49]K. Osvay, L. Stuhl, P. Varmazyar, T. Gilinḡer, Z. Elekes, A. Fenyvesi, P. K. Singh, S. Hussain, A. Börzsönyi, Z. Fülöp, T. Tajima, and G. Mourou, European Physical Journal Plus 139, 574 (2024).
- [50]J. D. Jackson, Classical Electrodynamics, 3rd ed. (Wiley, New York, 1999).
- [51]C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).
- [52]J. B. Khurgin, M. Clerici, and N. Kinsey, Laser & Photonics Reviews 15, 2000291 (2021).
- [53]P. Xie, W. Wang, and Y. Kivshar, Applied Physics Reviews 12, 023807 (2025).
- [54]H.-H. Huang, T. Nagashima, K. Kumagai, Y. Hayasaki, S. Juodkazis, V. T. Tikhonchuk, and K. Hatanaka, Light: Advanced Manufacturing 5, 30 (2024).
- [55]M. Ryu, S.-H. Ng, M. Han, V. Anand, T. Katkus, J. Vongsvivat, D. Appadoo, Y. Nishijima, S. Juodkazis, and J. Morikawa, Nanoscale Horizons 7, 1047 (2022).
- [56]R. Honda, M. Ryu, M. Moritake, A. Balčytis, V. Mizeikis, J. Vongs̄vivut, M. J. Tobin, D. Appadoo, J.-L. Li, S.-H. Ng, S. Juodkazis, and J. Morikawa, Nanomaterials 9, 732 (2019).
- [57]Y. Chai, H.-H. Fang, Z.-Z. Li, T.-W. Wang, S.-F. Liu, H.-R. Chen, S.-C. Li, X.-Y. Li, J.-M. Lyu, and H.-B. Sun, Light: Science & Applications 15, 297 (2026).
- [58]H. Iwase, S. Kokubo, S. Juodkazis, and H. Misawa, Opt. Express 17, 4388 (2009).
- [59]S. K. Saha, D. Wang, V. H. Nguyen, Y. Chang, J. S. Oakdale, and S.-C. Chen, Science 366, 105 (2019).
- [60]S. Gu, C. Mao, A. Guell Izard, S. Sadana, D. Terrel-Perez, M. Mettry-Yassa, W. Choi, W. Zhou, H. Yan, Z. Zhou, et al., Nature 648, 591 (2025).
- [61]T. Katkus, S.-H. Ng, H. Mu, N. H. A. Le, D. Stonyte, Z. Kha-jehsaeidimahabadi, G. Seniutinas, J. Baltrukonis, O. Ulcinas, M. Mikutis, V. Sabonis, Y. Nishijima, M. Rienacker, U. Roemer, J. Krueḡener, R. Peibst, S. John, and S. Juodkazis, Advanced Engineering Materials 26, 2400763 (2024).
- [62]M.-L. Hsieh, A. Kaiser, S. Bhattacharya, S. John, and S.-Y. Lin,
- [63]F. Haase, C. Hollemann, S. Schäfer, A. Merkle, M. Rienacker, J. Krueḡener, R. Brendel, and R. Peibst, Solar Energy Materials and Solar Cells 186, 184 (2018).
- [64]B. Henriques and A. F. Lasagni, Biomedical Materials & Devices 2, 566 (2024).
- [65]D. Moskal, J. Martan, and M. Honner, Micromachines 14, 1241 (2023).
- [66]S. Juodkazis, M. Sudzius, V. Mizeikis, H. Misawa, E. G. Gamaly, Y. Liu, O. A. Louchev, and K. Kitamura, Appl. Phys. Lett. 89, 062903 (2006).
- [67]P. Wu, X. Cao, W. Chu, Z. Chen, H. Yuan, R. Wang, S. Yao, S. Juodkazis, and W. Zhang, Opt. Lett. 48, 1379 (2023).
- [68]Y. V. Bykov, K. I. Rybakov, and V. E. Semenov, Journal of Physics D: Applied Physics 34, R55 (2001).
- [69]K. I. Rybakov, E. A. Olevsky, and E. V. Krikun, Journal of the American Ceramic Society 96, 1003 (2013).
- [70]P. Ravi and A. Nozariasbmarz, ACS Applied Materials & Interfaces 18, 29452 (2026).
- [71]C.-H. Tsai, C. P. Savant, M. J. Asadi, Y.-M. Lin, I. Santos, Y.-H. Hsu, J. Kowalski, L. Pelaz, W.-Y. Woon, C.-K. Lee, and J. C. M. Hwang, Appl. Phys. Lett. 121, 052103 (2022).
- [72]Y. V. Bykov, G. G. Denisov, A. G. Eremeev, and M. Y. Glyavin, IEEE Transactions on Microwave Theory and Techniques 66, 3162 (2018).
- [73]G. M. Lo Piccolo, M. Cannas, and S. Agnello, Materials 14, 7682 (2021).
- [74]L. Skuja, Journal of Non-Crystalline Solids 239, 16 (1998).
- [75]G. Pacchioni, L. Skuja, and D. L. Griscom, eds., Defects in SiO₂ and Related Dielectrics: Science and Technology (Springer, Dordrecht, 2000).
- [76]C. Zhang, S. Gao, M. Zhu, Z. Shao, L. Nie, H. Wang, Y. Jia, and B. Fu, Materials 18, 320 (2025).
- [77]X. W. Wang, R. Buivid̄as, F. Funabiki, P. R. Stoddart, H. Hosono, and S. Juodkazis, Applied Physics A 122, 194 (2016).
- [78]Y. Li, K. Song, M. Zhu, X. Li, Z. Zeng, K. Luo, Y. Jiang, Z. Zhang, C. Li, Y. Wang, B. Li, Z. Wang, Z. Zhang, and W. Hu, Science Advances 11, eady2349 (2025).
- [79]A. Shilov, “ASML’s High-NA chipmaking tool will cost $380 million,” Tom’s Hardware, https://www.tomshardware.com/tech-industry/manufacturing (2024).
- [80]J. Bhattacharyya and M. Packalen, “Stagnation and scientific incentives,” Tech. Rep. 26752 (National Bureau of Economic Research, 2020).
- [81]Grand View Research, “Photonic integrated circuit market size & share report, 2024–2030,” Market analysis report, https://www.grandviewresearch.com/industry-analysis/photonic-integrated-circuit-market (2024).
- [82]IEEE International Roadmap for Devices and Systems, “International roadmap for devices and systems (IRDS), 2024 edition,” IEEE, https://irds.ieee.org (2024).
- [83]B. Nagy, J. D. Farmer, Q. M. Bui, and J. E. Trancik, PLoS ONE 8, e52669 (2013).