Original scientific paper Informacije ^efMIDEM A Innrnal of Journal of Microelectronics, Electronic Components and Materials Vol. 44, No. 2 (2014), 119 - 125 Impact of Downscaling on Analog/RF Performance of sub-100nm GS-DG MOSFET P. K. Sahu, S. K. Mohapatra, K. P. Pradhan Nano Electronics Laboratory, Department of Electrical Engineering, National Institute of Technology (NIT), Rourkela, Odisha, India Abstract: This paper presents a systematic study to show the impact of channel length on the Analog/RF performances of gate stack (GS) silicon on insulator (SOI) architecture. The downscaling of channel length becomes the biggest challenge to maintain higher speed, low power and better electrostatic integrity for each generation. This investigation is done to find out the potential of the channel length in view of analog and RF performance measures of sub-100nm GS-double gate (DG) MOSFETs. The threshold voltage (Vth) is made constant by tuning the gate metal work function while downscale the channel length (L). The impact of channel length variation on subthreshold slope (SS), drain induced barrier lowering (DIBL), transconductance (gm), output conductance (gd), early voltage (VEA), transconductance generation factor (TGF), intrinsic gain (AV), cut-off frequency (fT), transconductance frequency product (TFP), gain frequency product (GFP) and gain transconductance frequency product (GTFP) are rigorously examined. It is shown that gate stack design results in higher cut-off frequency along with a broader analog 'sweet spot' in nanoscale MOSFETs thus offering better possibilities for analog/RF scaling below 50nm. For shorter gate length devices (L=30nm), the design results in an impressive 69.10% improvement in fT along with 36.31% enhancement in 'sweet spot' as compared to L=60 nm. The study generates an optimized channel length of L=40 nm for the designed device dimension in connection with the analog and RF performance for circuit design. Keywords: Gate Stack (GS), DG-MOSFETs, Metal Gate Technology, Analog/RF FOMs, Sweet Spot Vpliva pomanjševanja na analogne/RF lastnosti pod-100 nm GS-DG MOSFETa Izvleček: V članku je predstavljana sistematična študija vpliva dolžine kanala na analogne RF lastnosti arhitekture večplastnih vrat (GS) silicija na izolatorju (SOI). Največji izziv pri krajšanju kanala je ohranjanje visoke hitrosti, nizke moči in boljše elektrostatične celote. Raziskava odkriva potenciale dolžine kanala pri analognih in RF lastnostih. MOSFET-ov z dvojnimi vrati. Konstantnost pragovne napetosti pri krajšanju kanala se je ohranjala s spreminjanjem delovne funkcije kovinskih vrat. Raziskan je vpliv krajšanja kanala na podpragovni naklon (SS), ponorno vzbujano nižanje bariere, transkonduktanco (gm), izhodno prevodnost (gd), zgodnjo napetost (VEA), generacijski faktor transkonduktance (TGF), intrinsično ojačenje (AV), frekvenco reza (fT), produkt transconduktančne frekvence (TFP), produkt frekvence ojačenja (GFP) in produkt frekvenc ojačenja in transkonduktance (GTFP). Izkazalo se je, da večplastna vrata omogočajo višje frekvence reza skupaj s širšim območjem najboljšega delovanja v nanodimenzijskih MOSFET-ih, kar omogoča krčenje na 50 nm. Pri krajših dolžinah vrat (L=30nm) oblika izkazuje impresivno 69.10 % izboljšanje fT in 36.31 % izboljšanje območja najboljšega izplena v primerjavi z L = 60 nm. Izkazalo se je, da je optimalna dolžina kanala 40 nm. Ključne besede: večplastna vrata, DG-MOSFET, tehnologija kovinskih vrat, Analogni/RF FOM, območje najboljšega izplena * Corresponding Author's e-mail: kp2.etc@gmail.com 1 Introduction The use of low power and high frequency operated devices are having a high priority for future electronic applications. Silicon on Insulator (SOI) devices are excellent candidates as an alternative for the conventional bulk CMOS [1-2]. Advanced MOSFET structures such as ultra-thin body (UTB) SOI double gate (DG) MOSFET can be scaled more aggressively than bulk Si structure [3]. A double gate structure fabricated on SOI wafer has been utilized in CMOS technology due to their excellent scaling capability, outstanding Short Channel Effects (SCEs) immunity, high current drivability (Ion) and transconductance (gm) and lower leakage current (Ioff) as compared to the bulk MOSFETs [4-10]. Channel length scaling is limited by the ability to control off- state leakage current due to quantum tunneling and thermionic emission between the source and drain. The analog and RF circuits with a digital CMOS technology suffers from many challenges. More over the technology is optimized for digital design and the devices are characterized for current drive and gate delay. For System on Chip (SoC) and System in Package (SiP) applications, optimization of the devices becomes more challenging. SiP refers to a semiconductor device that incorporates multiple readily available chips into a single package; while SoC refers to a device that incorporates a system of differently functioning circuit blocks on a single silicon chip. So it is required to enhance the performance for digital and analog/RF circuit applications [11-18]. Major semiconductor companies such as IBM, RFMD, Honeywell, OKI, etc., have already produced several products for the telecommunication market based on SOI RF technologies. In this work, a study has been made with designing different device cases by varying channel length (L) keeping all other process parameters constant for high-k/ metal gate (HKMG) DG-MOSFETs. In our analysis, Analog/RF devices designed around the 'sweet spot' (compromise between power and speed performance, and linearity) signified by the peak of transconductance to current ratio (gm/ID) and cut-off frequency (fT) product, typically characterized by a gate voltage of 0.35 V. The device performance is investigated keeping in view for analog and RF circuit application. For the analysis, the threshold voltage (Vth) is maintained at a constant value for all cases by tuning the gate metal work function between 4.6 eV to 4.7 eV. The Vth is achieved a constant value of 0.2V for all channel length device cases. Section 2 describes the device description that includes all the dimensions, materials and doping concentrations of HKMG DG-MOSFET. Section 3 analyses the physics of the device using device numerical simulations and models activated for simulation. In section 4, various performance metrics of the device including the short channel effects (SCEs) like SS, DIBL, and the important analog and RF FOMs such as g^, g^, TGF, Vea, A,, Cg^, Cg^, fT, TFP, GFP, GTFP have been closely observed. Finally, the concluding remarks have been included in section 5. 2 Basic Structure For the purpose of study a planar DG SOI n-channel transistor is considered. For the simulation a symmetric device structure model has been designed. Four differ- ent channel lengths (L= 60nm, 50nm, 40nm, and 30nm) are chosen for the device to analyze the device perfor- mance. The schematic simulated structure is shown in Fig. 1, where the silicon channel is surrounded above and below by two layer of oxide, SiO2 and HfO2. Metal gate technology is attractive because, it eliminates the poly-Si gate depletion effect and the associated degradation in transistor performance. To achieve a constant Vth, the work function of the metal gate is tuned in between 4.6 eV to 4.7 eV while varying the channel length of the device.The doping profile for channel (p-type 1016 cm-3) and source, drain (n-type 1020 cm-3) are set. For a better comparison of analog/RF FOMs, the Vth is maintained at a constant value 0.2 V while varying the metal gate work function at VDS=0.1 V. Figure 1: Schematic structure of Double Gate N-MOS-FET 3 Device Simulation The 2-D numerical device simulator ATLAS is employed to simulate the planner DG-MOSFET with high-k/metal gate. According to International Technology Roadmap for Semiconductors (ITRS), the drain bias has been fixed at VDD = 1.0 V [19]. In order to study the analog performance, the simulation is performed with VDS = 0.5 V (which is half of the supply voltage i.e., VDD/2) [20-21] with variable VGS = 0 V to 1.0 V. Threshold voltage (Vth) is extracted using constant current (ID=10-6 A/ ^m) definition from the ID-VGS transfer characteristic at VDS=0.1 V. The threshold voltage is used to find out the gate overdrive voltage (VGT=VGS-Vth). The gate over drive voltage is an important property of amplifier circuit as it decides the region of operation. The increment of VGT increases drain current until saturation. Hence in our investigation all the analog and RF analysis has been done against VGT [22- 23]. For more accuracy, it is essential to account for the mobility degradation inside the inversion layers. The degradation normally occurs as a result of higher surface scattering near the semiconductor to insulator interface. Hence, during simulation, the inversion-layer Lombardi constant voltage and temperature (CVT) mobility model has been considered. The Shockley-Read-Hall (SRH) generation and recombination parameters simulate the leakage currents existing due to thermal generation are incorporated in the simulation. At thermal equilibrium, a semiconductor lattice obeys Fermi-Dirac statistics. The use of Boltz-mann statistics is normally justified in semiconductor device theory. However, Fermi-Dirac statistics are necessary to include for certain properties of very highly doped (degenerate) materials. The model Fermi-Dirac uses a Rational Chebyshev approximation that gives results close to the exact value. In the simulation all the junctions of the structure are considered to be of abrupt. Suitable empirical parameter ß is selected to calibrate the drift diffusion transport model [20]. The biasing conditions are granted as per the room temperature T=25 oC. Furthermore, we have chosen two numerical techniques Gummel and Newton to obtain the solutions [24]. 4 Results and Discussion In order to analyze the impact of channel length of the device on the performance, simulation is done for four different channel lengths as 60 nm, 50 nm, 40 nm, and 30 nm. f) (1 fl.2 0.4 0.6 a.8 Gate VoltaKe, (V) (a) (b) 0.2 0.4 0.6 0.8 Drain Voltage, (V) Figure 2: Drain current (ID) in both linear and log scale as a function of gate to source voltage (VGS) (b) Output conductance (gd) and drain current (ID) with respect to drain to source voltage (^J for different channel lengths. However, for analog and RF performance comparison among all the channel length cases, Vth kept constant (0.2 V) by adjusting the metal gate work function between 4.6 eV and 4.7 eV. The L-V transfer character- D GS istics both in linear and log scales have been shown in Fig. 2 (a) for different configurations at VDS=0.5 V. The leakage current (Ioff) is quit constant for three different channel length except L=30 nm. In the short channel device the Ioff increases due to random motion of charge carriers. As channel length decreases, it gives rise to high drain current because of the relation Id^1/L. However, from the log scale, the leakage current is also prominent for lower channel lengths.Drain current (ID) and output conductance (gd) against drain to source voltage (VDS) for different cases at VGS=0.5 V are presented in Fig.2 (b). As per the Fig.2 (b), the drain current is increasing with decrease in channel length which in turn makes; the gd high for lower L devices as gd=3lD/3VDS. As we know from the literatures that gain and early voltage are inversely proportional to output conductance, so the device having lower L gives higher gd which comprises lower gain and early voltage of the device. Transconductance generation factor (TGF=gm/ID) and transconductance (gm) as a function of gate over drive voltage (VGT=VGS-Vth) are presented in Fig.3 (a). From the figure, it is clear that as the channel length decreases the gm value is increasing because of high drain current. The gm/ID ratio demonstrates how efficiently the current is used to achieve a certain value of transconductance. The advantage of high transconductance-to-drain ratio is the realization of circuits operating at low supply voltage. As shown in the figure, gm/ID is maximized towards the subthreshold region of device operation. From the Fig.3 (a), it is clear that the structure having channel length 60 nm shows higher gm/ID ratio as compare to others and it decreases as channel length decreases. Fig.3 (b) shows the variation of the early voltage (VEA) and intrinsic gain (AV) as a function of gate over drive voltage (VGT) for different channel lengths. For better analog performance the VEA and AV should be as high as possible. An enormous improvement is observed in VEA for channel length L=60 nm as compared to others. The intrinsic gain of the device, which is a ratio of transconductance and output conductance for various channel lengths is plotted against gate voltage (VGS) for Vds=0.5 V is shown in Fig.3 (b). The intrinsic gain (AV=gm/ gd) is a valuable figure of merit for operational transconductance amplifier. From the figure, the device having channel length 60 nm gives highest gain from others and it decreases as the channel length decreases. Both the extracted values of SS and calculated values of DIBL for different channel lengths are tabulated in Table 1. 45-1 40- 35- 30- V 25- a t; E 20- M nT 15- Ü 10- H 5- 0- \ V„, = 0.5V - L=60 nm - - L=50 nm - L=40 nm - - L=30 nm 3.0x10"' H 3 2.5x10-' 3 o o 2.0x10"' 3 D. C O 1.5x10"' B 3 O 01 1.0x10"' Oq 5.0x10"' 2 0,0 0.2 0.4 \r(V) 0.6 0.8 (a) > ^20 M w 16- u I '' "o > 8 ^L-60 nm ^L=50nm —^ L=40 nm d' Vos = 0-5 v 50 3 40 2. 30 0.0 0.2 0.4 Vor (V) 0.6 0.8 o O p 20 ffQ (ff 10 Q. 03 (b) Figure 3: (a) Transconductance generation factor (TGF) and transconductance (gm) (b) Early voltage (VEA) and intrinsic gain (AV) as a function of gate over drive voltage (VGT) for different channel lengths. The DIBL calculation is performed for Vth at VDS=0.1 V and VDS=1.0 V. From the table, it is clear that tlie SS value increases as channel length decreases and it is high for channel length 30 nm. Similarly, the DIBL value also increases as channel length decreases and it gives a maximum value for channel length of 30 nm. These two parameters are very important for short channel effects, which should be minimized. The maximum values for gm, gd, VEA, AV, TGF are also tabulated in Table 1. ID increases for lower channel length devices which consequently increases gm values for devices having lower value of L. However, because of high gd, the VEA, and AV becomes lower as L decreases. Coming from L=60 nm to 30 nm, the DIBL and SS values are more prominent for lower channel length devices and also the TGF and Gain are decreases as L decreases. Figure 4: Gate to source capacitance (Cgs) and gate to drain capacitance (Cgd) as a function of gate over drive voltage (VGT) for different channel lengths. Fig. 4 shows the intrinsic capacitances (Cgs & Cgd) as a function of VGT. As shown in figure, the intrinsic capacitance parameters increase swiftly in the super threshold region. This is because of the increase in the fringing field lines emanating from the gate edges. The device having channel length 60nm shows higher values of intrinsic capacitances (both Cgs & Cgd) and it decreases as channel length decreases. From Fig.5 (a), the variations of cut off frequency (fT=gm/2n(Cgs+Cgd)) and gain transconductance frequency product (GTFP=AV*TGF*fT) can be observed with respect to VGT for different values of channel lengths. Here, the value of fT obtained for the device having low channel length is higher and it gradually decreases as the channel length decreases. fT is inversely proportional to the intrinsic capacitances (Cgs & Cgd). So, fT value is low due to high capacitance values for higher channel length devices. It is interesting to see that the device having Table 1: Electrostatic & Analog performances for different values of channel lengths D (m m SS (mV/dec-ade) Transconductance, gm (mS) Output Conductance, gd (^S) Early Voltage, VEA (V) Gain, AV (dB) TGF (V-1) Vds=0.5V L=60 20.71 62.31 2.81 12.94 24.130 46.742 43.637 L=50 21.12 63.12 2.93 16.57 20.737 44.969 42.840 L=40 37.24 65.95 3.05 24.53 15.787 41.921 40.639 L=30 49.88 68.35 3.18 54.84 8.499 35.458 34.669 Table 2: RF performances for different values of channel lengths Channel Length (nm) Cgs (fF) Cgd (fF) ft (GHz) GFP (GHz) TFP(GHz/V) GTFP (GHz/V) L=60 1.207 0.484 306.78 4.57*103 3.69*103 1.12*105 L=50 1.041 0.449 358.00 5.27*103 4.22*103 1.25*105 L=40 0.874 0.416 425.80 6.12*103 4.76*103 1.33*105 L=30 0.706 0.384 518.79 7.02*103 5.03*103 1.20*105 channel length L=40 nm shows a higher GTFP value as comparison to others. This is due to the reduction in peak electric field, lower output conductance of the device having channel length 40nm. The GTFP value is very low for L=30nm. 1.5x10" Vct(V) (a) 2.0x10 Vc,T(V) (b) Figure 5: (a) Cut off frequency (fT) and gain transconductance frequency product (GTFP) (b) Gain frequency product (GFP) and transconductance frequency product (TFP) as a function of gate over drive voltage (VGT) for different channel lengths. The product of gm/ID and fT represents a trade-off between power and bandwidth and is utilized in moderate to high speed designs. Fig.5 (b) gives the gain frequency product (GFP=AV*fT) and TFP against gate over drive voltage (VGT) for different values of channel lengths. From the figure, the value of GFP increases as channel length decreases and reaches utmost for the device having channel length 30nm. The same figure also shows transconductance frequency product (TFP) as a function of VGT for different values of channel lengths. Where the bandwidth is flexible and part of the overall optimization process, it is important to consider the product of gm/ID and fT as shown in Fig. 5 (b). This figure of merit exhibits a 'sweet spot' i.e., peak value for all device cases. Peak values for the product (TFP) increases from 3.69*103 GHzV-1 to 4.76*103 GHzV-1 when gate length (L) is reduced from 60nm to 30nm. The maximum values for TFP i.e., 'sweet point' is achieved at a gate voltage of 0.35 V which is 0.15 V more than from the threshold voltage. From the figure it is clear that the device having channel lengths 40nm and 30nm gives higher TFP values as comparison to others. All the extracted values for analog/RF FOMs are plotted in Table 2 for different values of channel lengths. It is clear from the Table that, while the gate length is reduced the RF FOMs like fT, GFP and TFP are also increased because of high drain current which results in higher gm values for shorter gate length devices. However, the improvement in GTFP, which is a unique and major FOM, with L downscaling reduces below 40 nm due to short channel effects. From Table 1, the AV values are obtained around 41.921 dB and 35.458 dB for 40nm and 30nm respectively whereas the gm/ID shifts lower to 34.669 V-1 for 30nm technology from 40.639 V-1 for 40nm technology. The GTFP is nothing but the product of three parameters i.e., TGF, gain and frequency, so the peak values are more for 40 nm technology as compare to 30 nm technology. 5 Conclusion A detailed study on the analog/RF performance measure of GS-DG MOSFETs has been presented. It has been demonstrated that analog/RF performance metrics can be significantly enhanced by down scaling the channel length in terms of (a) analog 'sweet spot' and (b) compromise between gain and cut-off frequency i.e., GFP. However, the improvement in GTFP with L down scaling continues up till 40nm where a peak value of 1.33*105 GHzV-1 is attained. Beyond 40 nm, peak GTFP value reduces to 1.20*105 GHzV"1 in 30 nm GS-DG MOS-FET due to short channel effects (peak gm/ID~34.669 V-1 and gain~35.458 dB, Table 1). So, from tiie results it is clear that all studied parameters are more sensitive to the channel length (L) of the device. The calculated and simulated results demonstrate that the optimized value of "L" will be 40 nm for the chosen device dimension. The values of the parameters like DIBL=37.24 mV/V, SS=65.95 mV/decade, Ion=1.912 mA, AV=41.92 dB, fT=425.80 GHz, GTFP=1.33*o1n05 GHz/V are achieved for L=40 nm. The other device parameters could also be properly chosen for further downscaling of the transistor. Due to lack of fabrication facilities, we can't validate our simulation results with the experimental results in literature. However, as our simulation is calibrated with experimental results, so we can give a formal assurance that degradation of studied parameters may not be expected when implementing MOSFETs in a real chip/ SoC. 6 References 1. D.J. Frank, R.H. Dennard, E. Nowak, P.M. Solomon, Y. Taur, H. Wong, "Device scaling limits of Si MOSFETs and their application dependencies', Proc. IEEE, Vol. 89 , No.3, pp. 259-288, Mar., 2001. 2. Leland Chang, Y-K Choi, D Ha, P Ranade, S Xiong, J Bokor, C Hu, And T-J King, "Extremely Scaled Silicon Nano-CMOS Devices", Proceedings of The IEEE, vol. 91, no. 11, pp. 1860 - 1873, Nov., 2003. 3. Isabelle Ferain, Cynthia A. Colinge, Jean-Pierre Col-inge, " Multigate transistors as the future of classical metal-oxide-semiconductor field-effect transistors', Nature, vol.-479, pp.310-316, Nov., 2011. 4. K. Suzuki, Y. Tosaka, T. Tanaka, H. Horie, Y. Arimoto, "Scaling theory of double-gate SOI MOSFET's", IEEE Transactions on Electron Devices, Vol. 40, No.12, pp.2326-2329, Dec., 1993. 5. Scott Thompson, Paul Packan, Mark Bohr, "MOS Scaling: Transistor Challenges for the 21st Century" Intel Technology Journal, vol. 2, issue 3, 1998. 6. Colinge J.P. "Multiple-gate SOI MOSFETs" Solid State Electronics, Vol. 48, No.6, pp.897-905, June, 2004. 7. C. Wann, K. Noda, T. Tanaka, M. Yoshida, and C. Hu, "A Comparative study of advanced MOSFET concepts", IEEE Transactions onElectron Devices, Vol. 43, pp. 1742, Oct. 1996. 8. H.-S. P. Wong "Beyond the conventional transistor" IBM J. Res. & Dev., Vol. 46, No. 2/3, Mar./May, 2002. 9. G. K. Cellera, Sorin Cristoloveanu, "Frontiers of sil-icon-on-insulator", Journal of Applied Physics, Vol., No.-9, pp.4955-4978, May, 2003. 10. Chaudhry, M.J. Kumar, "Controlling short-channel effects in deep submicron SOI MOSFETS for improved reliability: a review", IEEE Transaction on Device Material Reliability , Vol. 4, No.3, pp. 99-109, Mar, 2004. 11. Jean-Pierre Colinge, "Fully-Depleted SOI CMOS for Analog Application", IEEE Transaction on Electron Devices, Vol.-45, No.-5, pp.-1010-1016, May, 1998. 12. Behzad Razavi, "CMOS Technology Characteriza-tionfor Analog and RF Design", IEEE Journal Of SolidState Circuits, Vol.-34, No.-3, pp.-268-276, Mar., 1999. 13. Valeriya Kilchytska, A Neve, L Vancaillie, D Levacq, 5 Adriaensen, H V Meer, K De Meyer, C Raynaud, M Dehan, J P Raskin, D Flandre, "Influence of Device Engineering on the Analog and RF Performances of SOI MOSFETs", IEEE Transactions on Electron Devices, Vol. 50, No.3, pp. 577-588, Mar., 2003. 14. Lazaro, B. Iniguez, "RF and noise performance of double gate and single gate SOI"",Solid-State Electronics, Vol. 50, Issue 5, pp. 826-842, May, 2006. 15. Jean-Pierre Raskin, Tsung Ming Chung, Valeria Kilchytska, Dimitri Lederer, Denis Flandre, "Analog/ RF Performance of Multiple Gate SOIDevices: Wideband Simulations and Characterization""IEEE Transactions on Electron Devices, Vol. 53, No.5, pp. 1088-1095, May, 2006. 16. N. Mohankumar, B Syamal, C K Sarkar, "Influence of Channel and Gate Engineering on the Analog and RF Performance of DG MOSFETs", IEEE Transactions on Electron Devices , Vol. 57, No. 4, pp. 820826, Apr., 2010. 17. Angsuman Sarkar, Alok Kumar Das, Swapnadip De, Chandan Kumar Sarkar, "Effect of gate engineering in double gate MOSFETs for analog/RF applications", Microelectronics Journal, Vol. 43, pp.-873-882, July, 2012. 18. Rupendra Kumar Sharma, Matthias Bucher, "Device Design Engineering for Optimum Analog/ RF Performance of Nanoscale DG MOSFETs", IEEE Transactions on Nanotechnology, Vol.-11, no.-5, pp.-992-998, Sept., 2012. 19. The International Technology Roadmap for Semiconductors [Online]. Available: http://public.itrs. net/ reports.html. 20. K. P. Pradhan, S. K. Mohapatra, P. K. Sahu, D. K. Behera, '"Impact of High-k Gate Dielectric on Analog 6 RF Performance of Nanoscale DG-MOSFET", Microelectronics Journal , vol. 45, Issue 2, pp.-144-151, Feb., 2014. 21. Abhijit Mallik; Chattopadhyay, A., "Tunnel Field-Effect Transistors for Analog/Mixed-Signal System-on-Chip Applications,"" IEEE Transactions on Electron Devices, vol.59, no.4, pp.888,894, April, 2012. 22. N. Mohankumar, Binit Syamal, C.K. Sarkar, "Investigation of novel attributes of single halodual-ma-terial double gate MOSFETs for analog/RF applications'; Microelectronics Reliability,vol. 49, Issue 12,pp.-1491-1497, Dec., 2009. 23. Saurav Chakraborty, Abhijit Mallik, Chandan Kumar Sarkar, "Subthreshold Performance of DualMaterial Gate CMOS Devices and Circuits for Ul-tralow Power Analog/Mixed-Signal Applications", IEEE Transactions on Electron Devices , Vol. 55, No.3, pp. 827-832, Mar. 2008. 24. ATLAS manual: SILVACO Int. Santa Clara, 2008. Arrived: 18. 12. 2013 Accepted: 25. 03. 2014