Amid the different approaches to rigorously model the scattering of electromagnetic waves in sub-wavelength photolithography process, waveguide methods have been proven to be accurate and efficient. While totally different from the time domain methods such as FEM and TDFD, waveguide methods simultaneously compute diffraction of the incident plane waves with different incidence angles by a periodical dielectric structure. EM modes inside dielectric layers are solved for by decoupling the eigen-system and the electromagnetic field boundary conditions between each two adjacent layers are applied to stitch the modes and obtain the full scattering matrix. It is clear that the number of orders of the plane waves under consideration directly affects the accuracy of a waveguide simulator. However, in a 3-D simulation, simulation time and memory usage increase drastically with the increase of the number of orders. These limitations prevent waveguide methods from being applied to large layout patterns that require higher simulation orders. Since many cases under study in lithography process optimization and layout printability analysis are actually axial-symmetric about both x and y axes, it is possible to use this symmetry to simplify the calculation. Based on our definitions of the four fundamental groups of symmetric and anti-symmetric functions and the operations between these symmetric groups, we have made some fundamental discoveries of the waveguide propagating behaviors inside symmetric dielectric structures. In this paper, we will describe our new rigorous theory to decompose the propagating waves in a symmetric layer into four symmetric and anti-symmetric transmitting plane wave groups. Each group can be directly stitched to the next layer without losing symmetry or anti-symmetry. Furthermore, any form of incident waves can be reorganized into the above four plane wave groups, so can the reflective and refraction fields. We have shown that by using this field decomposition technique, we can compute scattering of the four smaller diffraction mode groups separately and the computational complexity will be greatly reduced as compared with the full mode simulation. We have developed a new version of the 3-D simulator, METRO, which incorporated this new method and achieved very high computational efficiency. For example, a 33 by 33-order, double precision simulation will only require 300M bytes memory and take less than 15 minutes without any loss of accuracy on a 1G processor compared with 1.7G memory usage and 10 hours simulating time on the same processor without symmetric simplification.
Printability challenges are already one of the key limits of new technology generations and the situation will become even worse with the delay in the introduction of 157nm lithography. To optimize the set-up of RET and verify their efficiency within a photo process window, it has become necessary to run detailed lithography simulations. While some 1-D test patterns can be simulated using the 2-D software tools, patterns such as line-ends, corners and metal islands should be evaluated by 3-dimensional simulators. Such a simulator can be calibrated using a small number of in-line measurements and then employed to analyze the “difficult to print patterns”. However, until recently the computational efficiency of such a 3-D tools has not been sufficient for practical applications. Since a lot of features under study are axial-symmetric, symmetry of the 3-D structures can be used to reduce computational complexity. To address this issue, we have developed a method to simulate three-dimensional axial-symmetrical structures. This method is based on the observation that scattering in and out of axial symmetrical structure can be classified into four kinds of scattering waves and each kind of scattering wave can be computed separately. Based on this observation, we have greatly reduced the run time and memory usage of the rigorous waveguide method without any loss of accuracy. According to our experiments, memory usage was reduced to one-forth and run time was shortened to 3~5% of the original simulation. (20-30 times speed up) In this paper, we will discuss applications of the new symmetrical simulator METROPOLE-3D to lithography simulation and layout analysis. We will discuss the fast calibration of METROPOLE-3D simulator. The simulated CD data was compared with experimental FEM wafer data to calibrate the parameters used in PEB and photoresist development. We are now able to analyze a lot of axial symmetrical dense and isolated 2-D structures with unprecedented speed. In this paper, wee will present a full lithography analysis example from a real product manufactured in the state-of-the-art lithography process.
Conventional lithography techniques have been losing their ability to easily support continuous shrinking of feature sizes, especially when the pattern half-pitch is <60 nm. EUV lithography is one of the leading contenders to replace these conventional techniques. Because the EUV mask structure is many times thicker than the illumination wavelength, scalar Fraunhofer diffraction calculations cannot describe the scattering of light from the EUV masks with enough accuracy. In this paper, we present a rigorous 3-D electromagnetic field simulator for EUV masks. The simulator is based on a 3-D waveguide method developed to calculate the characteristics of the light scattered from nonplanar EUV masks. A typical EUV mask contains as many as 80 reflective layers in addition to the absorbing layers, and we have developed a fast method to calculate the scattering matrix of the reflective layers. Also, based on the existing numerical techniques, we can describe the light scattering in the absorbing layers with complex index of refraction. The simulator has been recently modified to handle oblique illumination conditions, and this is the focus of our paper. Aerial images are calculated in the image plane of a typical EUV stepper, and a threshold resist model is used to predict the printed pattern size. We will first compare our work with published results on dense and isolated lines. Then, we will describe the results of our calculations for two-dimensional patterns (e.g., contacts and islands) under oblique illumination. The typical simulation time is less than 10 hours on a desktop workstation for two-dimensional patterns.
EUV lithography at 13nm wavelengths will require the use of reflective multi-layer mask substrates with a patterned absorbing top layer. As the height of the absorber is many times the wavelength of the exposure radiation, EUV masks are true topographical entities. Therefore, scalar Frauenhofer diffraction simulations are not completely valid and accurate. Rigorous solution of Maxwell's equations in the topography is necessary to accurately predict many of the lithographic effects of these masks. Additionally, the interactions of the non-vertical radiation reflectance from the multilayer stack with the absorber causes further non-intuitive lithographic results which need to be understood and optimized. To further the rigorous understanding of EUV mask effects, we have extended existing 2D and 3D rigorous simulators to model EUV lithography diffraction. These simulators, METRO and METROPOLE-3D respectively, use the waveguide method for fast and accurate computation of topographical mask structures. Modifications of these simulators enable an order of magnitude speedup in the calculation of results containing large numbers of EUV multilayers, improve mask dimension modeling accuracy and enable accurate modeling of the non-normal EUV source illumination through a range of pitch values. Rigorous illumination and multilayer defect analysis results will also be shown. Additional comparisons to recently published rigorous EUV simulations results and impressive runtime results on standard desktop workstations are presented. The results show that for EUV lithography the best focus position and the symmetry of the aerial image intensity profile with respect to best focus is affected by the pattern pitch.
The production of 70nm devices is projected for the year 2008. With this projection, optical lithography will become more challenging since as the device size goes down, the potential for introducing killer defects also increases dramatically. Wafer inspection will play a key role in controlling the defect mechanisms and keeping an acceptable yield for next generation VLSI manufacturing. Metrology tools for the new generation lithography will have the following features. First, the projection and collection lenses will have higher numerical apertures (NA) to obtain high-resolution images. Typically, the NA will be as high as 0.9. Second, the wavelength used for wafer inspection will be much smaller so that common wafer materials will become highly absorptive. Last, with the increased number of process steps, wafer inspection will need to provide information for more critical processes. Based on these features, accurate modeling of the next generation wafer inspection schemes is needed to aid in the characterization and optimization of the inspection tools. The simulation tool must be able to simulate inspection systems with high NA lens, DUV wavelengths, and highly absorptive wafer materials accurately and quickly. At Carnegie Mellon University, a simulator called METRO-3D was developed into a defect inspection simulator for DUV lithography processes. This simulator is able to successfully model various types of defect mechanisms. However, it has experienced occasional numerical instabilities, with discontinuous dielectric structures composed of highly absorptive materials. In order to provide a tool to simulate the wafer inspection scheme of SUB-70NM NODE LITHOGRAPHY, we have incorporated a new algorithm to model the wafer inspection system more accurately and robustly. Numerical experiments shows that the algorithm is capable of simulating topographies with discontinuous dielectric functions, yielding stable results even when the material is highly absorptive. To verify the accuracy of the simulator, several simulations were compared with both analytical models and results form other existing simulators. The results show good matches between METRO-3D and these well-established results. Finally, we performed simulations on industrial data and the results exemplified the ability of METRO-3D to model complex 3-D structures. In this paper we will present is efficient and stable EM solver and the results of the simulator applied to various sub-70nm node wafer inspection schemes.
In recent years, optical proximity correction (OPC) has become major technology to compensate the undesirable pattern distortions caused by systematic photolithography and etch variations. Model based OPC usually utilizes a scalar method to model diffraction limitations and coherence of optical system and empirical models to include etch variations. By calibrating the simulator to the experiential data from the focus exposure matrix (FEM) wafers, the lithography and etch errors can be predicted to some extent for 1D structures such as long lines and gaps. However, it has been difficult to predict the printability of such 2D patterns as line-ends, islands and small features with the simplified simulator. Furthermore, since the numerical aperture (NA) tends to increase in the state-of-the-art lithography and the wafer and mask topographies are becoming more complex, the optical simulator can no longer calculate aerial image correctly.
This paper presents the result of an experimental study to investigate the potential advantages of using passive vibration confinement over conventional active vibration control methods, as well as to investigate the benefits of using the two methods simultaneously. The general approach is to compare the result of actively controlling vibrations in a beam which sees various degrees of modal confinement. Vibration confinement is carried out passively, and the comparison is based on control effort required as well as vibration control performance achieved. To date, there has been a significant amount of work in the area of vibration confinement, or mode localization, but the focus has been primarily either 1) that it is an interesting phenomenon which exists in structures or 2) that it can be produced in structures through active, passive, or hybrid means to achieve some end such as vibration control. This paper presents an experimental follow-up to an earlier numerical study which directly compared confinement techniques to conventional active vibration control methods, and showed how confinement can be used to enhance conventional vibration control. Although not as dramatic, the result presented in this paper clearly support that study and show that passive vibration confinement can enhance active control through both performance and energy consumption.
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.