OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 15, Iss. 19 — Sep. 17, 2007
  • pp: 12539–12547
« Show journal navigation

Etching of photosensitive chalcogenide glasses:experiments and simulations

Raphi Dror, B. Sfez, Sh. Y. Goldin, and A. Cashingad  »View Author Affiliations


Optics Express, Vol. 15, Issue 19, pp. 12539-12547 (2007)
http://dx.doi.org/10.1364/OE.15.012539


View Full Text Article

Acrobat PDF (272 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We have developed a three-dimensional simulation algorithm based on fast marching method that mimics the etching behavior of chalcogenide photoresists, especially for maskless interference lithography. This lithography exposure is characterized by continuous variation of the exposure intensity inside the photoresist, without step like variation. Furthermore, the chalcogenide photoresist has a “gray-scale” behavior, without definite threshold. The resulting etching process is very sensitive to exposure dose and etching time. The optimal relations between these parameters are determined both theoretically and experimentally. A very good agreement between calculation and experimental results is shown, opening the door to complex nanostructures engineering.

© 2007 Optical Society of America

1. Introduction

The exposure of thin amorphous chalcogenide layers to light of suitable wavelength and intensity can result in significant structural changes due to their very unique physicochemical properties [1

1. K. Shimakawa, A. Kolobov, and S. R. Elliott, “Photoinduced effects and metastability in amorphous semiconductors and insulators,” Adv. Phys. 44, 475 (1995). [CrossRef]

4

4. G. Rosenblum, B. G. Sfez, Z. Kotler, V. Lyubin, and M. Klebanov, “Nonlinear optical effects in chalcogenide photoresists,” Appl. Phys. Lett. 75, 3249 (1999). [CrossRef]

]. One of the most interesting properties of amorphous chalcogenide thin layers is its chemical resistance to alkaline solvents that can be modified by light exposure [5

5. V. M. Lyubin, A. M. Sedikh, N. N. Smirnova, and V. P. Shilo, Microelectronica18, 523 (1989).

, 6

6. A. Arsh, M. Klebanov, V. Lyubin, L. Shapiro, A. Feigel, M. Veinger, and B. Sfez, “Glassy mAs2S3·nAs2Se3 photoresist films for interference laser lithography,” Opt. Mater. 26, 301–304 (2004). [CrossRef]

]. This phenomenon, together with near-atomic resolution capability and high transmission in the IR region enables them to be applied as a medium for the fabrication of nano-structured optical elements such as gratings [7

7. M. Vlcek, P. J. S. Ewen, and T. Wagner, “High efficiency diffraction gratings in As-S layers,” J. Non-Cryst. Solids 227230, 743 (1998). [CrossRef]

, 8

8. A. V. Stronski, M. Vlcek, A. Sklenar, P. E. Shepeljavi, S. A. Kostyukevich, and T. Wagner, “Application of As40S60-xSex layers for high-efficiency grating production,” J. Non-Cryst. Solids 266269, 973 (2000). [CrossRef]

] and filters. Additionally, chalcogenide glasses possess a very high refractive index, ranging from 2.3 to 3.5 which make them attractive materials for fabrication of photonic crystals [9

9. S. Wong, M. Deubel, F. Pérez-Willard, S. John, G. A. Ozin, M. Wegener, and G. von Freymann, “Direct Laser Writing of Three- Dimensional Photonic Crystals with a Complete Photonic Bandgap in Chalcogenide Glasses,” Adv. Mater. 18, 265–269 (2006). [CrossRef]

]. Recently our group has demonstrated double exposure two-dimensional photonic crystal and layer-by-layer three-dimensional (3D) photonic crystal (wood-pile) structures using maskless interference lithography [10

10. A. Feigel, M. Veinger, B. Sfez, A. Arsh, M. Klebanov, and V. Lyubin, “Two dimensional photonic band gap pattering in thin chalcogenide glassy films,” Thin Solid Films 488, 185–188 (2005). [CrossRef]

12

12. A. Feigel, M. Veinger, B. Sfez, A. Arsh, M. Klebanov, and V. Lyubin, “Three-dimensional simple cubic woodpile photonic crystals made from chalcogenide glasses,” Appl. Phys. Lett. 83, 4480 (2003). [CrossRef]

]. The interference lithography method is characterized by long range order periodic pattern, very high resolution, and infinite depth of focus. Therefore it is very attractive for fabrication of dense gratings and 2D and 3D photonic crystal.

In interference lithography, the exposure intensity inside the photoresist varies continuously, without sharp steps as in mask lithography. On the other hand, the etching rate of most photosensitive chalcogenide does not present the definite threshold characteristics of commercial organic photoresists, but is rather near to linear for lower exposure energies and saturates at high energies. Therefore it presents a gray-scale behavior. The combination of gray-scale intensity patterning and gray-scale etching behavior of the glass makes interference lithography patterning in chalcogenide material difficult to predict and therefore to engineer.

Currently there is no report of quantitative study of the etching process in chalcogenide glasses. Although simulation tools developed for virtual linear photosensitivity photoresist [17

17. J. A. Sethian and D. Adalsteinsson, “Overview of level set methods for etching, deposition, and lithography development,” IEEE Transactions on Semiconductor Devices 10, 167–184 (1997). [CrossRef]

] or commercial photoresist like SU-8 [13

13. R. C. Rumpf and E. G. Johnson, “Fully three-dimensional modeling of the fabrication and behavior of photonic crystals formed by holographic lithography,” J. Opt. Soc. Am. A 21, 1703–1713 (2004). [CrossRef]

] can be adapted, the specific etching of chalcogenide glasses at the sub-micron level has to be addressed specifically. The objective of this study is to present a comprehensive approach to the etching of chalcogenide glasses for sub-micron patterning.

In the first part of this paper we present the development of a dedicated algorithm based on the fast marching method for the simulation of the etching process. In the second part we describe the experimental conditions for the quantitative study. In the third part we present a detailed comparison between simulation and experimental results. We then propose a simple predictive model that is tested against experiments.

2. Developing of the etching model

2.1 Light interference pattern

The exposure model determines the distribution of absorbed energy within the resist. Four coplanar beams interfere: two primary beams coming from the laser and two reflections from the substrate. The primary sinusoidal interference pattern results from the contribution of two TE polarized coherent laser beams. This image is disturbed by absorption within the chalcogenide layer and by air and substrate reflections. The indexes of refraction and the coefficient of absorption involved are measured and, in the simulation, the complex Fresnel coefficient is calculated to determine the amplitude and phase of the reflected beams.

2.2 Formulation of Etching Model

Simulation of the etching process is a typical problem of front propagation. In this work we use the fast marching method developed by Sethian and Osher [14

14. J. A. Sethian, “A fast marching level set method for monotonically advancing fronts,” Proc. Nat. Acad. Sci. 93, 1591–1595 (1996). [CrossRef] [PubMed]

19

19. D. Adalsteinsson and J. A. Sethian, “A level set approach to a unified model for etching, deposition, and lithography II: three-dimensional simulations [integrated circuits],” J. Comp. Phys. 122, 348–366 (1995). [CrossRef]

]. Recently, implementation of the method to lithography in SU-8 was presented [20

20. R. C. Rumpf, P. Srinivasan, and E. G. Johnson, “Modeling the fabrication of nano-optical structures,” Proc. SPIE 6110, 611004 (2006). [CrossRef]

21

21. R. C. Rumpf and E. G. Johnson, “Comprehensive modeling of near-field nano-patterning,” Opt. Express 13, 7198 (2005). [CrossRef] [PubMed]

]. Instead of trying to determine the position of the solvent front at every time moment, we look for the time the front arrives at any given point. The problem then reduces to the solution of a partial differential equation as explained in the following.

In the case of the “wet etching” relevant here, the front propagates in its normal direction with a local etching time t0·γ(x,y,z). γ is a parameter without units which determines the relative etching time (the etching time of the illuminated area divided by the etching time of the non-illuminated area) and t0 is the etching time of a unit length non-illuminated sample. Applying γ on the illumination pattern I(x,y,z) provides the relative local etching time γ(x,y,z). The arrival time (T) is then governed by the Eikonal equation:

T(x,y,z)=t0·γ(x,y,z)=t0·γ(I(x,y,z))
(1)

The value of t0 is measured experimentally. The numerical scheme appropriate for the solution of this equation must be able to treat the possible cases when a sharp corner develops in the front, the front separates into a few pieces or different parts of the front merge. We use both the first-order and the second-order schemes. The first-order one is defined by the following finite-difference equation:

[[max(DijkxT,Dijk+xT,0)]2+[max(DijkyT,Dijk+yT,0)]2+[max(DijkzT,Dijk+zT,0)]2]12=t0γijk
(2)

where Dijk+xT=Ti+1,jkTijkxi+1,jkxijk and DijkxT=TijkTi1,jkxijkxi1,jk are the forward and the backward partial x-derivatives respectively; D±yijk T and D±zijk T are similarly defined for j and k indices respectively. The key point here is the correct choice of the right and the left derivatives. The second-order scheme is built on the same principles.

Figure 1 presents a movie showing the time-development simulation of the etching process for a sinusoidal illumination of a chalcogenide layer. The simulation grid resolution along the substrate normal direction is 9.375 nm. The simulation grid resolution along the grating periodicity direction is 21.25 nm. The all volume size is 1.7×1.7×1.7 µm3

The reflected beams are responsible for the modulations on the grating walls. The modulation pitch size is 100 nm (half wavelength inside the photoresist).

Fig. 1. (1,390KB) The movie simulates the etching of a dielectric grating written using interference of two laser beams.. [Media 1]

3. Experimental setup and photosensitivity measurement

3.1 Samples preparation and photosensitivity measurement

For the fabrication of thin chalcogenide films (750 nm thickness) we thermally evaporated glassy 5(As2S3):1(As2Se3) in a high-vacuum evaporation chamber at a pressure of 5×10-6 mbar, and with a rate, measured with a quartz microbalance, of approximately 2–3 nms-1. The absorption coefficient and index of refraction were measured using an Ellipsometer (V.A.S.E Series of J. A. Wollam Co.). At the writing wavelength 532 nm the index of refraction of 5(As2S3):1(As2Se3) composition is 2.79 and the absorption is approximately 10,000cm-1.

Each set contained 17 samples that were thermally evaporated together, for maximum uniformity and for reproducibility analysis. In principle, chalcogenide etching properties depend not only on the total illumination energy, but also on the illumination power [4

4. G. Rosenblum, B. G. Sfez, Z. Kotler, V. Lyubin, and M. Klebanov, “Nonlinear optical effects in chalcogenide photoresists,” Appl. Phys. Lett. 75, 3249 (1999). [CrossRef]

]. In order to avoid the use of additional parameter in this study, we kept the illumination time fixed at 200 sec, and we determined γ for different values of the illumination intensity. To measure γ, seven samples thermally evaporated together were used. Six of them were illuminated for 200 sec. The illumination intensity was varied and the total dose was in the 3–30 J/cm2 range. The etching time was measured for all the samples. Then, the ratio between the etching time of illuminated and non-illuminated samples was calculated. In Fig. 2 an experimental curve of γ values versus exposure dose is shown. This curve is slightly non-linear at low intensity and saturates at high intensity. In the simulation we use an analytic function that fits the experimental results. The best fit is given by

γ=1+a·In1+b·In;
(3)

With a=0.416, b=0.00716, n=1.92. The exposure dose units are J/cm2.

Fig. 2. Experimental results and fitted curve of relative etching time versus dose illumination for 5As2S3:1As2Se3 chalcogenide photoresist. The illumination intensity was varying up to 175mW/cm2 and the illumination time was kept constant to 200 sec.

3.2 Writing setup

The set-up used to write the gratings is based on the interference of two CW laser beams (frequency doubled Nd:YAG, COHERENT VERDI-5). The incoming beam is split by a polarizing beam splitter into two equal intensity beams that form the arms of the interferometer with a sample mounted at the point where the two beams interfere. A λ/2 plate is added to one of the arms to set TE polarization in both beams. The angle at which the two beams overlap determines the spacing of the fringe pattern and sets the grating pitch Λ=λ(2·sin(θ 2)), where λ=0.532 µm is the laser wavelength and θ is the angle between the arms. Each beam is spatially filtered and then collimated to a diameter of 20 mm (FWHM). The sample size is 8×8 mm2, smaller than the beam size, exposed only to the central flat part of the beam. The laser output power is limited to 900 mWatt (450 mWatt per beam) to prevent damage to the irises of the spatial filters.

The 0.532 µm writing wavelength of the laser is in the absorption tail of the 5(As2S3):1(As2Se3) chalcogenide material in order to enable effective photo-structural process as well as deep penetration of the laser light into the photoresist film. For example, 1 µm thick layer absorbs 60% of the light.

4. Experimental results and comparison with the simulation

The aim of this study is to experimentally validate the simulation algorithm, so that it can be used as a predictive tool for the fabrication of complex structures. Therefore, we have chosen two criteria for the comparison. The first criterion for comparison is the geometry of the structures. We compared the filling factor (the fraction of resist that remains after developing) and the shape (parallel walls or trapezium shape). The second criterion is the period of time during which the etched structure is stable inside the etching agent. For appropriate selection of the illumination intensity and material composition, the model predicts a stable period of time (dwell time) during which the structure remains almost unchanged during etching. After that, the walls get thinned until the structure breaks down.

For the comparative study, several sets of samples are illuminated and etched for different illumination energy doses and etching times. Then SEM cross sections of the etched samples are performed and analyzed. Finally we simulate the etching process for the same parameters using the function shown on Fig. 6. The following table summarizes the samples used in this work:

Group I:. Under exposure; non-stable etching process. The samples were used in the calibration of the simulation software. Group II: Stable etching process. The samples were used in the calibration of the simulation software. Group III: Stable etching process. These samples were added after calibration in order to validate the prediction of the stability zone.

table-icon
View This Table

A sensitive factor for the comparison of the simulation with the experiment is the unavoidable set-up mechanical vibrations. These vibrations smooth the actual light intensity pattern and create an effective background that reduces the maximum possible contrast. We model these mechanical vibrations by convoluting the laser intensity profile in the sample with a Gaussian kernel of standard deviation of 30 nm.

4.1 Shape and etching time comparison

In a first step, we focused on fully etched samples (Group II). Experimental and simulation results for one such sample are displayed in Fig. 3. We measured the wall thickness in the upper and lower surfaces, and compared the results with simulations. For a set of samples with total laser energy ranging from 12 J/cm2 to 22 J/cm2 the mismatch in the upper surface wall thickness was less then 10%. Comparing the lower surface walls thickness yields the same agreement.

Fig. 3. SEM image and simulated picture of a 5As2S3:1As2Se3 chalcogenide dielectrics grating cross section. The exposure dose is 17.5 J/cm2 and the etching time is 105 sec.

To study the time evolution of the etching process we used non-completely etched samples (Group I). For each experimental condition, we run a simulation and recorded the time at which the simulated structure best fit the experimental SEM cross-section (Fig. 4). This time defined a simulated etching time that could be directly compared with the experimental etching time. The variation between these values was less then 10%.

Fig. 4. SEM image and simulation of non-complete etching of 5As2S3:1As2Se3 chalcogenide dielectric grating crosses section. The exposure dose is 90 J/cm2 and the etching time is 70 sec.

5. Discussion

We then plotted tE and tB as a function of the dose (Fig. 6). These two graphs define a region of parameters, shaded in yellow on Fig. 6, for which the structure is stable. In this region, the dwell time tB - tE is positive and characterizes the structure stability: larger dwell time meaning better tolerance to parameter fluctuations. It can be seen that at about 30 J/cm2, the dwell time is maximal. In order to be less sensitive to parameters fluctuations, the etching time should be chosen in the vicinity of (tE+tB)/2.

Group I experimental samples (uncompleted etching) are plotted in this graph using square red dots. These dots are outside the yellow shaded region. Group II experimental samples (complete etching) are plotted using green dots and are located within the shaded region. Since group I and group II samples have been used for the simulation calibration of the mechanical vibrations, the good agreement between the simulation and the experiment is expected.

Fig. 6. Illumination doses and etching times used for preparing gratings in the 5As2S3:1As2Se3 chalcogenide composition. The area between the two lines is a region of stable etching (nonzero dwell time).

6. Conclusion

There is an increased interest in chalcogenide materials as high index, high resolution photoresists for the fabrication of complex optical nanostructures such as two or three-dimensional photonic crystals. Naïve approaches based on intensity threshold for the prediction of the final shape of the nanostructure miss the dramatic effects many process parameters have on the final structure geometry and etching stability. We have shown that there is a delicate interplay between these parameters and that proper structure design should take these effects into account.

In order to do so, we have developed a simulation tool based on a fast marching method and compared the simulation with experimental results. We found an excellent agreement between both, provided that we introduced empirically mechanical vibrations of about 30 nm amplitude. In addition, we defined a predictive chart based on the simulation that allows fast determination of the optimal illumination and etching parameters for obtaining a stable process.

This chart has been developed for the case of simple sinusoidal illumination using continuous wave laser. However, the presented method is general and can be adapted to different chalcogenide composition and different types of physical mechanism for the light absorption, such as single or multi photon absorption, continuous or pulsed illumination. In addition, it can be easily adapted to more complicated interference structures such as two or three-dimensional photonic crystals.

Acknowledgment

The authors express their gratitude to Prof. V. Lyubin and Dr. M. Klebanov from Physics Department, Ben-Gurion University, Beer-Sheva, Israel for providing the samples and for fruitful discussions.

References and links

1.

K. Shimakawa, A. Kolobov, and S. R. Elliott, “Photoinduced effects and metastability in amorphous semiconductors and insulators,” Adv. Phys. 44, 475 (1995). [CrossRef]

2.

K. TanakaA.V. Kolobov, ed.,Photo-Induced Metastability in Amorphous Semiconductors (Wiley, Weinheim, 2003) pp 69. [CrossRef]

3.

A. Ozols and K. Shvarts, “Photosensitivity of amorphous semiconductor As-S and As-Se films under CW, nanosecond and picosecond laser irradiation,” Cryst. Latt. Def. and Amorph. Mat. 17, 235–239 (1987).

4.

G. Rosenblum, B. G. Sfez, Z. Kotler, V. Lyubin, and M. Klebanov, “Nonlinear optical effects in chalcogenide photoresists,” Appl. Phys. Lett. 75, 3249 (1999). [CrossRef]

5.

V. M. Lyubin, A. M. Sedikh, N. N. Smirnova, and V. P. Shilo, Microelectronica18, 523 (1989).

6.

A. Arsh, M. Klebanov, V. Lyubin, L. Shapiro, A. Feigel, M. Veinger, and B. Sfez, “Glassy mAs2S3·nAs2Se3 photoresist films for interference laser lithography,” Opt. Mater. 26, 301–304 (2004). [CrossRef]

7.

M. Vlcek, P. J. S. Ewen, and T. Wagner, “High efficiency diffraction gratings in As-S layers,” J. Non-Cryst. Solids 227230, 743 (1998). [CrossRef]

8.

A. V. Stronski, M. Vlcek, A. Sklenar, P. E. Shepeljavi, S. A. Kostyukevich, and T. Wagner, “Application of As40S60-xSex layers for high-efficiency grating production,” J. Non-Cryst. Solids 266269, 973 (2000). [CrossRef]

9.

S. Wong, M. Deubel, F. Pérez-Willard, S. John, G. A. Ozin, M. Wegener, and G. von Freymann, “Direct Laser Writing of Three- Dimensional Photonic Crystals with a Complete Photonic Bandgap in Chalcogenide Glasses,” Adv. Mater. 18, 265–269 (2006). [CrossRef]

10.

A. Feigel, M. Veinger, B. Sfez, A. Arsh, M. Klebanov, and V. Lyubin, “Two dimensional photonic band gap pattering in thin chalcogenide glassy films,” Thin Solid Films 488, 185–188 (2005). [CrossRef]

11.

A. Feigel, Z. Kotler, B. Sfez, A. Arsh, M. Klebanov, and V. Lyubin, “Chalcogenide glass-based three-dimensional photonic crystals,” Appl. Phys. Lett. 77, 3221 (2000). [CrossRef]

12.

A. Feigel, M. Veinger, B. Sfez, A. Arsh, M. Klebanov, and V. Lyubin, “Three-dimensional simple cubic woodpile photonic crystals made from chalcogenide glasses,” Appl. Phys. Lett. 83, 4480 (2003). [CrossRef]

13.

R. C. Rumpf and E. G. Johnson, “Fully three-dimensional modeling of the fabrication and behavior of photonic crystals formed by holographic lithography,” J. Opt. Soc. Am. A 21, 1703–1713 (2004). [CrossRef]

14.

J. A. Sethian, “A fast marching level set method for monotonically advancing fronts,” Proc. Nat. Acad. Sci. 93, 1591–1595 (1996). [CrossRef] [PubMed]

15.

J. A. Sethian, Level Set Methods and Fast Marching Methods, (Cambridge Univ. Press, 2nd ed., 1999).

16.

S. Osher and R. Fedkiw, Level Set Methods and Dynamic Implicit Surfaces (Springer, 2003).

17.

J. A. Sethian and D. Adalsteinsson, “Overview of level set methods for etching, deposition, and lithography development,” IEEE Transactions on Semiconductor Devices 10, 167–184 (1997). [CrossRef]

18.

D. Adalsteinsson and J. A. Sethian, “A level set approach to a unified model for etching, deposition, and lithography. I. Algorithms and two-dimensional simulations,” J. Comp. Phys. 120, 128–144 (1995). [CrossRef]

19.

D. Adalsteinsson and J. A. Sethian, “A level set approach to a unified model for etching, deposition, and lithography II: three-dimensional simulations [integrated circuits],” J. Comp. Phys. 122, 348–366 (1995). [CrossRef]

20.

R. C. Rumpf, P. Srinivasan, and E. G. Johnson, “Modeling the fabrication of nano-optical structures,” Proc. SPIE 6110, 611004 (2006). [CrossRef]

21.

R. C. Rumpf and E. G. Johnson, “Comprehensive modeling of near-field nano-patterning,” Opt. Express 13, 7198 (2005). [CrossRef] [PubMed]

OCIS Codes
(050.7330) Diffraction and gratings : Volume gratings
(160.2900) Materials : Optical storage materials

ToC Category:
Materials

History
Original Manuscript: July 17, 2007
Revised Manuscript: August 9, 2007
Manuscript Accepted: August 13, 2007
Published: September 14, 2007

Citation
Raphi Dror, B. Sfez, Sh. Y. Goldin, and A. Cashingad, "Etching of photosensitive chalcogenide glasses:experiments and simulations," Opt. Express 15, 12539-12547 (2007)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-19-12539


Sort:  Year  |  Journal  |  Reset  

References

  1. K. Shimakawa, A. Kolobov, S. R. Elliott, "Photoinduced effects and metastability in amorphous semiconductors and insulators," Adv. Phys. 44, 475 (1995). [CrossRef]
  2. K. Tanaka, "Sub-Gap Photo-Induced Phenomena in Chalcogenide Glasses," in Photo-Induced Metastability in Amorphous Semiconductors, A.V. Kolobov, ed., (Wiley, Weinheim, 2003) pp 69. [CrossRef]
  3. A. Ozols and K. Shvarts, "Photosensitivity of amorphous semiconductor As-S and As-Se films under CW, nanosecond and picosecond laser irradiation," Cryst. Latt. Def. and Amorph. Mat. 17, 235-239 (1987).Q1
  4. G. Rosenblum, B. G. Sfez, Z. Kotler, V. Lyubin, and M. Klebanov, "Nonlinear optical effects in chalcogenide photoresists," Appl. Phys. Lett. 75, 3249 (1999). [CrossRef]
  5. V. M. Lyubin, A. M. Sedikh, N. N. Smirnova and V. P. Shilo, Microelectronica 18, 523 (1989).Q2
  6. A. Arsh, M. Klebanov, V. Lyubin, L. Shapiro, A. Feigel, M. Veinger, B. Sfez, "Glassy mAs2S3·nAs2Se3 photoresist films for interference laser lithography," Opt. Mater. 26,301-304 (2004). [CrossRef]
  7. M. Vlcek, P. J. S. Ewen and T. Wagner, "High efficiency diffraction gratings in As-S layers," J. Non-Cryst. Solids 227-230, 743 (1998). [CrossRef]
  8. A. V. Stronski, M. Vlcek, A. Sklenar, P. E. Shepeljavi, S. A. Kostyukevich, and T. Wagner, "Application of As40S60-xSex layers for high-efficiency grating production," J. Non-Cryst. Solids 266-269, 973 (2000). [CrossRef]
  9. S. Wong, M. Deubel, F. Pérez-Willard, S. John, G. A. Ozin, M. Wegener, and G. von Freymann, "Direct Laser Writing of Three- Dimensional Photonic Crystals with a Complete Photonic Bandgap in Chalcogenide Glasses," Adv. Mater. 18, 265-269 (2006). [CrossRef]
  10. A. Feigel, M. Veinger, B. Sfez, A. Arsh, M. Klebanov, V. Lyubin, "Two dimensional photonic band gap pattering in thin chalcogenide glassy films," Thin Solid Films 488,185-188 (2005). [CrossRef]
  11. A. Feigel, Z. Kotler, B. Sfez, A. Arsh, M. Klebanov and V. Lyubin, "Chalcogenide glass-based three-dimensional photonic crystals," Appl. Phys. Lett. 77, 3221 (2000). [CrossRef]
  12. A. Feigel, M. Veinger, B. Sfez, A. Arsh, M. Klebanov and V. Lyubin, "Three-dimensional simple cubic woodpile photonic crystals made from chalcogenide glasses," Appl. Phys. Lett. 83, 4480 (2003). [CrossRef]
  13. R. C. Rumpf and E. G. Johnson, "Fully three-dimensional modeling of the fabrication and behavior of photonic crystals formed by holographic lithography," J. Opt. Soc. Am. A 21, 1703-1713 (2004). [CrossRef]
  14. J. A. Sethian, "A fast marching level set method for monotonically advancing fronts," Proc. Nat. Acad. Sci. 93, 1591-1595 (1996). [CrossRef] [PubMed]
  15. J. A. Sethian, Level Set Methods and Fast Marching Methods, (Cambridge Univ. Press, 2nd ed., 1999).
  16. S. Osher and R. Fedkiw, Level Set Methods and Dynamic Implicit Surfaces (Springer, 2003).
  17. J. A. Sethian and D. Adalsteinsson, "Overview of level set methods for etching, deposition, and lithography development," IEEE Transactions on Semiconductor Devices 10, 167-184 (1997).Q3 [CrossRef]
  18. D. Adalsteinsson and J. A. Sethian, "A level set approach to a unified model for etching, deposition, and lithography. I. Algorithms and two-dimensional simulations," J. Comp. Phys. 120, 128-144 (1995).Q4 [CrossRef]
  19. D. Adalsteinsson and J. A. Sethian, "A level set approach to a unified model for etching, deposition, and lithography II: three-dimensional simulations [integrated circuits]," J. Comp. Phys. 122, 348-366 (1995).Q5 [CrossRef]
  20. R. C. Rumpf, P. Srinivasan, and E. G. Johnson, "Modeling the fabrication of nano-optical structures," Proc. SPIE 6110, 611004 (2006). [CrossRef]
  21. R. C. Rumpf and E. G. Johnson, "Comprehensive modeling of near-field nano-patterning," Opt. Express 13, 7198 (2005). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Supplementary Material


» Media 1: MOV (1426 KB)     

« Previous Article

OSA is a member of CrossRef.

CrossCheck Deposited