To evaluate the effectiveness of our proposal concept and spatiotemporal model, we conducted TPDOPL using a femtosecond pulse laser and a DMD (Fig. 3A, fig. S7). The DMD, with a megapixel-resolution projection layout of arbitrary features, is irradiated by a flat-top beam and focused onto a photoresist film on a cover glass using an oil-immersion objective lens (Nikon, 100×, NA 1.49). We chose a commercially available non-chemically amplified (non-CA) negative photoresist (AR-N 7520) because its degree of polymerization, driven by the stepwise photopolymerization mechanism, is easily controllable and quantifiable for NeTPA under few-photon irradiation. This resist has an absorption peak at 323 nm and an absorption cut-off wavelength of 353 nm (fig. S8), ensuring that only TPA occurs when using femtosecond pulse lasers at both 400 nm and 517 nm wavelengths. Utilizing this system, a uniform exposure of 80 × 100 µm2 can be achieved in a single exposure field, as demonstrated in fig. S9. Note that the number of excitable molecules in the photoresist by TPA should be less than the calculated NeTPA from the proposed model. The calculated NeTPA predicts the possible opportunity and distribution of eTPA under photon irradiation from the viewpoint of incident photons, but it ignores the molecular concentration, distribution, and quantum yield of TPA in the photoresist. Furthermore, the conversion efficiency from excited molecules by TPA to the practically initiated coupling reaction between molecules should be considered. The number of reaction sites of the photosensitive molecules is ultimately limited by the final absorbed NeTPA and their quantum yield to initiate the reaction.
We utilized incident light with an average power of 1 mW after passing through the objective lens as an illustrative example, specifically with parameters Nspp = 1250 (0.48 fJ/pulse pixel) and Npulse = 1 × 107. The distribution of NeTPA for a line composed of one pixel demonstrates a notable 24% reduction in full width at half maximum (FWHM) compared to the photon distribution (Fig. 3A). According to photopolymerization theory33,34, the relationship between the concentration of photosensitive molecules (M) excited by TPA and the photon distribution is nonlinear. As the reaction step is repeated, the molecular weight at the center of the exposure field increases exponentially due to the chemical cross-linking reaction. The concentration distribution of photosensitive molecules involved in the reaction is illustrated in Fig. 3B (Supplementary Information S3 for more details). The exponential increase in molecular weight leads to faster gelation at the exposure field's center, forming insoluble polymer networks more quickly than in the surroundings. Consequently, the superposition and coordination of optical and chemical nonlinearity can effectively reduce the feature size in TPDOPL under few-photon irradiation.
To investigate the effectiveness of the proposed spatiotemporal model for TPDOPL under few-photon irradiation, we fabricated separate lines using our TPDOPL system with a femtosecond laser wavelength (λ) of 517 nm and pulse width of 238 fs. Using a single-pixel DMD layout, a line with an average width of 41 nm and a minimum feature size of 28 nm (Fig. 3C) was achieved under the irradiation of a total incident photon number per pixel of 4.37 × 1011 (0.167 µJ) with accumulation Npulse of 8.5 × 107 pulses containing Nspp of 5.14 × 103 (1.97 fJ/pulse pixel). Correspondingly, we calculated the eTPA distribution using the same photon flux as the experimental result in Fig. 3C but only performed 8.5 × 103 pulses in simulation. The simulation result shown in Fig. 3D indicates that the eTPA distribution is concentrated in a central area of about 30 nm. This validates the effectiveness of our spatiotemporal model for predicting the feature size of TPDOPL under few-photon irradiation.
Photon irradiance density and the accumulated pulse numbers critically influence the line width of TPDOPL. By decreasing Nspp from 1.12×104 (4.30 fJ/pulse pixel) to 6.52 × 103 (2.51 fJ/pulse pixel), the average line width of the polymer line was reduced from 164 nm to 43 nm under the accumulation of 6 × 107 pulses, achieving a minimum feature size of 26 nm (1/20 λ), as shown in Fig. 3E. The NeTPA under different Nspp irradiations can be observed in fig. S10. The relationship between the polymer line width and photon irradiance density is depicted in Fig. 3F, indicating that the feature size can be reduced by decreasing the photon irradiance density. However, lower photon irradiance density may increase line roughness due to quantum noise, which can increase edge roughness for extremely fine lines (fig. S11). On the other hand, increasing the accumulation of pulses with a fixed photon flux density leads to a widening of the line width, as shown in Fig. 3G.
Another significant aspect pertains to periodic lines in photolithography, which determine the potential feature density achievable in device applications. Generally, the minimum distinguishable period between adjacent lines is dependent on the wavelength and determined by the equation HP (half pitch) = 0.5 λ/NA, following the Sparrow criterion35. When the design pattern period is less than the minimum resolvable distance between two lines, double patterning lithography (DPL) can overcome this problem36. For instance, at λ = 517 nm and NA = 1.45, this criterion yields an approximate value of 217 nm. We designed a line array using the DMD pixel period of 7.56 µm combined with 2 pixels on and 1 pixel off periodically (fig. S12A), corresponding to a period of 226.8 nm. Using irradiation conditions with Npulse = 6 × 107 and Nspp = 1.52 × 104 (5.84 fJ/pulse pixel), the lines were indistinguishable (fig. S12C). We efficiently utilized the flexibility of TPDOPL by using a DMD as a digital mask, enabling in-situ digital multiple exposures (iDME) to print dense features without being constrained by the diffraction limit. Exploiting DMD characteristics, two split layouts with a period of 2‘p’ are sequentially loaded in situ for double exposure, achieving an exposure result with a period of ‘p’, as depicted in Fig. 4A. Under twice alternating exposure of Npulse = 6 × 107 and Nspp = 8.53 × 103 (3.28 fJ/pulse pixel), we successfully achieved a dense line array with a period of 210 nm (HP ~ 0.3 λ/NA), a linewidth of 150 nm, and a gap spacing of 60 nm, as shown in Fig. 4B, surpassing the diffraction limit.
Taking advantage of TPDOPL-iDME, we can achieve distinguishable dense structure patterning. When the pitch is less than 5 pixels, a single exposure cannot meet the resolution consistent with the design pattern (fig. S13). Figure 4C shows a typical circuit layout selected from a commercial chip, including isolated and dense lines with a width of 3 or 7 pixels and intervals of 1 and 2 pixels between lines (fig. S14). We employ algorithms37 to strategically distribute polygons with interspacing distances below 2 pixels across distinct sub-masks, optimizing their arrangement for TPDOPL-iDME. SEM images show that direct single exposure causes indistinguishable results in dense line areas (Fig. 4D). By splitting this layout into two (Fig. 4E) and performing our TPDOPL-iDME approach, we successfully achieved the expected circuit patterning (Fig. 4F). The dense lines are clearly distinguished, and the periods agree well with the design. Furthermore, by optimizing exposure parameters and layout design for TPDOPL-iDME, line width, period, and gap distance can be controlled for finer and denser feature patterning.
Optical devices with curved and circular microstructures have been fabricated using TPDOPL, such as patterns including arrayed waveguide gratings and micro-ring resonators38. The radius of the ring affects the value of the free spectral range, and the gap or spacing between the guide and the ring affects the coupling ratio between the waveguide and the ring39. Through layout design and the TPDOPL-iDME method, we can fabricate micro-ring filters with varied radius pitches. The widths of the circular rings can be adjusted from 220 nm to 346 nm by increasing Npulse under the irradiation of Nspp = 8.53×103 (3.28 fJ/pulse pixel) (fig. S15 A). We patterned the line waveguides first, then fabricated circular rings with different diameters (Fig. 5A), leveraging TPDOPL-iDME. The gap distances between the line and circular rings can be finely adjusted from 66 nm to 480 nm (fig. S15B), optimizing the structures and improving the properties of photonic resonance devices.
The flexibility of TPDOPL-iDME allows us to create arbitrary patterns with various sizes, shapes, and densities, applicable not only in microelectronics and microphotonics but also in microfluidics40,41. Microfluidics in microbiology offer an in vitro platform for interactions among diverse cell types, enabling real-time observation and assessment of reaction processes42. We designed a rectangular module to substitute the cell chamber and a circular module to replace the cell secretion chamber, with channels of varied sizes to facilitate the addition and observation of multiple cell types and their reactions43. Figure 5B shows complex patterns of biological microfluidics fabricated by TPDOPL-iDME, where square cell incubators (3 × 3 µm²), rectangular cell chambers (2.8 × 6 µm²), and circular cell collectors with micrometer and sub-micrometer scales are connected by different channels with widths from 70 nm to 800 nm (Fig. 5B iii), effectively carrying and separating viruses of different sizes. Most biomolecular analytes are below microns in size44, especially foreign objects such as viruses45, which are usually 20–300 nm in size. Cross-scale biological microfluidics, from micrometer to nanometer, hold promise for providing research platforms for new diagnostic and therapeutic methods for viruses like the new coronavirus.