In atom beam lithography, a beam of kilo-electron-volt helium atoms illuminates a stencil mask and transmitted beamlets transfer the mask pattern to resist on a substrate. It shares the advantages of masked ion beam lithography but is immune to charging artifacts that limit resolution and pattern fidelity. This paper describes a high-brightness source of energetic He atoms, where He+ ions are extracted from a multicusp ion source, focused by two-stage accelerating optics, and neutralized by charge-transfer scattering in a differentially pumped, He-filled cell. Since scattering angles are extremely small, the straight line trajectories of scattered atoms are essentially tangent to the (possibly curved) trajectories of the parent ions. Space-charge repulsion prevents the ion beam crossing over; instead, it converges to a waist of minimum cross section before diverging further downstream. Atom trajectories produced by a cell placed in the region of intense space charge near the waist are strongly affected by the curvature of ion trajectories within the cell. The flaring of the ion beam due to space charge can be used to increase the width of the atom beam, although to the detriment of resolution. In this paper, the authors study a configuration where the cell is placed in the converging ion beam as far as practicable from the ion-beam waist. The atom beam then converges to a crossover, which becomes the virtual source seen by the mask. The source diameter and angular flux density initially increase with increasing cell pressure but saturate at higher pressures; the respective saturation values at 50 keV are 125 μm (2σ) and 8.7 × 1017 particles/s sr. Under these conditions, the beam diameter is ∼2.5 cm, 7 m from the source. A practical system for subnanometer printing is discussed with 0.2 nm (2σ) penumbral blur and 1.25 × 1013 particles/s cm2 flux density over a 1 cm circular field.

The helium ion microscope (HIM) has enabled the fabrication of a class of novel devices that rely on ∼30 keV He+ ion irradiation to directly modify substrates at a resolution far beyond the capability of resists. Examples include the first insulating tunnel junctions in YBa2Cu3O7 − δ (Ref. 1) by radiation damage to a single-crystal substrate and ultrafine graphene nanoribbons by surface ablation.2 Causer et al.3 have demonstrated chemical mixing in FePt3 at 15 keV without disrupting the crystal lattice, a critical development for applying direct substrate modification to nanomagnetic device technology. Moreover, molecular dynamic simulations suggest that ultralow energy, heavy ions (e.g., ∼60 eV Xe+) may enable the formation of semiconducting graphene nanoroads4 by selectively depopulating F-atoms on a fluographene sheet.5 The critical fluence required for these applications, 1016–1017 particle/cm2, is extremely high compared to resist-based approaches; this, coupled with a serial writing modality, severely limits the scale of circuits that can be fabricated by an HIM. Since beam current at the highest resolution (0.5 nm) is ∼1 pA, it would take ∼22 h (plus overhead) to write a 7 × 7 mm2 chip with 10% pattern density and 1 × 1016 particle/cm2 exposure fluence, and 220 h for 1 × 1017 particles/cm2. Moreover, the field size, 30 × 30 μm2 for subnanometer resolution, is extremely small due to a limited depth of field6 and field stitching capability is currently absent in HIM instruments. This paper describes a source of energetic He atoms for direct substrate modification through a stencil mask, a parallel exposure strategy with potentially much larger chip size and higher throughput than can be achieved by an HIM. Our goal is a source with the combination of virtual source size, angular flux density, and beam diameter required for system-level studies of these intriguing technologies.

We follow a classical approach7,8 where He+ ions are extracted from a plasma source and neutralized by charge-transfer scattering in a differentially pumped, high pressure cell. A downstream electrostatic deflector removes unscattered ions from the mixed beam emanating from the cell. The absence of space-charge repulsion in the neutral beam implies that the critical lithographic parameters, penumbral blur, flux density, and field size, depend only on the source—its diameter, angular flux density, and divergence angle-irrespective of the atomic particle flux through the system. It is, therefore, straightforward to extract the source parameters at any point along the beamline and to use that data to explore the design space of optimized systems for the applications discussed above.

Figure 1 illustrates the source concept. He+ ions are extracted from a decapole multicusp ion source,9 focused by two-stage accelerating optics, and neutralized by charge-transfer scattering10 in a high pressure, He-filled cell. An electrostatic deflector removes unscattered ions from the beam. The ions are accelerated to 3.1 keV in the first stage and to the final energy (typically 10–50 keV) in the second. Space charge, acting as a negative lens,11 prevents the ion beam from crossingover; instead, it converges to a waist of minimum cross section12 and diverges downstream. The first stage uses a repurposed lens from a Coates and Welter Cwikscan microscope; the second uses a three-element lens from a Varian ion implanter. The intermediate energy was chosen empirically to optimize ion beam uniformity downstream of the beam waist.

FIG. 1.

Neutral particle source schematic. Details are described in the text. To aid the eye, electrodes with like potentials have the same color.

FIG. 1.

Neutral particle source schematic. Details are described in the text. To aid the eye, electrodes with like potentials have the same color.

Close modal

Since charge-transfer scattering angles are extremely small (<0.2 mrad at 50 keV), atoms leave the point of neutralization on straight line trajectories that are essentially tangent to the, possibly curved, trajectories of the parent ions. These trajectories crossover, becoming the virtual source of the atom beam for downstream applications. The cell is positioned in the converging ion beam just downstream of the final electrode of the second stage lens. This minimizes space charge in the cell, which would cause the ion trajectories to curve, blurring the virtual atom source.

The charge-transfer cell is a 2.54 cm long cylinder with 3 mm circular differential pumping apertures. The flow of He gas to the cell is controlled by a needle valve. The pressure of the source chamber, measured by an ion gauge, was converted to cell pressure using the conductance of the apertures and the speed (1200 l/s) of the pumps.

The Faraday cup provides a near-absolute determination of flux density for singly charged ions. In lithography, they are used to determine the relationship between the thickness of resist remaining (TRR), after exposure and development, on ion fluence. Typically, this information is displayed as a Hurter and Driffield (H&D) curve,13 normalized TRR versus log10(fluence), where fluence is the integrated flux density over the duration of the exposure. This curve provides the important resist parameters including critical dose and contrast.

Faraday cups have also been used to monitor neutral atom flux density.14 In the experiments reported here, atoms pass through an aperture in a grounded enclosure, traverse a cylindrical electrode and impinge on a collector plate.15 The collector is connected to the input, a virtual ground, of a Keithley 414A picoammeter. For He+ ions, the cylindrical electrode is negatively biased (−90 V) to return secondary electrons to the collector. For atoms, it is positively biased (+90 V) to attract secondary electrons. Then, collector current Ic is equal to the secondary electron current generated by atom impact on the collector; Ic = eAγf, where e, A, γ, and f, are, respectively, the elementary charge, the area of the cup aperture, the secondary electron yield of the collector, and the atom flux density. To quantify the atom flux requires an accurate determination of γ. While γ is well known for clean metals, its value for contaminated surfaces can be highly variable. The collector of the present cup, which has not been cleaned in over 25 years, is coated by a brown, carbonaceous film due to ion/atom beam polymerization of residual hydrocarbons in the vacuum systems (2.0–5.0 × 10−7 Torr base pressure). It is likely that this coating has been converted to a terminal state of maximum graphitic order16 by the radiation. A method for measuring γ is discussed in the paragraph below.

Resist exposure processes are the same for He+ ions and He atoms of the same energy because the charge state of the incident particles converges within a few atomic layers to a value that depends only on the resist material.17,18 It follows that the fluence F required for a given value of normalized TRR is the same for ions and atoms. For neutral atoms, F = C/eAγ, where C is the secondary electron charge produced by neutral particle fluence F. Let T(log10(F)) be normalized TRR as a function of log10(F), the function represented by the H&D curve. For neutral particles, T(log10(F)) = T(log10(C/eAγ)) = T(log10(C/eA) − log10(γ)). Thus, the normalized TRR for neutral particles is T(log10(C/eA)), shifted by log10(γ).

Figure 2 shows the H&D curve for 50 keV He ions and normalized TRR versus log10(C/eA) for 50 keV He atoms for a plasma deposited resist.19 Best alignment occurs for γ = 0.99. The curves for γ = 0.94 and 1.04, which bound both the ion and best fit atom curves, illustrate the sensitivity of the atom curves to the value of γ. Lines between the data points are meant only to guide the eye. Data for atom fluences below 1.5 × 1014 particles/cm2 have been omitted due to systematic errors related to develop swelling of partially cross-linked films. We assumed for simplicity that γ is unity for the measurements of atom flux density and fluence reported in this paper.

FIG. 2.

Normalized resist thickness remaining after development as a function of log10(fluence) for 50 keV ions and atoms for secondary electron yields of 0.94, 0.99, and 1.04.

FIG. 2.

Normalized resist thickness remaining after development as a function of log10(fluence) for 50 keV ions and atoms for secondary electron yields of 0.94, 0.99, and 1.04.

Close modal

The flux density f of the ion/atom beam was measured using a scanning Faraday cup, with a 3 × 0.100 mm2 aperture, 7.0 m from the atom beam crossover. The angular flux density is given by fΩ=f/Ω, where Ω is the solid angle subtended by the cup aperture when viewed from the source.

The size and position of the crossover were determined by a pair of knife-edge experiments.20 Briefly, a Faraday cup is first used to profile the penumbral shadow of a razor blade at position x, assumed to be downstream from the virtual source. A cumulative Gaussian distribution is fitted to the profile and the standard deviation σp extracted. The diameter 2σv of the atom source is then given by 2σv=((xx0)/g)2σp, where x0 is the unknown position of the source and g is the distance from the knife-edge to the Faraday cup. A second measurement, at a different knife-edge position, enables 2σv and x0 to be determined.

The source size versus pressure data were smoothed for the calculation of brightness using a least square fit to a first order response function [σv=a(1ebP)+c, where a, b, c are the fitting parameters]. The standard deviations21 of these fits are 0.04 and 0.1 μm at 25 and 50 keV, respectively. The smoothed source size and the raw data for angular flux density were used to derive the brightness β=fΩ/2πσv2.22 

Figure 3 shows angular flux density versus Faraday cup position for 25 and 50 keV ion and atom beams. The charge-transfer cell pressure was 13 mTorr. The beam is wider for ions than for atoms due to space-charge repulsion near the beam waist. The charge-transfer efficiency, the ratio of total atom flux to total ion flux, is 85% at 25 keV and 65% at 50 keV, probably reflecting the difference in the total charge-transfer cross sections at these energies. The width of the uniform region in the center of the atom profiles is ∼2.5 cm for both energies; it is ∼2.7 cm, full-width-at-half-maximum. The angular flux density, source size (2σ), and brightness at 50 keV are 8.7 × 1017 particles/s sr, 124 μm, and 3.75 × 1021 particles/s sr cm2, respectively. At 25 keV, the angular flux density and source size are 3.7 × 1017 particles/s sr and 88 μm, respectively.

FIG. 3.

Angular flux density vs Faraday cup position 7.0 m downstream of the atom source for 25 and 50 keV He ion and atom beams. Charge-transfer cell pressures were 13 mTorr. The connecting lines between the data points are meant to guide the eye.

FIG. 3.

Angular flux density vs Faraday cup position 7.0 m downstream of the atom source for 25 and 50 keV He ion and atom beams. Charge-transfer cell pressures were 13 mTorr. The connecting lines between the data points are meant to guide the eye.

Close modal

Figure 4 shows (a) source size (2σv), (b) angular flux density fΩ at beam center, and (c) brightness as functions of cell pressure P for 25 and 50 keV helium atom beams. Angular flux density increases with increasing pressure due to the increased neutralization probability within the cell; it saturates at high cup pressure due to the establishment of a steady state between the neutralization and reionization processes. The increase in source size with increasing pressure is probably due to direct scattering of He-on-He;10 the reason for the saturation remains unclear. Note that brightness, which is included for completeness, is not proportional to the beam energy, as would be expected for ion beams.

FIG. 4.

(a) Source size, (b) angular flux density, and (c) brightness as functions of charge-transfer cell pressure at 25 and 50 keV. The trend lines though the source size data are least square fits to a first order response function. The logarithmic trend lines through the brightness data are meant only to guide the eye.

FIG. 4.

(a) Source size, (b) angular flux density, and (c) brightness as functions of charge-transfer cell pressure at 25 and 50 keV. The trend lines though the source size data are least square fits to a first order response function. The logarithmic trend lines through the brightness data are meant only to guide the eye.

Close modal

We assume for this model a 50 keV beam to minimize lateral scattering in the near-surface region of the substrate and a charge-transfer cell pressure of 13 mTorr to maximize charge-transfer efficiency. These are the same conditions used to obtain the data discussed in Sec. III A (beam profiles). Accordingly, penumbral blur 2σp, flux density f, and field size S for a short system with source-to-mask distance L and mask-to-substrate gap g are given, respectively, as follows:2σp=(g/L)2σv,f=fΩ/L2,S=S7(L/L7), where 2σv=124μm, fΩ=8.7×1017particles/ssr, S7=2.5cm, and L7=7m. We choose a gap g=5μm to minimize penumbral blur and limit diffraction. Though small, this gap can be easily maintained by clamping mica spacers between the mask and substrate—this is also an easy way to eliminate relative motion during exposure. Diffraction, which could become important in the low nanometer range of feature sizes with this gap, is poorly understood in this setting and will not be included in the model. Finally, exposure time t=F/f, where F is the fluence required to modify the substrate.1,2Figure 5 shows flux density and penumbral blur as functions of source-to-mask distance L. The scale at the top of the figure shows field size corresponding to the length scale below the figure. This density is the lower bound of the range 10161017particles/cm2, discussed in the introduction for direct substrate modification. Under these conditions, L = 2.8 m yields penumbral blur, exposure time, and field size of 0.2 nm, 14 min, and 1 cm, respectively. The other end of the range, 1017particles/cm2, L = 1.4 m, yields corresponding values of 0.5 nm, 35 min, and 0.5 cm.

FIG. 5.

Penumbral blur (2σ) and flux density as functions of distance from the source position for conditions described in the text.

FIG. 5.

Penumbral blur (2σ) and flux density as functions of distance from the source position for conditions described in the text.

Close modal

We have designed and characterized a 50 keV source of energetic He atoms for proximity lithography and explored its applicability to nanosystem prototyping based on resistless, direct substrate modification. We find that it can access the range of exposure densities up to 1017 particles/cm2 with penumbral blur below 0.5 nm, field diameters above 0.5 cm, and exposure times below 35 min; for 1016 particles/cm2, the corresponding values are 0.2 nm, 1 cm, and 14 min, respectively. While far from VLSI standards, throughput seems large enough to enable system-level studies of high-temperature superconducting and nanomagnetic device technologies. This paper addresses only the source issue for nanosystem prototyping by direct substrate modification. The fabrication of masks with the required resolution, line-edge roughness, and pattern density, as well as the potential impact of diffraction remains as a serious open issue. There is, however, a substantial body of knowledge on mask fabrication for integrated circuit manufacturing23 that can serve as a starting point for nanoscale mask development.

The research was supported by the National Science Foundation (NSF) under Award Nos. DMI-0521523 and ECS-0404308, the Texas Center for Superconductivity at the University of Houston, and the Cullen Foundation. The authors are grateful to Hongjie Guo for his early contributions to this project and to Michael R. Young for many helpful discussions and outstanding machine shop support over many years. The authors would also like to thank John Notte and Ray Hill of Carl Zeiss Microscopy for useful discussions of field size in the helium ion microscope and Horst Rogalla (University of Colorado, Boulder) for sharing his insight into the lithography requirements of nanosystem prototyping.

1.
Shane A.
Cybart
 et al,
Nat. Nanotechnol.
10
,
598
(
2015
).
2.
Ahmad N.
Abbas
 et al,
ACS Nano
8
,
1538
(
2014
).
3.
Grace L.
Causer
 et al,
Appl. Surf. Sci.
459
,
672
(
2018
).
4.
M. A.
Ribas
 et al,
Nano. Res.
4
,
143
(
2011
).
5.
O.
Lehtinen
 et al,
Phys. Rev. B
81
,
153401
(
2010
).
6.
R.
Flatabø
 et al,
Nanotechnology
29
,
275301
(
2018
).
7.
I.
Amdur
,
J. E.
Jordan
, and
S. O.
Colgate
,
J. Chem. Phys.
34
,
1525
(
1961
).
8.
D. J.
Rej
,
I.
Henins
,
R. I.
Fonck
, and
Y. J.
Kim
,
Rev. Sci. Instrum.
63
,
4934
(
1992
).
9.
C. M.
McKenna
,
J. E.
Wood
,
J. L.
Bartelt
,
R. D.
Olney
,
J. W.
Ward
, and
C. W.
Slayman
, U.S. patent 4,757,208 (12 July
1988
).
10.
J. H.
Newman
 et al,
J. Geophys. Res.
91
,
8947
(
1986
).
11.
D.
Gabor
,
The Electron Microscope
(
Hulton
,
London
,
1946
), p.
73
et seq.
12.
J. W.
Schwartz
,
RCA Rev.
18
(
1
),
3
(
1957
).
13.
F.
Hurter
and
V. C.
Driffield
,
J. Chem. Technol. Biotechnol.
9
,
455
(
1890
).
14.
D. W.
Fahey
,
L. D.
Schearer
, and
W. F.
Parks
,
Rev. Sci. Instrum.
49
,
503
(
1978
).
15.
C. W.
Slayman
,
J. L.
Bartelt,
C. M.
McKenna
, and
J. Y.
Chen
,
Opt. Eng.
22
,
208
(
1983
).
16.
T.
Venkatesan
,
R. C.
Dynes
,
B.
Wilkens
,
A. E.
White
,
J. M.
Gibson
, and
R.
Hamm
,
Nucl. Instrum. Methods B
1
,
599
(
1984
).
17.
J. F.
Ziegler
,
J. P.
Biersack
, and
M. D.
Ziegler
,
Nucl. Instrum. Methods Phys. Res. B
268
,
1818
(
2010
).
18.
M.
Gruntman
,
Rev. Sci. Instrum.
68
,
3617
(
1997
).
19.
Dhara
Parikh
,
Barry
Craver
,
Hatem
Nounu
,
Fu-on
Fong
, and
John C.
Wolfe
,
J. Microelectromech.
17
,
735
(
2008
).
20.
B. P.
Craver
and
J. C.
Wolfe
,
J. Phys. D Appl. Phys.
41
,
024007
(
2008
).
21.
NIST/SEMATECH e-Handbook of Statistical Methods, see http://www.itl.nist.gov/div898/handbook/.
22.
23.
Rainer
Kaesmaier
,
Hans
Löschner
,
Gerhard
Stengl
,
John C.
Wolfe
, and
Paul
Ruchhoeft
,
J. Vac. Sci. Technol. B
17
,
3091
(
1999
).