Abstract
We have systematically investigated the spatial and temporal dynamics of crystallization that occur in the phase-change material Ge2Sb2Te5 upon irradiation with an intense terahertz (THz) pulse. THz-pump-optical-probe spectroscopy revealed that Zener tunneling induces a nonlinear increase in the conductivity of the crystalline phase. This fact causes the large enhancement of electric field associated with the THz pulses only at the edge of the crystallized area. The electric field concentrating in this area causes a temperature increase via Joule heating, which in turn leads to nanometer-scale crystal growth parallel to the field and the formation of filamentary conductive domains across the sample.
Original language | English |
---|---|
Article number | 165702 |
Journal | Physical review letters |
Volume | 121 |
Issue number | 16 |
DOIs | |
Publication status | Published - 2018 Oct 19 |
Externally published | Yes |
ASJC Scopus subject areas
- Physics and Astronomy(all)
Access to Document
Fingerprint Dive into the research topics of 'Zener Tunneling Breakdown in Phase-Change Materials Revealed by Intense Terahertz Pulses'. Together they form a unique fingerprint.
Cite this
- APA
- Standard
- Harvard
- Vancouver
- Author
- BIBTEX
- RIS
Zener Tunneling Breakdown in Phase-Change Materials Revealed by Intense Terahertz Pulses. / Sanari, Yasuyuki; Tachizaki, Takehiro; Saito, Yuta; Makino, Kotaro; Fons, Paul; Kolobov, Alexander V.; Tominaga, Junji; Tanaka, Koichiro; Kanemitsu, Yoshihiko; Hase, Muneaki; Hirori, Hideki.
In: Physical review letters, Vol. 121, No. 16, 165702, 19.10.2018.Research output: Contribution to journal › Article › peer-review
}
TY - JOUR
T1 - Zener Tunneling Breakdown in Phase-Change Materials Revealed by Intense Terahertz Pulses
AU - Sanari, Yasuyuki
AU - Tachizaki, Takehiro
AU - Saito, Yuta
AU - Makino, Kotaro
AU - Fons, Paul
AU - Kolobov, Alexander V.
AU - Tominaga, Junji
AU - Tanaka, Koichiro
AU - Kanemitsu, Yoshihiko
AU - Hase, Muneaki
AU - Hirori, Hideki
N1 - Funding Information: This study was supported by PRESTO (No. JPMJPR1427) and CREST (No. JPMJCR14F1) Grants from JST and the KAKENHI Grant No. 17H06228 from JSPS. Funding Information: Sanari Yasuyuki 1 Tachizaki Takehiro 2,3 Saito Yuta 4 Makino Kotaro 4 Fons Paul 4 Kolobov Alexander V. 4 ,* Tominaga Junji 4 Tanaka Koichiro 3,5 Kanemitsu Yoshihiko 1 Hase Muneaki 6 ,† Hirori Hideki 1,3 ,‡ Institute for Chemical Research, 1 Kyoto University , Uji, Kyoto 611-0011, Japan Department of Optics and Imaging Science and Technology, 2 Tokai University , Hiratsuka, Kanagawa 259-1292, Japan Institute for Integrated Cell-Material Sciences (iCeMS), 3 Kyoto University , Kyoto, Kyoto 606-8501, Japan Nanoelectronics Research Institute, 4 National Institute of Advanced Industrial Science and Technology , Tsukuba Central 5, 1-1-1 Higashi, Tsukuba 305-8565, Japan Department of Physics, Graduate School of Science, 5 Kyoto University , Kyoto, Kyoto 606-8502, Japan Division of Applied Physics, Faculty of Pure and Applied Sciences, 6 University of Tsukuba , 1-1-1 Tennodai, Tsukuba 305-8573, Japan * Present address: Department of Physical Electronics, Herzen State Pedagogical University, St Petersburg 191186, Russia. † mhase@bk.tsukuba.ac.jp ‡ hirori@scl.kyoto-u.ac.jp 19 October 2018 19 October 2018 121 16 165702 29 June 2018 28 March 2018 Published by the American Physical Society 2018 authors Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. We have systematically investigated the spatial and temporal dynamics of crystallization that occur in the phase-change material Ge 2 Sb 2 Te 5 upon irradiation with an intense terahertz (THz) pulse. THz-pump–optical-probe spectroscopy revealed that Zener tunneling induces a nonlinear increase in the conductivity of the crystalline phase. This fact causes the large enhancement of electric field associated with the THz pulses only at the edge of the crystallized area. The electric field concentrating in this area causes a temperature increase via Joule heating, which in turn leads to nanometer-scale crystal growth parallel to the field and the formation of filamentary conductive domains across the sample. Precursory Research for Embryonic Science and Technology 10.13039/501100009023 JPMJPR1427 Core Research for Evolutional Science and Technology 10.13039/501100003382 JPMJCR14F1 Japan Science and Technology Agency 10.13039/501100002241 17H06228 Japan Society for the Promotion of Science 10.13039/501100001691 Ge 2 Sb 2 Te 5 (GST) is one of the best performing chalcogenide alloys for phase-change applications, and its use in optical recording media is due to its robust optical contrast upon optical excitation inducing a phase change between amorphous and crystalline phases [1–3] . Because the phase change also enables resistive switching, which can be induced by nanosecond electrical pulses, GST is also suitable for electrical nonvolatile memory applications [4,5] . It has been argued that GST sandwiched between two electrodes reaches its crystallization temperature through Joule heating. However, the electric field effects on resistive switching are still unclear, because many physical processes may simultaneously participate with both thermal and electric field effects. These processes range from simple atomic structural changes [2,6–11] to complex combinations of Joule heating-induced crystallization [4,12–14] , nanometer-scale crystal growth or amorphization [15,16] , nonlinear increases in electrical conductivity [17,18] , and formation of filamentary paths across the active material between the electrodes [19,20] . Additionally, the phase changes induced by subnanosecond duration electrical pulses lead to a dome-shaped crystalline area [21,22] , which is considered to form a result of isotropic heat diffusion governing the spatial expansion of the crystal area. The experimental conditions obscure the complex mechanisms of resistive switching [4,13,14] , where the effects of both heating and electrical fields contribute to the phase change. To understand the relevant carrier dynamics, it is necessary to separate the contributions of these components by performing a sophisticated experiment. Recent advances in generating intense terahertz (THz) pulses have yielded new strategies for manipulating the electron and lattice degrees of freedom on picosecond timescales [23,24] . The use of ultrashort THz pulses helps to suppress heat diffusion and may enable the study of crystallization mechanisms in which the lattice temperature exceeds the crystallization temperature on picosecond timescales. The present study used THz pulses of a few picoseconds to induce crystallization in amorphous GST. Pulses with fields up to a few MV / cm were generated through local field enhancement using a gold dipole antenna; the resulting fields were an order of magnitude higher than those in previous reports on electrical conductivity [25] . THz-pump–optical-probe spectroscopy combined with microscopy and micro-Raman measurements enabled a spatially resolved investigation of the effects of the electric field on the amorphous and crystalline GST areas within the antenna gap. The experimental configuration is shown in Fig. 1(a) [26] . THz pulses were generated with the tilted pulse-front technique using a LiNbO 3 crystal and the setup described in Refs. [27,28] . The amorphous GST sample (40 nm thick) was deposited by sputtering a GST alloy target onto a (100) Si substrate at room temperature. A thin ZnS - SiO 2 layer (20 nm) was grown on top of the film in the same sputtering chamber to prevent oxidation. Figure 1(b) shows the temporal profile of the incident THz electric field. The field was enhanced by gold antennas fabricated on the sample surface. Figure 1(c) plots finite difference time domain (FDTD) calculations of the spatial distribution of the enhanced electric field and its direction (white arrows) inside the sample located in the 5 - μ m gap between the antennas and at a depth of 10 nm from the interface between the ZnS - SiO 2 and GST layers. The gap of the present sample is 5 times larger than that used in a similar previous work [18] , which allows us to trace the crystal growth much easier. The color bar shows the ratio between the maximum electric field of the incident pulse E o and the electric field induced inside the sample E i . The enhancement factors were about five in the center of the gap and more than 15 near the antenna corners, producing E i values ranging from 1 to 3 MV / cm . 1 10.1103/PhysRevLett.121.165702.f1 FIG. 1. (a) Experimental setup. The Au antenna enhances the THz electric field over the 5 - μ m gap. (b) Measured temporal profile of the pulse. (c) Electric field enhancement in the gap. White arrows indicate the electric field directions. (d)–(g) Microscope images of spatial reflectivity changes ( Δ R / R ) after THz irradiation as a function of pulse number. The images were obtained by averaging over 450 pulses (225 pulses before and after the corresponding center pulse) and dividing the result with the data obtained without THz excitation. A THz pulse with a peak field of E o = 175 kV / cm (500 Hz repetition rate) was used to illuminate the sample, and the corresponding changes in reflectivity Δ R / R were monitored by an optical microscope equipped with a charge coupled device (CCD) [Figs. 1(d)–1(g) ]. The white areas in the images correspond to those with increased reflectivity in the visible light region. These figures show that changes begin to appear after about 10 000 pulses, and an altered region with an elongated shape develops during the subsequent 5000 pulses. Because this growth process coincides with the direction of the electric field in Fig. 1(c) , this phenomenon is due to the THz pulse irradiation [31] . We used Raman spectroscopy to identify the structure of the regions showing enhanced reflectivity. Figure 2(a) compares the spectrum of the region in Fig. 1(g) that became completely white (red line) with the spectra for amorphous GST (blue line) and crystalline GST (black line); the latter crystalline sample was prepared by optical-pulse irradiation (about 1 mJ / cm 2 , λ = 800 nm ) of the amorphous phase. The characteristic peaks of amorphous GST are at 130 and 148 cm − 1 ( A 1 modes), while those of crystalline GST are at 105 and 160 cm − 1 ( E g and A 1 g modes, respectively) [32,33] . The characteristics of the spectrum from the region altered by the THz pulses are identical to those of the crystal. Figure 2(b) shows the corresponding two-dimensional map of Fig. 1(g) , plotting the intensity ratio between the Raman peaks at 105 and 148 cm − 1 . The clear difference between the amorphous and crystalline phases means that the crystalline phase coincides with the region where the reflectivity change occurred. The growth speed in Figs. 1(d)–1(g) was estimated to be a few nanometers per pulse ( 5 μ m for 5000 pulses), and it proceeded along the field direction. 2 10.1103/PhysRevLett.121.165702.f2 FIG. 2. (a) Comparison of Raman spectra for point B in Fig. 1(g) (red line), crystalline GST (black line), and amorphous GST (blue line). Arrows indicate the mode positions; the curves are offset for clarity. (b) Spatial map of intensity ratio of Raman signals at (i) 105 and (ii) 148 cm − 1 . Since Au has a constant signal over a broad region, the ratio is about unity. To understand the processes occurring on an ultrafast timescale due to the illumination with THz pulses and the resulting elongated crystal growth observed in Figs. 1(d)–1(g) , the transient reflectivity change Δ r / r within the crystallized area was measured using a pump-probe technique. The probe light pulse had a wavelength of 800 nm and a beam spot diameter of ≈ 1 μ m [34] . The Δ r / r signal obtained from the crystalline phase at point B in Fig. 1(g) is shown in Fig. 3(a) for three different incident electric field intensities. To prevent changes in the sample structure, the measurements were carried out using field strengths ( E o ≤ 150 kV / cm ) well below those required for crystallization. Note that there was no difference between the signals obtained in measurements using the ramp-up or ramp-down in THz field, as shown in the inset of Fig. 3(a) ; i.e., the process was reversible. 3 10.1103/PhysRevLett.121.165702.f3 FIG. 3. (a) Time-resolved reflectivity change for point B [Fig. 1(g) ] for different THz electric field strengths (black, E o = 50 kV / cm ; blue, E o = 100 kV / cm ; red, E o = 150 kV / cm ) and optical fluence of 0.12 mJ / cm 2 (yellow). The inset shows the data for ramp-up or ramp-down in THz field. (b) THz intensity dependence of the peak reflectivity modulation Δ r max / r at the time delay of 25 ps in (a). No change was observed for the amorphous region because the electric field enhancement at the probe spot is significantly lower than the electric field enhancement near the antenna. While the latter reaches the crystallization temperature, the temperature of amorphous region at the probe spot increases only by about 3 K [35,38] . The scale of Δ T L is on the right. The data were measured by varying E o as follows: 150, 120, 100, 80, 50, 90, 110, 130, 60, 70, and 140 kV / cm . Because the data line up on a single curve, the sample condition was not altered during the measurements and the repeatability is confirmed. (c) Electric field enhancement in the 5 - μ m gap during the actual crystal growth indicated in Fig. 1(e) . The electric field is dramatically enhanced at the edge of the crystalline area with a conductivity of 1 × 10 3 S / cm . For the largest field value, 150 kV / cm [in Fig. 3(a) , red curve], the reflectivity dropped rapidly within the first few picoseconds and then recovered. This first component of the signal was linearly proportional to the square of the electric field (intensity) and thus shows the increase in absorption due to the electroabsorption effect [39] . A second drop in reflectivity with a delayed onset occurred at a delay of ≈ 25 ps , and the reflectively subsequently rose gently with a decay time of ≈ 500 ps (the heat diffusion length estimated from this trend was 4.5 nm [40] ). The slow component corresponds to an increase in lattice temperature due to electron-phonon coupling, suggesting that the intensity at t = 25 ps is proportional to the total temperature rise Δ T L . This rise time of 25 ps was determined from the electron-phonon coupling coefficient G , which is 1.7 × 10 9 W / ( cm 3 K ) according to Fig. 3(a) [44] . The temperature rise Δ T L can be expressed in terms of the electrical conductivity σ [the Joule heat σ E i 2 generated by the internal electric field E i ( ∝ E o ) ] and the specific heat of the lattice C L , which is 1.25 J / ( cm 3 K ) for the crystalline phase and 1.33 J / ( cm 3 K ) for the amorphous phase [45] : Δ T L = 1 C L ∫ σ ( t ) E i ( t ) 2 d t . (1) The relationship between Δ T L and Joule heating can be examined from that of the amplitude of the reflectivity change at 25 ps, Δ r max / r (which is proportional to Δ T L ), and E o 2 [Fig. 3(b) ]. If the electrical conductivity is constant, Eq. (1) predicts that Δ T L should increase linearly with E o 2 (dashed line), which is inconsistent with the experimental data. This suggests that the electrical conductivity of the crystalline phase increases upon strong THz excitation. Upon illumination with a THz pulse, the initial crystallization starts around the antenna edge solely through the Joule heat generated by the large electric field in this region, which was confirmed experimentally (crystallization usually started near the corners of the antenna in our experiment) [31,35] . Nonetheless, the subsequent stepwise elongated crystal growth parallel to the field cannot be explained by the heat distribution that is obtained from Fig. 1(c) . The increase in σ of the crystalline area probably originates from an increase in the electron density N . In particular, this suggests that Zener tunneling is responsible for the electron excitation in the crystalline phase; i.e., when a strong electric field E i is applied to a semiconductor with a small band gap, the band is modulated and electrons tunnel directly from the valence band into the conduction band. The number of electrons excited into the conduction band per unit volume and time can be expressed as [46–48] N Zener ( E i ( t ) ) = e 2 E i ( t ) 2 m μ 1 / 2 18 π ℏ 2 E G 1 / 2 exp ( − π m μ 1 / 2 E G 3 / 2 2 ℏ e | E i ( t ) | ) , (2) where ℏ is the Dirac constant, e the elementary charge, m μ the reduced mass [expressed as m μ = m e m h / ( m e + m h ) using the effective masses m e and m h of electrons and holes], and E G the band gap energy. We assumed that Eq. (2) holds for all time periods and determined the electron density in the crystal N for a homogeneous distribution of additional electrons throughout the crystal volume at all time steps of the FDTD calculation (see Supplemental Material [48] ). Here, each cell in the FDTD model employs the local value of E i in order to account for the distribution of the electric field enhancement in the crystalline phase induced by E o . The electrical conductivity σ can be determined by assuming that it is proportional to N in accordance with the Drude model. In turn, Δ T L can be calculated by inserting the obtained σ and E i at the probe point into Eq. (1) . Setting m e = 0.4 m [49] , m h = 0.2 m [50] ( m is the electron rest mass), and E G = 0.24 eV in Eqs. (1) and (2) , we obtained a curve that reproduces the experimental results well [red solid line in Fig. 3(b) ]. The relationship between Δ r max / r and N can be quantitatively examined using the THz-pump data with the signal obtained for 400-nm-pulse excitation (the yellow curve in Fig. 3(a) ; the calculation yields an optically excited electron density of N opt = 2 × 10 20 cm − 3 [57] ). Δ r max / r was 2.6 × 10 − 2 at 25 ps for E o = 150 kV / cm [Fig. 3(a) , red curve], which is 10 times larger than the Δ r max / r for optical excitation ( ≈ 2.8 × 10 − 3 ), implying that the carrier density generated by the THz pulse is on the order of 10 21 cm − 3 . We obtained N = 0.67 × 10 21 cm − 3 from a calculation based on Eq. (2) for E o = 150 kV / cm , which is quantitatively consistent with the experimental result ( 10 21 cm − 3 ). Thus, for the electric field range where crystal growth occurs ( E o = 175 kV / cm ), N is larger than 10 21 cm − 3 and, thus, σ ≥ 1 × 10 3 S / cm . To clarify the distribution of the effective electric field during crystal growth, we used the above value of σ = 1 × 10 3 S / cm for the crystalline region and performed a FDTD calculation. The results of the calculation [Fig. 3(c) ] show that a highly concentrated field exists at the upper edge of the crystalline area, with values about 20 times that of the electric field acting on the amorphous region. Since the enhancement at the upper crystal edge is strong enough to exceed the crystallization temperature ( Δ T L = 130 K [35] ), we infer that crystal growth occurs along the electric field directions [see Figs. 1(e) and 1(f) ] as a consequence of this concentration. Furthermore, if Zener tunneling is absent (i.e., for the case of σ = 100 S / cm ), the electric field enhancement factor at the edge of the crystalline area is only 8 times and Δ T L = 26 K , which suggests that the concentration of the electric field due to Zener tunneling induces the crystal growth in the region between the antennas. To verify the dimensionality of the crystallization in Fig. 1 , we compared the pulse number dependence of the crystal volume fraction f ( P ) with the Johnson-Mehl-Avrami-Kolmogorov (JMAK) model [59,60] . We parametrized f ( P ) with the reflectivity change Δ R / R at points A , B , and C in Fig. 1(g) and plotted the corresponding normalized crystal volume fraction f as a function of P in Fig. 4 . When crystallization starts after P i pulses at point i , the JMAK model leads to an Avrami equation of the form f ( P ) = 1 − exp [ - K ( P − P i ) n ] , with K and n being constant because the time for crystallization per pulse is constant ( P is the pulse number). Here, n is the Avrami constant from which one can estimate the dimensionality of crystal growth. If 1 < n < 2 , the crystal grows one-dimensionally, while for 2 < n < 3 , it grows two-dimensionally [61] . Crystallization begins at P A , P B , and P C , which define the time required for the crystalline nuclei to reach the critical radius [62] , and Δ R / R saturates at P f . f ( P ) was fitted by varying n ; the fitting confirmed that n i changes dramatically when a threshold, defined by P m , is reached (for P = P m , the crystalized area in the image connects the upper and lower Au antennas). The region before P m can be classified as one-dimensional growth, since growth proceeds with n < 2 , as illustrated in the inset of Fig. 4 . On the other hand, n a v g ( = 2.48 ) in the region after P m is consistent with crystallization driven by pulsed laser illumination ( n = 2.5 [61] ), indicating that it proceeds via isotropic heating of the entire crystal volume. 4 10.1103/PhysRevLett.121.165702.f4 FIG. 4. THz pulse number dependence of the crystal volume fraction f ( P ) obtained experimentally at the different locations marked as points A (red square), B (green square), and C (blue square) in Fig. 1(g) . The experimental f ( P ) is obtained by dividing the reflectivity Δ R ( P ) by the saturated maximum value: Δ R max = Δ R ( P = 20000 ) . The red, green, and blue solid curves are fits of the Avrami equation to the experimental values at points A , B , and C ranging from P A , P B , and P C to P m . The black curve is a fit using the averaged data measured at the points A , B , C for P > P m . P A = 1.1 × 10 4 , P B = 1.2 × 10 4 , P C = 1.4 × 10 4 , P av g = 1.5 × 10 4 , K A = 2.0 × 10 − 4 , K B = 6.6 × 10 − 5 , K C = 1.5 × 10 − 6 , K av g = 1.1 × 10 − 8 , n A = 0.90 , n B = 1.05 , n C = 1.67 , n av g = 2.48 . Inset: Mechanism of the crystal growth under applied electric field. In conclusion, the observed nonlinear increase in the conductivity of the crystalized region is a result of Zener tunneling; the tunneling leads to local heating which causes a one-dimensional crystal growth. Besides revealing the carrier dynamics in GST under the influence of electric fields that are relevant to crystal growth, our results indicate that the observed growth rate was a few nm per pulse, which provides a sophisticated method for miniaturizing memory devices down to the nanometer scale. Finally, we note that the role of the nucleation in the crystallization process is still a matter of debate. As such, the characterization of the atomic structure is essential to clarify the crystal growth mechanism under high electric fields. This study was supported by PRESTO (No. JPMJPR1427) and CREST (No. JPMJCR14F1) Grants from JST and the KAKENHI Grant No. 17H06228 from JSPS. We are grateful to Naoki Yokoyama and Hiroyuki Akinaga for valuable comments. [1] 1 S. R. Ovshinsky , Phys. Rev. Lett. 21 , 1450 ( 1968 ). PRLTAO 0031-9007 10.1103/PhysRevLett.21.1450 [2] 2 A. V. Kolobov , P. Fons , A. I. Frenkel , A. L. Ankudinov , J. Tominaga , and T. Uruga , Nat. Mater. 3 , 703 ( 2004 ). NMAACR 1476-1122 10.1038/nmat1215 [3] 3 M. Wuttig and N. Yamada , Nat. Mater. 6 , 824 ( 2007 ). NMAACR 1476-1122 10.1038/nmat2009 [4] 4 M. Boniardi , D. Ielmini , I. Tortorelli , A. Redaelli , A. Pirovano , M. Allegra , M. Magistretti , C. Bresolin , D. Erbetta , A. Modelli , E. Varesi , F. Pellizzer , A. L. Lacaita , and R. Bez , Solid State Electron. 58 , 11 ( 2011 ). SSELA5 0038-1101 10.1016/j.sse.2010.11.033 [5] 5 R. E. Simpson , P. Fons , A. V. Kolobov , T. Fukaya , M. Krbal , T. Yagi , and J. Tominaga , Nat. Nanotechnol. 6 , 501 ( 2011 ). NNAABX 1748-3387 10.1038/nnano.2011.96 [6] 6 N. Yamada , E. Ohno , K. Nishiuchi , N. Akahira , and M. Takao , J. Appl. Phys. 69 , 2849 ( 1991 ). JAPIAU 0021-8979 10.1063/1.348620 [7] 7 Z. Sun , J. Zhou , and R. Ahuja , Phys. Rev. Lett. 96 , 055507 ( 2006 ). PRLTAO 0031-9007 10.1103/PhysRevLett.96.055507 [8] 8 S.-W. Nam , H.-S. Chung , Y. C. Lo , L. Qi , J. Li , Y. Lu , A. T. C. Johnson , Y. Jung , P. Nukala , and R. Agarwal , Science 336 , 1561 ( 2012 ). SCIEAS 0036-8075 10.1126/science.1220119 [9] 9 T. H. Lee and S. R. Elliott , Phys. Rev. Lett. 107 , 145702 ( 2011 ). PRLTAO 0031-9007 10.1103/PhysRevLett.107.145702 [10] 10 I. Ronneberger , W. Zhang , H. Eshet , and R. Mazzarello , Adv. Funct. Mater. 25 , 6407 ( 2015 ). AFMDC6 1616-301X 10.1002/adfm.201500849 [11] 11 J. Kalikka , J. Akola , and R. O. Jones , Phys. Rev. B 94 , 134105 ( 2016 ). PRBMDO 2469-9950 10.1103/PhysRevB.94.134105 [12] 12 J. Orava , A. L. Greer , B. Gholipour , D. W. Hewak , and C. E. Smith , Nat. Mater. 11 , 279 ( 2012 ). NMAACR 1476-1122 10.1038/nmat3275 [13] 13 R. Jeyasingh , S. W. Fong , J. Lee , Z. Li , K.-W. Chang , D. Mantegazza , M. Asheghi , K. E. Goodson , and H.-S. P. Wong , Nano Lett. 14 , 3419 ( 2014 ). NALEFD 1530-6984 10.1021/nl500940z [14] 14 A. Sebastian , M. L. Gallo , and D. Krebs , Nat. Commun. 5 , 4314 ( 2014 ). NCAOBW 2041-1723 10.1038/ncomms5314 [15] 15 M. Salinga , E. Carria , A. Kaldenbach , M. Bornhöfft , J. Benke , J. Mayer , and M. Wuttig , Nat. Commun. 4 , 2371 ( 2013 ). NCAOBW 2041-1723 10.1038/ncomms3371 [16] 16 B. Casarin , A. Caretta , B. Chen , B. J. Kooi , R. Ciprian , F. Parmigiani , and M. Malvestuto , Nanoscale 10 , 16574 ( 2018 ). NANOHL 2040-3364 10.1039/C8NR04350G [17] 17 D. Ielmini and Y. Zhang , J. Appl. Phys. 102 , 054517 ( 2007 ). JAPIAU 0021-8979 10.1063/1.2773688 [18] 18 P. Zalden , M. J. Shu , F. Chen , X. Wu , Y. Zhu , H. Wen , S. Johnston , Z.-X. Shen , P. Landreman , M. Brongersma , S. W. Fong , H.-S. P. Wong , M.-J. Sher , P. Jost , M. Kaes , M. Salinga , A. von Hoegen , M. Wuttig , and A. M. Lindenberg , Phys. Rev. Lett. 117 , 067601 ( 2016 ). PRLTAO 0031-9007 10.1103/PhysRevLett.117.067601 [19] 19 R. Pandian , B. J. Kooi , G. Palasantzas , J. T. M. De Hosson , and A. Pauza , Adv. Mater. 19 , 4431 ( 2007 ). ADVMEW 0935-9648 10.1002/adma.200700904 [20] 20 E. Janod , J. Tranchant , B. Corraze , M. Querré , P. Stoliar , M. Rozenberg , T. Cren , D. Roditchev , V. Ta Phuoc , M.-P. Besland , and L. Cario , Adv. Funct. Mater. 25 , 6287 ( 2015 ). AFMDC6 1616-301X 10.1002/adfm.201500823 [21] 21 D. Loke , T. H. Lee , W. J. Wang , L. P. Shi , R. Zhao , Y. C. Yeo , T. C. Chong , and S. R. Elliott , Science 336 , 1566 ( 2012 ). SCIEAS 0036-8075 10.1126/science.1221561 [22] 22 F. Rao , K. Ding , Y. Zhou , Y. Zheng , M. Xia , S. Lv , Z. Song , S. Feng , I. Ronneberger , R. Mazzarello , W. Zhang , and E. Ma , Science 358 , 1423 ( 2017 ). SCIEAS 0036-8075 10.1126/science.aao3212 [23] 23 H. Hirori , K. Shinokita , M. Shirai , S. Tani , Y. Kadoya , and K. Tanaka , Nat. Commun. 2 , 594 ( 2011 ). NCAOBW 2041-1723 10.1038/ncomms1598 [24] 24 M. Liu , H. Y. Hwang , H. Tao , A. C. Strikwerda , K. Fan , G. R. Keiser , A. J. Sternbach , K. G. West , S. Kittiwatanakul , J. Lu , S. A. Wolf , F. G. Omenetto , X. Zhang , K. A. Nelson , and R. D. Averitt , Nature (London) 487 , 345 ( 2012 ). NATUAS 0028-0836 10.1038/nature11231 [25] 25 A. Pirovano , A. L. Lacaita , A. Benvenuti , F. Pellizzer , and R. Bez , IEEE Trans. Electron Devices 51 , 452 ( 2004 ). IETDAI 0018-9383 10.1109/TED.2003.823243 [26] 26 C. A. Werley , K. Fan , A. C. Strikwerda , S. M. Teo , X. Zhang , R. D. Averitt , and K. A. Nelson , Opt. Express 20 , 8551 ( 2012 ). OPEXFF 1094-4087 10.1364/OE.20.008551 [27] 27 H. Hirori , A. Doi , F. Blanchard , and K. Tanaka , Appl. Phys. Lett. 98 , 091106 ( 2011 ). APPLAB 0003-6951 10.1063/1.3560062 [28] 28 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.121.165702 for the properties of THz pulse (Sec. I), which includes Refs. [29,30]. [29] 29 L. V. Keldysh , J. Exp. Theor. Phys. 47 , 1945 ( 1964 ) JTPHES 1063-7761 29 [ L. V. Keldysh Sov. Phys. JETP 20 , 1307 ( 1965 )]. JTPHES 1063-7761 [30] 30 C. Lange , T. Maag , M. Hohenleutner , S. Baierl , O. Schubert , E. R. J. Edwards , D. Bougeard , G. Woltersdorf , and R. Huber , Phys. Rev. Lett. 113 , 227401 ( 2014 ). PRLTAO 0031-9007 10.1103/PhysRevLett.113.227401 [31] 31 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.121.165702 for the details about the initial crystallization (Sec. II). [32] 32 P. Němec , V. Nazabal , A. Moreac , J. Gutwirth , L. Beneš , and M. Frumar , Mater. Chem. Phys. 136 , 935 ( 2012 ). MCHPDR 0254-0584 10.1016/j.matchemphys.2012.08.024 [33] 33 V. Bragaglia , K. Holldack , J. E. Boschker , F. Arciprete , E. Zallo , T. Flissikowski , and R. Calarco , Sci. Rep. 6 , 28560 ( 2016 ). SRCEC3 2045-2322 10.1038/srep28560 [34] 34 A. T. Tarekegne , H. Hirori , K. Tanaka , K. Iwaszczuk , and P. U. Jepsen , New J. Phys. 19 , 123018 ( 2017 ). NJOPFM 1367-2630 10.1088/1367-2630/aa936b [35] The crystallization temperature is about Δ T L = 130 K for the initial temperature of 300 K [36] . For the case of E o = 175 kV / cm and an electric field enhancement of 18 times as indicated in Fig. 1(c) , a Δ T L of 130 K was obtained from Eq. (1) by using the conductivity of the amorphous phase. The electric conductivity in the amorphous region is usually described using the Poole-Frenkel (PF) model [19] , leading to an electrical conductivity σ PF ( E ) = σ 0 exp [ ( e 3 / 2 E 1 / 2 ) / ( k T π 1 / 2 ϵ 1 / 2 ) ] . Here, for the amorphous phase, σ 0 = 10 S / cm , ϵ = 30 ϵ 0 ( ϵ 0 is the permittivity of vacuum) [37] , and T = 300 K . [36] 36 S. Senkader and C. D. Wright , J. Appl. Phys. 95 , 504 ( 2004 ). JAPIAU 0021-8979 10.1063/1.1633984 [37] 37 F. Kadlec , C. Kadlec , and P. Kužel , Solid State Commun. 152 , 852 ( 2012 ). SSCOA4 0038-1098 10.1016/j.ssc.2012.02.018 [38] 38 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.121.165702 for the Δ r max / r in amorphous regions (Sec. IV). [39] 39 G. Weiser , U. Dersch , and P. Thomas , Philos. Mag. B 57 , 721 ( 1988 ). PMABDJ 1364-2812 10.1080/13642818808208489 [40] 40 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.121.165702 for the consideration of the thermal diffusion during THz pulse excitation (Sec. III), which includes Refs. [4,41–43]. [41] 41 J. P. Reifenberg , M. A. Panzer , S. Kim , A. M. Gibby , Y. Zhang , S. Wong , H.-S. P. Wong , E. Pop , and K. E. Goodson , Appl. Phys. Lett. 91 , 111904 ( 2007 ). APPLAB 0003-6951 10.1063/1.2784169 [42] 42 R. P. Sear , J. Phys. Condens. Matter 19 , 033101 ( 2007 ). JCOMEL 0953-8984 10.1088/0953-8984/19/3/033101 [43] 43 V. G. Karpov , Y. A. Kryukov , S. D. Savransky , and I. V. Karpov , Appl. Phys. Lett. 90 , 123504 ( 2007 ). APPLAB 0003-6951 10.1063/1.2715024 [44] 44 M. J. Shu , P. Zalden , F. Chen , B. Weems , I. Chatzakis , F. Xiong , R. Jeyasingh , M. C. Hoffmann , E. Pop , H.-S. P. Wong , M. Wuttig , and A. M. Lindenberg , Appl. Phys. Lett. 104 , 251907 ( 2014 ). APPLAB 0003-6951 10.1063/1.4884816 [45] 45 J.-L. Battaglia , A. Kusiak , V. Schick , A. Cappella , C. Wiemer , M. Longo , and E. Varesi , J. Appl. Phys. 107 , 044314 ( 2010 ). JAPIAU 0021-8979 10.1063/1.3284084 [46] 46 C. Zener , Proc. R. Soc. A 145 , 523 ( 1934 ). PRLAAZ 1364-5021 10.1098/rspa.1934.0116 [47] 47 E. O. Kane , J. Phys. Chem. Solids 12 , 181 ( 1960 ). JPCSAW 0022-3697 10.1016/0022-3697(60)90035-4 [48] 48 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.121.165702 for the details of Eq. (2) (Sec. IV), which includes Refs. [37,46,47,49–56]. [49] 49 L. Thamizhmani , A. K. Azad , J. Dai , and W. Zhang , Appl. Phys. Lett. 86 , 131111 ( 2005 ). APPLAB 0003-6951 10.1063/1.1896451 [50] 50 S. P. Mickan , K.-S. Lee , T.-M. Lu , E. Barnat , J. Munch , D. Abbott , and X.-C. Zhang , Proc. SPIE Int. Soc. Opt. Eng. 4591 , 197 ( 2001 ). PSISDG 0277-786X 10.1117/12.449149 [51] 51 Y. Zhao and D. R. Grischkowsky , IEEE Trans. Microwave Theory Tech. 55 , 656 ( 2007 ). IETMAB 0018-9480 10.1109/TMTT.2007.892798 [52] 52 N. Laman and D. Grischkowsky , Appl. Phys. Lett. 93 , 051105 ( 2008 ). APPLAB 0003-6951 10.1063/1.2968308 [53] 53 T. Siegrist , P. Jost , H. Volker , M. Woda , P. Merkelbach , C. Schlockermann , and M. Wuttig , Nat. Mater. 10 , 202 ( 2011 ). NMAACR 1476-1122 10.1038/nmat2934 [54] 54 A. Mendoza-Galván and J. González-Hernández , J. Appl. Phys. 87 , 760 ( 2000 ). JAPIAU 0021-8979 10.1063/1.371937 [55] 55 K. Shportko , S. Kremers , M. Woda , D. Lencer , J. Robertson , and M. Wuttig , Nat. Mater. 7 , 653 ( 2008 ). NMAACR 1476-1122 10.1038/nmat2226 [56] 56 J. Kim , J. Kim , and S.-H. Jhi , Phys. Rev. B 82 , 201312 ( 2010 ). PRBMDO 1098-0121 10.1103/PhysRevB.82.201312 [57] By using the optical fluence F , the reflectivity R , and the absorption coefficient α for the photon energy ℏ ω , the number of electrons N opt that are injected by the light can be expressed as N opt = F α ( 1 − R ) / ( ℏ ω ) [58] . For R and α we used the values 0.33 and 10.9 × 10 5 cm − 1 , respectively [55] . [58] 58 M. J. Shu , I. Chatzakis , Y. Kuo , P. Zalden , and A. M. Lindenberg , Appl. Phys. Lett. 102 , 201903 ( 2013 ). APPLAB 0003-6951 10.1063/1.4807731 [59] 59 W. A. Johnson and R. F. Mehl , Trans. AIME 135 , 416 ( 1939 ). [60] 60 M. Avrami , J. Chem. Phys. 7 , 1103 ( 1939 ). JCPSA6 0021-9606 10.1063/1.1750380 [61] 61 J. N. Hay , Br. Polym. J. 3 , 74 ( 1971 ). BPOJAB 0007-1641 10.1002/pi.4980030205 [62] 62 V. Weidenhof , I. Friedrich , S. Ziegler , and M. Wuttig , J. Appl. Phys. 89 , 3168 ( 2001 ). JAPIAU 0021-8979 10.1063/1.1351868
PY - 2018/10/19
Y1 - 2018/10/19
N2 - We have systematically investigated the spatial and temporal dynamics of crystallization that occur in the phase-change material Ge2Sb2Te5 upon irradiation with an intense terahertz (THz) pulse. THz-pump-optical-probe spectroscopy revealed that Zener tunneling induces a nonlinear increase in the conductivity of the crystalline phase. This fact causes the large enhancement of electric field associated with the THz pulses only at the edge of the crystallized area. The electric field concentrating in this area causes a temperature increase via Joule heating, which in turn leads to nanometer-scale crystal growth parallel to the field and the formation of filamentary conductive domains across the sample.
AB - We have systematically investigated the spatial and temporal dynamics of crystallization that occur in the phase-change material Ge2Sb2Te5 upon irradiation with an intense terahertz (THz) pulse. THz-pump-optical-probe spectroscopy revealed that Zener tunneling induces a nonlinear increase in the conductivity of the crystalline phase. This fact causes the large enhancement of electric field associated with the THz pulses only at the edge of the crystallized area. The electric field concentrating in this area causes a temperature increase via Joule heating, which in turn leads to nanometer-scale crystal growth parallel to the field and the formation of filamentary conductive domains across the sample.
UR - http://www.scopus.com/inward/record.url?scp=85055208798&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85055208798&partnerID=8YFLogxK
U2 - 10.1103/PhysRevLett.121.165702
DO - 10.1103/PhysRevLett.121.165702
M3 - Article
C2 - 30387634
AN - SCOPUS:85055208798
VL - 121
JO - Physical Review Letters
JF - Physical Review Letters
SN - 0031-9007
IS - 16
M1 - 165702
ER -