Time-resolved observation of band-gap shrinking and electron-lattice thermalization within X-ray excited gallium arsenide.

Femtosecond X-ray irradiation of solids excites energetic photoelectrons that thermalize on a timescale of a few hundred femtoseconds. The thermalized electrons exchange energy with the lattice and heat it up. Experiments with X-ray free-electron lasers have unveiled so far the details of the electronic thermalization. In this work we show that the data on transient optical reflectivity measured in GaAs irradiated with femtosecond X-ray pulses can be used to follow electron-lattice relaxation up to a few tens of picoseconds. With a dedicated theoretical framework, we explain the so far unexplained reflectivity overshooting as a result of band-gap shrinking. We also obtain predictions for a timescale of electron-lattice thermalization, initiated by conduction band electrons in the temperature regime of a few eVs. The conduction and valence band carriers were then strongly non-isothermal. The presented scheme is of general applicability and can stimulate further studies of relaxation within X-ray excited narrow band-gap semiconductors.

Intense ultrashort radiation available at the free-electron-laser (FEL) facilities all over the world12345 enables structural studies of matter with atomic resolution67, and, with a pump-probe scheme, investigation of transient states of matter with a few femtosecond resolution. A pump pulse (FEL, optical, THz) initiates a transition which is probed with another laser pulse delayed by a certain time with respect to the pump pulse. It was already found during the first solid-state experiments with FELs that an intense FEL laser pulse excites many electrons within the irradiated solid. This leads to a transient change of the optical properties of the target which follows the evolution of the electron density. In particular, this allows to estimate the timescale of the delay between pump and probe pulse. Therefore, in the first89 and later experiments1011121314 measuring transient optical properties of solids the detected changes have been used to design non-destructive time-delay diagnostic and pulse centeracterization tools for FELs.

However, this experimental scheme also opens possibilities to extract information about the physical processes following the relaxation of the laser-excited materials. Similar techniques have been used in two-color experiments with optical lasers to extract information on the electron-phonon coupling coefficient15 or to measure the temperature-dependent electron-lattice thermalization time161718. Transient changes of transmission within X-ray irradiated GaAs were reported in ref. 19 and interpreted as a signature of non-equilibrium relaxation effects within the irradiated solid.

In this study we revisit the results of the experiments measuring the transient optical reflectivity in GaAs irradiated with femtosecond pulses of soft and hard X-ray radiation820. GaAs is a direct band-gap semiconductor. Its band-gap width of  eV at room temperature21 is slightly smaller than the energy of the optical photons at  nm ( eV). Thus, in addition to the free carrier absorption, optical photons may also trigger electron excitation from the valence to the conduction band, which then becomes the dominant channel of the energy absorption2122. The rate of this interband excitation depends on the band occupations and on the band-gap width. As the band-gap width in semiconductors is a function of the lattice temperature, the interband transition rate can reflect the transient changes of the lattice temperature.

Electronic excitation and relaxation processes following the FEL irradiation of solids are well understood232425. The arriving FEL photons excite electrons from the valence to the conduction-band and, at sufficiently high photon energies, also from the atomic inner shells. In light elements the inner-shell electron excitation is followed by an Auger decay, resulting in the emission of another electron into the conduction-band. Energetic electrons within the conduction band can excite more electrons from the valence band through the electron impact ionization. The electrons also interact with each other, exchanging energy. This leads to a fast thermalization of the conduction band electrons. Its timescale depends on the photon energy and fluence. At a typical few tens of femtoseconds FEL pulse durations, the electron thermalization timescales were reported to be fs both by theory2425 and experiment820. Similar timescales were found in experiments with optical excitation262728. A recent study29 showed that the electron thermalization following FEL pulses is mainly determined by the secondary electron cascading, which timescale increases with the increasing photon energy.

However, thermalization between the valence and conduction band may take longer time than the carrier thermalization within the conduction band. The relation between these two timescales depends on the ratio of the effective masses of holes and electrons. In GaAs holes are almost times heavier than the conduction band electrons. This allows to expect that the absorbed energy (brought by X-ray photons) will be initially shared among the light conduction-band electrons and then slowly transferred to the heavy holes. Conduction and valence-band carriers then remain for some time strongly non-isothermal. Similar carrier non-isothermality has been observed in3031 and recently in32.

Already after the appearance of first photoelectrons within the conduction band, the energy exchange between free carriers and lattice starts. In semiconductors the timescale of this process is of the order of picoseconds, which is typically much longer than the electron thermalization time1718. While the temperature of the lattice increases, the band gap shrinks. Based on the experimental data, phenomenological fits were constructed to describe the band-gap shrinking as a function of the lattice temperature213334. We consider here the case when the densities of excited carriers are low and do not lead to an additional band-gap shrinking35. After the free carriers thermalize, recombination processes begin to contribute. In GaAs both radiative and non-radiative recombination occur36. However, the typical recombination times for GaAs are of the order of ns 36, which is much longer than the electron-lattice thermalization time17. Figure 1 shows the timescales of the predominant excitation and relaxation processes in X-ray irradiated GaAs.

Figure 1

Timescales for predominant excitation and relaxation processes in X-ray irradiated GaAs.

In this work we show that the experimental data on the transient reflectivity allow to follow in time the relaxation of FEL excited GaAs on timescales up to a few tens of picoseconds. We follow the ideas from the optical measurements described in refs 15,17,18 however, the photon wavelengths considered in this work are in X-ray regime. They heat the solid to much higher temperatures (∼few eVs) than those considered in1718. This excitation results in the centeracteristic ‘reflectivity overshooting’, i.e., after the relaxation the reflectivity reaches an equilibrium value higher than its initial value before the excitation. This effect has not been explained so far37.

Our proposed theoretical framework uses rate equations to describe the evolution of electron distribution as a function of time. The rate equations are coupled with the two-temperature model38, describing the electron-lattice equilibration. The Drude model is applied to follow the transient reflectivity as a function of free-carrier density. Here, this model is extended beyond the free-carrier absorption framework34 and it also includes the predominant contribution from interband transitions. Comparing this model to the available data, we identify the reflectivity overshooting as a result of the band-gap shrinkage and obtain predictions for the electron-lattice thermalization time as a function of pulse fluence and photon energy. If more experimental data were available, the presented scheme could be applied for a quantitative study of excitation and relaxation times in X-ray excited GaAs and in other narrow-band-gap semiconductors for which .

Results

In experiments with X-ray FELs820 relative changes of transient reflectivity within GaAs have been measured, , as a function of the time delay between the FEL pump pulse and the optical probe pulse. In the experiment with soft X-rays by Gahl et al.8, FEL photons had the energy of 40 eV. In the experiment by Krupin et al.20, photons of 800 eV energy were used. Various FEL fluences have been applied for pumping the material. The FEL pulse duration did not exceed a few tens of fs. Optical probe pulses were of low intensity so that they did not cause any damage in GaAs. Their wavelengths were: (i) 800 nm (∼1.55 eV) in20 and (ii) both 800 nm and 400 nm ( eV) in8. The error of fluence estimation in the earlier experiment8 performed in 2008 is expected to be much larger than in20 performed in 2012 as the pulse diagnostics methods had been significantly improved.

Both experiments observed a centeracteristic shape of transient reflectivity curve, . On the timescale of a few hundreds of femtoseconds since FEL irradiation, one can see an initial ultrafast drop and rise of optical reflectivity (Figs 2,3 in ref. 8, Fig. 5 in ref. 20). These changes are followed by a slower relaxation of reflectivity towards an equilibrium value on a timescale of tens of ps. The transient reflectivity then overshoots the initial value. This is in contrast to the behaviour of irradiated insulators (e.g., silicon dioxide10 and silicon nitride20). In what follows we prove that this effect can occur only for semiconductors with a band gap slightly smaller than the energy of the probing optical photon. In this case, the lattice temperature-sensitive electronic transitions between valence and conduction band triggered by the optical photons influence significantly the transient optical properties of the material. As the excited carrier densities (maximally of the order of a few permilles of the initial valence electron density) are much too low to affect the band gap35, the transient optical properties then only reflect the band-gap shrinking due to the increasing temperature of the lattice.

Figure 2

Relative change of transient reflectivity in GaAs as a function of time measured in: (a) the experiment by Krupin et al.8 at FEL photon energy of 800 eV and optical probe of 800 nm, (b) the experiment by Gahl et al.20 at FEL photon energy of 40 eV and optical probes of 800 nm and 400 nm.

Experimental values (points) and theory predictions (lines) are compared for: (a) various FEL fluences: and 40 mJ/cm2, and (b) various wavelengths of probe pulse. The theory results in (b) are obtained for a FEL pulse fluence of  mJ/cm2.

With our modeling tool we performed simulations of the transient reflectivity changes within GaAs for the experimental conditions as described in refs 8,20. Our tool can describe optical properties within GaAs bulk, far from its surface. Therefore, all emitted electrons are assumed to stay within the material, i.e., they are emitted to the conduction band (below the continuum level of the material). The centerge neutrality is then preserved.

First, we estimated the final electron-hole density after FEL irradiation, knowing the pulse fluence and the photoabsorption cross section at a given photon energy. We then followed the increase of the electron-hole density until the maximal density was reached. Such behaviour is typical for FEL irradiated semiconductors as described, e.g., in2425.

After the maximal electron-hole density is reached, the system starts to relax. Electron-lattice thermalization and electron-hole recombination are the predominant relaxation channels. In GaAs the latter one contributes less significantly, as the typical recombination timescales are of the order of ns36. Electron-lattice thermalization in GaAs is expected to act on ps time scales1718. The resulting lattice heating leads to a shrinkage of the band gap39, when compared with its initial width of eV at 300 K. In what follows we demonstrate that this affects the interband photoabsorption, which results in the experimentally observed ‘overshooting’ of the transient reflectivity.

The extended Drude model including the interband contribution is used to calculate the transient reflectivity change from the dielectric function, , which is the function of optical coefficients, and . See Supplemental Material at [URL will be inserted by publisher] for more details on the model. The dielectric function is parametrized as in ref. 34

where describes all contributions to the dielectric function beyond the free-carrier absorption. Here, , where describes the contribution from the transition between valence and conduction bands, using the interband absorption coefficient for a direct interband transition, , parametrized as in ref. 34 (Eq. (5.31) therein). The absorption coefficient scales with the photon energy , as   Band-gap shrinking is described with the phenomenological relation from39. The interband absorption coefficient also contains the matrix element, , which couples states with the same electron wave vector in the valence and conduction bands. We parametrize it, using the measured absorption coefficients for GaAs from ref. 22 (Fig. 3 therein). The time is the electron (hole) collision time and the frequency is the plasma frequency for electrons (holes). Carrier density determines the plasma frequency of the material . The frequency is the photon frequency and is its corresponding wavelength.

The average electron collision time is fitted in order to match the minimum of curve. The accuracy of the fit depends on the accuracy with which we can estimate the minimum reflectivity. This is determined by the time resolution of the data. An improvement of the experimental temporal resolution would improve the accuracy of all model fits performed. The average hole collision time can be estimated from the electron one, using mass scaling relation for electron and hole collision frequencies40. The initial value of the lattice temperature is 300 K. Lattice temperature does not change much during the first 100–200 fs after the FEL irradiation (during electron thermalization) and so the band-gap width does not change either.

Rate equations describe the changes of electron-hole density due to ionization and recombination processes. Heat capacities of the free-electron gas from ref. 40 and of the lattice from ref. 21 are used in the temperature equation, describing the exchange of thermal energy between free electrons and lattice. Fitted (i.e., iteratively adjusted) parameters are: (i) the free-electron temperature at the minimum of curve, , and (ii) the thermalization time, . Plots (Fig. 2) show the theoretical calculations (solid lines) compared to the experimental data (points).

The ultrafast drop and rise of optical reflectivity is due to the production of electron-hole pairs in the irradiated semiconductors as described in ref. 10. The created free carriers contribute to the optical properties of the material. Their contribution can be described by the Drude model.

The theoretical calculations shown in Fig. 2a were performed with experimentally measured fluence values. The Drude model could not be applied at the experimentally measured fluence value by Gahl et al. that lead to an overcritical electron density. However, at the pulse fluences close and larger than the one leading to the critical electron density, the minimum of curve does not change (Fig. 3b in8). We then obtained results presented in Fig. 2b, using a fluence, corresponding to a subcritical electron density, at which the Drude model is still applicable. The predictions for 400 nm are obtained for the same FEL pulse parameters. Note that they cannot be directly compared to the experimental curve from Fig. 2b in8: according to Fig. 3a therein, it corresponds to a different set of FEL pump parameters.

At a 400 nm probe-pulse the interband absorption coefficient is about 100 times larger than at 800 nm2122. The optical coefficients are sensitive to the value of absorption coefficient which, according to various estimations, may differ by a factor of 2–32122. This strongly influences the accuracy of theory predictions (Fig. 2). Reflectivity overshooting does not show up for 400 nm as the large absorption coefficient is then not sensitive to the small change caused by energy shift due to the band-gap shrinking ().

As we have discussed above, the centeracteristic time of the overshooting directly follows the timescale of the band-gap shrinking. The experimental data on reflectivity change can then be used to determine the electron-lattice thermalization time. Experimental studies of electron-phonon coupling and electron-lattice thermalization time have been extensively performed in metals154142434445, semiconductors171846, superconductors47, and warm-dense matter484950. Theoretical studies are frequently performed within the framework of two-temperature model38435152535455. New developments can be found in565758.

Here, we apply our model to the experimental data to extract information on the electron-lattice thermalization time in GaAs. In Table 1 we list the parameters iteratively adjusted to obtain the predictions in Fig. 2a, i.e., the thermalization time and the free electron temperature at the minimum of curve . To clarify, our numerical scheme calculates transient reflectivity, adjusting iteratively: (i) the thermalization time, and (ii) the free electron temperature at the minimum of curve, so as to obtain the best agreement of our prediction with the experimental curve for times after the minimum. The average relative distance, , between the model predictions and experimental data is the corresponding error measure. It is defined as , where and N is the number of the available experimental points. The values of thermalization time and electron temperature presented in Table 1 were obtained for the smallest achieved which was % for all presented parameter values. These values lay within the range reported in ref. 17. The product of collision time and photon frequency is 1 for  nm and 2 for  nm. Note the increase of with fluence, and the corresponding decrease of the thermalization time with the electron temperature17.

Table 1

Parameters used to obtain predictions in Fig. 2a: thermalization time () and the free electron temperature at the minimum of curve ().

40	2.0	2.8
20	2.5	2.2
10	3.0	1.6

The values of and for Fig. 2b (experiment by Gahl et al.) at the pulse fluence of mJ/cm2 are ps and  eV, respectively, for both  nm and  nm.

Discussion

The estimated electron temperature at the end of the electron thermalization (the minimum of curve) is much higher ( eV) than the temperature that we would obtain assuming a full thermalization of valence and conduction bands at this time ( eV). As discussed above, this results from the delayed thermalization between electron and hole systems due to their large mass difference (). A crude estimation of the average kinetic energy of a free electron within the conduction band is , where is the average pair creation energy21, yields a much higher average electron temperature of eV. This is a limiting case for which we assumed that the created hole has a negligible kinetic energy. As electron-lattice thermalization starts before the electron thermalization is completed, the estimated should be interpreted as kinetic temperature which is slightly larger than . We consider the high electron temperature as an indication for a strong non-isothermality of conduction and valence-band carriers that maintains up to a few hundred femtoseconds since exposure to the FEL pulse and even after the carriers within each band are already thermalized. To compare, in31 the heavy hole thermalization time was found to be fs within GaAs at room temperature. In30 the equilibration of electron and hole temperatures was predicted to occur within 10 ps since the exposure to an optical pulse in a photoexcited bulk semiconductor ().

In addition to the processes discussed above, thermal diffusion of electrons from the interaction region can influence the thermalization process. During the diffusion a fraction of electrons is slowly leaving the interaction region59, carrying a fraction of the pulse-absorbed energy outside. The contribution of diffusion is related to the attenuation length of photons which defines the size of the interaction region. For 40 eV photons, it is microns, for eV photons, it is microns in GaAs60. Created electrons can then escape to the radiation-unaffected part of GaAs bulk, ‘disappearing’ from the interaction region. We can therefore expect that the diffusion should affect much more the data obtained with eV photons than those obtained with eV photons. However, in both cases the estimated thermalization times and initial electron temperatures are of the same order, indicating that the effect of thermal diffusion is negligible for our data analysis.

To sum up, in this study we presented a theoretical model that follows the relaxation of X-ray-laser excited GaAs on a few tens of picoseconds timescale. In particular, it explains the reflectivity overshooting observed in820 as an effect of thermal band-gap shrinking during lattice heating. The model includes interband transitions and uses rate and two-temperature equations to follow the relaxation of FEL excited GaAs. Model results show that the reflectivity overshooting is a signature of electron-lattice thermalization due to electron-phonon coupling and can be used to determine the electron-lattice thermalization timescale. The presented scheme is of general applicability. We expect it to inspire dedicated, quantitative studies of relaxation times in X-ray excited GaAs and other narrow band-gap semiconductors for which the thermalization and recombination times fulfill the condition .

Methods

The applied theoretical model uses rate equations to describe the evolution of free-carrier densities within irradiated GaAs bulk as a function of time. The rate equations are coupled with the two-temperature model38, describing the electron-lattice equilibration. The Drude model is applied to follow the transient reflectivity as a function of carrier densities. Here, this model is extended beyond the free-carrier absorption framework34, and it also includes the predominant contribution from interband transitions. More details on the model are given in the Supplemental Material attached at [URL will be inserted by publisher]. Results of the model are compared to the existing experimental data on transient optical properties of GaAs taken from refs 8,20.

Additional Information

How to cite this article: Ziaja, B. et al. Time-resolved observation of band-gap shrinking and electron-lattice thermalization within X-ray excited gallium arsenide. Sci. Rep. 5, 18068; doi: 10.1038/srep18068 (2015).