SciELO - Scientific Electronic Library Online

 
vol.65 número1Characterization of high mobility inverted coplanar Zinc Nitride thin-film transistorsCO2 adsorption on a modified graphite surface with sodium dodecyl sulfate surfactants: a molecular dynamics study índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.65 no.1 México ene./feb. 2019  Epub 09-Nov-2019

 

Investigación

Effect of the paramagnetic to spin-glass phase transition on the fundamental absorption edge of MnIn2Se4 magnetic semiconducting compound

V. Sagredoa 

T. E. Torresb 

G.E. Delgadoc 

C. Rincónd 

aLaboratorio de Magnetismo, Departamento de Física, Facultad de Ciencias, Universidad de Los Andes, Merida 5101, Venezuela.

bInstituto de Nanociencia de Aragon & Laboratorio de Microscopias Avanzadas, Universidad de Zaragoza, C/Mariano Esquillor s/n, CP 50018, Zaragoza, Spain.

cLaboratorio de Cristalografía, Departamento de Química, Facultad de Ciencias, Universidad de Los Andes, Merida 5101, Venezuela.

dCentro de Estudios de Semiconductores, Departamento de Física, Facultad de Ciencias, Universidad de Los Andes, Merida 5101, Venezuela.


Abstract

The temperature dependence of the indirect and direct fundamental band gaps of the layered compound MnIn2Se4, that crystallizes in a rhombohedral defect structure with space group R3-mD3d5, was studied by optical absorption spectra. The data were analyzed in terms of current theoretical models that take into account the contribution of antiferromagnetic exchange interaction between spins to the shift of the energy gap E G with T in the vicinity of the critical paramagnetic to spin-glass phase transition. It was established that short- and long-range effect spin correlations dominate the contribution to this shift in the critical region below about 20 K, near the spin-glass freezing temperature T f , and noncritical one, between about 70 and 160 K, far from T f , respectively. An intermediate temperature region, compatible with the behaviour expected for a cluster-glass transition where a gradual freezing of the magnetic moments occurs, was also observed.

Keywords: Semiconductors; electronic band structure; exchange and superexchange; order-disorder effects

PACS: 71.20.Nr; 75.10.Nr

Introduction

The family of ordered-vacancy compounds (OVC’s) with the general formula MnB2IIIX4VI, with B: Al, Ga or In; and X: S, Se, or Te, where ⃞ symbolizes a cation vacancy, or simply MnB2IIIX4VI, has received considerable interest in the last three decades [1-12]. These magnetic compounds, which have three metal cations (one Mn and two B’s) around each anion site (X) while the fourth empty site forms an ordered array of vacancies, could exist in different crystal structures and have shown several notable physical and chemical properties associated with the crystallographic ordered array of these vacancies. In addition, the presence in these OVC’s of magnetic Mn+2 ions make them of great interest from the point of view of their magneto-optical properties. Thus, magnetic properties of MnIn2Te4, MnGa2Se4, MnIn2Se4; which crystallize in tetragonal defect-chalcopyrite (space group I4-2m) [2,3], defect-thiogallate (s. g. I4-) [7], and cubic spinel (s. g. Fd3-m) [12] structures at ambient conditions, respectively, and its alloys, have been studied in detail for several authors [2-8,13-15]. MnIn2Se4, a member of this family which crystallizes in a layered rhombohedral structure (s. g. R3-m) [16,17], has also been studied with considerable interest. Measurements of magnetic susceptibility as a function of temperature show antiferromagnetic behavior with paramagnetic Curie temperature of about −96 K and the formation of spin-glass-like phase below about 3.6 K [6,14,15]. The occurrence of this spin-glass transition is confirmed by ac in-phase magnetic susceptibility measurements close to the spin freezing transition [14]. Here, the typical downward shift in the freezing temperature T f , when the frequency is decreased, was observed. The zero-frequency transition temperature is T f = (3.58 ± 0.02) K. On the other hand, analysis of its optical absorption coefficient spectra [1,18] reveals that the lowest-energy-gap of this compound is indirect between parabolic bands that vary from 1.55 to 1.43 eV in the temperature range from 10 K to 295 K [18]. This indirect band-to-band transition arises from the uppermost Γ4(z) valence band at Γ point, to the conduction band minimum located at a lower-symmetry point different from Γ [12]. Additionally [18], two direct band-to-band transitions, beginning at 1.72 and 1.85 eV, at 295K, and at 1.82 and 1.96 eV, at 10 K, from the uppermost Γ4(z) heavy-hole and middle Γ5(z) light-hole bands, to the Γ1(s) conduction band, were observed in the high energy region. At energies around 2.2 eV, an internal transitions from 6A1 ground state to the 4 T 1 d-excited level of Mn+2 ion, was also observed from photoluminescence spectra.

A semiconductor property which is affected by the presence of Mn+2, especially at low temperatures, is the band gap E G . Theoretical studies have shown [19,20] that nearby a critical point, as for example the spin glass transition temperature T f , a magnetic contribution to the variation E G with temperature, due to the exchange interaction between the d electrons of Mn+2 ions and the electrons in the valence and conduction bands, should occur. For antiferromagnetic exchange interaction, as is the case of MnIn2Se4, an increase in E G around T f is expected [19]. In spite of these advances, to the best of our knowledge, a detailed study of the temperature dependence of E G for semiconductors where a paramagnetic to spin-glass transition occurs, has not been reported. This is of interest in relation to whether it is a true phase transition from a thermodynamic point of view or a gradual freezing of the magnetic moments. For this reason, in the present work we present a study of the temperature dependence of the fundamental band gaps in MnIn2Se4. This dependence is discussed in terms of the theoretical models reported in the literature that analyze the magnetic contribution to the shift of E G with T.

2. Crystal growth and experimental details

Bulk single-crystals of MnIn2Se4 were grown by chemical vapor transport technique from polycrystalline MnIn2Se4 which was synthesized from stoichiometric amounts of 5 N pure Mn, In, and Se, in an evacuated and closed ampoule. This was carried out by gradual heating of the mixture, subsequent melting, and gradual cooling through the freezing point. For the crystal growth, a 20 cm long and 1.5 cm inner diameter quartz stationary ampoule loaded with 1.5 g charge of the pre-synthesized material was used. I 2, with a concentration of 4 mg/cm3 of ampoule volume, was employed as transport agent in the reaction. The transport reaction was carried out in a two-temperature zone furnace during seven days between 800oC, in the highest-temperature region of the ampoule where the source pre-reacted material is initially placed, and 750oC, in the lowest-temperature region where the obtained bulk crystals are deposited. As-grown bulk crystals have typical dimensions of about 20-30 mm2 and show a dark color.

A SIEMENS D-5005 diffractometer, operated at 30 kV and 15 mA, with a Cu-target (λ−Cu Kα = 0.15418 nm) tube, a graphite monochromator and a scintillation detector, was used to register the powder X-ray diffraction pattern at room temperature. The intensity data were recorded in a θ/θ reflection mode from 10 to 100o with a step size of 0.02o and a scan speed of 25 s/step. Quartz was used as an external standard.

The absorption coefficient spectra α was obtained from the optical transmittance T r spectra through the relation [21] T r = [(1 − R)2 exp(−αx)]/[1 − R 2 exp(−2αx)], where x is the thickness sample and R ≈ 0.223 [18], is the reflectivity. T r spectra were measured by using an Acton Research Spectra P-500 spectrometer in the energy range from 1.2 to 2.4 eV. For the measurements of these spectra at various temperatures from 10 to 295 K, the samples were mounted on a cold finger of a Janis closed cycle refrigerator.

3. Powder X-ray diffraction analysis

To determine the crystal structure and unit cell parameters of MnIn2Se4, powder X-ray diffraction analysis was carried out. The structure was refined by means of the Rietveld method using the Fullprof program. For the refinement, atomic coordinates of MnIn2Se4 single-crystal reported by Range et al. [17], where used in the analysis. The Rietveld plot is shown in Fig. 1. The final refinement converged to the profile agreement factors: R p = 9.9 %, R wp = 10.5 %, R exp = 9.1%, S = 1.2. MnIn2Se4 crystallizes in the centrosymmetric space group R3-mD3d5, Z = 3, with unit cell parameters: α = 0.4055(1), c = 3.9467(2) nm, c/α = 9.733, and V = 0.5620(2) nm3. These values are in good agreement with those reported in [16,17] and also shown in Table I.

FIGURE 1 Observed (circles), calculated (solid line), and difference plot of the final Rietveld refinement of MnIn2Se4

TABLE I Unit cell parameters α, c, and c/α, reported in the literature for MnIn2Se4. 

Compound Structure α (nm) c (nm) c/α Ref
Rhomboedral 0.4051(1) 3.9464(2) 9.74 [16]
MnIn2Se4 (space group: 0.405(6) 3.94(9) 9.73 [17]
R3-m) 0.4055(1) 3.9467(2) 9.733 Present work

The crystallographic analysis confirms that MnIn2Se4 belongs to a rhombohedral defect structure with tetrahedral and octahedral cation-to-anion coordination, which can be expressed by the chemical formula [M]oct([M]tetr)2Se4, where [M] means either Mn and In cations distributed in a disordered arrangement. [M]oct being the disordered cations occupying the octahedral sites, six-fold coordinated by Se anions, while [M]tetr being these cations occupying the tetrahedral sites, four-fold coordinated by Se anions.

The unit cell diagram of MnIn2Se4 is shown in Fig. 2. A close-packed Se layer, with Mn and In cations at random in octahedral and tetrahedral voids, leaving the interlayers between h-stacked layers vacant, can be observed in this figure.

FIGURE 2 The rhombohedral structure of MnIn2Se4 with octahedral (purple) and tetrahedral (green) polyhedra. Mn, In, Se atoms and cation vacancies, are represented by green, orange, yellow and white spheres, respectively. 

4. Temperature dependence of the band gap energy

The nature of the fundamental absorption edge in MnIn2Se4 was examined by us in a previous article [18] by studying the optical absorption coefficient spectra α c of this compound at several temperatures between 10 and 295 K. It was found, at values of α c lower than about 2 × 103 cm−1, that the lowest energy gap is indirect. On the other hand, for values of α c greater than 2 × 103 cm−1, the absorption coefficient spectra show the presence of two direct band-to-band transitions at higher energies. These two processes, with characteristic energies EGDI and EGDII, correspond [18] to direct transitions from the uppermost Γ4(z) and the middle Γ5(z) valence bands, to the Γ1(s) conduction band of MnIn2Se4, respectively.

The shape of E G vs. T curve in non-magnetic semiconductors is approximately flat at low temperatures, but smoothly changes to a linear decrease in the high temperature region. However, as can be seen in Fig. 3, where the temperature variation of E GI , EGDI and EGDII energy gaps of MnIn2Se4 in the form of solid circles are plotted, this is not observed in the present case where going to low temperatures, the curves also shown a nearly linear variation with T.

FIGURE 3 Temperature variation of E GI , EGDI and EGDII energy gaps of MnIn2Se4. Solid circles (•) and broken lines represent the experimental points and the fit to the experimental curves in the high-temperature region (T > 150 K) by using Eq. (5), respectively. 

Such a behavior at low T, which has also been observed in other materials having antiferromagnetic exchange [22,24], is also consistent with the theoretical study [19] about the effects of spin fluctuations on the shift of the E G with temperature around a magnetic critical temperature T C . According to this study, for an antiferromagnetic exchange interaction between spins, a magnetic contribution is expected that should cause an increase ∆E GM in the band gap of the semiconductor in the vicinity of the spin-glass transition temperature T f .

To describe the behavior of E G near the critical freezing temperature T f in terms of a power law, the reduced temperature difference τ = |TT f |/T f , which is zero at the phase transition, is introduced in the analysis of the data [19]. For a semiconductor with antiferromagnetic exchange interaction, the theory indicates that ∆E GM satisfies a relation of the form

(d/dτ)ΔEGM=-Aτ-λ (1)

where λ is the critical exponent of the band-gap variation ∆E GM and A is a parameter that measures the strength of the total spin interaction on this variation.

Two different temperature regions should be taken into account in this analysis. The critical region, near the spinglass freezing temperature T f , where the main contribution to ∆E GM vs. T arises from short-range spin correlations, and the non-critical one, far from T f , (T > Tf) where longer-range effects turn out to be important. Kasuya and Kondo [25], when studying the anomalies of the electrical resistivity in antiferromagnetic materials that exhibit orderdisorder second-order phase transitions, have found that in the critical region λ is equal to the critical index of the specific heat magnetization α, which has a value close to zero. On the other hand, for τ outside the critical region (at T > T f ), Alexander et al. [19] predict that the value of exponent λ should be 1/2. Analyses of experimental E G vs. T data for various diluted magnetic semiconductors showing antiferromagnetic behaviour [22-24] confirm the predictions of both these models.

In order to obtain a theoretical expression for ∆E GM that can be compared with the data, Eq. (1) should be integrated.

ΔEGM(τ)=-Aτ1-λ/(1-λ)+c (2)

then, if the value of ∆E GM at T f is expressed as ∆E f , the constant of integration C in Eq. (2) can be removed and this equation can be written as:

ΔEGM(τ)-ΔEf=-Aτ1-λ/(1-λ) (3)

-ln[ΔEGM(τ)-ΔEf]=(1-λ)lnτ+ln[A/(1-λ)]. (4)

Thus, values of the magnetic shift ∆E GM determined from the experimental data can be analyzed in terms of Eq. (4). To calculate ∆E GM it is necessary to know the values of E G in the absence of any magnetic effects. These can be determined by extrapolating to low temperatures the E G (T) curves at high temperatures where the magnetic effects are unimportant. This may be done by using the expression for the E G vs T variation based on Bose-Einstein phonon model, given by [26]

EG(T)=EB-ab{1+2[exp(Φ/T)-1]-1}, (5)

where E B α b = E G (0) is the value of the energy gap at 0 K, and Φ an average phonon related to the Debye temperature.

The difference between the measured curves and those extrapolated by using Eq. (5) is thus attributed to the magnetic contribution∆E GM . Hence, the data shown in Fig. 3 have to be fitted to Eq. (5) in the high-temperature region (T ≥ 150 K). Thus, in order to reduce the number of the free parameters to be determined from the data, i.e. E B , α B , and Φ, to only two, we have taken Φ = 200 K, which is the Debye temperature for MnIn2Se4 estimated from its melting point (T M ≈ 1178 K) by using the classical Lindemann expression, which relates Φ to the melting point T M through the expression Φ ∼ TM 1/2 [27]. Hence, using Φ ≈ 200 K as a fixed parameter, Eq. (5) was fitted to the experimental E GI (T), EGDI (T), and EGDII (T) data of Fig. 3, in the temperature range from 150 to 300 K. Values of E G (0) and a B , obtained from this fit are also given in Table II. Theoretical extrapolated E GIE (T), EGDEI (T), and EGDEII (T) curves, obtained from Eq. (6), are also shown in Fig. 3 by broken lines. The magnetic shift can thus be obtained as ∆E GM (T) = E G (T)−E GE (T), for these three curves. The value of ∆E G at T f ≈ 3.6 K for each of the curves in Fig. 3, as obtained directly from these curves, was 34, 39 and 33 meV, for E GI (T), EGDI(T), and EGDII(T), respectively.

TABLE II Values of parameters E G (0) and a B , obtained from the fit of E GI (T), EGDI(T) and EGDII(T) data of Fig. 3 to Eq. (5) using Φ ≈ 200 K as a fixed parameter, and λ c , λ nc , A c and A nc , obtained from the fit of ∆E GM (τ) data to Eq. (4) in the critical and noncritical regions. 

Optical
transition
E G (0)
(eV)
αb
(eV)
λc λnc A c
(meV)
A nc
(meV)
EGI 1:516 ± 0:001 0:044 ± 0:06 0.02 0.46 1.3 2.0
EGDI 1:775 ± 0:001 0:031 ± 0:001 0.08 0.49 1.6 2.7
EGDII 1:925 ± 0:005 0:035 ± 0:001 0.09 0.52 1.4 2.5

The temperature variation of ΔEGM for the three band-to-band transitions in the form of a plot of ln(ΔEGM-ΔEG) vs. ln τ, in the temperature range where the magnetic contribution to EG vs. T is appreciable, is shown a broken line, in Figs. 4, 5a, and 5b, for EGI, EGDI, and EGDII, respectively. It can be seen that in the low τ region all the curves show a linear behaviour, in good agreement with Eq. (4). The slope (1-λ) in this temperature region is found to be 0.98, 0.92, and 0.91 (±0.02), for E GI , EGDI, and EGDII curves. These give λ c ≈ 0.02, 0.08 and 0.09 for these curves, respectively. This is also shown in Table II. These values are in good agreement with theoretical analysis that predicts λ ≈ 0 in the critical region [25]. In addition, in the high-τ region, the slope of the line is found to be 0.54, 0.51, and 0.48 (±0.02). This gives λ nc ≈ 0.46, 0.49, and 0.52 (±0.02) for these curves, respectively, which is also in agreement with the theory that predicts λ nc ≈ 0.5 in this region [19].

FIGURE 4 Temperature variation of ∆E GM for the indirect bandto-band transition E GI observed in MnIn2Se4 in the form of a plot of ln(∆E GM − ∆E f ) vs. lnτ, in the temperature range where the magnetic contribution to E G vs. T is appreciable (T < 150 K) (broken line). The slopes in both low and high-τ regions are indicated by bold continuous straight lines. 

FIGURE 5 Temperature variation of ∆E GM for the indirect band-to-band transition GGDI and GGDII observed in MnIn2Se4 in the form of a plot of ln(∆E GM − ∆E f ) vs. lnτ, in the temperature range where the magnetic contribution to E G vs. T is appreciable (T < 150 K) (broken line). The slopes in both low and high-t regions are indicated by bold continuous straight lines  

However, it is also observed that an intermediate temperature region exists where λ which was obtained from the slope of the broken lines at each temperature in Fig. 5, gradually increases between these two values. This is shown in Fig. 6 where the variation of λ with lnτ is plotted for the EGDI(T) optical transition. Values of λ for two alloys of Cd x Zn y Mn z Te pseudo-ternary system, a material for which a transition from antiferromagnetic-to-spin glass phase has been observed [34], are also shown in this figure for comparison.

FIGURE 6 Variation of λ¸ obtained from the slope of broken lines in Fig. 5, with ln τ in the critical and noncritical regions for the direct transition EGDI in MnIn2Se4 (black line). Values of λ¸ for the lowest band gap of Cd0:3Zn0:3Mn0:4Te (blue line) and Cd0:25Zn0:25Mn0:5Te (red line) CdxZnyMnzTe pseudo-ternary alloys, taken from [22] are also shown. 

The existence of this intermediate region, which does not appear in ideal antiferromagnetic compounds where an abrupt change when passing from critical to noncritical regions is observed [22-24], indicates that in the present case a gradual change is occurs when MnIn2Se4 pass from the spinglass phase to the paramagnetic one. This fact reveals that it is not a true phase transition in the thermodynamic sense but a gradual freezing of the magnetic moments.

Thus, it is found that in MnIn2Se4 below about 20 K, short-range spin correlations give the main contribution to ΔEGM vs. T, while longer-range effects becomes dominant at temperatures between about 70 to 150 K. In the intermediate temperature region, where λ¸ varies from about 0.02 to 0.40, both these two spin correlation effects should be taken into account in the analysis of magneto-optical properties of this compound.

Values of the parameters A in the critical and non-critical regions, obtained from the fit are also shown in Table II. These values of A c and A nc , are in good agreement with those reported for antiferromagnetic semiconductors [22,24] that vary from 0.2 to 28.0 and 4.8 to 18.0 meV, respectively.

5. Conclusion

The temperature dependence of the indirect and direct fundamental band gaps of single crystals on the layered compound MnIn2Se4, that crystallize in the rhombohedral defect structure with space group R3-mD3d5, were studied by optical absorption spectra. The data were analyzed in terms of the theoretical models by Alexander et al. [19] and Kasuya and Kondo [25] that takes into account the magnetic contribution to the shift of E G with T. It is suggested that short-range effect spin correlations dominate the contribution to the shift below about 20 K, in the critical region of the paramagnetic to spin-glass phase transition, whereas longer-range effect spin correlations control the contribution to this shift in the temperature range from about 70 to 150 K, in the noncritical region. An intermediate temperature region, where a gradual freezing of the magnetic moments occurs, was also observed. This probably indicates that, from a thermodynamic point of view, the spin-glass phenomenon is a gradual freezing of the magnetic moments instead of a true phase second order transition.

Acknowledgments

This work was supported by the Consejo de Desarrollo Científico, Humanístico, Tecnologico y de las Artes (CD-CHTA) de la Universidad de Los Andes, Merida-Venezuela, under Project No. C-1658-09-65-A.

References

1. H. Neumann, C. Bellabarba, A. Khan, V. Riede, Cryst. Res. Technol 21 (1986) K21. [ Links ]

2. K. J. Range, and Hubner, Z. Naturf. B 30 (1975) 145. [ Links ]

3. G. E. Delgado, C. Chacon, J. M. Delgado, V. Sagredo, Phys. Status. Solidi A 134 (1992) 61. [ Links ]

4. M. Cannas, L. Garbato, A. G. Lehmann, N. Lampis, F. Ledda, Cryst. Res. Technol . 33 (1998) 417. [ Links ]

5. V. Sagredo, E. ter Haar, G. Attolini, Physica B 320 (2002) 407. [ Links ]

6. R. Tovar et al ., Mater. Res. Bull. 37 (2002) 1011. [ Links ]

7. N. N. Niftiev, M. A. Alidzhanov, O. B. Tagiev, F. M. Mamedov, M. B. Muradov, Semiconductors, 38 (2004) 531-532. [ Links ]

8. V. Sagredo, M. C. Moron, L. Betancourt, G. E. Delgado, J. Magn. Magn. Mater. 312 (2007) 294. [ Links ]

9. F. J. Manjon et al ., Phys. Status Solidi. B 244 (2007) 229. [ Links ]

10. M. C. Menard et al ., Chem. Mater. 23 (2011) 3086. [ Links ]

11. B. G. Tagiev, et al., Semiconductors, 46 (2012) 705. [ Links ]

12. S. H. You, K. J. Hong, T. S. Jeong, K. Y. Lim, C. J. Youn, J. Cryst. Growth 404 (2014) 116. [ Links ]

13 M. Quintero, M. Morocoima, A. Rivero, P. Bocaranda, J. C. Woolley, J. Phys. Chem. Sol. 58 (1997) 491. [ Links ]

14. J. C. M. Ochoa et al, J. Magn. Magn. Mater, 272-276 (2004) 1308. [ Links ]

15. J. C. Mantilla et al., J. Phys. Condens. Mat. 17 (2005) 2755. [ Links ]

16. G. Doll, M. Ch. Lux-Steiner, Ch. Klock, J.R. Baumann, E. Bucher, J. Cryst. Growth 104 (1990) 593. [ Links ]

17. K.-J. Range, U. Klement, G. Doll, E. Bucher, J. R. Baumann, Z. Naturf. 46b (1991) 1122. [ Links ]

18. C. Rincon, T. E. Torrres, V. Sagredo, S. J. Jiménez-Sandoval, E. Mares-Jacinto, Physica B 477 (2015) 123. [ Links ]

19. S. Alexander, J. S. Helman, J. Balberg, Phys. Rev. B 13 (1976) 304. [ Links ]

20. R. B. Bylsma, W. M. Becker, J. Kossout, U. Debska, D. Yoder-Short , Phys. Rev B 33 (1986) 8207. [ Links ]

21. J. I. Pankove, Optical Properties of Semiconductors. (Dover Pub., Inc. New York, 1971). [ Links ]

22. T. Donofrio, G. Lamarche, J. C. Woolley, J. Appl. Phys. 57 (1985) 1932. [ Links ]

23. M. Quintero, B. D. Marks, J. C. Woolley, J. Appl. Phys. 66 (1989) 2402. [ Links ]

24. X. L. Chen, A. M. Lamarche, G. Lamarche, J. C. Woolley, J. Phys. Condens. Mat. 5 (1993) 7143. [ Links ]

25. T. Kasuya and A. Kondo, Solid State Commun. 14 (1974) 249. [ Links ]

26. L. Viña, S. Logothetidis, M. Cardona, Phys. Rev. B 30 (1984)1979. [ Links ]

27. L. E. Shelimova and S. K. Plachkova, Phys. Status Sol. A 104 (1987) 679. [ Links ]

28. T. Donofrio, G. Lamarche, and J. C. Woolley, J. Appl. Phys. 57 (1985) 1937-1940. [ Links ]

Received: August 01, 2018; Accepted: September 21, 2018

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License