research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Evaluation of the X-ray/EUV Nanolithography Facility at AS through wavefront propagation simulations

crossmark logo

aDepartment of Mathematical and Physical Sciences, School of Computing, Engineering and Mathematical Sciences, La Trobe University, Bundoora, Victoria 3086, Australia, and bAustralian Synchrotron, Australian Nuclear Science and Technology Organisation (ANSTO), Clayton, Victoria 3168, Australia
*Correspondence e-mail: j.knappett@latrobe.edu.au, g.vanriessen@latrobe.edu.au

Edited by M. Yabashi, RIKEN SPring-8 Center, Japan (Received 16 November 2023; accepted 18 March 2024; online 17 April 2024)

Synchrotron light sources can provide the required spatial coherence, stability and control to support the development of advanced lithography at the extreme ultraviolet and soft X-ray wavelengths that are relevant to current and future fabricating technologies. Here an evaluation of the optical performance of the soft X-ray (SXR) beamline of the Australian Synchrotron (AS) and its suitability for developing interference lithography using radiation in the 91.8 eV (13.5 nm) to 300 eV (4.13 nm) range are presented. A comprehensive physical optics model of the APPLE-II undulator source and SXR beamline was constructed to simulate the properties of the illumination at the proposed location of a photomask, as a function of photon energy, collimation and monochromator parameters. The model is validated using a combination of experimental measurements of the photon intensity distribution of the undulator harmonics. It is shown that the undulator harmonics intensity ratio can be accurately measured using an imaging detector and controlled using beamline optics. Finally, the photomask geometric constraints and achievable performance for the limiting case of fully spatially coherent illumination are evaluated.

1. Introduction

As of 2019, the semiconductor device manufacturing industry has adopted lithography technology utilizing extreme ultraviolet (EUV) radiation of wavelength λ = 13.5 nm, corresponding to a photon energy of 91.8 eV, for high-volume manufacturing (Fomenkov, 2019[Fomenkov, I. (2019). Proceedings of the 2019 ASML Source Workshop, 5 November 2019, Amsterdam, The Netherlands, p. S1.]). To keep up with the demands of device scaling predicted by Moore's law, a future transition to 6.7 nm wavelength (185 eV) sources is anticipated. The availability of laser-pulsed plasma light sources (Otsuka et al., 2012[Otsuka, T., Li, B., O'Gorman, C., Cummins, T., Kilbane, D., Higashiguchi, T., Yugami, N., Jiang, W., Endo, A., Dunne, P. & O'Sullivan, G. (2012). Proc. SPIE, 8322, 832214.]) and multilayer optics (Uzoma et al., 2021[Uzoma, P. C., Shabbir, S., Hu, H., Okonkwo, P. C. & Penkov, O. V. (2021). Nanomaterials, 11, 2782.]) at this wavelength make it a particularly promising candidate for future lithography. A shift to shorter wavelengths brings with it significant challenges, including an increased significance of stochastic effects for higher resolution patterning (De Bisschop, 2017[De Bisschop, P. (2017). J. Micro/Nanolith. MEMS MOEMS, 16, 041013.]), and a need for understanding the wavelength-dependent effects of mask defects (Goldberg & Mochi, 2010[Goldberg, K. A. & Mochi, I. (2010). J. Vac. Sci. Technol. B, 28, C6E1-C6E10.]). Interference lithography (IL) using synchrotron radiation has recently emerged as a powerful tool for understanding the challenges for future lithographic processes, including photoresist performance from EUV (Mojarad et al., 2015a[Mojarad, N., Gobrecht, J. & Ekinci, Y. (2015a). Sci. Rep. 5, 9235.]) to soft X-ray (SXR) wavelengths (2.5 nm) (Mojarad et al., 2021[Mojarad, N., Kazazis, D. & Ekinci, Y. (2021). J. Vac. Sci. Technol. B, 39, 042601.]). Synchrotron sources are ideal for the development of lithographic technology as they allow for the control of properties such as flux, coherence and polarization, and can cover the entire energy range from EUV to SXR. Partially coherent wavefront propagation simulation can provide critical insight into the design of optical systems and EUV/SXR-IL process development using existing synchrotron radiation sources.

For high-resolution patterning using IL, it is critical to achieve high contrast in the aerial image formed where beams diffracted from two or more gratings interfere. This requires that light sources used in IL provide high intensity, spatial coherence and stability (Mojarad et al., 2015b[Mojarad, N., Gobrecht, J. & Ekinci, Y. (2015b). Microelectron. Eng. 143, 55-63.]). The necessary spatial coherence length to ensure high-contrast aerial images is determined by the size of the desired exposure area (Solak et al., 2003[Solak, H., David, C., Gobrecht, J., Golovkina, V., Cerrina, F., Kim, S. & Nealey, P. (2003). Microelectron. Eng. 67-68, 56-62.]), which is, in turn, defined by the maximum separation between points in a set of two or more radially arranged gratings. Fully coherent illuminaton of the grating set typically requires a coherence length three times the linear dimensions of the desired exposure area (Solak et al., 2002[Solak, H. H., David, C., Gobrecht, J., Wang, L. & Cerrina, F. (2002). Microelectron. Eng. 61-62, 77-82.]). Typical grating masks used for EUV-IL require a lateral coherence length at the mask plane of between 150 µm and 1.2 mm (Ekinci et al., 2014[Ekinci, Y., Vockenhuber, M., Mojarad, N. & Fan, D. (2014). Proc. SPIE, 9048, 904804.]; Meng et al., 2021[Meng, X., Yu, H., Wang, Y., Ren, J., Xue, C., Yang, S., Guo, Z., Zhao, J., Wu, Y. & Tai, R. (2021). J. Synchrotron Rad. 28, 902-909.]), although small area patterning using EUV-IL has been reported using a coherence length of just 43.2 µm (Sahoo et al., 2023[Sahoo, A. K., Chen, P., Lin, C., Liu, R., Lin, B., Kao, T., Chiu, P., Huang, T., Lai, W., Wang, J., Lee, Y. & Kuan, C. (2023). Micro Nano Eng. 20, 100215.]).

When evaluating the performance of a light source for interference lithography, the relevant quantity that determines throughput is the exposure time required to provide the radiation dose-on-wafer Dw(λ) to transfer a pattern to a film of photosensitive resist with particular sensitivity at a particular wavelength λ. This, in turn, depends on the coherent intensity at the mask and mask efficiency. Practically, for studies of lithographic processes at high resolution, through­put is of lower priority than lithographic quality, and the exposure time determines the required mechanical stability between the mask and wafer.

Photoresist sensitivity is defined by the dose D0(λ) required to remove 50% of the resist thickness during development for a particular wavelength (Ekinci et al., 2013[Ekinci, Y., Vockenhuber, M., Hojeij, M., Wang, L. & Mojarad, N. (2013). Proc. SPIE, 8679, 867910.]; Mojarad et al., 2015a[Mojarad, N., Gobrecht, J. & Ekinci, Y. (2015a). Sci. Rep. 5, 9235.]). Resist sensitivity generally decreases for shorter wavelengths, i.e., D0 ∝ 1/λ. Combined with the requirement for higher doses for patterning thinner resists and smaller feature sizes, the scaling of source power requirements with industry device scaling targets is a major challenge (Van Schoot, 2021[Van Schoot, J. (2021). Proc. SPIE, 11609, 1160905.]; Levinson, 2022[Levinson, H. J. (2022). Jpn. J. Appl. Phys. 61, SD0803.]). For the present work, it is convenient to define the dose-on-mask, Dm, required to achieve a target value of D0, taking into account diffraction efficiency, η(λ), of the grating mask and the transmission of the substrate, T(λ), at the relevant wavelength,

[D_{\rm m} = {{D_{0}} \over {\eta TN^{\,3}}}, \eqno(1)]

where N is the number of gratings contained in the mask and the power of three is due to the definition of the efficiency of an IL grating, ηIL = Nη (Mojarad et al., 2015b[Mojarad, N., Gobrecht, J. & Ekinci, Y. (2015b). Microelectron. Eng. 143, 55-63.]), and periodicity in the aerial image intensity. The wavelength dependence has been suppressed for simplicity. While diffraction efficiencies of up to 28% have been reported for low-resolution IL gratings (Braig et al., 2011[Braig, C., Käsebier, T., Kley, E.-B. & Tünnermann, A. (2011). Opt. Express, 19, 14008-14017.]), Wang et al. (2021[Wang, X., Kazazis, D., Tseng, L.-T., Robinson, A. P. & Ekinci, Y. (2021). Nanotechnology, 33, 065301.]) reported diffraction efficiencies for a bilayer of hydrogen silsesquioxane and spin-on-carbon of η > 6% for gratings with half-pitch down to 12 nm at EUV wavelengths.

Organic chemically amplified resists (CARs) have been used for performance comparisons of EUV and SXR, as they offer high-resolution patterning in both wavelength ranges (Mojarad et al., 2013[Mojarad, N., Vockenhuber, M., Wang, L., Terhalle, B. & Ekinci, Y. (2013). Proc. SPIE, 8679, 867924.]). Mojarad et al. (2015a[Mojarad, N., Gobrecht, J. & Ekinci, Y. (2015a). Sci. Rep. 5, 9235.]) reported CAR sensitivities of D0 (13.5 nm) = 11.6 mJ cm−2 and D0 (6.5 nm) = 33.2 mJ cm−2. The 2021 International Roadmap for Devices and Systems has predicted that D0 > 80 mJ cm−2 will be required for 10 nm half-pitch patterning (Levinson, 2022[Levinson, H. J. (2022). Jpn. J. Appl. Phys. 61, SD0803.]). Assuming a SiO2 substrate of 40 nm thickness and a four-grating mask with η = 6% efficiency, the corresponding minimum dose on mask for EUV-IL can be estimated to be Dm ≃ 11 mJ cm−2.

In this work, partially coherent wavefront propagation was used to simulate the complex wavefield at the mask plane of an EUV-IL instrument being developed at La Trobe University for one of two branches of the SXR beamline of the Australian Synchrotron (AS). Simulations were performed primarily at 90.44 eV and 184.76 eV photon energy and the relevant characteristics of the beam were compared with direct measurements of the intensity taken at the beamline with photon energy ranging from 90 eV to 300 eV. We demonstrate, by comparison with photoelectron spectroscopy measurements, that the ratio of undulator harmonic intensity can be measured using an imaging detector. We discuss how the higher harmonic intensity scales with decreasing photon energy and is affected by collimating slits also used to control spatial coherence. The beamline's suitability for EUV-IL is evaluated with respect to the compromise between the grating area that can be coherently illuminated and the intensity of the illumination.

2. Theory and modelling

2.1. Interference lithography

Light incident on a grating diffracts at multiple angles, θm, corresponding to orders, m, governed by the grating equation. For a monochromatic beam of wavelength λ, incident at angle θi, on a grating with period, pG, the angle of diffraction of the mth order beam is (Hecht, 1987[Hecht, E. (1987). Optics, 2nd ed. Addison-Wesley.])

[\sin\theta_{m} = {{m\lambda}\over{p_{\rm{G}}}}+\sin\theta_{\rm i}. \eqno(2)]

As shown in Fig. 1[link], the mth-order beams from each grating will interfere at the aerial image plane at a distance zm from the mask, given by

[z_{m} = {{d}\over{2\tan\theta_{m}}} = {{ d\,\left(\,p_{\rm{G}}^{2}-m^{2}\lambda^{2}\right)^{1/2} }\over{ 2m \lambda }}. \eqno(3)]

For the simplest IL setup, consisting of two plane waves intersecting at a 2θm angle, the aerial image will have a period p given by

[p = {{\lambda}\over{2\sin\theta_{m}}}. \eqno(4)]

The standard measure of the resolution of an aerial image is the half-pitch, HP = p/2. For a grating illuminated at normal incidence, equations (2)[link] and (4)[link] can be combined to show that

[{\rm{HP}} = {{\lambda}\over{4\sin\theta_{m}}} = {{p_{\rm{G}}}\over{4m}}. \eqno(5)]

The ultimate resolution for a first-order aerial image is therefore pG/4 (Karim et al., 2015[Karim, W., Tschupp, S. A., Oezaslan, M., Schmidt, T. J., Gobrecht, J., van Bokhoven, J. A. & Ekinci, Y. (2015). Nanoscale, 7, 7386-7393.]). Equation (5)[link] implies that sub-10 nm patterning by IL requires a grating mask with pG < 40 nm, which can be readily fabricated using electron beam lithography (Vieu et al., 2000[Vieu, C., Carcenac, F., Pépin, A., Chen, Y., Mejias, M., Lebib, A., Manin-Ferlazzo, L., Couraud, L. & Launois, H. (2000). Appl. Surf. Sci. 164, 111-117.]). It also shows the achromatic nature of IL, as the right side of equation (5)[link] is independent of λ.

[Figure 1]
Figure 1
Schematic of an interference lithography setup using a binary grating mask. The first-order (m = ±1) beams form an aerial image at the image plane, a distance z1 from the mask, which is then transferred onto a photoresist. The zero-order (m = 0) and second-order (m = 2) diffracted beams are also shown.

2.2. Source model

Synchrotron Radiation Workshop (SRW) (Chubar & Elleaume, 1998[Chubar, O. & Elleaume, P. (1998). Proceedings of the Sixth European Particle Accelerator Conference (EPAC'98), 22-26 June 1998, Stockholm, Sweden, pp. 1177-1179. THP01G.]) was used to construct a model of branch B of the SXR beamline at the AS, shown in Fig. 2[link] (Knappett, 2021[Knappett, J. (2021). Wavefront Propagation Simulation of X-ray Interference Nanolithography. Master's thesis, La Trobe University, Australia.]). A new model discussed in the present work includes the proposed EUV-IL endstation, permitting evaluation of the influence of all beamline parameters on the properties of the aerial image. SRW is commonly used to model SXR/EUV synchrotron beamlines and is capable of fully or partially coherent propagation (He et al., 2020[He, A., Chubar, O., Dvorak, J., Hu, W., Hulbert, S., Walter, A., Wilkins, S. & Salort, C. (2020). Proc. SPIE, 11493, 1149307.]; Meng et al., 2021[Meng, X., Yu, H., Wang, Y., Ren, J., Xue, C., Yang, S., Guo, Z., Zhao, J., Wu, Y. & Tai, R. (2021). J. Synchrotron Rad. 28, 902-909.]).

[Figure 2]
Figure 2
Schematic illustration of the SXR beamline. The distance of each element from the source is indicated (not to scale). Proposed IL optics and an imaging detector are included in branch B. Branch A was used to obtain measurements of the undulator harmonic intensity using photoelectron spectroscopy.
2.2.1. Storage ring and electron beam

The parameters used to model the electron beam and storage ring were obtained from Wootton et al. (2012[Wootton, K. P., Boland, M. J., Dowd, R., Tan, Y. E., Cowie, B. C. C., Papaphilippou, Y., Taylor, G. N. & Rassool, R. P. (2012). Phys. Rev. Lett. 109, 194801.], 2013[Wootton, K. P., Rassool, R. P., Boland, M. J. & Cowie, B. C. C. (2013). Proceedings of the 2nd International Beam Instrumentation Conference (IBIC2013), 16-19 September 2013, Oxford, UK, pp. 546-549. TUPF19.]) and are listed in Table 1[link].

Table 1
Key values for the Australian Synchrotron storage ring; parameters obtained from Wootton et al. (2012[Wootton, K. P., Boland, M. J., Dowd, R., Tan, Y. E., Cowie, B. C. C., Papaphilippou, Y., Taylor, G. N. & Rassool, R. P. (2012). Phys. Rev. Lett. 109, 194801.])

  Parameter Value
E0 Beam energy 3.01 GeV
I0 Current 0.2 A
σE RMS energy spread 0.1021%
εx Horizontal emittance 10 nm rad
εy Vertical emittance 0.009 nm rad
βx Horizontal beta function 9 m
βy Vertical beta function 3 m
2.2.2. APPLE-II undulator

The SXR beamline utilizes an elliptically polarizing APPLE-II undulator source. This undulator was modelled as ideal, with magnetic fields defined by the undulator length, L, magnetic period, λu, and number of periods, Nu. The undulator parameters were obtained from Wootton et al. (2012[Wootton, K. P., Boland, M. J., Dowd, R., Tan, Y. E., Cowie, B. C. C., Papaphilippou, Y., Taylor, G. N. & Rassool, R. P. (2012). Phys. Rev. Lett. 109, 194801.], 2013[Wootton, K. P., Rassool, R. P., Boland, M. J. & Cowie, B. C. C. (2013). Proceedings of the 2nd International Beam Instrumentation Conference (IBIC2013), 16-19 September 2013, Oxford, UK, pp. 546-549. TUPF19.]) and are shown in Table 2[link].

Table 2
Key values for the APPLE-II undulator used for the SXR beamline at the Australian Synchrotron

Ku, which determines Bu, was adjusted to produce a fundamental harmonic at 184.76 eV. Other values obtained from Wootton et al. (2012[Wootton, K. P., Boland, M. J., Dowd, R., Tan, Y. E., Cowie, B. C. C., Papaphilippou, Y., Taylor, G. N. & Rassool, R. P. (2012). Phys. Rev. Lett. 109, 194801.]).

  Parameter Value
Lu Undulator length 1.875 m
λu Undulator period length 75 mm
Nu Number of undulator periods 25
Bu Peak magnetic field 0.4611 T
Ku Deflection parameter 3.230

2.3. Beamline model

A schematic of the optical elements in use for each branch of the SXR beamline is shown in Fig. 2[link]. The key optical components include white-beam slits (WBS), a toroidal mirror, a planar grating monochromator (PGM) and a cylindrical mirror that can be rotated to direct the beam to branch A or B. The schematic shows the proposed location of IL optics in branch B, along with an existing EUV/SXR area detector. For typical operating conditions of branch B, the WBS are wide open (not touching the beam) and the bandwidth of the beam is defined by the PGM and secondary source aperture (SSA). The PGM disperses the monochromatic components of the incident beam at angles dependent on the photon energy as defined in equation (2)[link], the beam is then focused at the SSA, the size of which defines the resolving power of the beam as (λ/dλ). All relevant parameters are listed in Table 3[link].

Table 3
Typical geometry of the optical elements used in the model of the beamline and IL optics

The defining characteristics of a mirror are given as the major radius, R, and minor radius, ρ. The defining characteristic of a grating is the periodicity, and of a mask is the grating size, G. Aperture sizes shown were used for all simulations and measurements shown in this work unless stated otherwise.

Element name Element type Propagation distance (m) Distance from source (m) Dimensions (mm) Defining characteristics Incident angle (°)
White-beam slits Aperture 14.30 14.30 4 × 3 N/A 90
Toroidal mirror Mirror 0.82 15.12 420  ×  30 R = 6668 m, ρ = 5.262 m 1
Planar mirror Mirror 3.38 18.50 460 × 50 N/A 2
Grating Grating 0.20 18.70 150 × 20 250 lines mm−1 1.173
Exit aperture Aperture 1.00 19.70 10 × 20 N/A 90
Cylindrical mirror Mirror 0.60 20.30 240 × 40 R = 100 km, ρ = 0.2443 m 1.5
Secondary source aperture Aperture 6.50 26.80 0.025 × 0.025 N/A 90
Grating mask Mask 9.50 36.30 0.15 × 0.15 G = 50 µm 90

2.4. Partially coherent simulations

SRW implements fully spatially coherent wavefront propagation through the propagation of monochromatic radiation emitted by a single relativistic electron passing through an undulator. The wavefield is represented by the monochromatic, transverse components of the electric field Eλ and paraxial propagation is performed through fast Fourier transforms (Chubar & Elleaume, 1998[Chubar, O. & Elleaume, P. (1998). Proceedings of the Sixth European Particle Accelerator Conference (EPAC'98), 22-26 June 1998, Stockholm, Sweden, pp. 1177-1179. THP01G.]). Taking x and y as the horizontal and vertical directions, respectively, and z as the direction along the beam axis, the propagated wavefield, Eλ(x, y, z) from each wavefield component at z = 0, Eλ(x, y, z = 0), can be written as

[\eqalignno{ E_{\perp\lambda}(x,y,z) = {}& -{{ik\exp(ikz)}\over{2\pi{z}}} \int\!\!\!\int\limits_{-\infty}^{\infty} E_{\perp\lambda}(x^{\prime},y^{\prime},z=0) \cr & \times \exp\left\{ {{ik}\over{2z}}\left[(x-x^{\prime})^{2}+(y-y^{\prime})^{2}\right] \right\} \,{\rm{d}}x^{\prime}\,{\rm{d}}y^{\prime}, & (6) }]

where k = (kx 2+ky 2+kz 2)1/2 is the spatial frequency of the plane wave component over an area of phase space dkxdky. In all simulations shown in this work, only a single monochromatic component is considered, where the spatial frequency k = 2π/λ, and a quadratic approximation to the exponential phase term in equation (6)[link] (Chubar & Celestre, 2019[Chubar, O. & Celestre, R. (2019). Opt. Express, 27, 28750-28759.]) was used where appropriate.

Partially coherent monochromatic propagation is achieved by randomly sampling the phase space of the entire electron beam using the Monte Carlo method (Laundy et al., 2014[Laundy, D., Alcock, S. G., Alianelli, L., Sutter, J. P., Sawhney, K. J. S. & Chubar, O. (2014). Proc. SPIE, 9209, 920903.]). The number of electrons needed to accurately represent the electron beam increases with the electron beam emittance (Chubar et al., 2016[Chubar, O., Chu, Y. S., Huang, X., Kalbfleisch, S., Yan, H., Shaftan, T., Wang, G., Cai, Y. Q., Suvorov, A., Fluerasu, A., Wiegart, L., Chen-Wiegart, Y., Thieme, J., Williams, G., Idir, M., Tanabe, T., Zschack, P. & Shen, Q. (2016). AIP Conf. Proc. 1741, 040002.]). For all simulations shown in this work, 4000 electrons were sufficient to accurately model the monochromatic beam (Knappett, 2021[Knappett, J. (2021). Wavefront Propagation Simulation of X-ray Interference Nanolithography. Master's thesis, La Trobe University, Australia.]). The electric field of each electron is then propagated coherently using equation (6)[link] and the partially coherent electric field is calculated as the average of the electric field distributions from each electron, seeded over the phase space occupied by the beam (Chubar, 2014[Chubar, O. (2014). Proc. SPIE, 9209, 920907.]).

3. Simulation results

3.1. Beam profile

The beamline model was used to propagate a partially coherent wavefield from the source through to the grating mask plane at λ = 4.96, 6.7, 9.2 and 13.5 nm. The intensity was modelled from the complex electric field as I(x, y) = |E(x, y)|2, with the pixel size chosen such that it satisfies the Nyquist–Shannon sampling theorem to ensure satisfactory representation of phase (Shannon, 1949[Shannon, C. (1949). Proc. Inst. Radio Eng. (IRE), 37, 10-21.]). The full width at half-maximum (FWHM) of the intensity profile at the mask plane was calculated for each wavelength.

3.2. Flux and coherence at the mask plane

The total flux and spatial coherence of the wavefield at the mask plane was evaluated after partially coherent propagation through the beamline model. The total flux, Φ, was calculated in units of photons per second per 0.1% bandwidth from the propagated complex electric field E(x, y) by

[\eqalignno{ \Phi\ \big[{\rm{photons}}\ {\rm{s}}^{-1}\ & (0.1\%\,\, {\rm{bandwidth}})^{-1}\big] \cr & =\, {\rm{dx}}\,{\rm{d}}y\times10^{-4}\sum\limits_{x}\sum\limits_{y}\big|E(x,y)\big|^{2}, & (7) }]

where dx and dy are the pixel size of E in x and y, respectively. The flux was corrected for the reflectivity of each mirror and efficiency of the grating used in the PGM, shown in Table 3[link]. Since the reflection grating efficiency at EUV/SXR wavelengths is not accurately known, a constant efficiency of 10% was assumed for all simulations, with a third-order efficiency of 0.1%.

The spatial coherence length lc at the mask plane was calculated by first computing the one-dimensional mutual intensity J using

[J\left(x_{1},x_{2}\right) = \left\langle E^{\,*}(x_{1}) \, E(x_{2}) \right\rangle, \eqno(8)]

where x1 and x2 are different sets of points in the central horizontal axis (x, y = 0), so that when x1 = x2, J reduces to the horizontal intensity profile. The degree of coherence between any point along (x, y = 0) and the central point (x = 0, y = 0) was then calculated following Meng et al. (2021[Meng, X., Yu, H., Wang, Y., Ren, J., Xue, C., Yang, S., Guo, Z., Zhao, J., Wu, Y. & Tai, R. (2021). J. Synchrotron Rad. 28, 902-909.]) using

[\gamma = {{ J(x,0) }\over{ \left[\,\left\langle|E(x)|^{2}\right\rangle \,\left\langle|E(0)|^{2}\right\rangle\,\right]^{1/2} }}. \eqno(9)]

The coherence length lc(x, y) can then be determined as the distance from the central point at which γ < 0.8. Repeated propagations were undertaken for different horizontal SSA sizes, ranging from 25 µm to 350 µm. The results of the simulations are shown in Fig. 3[link], which shows that the total flux at the mask plane increases linearly with horizontal SSA size from 2.7 × 109 photon s−1 to 2.9 × 1010 photons s−1, while the horizontal coherence length at the mask plane decreases inversely proportional to the slit width, from ∼1.17 mm to ∼125 µm. Previous unpublished measurements by the authors using a Young's double-slit array with a maximum slit separation of 30 µm showed no significant loss of visibility with any horizontal SSA size, indicating that the horizontal coherence length at 6.7 nm wavelength is always greater than 30 µm. The finite horizontal coherence may impose a practical limit on the horizontal extent of an IL grating mask as it will generally be smaller than the vertical coherence length due to the asymmetric source emittance (Table 1[link]).

[Figure 3]
Figure 3
(a) Total photon flux (red), horizontal coherence length (blue), and (b) resolving power at the mask plane of beamline branch B for different SSA sizes. The resolving power is shown as a function of SSA height because the energy resolution of the PGM is defined only by the vertical SSA size. The grey dotted lines correspond to the maximum SSA size that provides coherent illumination of a representative grating mask (see Fig. 7[link]).

The bandwidth was measured using X-ray photoelectron spectroscopy (XPS) of the Au Fermi edge. The Fermi edge broadening was measured with different vertical SSA sizes and the energy bandwidth of the fundamental undulator harmonic at λ = 6.7 nm was estimated for each SSA size. The results are shown in Fig. 3[link].

4. Experimental measurements

Measurements of the intensity profile of the beam were taken at branch B of the beamline (Fig. 2[link]), using an AXIS-SXR sCMOS detector with an EUV-Enhanced GSENSE400 BSI sensor (Harada et al., 2019[Harada, T., Teranishi, N., Watanabe, T., Zhou, Q., Bogaerts, J. & Wang, X. (2019). Appl. Phys. Expr. 13, 016502.]). The detector was situated z = 0.5 m from the beam-defining aperture (BDA) plane. At the BDA plane, a 4.064 µm-thick Ultralene (C3H6) filter was positioned in the beam path to attenuate the beam and avoid over-saturating the detector. For each measurement, the WBS were set to 0.84 mm × 1.0 mm (x, y), and the SSA was set to 25 µm × 25 µm (x, y). The PGM fixed focus constant (Cff) was set to 1.4 instead of the standard operating value of 2, as the higher harmonic content (ζ) is reduced for small Cff values (Kleemann et al., 1997[Kleemann, B. H., Gatzke, J., Jung, C. & Nelles, B. (1997). Proc. SPIE, 3150, 137-147.]); all other parameters were as shown in Table 3[link]. These small slit sizes mean all results presented in this section show less flux and intensity than can be expected with beamline parameters optimized for high flux. The photon energy was varied from 90 eV to 270 eV in 10 eV steps, with extra steps at 92 eV (13.5 nm) and 185 eV (6.7 nm). One hundred exposures were taken at each photon energy which were processed in sets of ten. Horizontal and vertical profiles across the summed intensity in each set were fitted with two Gaussians, as shown in Fig. 4[link], which represent the fundamental (harmonic n = 1) intensity profile and the intensity profile of the third (n = 3) harmonic. For the fitted Gaussian representing the fundamental, the FWHM was calculated (Fig. 5[link]).

[Figure 4]
Figure 4
The measured horizontal intensity profile, 50 cm from the BDA with a fundamental photon energy of 185 eV and Cff = 2. The horizontal line profile through the beam, overlaid with a fitting of two Gaussians representing the fundamental (harmonic n = 1) and the third harmonic (n = 3).
[Figure 5]
Figure 5
The intensity FWHM measured in the horizontal (red) and vertical (black) direction for photon energies from 90 eV to 270 eV and Cff = 1.4. The calculated FWHM of intensity profiles generated through partially coherent wavefront propagation through the beamline model are also shown for photon energies of 90.44 eV, 135 eV, 184.76 eV and 250 eV.

The total photon flux of each harmonic, Φn, was also calculated, as well as the intensity over a four-grating mask area, [I^{\,{\rm{G}}}_{n}]. The four-grating mask, typical of those used to produce two-dimensional patterns by IL (Mojarad et al., 2015b[Mojarad, N., Gobrecht, J. & Ekinci, Y. (2015b). Microelectron. Eng. 143, 55-63.], 2021[Mojarad, N., Kazazis, D. & Ekinci, Y. (2021). J. Vac. Sci. Technol. B, 39, 042601.]), consists of 50 µm × 50 µm gratings, arranged radially around the beam centre, with a 50 µm × 50 µm central beam stop. A schematic of the grating mask area over the two-dimensional beam intensity is shown later in Fig. 7. Errors for each fitted parameter have been taken as the standard deviation of the fitted value for each set of ten processed images. For each experimentally measured beam parameter discussed henceforth, equivalent parameters extracted from the partially coherent simulated intensity at 90.44, 135, 184.76 and 250 eV have been included for comparison.

The total flux was measured again as a function of photon energy using the AXUV100 photodiode at branch A (Fig. 2[link]). Photodiode measurements were adjusted to account for the reflectivity of the extra mirror used in branch A. Photodiode measurements were taken with an SSA size of 14000 µm horizontal (i.e. wide open) × 20 µm vertical (branch A), while intensity measurements were taken with an SSA of 25 µm × 25 µm (branch B). The difference in horizontal SSA size between branch A and branch B leads to a 100-fold increase in flux at the photodiode, which has been accounted for by scaling the photodiode flux values in Fig. 7 accordingly.

Measurements of the harmonic content of the beam, ζ, were taken using XPS of Au 4f7/2 peaks corresponding to different undulator harmonics. Through analysis of the relative peak areas, an estimate of ζ was calculated for fundamental photon energies from 130 eV to 350 eV in 20 eV steps,

[\zeta = {{\Phi_{n\,\gt\,1}} \over {\Phi_{n}}}, \eqno(10)]

where Φn>1 is the photon flux contained in higher undulator harmonics and [\sum_{n}\Phi_{n}] is the total flux.

5. Experimental results

5.1. Beam profile

The FWHM of the fundamental (n = 1) intensity obtained from the fit to the measured intensity distribution at the grating mask plane is shown in Fig. 5[link] for photon energies 90 eV to 270 eV. Comparison with simulated intensity FWHM shows a similar dependence on photon energy. However, at low photon energies the comparison between the two methods shows less agreement, which may be attributed to uncertainties resulting primarily from the low signal-to-noise for the intensity measurements taken with low photon energy.

5.2. Higher-harmonic content

Measurements were taken by XPS with Cff values of 2 and 1.4. The total ζ was found to significantly reduce at Cff = 1.4 for all photon energies measured as expected (Kleemann et al., 1997[Kleemann, B. H., Gatzke, J., Jung, C. & Nelles, B. (1997). Proc. SPIE, 3150, 137-147.]), with a reduction from ∼10% harmonic content at Cff = 2 to ∼2% at Cff = 1.4 for a fundamental photon energy of 185 eV. The results of the harmonic content measurements are shown in Fig. 6[link] for Cff = 1.4 and compared with the harmonic contamination given by the fitting of Gaussians to the direct intensity measurements. The XPS measurements show good agreement with the direct beam measurements when comparing total intensity, and the simulation results show close agreement to both measurements.

[Figure 6]
Figure 6
Higher harmonic content (ζ) measurements using direct beam measurements (black) and XPS measurements (red). A single calculation from simulation has been included at 184.76 eV (green).

5.3. Flux and intensity at the mask plane

The total flux as a function of photon energy is shown in Fig. 7[link] as measured by three different methods: the AXUV100 photodiode at branch A beamline (Fig. 2[link]), the fundamental fitting of the direct beam intensity measurements taken at branch B, and simulated wavefront propagation.

[Figure 7]
Figure 7
(a) The simulated and measured total flux (Φn) and (b) the total intensity over the area of a four-grating mask ([I_{n}^{\,{\rm{G}}}]) consisting of 50 µm × 50 µm gratings arranged symmetrically with a 50 µm gap in the centre (shown in inset, where the gratings are shown as black squares, the n = 1 harmonic in blue, and the n = 3 harmonic in green). Flux and intensity measurements were taken at branch B, while flux measurements were also taken on branch A using an AXUV100 photodiode. Cff = 1.4 for all measurements shown.

The intensity measured over a typical IL-grating area, [I^{\,{\rm{G}}}_{n}], described earlier and shown in Fig. 7[link], is also included for each photon energy.

The simulated results show significantly greater flux and intensity over the grating area compared with the direct intensity measurements; however, they show close agreement with the flux measurements using the photodiode. Fig. 7[link] also shows that the intensity-on-mask of n = 1 at 185 eV is [I^{\,{\rm{G}}}_{1}] ≃ 2 × 108 photons s−1 cm−2, which gives a dose-on-mask of just Dm ≃ 6 × 10−6 mJ s−1 cm−2. However, recall that the intensity measurements shown in Fig. 7[link] were made with `low flux' beamline settings, using small WBS and SSA sizes to avoid over-saturation of the detector. Using `high flux' settings, the intensity at the mask plane is expected to be closer to the simulated results shown in Fig. 3[link]. Hence, the results shown in Fig. 7[link] serve mainly to confirm that the simulation can accurately predict intensity at the mask plane.

To evaluate the beamline's suitability for IL, partially coherent wavefront propagation was used to calculate Φ1 and [I^{\,{\rm{G}}}_{1}] as a function of SSA size for a photon energy of 185 eV. The SSA height and width were kept equal and varied from 200 µm to 1500 µm in 100 µm steps. [I^{\,{\rm{G}}}_{1}] was found to reach a maximum of 5.85 × 1013 photons s−1 cm−2 for a 300 µm × 300 µm SSA. For the same representative grating with efficiency η = 6% illustrated in Fig. 7[link], this value of [I^{\,{\rm{G}}}_{1}] corresponds to a value of Dw equal to 5.46 mJ s−1 cm−2. For Dm = 11 mJ cm−2, the minimum exposure time is then 1.9 s. The coherence length at the mask plane for these settings is 137 µm, which would allow for the coherent illumination of a mask containing 45 µm gratings.

6. Conclusion

Partially coherent wavefront propagation simulations have been used to accurately model the SXR beamline at the Australian Synchrotron. The model was shown to give a quantitative representation of photon flux that agrees with direct experimental measurements using an AXUV100 photodiode. However, the simulated flux and the photodiode flux measurements did not show agreement with the direct intensity measurements. This could be due to an error in the thickness and density of the Ultralene filter used to attenuate the beam, which has been found to vary in thickness by 10% from the nominal value of 4 µm (Surowka et al., 2020[Surowka, A. D., Birarda, G., Szczerbowska-Boruchowska, M., Cestelli-Guidi, M., Ziomber-Lisiak, A. & Vaccari, L. (2020). Anal. Chim. Acta, 1103, 143-155.]). The discrepancy could also have come from errors in simulation due to the assumption of ideal beamline elements, reflectivities and diffraction efficiencies, or a possible misalignment of the SSA in branch B of the SXR beamline, causing a drop in flux at the detector plane.

A method for higher harmonic suppression by adjusting Cff of the PGM from 2 to 1.4 was found to reduce total harmonic contamination at the mask plane from ∼10% to ∼2%.

The beamline was shown to possess adequate power and coherence for IL at EUV/SXR wavelengths when operating with a WBS of 4 mm × 3 mm and an SSA of 300 µm × 300 µm, allowing for 1.9 s patterning of a 45 µm maximum grating size, which is compatible with that used at other EUV-IL facilities (Sahoo et al., 2023[Sahoo, A. K., Chen, P., Lin, C., Liu, R., Lin, B., Kao, T., Chiu, P., Huang, T., Lai, W., Wang, J., Lee, Y. & Kuan, C. (2023). Micro Nano Eng. 20, 100215.]). The exposure time is sufficiently short for maintaining the nanometre mechanical stability required for studying lithographic processes at high spatial resolution. Greater power is possible at the cost of coherence — and therefore grating size — with experimental evidence that the transverse coherence length exceeds 30 µm for any choice of SSA size.

Acknowledgements

Synchrotron experiments were carried out at the SXR beamline of the Australian Synchrotron, part of ANSTO. The authors would like to thank Trey Guest for his early contributions to beamline modeling and helpful discussions. Open access publishing was facilitated by La Trobe University, as part of the Wiley–La Trobe University agreement via the Council of Australian University Librarians.

Funding information

This research was supported by the Australian Research Council (ARC) Linkage Infrastructure, Equipment and Facilities (LIEF) scheme, which funded part of the equipment used in this work. Two of the authors (JBMK, BH) were supported by an Australian Government Research Training Program (RTP) Scholarship, and one (JBMK) by an AINSE Ltd Postgraduate Research Award (PGRA).

References

First citationBraig, C., Käsebier, T., Kley, E.-B. & Tünnermann, A. (2011). Opt. Express, 19, 14008–14017.  CrossRef CAS PubMed Google Scholar
First citationChubar, O. (2014). Proc. SPIE, 9209, 920907.  Google Scholar
First citationChubar, O. & Celestre, R. (2019). Opt. Express, 27, 28750–28759.  Web of Science CrossRef PubMed Google Scholar
First citationChubar, O., Chu, Y. S., Huang, X., Kalbfleisch, S., Yan, H., Shaftan, T., Wang, G., Cai, Y. Q., Suvorov, A., Fluerasu, A., Wiegart, L., Chen-Wiegart, Y., Thieme, J., Williams, G., Idir, M., Tanabe, T., Zschack, P. & Shen, Q. (2016). AIP Conf. Proc. 1741, 040002.  Google Scholar
First citationChubar, O. & Elleaume, P. (1998). Proceedings of the Sixth European Particle Accelerator Conference (EPAC'98), 22–26 June 1998, Stockholm, Sweden, pp. 1177–1179. THP01G.  Google Scholar
First citationDe Bisschop, P. (2017). J. Micro/Nanolith. MEMS MOEMS, 16, 041013.  Google Scholar
First citationEkinci, Y., Vockenhuber, M., Hojeij, M., Wang, L. & Mojarad, N. (2013). Proc. SPIE, 8679, 867910.  CrossRef Google Scholar
First citationEkinci, Y., Vockenhuber, M., Mojarad, N. & Fan, D. (2014). Proc. SPIE, 9048, 904804.  Google Scholar
First citationFomenkov, I. (2019). Proceedings of the 2019 ASML Source Workshop, 5 November 2019, Amsterdam, The Netherlands, p. S1.  Google Scholar
First citationGoldberg, K. A. & Mochi, I. (2010). J. Vac. Sci. Technol. B, 28, C6E1–C6E10.  CrossRef CAS Google Scholar
First citationHarada, T., Teranishi, N., Watanabe, T., Zhou, Q., Bogaerts, J. & Wang, X. (2019). Appl. Phys. Expr. 13, 016502.  CrossRef Google Scholar
First citationHe, A., Chubar, O., Dvorak, J., Hu, W., Hulbert, S., Walter, A., Wilkins, S. & Salort, C. (2020). Proc. SPIE, 11493, 1149307.  Google Scholar
First citationHecht, E. (1987). Optics, 2nd ed. Addison-Wesley.  Google Scholar
First citationKarim, W., Tschupp, S. A., Oezaslan, M., Schmidt, T. J., Gobrecht, J., van Bokhoven, J. A. & Ekinci, Y. (2015). Nanoscale, 7, 7386–7393.  CrossRef CAS PubMed Google Scholar
First citationKleemann, B. H., Gatzke, J., Jung, C. & Nelles, B. (1997). Proc. SPIE, 3150, 137–147.  CrossRef Google Scholar
First citationKnappett, J. (2021). Wavefront Propagation Simulation of X-ray Interference Nanolithography. Master's thesis, La Trobe University, Australia.  Google Scholar
First citationLaundy, D., Alcock, S. G., Alianelli, L., Sutter, J. P., Sawhney, K. J. S. & Chubar, O. (2014). Proc. SPIE, 9209, 920903.  Google Scholar
First citationLevinson, H. J. (2022). Jpn. J. Appl. Phys. 61, SD0803.  CrossRef Google Scholar
First citationMeng, X., Yu, H., Wang, Y., Ren, J., Xue, C., Yang, S., Guo, Z., Zhao, J., Wu, Y. & Tai, R. (2021). J. Synchrotron Rad. 28, 902–909.  Web of Science CrossRef IUCr Journals Google Scholar
First citationMojarad, N., Gobrecht, J. & Ekinci, Y. (2015a). Sci. Rep. 5, 9235.  CrossRef PubMed Google Scholar
First citationMojarad, N., Gobrecht, J. & Ekinci, Y. (2015b). Microelectron. Eng. 143, 55–63.  CrossRef CAS Google Scholar
First citationMojarad, N., Kazazis, D. & Ekinci, Y. (2021). J. Vac. Sci. Technol. B, 39, 042601.  CrossRef Google Scholar
First citationMojarad, N., Vockenhuber, M., Wang, L., Terhalle, B. & Ekinci, Y. (2013). Proc. SPIE, 8679, 867924.  CrossRef Google Scholar
First citationOtsuka, T., Li, B., O'Gorman, C., Cummins, T., Kilbane, D., Higashiguchi, T., Yugami, N., Jiang, W., Endo, A., Dunne, P. & O'Sullivan, G. (2012). Proc. SPIE, 8322, 832214.  CrossRef Google Scholar
First citationSahoo, A. K., Chen, P., Lin, C., Liu, R., Lin, B., Kao, T., Chiu, P., Huang, T., Lai, W., Wang, J., Lee, Y. & Kuan, C. (2023). Micro Nano Eng. 20, 100215.  CrossRef Google Scholar
First citationShannon, C. (1949). Proc. Inst. Radio Eng. (IRE), 37, 10–21.  Google Scholar
First citationSolak, H., David, C., Gobrecht, J., Golovkina, V., Cerrina, F., Kim, S. & Nealey, P. (2003). Microelectron. Eng. 67–68, 56–62.  CrossRef CAS Google Scholar
First citationSolak, H. H., David, C., Gobrecht, J., Wang, L. & Cerrina, F. (2002). Microelectron. Eng. 61–62, 77–82.  Web of Science CrossRef CAS Google Scholar
First citationSurowka, A. D., Birarda, G., Szczerbowska-Boruchowska, M., Cestelli-Guidi, M., Ziomber-Lisiak, A. & Vaccari, L. (2020). Anal. Chim. Acta, 1103, 143–155.  Web of Science CrossRef CAS PubMed Google Scholar
First citationUzoma, P. C., Shabbir, S., Hu, H., Okonkwo, P. C. & Penkov, O. V. (2021). Nanomaterials, 11, 2782.  CrossRef PubMed Google Scholar
First citationVan Schoot, J. (2021). Proc. SPIE, 11609, 1160905.  Google Scholar
First citationVieu, C., Carcenac, F., Pépin, A., Chen, Y., Mejias, M., Lebib, A., Manin-Ferlazzo, L., Couraud, L. & Launois, H. (2000). Appl. Surf. Sci. 164, 111–117.  CrossRef CAS Google Scholar
First citationWang, X., Kazazis, D., Tseng, L.-T., Robinson, A. P. & Ekinci, Y. (2021). Nanotechnology, 33, 065301.  CrossRef Google Scholar
First citationWootton, K. P., Boland, M. J., Dowd, R., Tan, Y. E., Cowie, B. C. C., Papaphilippou, Y., Taylor, G. N. & Rassool, R. P. (2012). Phys. Rev. Lett. 109, 194801.  CrossRef PubMed Google Scholar
First citationWootton, K. P., Rassool, R. P., Boland, M. J. & Cowie, B. C. C. (2013). Proceedings of the 2nd International Beam Instrumentation Conference (IBIC2013), 16–19 September 2013, Oxford, UK, pp. 546–549. TUPF19.  Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds