As complex NTD resist behavior significantly impacts AF wafer printing , there is a need to better model NTD AF printing. We present our work to further enhance NTD compact modeling accuracy for AF printing prediction by physics-based and data-based methods. The physics-based enhancements are derived from improvements in behavioral mechanisms in exposed and partially exposed NTD materials. The data-based enhancements are derived from learning methodologies developed for predicting lithography hot-spots at the limits of process control. Both types of enhancements are needed to predict fine changes in imaging and resist behavior where traditional compact models break down.
In photolithography, we need accurate models as computation engine for optical proximity correction (OPC). Traditional OPC modeling consists of a series of components for photo mask, optical exposure system, and resist materials. These models are trained using compact model forms based on wafer-level critical dimension (CD) or edge placement error (EPE) measurements. In recent years, advancements in neural networks and machine learning have had significant advancements. In this work, we evaluated advanced neural network-based resist models on a Tensor Flow machine learning platform. This work describes resist and optical response of machine learning (ML) model through process window to achieve improved model representation of lithography process. Using ML OPC vias mask as an example, we will show improved accuracy through dose and focus process conditions and verify model accuracy with physical hardware data. Also, we will compare multiple neural network-based modeling approaches, investigate the ML models’ impacts on OPC correction and verification recipes, and dataprep runtime. The machine learning based OPC with ML model and best practice will be implemented in cloud production environment.
We provide background on differences between traditional and machine learning modeling. We then discuss how these differences impact the different validation needs of traditional and machine learning OPC compact models. We then provide multiple diverse examples of how machine learning OPC compact validation modeling can be appropriately validated both for modeling-specific production requirements such as model signal/contour accuracy, predictiveness, coverage and stability; and also general OPC mask synthesis requirements such as OPC/ILT stability, convergence, etc. Finally we conclude with thoughts on how machine learning modeling methods and their required validation methods are likely to evolve for future technology nodes.
Due to the semiconductor industry’s ever increasing need for finer resolution and improved critical dimension (CD) control, negative tone development (NTD) photoresists (resists) have been adopted for several advanced applications in lithographic patterning. NTD resists enable brightfield imaging by using an organic solvent developer to penetrate and remove the unexposed regions of the resist [1]. For certain critical patterning layers, such as metal trenches and vias, NTD resists are able to provide better resist imaging quality compared to the previous positive tone development (PTD) resist process. However, there are several additional engineering difficulties which must be addressed for an NTD resist process. Specifically, NTD resists have low contrast organic solvent development and in an NTD process the material remaining on the wafer substrate is exposed resist which has been substantially transformed both chemically and mechanically. Therefore, the remaining exposed resist shows significantly more complex physical behavior than the remaining PTD resist and these behaviors require substantial improvement in an OPC (compact) model’s physical modeling accuracy in order to match wafer data and trends [2,3]. Additionally, these more complex resist behaviors place further requirements on the physical validation of OPC modeling inference. In this paper, we present results of our work to understand and improve the optimization and physical validation of physics-based NTD compact modeling flows by utilizing new methods for analysis and automation. We utilize a complete compact model flow containing physics-based resist model forms for chemically amplified resist (CAR) exposure, CAR reaction-diffusion, resist top-loss due to exposure combined with post-exposure bake (PEB), low contrast organic solvent development of resist, and mechanical deformation effects in multiple process steps. We present solid evidence that this physically-based flow has been validated for accuracy and predictability by comparing it to several experimental NTD datasets and to results of rigorous 3D lithography simulation models which were trained to fit other experimental NTD data. We additionally compared key physics-based model forms from the compact model to the more complex full time-based moving surface NTD models of the rigorous 3D simulation. We next analyzed the key physics-based compact model forms for sensitivity to input testpattern type, layout and mask dimension (e.g., linearity and MEEF), traditional dose-focus variations, as well as systematic and random noise in CD metrology. We present the results of this study and make recommendations for minimum testpattern and overall process space data to include in NTD compact model datasets. We also present flow benefits obtained from automating different validation tests including the usefulness of employing rigorous lithography simulation NTD results early in the compact modeling flow to improve overall model quality. [1] S-H. Lee, et all. Understanding dissolution behavior of 193-nm photoresists in organic solvent developers.
As feature resolution and process variations continue to shrink for new nodes of both DUV and EUV lithography, the density and number of devices on advanced semiconductor masks continue to increase rapidly. These advances cause significantly increased pressure on the accuracy and efficiency of OPC and assist feature (AF) optimization methods for each subsequent process technology. Several publications and industry presentations have discussed the use of neural networks or other machine learning (or even deep learning) to provide improvements in efficiency for OPC main feature optimization or AF placement. However, these two mask synthesis steps are not independent. OPC affects AF optimum position and size; and AF position and size both affect the final optimum OPC main feature correction. A challenging example of these interactions is the need for OPC and AF methods to be aware of potential AF wafer printing. AF printing on the wafer can lead to catastrophic device failure. If an AF is at risk of printing in photoresist, both the OPC and the size (and potentially the position) of the AF need to be modified accurately and efficiently. Recent advancements in lithography utilizing negative tone develop (NTD) photoresists (resists) with strong physical shrink effects also further increase the difficulty of accurately modeling AF printing. In this paper, we present results of our work to explore the requirements, the issues and the overall potential for developing robust, accurate and fast integrated machine learning methods to optimize OPC and AFs.
The model accuracy for Optical Proximity Correction (OPC) and Inverse Lithography Technology (ILT) applications is a critical factor for patterning success in advanced technology nodes. One difficult challenge has been the accurate and fast simulation of Negative Tone Development (NTD) photoresist processes. It has been widely observed that CD measurements, top-down SEM contour images and X-section or AFM resist side wall profile measurements cannot be adequately predicted by conventional lithography process models, e.g., [1]. Therefore, conventional OPC models were often unable to meet the demanding accuracy requirements of advanced logic or memory manufacturing. A key to achieving the demanding model accuracy requirements in NTD photoresist processes is to consider the photoresist shrinkage effects both in the Post-Exposure-Bake (PEB) step and photoresist development process step. Starting from continuum elastic mechanics, e.g., [2], we have developed a fast and accurate full 3D compact shrinkage model and validated its accuracy and usefulness vs. experimental results for several advanced processes and vs. rigorous simulation using a full physical lattice model. The compact model captures the significant photoresist shrinkage and deformation behaviors such as surface topography, resist sidewall angle (SWA) and layout pattern dependency [3], with much faster runtime capable of use in OPC and ILT mask data optimization. The speed and flexibility of the model are such that it can also be applied to help increase accuracy of simulation for some complex physical behaviors seen in other photoresist processes such as EUV and positive tone develop (PTD) photoresist.
The optical proximity correction (OPC) model and post-OPC verification that takes the developed photoresist (PR) 3D profile into account is needed in the advanced 2Xnm node. The etch process hotspots caused by poor resist profile may not be fully identified during the lithography inspection but will only be observed after the subsequent etch process. A complete mask correction that targets to final etch CD requires not only a lithography R3D profile model but also a etch process compact model. The drawback of existing etch model is to treat the etch CD bias as a function of visibility and pattern density which do not contain the information of resist profile. One important factor to affect the etch CD is the PR lateral erosion during the etch process due to non-vertical PR side wall angle (SWA) and anisotropy of etch plasma source. A simple example is in transferring patterns from PR layer to thin hard mask (HM) layer, which is frequently used in the double pattern (DPT) process. The PR lateral erosion contributes an extra HM etch CD bias which is deviated from PR CD defined by lithography process. This CD bias is found to have a nontrivial dependency on the PR profile and cannot be described by the pattern density or visibility. In this report, we study the etch CD variation to resist SWA under various etch conditions. Physical effects during etch process such as plasma ion reflection and source anisotropy, which modify the local etch rate, are taken into considerations in simulation. The virtual data are generated by Synopsys TCAD tool Sentaurus Topography 3D using Monte Carlo engine. A simple geometry compact model is applied first to explain the behavior of virtual data, however, it works to some extent but lacks accuracy when plasma ion reflection comes into play. A modified version is proposed, for the first time, by including the effects of plasma ion reflection and source anisotropy. The new compact model fits the nonlinear etch CD bias very well for a wide range of resist SWAs from 65 to 90 degrees, which covers the resist profile diversities in most real situations. This result offers a potential application for both resist profile aware and etch process aware mask correction model in the mask synthesis flow.
The Optical Proximity Correction (OPC) model, a key to process yield in the mask synthesis flow, is getting more and more complicated and challenging at advanced technology nodes (1X nm). To achieve accurate critical dimension (CD) prediction and model robustness on varied designed patterns, a rigorously tuned compact model (RTCM) [1] that takes the photoresist chemical effects into considerations is strongly desired. A lithography process consists of three main stages: Exposure, Post-Exposure Bake (PEB), and Photoresist Development. Each stage is characterized by its fundamental physics or chemistry that governs the process of illumination induced photo-acid generation, thermally activated chemical reaction-diffusion, and developer dependent photoresist dissolution, respectively. The final resist profile is determined by the process details of all these stages directly or indirectly. For an ideal resist that the development contrast approaches infinity, resist development is aptly represented by a threshold model applied to the PEB latent image (acid or inhibitor concentration). So the quality of OPC modeling is largely determined by the fidelity of PEB latent image [2,3]. However, for some types of resist and developer used in Negative Tone Development (NTD), the development contrast shows a long tail without a sharp transition. For such low-contrast resist, the developed resist profile is no longer described well by the equilevel surface of PEB latent image. Going beyond the threshold approximation, we start from the fundamental equations of resist development physics and analyze the time evolution of development front that determines the resist profile. In this paper, a new compact model is derived to catch the main physics in resist Development, which is also simple and computationally efficient to suit for OPC applications. Comparison with S-LITHO rigorous solutions and real-wafer experiments with 1D and 2D test patterns have showed that the new compact model, with fewer free parameters, provides better CD prediction than the existing empirical lumped parameter models for low-contrast resists. The new physical compact model offers a more accurate and extendable solution for OPC modeling at the 10nm node and onward.
3D Resist profile aware OPC has becoming increasingly important to address hot spots generated at etch processes
due to the mass occurrence of non-ideal resist profile in 28nm technology node and beyond. It is therefore critical to
build compact models capable of 3D simulation for OPC applications. A straightforward and simple approach is to
build individual 2D models at different image depths either based on actual wafer measurement data or virtual
simulation data from rigorous lithography simulators. Individual models at interested heights can be used by
downstream OPC/LRC tools to account for 3D resist profile effects. However, the relevant image depths need be
predetermined due to the discontinuous nature of the methodology itself. Furthermore, the physical commonality
among the individual 2D models may deviate from each other as well during the separate calibration processes. To
overcome the drawbacks, efforts are made in this paper to compute the whole bulk image using Hopkins equation in
one shot. The bulk image is then used to build 3D resist models. This approach also opens the feasibility of
including resist interface effects (for example, top or bottom out-diffusion), which are important to resist profile
formation, into a compact 3D resist model. The interface effects calculations are merged into the bulk image
Hopkins equation. Simulation experiments are conducted to demonstrate that resist profile heavily rely on interface
conditions. Our experimental results show that those interface effects can be accurately simulated with reference to
rigorous simulation results. In modeling reality, such a 3D resist model can be calibrated with data from discrete
image planes but can be used at arbitrary interpolated planes. One obvious advantage of this 3D resist model
approach is that the 3D model is more physically represented by a common set of resist parameters (in contrast to
the individual model approach) for 3D resist profile simulation. A full model calibration test is conducted on a
virtual lithography process. It is demonstrated that 3D resist profile of the process can be precisely captured by this
method. It is shown that the resist model can be carried to a different lithography process with same resist setup but
a different illumination source without model any accuracy degradation. In an additional test, the model is used to
demonstrate the capability of resist 3D profile correction by ILT.
A single compact resist model capable of predicting 3D resist profile is strongly demanded for the advanced technology
nodes to avoid the potential hotspots due to imperfect resist pattern shape and its lack of resistance in the subsequent
etch process. In this work, we propose a resist 3D (R3D) compact model that takes acidz-diffusion effect into account.
The chemical reaction between acid and base along z-direction is treated as second order effect that is absorbed into the
anisotropic diffusion length as a fitting parameter. Meanwhile, the resist model in the x-y wafer plane is still kept in
general by applying the compact solution of 2D reaction-diffusion equation. In order to have the 2D contour
predictability at arbitrary resist height, calibration from entire 3D data (CDs at several heights) areconducted
simultaneously witha single cost function so that the R3D compact model is described by a common set of resist free
parameters and threshold for all resist heights. With the low energy approximation, the acid z-diffusion effect is
equivalent to a z-diffused TCC that takes the form of linear combination of pure optical TCCs sampled at discrete
image-depth which can be pre-calculated. With this benefit, the R3D compact model offers a more physical approach but
adds no runtime concern on the OPC and verification applications. The predicted resist cross-section profiles from our
test patterns are compared those computed with rigorous lithography simulator SLITHO and show very good matching
results between them. The demonstration of the AF printability check from the predicted cross-section profile at AF
indicates the success of our R3D compact model.
Access to the requested content is limited to institutions that have purchased or subscribe to SPIE eBooks.
You are receiving this notice because your organization may not have SPIE eBooks access.*
*Shibboleth/Open Athens users─please
sign in
to access your institution's subscriptions.
To obtain this item, you may purchase the complete book in print or electronic format on
SPIE.org.
INSTITUTIONAL Select your institution to access the SPIE Digital Library.
PERSONAL Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.