C. R. Physique ••• (••••)

JÈûM&âtm ELSEVIER

Contents lists available at ScienceDirect

Comptes Rendus Physique

www.sciencedirect.com

Thermoelectric mesoscopic phenomena / Phénomènes thermoélectriques mésoscopiques

Reduction of phonon mean free path: From low-temperature physics to room temperature applications in thermoelectricity

Reduction du libre parcours moyen des phonons : de la physique des basses températures aux applications en thermoélectricité à l'ambiante

Olivier Bourgeoisa b *, Dimitri Tainoffa b, Adib Tavakolia,b, Yanqing Liua b, Christophe Blanca b, Mustapha Boukharic, André Barskic, Emmanuel Hadjic

a Institut Néel, CNRS, 25, avenue des Martyrs, 38042 Grenoble, France b Université Grenoble Alpes, Inst Néel, 38042 Grenoble, France

c Institut Nanosciences et Cryogénie, SP2M, CEA-UJF, 17, rue des Martyrs, 38054 Grenoble, France

A R T I C L E I N F 0 A B S T R A C T

It has been proposed for a long time now that the reduction of the thermal conductivity by reducing the phonon mean free path is one of the best way to improve the current performance of thermoelectrics. By measuring the thermal conductance and thermal conductivity of nanowires and thin films, we show different ways of increasing the phonon scattering from low-temperature up to room-temperature experiments. It is shown that playing with the geometry (constriction, periodic structures, nano-inclusions), from the ballistic to the diffusive limit, the phonon thermal transport can be severely altered in single crystalline semiconducting structures; the phonon mean free path is in consequence reduced. The diverse implications on thermoelectric properties will be eventually discussed.

© 2016 Académie des sciences. Published by Elsevier Masson SAS. This is an open access

article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

R É S U M É

Il a été proposé depuis longtemps maintenant que diminuer la conductivité thermique en réduisant le libre parcours moyen des phonons est une des meilleures façons d'améliorer les performances actuelles des matériaux thermoélectriques. En mesurant la conductance thermique et la conductivité thermique de nanofils et de couches minces, nous montrons différentes manières d'augmenter la diffusion des phonons des basses températures jusqu'à la température ambiante. Il est montré qu'en jouant sur la géométrie (constriction, structures périodiques, nano-inclusions), à partir de la limite balistique jusqu'à la limite diffusive, le transport thermique phononique peut être altéré de façon significative dans des structures monocristallines semiconductrices, le libre parcours moyen

Article history: Available online xxxx

Keywords: Phonon

Thermal transport

Thermoelectricity

Semiconductor

Nanowire

Inclusions

Mots-clés : Phonon

Transport thermique

Thermoélectricité

Semiconducteur

Nanofil

Inclusions

* Corresponding author.

E-mail address: olivier.bourgeois@neel.cnrs.fr (O. Bourgeois).

http://dx.doi.org/10.1016/jxrhy.2016.08.008

1631-0705/© 2016 Académie des sciences. Published by Elsevier Masson SAS. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

2 O. Bourgeois etal./C. R. Physique ••• (••••) •••-•••

des phonons étant réduit en conséquence. Enfin, les diverses implications sur les propriétés thermoélectriques seront discutées. © 2016 Académie des sciences. Published by Elsevier Masson SAS. This is an open access

article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

1. Introduction

In the two last decades, materials structured at the nanoscale to have low thermal conductivity have attracted a lot of attention, due to their high-potential thermoelectric applications [1-4]. Indeed, as proposed by Dresselhaus et al. in 1993, nanostructuration can significantly improve the thermoelectric efficiency [5,6]. This efficiency of thermoelectric materials is given by the dimensionless figure of merit ZT = S2Ta/k, where S is the Seebeck coefficient, a the electrical conductivity, and k the thermal conductivity. If we focus only on the thermal properties, ZT will be increased by reducing the phonon thermal transport. This reflects the fact that the heat source and the cold source must be as thermally isolated as possible to prevent too much heat leaking, as that will degrade the overall performances of the thermoelectric modules.

For electrons, the thermal and the electronic properties are coupled through the Wiedemann-Franz law. Thus, the only way of reducing the thermal transport is to reduce the phonon thermal conductivity. Two major directions are available to act on the phonon properties: incoherent and coherent effects. On the one hand, incoherent effects are ruled by the phonon mean free path (MFP), then the reduction is obtained by adding scattering centers, by reducing the size, having rough surface or by adding nano-inclusions. On the other hand, coherent effects will manifest themselves through the reduction of the group velocity (band flattening), opening of band gaps or having destructive interferences. Two significant length scales govern phonon physics, involving incoherent or coherent effects: the phonon mean free path (A) and the dominant phonon wavelength given by:

^dom — '

2.82kB T (1)

where h the Planck constant, kB the Boltzmann constant and vs is the speed of sound. This length is indeed temperature dependent, increasing as the temperature decreases (see Fig. 1).

Here, we will focus our attention on the experimental attempts to reduce the phonon mean free path in single-crystal semiconducting materials by playing with the geometry of its surface or with its internal structure. Other means based on coherent effects will be discussed at the end of the paper. There are actually two complementary approaches to realizing an efficient phonon engineering at the nanoscale: the top-down one based on e-beam lithography to make small structures like nanowire and membrane, and the second one, the bottom-up, generally based on specific growth of nanostuctured materials (bulk or nanowire) using for instance molecular beam epitaxy (MBE).

These two important limits will be experimentally addressed in this paper summarizing the recent findings of our group: the low-temperature one set by boundary scattering limit where both important lengths (A and Xdom) play a role and the room temperature limit, where for the time being only the mean free path can be engineered efficiently.

2. Towards the ballistic limit: thermal transport in Si nanowires

In low-dimensional systems, the surface scattering of phonons plays a crucial role in the understanding of thermal transport [7-9]. The so-called Casimir regime occurs when the bulk phonon MFP is much bigger than the actual size of the

Phonon spectral density _300K

200K 100K

Fig. 1. Phonon spectral density as a function of the wavelength for different equilibrium temperatures of a black body. The dominant phonon wavelength for each temperature corresponds to the maximum of the curves. Typically, the dominant phonon wavelength is of the order of 1 nm at 300 K and 100 nm around 1 K.

O. Bourgeois et al. / C. R. Physique ••• (••••) •

relevant length scales transport regime

models

relevant formulae

, Diffusive regime

Xdom{T)<Aph (T)< L (incohereënt)

Kinetic model

K=-CvA

*d„m (T) < L ~ Aph < Aph-bulk Casimir regime (incoherent) Casimir model = ß^CasT'

*dom (T) ~ L < Aph Ziman regime (partly coherent) Ziman model A ph ~ , ACas I- P

L < *dom (T) < Aph Ballistic regime (coherent) Landauer model T= 1 l+3L/4Aeff

Fig. 2. Classification of the different regimes of phonon transport in nanoscaled systems at low enough temperature, typically below 30 K for silicon. L is the sample length, Aph is the inelastic phonon mean free path, and !jom the dominant phonon wavelength. Concerning the Landauer limit, the parameter T denotes the transmission coefficient between a nanowire and its heat reservoir. Other important lengths are explained in the text.

nanosystem under study. In that case, the phonon MFP is only set by the scattering on the surfaces. This limit is reached at low enough temperature when other scattering mechanisms do not dominate (phonon-phonon interaction, electron-phonon interaction, scattering on impurities etc.). As a result, effects linked to surface roughness have a significant impact on phonon transport [10,11]. These effects can be simple incoherent scattering when the phonon wavelength is much bigger than the mean asperity or, as the temperature is lowered, may become wavelength dependent. Indeed, if the wavelength is bigger than the mean roughness, specular reflections of phonons may occur on the surface. If elastic scattering is involved, the energy is conserved, leading to an increase in the MFP.

2.1. Normal phonon transport in Si nanowires at low temperature

In order to understand the thermal transport at low temperature, a summary of the Casimir and Ziman models is necessary. In the Casimir model, [7,9], the thermal conductance is given by:

K Cas = 3.2 X 103

2 k4 kB

5h3 v 3

e x w Aca: L

T3 = ßca

where kB is the Boltzmann constant, h the Planck constant and Vs refers to the sound velocity. Acas = 1. 12e x w is the Casimir MFP of the phonons in a nanowire having a rectangular section, e refers to the thickness and w to the average width of the nanowire; L is its length. This formula is obtained when the mean free path of phonons is only limited by the cross-section of the wire, i.e. by boundary scattering. Indeed, in this limit, each phonon hitting the rough surface is assumed to be absorbed an re-emitted in all direction and energy following the black body radiation law (see Fig. 1). In other words, in this Casimir limit, the surface is considered as infinitely rough whatever the phonon wavelength.

As the temperature is lowered, the dominant phonon wavelength actually grows (see Eq. (1)) and becomes eventually larger than the roughness amplitude; then specular reflections become more probable and cannot be neglected like in the Casimir model. To take this trend into account, the Casimir mean free path should be replaced in Eq. (2) by an effective MFP Aeff, as proposed by Ziman and colleagues [8,9]:

Aeff = 1-P ACa

where p is the parameter describing the probability for a phonon to be specularly reflected at the surfaces. The value of this parameter will depend on both the wavelength of the phonons and the roughness of the wire surface. The Casimir purely diffusive case is given by p = 0 and the purely specular one, when p = 1, will be referred to as the ballistic limit; the MFP is diverging to infinite. The Ziman model takes into account a value for p that is different from zero; as the temperature is lowered, the phonon mean free path can be bigger than the section of the nanowire. The different phonon transport regimes and their associated models are summarized in the Table given in Fig. 2.

As a conclusion, the thermal transport in silicon nanowires having a section of 100 nm will be well described by the Casimir model, and then at lower temperature by the Ziman model, as it will be seen in the next section. Indeed, the length L is smaller than the bulk phonon mean free path Aph-bulk in Si, and on the same order of magnitude as the dominant phonon wavelength. At even lower temperature, the fully ballistic limit should be described by a Landauer-type model where the thermal transport is given by transmission coefficients between the nanowire and the heat bath. This goes clearly beyond our focus in this paper; there is currently no experimental data in this regime.

4 O. Bourgeois et al. / C. R. Physique ••• (••••)

Fig. 3. A SEM images of: a) a suspended nanowire linked to the two pads where the electrical connections permit the 3m measurements; b) a straight nanowire c) corresponds to the side view of a S-shaped nanowire.

3.0 -T-,-,-,-T-,-,-1-T-2

0.0 I—-1-.-1-1-1-.-1-.-

0 1 2 3 4 5

Fig. 4. Thermal conductance of 5 |im long nanowires normalized to four times the universal value of thermal conductance versus temperature. The thermal conductance of the straight nanowire is much bigger than the thermal conductance of the S-shape nanowires. The presence of the serpentine section stops the transmission of ballistic phonons, reducing the overall phonon mean free path.

2.2. Phonon blocking with S-shaped chicanes in Si nanowires

We will illustrate in the following section these different scatterings and phonon transport regimes in three situations: firstly the thermal transport in a normal straight silicon nanowire, secondly in a Si nanowire having S-shaped chicanes, and thirdly in Si nanowires with engineered surfaces. The nanowire thermal conductance has been measured using longitudinal 3« technique with a highly sensitive thermometry based on niobium nitride (NbN) [12,10,11,13]. The Si nanowires are fabricated from silicon on insulator (SOI) substrates by e-beam lithography. The total length of each structure has been purposely set to 10 jam in order to easily compare the thermal measurement of all kinds of nanowires (see Fig. 3). The section of the nanowires is 100 nm by 200 nm. The central symmetry of the structures is necessary for the 3« method employed to measure the thermal conductance [12].

The temperature profile looks like a parabola with the hottest point located at the center of the nanowire. The heat generated at each points of the suspended structure flows from the center to the heat bath on both sides. The thermal conductance is generally measured using a power of a few tens of femtowatts dissipated in the niobium nitride transducer deposited on top of the nanowires. This power creates on average a temperature gradient smaller than 1 mK. The values of the measured thermal conductance are around a few universal value of thermal conductance K0 = n2k2T/3h ~ 10-12T W/K2; a value used to normalized the data.

Equation (2) gives the correct order of magnitude for the thermal conductance, since a K value close to 10-11 W/K is predicted at 3 K and a value of 1.5 x 10-11 W/K is measured on a straight nanowire. For the explanation of the temperature variation of the thermal conductance, it is needed to take into account the contribution of specular reflection by introducing the Ziman mean free path, and hence the parameter p. This has been described in details in the past [10,11,14]. We would like to highlight here that the contribution of ballistic phonons to the thermal transport is significant. We have proposed a solution for thermal management that will not perturb diffusive transport but will affect only the ballistic one. Specific suspended nanowires have been made, including a double-bend chicane structures 400 nm in length; the total extended length of the S-shaped nanowire itself being identical to the one of the straight nanowire (see Fig. 3c and Fig. 4).

In recent works [14,15], we have shown that the S-shaped nanowires have a much lower thermal conductance than the straight one, as it can be seen in Fig. 4. If transport were only diffusive, this constriction would have no effect at all on thermal conductance. Actually, the severe reduction observed (reduction by more than 30 %) demonstrates two significant facts: first, this illustrates the contribution of ballistic phonons to thermal transport and second, including a chicane in a thermal conductor is a very efficient way of reducing thermal transport. Phonons having long mean free path (bigger than the Casimir mean free path Acas) will be blocked and reflected by the chicane acting like a mirror for the long-wavelength acoustic plane waves. Conversely, short-wavelength phonons that are more scattered by the nanowire surfaces will not be much affected by the chicane, their mean free path being already set by the section of the nanowire like in the Casimir

O. Bourgeois etal./C. R. Physique ••• (••••) •••-••• 5

■0.6 -0.4 -0.2 0.0 0.2 0.4 0.6

Fig. 5. Normalized thermal conductance versus temperature for a straight nanowire (Straight) and four corrugated nanowires (C1 to C4) in a log-log plot. The normalization coefficient /¡Cas corresponds to the Casimir limit value of the thermal conductance at 1 K. The corrugated nanowire Q can have much lower thermal conductance that smooth nanowire. The phonon MFP is reduced severely, even below the Casimir limit. As an illustration, in insert we compare the average phonon mean free path in the straight nanowire (in yellow), in the Casimir limit (in red), and as measured in corrugated nanowire (in blue). The corrugations reduce significantly the phonon mean free path.

regime. This concept has been successfully used to thermally isolate a suspended platform dedicated to sensitive thermal measurements [16].

2.3. Reduced thermal transport in corrugated nanowires

In the previous section, we have seen how a macroscopic constriction can act as a phonon filter. In this section, we will illustrate how phonon backscattering or multiple scatterings can be introduced by using nano-engineered surfaces. In this respect, phononic crystals, defined as a regular modulation of the geometry of a phonon waveguide, have attracted significant attention in the past decade as a promising means to reduce phonon transport. In the present samples, periodic nanostructurations are used (see Fig. 5) to induce phonon backscattering and reduce incoherently the thermal conductance, not coherently like in a real thermocrystal. We have reported that in that kind of engineered corrugated nanowires, phonon transport can be decreased even below the Casimir limit [15,17].

In order to compare thermal transport between various geometries of nanowires, the thermal conductances have to be normalized to a geometrical factor called the Casimir coefficient ^Cas, according to Eqs. (2) and (3). Fig. 5 reports the thermal conductance plotted as K/^Cas.

This normalization yields the temperature power law and allows a direct comparison between the thermal transport of each nanowire. Indeed, from Eqs. (2) and (3), the normalized conductance can be expressed as:

log(K/fas) = a + b log(T) (4)

(K/fas) = Tb (5)

From these two equations, we can extract the power-law factor b of the temperature variation of the thermal conductance, the value of the MFP Aeff and the parameter p [15,17].

The thermal conductance of the corrugated nanowires is clearly below the thermal conductance of straight nanowires, even if the power law has not been changed; it is very close to 3 (b ^ 2.6). The presence of artificial periodic nanostruc-turation increases phonon scattering (backscattering, phonon trapping etc.). Surface roughness is not changed, but the mean free path is clearly reduced. This can be depicted by ray-tracing Monte Carlo numerical simulations. Indeed, the multiple scattering of ballistic phonons can be seen through a corpuscular picture explaining the fact that the phonons are somehow trapped in the core of the nanowire. As a consequence, the phonon MFP that can be even smaller than the Casimir limit, i.e. smaller than the mean diameter of the nanowire. This is characterized by the necessity of putting a negative p parameter to fit the data. Even if this special result does not involve phonon band engineering like in a real phononic crystal, it gives however a fantastic way of turning a single-crystal semiconductor into a phonon glass.

3. In the diffusive limit: nano-inclusions limit thermal conductivity in semiconductors

After having shown in details the different ways of affecting the thermal transport in the presence of ballistic phonons, we will focus now on room-temperature effects where the diffusive regime is governing the thermal conductivity. We will show that engineering the phonon thermal properties is also possible with a semiconductor at room temperature (RT). As

O. Bourgeois et al. / C. R. Physique ••• (••••) •

Fig. 6. TEM image of the germanium matrix with GeMn spherical nano-inclusions at 8 % of Mn. The presence of numerous inclusions will significantly increase the phonon scattering without perturbing the electronic transport. The scale bar represents 50 nm. The distance between the nano-inclusions is indeed smaller than the average phonon mean free path.

compared to the ballistic limits, working at RT in the diffusive limit will have significant implications for real applications in managing heating, cooling, or harvesting energy.

The great challenges mostly rely on our capacity to improve the performances of semiconductor materials. As, on the one hand, a semiconductor is by far too thermally conductive (around 100 W/mK), on the other hand its Seebeck coefficient is generally quite high (above 200 uV/K). Then, obtaining an efficient thermoelectric material having a high ZT based on a semiconductor is possible by decreasing significantly its thermal conductivity.

The increase in the thermoelectric efficiency is feasible thanks to the difference of mean free path existing between phonons (~100 nm at 300 K) and electrons (~1 nm at 300 K). By introducing in a given material a structural disorder at the nanometer scale, it is possible to induce phonon diffusion without affecting charge carriers, justifying the innovative concept of electron crystal-phonon glass material.

Indeed, in a semiconducting matrix, phonons contribute the most to heat transport, having an average mean free path of ~100 nm at room temperature and a dominant wavelength of 1 nm at RT. With the objective of manipulating and controlling phonon transport in a SC, the engineering of materials can then follow the two different approaches mentioned in the introduction: either by playing with phonon wavelength [18-21] or with phonon mean free path [17,20,21]. At room temperature, the most efficient effect will be the reduction in the phonon mean free path using a geometrical scattering approach. Indeed, as it has been recently experimentally demonstrated, the phonon wavelength is far two small at RT for the nanostructuration to create sizable effects [20,21].

The best large-scale applicability is expected from bulk 3D SC material that has been engineered at the low dimension for thermal management. One of the most important concepts is the electron crystal-phonon glass, being an inclusion-embedded semiconducting system, which shows great advantages and potentials [22]. For crystalline semiconductors, the ideal model would be a defect-free highly crystalline matrix containing nanoscaled inclusions with a wide diameter distribution. It can be theoretically predicted that as the inter-distance of the inclusions is comparable to the phonon mean free path in the matrix, the nano-inclusions can prohibit efficiently phonon transport. Indeed, by introducing more phonon scattering processes, we expect to get a germanium based electron crystal-phonon glass material.

The materials of interest here are thin films of a single-crystalline germanium matrix embedded with randomly distributed Ge3Mn5 nano-inclusions grown by epitaxy on a germanium substrate (see Fig. 6). The nano-inclusions are found to be nearly spherical with a diameter distribution from 5 to 50 nm, depending on the growth parameters. Furthermore, manganese atoms act as p-type dopant in the Ge matrix, which ensures a good electrical conductivity of the material.

The experimental measurements of the thermal conductivities of the Ge:Mn materials have been done using a highly sensitive differential 3w technique [23]. Experimental results of 3w measurements revealed strong reductions in the thermal conductivity, compared to the Ge bulk at room temperature, for Ge:Mn thin films containing different Mn concentrations.

The major experimental finding is a reduction by more than a factor of ten with respect to the bulk value when there are 8 % of Mn in the Ge matrix. A thermal conductivity below 8 W/mK has been measured, which has to be compared to the bulk thermal conductivity, i.e. 65 W/mK. The reduction is fully ascribed to the presence of evenly distributed phonon scatterers that significantly reduce the phonon mean free path.

4. Conclusions

Different concepts of phonon thermal transport reduction have been shown experimentally, from the ballistic limit to the diffusive limit. These new paradigms will be fruitfully applied for the development of new materials or devices in thermoelectricity and energy harvesting in general. It seems that now, we are facing a kind of intrinsic limit in that reduction by non-coherent effects, a lower limit close to the value of the alloy's thermal conductivity. Unfortunately, engineering the

O. Bourgeois etal./C. R. Physique ••• (••••) •••-••• 7

phonon band by opening gaps like in a phononic crystal seems inappropriate for room-temperature thermal management, because phonon wavelengths are far too small as compared to the length accessible by current fabrication methods.

Other technological solutions may be found in playing also on the electrical properties of semiconductor-based systems to improve the Seebeck coefficient. Indeed, as proposed by Hicks and Dresselhaus [5], it could be possible to increase the Seebeck coefficient by playing with the electronic density of state at the Fermi level. Then, engineering the electronic properties by various means (nanostructuration, adding quantum dots, etc.) seems absolutely necessary if one wants to boost the performances of future thermoelectric materials.

Acknowledgements

We acknowledge technical supports from Nanofab, the Cryogenic and the Electronic facilities and the "Pôle Capteur thermométrique et Calorimétrie" of the Institut Néel for these experiments. Funding for this project was provided by a grant from the "Région Rhône-Alpes" (Cible), by the "Agence nationale de la recherche" (ANR) through the project QNM No. 0404 01 and Mesophon grant No. 15-CE30-0019 by the European Projects: MicroKelvin FP7 low-temperature infrastructure, Grant No. 228464, and MERGING Grant No. 309150.

References

[1] D.G. Cahill, W.K. Ford, K.E. Goodson, G.D. Mahan, A. Majumdar, H.J. Maris, R. Merlin, S.R. Phillpot, J. Appl. Phys. 93 (2003) 793.

[2] A. Shakouri, Annu. Rev. Mater. Res. 41 (2011) 399.

[3] D.G. Cahill, P.V. Braun, G. Chen, D.R. Clarke, S.H. Fan, K.E. Goodson, P. Keblinski, W.P. King, G.D. Mahan, A. Majumdar, H.J. Maris, R. Merlin, S.R. Phillpot, E. Pop, L. Shi, J. Appl. Phys. 1 (2014) 011305.

[4] S. LeBlanc, S.K. Yee, M.L. Scullin, C. Dames, K.E. Goodson, Renew. Sustain. Energy Rev. 32 (2014) 313.

[5] L.D. Hicks, M.S. Dresselhaus, Phys. Rev. B 47 (1993) 16631.

[6] M.S. Dresselhaus, G. Chen, M.Y. Tang, R.G. Yang, H. Lee, D.Z. Wang, R.F. Ren, J.P. Fleurial, P. Gogna, Adv. Mater. 19 (2007) 1.

[7] H.B.G. Casimir, Physica 5 (6) (1938).

[8] J.M. Ziman, Electrons and Phonons, Clarendon Press, Oxford, UK, 2001.

[9] R. Berman, E.L. Foster, J.M. Ziman, Proc. R. Soc. Lond. Ser. A 231 (1955) 130.

[10] J.S. Heron, T. Fournier, N. Mingo, O. Bourgeois, Nano Lett. 9 (2009) 1861.

[11] J.S. Heron, T. Fournier, N. Mingo, O. Bourgeois, Nano Lett. 10 (2010) 2288.

[12] O. Bourgeois, T. Fournier, J. Chaussy, J. Appl. Phys. 101 (2007) 016104.

[13] O. Bourgeois, E. André, C. Macovei, J. Chaussy, Rev. Sci. Instrum. 77 (2006) 126108.

[14] J.S. Heron, C. Bera, T. Fournier, N. Mingo, O. Bourgeois, Phys. Rev. B 82 (2010) 155458.

[15] C. Blanc, J.-S. Heron, T. Fournier, O. Bourgeois, Appl. Phys. Lett. 105 (2014) 043106.

[16] S. Alaie, D.F. Goettler, K. Abbas, M.F. Su, C.M. Reinke, I. El-Kady, Z.C. Leseman, Rev. Sci. Instrum. 84 (2013) 105003.

[17] C. Blanc, A. Rajabpour, S. Volz, T. Fournier, O. Bourgeois, Appl. Phys. Lett. 103 (2013) 043109.

[18] N. Zen, T.A. Puurtinen, T.J. Isotalo, S. Chaudhuri, I.J. Maasilta, Nat. Commun. 5 (2014) 3435.

[19] M. Nomura, J. Nakagawa, Y. Kage, J. Maire, D. Moser, O. Paul, Appl. Phys. Lett. 106 (2015) 143102.

[20] J. Maire, R. Anufriev, H. Han, S. Volz, M. Nomura, arXiv:1508.04574.

[21] M.R. Wagner, B. Graczykowski, J.S. Reparaz, A. El Sachat, M. Sledzinska, F. Alzina, C.M. Sotomayor Torres, arXiv:1511.07398.

[22] W. Kim, J. Zide, A. Gossard, D. Klenov, S. Stemmer, A. Shakouri, A. Majumdar, Phys. Rev. Lett. 96 (2006) 045901.

[23] Y.Q. Liu, D. Tainoff, M. Boukhari, J. Richard, A. Barski, P. Bayle-Guillemaud, E. Hadji, O. Bourgeois, IOP Conf. Ser., Mater. Sci. Eng. 68 (2014) 012005.