JPCS 3
JPCS 3
JPCS 3
A R T I C L E I N F O A B S T R A C T
Keywords: Several quaternary nanocomposites, doped with two transition metal oxides, of the generalised composition
Quaternary nanocomposites formula 0.4ZnO–0.1P2O5–0.5[xV2O5–(1-x) MoO3] were prepared via the melt quenching process. We investi
X-ray diffraction gated the microstructural and conductivity properties of two ternary systems and four quaternary systems to
Optical bandgap
determine the better semiconductor nanocomposite system for wider application purposes. The X-ray diffraction
Dc and ac conductivity
CBH and NSPT model
results showed the presence of superposed nanocrystallites within amorphous networks or matrices. Analysis of
ultraviolet–visible absorption spectra showed that optical bandgap energy (Eopt) values varied with V2O5 con
centration, and there was an inverse relationship between Eopt and average nanocrystallite size. The responsible
mechanisms of ac conductivity were examined using the model of Jonscher’s universal power-law and
Almond–West formalism, and it was found that ac conductivity increased as temperature rose, exhibiting sem
iconducting features. All the ternary (two examples) and quaternary (four examples) nanocomposite systems
established non-linearity in dc conductivity, caused by the small polaron hopping process with dissimilar acti
vation energies at different temperature regions. The power-law exponent (s) of Jonscher’s universal power-law
revealed that the ac conductivity mechanism could be described by correlated barrier-hopping (CBH) or non-
overlapping small polaron tunnelling (NSPT) models for different samples due to major structural alterations.
We applied modified CBH and NSPT models (as s > 0.8) to get realistic values of different fitting parameters, and
in this process, we also obtained the values of ideal or hypothetical glass transition temperature from the
theoretical models. The different activation energies associated with ac conductivity and for the small polaron
migration process diminished with an increment in conductivity. Scaled ac conductivity spectra demonstrated
that the conductivity relaxation process relied upon the composite structural features and did not depend on
temperature. These materials can be used for such applications as gas sensors, even at higher temperatures due to
their semiconducting nature and the different valence states of transition metal ions.
* Corresponding author.
E-mail address: anindyasundardas03@gmail.com (A.S. Das).
https://doi.org/10.1016/j.jpcs.2020.109505
Received 13 November 2019; Received in revised form 10 April 2020; Accepted 11 April 2020
Available online 21 April 2020
0022-3697/© 2020 Elsevier Ltd. All rights reserved.
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
reported [1–3]. The MoO3-doped phosphate glasses have electro-optical compare variations in the structural and conductivity properties among
applications caused by electrochromism properties and higher ionic all the as-quenched ternary and quaternary glass compositions to show
conductivity [4,5]. The glass formation and structural features of which system has better conductivity properties for diverse types of
ZnO–P2O5–MoO3 glasses have been reported [6]. Furthermore, the in applications.
teractions of molybdenum (Mo) in phosphate glasses are advantageous
to manage radioactive wastes [6]. The electrical transport mechanisms 2. Experiment section
in ZnO–MoO3–P2O5 glasses have been described by the small polaron
hopping (SPH) process [7]. Three bridging oxygen atoms and one 2.1. Materials and sample preparation
double-bonded oxygen atom make tetrahedral units of PO4, which is the
building block P2O5-based glasses [8]. In the absence of cations, the Four quaternary and two ternary glass nanocomposite samples of the
development of different phosphate anions takes place due to the linking common formula 0.4ZnO–0.1P2O5–0.5 [xV2O5–(1-x) MoO3], with the
of covalent bridging oxygen ions with PO4 tetrahedral units. The pres variation of x in steps of 0.2, were synthesised via melt quench tech
ence of different valence states and also different co-ordinations, which nique. Suitable volumes of high-purity chemicals ZnO (Loba Chemie:
form corner-sharing VO4/VO5 polyhedral units, makes V2O5 a good 99% purity), V2O5 (Loba Chemie: 98% purity), MoO3 (Loba Chemie:
glass network former, and V2O5 also performs as a better glass network 99% purity), and P2O5 (Merck Specialities Private Limited: 98% purity)
modifier when combined with the phosphate-based glass nano were carefully mixed in ratios according to the composition formula of
composites. Transition metal ions (TMIs) are widely used as glass the individual system. Then the mixtures were heated to the melting
composite materials since their outer d-electron orbital functions have a point of the compositions in an electric furnace in atmospheric air
wide radial distribution and owing to their profoundly sensitive condition within a temperature range of 790� C-945� C depending on the
response to the variations in neighbouring cations [9–11]. As the longer molar ratio of the composition. Ternary bulk glassy systems, i.e. samples
covalent chains shorten and local structure reforms, the phosphate with x ¼ 0.0 and 1.0, melted at 945� C and 810� C, respectively. Whereas,
network depolymerises, evolving efficiently to a refined vanadium the quaternary bulk glass compositions (x ¼ 0.2 to 0.8) melted within
phosphate network with isolated PO4 tetrahedral units [11–13]. The the temperature range from 885� C to 785� C with an increase in V2O5
literature shows that the inclusion of TMIs enhances the chemical content. We continuously monitored the melting process in every step
durability and stability of glassy systems by substituting bridging P–O–P and, after the melting of the oxide materials, we checked the fluidity.
bonds in phosphate glasses with additional moistness resistant P–O–TM Then we kept the melt at the melting temperature for 20 min to
bonds [11–13]. homogenise and equilibrate the melt and to achieve better fluidity for
The electrical transport properties of TMI-doped phosphate glasses pouring. Then the melt was transferred onto highly polished aluminium
can be enhanced as mobile charge carriers or small polarons start the plate, and quickly quenched by pressing the melt in a technique using
hopping process from lower to higher valence states. Glasses doped with another highly polished aluminium plate at room temperature. Conse
two TMIs may attain super-exchange of electrons to a great extent be quently, small glassy flakes of thickness ~0.2–0.5 mm were obtained.
tween the redox ion pairs, known as hopping, leading to the non-linear Then as-quenched glassy flakes were ground into a fine pure powder for
properties of activation energy [12,13]. Such non-linearity occurs owing XRD and UV–Vis absorption measurements.
to the mixed transition metal-ion effect (MTIE) [14–16]. In glassy sys
tems with mixed TMIs, the electrons or small polarons also hop between 2.2. Volumetric mass density measurement
diverse TMI sites, for example, Mn2þ to Fe3þ and Fe2þ to Mn3þ in
Fe2O3–MnO–TeO2 glasses [12]. After doping with two TMOs, the con Using the renowned Archimedes principle, volumetric mass density
ductivity of the Fe2O3–P2O5–V2O5 glassy system reduces owing to the (ρ) values were measured by immersing samples in acetone because
Anderson localisation effect [13]. After adding MoO3 to the glassy sys acetone is non-reactive with the samples being examined. The ρ is
tem, Mo ions exist in the stable Mo5þ (donor level) and Mo6þ (acceptor considered a significant parameter, and an influential tool to study the
level) valence states. The glass matrix becomes more compact as Mo6þ structural modifications occurring within the polycrystalline or amor
ions reside in the network-forming positions with MoO4 tetrahedral phous material. The value of ρ of all bulk glass samples was measured
units. Contrariwise, in Mo5þO3 complexes, Mo5þ ions are positioned in using the following relationship [21]:
place-modifying positions [17]. In general, the coexistence of V4þ � �
(donor level) and V5þ (acceptor level) states of V ions are observed in ρ¼
Wair
� ρacetone (1)
vanadate-based glasses. The octahedral V4þO6 and VO4 tetrahedral sites Wair Wacetone
are formed due to the presence of V4þ ions, while V5þ ions create VO5
where ρacetone is the density of acetone and Wacetone and Wair are the
pyramidal units [18,19]. The distribution of charge carrier jump dis
weights of the bulk glassy sample immersed in acetone and suspended in
tances or the energy differences of electrons positioned at the unrelated
air, respectively. Average ρ was determined after measuring ρ of each
potential wells due to the presence of charged impurity atoms, in
sample quite a few times. The value of ρ provides some basic structural
fluences the electrical conduction process [20]. Additionally, the vari
features such as molar volume (VM) and the ionic packing ratio. The VM
ation in the electrical transport mechanism and the activation energy
of each composite was computed from the values of total molecular
necessary for the conductivity mechanism is instigated by the
weight (M) of the composition and ρ according to the following rela
thermally-activated structural vibrations or due to differences in oxygen
tionship [22]:
coordination [20].
The purpose of this work is initially to synthesise four quaternary and X xi � MWi
VM ¼ (2)
two ternary glassy systems of composition 0.4ZnO–0.1P2O5–0.5 ρ
[xV2O5–(1-x) MoO3] by using the melt quench technique. Here, we
study the electrical transport mechanism of the quaternary systems where MWi is the ith oxide molecular weight, xi is the molar fraction of
using the property of MTIE due to the presence of two TMOs: V2O5 and the ith oxide, and ρ is the determined value of the average density of each
MoO3. Furthermore, one objective of this work is to explore the micro glass composition.
structural and optical absorption properties of all nanocomposite sys
tems using X-ray diffraction (XRD) and ultraviolet–visible (UV–Vis) 2.3. Characterisation techniques
absorption techniques. A second objective is to investigate the probable
dc and ac electrical conduction mechanisms and to determine the For structural characterisation, an X-ray diffractometer (Rigaku
different activation energies required for the processes. Finally, we TTRAX-III) was utilised, and the chosen wavelength of CuKα X-ray
2
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
radiation was 1.5418 Å. The scanning speed was 4� per minute in steps 3. Results and discussion
of 0.02� , while Bragg’s angle (2θ) ranged within 10–80� . According to
the Gaussian distributions [23], the microstructural properties of the 3.1. Analysis of X-ray diffractograms
samples were inspected using the line profile analysis (LPA) method.
A PerkinElmer Lambda-750 UV–Vis spectrophotometer was used in Fig. 1a presents the X-ray diffractograms of ternary (two examples)
absorbance mode to obtain the optical absorption spectra of the nano and quaternary (four examples) finely powdered nanocomposite sam
composites at 300 K in a wide range of wavelengths (190–700 nm). The ples. All samples clearly exhibit a definite amount of amorphous phase,
powdered samples were dissolved in ethyl alcohol to form colloidal which is the feature of short-range order. The amount of amorphous
suspensions for the absorption procedure and homogenised using an phase present within the samples can be calculated from the relative
ultrasonic sonicator (Takashi SK-500F). percent crystallinity of each nanophase. Relative percent crystallinity
can be determined by deconvolution of the amorphous state and crys
2.4. Electrical conductivity measurement tallinity or by the ratio of the crystalline peak intensity to the sum of the
amorphous and crystalline peak intensities, mathematically expressed
For electrical measurements, highly conductive and high- by the following relationship [25]:
temperature resistant silver paste was placed on either side of the bulk
glassy flakes to function as the electrodes. The sample temperature Relative percent crystallinity ¼ Icrystalline/(Icrystalline þ Iamorphous) (4)
within the furnace was recorded by a thermocouple, located nearby the Furthermore, according to Ruland [26], the ratio of the total area
sample surface. The thermocouple was connected to a temperature under any crystalline peak to the total area under the total diffraction
control unit to measure the operating temperature of samples with �5% pattern (crystalline plus amorphous peaks) provides the value of relative
error. Alternating current (ac) electrical measurements were performed percent crystallinity as expressed below [27,28]:
using an automatic LCR Hi-Tester (Model No: 3532–50; Hioki) in the
frequency window of 42 to 5 MHz. The values of parallel capacitance Relative percent crystallinity ð%Þ ¼
Area under the crystalline peak
(CP), dielectric loss tangent (tan δ), and conductance (G) were recorded Total area under the diffraction curve
and the experimental value of total electrical conductivity (σtotal) was ðcrystalline þ amorphousÞ
calculated using the following relationship [24]: � 100
(5)
Conductance ðG in siemensÞ � thickness ðt in cmÞ
σtotal ¼ � (3)
Area A in cm2 The values of relative percent crystallinity of each nanophase are
shown in Table 1, illustrating the development of a few nanocrystallites
The accuracy of the temperature measurements was �3 K, and within the amorphous structure of the glass network or matrices.
conductance measurements had a maximum error of �2 μS. The relative percent crystallinity values are in agreement with
Fig. 1a, which also confirms the superposition of a few nanocrystallites
over the amorphous glass network irrespective of the composition
molarity, which is characteristic of long-range order. The nanophases of
each distinct diffraction angle peaks have been identified and indexed
through laborious evaluation from existing literature data. The
Fig. 1. (a) X-ray diffractograms of all bulk nanocomposites with variations in x and (b) the composition-dependent average crystallite size of all bulk samples with
error bars, and the dotted lines are guidelines.
3
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
Table 1
Values of diffraction peak angle (2θ), Corrected Full width at half maxima (FWHM), nanocrystallite size using Scherrer equation (dC), nanocrystallite phases and
corresponding Miller indices (h, k, and l), and percentage of crystallinity of 0.4ZnO–0.1P2O5–0.5 [xV2O5–(1-x) MoO3] glass nanocomposites. Ranges of errors are
indicated within the parentheses.
x (mol %) 2θ (degree) FWHM (βcorrect) dC (nm) Nanocrystallite phase h k l Percentage
(�0.03) (�0.0035) (�1.05) Crystallinity (%) (�0.25)
development of superimposed ZnP4O11 [29], ZnMoO4 [30], MoP3O9 nanocomposites (x ¼ 0.2 to 0.8) are seen in Fig. 1a and Table 1.
[31], Zn2P2O7 [32], Zn2.5MoVO8 (ICDD-50-1894), MoV2O8 [33], Although samples of x ¼ 0.0 (ZnO–P2O5–MoO3 ternary system) and x ¼
Zn3MoV2O11 [34], V4[P2O7]3 [35], VP3O9 [36], Zn2V2O7 [37], VOPO4 1.0 (ZnO–P2O5–V2O5 ternary system) reveal the presence of a definite
[38], and V2O5 [39] nanophases within the amorphous glassy networks amount of amorphous phase (Table 1), the number of developed nano
or matrices have been predicted and confirmed. The LPA method is used crystallites exceeds those in the quaternary systems except for x ¼ 0.8,
to analyse the microstructural properties following Gaussian distribu and these samples are considered polycrystalline due to the absence of a
tions [40]. The broadening effect of the instrument is corrected by broad hump. For quaternary systems with values of x ¼ 0.2, 0.4, and 0.6,
considering the βinstrumental value of standard crystalline silicon. The a broad hump is visible around 2θ ¼ 27� , which confirms the total
corrected full width at half maximum (FWHM) (βcorrect) values of the amorphous structure of the samples with very few nanocrystallites. The
diffraction peaks are the result of subtraction of the observed broad evolution of total amorphous structure with the presence of a broad
ening, βobserved, and instrumental broadening, βinstrumental, effects. Thus, hump appears for x ¼ 0.2 and then increases progressively with rising
the value of the Gaussian distribution function, βcorrect, used in the cal V2O5 concentration up to x ¼ 0.6; for x ¼ 0.8, the broad hump disap
culations can be expressed as follows [41]: pears due to major structural variations, thus representing a poly
�0:5 crystalline sample. The presence of superposed nanocrystallites inside
βcorrect ¼ β2observed β2instrumental (6) the amorphous glass networks or matrices alters the structural features
Using the well-known Scherrer equation, the size of each nano and also electrical conductivity.
crystallite (dC) is determined by the following equation [41,42]: Fig. 1b illustrates average size of nanocrystallites with variation in
V2O5 concentration (x). For x ¼ 0.0, the average size of nanocrystallites
dC ¼
0:89 λ
(7) is highest and followed by a decrement in the average size of nano
βcorrect cos θ crystallites with inclusion of V2O5 up to x ¼ 0.6 due to formation of the
total amorphous nature of the corresponding composites. Average dC
where λ is Cu-Kα X-ray radiation wavelength (1.54184 Å). The super then increases as samples (x ¼ 0.8 and 1.0) develop more crystallinity
posed nanophases with their peak θ, βcorrect, values of the equivalent than the pure amorphous state.
Miller index, and each dC are presented in Table 1. The distinct varia
tions of the XRD patterns of ternary (x ¼ 0.0 and 1.0) and quaternary
4
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
3.2. Analysis of ρ and VM separation <dMo–Mo> [22,45]. For that, the volume of 1 mol of V (VVM)
and Mo (VMoM ) within the given structure can be computed using the
Fig. 2 presents the differences in ρ (Eq. (1)) and VM (Eq. (2)) values following expressions [22,45]:
with variations in V2O5 concentration (x) and shows the usual inverse
VM VM
relation of VM and ρ [22]. According to the composition formula, the VVM ¼ and VMo
M ¼ (8)
2ð1 XV Þ ð1 X Mo Þ
molecular weight of each sample increases with increasing V2O5 con
ion ion
centration (x) (Table 2) because of the higher molecular weight of V2O5 where XV-ion and XMo-ion are the molar fractions of V2O5 and MoO3,
(181.88 g mol 1) compared to MoO3 (143.94 g mol 1). Thus, incorpo respectively. The <dV–V> and <dMo–Mo> (Table 2) are evaluated using
ration of V2O5 in the glass matrix acts as a modifier, causing an incre the following relationships [22,45]:
ment in non-bridging oxygen ions (NBOs), which escalates the breaking
� �1=3 � �1=3
of covalent bonds within the network. For x ¼ 0.0 (ternary system of VVM VMo
composition ZnO–P2O5–MoO3), the value of ρ is higher than x ¼ 1.0 < dV V > ¼ and < dMo Mo > ¼ M
(9)
NA NA
(ternary system of composition ZnO–P2O5–V2O5) due to MoO6 poly
hedra group formation [43]. where NA is Avogadro’s number (6.0228 � 1023 mol 1). The <dV–V>
The higher concentration of MoO3 typically reorganises the structure value increases and <dMo–Mo> value reduces progressively with incre
of the glass network or matrix by creating additional oxygen ions, ment in V2O5 concentration (x) (Table 2). Therefore, due to the cumu
resulting in the changeover of MoO4 to Mo2O7 or MoO6 structural units lative effect of inclusion of V2O5 strongly expanding the quaternary glass
[43]. During the development of quaternary glassy systems with the network, this directly supports the decline in ρ and the increasing trend
inclusion of V2O5, the ρ value clearly declines up to x ¼ 0.8 (Fig. 2), of VM for as-quenched nanocomposites. Such nature of ρ and VM may
revealing a loosely packed amorphous glass network, and ρ increases arise from the formation of excess NBO bonds, thus expanding the glass
again for x ¼ 1.0 (ZnO–P2O5–V2O5, ternary system). The decrement of ρ network structure [22,45]. Using the concept of chemical bond
value with V2O5 concentration rise is caused by a number of factors, approach in explaining VM and glass ρ behaviour, the oxygen molar
firstly, introduced V ions reside in interstitial sites, thus reducing the volume (VO) and oxygen-packing density (OPD) of each nanocomposite
vacancies inside the glass matrix [44]. Furthermore, the conversion of are computed using the following ratios [45]:
six- to four-fold coordination states of V and Mo ions also reduces the ρ � ρ �
VM
value [44]. For x ¼ 1.0 (ZnO–P2O5–V2O5, ternary system), the ρ value VO ¼ P and OPD ¼ 1000 C (10)
x O M
increases considerably due to a slight decrease in volume due to the n n n
5
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
Table 2
Molecular weight (M) of each composition, vanadium volume (VVM), average vanadium–vanadium separation <dV–V>, molybdenum volume (VMo M ), average molyb
denum–molybdenum separation <dMo–Mo>, molar volume of oxygen (VO), oxygen-packing density (OPD), V ion concentration (NV-ions), and Mo ion concentration
(NMo-ions) of all the glass samples for different values of x. For the calculation, we considered the valence number of V as five and Mo as six. Ranges of errors are
indicated within the parentheses.
x (mol Molecular Weight VVM (g <dV–V> VMo
M (g <dMo–Mo> VO (cm3 OPD (mol NV-ions ( � 1022 ions/ NMo-ions ( � 1022 ions/
%) (M) cm–3) (Å) cm–3) (Å) mol 1) l 1) cm3) cm3)
(g mol 1) (�0.20) (�0.04) (�0.20) (�0.04) (�0.20) (�0.05) (�0.25) (�0.25)
(�0.02)
6
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
of a number of defect states and disorder forms localised or extended 0.0, 0.4, 0.8, and 1.0 are shown in Fig. 6a–d, respectively. Qualitatively
states inside the optical bandgap or mobility gap. In the non-linearity all the other samples also exhibit similar σtotal behaviour. The σtotal of all
plot of the absorption coefficient (ln α), near the optical bandgap bulk nanocomposites increases non-linearly with equal increment in
edge, an exponentially decaying band tail is detected that is expressed as temperature, confirming semiconducting behaviour. Individually iso
Urbach energy (EU) or band tail energy. The EU value can be estimated therms illustrate two typical features, i.e. a plateau or a flat dc con
using the following Urbach empirical rule [22,50]: ductivity (σdc) region, which is nearly frequency independent owing to
� � the effect of space charge polarisation and the conduction process is the
hν
α ¼ α0 exp (14) consequence of small polarons or charge carrier hopping mechanism.
EU
The second feature is the dispersion that varies almost linearly at the
higher frequencies. The range of frequency-independent conductivity
where α0 is a constant term when (hν/EU) → 0. Generally, EU refers to the
increases to a higher frequency with temperature rise and dispersion
amount of defect, extended, or localised states present in the vicinity of
regions nearly coincide with each other at higher frequencies due to the
Eopt [51]. The value of the reciprocal of the slope of the linear region of
mixed electronic and ionic conductivity nature of the systems (Fig. 6a
the plot of ln (α) versus hν provides the EU value, as demonstrated by
and b). In contrast, for x ¼ 0.8 and 1.0, the range of frequency inde
solid lines in Fig. 5a. The composition-dependent values of EU with the
pendent and dispersive regions is almost the same for the entire tem
variation of V2O5 molar ratio (x) are displayed in Fig. 5b. The EU value
perature range (Fig. 6c and d), i.e. the crossover frequency does not alter
drops up to x ¼ 0.6 and then rises again for x ¼ 0.8 and 1.0. Thus, it can
much, which confirms that conductivity is electronic or polaronic in
be concluded that for the present nanocomposite systems, the EU values
nature. The dispersal behaviour of frequency-dependent conductivity
decline as the degree of amorphous state rises. Comparative study of Eopt
originates either due to tunnelling or hopping mechanisms of small
and EU, as shown in the inset of Fig. 5b, reveals the inverse relationship
polarons or mobile charge carriers through or among the defect or
between Eopt and EU, i.e. as Eopt increases, EU falls, and dC also declines
localised states existing within the range of Eopt or owing to the ionic
(Fig. 1a). The XRD data show the existence of low crystallinity or
motion of mobile ions inside the glassy matrices [52,53]. The unified
amorphous nature of all the nanocomposites (Fig. 1a), which is
forward and backward motion of charge carriers or small polarons
confirmed by the optical absorbance spectra.
during hopping or tunnelling mechanisms causes the dispersive nature
of σtotal spectra. During the hopping process, small polarons or charge
3.4. Study of electrical conductivity carriers generally hop or jump from one localised or defect state to
another, which requires less energy than Eopt, thus resulting in the
Characteristic plots of total ac conductivity (σtotal) isotherms for x ¼ reduction of activation energy with a rise in frequency. All σtotal plots
(Fig. 6) present characteristics of the electronic or small polaronic
transport mechanism between sites with randomly separated space and
energy. The σtotal plots disclose that σtotal is the sum of both ac con
ductivity (σac (ω)), which is frequency-dependent but weakly
temperature-dependent and the temperature-dependent and
frequency-independent σdc. The σtotal can be examined by the following
Almond–West formalism [52,53]:
� �
ω n
σtotal ¼ σdc 1 þ (15)
ωH
7
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
Fig. 6. (a) Temperature- and frequency-dependent total ac conductivity (σtotal) spectra of the glass compositions for (a) x ¼ 0.0, (b) x ¼ 0.4, (c) x ¼ 0.8, and (d) x ¼
1.0. Black solid lines indicate best non-linear curve fitting data using Almond–West formalism.
Table 3
The average values of frequency exponent (n), composition-dependent values of ac conductivity activation energy (Eac), dc conductivity activation energy (Edc), the
radius of small polarons and density of defect states at the high-temperature region with errors within the parentheses.
x (mol %) Average value of n Eac (eV) Dc conductivity activation energy (Edc) (eV) Small polaron radius Density of defect states at high temperature
(�0.05) N(EF) ( � 1021 eV 1cm 3)
Elow (eV) Ehigh (eV) RP-V (Å) RP-Mo (Å) N(EF) of V ions N(EF) of Mo ions
(�0.02) (�0.02) (�0.02) (�0.02)
increment in V2O5 concentration, but a drop in the rate of crossover or independent of frequency.
hopping frequency is also detected as observed for x ¼ 0.2 (Fig. 7). Hence, it can be asserted that the reciprocal temperature response of
Using the fitting of Almond–West formalism with the experimental σac directs that the conduction process is due to space charge conduc
data, we obtain σdc and σac is calculated as follows: σac ¼ σtotal – σdc. The tion. The calculation of ac conductivity activation energy (Eac) for each
σac dependency with the variation of frequency of reciprocal tempera nanocomposite at various frequencies is carried out using the following
ture for x ¼ 0.4 reveals that σac increases consistently with an increment empirical relationship [55,56]:
in frequency and temperature (Fig. 8a). The σac usually depends on the � �
Eac
structural symmetry, grain and crystallite size, and distribution of grains σac ¼ σ0 exp (16)
KB T
within the glassy network [55]. In bulk polycrystalline or amorphous
materials, the σac mechanism is either due to the long-range hopping
where KB is Boltzmann constant (¼ 8.617 � 10 5 eV K 1) and σ0 is the
mechanism of charge carriers or small polarons or to the localised
value of ac conductivity at very high temperature when (Eac/KBT) → 0.
transportation resulting from oxygen vacancies [55]. The slopes of the
The Eac can be determined by the slope of the linear relationship, as
Arrhenius plot reduce with increasing frequency, and σac increases with
revealed by solid lines in Fig. 8a and b and Eq. (16). The estimated Eac
a rise in frequency, indicating the dispersive response of σac (Fig. 8a).
values are shown in Fig. 8a for different frequencies and in Table 3 for
This is the result of the possible release of space charges at higher fre
different composites of varying values of molar ratio (x). The Eac value
quencies, even at lower temperatures. At a higher temperature for
reduces as frequency rises (Fig. 8a), i.e. less energy is required to counter
different frequencies, σac does not differ much because at a higher
the oscillations of excited charge carriers at high frequencies due to
temperature the space charges release and recombine almost at the same
improved jumping or hopping of electrons between neighbouring
rate, therefore, the σac of the composite samples becomes almost
localised states [56]. Frequency-dependent Eac values show that for
8
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
(0.3 eV) [57]. A comparison of Eac for different compositions with the
variation in the molar ratio of V2O5 concentration (x) at a specific fre
quency (3 kHz) is presented in Fig. 8b. The Eac values of different
composites are presented in Table 3, indicating the inverse relationship
of conductivity and Eac. The Eac value for x ¼ 0.0 and 0.2 is nearby the
activation energy of the second ionisation of oxygen vacancies (0.7 eV;
Table 3) [57]. Relatively, Eac values for x ¼ 0.4–1.0 are in the range of
the first ionisation of oxygen vacancies (0.3 eV) [57]. With increasing
V2O5 concentration, Eac value falls as the inclusion of V2O5 induces an
increment in the number of glass network defects as the amorphous
phase of the samples increases (Fig. 1a) [55]. The diminishing nature of
Eac can also be ascribed to electrons trapped on the defects [55].
The variations of crossover or hopping frequency (ωH) of charge
carriers or small polarons with the reciprocal of temperature, which
obeys the Arrhenius equation, are shown in Fig. 9. The values of the
small polaron migration or crossover frequency activation energy (EH)
can be realised using the least-squares straight-line fit of the data as
shown by the solid lines in Fig. 9 using the following Arrhenius equation
[58]:
� �
Fig. 7. Total ac conductivity (σtotal) spectra of all nanocomposite samples with ωH ¼ ωe exp
EH
(17)
the variation of V2O5 molar ratio (x) at 653 K. KB T
frequency up to 1 MHz, Eac > 1 (Fig. 8a), which describes the σac
mechanism on account of long-range hopping of electrons or small po Fig. 9. Plot of the reciprocal temperature dependence of hopping or crossover
larons. Whereas, in the high frequency (�3 MHz) region, Eac value is frequency of all nanocomposites. Solid lines represent Arrhenius linear fits with
nearby the activation energy of the first ionisation of oxygen vacancies calculated activation energy values.
9
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
well, electrons are trapped, forming localised sites, and electrons or governed by the hopping process between similar TMI sites (between
small polarons can hop between different valence states of similar or V4þ and V5þ and between Mo5þ and Mo6þ) as well as unlike TMI sites,
dissimilar TMIs. such as between V5þ and Mo6þ and between V4þ and Mo5þ [14–16]. The
The σdc data, obtained from Almond–West formalism, strongly maximum conductivity for x ¼ 0.8 is due to the diminishing effect of
depend on the temperature; and the non-linearity of σdc is due to the Anderson localisation instigated by an increment in disorder states when
involvement of dissimilar conduction mechanisms [59,60]. The two TMIs are mixed, promoting electronic or polaronic motion, thus
temperature-dependent σdc plot is demonstrated in Fig. 10, and σdc can causing an enhanced σdc. Additionally, the formation of Zn3MoV2O11
be expressed by the Arrhenius equation [59,60]: nanocrystallites greatly raises σdc owing to a higher degree of MTIE.
� � � � Significant lattice distortion is a criterion for SPH to occur, and this
E high E low
σdc ¼ σ high exp þ σ low exp (18) process is usually a multi-phonon hopping process. In TMO-doped bulk
KB T KB T
glass, the SPH model describes the σdc using the following relationship
Two terms appear in Eq. (18) due to the two different σdc mecha [61,62]:
nisms in dissimilar temperature regions [60]. In Eq. (18), the first term � � � �
ν e2 W
defines σdc in the high-temperature region due to the band conduction σdc ¼ C ð1 CÞexpð 2αRÞexp (19)
KB T R KB T
process through the extended states, developed inside the mobility gap
as a result of the amorphous nature of the nanocomposites. The second
where W is activation energy, ν is optical phonon frequency, T is abso
term of Eq. (18) refers to the conduction process arising due to the
lute temperature (in Kelvin), and c is the redox ratio of TMIs [(reduced
hopping of small polarons or electrons via localised or extended or
TMI)/(total TMI)]. Here R implies the hopping distance concerning the
defect states at a lower temperature range or due to the tunnelling
available localised centres, and α represents the wave function of small
mechanism between unoccupied sites of the nearest neighbouring sites.
polarons at a localised centre decay as exp ( αR). The small polaron
The pre-exponential factors σhigh and σlow are strongly composition
radius (RP) can be predicted, knowing the value of the inter-nuclear
dependent [59,60].
distance of V and Mo ions using the following formula [22]:
However, values of activation energies associated with dc conduc � �
tivity at high (Ehigh) and low (Elow) temperature regions are obtained RV = RMo �π�1=3
RPðV=MoÞ ¼ (20)
from the gradient of the fitted lines indicated by the solid lines in Fig. 10. 2 6
The estimated values of Ehigh and Elow are given in Table 3. The Elow
values decrease with rising V2O5 concentration, and the lowest value of with a rise in V2O5 concentration, the estimated value of the small
Elow is for x ¼ 0.8, which reveals the highest conductivity (Fig. 10). The polaron radius of V ions (RP-V) reduces and the small polaron radius of
Ehigh value reduces as the atomic bonding disorder among the adjacent Mo ions (RP-Mo) increases (Table 3). According to the concept of small
sites increases, which enhances the density of tail and extended states polaron radius, σdc increases as the radii of small polarons reduces,
within the region of Eopt [59,60]. With the addition of V2O5 in the therefore, it can be concluded that V ions or V2O5 act as glass network
ternary system of the composition ZnO–P2O5–MoO3, the Ehigh value modifiers, and assist in enhancing conductivity for our present systems.
increases for x ¼ 0.2 due to substantial structural changes or configu The Fermi level defect states density [N(EF)] can be anticipated by
ration disorder, materialising because of the mixing of two TMIs [11] applying the assessed values of RMo and RV using the following expres
and the development of new crystallites over the amorphous matrix. sion [22]:
Initially, for x ¼ 0.2, the σdc is not much altered at the low temperature 3
region, but the slope changes noticeably for higher temperatures NðEF Þ ¼ (21)
4 π ðRMo =RV Þ3 Ehigh
(Fig. 10) owing to the interference introduced by Mo and V ions to
electronic motion. With further increasing V2O5 content, σdc increases The approximated and reasonable N(EF) values for localised states in
due to the MTIE, where the small polaron transport mechanism is the high-temperature region are noted in Table 3. With rising V2O5
concentration (x), the N(EF) of V ions advances more than for Mo ions,
which affirms the glass network modifier function of V2O5. The N(EF)
produces higher quantities of localised states inside the optical bandgap
or mobility gap, supporting the σdc hopping or diffusion process.
The frequency-dependent ac conductivity dispersion region is usu
ally examined using the Jonscher universal power-law equation as fol
lows [54]:
σðωÞ ¼ A ωs (22)
10
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
The modified NSPT model parameters, WH, τH, and T0, are obtained
Fig. 11. Variation of power-law exponent (s) as a function of temperature (T) by non-linear curve fit of the s–T plot with the modified NSPT model (Eq.
of all nanocomposites with different values of x. Solid lines are the non-linear (26)) at a constant f ¼ 1 kHz; the as-obtained values of WH, τH, and T0
best-fit data of (a) modified CBH model and (b) modified NSPT model. drop with an increase in V2O5 concentration (x) (Table 4). Therefore, the
increment in ac conductivity according to the NSPT model is the effect of
temperature (Fig. 11), which efficiently shows the most appropriate the reduction in values of WH and τH.
model for the ac conduction mechanism. Scaling property is used to identify the dependence of the conduc
The values of s decrease with rising temperature for x ¼ 0.2, 0.4, 0.6, tivity relaxation process on temperature and frequency as the scaled ac
and 1.0 (Fig. 11a), indicative that a correlated barrier-hopping (CBH) conductivity spectra merge into a solitary master curve [69,70]. The
mechanism is relevant for ac conductivity [66,67]. The s values increase frequency (X-axis) and ac conductivity (Y-axis) of each ac conductivity
with temperature increment (Fig. 11b), suggesting that a spectra are scaled or normalised by dividing by the values of ωH and σdc
non-overlapping small polaron tunnelling (NSPT) process is applicable individually [69,70]. The scaled spectra of a particular composition at
for ac conductivity for x ¼ 0.0 and 0.8 [63,68]. different temperatures are called temperature scaling spectra, as
In the CBH process, it is presumed that small polarons or charge depicted in Fig. 12a for x ¼ 0.8. For the present complicated
carriers do not interact with each other and are located in the crystallite
potential well, permitting nearest-neighbour jumps only [66,67]. The
small polaron jump or hopping probabilities with time can be presumed Table 4
as symmetric or asymmetric, i.e. back and forth hopping or jumping All values of fitting parameters of modified CBH model (for x ¼ 0.2, 0.4, 0.6, and
processes among the localised or defect sites or among localised states to 1.0) and the modified NSPT model (for x ¼ 0.0 and 0.8). Errors are indicated in
defect states. Due to the effect of thermal energy, the vibrational fre parentheses.
quency of electrons or small polarons at different well-defined sites Modified CBH model Modified NSPT model
enhances, leading to escalation of the SPH process. The localised or x (mol
WM (eV) τ0 (s) T0 WH (eV) τH (s) T0
%)
defect states also account for the ac conductivity process in response to (�0.004) (�0.05 � (K) (�0.014) (�0.05 � (K)
their closer proximity, in that way providing the distribution of different 10 –3) (�2) 10 –3) (�2)
heights of energy barrier separating them. It is predicted that separation 0.0 – – – 0.185 2.04 � 632
between two centres of a localised or defect state is higher than the 10 6
radius of the small polaron wave function of each site. According to the 0.2 0.089 1.38 � 604 – – –
10 3
CBH model, the relation of different parameters such as s, Debye
0.4 0.053 3.21 � 564 – – –
relaxation time (τ0), and temperature (T) is described by the following 10 3
expression [66,67]: 0.6 0.044 9.11 � 532 – – –
10 4
6 KB T 0.8 0.113 3.68 � 484
s¼ 1 (23) – – –
WM þ KB T ln ðωτ0 Þ 10 9
1.0 0.034 2.01 � 440 – – –
With a rise in temperature, the s value drops in Eq. (23), where WM is 10 5
11
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
Data availability
Fig. 12. Normalised plot of (a) temperature scaling spectra of ac conductivity Declaration of competing interest
for x ¼ 0.8 and (b) composition scaling spectra of ac conductivity at 693 K for
all nanocomposite samples with the variation in V2O5 concentration (x). The authors declare that they have no known competing financial
interests or personal relationships that could have appeared to influence
(quaternary) glassy systems, temperature scaling spectra comply with the work reported in this paper.
the time-temperature superposition (TTS) principle because all the
scaled spectra at different temperatures overlap almost perfectly [69, CRediT authorship contribution statement
71]. A wide distribution of electron trapping sites forms due to the ex
istence of structural disorder, with unlike energy barrier heights, Dipankar Biswas: Investigation, Methodology, Software, Valida
creating a wide-ranging distribution of back and forth hopping or tion. Anindya Sundar Das: Conceptualization, Visualization, Supervi
tunnelling rates that maintain the TTS principle [71]. The scaling pro sion, Writing - original draft, Writing - review & editing. Rittwick
cess for all the composites at a precise temperature is designated as the Mondal: Validation, Formal analysis. Anindita Banerjee: Software,
composition scaling spectra (Fig. 12b), which indicates that composition Data curation. Abhigyan Dutta: Resources, Supervision. Soumyajyoti
scaling spectra are not suitably merged into a single unique curve. Kabi: Resources, Funding acquisition. Debasish Roy: Resources. Loi
Conductivity scaling phenomena (i.e. temperature scaling and compo tongbam Surajkumar Singh: Project administration.
sition scaling) reveal that the conductivity relaxation process of mobile
charge carriers or small polarons does not change with variations of Acknowledgements
temperature but relies on the structural features of each composition, i.
e. on the molar ratio of V2O5 concentration (x). The authors wish to acknowledge Mr. Anish Karmahapatra and Ms.
Soma Roy of SINP, Kolkata, for their assistance in producing the XRD
4. Summary and conclusions patterns and UV–Vis absorption spectra, respectively. Dr. Abhigyan
Dutta thankfully acknowledges the financial support of SERB (Govt. of
Two ternary and four quaternary nanocomposite samples of India) through Grant no. EMR/2017/000325.
0.4ZnO–0.1P2O5–0.5[xV2O5–(1-x) MoO3] composition, where x varies
in steps of 0.2, were synthesised via the melt quenching process. An Appendix A. Supplementary data
inverse proportional relationship between ρ and VM was recognised. The
XRD analysis confirmed the development of ZnP4O11, ZnMoO4, MoP3O9, Supplementary data to this article can be found online at https://doi.
Zn2P2O7, Zn2.5MoVO8, MoV2O8, Zn3MoV2O11, V4[P2O7]3, VP3O9, org/10.1016/j.jpcs.2020.109505.
Zn2V2O7, VOPO4, and V2O5 nanocrystallites, superposing within the
amorphous glass networks. The dC and Eopt were in an inverse rela
tionship, whereas EU values were consistent with nanocrystallite size.
The formation of additional defect or localised states within the optical
12
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
13
D. Biswas et al. Journal of Physics and Chemistry of Solids 144 (2020) 109505
[68] S.R. Elliott, A theory of ac conduction in chalcogenide glasses, Philos. Mag. 36 [70] R. Mondal, D. Biswas, D. Mandal, A.S. Das, S. Kabi, Consequence of the heat-
(1978) 1291–1304. treatment on the ionic conductivity of silver-tellurite glass nanocomposites, J. Non-
[69] A. Ghosh, A. Pan, Scaling of the conductivity spectra in ionic glasses: dependence Cryst. Solids 506 (2019) 51–57.
on the structure, Phys. Rev. Lett. 84 (2000) 2188–2190. [71] B. Roling, A. Happe, K. Funke, M.D. Ingram, Carrier concentrations and relaxation
spectroscopy: new information from scaling properties of conductivity spectra in
ionically conducting glasses, Phys. Rev. Lett. 78 (1997) 2160–2163.
14