1.IntroductionAs the size of integrated circuit elements approaches the molecular scale, the photoresist films used to form those elements can no longer be approximated as a continuum material. This is particularly true for organic and hybrid chemically amplified resists exposed using extreme ultraviolet (EUV) radiation. A few nanometer-scale printed features will contain only a small number of chemical moieties, including a photoacid generator (PAG), base quencher (Q), and the base polymer. When combined with the small number of EUV photons typically employed to create the feature pattern in the film, these moieties lead to statistical variations of the chemical composition within the exposed and post-expose baked feature that in turn cause readily apparent deviations from the intended ideal shape. Two ways in which this statistical behavior is displayed are line-edge roughness (LER),1 which measures how the edge of a printed feature deviates from an ideal straight line, and stochastic printing failures,2 where a feature is missing from an extended array in a rare, random way. Consideration of the factors that control image quality in chemically amplified (CA) DUV photoresists identified a resolution—line edge roughness—sensitivity (R-L-S) tradeoff.3 Because of the statistical variations described above, it is not obvious that how the tradeoff operates is the same in EUV lithography. The use of ionizing radiation EUV resists adds a significant level of complexity and an additional source of inhomogeneity to the chemical composition of the exposed resist. The small number of EUV photons required to print a feature leads to a variation in the number of photons absorbed from one feature to the next. This variation, shot noise, follows Poisson statistics and leads to a significant difference in the extent of photolysis within a collection of features. The resulting formation of highly reactive species in the films such as positive and negative ions and free radicals is likely to be an additional factor that degrades resolution and LER but is not necessarily directly connected to sensitivity in a simple way. It is evident that improved control over resist composition variations and the attendant radiation chemistry may offer a means of reducing LER and improving resolution while maintaining a desired sensitivity. There are several aspects to consider for the current organic positive tone resist families. First, the components of an amorphous resist film are randomly located within the solid solution and their position will vary according to Poisson statistics. This aspect of the resist structure convolves with shot noise during EUV exposure. The extent of conversion of the matrix polymer to a soluble form depends on the kinetics of the acid-catalyzed deprotection reaction, which depends in turn on the details of the local radiation-induced chemical reactions, especially ones that generate acid. Second, aggregation of PAGs and Q compounds the complexity of the chemistry. Third, the polymer characteristics add another dimension of statistical variability. For a typical resist polymer prepared by free radical polymerization of two or more monomer species, the final product will contain a statistical distribution of chain lengths and compositions. For positive tone systems, the solubilization of an individual polymer chain in an aqueous base solution depends on the removal of hydrophobic groups by acid-catalyzed chemistry until a threshold of the ionizable deprotection products is reached.4 With chains of different lengths, composition, and local acid catalysis conditions, each polymer molecule will have a different threshold where solubility is attained. Several approaches to improve photoresist materials homogeneity have been pursued, for example, the use of polymerization techniques that provide narrow molecular weight distributions,5–7 introduction of polymer-bound PAG and quencher,8,9 use of molecular glasses,10 and more recently defined-sequence polymers (DSP) with uniform lengths. Synthetic strategies to prepare DSP materials using methods such as reversible deactivation radical polymerization and solid-phase synthesis have been reviewed.11 Preparation and applications of peptoid-based DSP materials to hierarchical nanomaterials for use in biomedical systems have also been described.12 The potential for polypeptides to improve control over the photoresist film for DUV and EUV exposures has recently been brought forward in studies describing their synthesis and imaging.13,14 DSP-based resist materials offer a way to mitigate polymer-level imperfections through precise placement of resist moieties and PAGs. These structural characteristics can potentially provide further benefits within the amorphous photoresist films themselves if they influence the degree of randomness of photoresist components. A reduction in the random distribution of reactive moieties within the films may provide a way to control the details of inevitable radiation-induced chemical reactions occurring during EUV exposure. In this study, we investigate how imaging of a series of model CA resist films made of DSP materials compares with imaging of random copolymers of the same composition using a computational polymer packing model combined with a simple exposure-deprotection model. The simulations provide some initial steps toward an understanding of the characteristics of DSP-based CA photoresist films and the developable images they form within an amorphous layer. 2.Materials and MethodsThe goal of this study is to calibrate thinking about how controlled polymer sequence and length, resist component aggregation, and polymer-bound resist components affect the overall spatial homogeneity of an amorphous CA-positive tone photoresist film and image quality. Similar considerations will apply to negative tone photoresists. Our primary metric for homogeneity is the spacing between the various monomers, so we employ a coarse-grained model of the resist film with a simplified representation of polymer chains as linear strings of spherical monomers. Monomers are connected by bonds with fixed bond lengths and bond angles, while the torsional angles along the chain are unconstrained. An ensemble of chains is randomly packed to form a dense amorphous glassy film. This representation of the polymers enables large film volumes to be calculated to improve precision of the statistical analysis of monomer distributions. The volume dimensions used in this work are . 2.1.Polymer Sequences and Film CompositionsTo evaluate the influence of polymer chain structure and location of PAG and Q on imaging, we examine several distinct types of polymer sequences within photoresist films of the same composition. Two are random copolymers where the PAG and Q are either bound or unbound to the photoresist polymer. The third is a series of DSPs with bound PAG and Q and one DSP with unbound PAG and Q. We also include a nonpolymeric case, which is a random distribution of monomers unconstrained by chain links. The base CA photoresist film composition selected for evaluation is representative of a conventional positive tone EUV resist of the ESCAP-class15 and has relatively well-understood chemistry stemming from deep UV16 and EUV photolithography studies.17 It includes a copolymer of 4-hydroxystyrene (HOST), tert-butyl methacrylate(TBMA), and styrene (STYR). The composition also includes a PAG and a Q, with both species present either in polymer-bound form or as free molecules. In these calculations, PAG and Q are not chemically structured, only represented as spherical monomers. Their function is assumed to be similar to that of a perfluorosulfonic acid and an amine for the purposes of estimating developable images as described in Sec. 2.2.2 below. PAG and Q in free molecule form may be random distributions of individual molecules or randomly sized aggregates with a mean cluster size of three molecules. In all cases, the mole ratios are HOST:TBMA:STYR:PAG:Q = 50:30:6.7:10:3.3. Table 1 summarizes the model compositions examined in this study. Table 1Combinations of polymer-PAG-quencher, with variations in sequence type, chain length, and polydispersity.
Table 2 presents the DSP sequences utilized in this work. Each is a 30-mer with distinct placements of bound PAG and quencher. The 90-mers are three repetitions of the specified 30-mer sequence. The sequences are not exhaustive but are designed to assess the effects of proximity on a chain. In each 30-mer, the first and third sets of 10 are identical, while the second, which contains a quencher but no styrene, differs from them. Sequence A serves as a reference, and results from it are the main ones reported here. Sequences B-E vary the position of the quencher relative to a PAG and of the TBMA relative to the PAG in several combinations. It was found that in general, the statistical analysis results (see Sec. 2.2) for polymer film structures using sequences B-E in the present work are nearly identical to those for sequence A, except for PAG-Q nearest neighbor distances. Table 2Defined sequence polymer types.
2.2.Computational Methods2.2.1.Polymer packing and film formation-polyscopeIn this work, we introduce Polyscope, an open-source code18 to construct films from polymer chains based on the polymer chain growth algorithm of Mueller et al. (PolyGrow and Embed, Ph.D Thesis, Chapters 3 and 4).19 PolyGrow treats chains as simple ball-and-stick assemblies, and builds each chain by choosing a random starting point in the film volume and adding monomers one by one with a randomly chosen torsional angle while avoiding overlap with already-placed chains. For this, monomers are treated as spheres with a user-specified amount of allowable overlap during packing. As the film volume becomes densely packed, the procedure performs a local relaxation of chains near the chain end under construction to facilitate addition of more monomers. Polyscope implements this chain growth and packing algorithm to form a representative photoresist film slice to which periodic boundary conditions are applied on all faces. Representative film slices constructed in this work are shown in Fig. 1. Polyscope calculates the spatial distribution of polymer moieties within the film slice as a function of user-supplied key parameters, including polymer chain length, polydispersity, and random or defined sequencing. The code package enables detailed control over the film’s chemical composition. An arbitrary number of different types of chemical building blocks can be specified. The building blocks can be monomer molecules of varying structures from which polymer chains are constructed and additive molecules such as PAG and Q. The sequence of building blocks in the polymer chains is specified to be random or a specific fixed sequence that is replicated in every chain. For random copolymers, all reactivity ratios (i.e., probability of a monomer attaching to a particular chain end building block) are assumed to be unity. PAG and Q can be specified to be part of the polymer chain (bound state) or dissolved in the photoresist film as an additive. Additives may be placed in the film as a random distribution of single molecules or aggregates of user-specified cluster size. Aside from aggregate formation, we do not include in the model other potential interactions between building blocks, such as hydrogen bonding or stacking. At present, the roles such interactions might play in influencing the structure of spin-cast films are not well characterized and are outside the scope of the present study. A typical film slice simulation contains 30,000 building block units, each 1 nm in diameter.20 These building block units are apportioned as monomers and additives of different chemical compositions according to the user’s specifications. The packing density is 1 building block unit per cubic nm. The center-to-center distance of adjacent monomers in a chain is set to be 0.9 nm. The bond angle between adjacent monomers is 70.5 deg. During packing, a maximum 10% overlap of building block units is allowed, so the smallest possible center-to-center distance between building block units is 0.9 nm. After the construction of a film slice is complete, PolyScope outputs a text file with the spatial location and chemical type of each monomer in each polymer chain. Each monomer is annotated with an index that marks which chain contains the monomer. Analogous location and chemical information for each nonpolymeric additive molecule are also contained in the text file. For a typical packing simulation (chain length =30 monomers, polydispersity =1.5), we calculate an average radius of gyration and an average end-to-end distance of 7.0 nm. These values are comparable to those obtained by Patsis et al. for molecular model calculations of a linear resist polymer of 30 monomers with monomer dimensions of 1 nm (2.7 and 6.3 nm, respectively).21,22 2.2.2.Statistical analysisStatistical analyses of the simulation outputs are performed using the program , an open-source software environment for statistical computing and graphics (version 4.3.2)23 with the RStudio integrated development environment (version 2023.12.1.402).24 The spatstat library for (version 3.0–7) provides three-dimensional spatial statistical functions.25 Additional statistical analyses (histograms, median values) are performed using the QtiPlot data analysis tool (version 1.2.1).26 A key metric of interest is the films’ compositional uniformity. Figure 1 does not reveal any visibly obvious differences, so to assess this, we use the pair correlation function .25 The function counts the average number of neighbors falling on the surface of a sphere of radius r for all chemical building blocks. This provides a measure of the probability of observing a chemical building block of a specific type separated by a distance from a reference chemical building block, divided by the corresponding probability for a random Poisson spatial distribution. Hence, is sometimes described as indicating “excess probability”; a value greater than 1.0 indicates clustering while a value less than 1.0 indicates a degree of regularity or non-randomness. 2.2.3.Modeling film exposure and deprotectionTo examine the effects of EUV exposure on resist image spatial statistics using these slices, we apply a simple blanket exposure algorithm to the film compositions described above. For each composition, an energy dose is applied that in the end produces an overall partial deprotection of 50% of the TBMA-protecting groups in the film slice. The 50% value was chosen to represent the extent of deprotection at the resist image line edge. This is where the film undergoes an abrupt change in dissolution rate and where stochastic effects are manifested. The general exposure/deprotection protocol is as follows. To model EUV exposure, a photon absorption site is randomly selected in the film. A volume centered on the absorption site is defined, with a radius equal to the electron blur length . It should be noted that this volume does not necessarily contain a single polymer chain, but parts of multiple chains. Within that volume, there is a number of PAG molecules given by the PAG concentration. An average of PAGs (where is the quantum yield, acid molecules produced/photon absorbed) is randomly selected within the volume and converted to photogenerated acid. The absorption site and PAG selection steps are repeated until a specified number of acid molecules sufficient to produce 50% deprotection is produced. Following the exposure process, a volume centered on a randomly selected photogenerated acid with a radius equal to the acid diffusion length is defined. Within this volume, all chemical building blocks are selected in random order and two types of reaction steps are allowed until no further reaction steps can occur. If the selected chemical building block is a TBMA with a protecting group, then the TBMA is converted to its deprotected product methacrylic acid (MAA). This step is allowed as long as there is TBMA remaining in the volume. If the selected chemical building block is a Q, then the acid is neutralized and no further deprotection can take place. The acid selection and reaction step processes are repeated until every photogenerated acid has been addressed. For this study, we use an electron blur value of 4.0 nm27,28 defining a volume that contains building blocks. We assume a value of 4.0 that is each photon absorbed by a PAG produces on average four acids within the blur volume. We examine two acid diffusion lengths for the deprotection step, with radius values of 2.0 and 5.0 nm. 3.Results and DiscussionThe film slice compositions shown in Tables 1 and 2 are analyzed to compare the spatial distributions of key components—TBMA with acid-labile protecting groups, PAG and Q—within a random and a DSP of the same composition. The simple exposure and deprotection steps are then used to evaluate the impact of these spatial distributions on formation of a developable image. 3.1.Spatial Distribution of PAGOne might anticipate that the spatial distribution of groups deprotected by acid would depend in part on the spatial distribution of PAG in the photoresist. Figure 2 illustrates typical pair correlation functions for a selected range of different model compositions, where each composition contains identical numbers of the different chemical building blocks but in increasingly structured form. Figure 2(a) also includes, for reference, obtained for a random point cloud distribution of discrete points (black dotted lines), showing a completely uniform distribution, as well as for the case where the building blocks are present in the correct proportion but in non-polymerized form, the composition NonPoly in Table 1 (dark blue curve). The random non-polymerized PAG-PAG function shows low values inside a radius of 0.9 to 1 nm. This is a consequence of excluded volume—the centers of two building blocks cannot be closer than the sum of their radii (minus the allowed 10% overlap). Above 1 nm radius, the curve overlays the point cloud curve, with some periodic structure that arises because the finite size of the building blocks excludes some of the center-to-center distances. When the monomer building blocks (HOST, TBMA, and STYR) are polymerized as a random copolymer, and the additive building blocks (PAG and Q) are present as free individual molecules (the composition R-FP-FQ-90-1.0 in Table 1), the function (dark cyan curve) is little changed from the non-polymerized case. The orange curve depicts the function that results when the free PAG molecules are aggregated to a mean value of three molecules/cluster (composition R-AP-FQ-90-1.0). As expected, the probability that another PAG is the nearest neighbor (at a distance) is high compared to the unaggregated case. If the PAG and Q additives are polymer-bound, i.e., incorporated into the random copolymer (composition R-BP-BQ-90-1.0), the function (green curve) is very similar to both the nonpolymeric composition and the case where the additives are free individual molecules. These trends are in contrast to that for a polymer with a defined sequence (composition SA-BP-BQ-90-1.0 in Table 1), which produces the function shown as a violet line. The DSP pair correlation function at distances less than a 3-nm radius falls below that of a random distribution, indicating local order. All defined sequences with bound PAG and bound Q listed in Table 2 yield essentially identical functions. However, a DSP composition where PAG and Q are free molecules (composition SF-FP-FQ-90-1.0), gives a function (dashed tan curve) equivalent to the random polymer cases. Figure 2(b) shows functions for the same cases as Fig. 2(a) with a mean chain length of 30 monomers and a polydispersity of 1.5. The PAG-PAG pair correlation functions in Figs. 2(a) and 2(b) are quite similar to each other, indicating that chain length and polydispersity do not play a role in determining them. It has been implied that PAG incorporation in a resist polymer improves the photoresist’s spatial uniformity.29,30 Our calculations show that incorporation of the PAG functionality into a random polymer chain (a polymer-bound PAG) as described in Table 1 does not intrinsically impart a more regular spatial distribution of PAG functionalities in the film. While binding PAG to a polymer likely reduces PAG aggregation/segregation compared with a conventional formulation that contains a poorly soluble free PAG,31 our calculations indicate that there is no improvement to the uniformity of the PAG spatial distribution if it is polymer-bound compared with the case where the PAG is fully soluble (unaggregated) in the solid polymer. In the DSP case, the traces in Fig. 2 show that the probability of one PAG having another PAG as the nearest neighbor is decreased compared to random polymers with free and bound PAGs. In the random case, all of a PAG’s nearest neighbor sites have an equal probability of containing another PAG. In the defined sequence case, however, the monomers vicinal to a PAG along the chain are strictly excluded by design from being another PAG, decreasing the overall likelihood of a PAG having another PAG as the nearest neighbor. Presumably, the ideal spatial distribution in the film of each of the resist components would be regularly spaced, which would minimize contributions of compositional inhomogeneity to resist stochastics. We can examine how closely the various compositions tested approach ideality by assessing the spacing between chemical building blocks of a given type. Figure 3 shows histograms of PAG-PAG nearest neighbor distances for the cases in Fig. 2(a). Figures 3(a) and 3(b) show that for a random copolymer with bound and free PAG, the distance distributions are nearly identical. A large proportion of the PAG population (36% for bound and 35% for unbound) is within 1 nm of another PAG. The primary difference is that in the bound PAG case [Fig. 3(a)], there is a spike at 1.46 nm due to geminal PAGs in the chain, that is, PAGs geminally bonded to the same building block. Figure 3(c) shows the histogram for the defined sequence case with free PAG and free Q. Its form is similar to the random polymer cases; 32% of the PAG population is within 1 nm of another PAG. Figure 3(d) shows the histogram for the defined sequence case with bound PAG and bound Q, which shows that there is still a significant proportion (20%) of PAGs within 1 nm of another PAG, smaller than the random copolymer cases. Figure 3(e) displays the histogram for the case where free PAG is aggregated within a random copolymer film; in this case, essentially all PAG molecules (94%) have another PAG within 1 nm. For comparison, Fig. 3(f) shows histograms for near-ideal spatial distributions of PAG where regular, uniform spacing is enforced. For the formulation used in this work (3000 PAGs distributed in a 30,000 nm3 volume), we calculate that precisely regularly spaced PAG molecules would be 2.15 nm apart. As shown in the coral curve of Fig. 3(f), when the positions of the PAG molecules are allowed to deviate randomly from the ideal uniform spacing coordinates by a small amount (), a nearest-neighbor distribution centered on the precise value is produced. In the purple curve in Fig. 3(f), those positions are allowed to deviate randomly by a larger amount (), which not only broadens the nearest-neighbor distribution but shifts it to smaller separation distances because of the significant proportion of PAGs that are closer to each other than in the uniform spacing case. A comparison of Fig. 3(f) to Figs. 3(a)–3(d) highlights the important finding from this work concerning the impact that chain packing alone has on the spatial distribution of resist components. The simulations show that the random chain packing that is characteristic of glassy amorphous polymer films leads to random, nonuniform spatial distributions of components, with a significant fraction of the component population located adjacent to another of the same component. This is found for both random sequence and defined sequence polymers and indicates that it is highly unlikely that a way could be found to create polymers that can confer significant spatial uniformity within the amorphous structure. 3.2.Spatial Distribution of Deprotectable TBMAAn analysis of the spatial distribution of TBMA, similar to that performed for the PAGs, leads to similar conclusions. Figure 4(a) illustrates pair correlation functions for the five model compositions examined in Fig. 2(a). All cases with random copolymers yield essentially identical functions that are consistent with a random spatial distribution of TBMA in the film volume. Both defined sequence cases deviate from the random cases in a manner analogous to that found for the bound PAG: the probability of one TBMA unit having another TBMA unit as a nearest neighbor is decreased compared to the other cases. All defined sequences in the present study (Table 2) were designed to avoid placement of TBMA monomers consecutively on the polymer chain, so this result is not surprising for the reasons discussed above for the PAG spatial distribution. The TBMA-TBMA pair correlation functions for cases where the mean chain length is 30 monomers with a polydispersity of 1.5, shown in Fig. 4(b), are virtually identical to those in Fig. 4(a), again demonstrating that chain length and polydispersity do not play a significant role in determining them. Figure 5 shows histograms of the nearest TBMA-TBMA neighbor distance for the cases in Fig. 4(a). As found in Fig. 3 for the PAG distributions, a characteristic of a random spatial distribution is that a significant portion of the TBMA population has another TBMA immediately adjacent to it. In all cases where the polymer sequence is random, 74% of the TBMA has another TBMA within 1 nm distance. With the defined sequence cases, where TBMA is not immediately adjacent to another TBMA on the chain, this portion is reduced to 50% (free PAG/Q case) and 54% (bound PAG/Q case). 3.3.Spatial Correlations of PAG-Quencher PlacementsThe proximities of PAGs and Q in various film compositions may influence the imaging characteristics of the compositions, especially photoresist sensitivity and the clumpiness of deprotected polymer moieties. Figure 6 shows histograms of nearest PAG-Q neighbor distances for several random and DSP structures, all for 90 mers with a polydispersity of 1.0. Figures 6(a) and 6(d) show that bound PAG/Q have a spike at the geminal distance of 1.5nm, as expected from Fig. 5(a), because the two species are held at specific distances by the chains. The distributions in Figs. 6(d) and 6(f) show that if PAG and Q are separated by more than one monomer on a DSP, their spatial distribution becomes very similar to that of a DSP with unbound PAG/Q and all of the random polymer cases [Figs. 6(a), 6(b) and 6(c)]. This indicates that a DSP with bound PAG and Q can only influence the PAG-Q spatial distributions if the PAG and Q are placed at a specific location close to one another on the polymer chain. 3.4.Impact of Polymer Film Structure on Exposure and DeprotectionThe EUV exposure and deprotection process substantially alters the spatial statistics of components within the resist film by changing their chemical identities. An individual photon absorption at a random location in the film is followed by the production of a cascade of ionization events that produce electrons and diverse reactive intermediates localized near the absorption site, within a distance on the order of the electron blur length .29 In a positive tone photoresist, these products interact with resist components, ultimately forming some quantity of photogenerated acids near the original absorption site. Upon heating, these acids affect the deprotection of the TBMA protecting group in a volume whose dimension is on the order of the acid diffusion length . The net effect is the addition of an irregular “clumpy” compositional texture to the existing structure of the film, the spatial characteristics of which depend on these two lengths, the proximity of a quencher, the number of absorbed photons, and the efficiencies of the acid generation and deprotection processes. For regions near the center of a resist pattern feature, where the exposure dose is relatively high and the initial acid distribution is dense, the deprotected volumes formed by those acids will overlap. Accordingly, those regions will be more uniformly deprotected and more likely to be fully soluble in the developer. For unexposed regions far from the nominal line edge, where the dose is near zero, any deprotected volume will be sparse, isolated, and unlikely to dissolve because the number of ionizable groups has not reached the critical threshold (assumed to be 50% in this work). In the median dose range, near the line edge, the film’s composition after deprotection transitions from insoluble to soluble. In this region, the irregular spatial composition is consequential and its non-uniform solubility will be manifested as edge roughness. Of particular interest in this study are the spatial distributions of a deprotection product MAA, formed from TBMA. Figure 7 compares functions for both acid diffusion length values of 2.0 and 5.0 nm. All random copolymers have a TBMA-TBMA correlation function consistent with random spatial distribution, with a decreased probability of adjacent TBMA groups in the defined sequence (see Fig. 4). Figure 7 shows that MAA produced from TBMA by the exposure and deprotection process leads to an MAA-MAA pair correlation function with values above 1 in all cases, consistent with the formation of clusters. This is to be expected since each acid can deprotect multiple TBMAs until it is quenched or no TBMA remains within its diffusion volume. This is found for both random and defined sequence polymers. As shown in Fig. 7(a), when the acid diffusion length is short all cases with random copolymers show essentially identical curves. For the defined sequence cases, however, is in general lower than for the random polymer cases at MAA-MAA distances below 2 nm. For the longer diffusion length [Fig. 7(b)], the pair correlations separate into two groups. The DSP and random sequence cases with free PAG and/or free Q have very similar MAA-MAA pair correlation functions above a 2 nm radius, with the DSP case being more narrowly dispersed below that distance, as expected. The cases with bound PAG and Q, on the other hand, are more similar to those in Fig. 7(a), indicating that acid diffusion length has less impact on the MAA-MAA nearest neighbor distances than when PAG and Q are unbound. The MAA distributions as a function of photoresist structure indicate that the spatial arrangement of resist components has somewhat less of an influence on the “clumpiness” of the deprotected regions if the acid diffusion length during post-exposure processing is short. If it is relatively long, on the other hand, it is possible that differences can be pronounced. This possibility is examined by considering the photoresist’s solubility chain by chain. Resist polymer solubility in aqueous base is controlled by the fraction of ionizable groups on the polymer chain.4 In all the compositions examined here, both HOST and MAA groups are ionizable. Ensembles of chains that comprise the random copolymer film slices include a broad range of different compositions with different fractions of ionizable HOST groups. Because the TBMA fractions have a similarly broad range across the chains, and exposure and deprotection events vary randomly, the ionizable MAA groups that form further alter each chain’s solubility in complex ways. By contrast, the defined sequence polymer is a single molecular structure. When assembled into a film, is each DSP chain uniform in composition after exposure at a given local photon dose and deprotection? The pair correlation functions do not provide this information because they analyze only the overall packing within the film. If individual chains are uniform, one can hypothesize that resist stochastic effects especially at the line edges would be attenuated. This hypothesis is tested by an additional set of calculations examining solubility after exposure using simple assumptions. We examine chain-by-chain statistics as follows. We analyze the composition of each ensemble of chains comprising a given resist composition both before and after exposure and deprotection. For an unexposed film, the properties of interest are the number of acid-labile TBMA groups and ionizable HOST groups in each chain. For an exposed film, in each case, we find the distribution of the numbers of MAA groups formed by deprotection of TBMA and the distribution of numbers of total ionizable groups (MAA + HOST) for each chain. We have performed this analysis for all cases and found that the results are well-represented by four particular cases: a random copolymer of polydispersity 1.0, with bound PAG and bound Q or unbound PAG and Q, and a defined sequence polymer of polydispersity 1.0, with bound PAG and bound Q or unbound PAG and Q. Figures 8(a)–8(c) show results for the random copolymer case with unbound PAG and Q (composition R-FP-FQ-30-1.0). Figure 8(a) presents histograms of the numbers of TBMA and HOST in the chains prior to exposure. Both have very broad distributions, centered on the expected values of 0.346 for TBMA and 0.577 for HOST for a polymer of overall composition HOST:TBMA:STYR=50:30:6.7. Figure 8(b) shows histograms of chain compositions after exposure and deprotection for . Although the TBMA distribution is centered at 0.346 (probability normalized to 1), the resulting MAA distribution is centered at 0.167 at the applied dose reflecting the total extent of deprotection being limited to 50%. The final distribution of ionizable groups, MAA + HOST, which controls the chains’ solubility, is approximately Gaussian and broadened compared with the initial HOST distribution. Only chains with can dissolve. Figure 8(c) shows histograms of the chain compositions after exposure and deprotection for . With the longer diffusion length, there emerges a significant probability that there will be chains with no MAA. This result can be understood by considering that in our model a larger diffusion length will lead to a longer catalytic chain, that is, more deprotection events per acid formed, because the accessible film volume per acid-forming event is larger. Thus, fewer acids and fewer exposure events are required to meet an overall deprotection of 50% in the film. Fewer exposure events will lead to a more “clumpy” spatial character as indicated by the significant number of chains that contain no MAA (19% to 33% of the chains depending on the film composition). Figures 8(d)–8(f) present the same analyses for random copolymer with bound PAG and Q (composition R-BP-BQ-30-1.0). The distributions are qualitatively the same as in Figures 8(a)–8(c), with shifts in all the distributions toward smaller median values (e.g., TBMA centered at 0.30 and HOST centered at 0.5, the values expected for a polymer of overall composition HOST:TBMA:STYR:PAG:Q = 50:30:6.7:10:3.3). Figures 8(g)–8(i) show the analysis of the defined sequence polymer case with bound PAG and Q (composition SA-BP-BQ-30-1.0). Figure 8(g) reflects the homogeneous composition of the ensemble of chains. Once subjected to exposure and deprotection, the random packing of these identical chains and the random characteristics of exposure leads to a broad distribution in chain compositions in the final film. Figure 8(h) shows the resulting distribution for a diffusion length . As found for the random copolymer cases, the MAA and MAA + HOST distributions are roughly Gaussian although with the defined sequence polymer, their widths are somewhat narrower. Figure 8(i) shows histograms of the chain compositions after exposure/deprotection where . As in the random copolymer case, there are a significant number of chains (fraction) where no MAA is formed, suggesting that acids are immediately quenched in the sequence A case (Table 2) even though PAGs and Qs are not adjacent on the chain, they are close enough to interact within the film (Fig. 7). However, there is large number of chains (fraction) that have all TBMA groups in the volume converted to MAA, creating a bimodal distribution. This means that at the assumed dose and diffusion length, when the acids formed can access a larger polymer film volume, TBMA is more likely to be fully deprotected. By starting with a uniform chain composition, rather than a very complex random composition distribution, and allowing each acid to react over a large volume, it is more likely that the deprotection process will go to completion. This bimodal form is replicated in the total fraction of ionizable groups per chain because each chain has an identical number of HOST groups. The MAA + HOST histogram is simply displaced to the right by a fixed amount equal to the fraction of HOST. These distributions can be compared with a DSP with free PAG and Q (composition SF-BP-BQ-30-1.0), Figures 8(j)–8(l). The ionizable fraction distributions in Fig. 8(b) are narrower than in Fig. 8(k), but their median values are very similar, indicating that it is the state of the PAG and Q, not the polymer structure, that most strongly influences solubility when the acid diffusion distance is 2 nm. At 5 nm, Fig. 8(l) is intermediate between the distributions found in Figs. 8(c) and 8(i). These findings raise the interesting question of how the combination of acid blur and defined sequences affect local solubility at line edges and hence line edge roughness, as well as what the resulting tradeoffs between resolution and line edge roughness might be. The present work was only performed for one dose, chosen to form 50% ionizable groups on average. These studies are currently being extended to examine a range of doses to evaluate the R-L-S tradeoff for DSPs. 4.ConclusionsIn this paper, we report a computational study of the internal structure of polymer slices representing photoresists modeled on the ESCAP composition, containing TBMA, styrene, and HOST on the polymer chains, and PAG and quencher either mixed in with the polymer chains or bound to them. The polymer chains are either random or with a defined sequence, with variations in their length and polydispersity. After assembly of the polymer slices, statistical analyses of the spatial correlations from PAG-PAG, TBMA-TBMA, and PAG-Q have been made. Comparison of the internal structures of random and defined sequence polymers of the same compositions show differences in the nearest neighbor identities (within 1–2 nm), but no differences beyond that distance that would indicate that DSPs have long-range order. These correlations are not affected by chain length or polydispersity. Imaging of these polymer slices using a simple exposure/deprotection model shows that the numbers of soluble chains formed using a 2-nm acid diffusion length have very similar distributions for both DSP and random copolymers with and without bound PAG and Q. Only one case was found to have distinctive solubility outcomes, DSP with either unbound or bound PAG and Quencher with an acid diffusion length of 5 nm. A remarkable bimodal distribution emerges in this case whose physical origin is not clearly related to spatial distributions of the resist components. How this finding depends on dose and how it affects image quality are currently under study. Code and Data AvailabilityThe data presented in this paper are publicly available in Zenodo, DOI: 10.5281/zenodo.13765732. The archived version of the codes described in this paper can be freely accessed through the Zenodo link as well as GitHub ( https://github.com/hnsbrg/PolyScope). AcknowledgmentsThis paper is based upon work by Columbia Hill Technical Consulting (WDH) and by the Center for High Precision Patterning Science (CHiPPS), an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Sciences (FAH). The authors are grateful to Dr. Matthias Mueller-Fischer for permission to use his algorithm from his PhD dissertation and to Dr. Greg Wallraff for helpful discussions. ReferencesG. M. Gallatin, F. A. Houle and J. L. Cobb,
“Statistical limitations of printing 50 and 80 nm contact holes by EUV lithography,”
J. Vac. Sci. Technol. B, 21
(6), 3172
–3176 https://doi.org/10.1116/1.1629294 JVTBD9 1071-1023
(2003).
Google Scholar
W. D. Hinsberg et al.,
“Contribution of EUV resist counting statistics to stochastic printing failures,”
J. Micro-Nanopatterning Mater. Metrol., 20
(1), 014603 https://doi.org/10.1117/1.Jmm.20.1.014603
(2021).
Google Scholar
G. M. Gallatin,
“Resist blur and line edge roughness,”
Proc. SPIE, 5754 38
–52 https://doi.org/10.1117/12.607233 PSISDG 0277-786X
(2005).
Google Scholar
P. C. Tsiartas et al.,
“The mechanism of phenolic polymer dissolution: a new perspective,”
Macromolecules, 30
(16), 4656
–4664 https://doi.org/10.1021/ma9707594 MAMOBX 0024-9297
(1997).
Google Scholar
S. Fleischmann and V. Percec,
“Synthesis of well-defined photoresist materials by SET-LRP,”
J. Polym. Sci. Pol. Chem., 48
(10), 2251
–2255 https://doi.org/10.1002/pola.24007
(2010).
Google Scholar
Y. Guo et al.,
“Terpolymerization of styrenic photoresist polymers: effect of RAFT polymerization on the compositional heterogeneity,”
Macromolecules, 48
(11), 3438
–3448 https://doi.org/10.1021/acs.macromol.5b00683 MAMOBX 0024-9297
(2015).
Google Scholar
G. G. Barclay et al.,
“Narrow polydispersity polymers for microlithography: synthesis and properties,”
Proc. SPIE, 2724 249
–260 https://doi.org/10.1117/12.241823
(1996).
Google Scholar
M. X. Wang et al.,
“Novel anionic photoacid generators (PAGs) and corresponding PAG bound polymers for sub-50 nm EUV lithography,”
J. Mater. Chem., 17
(17), 1699
–1706 https://doi.org/10.1039/b617133h JMACEP 0959-9428
(2007).
Google Scholar
D. N. Tuan, H. Yamamoto and S. Tagawa,
“Study on resist performance of polymer-bound and polymer-blended photo-acid generators,”
Jpn. J. Appl. Phys., 51
(8), 086503 https://doi.org/10.1143/JJAP.51.086503
(2012).
Google Scholar
K. Young-Gil et al.,
“A positive-working alkaline developable photoresist based on partially tert-Boc-protected calix[4]resorcinarene and a photoacid generator,”
J. Mater. Chem., 12
(1), 53
–57 https://doi.org/10.1039/b106278f JMACEP 0959-9428
(2002).
Google Scholar
J. De Neve et al.,
“Sequence-definition from controlled polymerization: the next generation of materials,”
Polym. Chem., 9
(38), 4692
–4705 https://doi.org/10.1039/C8PY01190G
(2018).
Google Scholar
Z. L. Li et al.,
“Hierarchical nanomaterials assembled from peptoids and other sequence-defined synthetic polymers,”
Chem. Rev., 121
(22), 14031
–14087 https://doi.org/10.1021/acs.chemrev.1c00024 CHREAY 0009-2665
(2021).
Google Scholar
F. Kaefer et al.,
“Controlled sequence photoresists from polypeptoids,”
J. Photopolym. Sci. Technol., 35
(1), 29
–33 https://doi.org/10.2494/photopolymer.35.29 JSTEEW 0914-9244
(2022).
Google Scholar
F. Käfer et al.,
“Polypeptoids: exploring the power of sequence control in a photoresist for extreme-ultraviolet lithography,”
Adv. Mater. Technol., 8
(23), 2301104 https://doi.org/10.1002/admt.202301104
(2023).
Google Scholar
H. Ito et al.,
“Environmentally stable chemical amplification positive resist; principle, chemistry, contamination resistance and lithographic feasibility,”
J. Photopolym. Sci. Technol., 7 433
–447 https://doi.org/10.2494/photopolymer.7.433 JSTEEW 0914-9244
(1994).
Google Scholar
G. M. Wallraff and W. D. Hinsberg,
“Lithographic imaging techniques for the formation of nanoscopic features,”
Chem. Rev., 99
(7), 1801
–1821 https://doi.org/10.1021/cr980003i CHREAY 0009-2665
(1999).
Google Scholar
T. Manouras and P. Argitis,
“High sensitivity resists for EUV lithography: a review of material design strategies and performance results,”
Nanomaterials, 10
(8), 1593 https://doi.org/10.3390/nano10081593
(2020).
Google Scholar
M. Mueller, The structure of dense polymer systems: geometry, algorithms, software, 117 Eidgenossische Technische Hochschule, Zurich
(1999). Google Scholar
G. M. Schmid et al.,
“Understanding molecular level effects during post exposure processing,”
Proc. SPIE, 4345 1037
–1047 https://doi.org/10.1117/12.436829 PSISDG 0277-786X
(2001).
Google Scholar
G. P. Patsis, V. Constantoudis and E. Gogolides,
“Effects of photoresist polymer molecular weight on line-edge roughness and its metrology probed with Monte Carlo simulations,”
Microelectron. Eng., 75
(3), 297
–308 https://doi.org/10.1016/j.mee.2004.06.005 MIENEF 0167-9317
(2004).
Google Scholar
G. P. Patsis et al.,
“Simulation of roughness in chemically amplified resists using percolation theory,”
J. Vac. Sci. Technol. B, 17
(6), 3367
–3370 https://doi.org/10.1116/1.591012 JVTBD9 1071-1023
(1999).
Google Scholar
R. C. Team, R: a language and environment for statistical computing, R Foundation for Statistical Computing, Vienna
(2023). Google Scholar
P. Team, RStudio: integrated development environment for R, Posit Software, PBS, Boston, Massachusetts
(2024). Google Scholar
A. Baddeley, E. Rubak and R. Turner, Spatial point patterns: methodology and applications, R. Chapman and Hall/CRC Press, London
(2015). Google Scholar
, “QtiPlot – data analysis and scientific visualisation,”
https://www.qtiplot.com/
(2023).
Google Scholar
P. de Schepper et al.,
“XAS photoresists electron/quantum yields study with synchrotron light,”
Proc. SPIE, 9425 942507 https://doi.org/10.1117/12.2085951 PSISDG 0277-786X
(2015).
Google Scholar
T. Kozawa and S. Tagawa,
“Radiation chemistry in chemically amplified resists,”
Jpn. J. Appl. Phys., 49
(3), 030001 https://doi.org/10.1143/JJAP.49.030001
(2010).
Google Scholar
U. Okoroanyanwu, Chemistry and Lithography, Second Edition, Vol. 2: Chemistry in Lithography, SPIE Press, Bellingham, Washington
(2023). Google Scholar
T. Watanabe et al.,
“CA resist with side chain PAG group for EUV resist,”
J. Photopolym. Sci. Technol., 19
(4), 521
–524 https://doi.org/10.2494/photopolymer.19.521 JSTEEW 0914-9244
(2006).
Google Scholar
M. X. Wang et al.,
“Novel polymeric anionic photoacid generators (PAGs) and corresponding polymers for 193 nm lithography,”
J. Mater. Chem., 16
(37), 3701
–3707 https://doi.org/10.1039/b607918k JMACEP 0959-9428
(2006).
Google Scholar
BiographyWilliam D. Hinsberg received his doctorate in chemistry at the California Institute of Technology, followed by postdoctoral work at Stanford University. He joined IBM Corporation in 1982, working in thin-film process chemistry for magnetic recording devices and the chemistry of new photoresist materials. He has numerous awards, including an IBM Corporate Award for development of stochastic kinetics simulation methods. He has published more than 150 scientific papers and is an inventor of more than 30 issued US patents. Currently, he is founder and president of Columbia Hill Technical Consulting. Frances A. Houle received her doctorate in chemistry at the California Institute of Technology. She has held appointments at the IBM Research Division and InVisage Technologies and is now a senior scientist at Lawrence Berkeley National Laboratory. Her current research focuses on chemical processes in EUV photolithography and solar photochemistry for energy conversion. She has received numerous awards including an IBM Corporate Award and the AVS John A. Thornton Memorial Award and Lecture. She has published more than 170 scientific papers and has 28 issued US patents. ORCID: 0000-0001-5571-2548 |