3. Studies of the Nonlinear Kinetic Equation
The phenomenological model of the ions diffusion transfer in the dielectric (protons in HBC) in the polarizing electric field is based on a system of the nonlinear equations of the Fokker–Planck and Poisson type (expressions (11), (14) of [
3], p. 13)
considering initial and boundary conditions ([
3], p. 13)
In (2)–(6), the following designations are adopted:
is an ion excess concentration over its steady-state concentration
;
,
,
,
,
,
, where
,
are diffusion and mobility factors (coefficients) for ions (protons in HBC) ([
1], p. 44, p. 75; [
4], p. 156);
—high-frequency dielectric constant of crystal; d—dielectric thickness ([
3], pp. 12–13).
We also introduce an additional dimensionless small parameter
< 1. It is important from the point of view of the further expansion of solutions to Equation (2) in series in the powers of the parameter
([
3], p. 16), which, in turn, determined the interaction parameter of relaxation modes, at the base frequency
[
3].
The authors of [
3] have built their solutions to the equation systems (2)–(6) by methods of the perturbation theory using power series ([
3], p. 12)
in the third approximation by parameter
. In [
3], based on the functions
,
,
by the method of mathematical induction the recurrence formulas have been devised for calculating relaxation modes
,
—volume density of charge
, calculated in the k-th approximation of the perturbation theory, at the frequencies
,
of the variable field, respectively, (expressions (18), (19) from [
3], p. 15). On this basis, the first two frequency harmonics
,
of the function
are calculated. They are expressions (20), (21) of ([
3], c.16) and the frequency harmonics of polarization ([
3], p. 16)
In (8),
is the dimensionless relaxation time for n-th relaxation mode, where, according to expressions
,
,
, the formula is built for the dimension relaxation time
determined with the help of
—diffusion relaxation time for n-th mode,
,
—Maxwell relaxation time ([
1], p. 44, [
3], p. 15, [
4], p. 60). In this case, the ratio of relaxation times allows us to construct the impotent equality of the semi-classical kinetic theory of the relaxation processes
, where the critical temperature separating temperature diffusion zones and Maxwell relaxation
may be calculated. Here we use the ratio for the diffusion and mobility factors
, with the expressions for the zeroth and linear approximations of the perturbation theory in terms of the small dimensionless parameter
for the rate of the probability of the transition of the relaxer through the potential barrier. Therefore, according to (1)
we can write the transcendental equations for critical temperature
Here, is determined from equality , another critical temperature separating the regions (areas) of tunneling (quantum), when , , and thermally activated (classical), when , , transitions of relaxers through a potential barrier.
An attempt to apply a similar method to the calculation of the function
in the approximation of k at the frequency
gives such bulky expressions
,
, etc., that the conclusion of the recurring formula
requires a fundamentally new, more general approach. The first and simplest member of this series is the expression
formulated ([
3], p. 17). Let us move on to the implementation of the general analytical method. Then, substituting the series (7) into the system (2)–(6) gives
Here
,
([
3], p. 13). In all subsequent approximations
We expand
into the Fourier series according to the orthogonal functions
, on the segment
:
The function
has the meaning of the relaxation mode of the n-th order k-th approximation of the perturbation theory,
is the complex amplitude of the relaxation mode
. According to
staking into account (13), rewrite (10) in the form of an operator equation
where
. Integrating (15) with
, we get
Using additional decomposition
and denoting
, we convert (16) to the recurring formula
Since the generation of relaxation modes with complex amplitudes
begins with the first order of perturbation theory [
1]
from the first term of (17) it is not difficult to see
where we get the recurring formula
To develop a recurring expression for the complex amplitudes
of higher orders (multiples of frequency
), we rewrite the expression (17) as
In Equation (18) uses integral operators
As the generation of the relaxation modes with the amplitudes
begins from the second order of the perturbation theory, according to (18)
In (18a) uses integral operators
According to (17a)
we get the recurrent relation from (18a)
where
.
As the generation of the relaxation modes with the amplitudes
begins from the third order of the perturbation theory, according to (18)
In (18c), integral operators
have been introduced.
We get the contribution to amplitude
from the consequent recursion (18c), taking into account (18b)
where
Substituting of (18d) from (18c) into addend gives
where the notation
is introduced. It is used in numerical calculations as
From (18c) it directly follows that
only for the odd modes n. We consider the contribution of the amplitude
from the third component of recursion (18c). For this purpose, from (18b) we calculate
The contribution from product of will be equal to zero, since for any p values.
Non-zero contribution gives only the product
. After that, from (18c)
The total contribution to the amplitude
from the terms
and
gives
Expression (18g), by virtue of (18c), is nonzero only for odd modes, therefore, it is further necessary to take into account the contribution to the amplitude
from the first term of the recurring expression (18c). Each expression (18g) of the g-th approximation, substituted in subsequent approximations in the first term (18c), provides the following contributions to the amplitude
:
In (18h), the sum of the series
, by virtue of (d1), reduces to the form
In (18h), the sum of the series is also introduced.
The total contribution of expression (18h) to the complex amplitude in the approximation is determined by summing the elements by g from g = 3 to g = k − 1: .
The complete expression for the complex amplitude
becomes
Further, representing (7) in the form of
,
, we will write the expansion of the frequency harmonic number
r for the function
in power series by the powers of the parameter
The decomposition components (19) are represented by the k-th order relaxation modes of the perturbation theory generated at
frequency, respectively
Expression (20), by virtue of (14), has the meaning of superposition of the relaxation modes of the n-th order
with the amplitudes calculated according to the recurring formula (18), taking into account (17a), (18b), (18j). Then, from (19), (20) we obtain
When
r = 1,
r = 2, from (21) we obtain
and
. Taking (17a), (18b) into account, we confirm the expressions (20), (21) from [
1,
7].
When
r = 3, according to
and taking account of (18j), we obtain
From where, after calculating the sums of the series, we get
Based on the expression [
1]
taking account of (23)
Expression (25), in addition to (8), confirms that odd r relaxation modes , as well as , give nonzero contribution to polarization. Manifestation of this pattern at higher frequencies ,,…, is obvious.
4. The Complex Dielectric Permittivity
The dielectric polarization, in the general case, calculated as
In (26)
is a complex function defined on set
of parameters of relaxers and temperature
, which has meaning of
order component from decomposition of the polarization
into a series according to odd degrees of intensity of the polarizing field
. Here
. The frequency harmonic of the polarization, of order
, according to (7),
, taking into account (20), becomes
Comparing the formula
([
4], p. 151) with the expression (26), let us write the complex dielectric susceptibility (CDS) in the form of
In (28) .
Taking the relaxation complex polarization as
([
4], p. 151), we find the complex dielectric permittivity (CDP)
Substitution (28) into (29) gives
In (30)
; limiting the approximation of
,
we find the CDP with accuracy to the quadratic term along the field
Substituting (8) and (25) in expressions of the form
we obtain
where relaxation parameters are introduced that characterize the first frequency harmonic of the polarization (a multiple of the first power of the alternating field (8))
The expression
formally coincides with the formulae (2) of ([
1], p. 44), which indicates the validity of the “zero” field of approximation and consistent with the expressions obtained in [
3] for CDP ([
3], p. 17).
Also, based on (25), we obtain a multiple of the second power of the alternating field component of the complex permittivity
In (34) the following notations are adopted:
Based on (32), taking (33) and (34) into account, we obtain
The abridged notations are adopted in (35):
Separating the real and imaginary parts in (35), we get
It is easy to see that the first terms in the right parts of (36), (37) also formally coincide with the CDP components calculated in the “zero” field approximation ([
1], p. 44; [
3], p. 17).
The second terms in the right parts of (36), (37), respectively, reflect the influence of the polarizing field on the CDP spectra, in the regions of both strong fields (10–1000 MV/m) and ultra-high temperatures (550–1500 K)—high-temperature volumetric-charge polarization; and weak fields (100 kV/m–1000 kV/m) and ultra-low temperatures (1–10 K)—low-temperature quantum tunnel polarization.
The functions
,
reflect the non-linear effects of the third order on the field frequency associated with the interaction of the relaxation modes generated at frequencies
,
,
with the corresponding relaxation times
. An additional relaxation mode generated at the frequency of the
field with relaxation time
appears in functions
,
. At the same time, the parameters
,
reflect the non-linear effects of the first order in frequency, with relaxation times
. Parameters
,
of (18h), (23) are a measure of the effect of the non-linear effects of the third order (interaction between the relaxation modes at frequencies
,
,
and corresponding relaxation times
) on the polarization. So, provided that the parameters
,
,
are small, the expression (23) is determined by the relaxation mode
, corresponding to the third order of the perturbation theory (k = 3). This is the simplest (by the structure of complex amplitudes
) of the set of modes
generated by the frequency
, which is consistent with ([
3], p. 17).
Numerical calculation according to the formulas (36), (37) will allow analysis in more detail of the influence of the non-linear effects of the third order on the field frequency, on polarization, depending on the values of the parameters
,
,
,
([
4], pp. 150–152) and the thickness of the dielectric.
The mathematical model developed in this paper is suitable for theoretical studies, previously experimentally installed in crystals of zirconia alumina ceramics (ZAC) (ZrO
2-Y
2O
3)-n (Al
2O
3), at the temperature of T = 1250 K and frequency of 1 kHz, abnormally high dielectric permittivity ε = 5
10
6 [
5]. The nonlinearities of the first and third order can be considered as the reason for the semi-ferroelectric effect [
6] associated with the rearrangement of the oxygen subarray in the ZAC, near the phase transition point.
The areas closest to electrical and radio engineering in the practical application of this model are laser and fiber optic technology: femtosecond lasers, optical light manipulators, incident radiation direction detectors, differential phase and frequency analyzers, etc. [
8,
9,
10,
11].
5. Effect of Temperature on Theoretical Spectra of Complex Dielectric Permittivity
Further study of the analytical dependencies (36), (37) is a separate theoretical task, the solution of which, together with expressions (33a), (34a)–(34h), will allow us to establish in more detail the influence of the nonlinear effects on dielectric loss tangent spectra and thermos-stimulated polarization and depolarization currents. It depends both on the patterns as macroscopic relaxation processes occurring throughout the dielectric and related to polarizing parameters (external) electric field (EMF amplitude and frequency) and temperature and microscopic processes associated with transitions of main charge carriers (ions of both charge signs in the general case, and protons in HBC) through a potential barrier. The parameters of the potential barrier (reflecting the molecular mechanism of interaction of the conduction ions with the crystal lattice, and protons with an anionic sublattice in HBC) significantly determine the efficiency of both classical (thermally activated) and quantum tunnel transitions of microparticles (relaxers) through the potential barrier. Within the framework of the quasi-classical model, they determine the influence of classical effects on the nonlinear volume–charge polarization in the field of sufficiently high temperatures (250–550 K) and quantum effects on tunneling polarization manifested in HBC due to the quantum relaxation movement of protons in the region of low (50–100 K) and ultra-low (4–25 K) temperatures.
The analytical expressions for the relaxation coefficients (33a) calculated at the base frequency of the alternating electric field in the areas of diffusion and Maxwell relaxation, away from the critical temperature , should differ significantly in the type of functional dependence of these parameters on temperature. In this case, we will analyze the dispersion expressions (36), (37) taking into account the temperature dependencies for the complete relaxation time ,.
We convert the “zero” component of the complex dielectric permittivity (CDP) by the polarizing field from (32), taking into account (33), into
. From where, separating the real and imaginary components, we obtain the dispersion ratios for the CDP components calculated in the nonlinear approximation of the perturbation theory (up to infinite approximation by the small dimensionless parameter
) at the main frequency
of the external polarizing alternating electric field, taking into account the effects manifested, within the phenomenological model of the relaxation polarization in the interaction of the first-order relaxation modes by the field frequency
In (38), the relaxation coefficients (33a) are a converging infinite functional series, obtained as a result of solutions of the Fokker–Planck Equation (2) together with the Poisson Equation (3) in an infinite approximation of the perturbation theory by a small dimensionless parameter
according to Equation (7), at the fundamental frequency of the alternating electric field and thereby are the consequences of sufficiently stringent solutions to the generalized nonlinear kinetic Equation (2). They allow us to calculate the spatially inhomogeneous bulk density of the excess electric charge in arbitrary approximation (multiplicity
r) on frequency harmonics (21). In particular, cases that are used to calculate the values of polarization of the dielectric at the first two frequencies (r= = 1, r = 3) from expressions (8) and (25) are important for comparison with the experiment. In this case, temperature dependencies of the relaxation parameters (33a) corresponding to the first frequency harmonic of polarization (8) are similar to temperature dependencies for the simplest forms of the relaxation parameters (34h), corresponding to the third frequency harmonic of polarization (25). This allows us, with a certain degree of accuracy, to judge the significant influence of the temperature range on the mechanism of the nonlinear relaxation processes when generating polarization of the dielectric in the wide range of field parameters. It is obvious that in the region of weak fields (0.1–1 MV/m) the non-linearities will appear due to low-temperature (10–50 K) polarization effects (quantum diffusion polarization), and in the region of strong fields (10–1000 MV/m) the non-linearities will appear due to high-temperature (550–1500 K) polarization effects (nonlinear bulk-charge polarization of the mixed type) [
7,
12,
13,
14,
15,
16].
In the field of Maxwell relaxation
, when
, taking
, due to the identity
, from (33a)
convert (38) into
Here , .
Tangent of dielectric loss angle
[
1], according to (39), takes asymptotic form
Substituting
[
1,
2] in (40) and comparing the obtained result with the known expression for conductivity loss
[
2], we find the coefficient of
ion (in HBC-proton) relaxation electrical conductivity
In the field of the diffusion relaxation
, the parameters (33a) converted into
at the area
, are reduced to equalities
In the region of ultra-low frequencies, when
, based on (42), taking into account the identity
, passing to the approximate equations of the form
,
and using formal expressions of the form (38)
we obtain
Further, in the diffusion relaxation
, when the values of the Debye shielding radius are much larger than the crystal thickness
[
1,
2], taking
, from (44) we obtain an expression for the real component of the CDP corresponding to the classical law for dipole polarization
Another criterion of validity of formulae (8) and (25) of the nonlinear approximation of the perturbation theory investigated above is their ultimate transition to results of the linear kinetic theory of the relaxation polarization [
7,
12,
13,
14,
15,
16]. To this end, we simplify the expression for the polarization (8).
Taking the condition
in (8), by virtue of
we obtain the expression
corresponding at the fundamental frequency of the field
, the first approximation by parameter
[
17], where
. Considering
[
1], calculating
;
, we find
In the case of
,
, by substituting
,
in (46), we obtain the expressions (39) approximated in the first approximation by the parameter
From this, by virtue of
,
, entering the designation
, we have
which is consistent with the results of the linear theory of dielectric losses [
1,
2] in the temperature region
.
In the case of
,
, by rewriting (46) while taking into account (42), we have
From here, by virtue of
,
,
, entering the designation
, we have
which is also consistent [
2] in the temperature region
.
In the theoretical studies of the solution schemes of the nonlinear kinetic Equation (2) and its annexes to the issues of studying the mechanisms of the relaxation polarization described by formulas (8), (25), (21), it is impossible not to consider the question of the generalized calculation of the relaxation parameters (33a) for the wide range of field and temperature parameters.
Introducing dimensionless variables
,
, we present the generalized expressions (33a) in the form of
It is convenient for further analysis of the effects of various kinds of the polarization nonlinearities (interactions of the multiple relaxation modes of both the same frequency and different field frequencies; quantum microscopic effects associated with proton tunneling in HBC; the phenomena of nonlinear volume–charge polarization, etc.) and temperature per dispersion expressions (38) with numerical calculation of frequency–temperature spectra of the complex dielectric permittivity (CDP).
When calculating by formulae
it becomes possible, regardless of the properties of the parameters
, , in a certain temperature range (see particular cases
,
, from (42) and
,
), to convert the sum of infinite series to an explicit analytical form. So, calculated in the functions of temperature and frequency of an alternating polarizing field, infinite sums of converging series
taking into account identities
according to
take the form
Special cases have application significance. In the field of the
diffusion relaxation, at temperatures far from critical
, in a wide range of the variable field frequencies
, when the conditions
,
,
are met, taking the variable
and moving to (49) and, respectively, in (50), (51), in the area of ultra-low temperatures
, to the limits of
,
, we write approximate expressions
In the ultra-low frequency region, according to
, taking in expressions (52), (53)
and, moving to the limits at
, we obtain the expressions
,
. They can act as a criterion for the validity of the results of previously performed calculations when moving from equalities (49) to (50), (51) and from (50), (51) to (52), (53), since the expressions for the
,
limits completely coincide with the results of the direct calculations of the sums of infinite series resulting from their (49) while meeting the formal conditions
,
when
,
. Moreover, it is not difficult to see that built formally, on the basis of (50), (51), for reasons of the smallness of a parameter
expressions (52), (53), in physical relation, fully satisfying the conditions of the diffusion relaxation
, coincide with the parameters
,
from the equalities (42). Therefore, after substitution (52), (53) in expressions (43), the generalized dispersion expressions can be obtained that are suitable for calculating the components of the complex dielectric permittivity (CDP) in the field of diffusion relaxation at arbitrary frequencies of a variable field according to generalized formulae (43). Then, based on the partial spectral expressions (46), approximate formulas of the linear approximation of perturbation theory (46.4) for the CDP components in the diffusion relaxation region
, far from the critical temperature
but not at absolute zero, according to (52), (53), are reduced to
a transcendent form
It is not difficult to see that the expressions (54), (55), in comparison with (46d), (46e), are more obvious functional dependencies of the CDP components on the frequency of the alternating field (against the background of a weak influence of temperature). They are also more convenient for numerical calculations of the theoretical frequency spectra of the CDP and analysis of the effects of field and temperature parameters on the tangent spectra of the dielectric loss angle in the diffusion relaxation region, far from the critical temperature , at the main field frequency, in the linear approximation of perturbation theory.
The study of the laws of dielectric relaxation at arbitrary temperatures (including close to critical
) when
, in the low frequency range of the variable field
, when in expressions (49), and accordingly in (50), (51), at the frequency of the field tending zero, the conditions
,
,
,
can be assumed, is an important question for theory. From there, entering the variable
and, going to (49) and (50), (51), to limits
,
we will write approximate expressions
In this case, the analytical expressions (56), (57) are functions only of temperature.
Also note that for the expressions (56), (57), as for the expressions (52), (53), the criteria of reliability of results obtained during the formula transformations (49)–(50), (51) and (50), (51) to the form (56), (57) are met. This is confirmed by the transitions in formulae (56), (57) to the limits , , which completely coincide with the results of direct calculations of the sums of infinite series resulting from (49) when the formal conditions , are met simultaneously, when , .
Note that, although the expressions (56), (57), when performing inequality
, in a certain relation satisfy the conditions of the diffusion relaxation
, in the general case (by virtue of
), do not coincide with the parameters
,
of equalities (42), because formulae (56), (57) are constructed primarily for the low frequency range of the variable field (
). This means that, for a mixed type of relaxation (near the critical temperature), when calculating the components of the complex dielectric permittivity (CDP), it is advisable to substitute equations (56), (57) in expressions (38). Then, based on the partial spectral expressions (46), the approximate formulas of the linear approximation of the perturbation theory for the CDP components, at ultra-low frequencies (
), according to (56), (57) and (48), are reduced to the form
The conversion of formulae (48) for the case
makes them more convenient for comparison with the experiment. To do this, rewrite (48) using (56), (57), entering the variable
. Then
On the other hand, from (47), with
, we get the same result as in (60a)
Here .
It is not difficult to see that the identity (60.1) makes it possible to significantly simplify temperature dependencies in generalized expressions (38)
Expression (61) is suitable for investigating the temperature dependencies of the real component of the complex dielectric permittivity
over the wide temperature range at the variable field frequency of zero (stationary relaxation polarization). It results from a generalized variance expression (38) constructed at the base frequency of the variable field in an infinite approximation of the theory of perturbations over a small parameter. It contrasts the expression (57), constructed on the basis of linear expressions (46) in the first approximation of the theory of perturbations by the parameter. This expression is more limited in terms of physical completeness of coverage of the nonlinear polarization effects associated with the effects of the interactions of multiple relaxation modes of the basic field frequency on the dielectric relaxation mechanism in the low and ultra-low temperature region. Thus, the expression (61), which is generalized to the wider range of the physical nonlinearities, is more effective in detecting the effects of quantum phenomena (due to tunnel transitions of relaxers (protons) between anionic sublattice ions) on the mechanism of low-temperature migration polarization (50–100 K) in proton semiconductors and dielectrics and their nanoscale layers (1–10 nm) at ultra-low temperatures (1–10 K). Within the framework of this scientific paper, we do not investigate the temperature relationships of the statistically averaged quantum transparency coefficient of a potential barrier for various configurations of the energy spectrum (discrete degenerate or non-degenerate; quasi-continuous) tunneling protons depending on the parameters of the undisturbed potential crystal lattice field (hydrogen bond field). Detailed studies of this issue carried out in the works [
1,
2,
4,
14,
15] showed that the nature and properties of the quantum transparency coefficient for protons in the hydrogen-bonded crystals (HBC) are most qualitatively disclosed for the discrete degenerate model (split into sublevels within individual energy zones, as a result of the effect of neighboring potential pits on a proton localized in the region of a given potential pit). The numerical values of the statistically averaged quantum transparency coefficient are significantly dependent on the parameters of the potential barrier (especially from its height, which makes sense of the activation energy), crystal thickness and temperature range. Thus, in [
16], it was established that accounting for splits of proton energy levels in the HBC on sublevels (in an unperturbed potential field) with the proton activation energy of 0.01 eV leads to a shift in the temperature maximum of the statistically averaged transparency coefficient from 55 K with an amplitude of 0.12 (in the absence of energy spectrum degeneration) to 150 K with an amplitude of 0.14 (in the presence of energy spectrum degeneration). Ref. [
14] shows the effects of the thickness of the crystal layer on quantum transparency and dielectric constant of the crystal of the HBC type at nanoscale thicknesses of the sample, when, in layers with a thickness of 10 nm, in the case of the zone structure of the energy spectrum of protons with an activation energy of the proton of 0.01 eV, there is an increase in the amplitude of the quantum transparency coefficient to 0.98–0.99 (proton superconductivity) at a maximum temperature of 10–15 K, and the amplitude of the stationary dielectric permittivity
calculated according to formula (61) reaches abnormally high values (2.5–3.0 million) [
16]. This indicates the formation of the ferroelectric state in HBC due to quantum tunnel movement of protons in hydrogen sublattice of HBC nanofilms in the region of liquid helium temperatures.
The effects of the quantum effects on the frequency spectra of the complex dielectric permittivity should be carried out according to the generalized dispersion formulas (38), taking into account the most common expressions for relaxation coefficients in the form (50), (51). This will allow taking into account the influence of all the above-mentioned parameters of the crystal lattice structure and the parameters of the energy spectrum of protons, crystal thickness and temperature on the parameters of maxima of the CDP by frequency and amplitude. The advantages of expressions (38), compared to the particular linear equations (46) are that formulae (39) are constructed, as noted above, at the base frequency of the alternating polarizing field, everything up to infinite is taken into account, approximation, by a small dimensionless parameter of the perturbation theory when solving the Fokker–Planck Equation (2) and, as a consequence, the results of the effects of the spatially inhomogeneous polarization-induced electric field in the dielectric. This is very important in studies of the effects of the volume–charge distribution of relaxers in the mixed-type relaxation region, manifested in the region of sufficiently high temperatures (250–450 K), when, against the background of the dominant mechanism of thermally activated proton transitions [
1,
2,
15,
16], proton tunneling continues to make a significant contribution to the kinetics of nonlinear thermally activated depolarization [
1,
2]. As noted in [
3,
4], the polarization nonlinearities associated with the formation of charge–volume polarization in dielectrics with the ion-molecular type of chemical bond and manifesting in the region of ultra-high temperatures (550–1500 K) and strong electric fields (10–100 MV/m) [
7,
14], result from a mixed type of relaxer movements (protons) activated both by classical (due to thermal motion and the interaction of protons with the crystal lattice) transitions and quantum tunnel transitions.
Numerical studies and analysis of the properties of the frequency spectra of the complex dielectric constant, performed according to the generalized dispersion expressions (39) and (50), (51), are a separate task, the results of which will be presented in a future work.
6. Properties of Theoretical Spectra of Complex Dielectric Permittivity in the Region of Maxwellian Relaxation
From the point of view of organizing and planning experiments to study the electrical properties of circuit elements of control and measuring and radio-electronic equipment used in the electric power industry, insulating and cable technology and industrial electronics, it is relevant to study the frequency–temperature spectra of the dielectric loss tangent in the temperature range T = 150–550 K [
2], when the main contribution to the ion-relaxation polarization in dielectrics is made by the Maxwellian relaxation of the main charge carriers (in ionic dielectrics-interstitial ions and in HBC-protons localized on hydrogen bonds).
In the area of Maxwellian relaxation at temperatures higher than critical
, when
[
1], by approximating parameters (33a)
with respect to the dimensionless parameter
, we bring the expressions (38)
to the form
Here .
The tangent of the dielectric loss angle at the fundamental frequency of the field
is
Figure 1 and
Figure 2 show the dependencies graphs
,
for various values of the parameter
, calculated using a computer program, according to formulas (64a), (64b), (65). With this representation, relaxation processes that have the same values
, are depicted on the graph
as one curve. In accordance with
Figure 1, an increase in the parameter
leads to an increase in the maximum and to a shift in its position to the area of small values, while at
the maximum position changes slightly with increasing
and remains approximately identically and equal to
. With small values of the parameter
the maximum position
determined by the criterion
or
. The curves
in the area of large values
are nearly congruent.
In conclusion, it should be pointed out that the anomalously high dielectric permittivities ε = 5 million found experimentally in [
5], in samples of corundum–zirconium ceramics (CZC), with an alternating field frequency of 1 kHz, at the point T = 1250 K, can be theoretically explained and further investigated at a higher analytical level using expressions (64a), (64b), (65) that indicate the correspondence of the formulated theoretical methodology to the experimental regularities that manifest themselves in high-temperature ionic superconductors near the second-order phase transition temperature (quasi-ferroelectric effect), which is topically for the design of the theoretical methods for forecasting the nonlinear electrophysical properties of hydrogen-bond ferroelectrics (KDP, DKDP) used in laser technology as regulators of electromagnetic radiation parameters and electric gates [
1,
2,
3,
4].
In this regard, promising are the designs of high-speed nonvolatile memory devices based on thin films of ferroelectrics with a rectangular hysteresis loop (RHL), characterized by anomalously long times remanent polarization relaxation, increased mechanical performance and thermal stability.
For the practical application of the calculation formulas (64a), (64b), (65), it is convenient to express them in terms of the static permittivity
and permeability at high frequencies
. In order to obtain such phenomenological relations, we will use the expression obtained from (61) in the Maxwellian relaxation region, when
and, at
, by virtue of
, we have
, where
. Then, based on (64a), (64b), (65) we have
In formal form, the dispersion relations (66) resemble the classical expressions for the Debye dispersion [
3], but differ from them in coefficients, which leads to other expressions for determining the maxima of the functions
и
, as well as their values.
The Debye expressions for the complex permittivity give the maximum value of the function
, reduced to
, compared with (66), but under the condition that is typical for dielectrics with a large depth of dispersion, when the magnitudes of the maxima will be approximately the same. The Debye relations [
3] for determining the position of the maxima
and
are several
times higher than the corresponding values calculated by formulas (66), (67). The Debye relations [
3] for determining the position of the maxima of functions
,
are
several times higher than, the corresponding values calculated by the formulas (66), (67). Thus, Debye expressions, which are used quite often to determine the parameters of defects during space charge relaxation, lead to a significant error in determining the relaxation time, increasing it by a factor of
.
8. Conclusions (with Elements of Results Analysis)
1. By the methods of the quasi-classical kinetic theory of the relaxation polarization and conductivity in crystals with ion-molecular type of chemical bonds (proton–relaxation polarization and proton conduction in HBC), taking into account the quantum effects caused by tunnel transitions of the main charge carriers or relaxers (first of all, protons, in HBC) on chemical bonds, for the model of a one-dimensional potential crystalline field, the transcendental Equation (8c) was built. (8d) allows the calculation of the critical temperature separating the temperature regions (zones) of diffusion and Maxwell relaxation and establishes the effects of molecular parameters of relaxers (lattice constant a; potential barrier width ; activation energy (potential barrier height) ; equilibrium concentration of relaxers (in the state of thermodynamic equilibrium in the absence of an external disturbing field); natural frequency of relaxer oscillations in an unperturbed potential well ). It has been found that the critical temperature for dielectric relaxation (dependent, in addition to molecular parameters, also on the thickness of the crystal) is associated with another previously calculated critical temperature for microscopic transitions of the relaxers over the potential barrier , which separates the areas of tunnel and thermally activated relaxation by ions (protons in HBC).
2. The analysis of basic positions of the physical and mathematical model of the nonlinear relaxation polarization in the hydrogen-bonded crystal (HBC) was performed [
1,
2,
3,
4]. It was established that the existing methods [
1,
3] are incorrect when calculating the polarization of the dielectric at the frequencies of the field, which are multiplied by fundamental frequency
, starting with the second frequency harmonic
. An improvement in the methodology for solving the system of equations of the basic phenomenological model (expressions (2)–(7)), in relation to the regions of abnormally high non-linearities (1–10 K, 100–1000 kV/m; 550–1500 K, 10–1000 MV/m) is required ([
1], p. 43; [
2], p. 74).
3. Generalized (in the k-th approximation of the perturbation theory, by parameter ) analytical solutions of the nonlinear system of the Fokker–Planck and Poisson equations with blocking electrodes at the boundaries of the crystal were obtained. Recurrence expression (18) was constructed to calculate the complex amplitudes of the relaxation modes generated in the stationary polarization mode at the frequency (20). Volumetric charge density was calculated with accuracy to the third frequency harmonic of the variable field (22), (23). In this case, is the first, is the second and is the third field approximation.
4. The generalized (in case of arbitrary approximation by the frequency of variable field
) non-linear by the field
expressions for complex dielectric permittivity (CDP)
calculated at multiple frequencies
were constructed in the form of decomposition into a series by even frequency harmonics of variable field
. Here
. A particular case at
was examined in [
1,
3]. The analytical expressions for the complex dielectric permittivity (CDP) were obtained with accuracy
on the set
in a quadratic field approximation. They reflect the non-linear polarization effects of the third order in field frequency, which can serve as a theoretical basis for creating microwave generators by tripling the frequencies of the radio range.
5. Based on the expression for dielectric polarization at the fundamental frequency of the alternating polarizing (external) field (see formula (8)) the generalized nonlinear analytic expressions are constructed (38). They allow us to examine and numerically calculate theoretical frequency and temperature spectra of the complex dielectric constant of the dielectric at the main frequency of the polarizing field, taking into account the nonlinear polarization effects caused by the effects on polarization of the spatially heterogeneous electric field induced in the crystal during volume–charge polarization (in the region of high temperatures) or due to the tunnel movement of ions (primarily protons), in low-temperature quantum polarization (in the low-temperature region) in the hydrogen-bonded crystals (HBC). Formally, according to the spectral expressions (38), the nonlinear polarization effects in the dielectric appear in the form of additional terms reflecting the interactions the of relaxation modes with different (by the Fourier order number of the harmonic of the bulk charge density) relaxation times at the fundamental frequency of the external field . In this paper, the relaxation , parameters are introduced for the first time. These parameters are dimensionless values written in the form of converging number series in the functions of the frequency of the alternating polarizing field and temperature (33a) and, further, converted to analytical expressions (48) calculated together with (50), (51) in the functions of dimensionless parameters , . Further studies of the expressions (50), (51), for particular cases of diffusion relaxation over the wide range of field frequencies (, ) in the form of (52), (53), and for the case of mixed type relaxation in the low frequency region of the field () in the form of (56), (57) confirms both a high degree of mathematical accuracy in justifying the expressions (51), (52), and the physical meaning laid down in the preceding equalities (33a).
6. In parallel with the results indicated in paragraph 5 of the conclusions, based on the expression
, we introduce the Debye screening length
. This allows dividing the temperature regions (zones) of relaxation of the diffusion type
, when the relaxation parameters (33a) are written as
and the Maxwell type
, when the relaxation parameters (33a) are recorded differently
Additional criteria of the validity of the dispersion expressions (38) and related relaxation parameters (33a) are established. So, for the case of diffusion relaxation at temperatures much lower than the critical
, according to the limit expressions (42), in the low frequency region of the field
, when equals (42) are approximated to the form
,
, polarization goes into a quasi-stationary mode. According to the dispersion expressions (43), which are the same dispersion expressions in (38) written formally, for the diffusion type case
relaxation, the real and imaginary components of the dielectric constant are converted to a form corresponding to classical dipole polarization
The only clarification is that here the relaxation time ratio
is a more complex function of temperature T than in the classical model. It takes into account not only thermally activated (classical) transitions of particles (relaxers)
, but also tunneling of particles (relaxers) through a potential barrier (see item 1 of the conclusions)
which is reflected in the expressions (8a), (8b) for the transition probability rates of the relaxers
,
. The values included in the right part of the formula for the parameter
are disclosed in the item 1 of the conclusions and in
Section 1 and
2 of this paper. Here, a constant
with a temperature dimension is adopted. In particular, in the case of the dominant contribution of thermally transitions of relaxers in the process of polarization formation
, it is mentioned in item 1 of the conclusions by the critical temperature
separating the areas of diffusion
and Maxwell relaxation
.
7. It was found that for the case of Maxwell relaxation at temperatures much higher than critical
, the above partial
formulas for relaxation coefficients
,
take the limit form of
In this case, the dispersion expressions (38) are converted to the analytical form (39) corresponding to the generalized spectra of the complex dielectric constant of the crystal, calculated by the methods of the quasi-classical kinetic theory, on the fundamental frequency harmonic polarization (8) formed at temperatures much higher than the critical relaxation temperature
The expression constructed on the fundamental frequency harmonic of polarization (8) for the tangent of the dielectric loss angle in the temperature region
is reduced to the form
corresponding to the conduction losses caused by macroscopic effects associated with the through movement of ion relaxers (ion conductivity currents) in the dielectric in the region of ultra-high temperatures (550–1550 K). This is consistent with experimental data on measuring the temperature spectra of the current density of the thermally activated depolarization (TSCD) in the hydrogen-bonded crystals in the vicinity of an ultra-high temperature maximum (450–550 K). Note that the application of the quasi-classical (corresponding to the system of nonlinear equations (2), (3)) formulas for Maxwell’s relaxation time
, where
calculated in linear approximation by the parameter
, the ion mobility coefficient [
61,
62] gives the expression for the tangent of the dielectric loss angle in the temperature region much higher than the critical relaxation temperature
to the form
(see formula (41)). Here, the value
makes sense of the conductivity coefficient for ions calculated as a function of the temperature and crystal structure parameters (see listing of these parameters in the item 6 of the conclusions)
8. From the results noted in the conclusions in item 6, it was found that in the low-frequency range
, during diffusion relaxation at temperatures near the critical
, while maintaining the above approximation
, when
, the real component of the dielectric constant reaches abnormally high values of
and in dielectrics with the complex crystal lattice structure (mica, talc, vermiculites, ceramics, etc.) there is a quasi-ferroelectric effect established earlier experimentally for samples of corundo–zirconium ceramics [
13].
In the case of the dominant contribution of thermally activated relaxer transitions , when , by virtue of , the point of phase transition of the crystal to the quasi-magnetic electric state is determined from the expression .
In general, when , by virtue of , by analogy with (8c), (8d), a transcendent equation generalized to the wide temperature range is constructed to determine the point of the phase transition of the crystal to the quasi-hydroelectric state . After transformations, this expression takes the form . Comparison of numerical values of temperatures and is established from the solution of Equation (8c) or (8d).
9. The theoretical scheme has been developed for a more stringent, in comparison with item 7 of the conclusions, analysis of the behavior of the real component of the dielectric constant in the low-frequency range , with a mixed dielectric relaxation mechanism () when the relaxation parameters (33a), (47), (48) are converted to the form (60a), (60b). This, as in the case of (42), also allows us to proceed to the study of the quasi-stationary polarization mode, according to the condition , when , . Variable has been entered. The first of the dispersion expressions (38) is brought to a stationary form, which made it possible to present the generalized analytical expression for static dielectric constant in the form convenient for comparing the results of the theory and experiment . This is a relevant result for this field of research, allowing us to identify, at the theoretical level, in the wide temperature range (0–1500 K), the influence of the microscopic and macroscopic mechanisms of relaxation processes in dielectric (diffusion and Maxwell relaxation; tunnel and thermally activated relaxation; electret effect; quasi-ferroelectric effect, etc.) on the nonlinear polarization effects associated primarily with the interactions of relaxation modes of the volumetric charge at the fixed frequency of an alternating field. In this matter, we limited ourselves to the fundamental frequency harmonic of polarization in the form of (8) in the framework of this paper. The study of the non-linear relaxation processes on the frequency harmonics of the polarization of the higher orders of magnitude, at least from the third (see formula (25)), etc., is the subject of a separate research and will be carried out in subsequent works.
10. The results obtained in this paper relate primarily to the category of applied theoretical studies in the field of dielectric physics and condensed state physics. They are aimed at improving existing and developing additional schemes for the theoretical analysis of nonlinear kinetic phenomena associated with the relaxation transfer of the most mobile charge carriers (relaxers) in the crystal structure of dielectrics with an ion-molecular type of chemical bonds (ceramics, layered silicates, crystalline hydrates, perovskites, vermiculites, clay minerals, halloysite, etc.) when polarizing dielectrics in variable external fields. The results obtained in this paper described in claims one to eight of the conclusions are relevant and scientifically significant for a number of fields of modern science. In particular, we could point out their potential application in the development of theoretical ideas about quantum mechanisms for the formation of a ferroelectric state due to the tunneling of protons in the hydrogen sublattice of KDP, DKDP crystals, which are of significant practical importance for laser technology (regulators of radiation parameters), nonlinear optics and optoelectronics (electro-optical sensors, the sensors of the deformations of hard rocks and building materials, the parts of machines and mechanisms, etc.). They are also significant for information technologies (thin-film ferroelectric elements of microcircuits for fast-acting non-volatile memory devices with abnormally long residual polarization retention time (up to 10 years)). A description of the prospects for the practical application of the results obtained for various branches of modern technology and technology is made in
Section 7 of this paper.