- Nano Express
- Open Access
- Published:

# Tailoring of Seebeck coefficient with surface roughness effects in silicon sub-50-nm films

*Nanoscale Research Letters*
**volume 7**, Article number: 169 (2012)

## Abstract

The effect of surface roughness on the Seebeck coefficient in the sub-50-nm scale silicon ultra thin films is investigated theoretically using nonequilibrium Green's function formalism. For systematic studies, the surface roughness is modelled by varying thickness periodically with square wave profile characterized by two parameters: amplitude (*A*_{0}) and wavelength (*λ*). Since high Seebeck coefficient is obtained if the temperature difference between the ends of device produces higher currents and higher induced voltages, we investigate how the generated current and induced voltage is affected with increasing *A*_{0} and *λ*. The theoretical investigations show that pseudoperiodicity of the device structure gives rise to two effects: firstly the threshold energy at which the transmission of current starts is shifted towards higher energy sides and secondly transmission spectra of current possess pseudobands and pseudogaps. The width of the pseudobands and their occupancies determine the total generated current. It is found that current decreases with increasing *A*_{0} but shows a complicated trend with *λ*. The trends of threshold energy determine the trends of Seebeck voltage with roughness parameters. The increase in threshold energy makes the current flow in higher energy levels. Thus, the Seebeck voltage, i.e. voltage required to nullify this current, increases. Increase in Seebeck voltage results in increase in Seebeck coefficient. We find that threshold energy increases with increasing *A*_{0} and frequency (1/*λ*). Hence, Seebeck voltage and Seebeck coefficient increase vice versa. It is observed that Seebeck coefficient is tuneable with surface roughness parameters.

## 1. Introduction

It is reported that when the two ends of the thermoelectric material are kept at different temperatures, the voltage is generated [1]. This effect is called as the Seebeck effect. This effect has useful applications as it shows that heat can be converted to electricity and vice versa. In the last few years due to quest for finding, the highly efficient ways to produce and use energy there has been an increased interest in the thermoelectric performance of a device [2]. Moreover, thermoelectric energy conversion paves the way for utilization of waste heat generated in the devices into useful electric power. The extraction of waste heat is an important issue, especially nowadays, when the devices are reaching the sub-50-nm scale. Thus, intensive research has been reported for the effect of dimensionality on the thermoelectric phenomenon [3–5]. As the dimensions reduce down to sub-50-nm scale, the surface roughness plays a crucial role in determining the performance of the device. In this scale, the spatial extent of variations due to roughness is comparable to the device dimensions [6]. Hence, the investigation on the effect of surface roughness is important and crucial [7–9]. In this article, we have made a systematic study of the surface roughness effects on Seebeck coefficient. The surface roughness has been modelled by varying the thickness of the film periodically with a square wave profile as shown in Figure 1. The roughness is characterized by two parameters: amplitude (*A*_{0}) and wavelength (*λ*) of the square wave [10]. These parameters are analogous to the root mean square roughness and the roughness autocorrelation length which are commonly used to determine the surface roughness morphology [11]. We study the effect of roughness with increasing *A*_{0} and *λ* on the Seebeck effect. In other words, we investigate how does the induced voltage, *V*_{Seebeck}, for a fixed temperature difference between the two ends of a device, is affected due to increasing *A*_{0} and *λ* of the surface roughness of the device. To clearly bring the role of surface roughness the electron transport is considered to be coherent inside the device. The phase coherent transport is possible in the devices with the dimensions in the sub-50-nm scale [12].

When the two ends of the device are subjected to a temperature difference, initially, there is a net flow of current *I*_{net} from one end of the device to the other as shown in Figure 2a. This is so because, due to temperature difference, there is a different distribution of electrons in different energy states at the two ends. Hence, the number of electrons moving from the cold end to the hot end in an energy level '*ε*' is not the same as that of from hot end to cold end. This results in net flow of electrons in an energy level '*ε*'. The total current *I*_{net} is obtained from adding the contributions from all the energy levels. This net flow of current results in an accumulation of electrons at one end which builds an electric field opposing the incoming electrons as shown in Figure 2b. Further accumulation of charge stops when no more electrons are able to make it due to the building of electric field. The voltage developed due to the electric field is called Seebeck voltage *V*_{Seebeck}. Seebeck coefficient is obtained by dividing the Seebeck voltage by temperature difference [13]. As shown in Figure 2c, Seebeck voltage can be measured by calculating the amount of voltage required '*V*_{applied}' to nullify the current '*I*_{net}' due to temperature gradient. In this article, we focus on two important issues: firstly we investigate how does, for the same temperature difference between the two ends, *I*_{net} changes with *A*_{0} and *λ*. Secondly, it has been observed that for some roughness parameters, though the *I*_{net} decreases but the voltage applied to nullify it, *V*_{Seebeck}, increases. We investigate the reasons behind this behaviour.

## 2. Theory

The expression for the current has been calculated using nonequilibrium Green's function formalism (NEGF) [14]. For the devices in the sub-50-nm scale quantum effects become important and quantum transport model is required for description of transport phenomenon. In this regards, NEGF formalism, in which quantum effects are inherent, provides an efficient framework to model the electron transport in thermoelectric devices. A brief description is given below and detailed development can be found in [15].

The isolated channel and its energy levels are described by Schrodinger equation using Hamiltonian H, potential U and energy eigenvalues of the electron *ε*_{
α
}

But when the channel gets connected to the contacts it becomes an open system and the electrons can come and go out of the channel. To represent such an open system, two extra terms get introduced in the Schrodinger equation. The Schrodinger equation becomes

∑_{1} and ∑_{2}, called self-energy, represent the coupling between the channel and the two contacts. ∑'*s* are a complex quantity. Its imaginary part signifies that the electrons can decay out of the channel.

The real part of ∑'*s* just lifts up the energy eigenvalues of the channel. The second term that gets introduced in the Schrodinger equation is {*S*}, it is called the source term and is like an actual source driving the system.

Solution of the Schrodinger equation is given by

where

[*G*] is called the Green's function of the system. In the energy domain, Green function gives the energy eigenvalues for the eigenstates that are occupied in response to the applied impulse.

*A*_{1}(*E*) and *A*_{2}(*E*) represent density of states due to coupling of the channel with contacts 1 and 2, respectively. The probability of filling of these states is determined by the Fermi Dirac distribution function, *f*^{1}(*E*) and *f*^{2}(*E*), of the contacts 1 and 2, which are at temperature *T* and *T+* Δ*T* respectively.

Carrier concentration is given by:

The net inflow of current in channel at contact 1 is the same as the net outflow of current at contact 2. Net current flow in channel at contact 1 is the difference between the inflow and outflow of current at contact 1.

If *f*^{1}(*E*) = 1 and *f*^{2}(*E*) = 0 then *T*(*E*) represents the maximum current flow in the channel. For the electron transport in *x*-direction with confinement present in *z*-direction, the expression for I(E), after integrating over the energies along *y*-direction, becomes

where *n*_{
z
} is the subband index.

## 3. Results and discussion

The information about how does *I*_{net}, i.e. the current due to the temperature gradient, and, *V*_{Seebeck}, i.e. the voltage required to nullify it, depend on *A*_{0} and *λ*, is obtained by plotting the current voltage (*I*-*V*) characteristics for devices with different surface roughness parameters. Figure 3 gives the *I*-*V* characteristics for (a) different amplitudes: *A*_{0}'s = 0.1 nm, 0.3 nm and 0.5 nm and *λ* = 2.5 nm and (b) different wavelengths: *λ*'s = 2.5 nm, 3.3 nm and 5 nm and *A*_{0} = 0.5 nm. The plot for smooth surface is given in Figure 3a for comparison. The temperature gradient of Δ*T* = 20 K is applied across the device. The graphs are obtained for the devices in the sub-50-nm scale by considering the parameters: length of the channel *L*_{
c
} = 10 nm and average thickness *T*_{
b
} = 2 nm. As *I*_{net} is the current generated only due to temperature gradient with *V*_{applied} = 0, it is given by the intercept with the *y*-axis. It is observed from Figure 3a, b that *I*_{net} decreases with increasing *A*_{0} but with *λ* there is no definite trend. *V*_{Seebeck} is the voltage applied to make the current due to temperature gradient equal to zero and is obtained by the intercept with the *x*-axis. It is observed from Figure 3a, b that *V*_{Seebeck} increases with increasing *A*_{0} and 1/*λ*.

The above trends of *I*_{net} and *V*_{Seebeck} with the roughness *A*_{0} and 1/*λ* can be understood by examining the current flowing through each energy level *I*(*E*). As seen from Equation 7, *I*(*E*) is given as the product of transmission *T*(*E*) and $\left({f}_{1\mathrm{D}}^{\mathsf{\text{cold}}}\left(E\right)-{f}_{1\mathrm{D}}^{\mathsf{\text{hot}}}\left(E\right)\right)$, i.e. difference in the occupancy of cold ${f}_{1D}^{\mathsf{\text{cold}}}\left(E\right)$ and hot junctions ${f}_{1D}^{\mathsf{\text{hot}}}\left(E\right)$. The transmission *T*(*E*) gives the maximum current through the energy level '*E*' and depends on the density of states *A*(*E*) and velocity, i.e. rate of flow of electrons in and out of that energy level. As such the geometry of the device does not affect the occupancy and hence $\left({f}_{1\mathrm{D}}^{\mathsf{\text{cold}}}\left(E\right)-{f}_{1\mathrm{D}}^{\mathsf{\text{hot}}}\left(E\right)\right)$ but it affects the transmission *T*(*E*) of electrons from energy level *E* through density of states *A*(*E*). Thus, the trends of *I*_{net} and *V*_{Seebeck} with *A*_{0} and *λ* can be understood by examining how the transmission gets affected by roughness.

### 3.1. Effect of roughness on transmission spectra

Figure 4 gives the plots of *T*(*E*) for (a) roughness wavelength: *λ*'s = 2.5 nm, 3.3 nm and 5 nm with *A*_{0} = 0.5 nm and (b) roughness amplitude: *A*_{0}'s = 0.1 nm, 0.3 nm and 0.5 nm with *λ* = 2.5 nm. The following trends are observed: (i) threshold energy, i.e. the 'onset' energy at which the transmission starts, increases with increasing *A*_{0} and 1*/λ*. (ii) There are regions of psuedobands and pseudogaps. The width of the pseudobands decreases with increasing *A*_{0} and *λ*. Within the pseudobands, there are peaks in transmission at some energy. The phenomenon of increase in threshold energy has also been observed experimentally in quantum wells with thickness < 4 nm [16]. We would like to mention here that it has been observed only for thicknesses below 4 nm because surface roughness effects become prominent for ultra thin films. The criterion for the thickness below which surface roughness effects become prominent is discussed in [16]. So, our discussions can be experimentally validated only for ultra thin films. We will analyse the implications of the above trends on *I*_{net} and *V*_{Seebeck} but before that we discuss in brief the presence of the above trends.

The increase in threshold energy occurs due to the varying thickness seen by the electron while crossing the channel. The regions of the channel where the thickness is large, the confinement is small and the regions where the thickness is small, the confinement is high as seen from Figure 1. Since the lowest value of *ε*_{
z
} , which corresponds to the subband energy, depends on the confinement it alternates with the same periodicity as the thickness as shown in Figure 5. The value of *ε*_{
z
} , i.e. subband energy is large for high confinement and vice versa. The subband energy corresponds to the potential energy seen by the electron. The electron crossing the channel experiences this varying potential energy landscape which was not the case when the surface was smooth. As the electron crossing the channel moves from the region of low potential energy to the high potential energy it feels an additional confinement effect which was not present for smooth surface. This additional confinement effect causes the increase in lowest energy which the electron can take. Thus, the threshold energy, which corresponds to the lowest energy eigenvalue, increases. The confinement effect will be more for increasing *A*_{0} as the electron while moving from low potential energy to high potential energy region sees higher barrier heights. Thus, the lowest energy eigenvalue, i.e. threshold energy, increases with increasing *A*_{0}. The increase in threshold energy with increasing 1/*λ* can be explained on the similar lines. The longer *λ* with the same *A*_{0} corresponds to the wider regions of the low potential energy with the same barrier height. This means the confinement effect reduces with increasing *λ*. Hence, the threshold energy decreases with increasing *λ*. In other words threshold energy increases with increasing 1/*λ*.

The presence of pseudobands and pseudogaps, a feature of only rough surfaces, is also due to alternate regions of low and high potential energy seen by the electron while crossing the channel. The energy eigenvalues for an electron in such a periodic potential energy profile are given by Kronig Penny model. The solution of the Schrodinger equation for such a potential energy profile shows the presence of energy bands and energy gaps. The investigations on the bandwidth (BW) of the pseudobands with respect to roughness parameters is important as the width of the pseudoband determines the area under the transmission curve and play an important role in determining the total current *I*_{net}. It is seen from Figure 4a, b that BW of the pseudobands decreases with increasing *A*_{0} and *λ*. These trends of BW on roughness parameters can be understood by looking, again, at the Kronig Penny model. According to Kronig Penny model, if the barrier heights and widths of the periodic potential profile are high, then the BW's are small. This is so because higher barrier heights and widths give smaller tunnelling probability. As the tunnelling probability becomes smaller the BW reduces. Along the same lines increasing *A*_{0}, which corresponds to higher barrier heights, and increasing *λ*, which corresponds to larger barrier widths of the pseudoperiodic potential profile, would give rise to smaller BW's. Thus BW decreases with increasing *A*_{0} and *λ*.

### 3.2. Effect of roughness on total current *I*_{net}

As *I*_{net} represents the current only due to temperature gradient with *V*_{applied} = 0, it is given by the intercept with the *y*-axis of *I*-*V* characteristics. It is observed from Figure 3a that *I*_{net} decreases with increasing amplitude but Figure 3b shows that there is no definite trend with *λ*. These features can be explained by considering that total current depends on the bandwidth of the pseudobands and their occupancy. For increasing *A*_{0} as seen from Figure 4b, threshold energy increases and BW of the pseudobands decreases. The smaller bandwidths along with their presence at higher energy sides result in the decrease in occupancy of these bands. This, in turn, decreases the total no. of carriers contributing to the current. Thus, *I*_{net} decreases with increasing *A*_{0}. The trends of *I*_{net} with *λ* are complicated. As seen from Figure 4a though the pseudobands are present at lower energy sides with increasing *λ* but their BW decreases. Thus, on one hand their presence at lower energy sides increases the occupancy on the other hand smaller bandwidths decrease the total occupancy which is obtained by summing the occupancy from all energy levels. These two competing features complicate the trends of *I*_{net} with *λ*. Figure 4a shows that BW for *λ* = 3.3 nm is smaller than *λ* = 2.5 nm but since the band is present at lower energy sides its occupancy is more. This results in more current for *λ* = 3.3 nm than for *λ* = 2.5 nm. For *λ* = 5 nm the BW is so small that the current is the least.

### 3.3. Effect of roughness on Seebeck voltage *V*_{Seebeck}and Seebeck coefficient

As already mentioned *V*_{Seebeck} increases with increasing *A*_{0} and 1/*λ*. This trend can be understood by examining the circumstances for which *I*_{total} (= *I*_{net}) becomes zero. *I*_{total} is the sum of current from all energy levels, i.e. *I* = ∫*I*(*E*)*dE*. Figure 6c gives *I*(*E*) versus *E* for *λ* = 2.5 nm and *A*_{
0
} = 0.5 nm. It is seen that below *E*_{null} the current is +ve and above *E*_{null} the current is -ve. Hence, the total current obtained by the sum of current from all energy levels is small. We need to check when does it become zero. As seen from Equation (7), *I*(*E*) is product of transmission *T*(*E*) and $\left({f}_{1\mathrm{D}}^{\mathsf{\text{cold}}}-{f}_{1\mathrm{D}}^{\mathsf{\text{hot}}}\right)$, i.e. difference in the occupancy of cold ${f}_{1\mathrm{D}}^{\mathsf{\text{cold}}}$ and hot junctions ${f}_{1\mathrm{D}}^{\mathsf{\text{hot}}}$. Since *T*(*E*) ≥ 0 for all energies therefore current changes sign whenever $\left({f}_{1\mathrm{D}}^{\mathsf{\text{cold}}}-{f}_{1\mathrm{D}}^{\mathsf{\text{hot}}}\right)$ changes sign as seen from Figure 6b. For those energy levels in which the number of electrons at the cold side are more than the number of electrons at the hot side, the net current flows from cold to hot and $\left({f}_{1\mathrm{D}}^{\mathsf{\text{cold}}}-{f}_{1\mathrm{D}}^{\mathsf{\text{hot}}}\right)$ is +ve and *I*(*E*) is +ve. These are the energy levels below *E*_{null}. It is vice versa for the energy levels above *E*_{null}. At *E*_{null}, $\left({f}_{1\mathrm{D}}^{\mathsf{\text{cold}}}-{f}_{1\mathrm{D}}^{\mathsf{\text{hot}}}\right)$ = 0. This *E*_{null} is a function of *V*_{applied} as occupancies of energy states at cold side, ${f}_{1\mathrm{D}}^{\mathsf{\text{cold}}}$ and hot side, ${f}_{1\mathrm{D}}^{\mathsf{\text{hot}}}$, depend on *V*_{applied}. *E*_{null} shifts towards higher energy sides with increasing *V*_{applied}. We apply voltage to shift *E*_{null} to such an extent that total current coming from energy levels below *E*_{null} balances the total current coming from energy levels above *E*_{null} so that the total current becomes zero.

To understand how does *E*_{null} shifts towards higher energy states with increasing *V*_{applied}, we observe that as the threshold energy increases the electrons occupy higher energy states and the current contributions come from higher energy levels. The net current *I*_{net}, which flows from hot side to cold side is also made of electrons flowing in higher energy levels. We need to apply voltage, *V*_{Seebeck}, to nullify this current *I*_{net}. Thus, we need to increase the potential of the cold side so that electrons at the cold side occupy higher energy states to balance this current. Thus, as the threshold energy increases we need to apply higher *V*_{Seebeck} to balance the net current. Also, since the occupancy of higher energy states at the cold side increases with increasing *V*_{Seebeck} the *E*_{null}, for which $\left({f}_{1\mathrm{D}}^{\mathsf{\text{cold}}}-{f}_{1\mathrm{D}}^{\mathsf{\text{hot}}}\right)=0$, shifts towards higher energy sides. Since the threshold energy increases with increasing *A*_{0} and 1/*λ*, hence Seebeck voltage *V*_{Seebeck} increases with the increasing roughness amplitude and frequency. Since Seebeck coefficient is obtained by taking the ratio of *V*_{Seebeck} with temperature difference Δ*T*, i.e. *S* = *V*_{seeback}/Δ*T*, therefore the increase of *V*_{Seebeck} for fixed Δ*T* causes the increase in Seebeck coefficient with the similar trends. It is seen from Tables 1 and 2 that Seebeck coefficient increases with increasing *A*_{0} and 1/*λ*. The dependence of Seebeck coefficient on roughness parameters suggest that it can be tailored by choosing the appropriate roughness parameters.

For the devices if *A*_{0} is increased at constant *λ* it is observed that though *I*_{net} decreases but still *V*_{Seebeck} increases. This happens because BW decreases and threshold energy increases. Smaller BW's at higher energy sides reduce occupancy and hence *I*_{net}. Though the current reduces, it is made up of electrons present in higher energy levels and as already explained higher voltages are required to nullify the current from higher energy states.

The discussions show that it is the increase in threshold energy which causes the current to flow in the higher energy levels and hence results in increase in Seebeck coefficient. This result implies that any physical geometry whether periodic or aperiodic which results in an increase in threshold energy will show an increase in Seebeck coefficient. It also implies that if the threshold energy dependence on roughness parameters become feeble then Seebeck coefficient will saturate, i.e. it will not change much with the change in roughness parameters. The detailed investigations on the above two implications are in progress.

## 4. Conclusions

The surface roughness causes the net current *I*_{net}, produced due to the temperature gradient, to decrease with increasing roughness *A*_{0} but trends with *λ* are complicated. The voltage *V*_{Seebeck}, required to nullify this *I*_{net} shows a definite trend. It increases with increasing *A*_{0} and 1/*λ*. Increase in *V*_{Seebeck} results in increase in Seebeck coefficient. The reasons for the above trends have been attributed to the varying periodic thickness seen by the electron crossing the channel. The periodically varying thickness corresponds to varying confinement and hence periodically varying subband energy. The subband energy corresponds to the potential energy seen by the electron. For the electron moving in a pseudoperiodic potential, the transmission spectra shows the following features: (1) The threshold energy at which the transmission of current starts is shifted towards higher energy sides. It increases with increasing *A*_{0} and 1/*λ*. (2) There are regions of pseudobands and pseudogaps. The BW of the pseudobands decreases with increasing *A*_{0} and *λ*.

The above features of transmission spectra result in following trends for *I*_{net} and *V*_{Seebeck}: (1) It is observed that BW decreases and shift towards the higher energy sides for increasing *A*_{0}. Both these features result in reduced occupancy of the band and hence total carrier concentration. Since *I*_{net} depends on the total carrier concentration thus *I*_{net} reduces with increasing *A*_{0}. (2) It is observed that pseudobands are present at lower energy sides with increasing *λ* but their BW decreases. Thus, on one hand the presence of bands on the lower energy side increases the occupancy and on the other hand the smaller bandwidth reduces the total occupancy. These two competing features give rise to complicated trends of *I*_{net} with *λ*. (3) The trends of threshold energy determine the trends of *V*_{Seebeck} with roughness parameters. The increase in threshold energy results in the electrons to occupy higher energy states. Thus the net current *I*_{net} flowing from hot to cold end is made up from higher energy levels. To nullify this current, the occupancy of higher energy states at the cold end is increased by applying higher voltages. Thus higher threshold energy requires higher applied voltages, *V*_{Seebeck}, to nullify the current. Since threshold energy increases with increasing *A*_{0} and 1/*λ*, hence *V*_{Seebeck} and Seebeck coefficient increases vice versa. (4) The dependence of Seebeck coefficient on roughness parameters suggest that it can be tailored by choosing the appropriate roughness parameters.

## References

- 1.
Jang M, Park Y, Jun M: The characteristics of Seebeck coefficient in silicon nanowires manufactured by CMOS compatible process.

*Nanoscale Res Lett*2010, 5: 1654–1657. 10.1007/s11671-010-9690-2 - 2.
Shakouri A: Recent development in semiconductor thermoelectric physics and materials.

*Annu Rev Mater Res*2011, 41: 399–431. 10.1146/annurev-matsci-062910-100445 - 3.
Hicks LD, Dresselhaus MS: Effect of quantum-well structure on thermoelectric figure of merit.

*Phys Rev B*1993, 47: 12727–12731. 10.1103/PhysRevB.47.12727 - 4.
Kim R, Datta S, Lundstorm MS: Influence of dimensionality on thermoelectric device performance.

*J Appl Phys*2009, 105: 034506–034511. 10.1063/1.3074347 - 5.
Demchenko DO, Heinz PD, Lee B: Determining factors of thermoelectric properties of semiconductor nanowires.

*Nanoscale Res Lett*2011, 6: 502. 10.1186/1556-276X-6-502 - 6.
Andrei P, Mayergoyz I: Quantum mechanical effects on random oxide thickness and random doping induced fluctuations in ultrasmall semiconductor devices.

*J Appl Phys*2003, 94: 7163–7172. 10.1063/1.1625084 - 7.
Markussen T, Jauho AP, Brandbyge M: Surface decorated silicon nanowires: a route to high-ZT thermoelectrics.

*Phys Rev Lett*2009, 103: 055502–055505. - 8.
Hochbaum AI, Chen RK, Delgado RD, Liang WJ, Garnett EC, Najarian M, Majumdar A, Yang PD: Enhanced thermoelectric performance of rough silicon nanowires.

*Nature*2008, 451: 163–168. 10.1038/nature06381 - 9.
Zianni X: Eff

*icient*thermoelectric energy conversion on quasi-localized electron states in diameter modulated nanowires.*Nanoscale Res Lett*2011, 6: 286. 10.1186/1556-276X-6-286 - 10.
Low T, Ansari D: Signature of quantum transport through two dimensional structures with correlated and anticorrelated interfaces.

*Phys Rev B*2008, 78: 165301–165309. - 11.
Goodnick SM, Ferry DK, Wilmsen CW, Liliental Z, Fathy D, Krivanek OL: Surface roughness at the Si(100)-SiO

_{2}interface.*Phys Rev B*1985, 32: 8171–8186. 10.1103/PhysRevB.32.8171 - 12.
Schliemann A, Worschech L, Forchel A, Curatola G, Iannaccone G: Universal signature of ballistic transport in nanoscale field effect transistors.

*Tech Dig Int Electron Devices Meet*2004, 1039–1042. - 13.
Bulusu A, Walker DG: Quantum modelling of thermoelectric properties of Si/Ge/Si superlattices.

*IEEE Trans Electron Dev*2008, 55: 423–429. - 14.
Wu Y, Childs PA: Conductance of Graphene nanoribbon junctions and the tight binding model.

*Nanoscale Res Lett*2011, 6: 62. - 15.
Datta S:

*Quantum Transport: Atom to Transistor*. New York: Cambridge University Press; 2005. - 16.
Uchida K, Watanabe H, Kinoshita A, Koga J, Numata T, Takagi SI: Experimental study on carrier transport mechanism in ultrathin-body SOI n and p-MOSFETs with SOI thickness less than 5 nm.

*Tech Dig-Int Electron Devices Meet*2002, 47–50.

## Acknowledgements

SN acknowledges the financial support from Department of Atomic Energy, BRNS, Mumbai, India (grant no.2010/34/47/BRNS/2317). MK thanks Prof. Arun Kumar (IIT Delhi) for help and support during the preparation of manuscript.

## Author information

### Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors' contributions

MK carried out the NEGF calculations and implemented the code. AB and SN conceived of the study. AB participated in the NEGF calculations and implementation of the code. SN, with the perspective of designing the device in his lab, participated in the discussions of theoretical results and provided the information of the related experimental works. All authors read and approved the final manuscript.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

**Open Access**
This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (
https://creativecommons.org/licenses/by/2.0
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Kumar, M., Bagga, A. & Neeleshwar, S. Tailoring of Seebeck coefficient with surface roughness effects in silicon sub-50-nm films.
*Nanoscale Res Lett* **7, **169 (2012). https://doi.org/10.1186/1556-276X-7-169

Received:

Accepted:

Published:

### Keywords

- Threshold Energy
- Seebeck Coefficient
- Roughness Parameter
- Schrodinger Equation
- Cold Side