Multibandgap quantum well wafers by IR laser quantum

Multibandgap quantum well wafers by IR laser quantum well intermixing: simulation of the lateral resolution of the process O. Voznyy, R. Stanowski, J. J. Dubowski Department of Electrical and Computer Engineering Research Center for Nanofabrication and Nanocharacterization Université de Sherbrooke, Québec J 1 K 2 R 1 Canada

Outline 1. Motivation 2. Modeling heat distribution and photoluminescence (PL) in QW wafers 3. Temperature profiles induced in In. Ga. As/In. Ga. As. P wafers by moving laser beam 4. PL shift profiles 5. Summary 2

Motivation > Quantum well intermixing Multibandgap materials are needed for creation of photonic integrated circuits (lasers, modulators, waveguides, multicolor detectors etc. fabricated on same wafer) E 0 E 1 E 2 E 3 Quantum well intermixing (QWI) – interdiffusion of wells and barriers resulting in the change of the well width, potential barrier height and energy of confined states. 3

Motivation > Current state of the problem Simulations [1] predict transition region ~300μm using CW Nd: YAG laser irradiation (photoabsorbtion induced disordering) with a shadow mask [1]. Also, pulsed laser IR disordering (2 -step process) has been proposed (~2μm transtion region possible [2]). Our aim is to investigate Laser-RTA (annealing with a moving CW laser beam) as a flexible (1 -step process) and potentially cost-effective technique. [1] A. Mc. Kee, et. al. , IEEE J. Quantum Electron. , vol. 33, pp. 45– 55, Jan. 1997. [2] B. S. Ooi, et. al. IEEE J. Quantum Electron. , vol. 40, pp. 481– 490, May 2004 4

Motivation > Moving laser beam In previous work [3] array of 12 lines of intermixed Ga. As/Al. Ga. As QW material was successfully written with 5 cm/s, 0. 7 mm CW Nd: YAG laser beam in a 14 mm x 6 mm sample. This approach has the potential to write complex patterns of intermixed material. [3] J. J. Dubowski, et. al. , Proc. SPIE, 5339, (2004). Quantum well PL peak position measured across the sample irradiated with a fast scanning laser beam that was used to generate a 12 -line pattern. 5

Computation details > Finite Element Method simulations To find heat distribution in a wafer we used FEMLAB commercial software. Geometry is divided into small mesh elements with their own PDE parameters. Then the resulting system of PDEs is solved. Heat transfer PDE: Subdomain equation: Q - (k T) = Cp( T/ t) Boundary equation: k T=q 0 + h(Tinf – T) + εσ(Tamb 4 – T 4) For correct results temperature dependent thermal conductivity k and optical absorption α should be taken into account. 6

Computation details > Finding PLshift(D) 1. Take diffusion coefficient as parameter 2. Find concentration profile for given D and time 3. Find energy profiles for electrons and holes (take into account bandgaps, band offsets, bandgap bowing) 4. Solve Schrödinger equation, find energy levels and PL 5. Approximate results as some function D(PL shift) If T(t)=const (like with RTA): LD = – diffusion length. Otherwise one needs to solve numerically D assumed to be the same for different atomic species. 7

Computation details > Finding D(T) and PLshift(T, t) Compare simulations and experimental PLshift(Tanneal) data for the same annealing time, find D(Tanneal) Build Arrhenius plot ln. D(1/k. T) and find parameters for D=D 0 exp(-EA/k. T) Now we can find PL shift for any T and time. 8

Computation details > Laser power density and surface damage To achieve T needed for intermixing, different power needed for different beam diameters. For small diameters <0. 5 mm power densities become higher than surface damage threshold (>30 W/mm 2). 270 W/mm 2 700 W/mm 2 1500 W/mm 2 Needed power density can be reduced using background heating. 9

Computation details > Power density for moving beam With laser fast scanning (Laser-RTA) we can heat samples to same temperatures, with smaller beam diameters and avoid surface damage. Power needed to heat the wafer to TQWI increases a little, but fluence drops down significantly (shorter dwell time). TQWI 10

Computation details > Depth dependence For small beam diameters T drops down with depth very fast. In. P is transparent to Nd: YAG wavelength at RT, but Eg(In. P)=1. 165 e. V at 500°C, and α=104 -106 cm-1 at higher T. Thus, all the energy is absorbed on the surface and goes inside only by heat conduction. d=100μm d=12μm Depth, μm 100 50 50 100 Lateral, μm 0 50 100 Lateral, μm 11

Temperature profiles > Scanning speed and bg heating • For small samples slower speed results in raise of background temperature. • For big wafers heat dissipates faster and temperature profiles don’t depend on scanning speed (laser power is adjusted to achieve same Tsurface). • Background heating helps to achieve wanted T. 12

PL shift profiles > Temporal T behavior during scan To calculate PL shift profile for moving beam we need: • calculate concentration and energy profiles using given T(t) and D(T) at different distances from line center, • solve Schrödinger equation and find PL shift. 13

PL shift profiles > PL shift profile for moving beam Due to varying T(t), PL shift profile for moving beam differs from that of stationary beam, although temperature profiles are the same. PL shift profile shape doesn’t depend on Tmax. Higher temperatures reduce processing time significantly. 14

PL shift resolution and processing time • Processing time for 100 nm PL shift along one 2 -inch line assuming Tmax=1073 K (which requires 90 s to get the same PL shift with RTA). • Practical applications will require shifts < 50 nm. 15

Summary 1. 2. 3. 4. 4. 5. Irradiation with the moving CW Nd: YAG laser been has been investigated for selective area writing of the QWI material. For large size wafers (2 inch) temperature profiles don’t depend on scanning speed (assuming that beam power is adjusted to achieve the same Tmax). Processing time to achieve targeted PL (badgap) shifts depends on beam diameter and Tmax. To achieve reasonable processing time without loss in resolution a) QWs should be very close to surface, b) Tmax should be as high as allowed by material decomposition temperature Background heating can be used to further decrease processing time (especially for deep QWs) but decreasing also resolution. Lateral PL shift resolution of 5μm is feasible (In. Ga. As/In. Ga. As. P QW material system) with the 12μm beam Laser-RTA. Support Natural Sciences and Engineering Research Council of Canada (NSERC) Canada Research Chair (CRC) Program 16