1. Introduction
In recent years, the intensive development of nanotechnologies has provided a significant stimulus to clarify patterns of the size effect in systems with ultradisperse particles of different topologies and dimensions. As an example of such systems, consider composite materials that are polar dielectrics embedded in porous structures having a branched network of through pores. Here, the embedded substances form either a system of isolated particles or a complex dendritic structure determined by the size and topology of the pores in original matrix, as well as by the surface tension, wettability, etc. [
1].
The effect of “restricted geometry” manifests differently in the case of various ferroelectrics. It is assumed [
2] that a decrease in the sizes of the ferroelectric crystal leads to an increase of the surface energy part in the total energy of the particle, an increase in the energy of the depolarizing field, and the energy associated with mechanical and chemical interactions between the ferroelectric particles and the substrate or matrix.
A combination of these mechanisms can result in both a decrease and an increase in the temperature of the structural phase transition (T
C). For example, sodium nitrite embedded in a porous glass matrix with nanosize through pores demonstrates shifting of ferroelectric (FE) phase transition in the low-temperature direction [
3]. A similar situation occurs in the case of ammonium hydrogen sulfate, in which the temperatures of the both ferroelectric and antiferroelectric (AFE) phase transitions are decreased [
4,
5,
6,
7].
Alongside this, potassium dihydrogen phosphate [
8,
9,
10,
11,
12,
13], Rochelle salt [
14,
15], and triglycine sulfate [
14,
16,
17,
18,
19], in contrast, are characterized by an increase of ferroelectric phase transition temperature under “restricted geometry” conditions. However the data on the magnitude of T
C shifting as a function of the matrix mean that the pore diameter varied significantly for different researchers.
It was suggested [
8] that the rise of the Curie temperature in KDP inclusions is the result of the tensile deformations which appear due to differences in the temperature coefficients of linear expansion (TKL) of the embedded substance and matrix material. This assumption was confirmed by a comparative analysis of shifting of the phase transition temperature in isomorphous crystals of the ferroelectric potassium dihydrogen phosphate and the antiferroelectric ammonium dihydrogen phosphate with the different baric coefficients (γ = dT
C/dP, where P is a hydrostatic pressure) and the different temperature coefficients of linear expansion [
9].
However the expansion of the temperature range, within which the ferroelectric phase exists in the embedded nanoparticles of triglycine sulfate and Rochelle salt [
14,
15,
16,
17,
18,
19] cannot be related to the action of mechanical stresses. The phase transition temperatures of these salts do not differ so much from the temperature at which they were introduced into the matrix.
These circumstances indicate a strong “non mechanical” interaction between the matrix and the embedded material. However, mechanisms of such interaction have not yet been clarified.
This work is a comparative study of the dielectric properties of composite materials obtained by embedding the related H – bonded crystals potassium dihydrogen phosphate (KH2PO4 - KDP), cesium dihydrophosphate (CsH2PO4 - CDP), and ammonium dihydrogen phosphate (NH4H2PO4 - ADP) in porous matrixes with different average pore diameters. The KDP, CDP, and ADP crystals have a similar chemical composition, but differ in the baric coefficients and in the temperature coefficients of linear expansion.
2. Samples Preparation and Validation
The KDP, CDP, and ADP polycrystalline samples were prepared by compacting appropriate powders in the form of disk with a diameter of 10 and a thickness of 1 mm. The grain sizes were within 0.1–0.4 mm.
Composite materials were fabricated by embedding corresponded salts into glass matrices with average diameter of through pores equal to approximately 320, 46, and 7 nm (the pore sizes were determined by the SEM method).
The samples of the glass matrices were cut off in the form of flat plates of about 10 × 10 × 0.7 mm. The porosity of used glasses matrix, which ranged from 40% to 55%, was determined by the relative mass decrement method during the preparation of the samples. (The technology of porous glasses and methods of their certification are presented in references [
7,
20]).
Before filling, the glass matrices were annealed and weighed. Then they were immersed in a saturated aqueous solution with corresponding quantities of salt at a temperature of 90 °C for four hours.
After the filling, the samples were pulled out and dried at room temperature for 24 h. Thereafter, they were subjected to 6-h thermal treatment at a temperature of 150 °C. Then, the surface layer was mechanically removed. The share of the substance embedded in the matrix had around a 0.2 volume fraction.
The samples were labeled as follows: KDP-PGX, ADP-PGX, and CDP-PGX, where X was the average matrix pore diameter in nanometers.
The X-ray diffraction patterns for the PG-320 glass matrix, the original KDP powder, and the KDP-PG320 composite, obtained at room temperature using a x-ray diffractometer (Model-DRON-3) with Bragg–Brentano geometry, Cu K
α1-radiation, are shown in
Figure 1 a–c correspondingly.
One can see the gallo in the case of a glass matrix, which is characteristic for amorphous materials. Lines of intensity (I) in
Figure 1b correspond to the reference data for potassium dihydrogen phosphate. The X-ray pattern obtained for the composite KDP-PG320 (
Figure 1c) is a superposition of the two previous I(2θ) dependences. Thus, the crystal structure of the bulk and embedded material is identical, and there are no extraneous phases.
For experiments, a conductive silver-containing paste was applied on large facets of the samples, as electrodes. The samples were placed in a cryostat, containing an evacuated measuring cell. Before each experiment, the samples were annealed about a one hour at a temperature of ≈ 100 °C to remove the adsorbed water. Its content in the sample was controlled qualitatively by measuring dielectric permittivity and the dielectric losses tangent (tgδ) during thermal treatment. The cessation of ε’ and tgδ variations with time indicated the removal of most of the water.
During the experiment, the temperature in the cryostat varied within 85–300 K, and its error measurement did not exceed ± 0.5 K. Measurements of the dielectric permittivity (ε’) were carried out in a slow heating mode with a rate of about 0.5 K/min using the E7-20 immittance meter at the frequency f = 1 kHz.
3. Experimental Results
The temperature dependences of dielectric permittivity of ammonium dihydrogen phosphate, potassium dihydrogen phosphate, and cesium dihydrogen phosphate, as well as composites based on them, are shown in
Figure 2,
Figure 3 and
Figure 4. One can see that dielectric permittivity of the investigated polycrystalline samples ADP, KDP, and CDP is reduced in comparison with known references data [
21]. This is due to significant porosity of the prepared samples, which, however, does not affect the corresponding Curie temperatures.
In the case of ammonium dihydrogen phosphate, an abrupt change of ε’ is observed near the T
C ≈ 150 K due to the first-order AFE phase transition (
Figure 2a).
For ADP-PG composites, as well as for bulk ADP, antiferroelectric phase transition was accompanied by step-like anomalies of the dielectric response (
Figure 2b–d). These anomalies were visibly diffused as the average pore diameter of the matrixes is decreased. Apparently, this is a consequence of the phase transition diffuseness in embedded ADP particles.
It is known [
21] that the appearance of an order parameter at FE or AFE phase transition leads to a decrease in dielectric permittivity. We used this to determinate T
C, which was found as the temperature below for which a decrease in ε’ was observed.
The temperature of AFE phase transition in ADP nano-particles, included in the composite, is slightly higher than T
C of a bulk ammonium dihydrogen phosphate (
Table 1).
The ε’(T) dependencies for potassium dihydrogen phosphate and cesium dihydrogen phosphate polycrystalline samples, as well as composites based on them, are shown in
Figure 3 and
Figure 4, respectively.
In the vicinity of the FE phase transition in the bulk KDP (T
C ≈ 122 K) and CDP (T
C ≈ 154 K) samples, distinct maxima of the dielectric permittivity were observed. Their position in the temperature axis is in good agreement with references [
21,
22].
In the case of composite materials, ε’ maxima were smeared and shifted in the high-temperature direction with decreasing of the mean pore diameter of the glass matrices (
Table 1). Obviously this is a result of increases in temperature and a smearing of the FE phase transitions in particles of embedded salts. The KDP-PG46 composite is the one exception here, because it has the same phase transition temperature as the bulk KDP.
Comparing the values of TC shifting for ADP, KDP, and CDP crystals under restricted geometry conditions, we can state that there was an increase in the phase transition temperature in embedded particles. The largest shifting of TC was found for the ferroelectric KDP, and the smallest one was found for the antiferroelectric ADP.
4. Discussion
An influence of the pore size on the phase transition temperature in KDP and ADP particles, embedded in porous matrices, has already been discussed [
8,
9,
10,
11,
13]. However, experimental results obtained by different authors contain some contradictions. For example, the experimental study of the dielectric properties of KDP in porous structures with pore dimensions D = 7–100 nm demonstrate that a decrease in the pore diameter leads to a monotonous increase in the T
C. At the same time, the nonmonotonic behavior of the T
C(D) dependence was revealed in references [
10,
11] on the basis of dielectric and dilatometric measurement results.
4.1. “Pressure Effect”
The results of our experiments showed an increase in the T
C in ADP, CDP, and KDP particles embedded in porous glass matrices with D ≈ 320 nm (
Table 1).
Following the authors of reference [
8], we assumed that the observed increase in the phase transition temperatures is a consequence of mechanical tensile stresses in the embedded particles. These stresses are the result of differences between the temperature expansion coefficients (TKL) of the matrix and the embedded material, and arise during the cooling of the sample from the temperature, at which the filling matrices was carried out (≈ 370 K), to the T
C. We will call this phenomenon the “pressure effect”.
The temperature coefficient of linear expansion for the quartz glass was estimated by the value α
SiO2 < 10
−6 K
−1 [
23], whereas the coefficients α
i for ADP in the temperature range of 203–407 K were: α
1 = α
2 ≈ 34.0–39.3 × 10
−6 K
−1 and α
3 ≈ 1.9–5.3 × 10
−6 K
−1 [
24]. For KDP crystal in the temperature range of 123–363 K, the temperature coefficients were: α
1 = α
2 ≈ 20–26.6 × 10
−6 K
−1 and α
3 ≈ 34.3–44.6 × 10
−6 K
−1 [
24]. And for cesium dihydrogen phosphate, the same values are: α
1 ≈ 31 × 10
−6, α
2 ≈ 84 × 10
−6, and α
3 ≈ −25 × 10
−6 K
−1 [
25].
Silicon dioxide has a lower TKL value in comparison with ADP, CDP, and KDP, so the tensile stresses arise in embedded particles during cooling of the composite sample, and due to the compressive stresses in the matrix.
If pores of the matrix are partially filled, and the thickness of the “layer” of embedded substance is much smaller than the thickness of the interchannel space of the porous matrix, then the deformation of the interchannel space can be neglected. Under these conditions, the effective elastic stresses in particles of the embedded material may be estimated using the following expression:
where α
av = (α
1 + α
2 + α
3)/3 is a mean TKL value, c
av = (c
11 + c
22 + c
33)/3 is a mean value of elastic modulus of the embedded substance, and ΔT is the difference between the embedding temperature and the T
C.
We assume that the tensile stress σ
ef affects on the phase transition in the same way as the hydrostatic pressure. Taking this into account, we can estimate the value of the phase transition shifting as:
where γ ≈ −8.5 × 10
−8, −4.5 × 10
−8, and −3.4 × 10
−8 K/Pa for CDP, KDP, and ADP particles, respectively [
21,
22].
Using the reference data for ADP, KDP [
24], and CDP [
26,
27] crystals, one can find the corresponding average values of the elastic modulus c
av = 5.5 × 10
10, 7 × 10
10 and 4 × 10
10 Pa, the mean stress values σ
ef ≈ 2.06 × 10
8; 2.2 × 10
8 and 1.73 × 10
8 Pa, and to estimate the shifting values ΔT
C ≈ 7, 10, and 14.7 K for ADP, KDP, and CDP, respectively.
Thus, the obtained ΔTC estimates are of the same order of magnitude as the values obtained from the experiment.
However, we should note that if thicknesses of embedded matter layers and the interporous space are comparable, then there will be noticeable deformations in the glass matrix, which will lead to a decrease in the elastic stresses in particles of the embedded material. This is most likely to occur for composites with a high content of embedded matter. The results of dilatometric measurements for composites ADP-PG160 and ADP-PG23 [
11] have revealed a significant deformation of glass matrix in the AFE phase transition region, which illustrates this assumption.
Among the porous matrices used in these experiments, the PG46 matrices are markedly more brittle. This indicates thin interporous partitions in them. It is clear that that this circumstance led to a very small T
C shifting in the KDP-PG46 composite in comparison with the bulk material (
Table 1).
4.2. Transversal Tunneling Ising Model
A strong “pressure effect” in crystals of potassium dihydrogen phosphate family is related to the mechanism of the ferroelectric phase transition, initiated by the ordering of protons in two-minimal O-H … O hydrogen bonds [
21]. Proton tunneling prevents their ordering and therefore significantly affects the temperature of the phase transition.
The dependence of T
C on tunneling frequency Ω, which is related to proton tunneling, can be described within the framework of Ising model in a transverse tunneling field, developed for a bulk KDP crystal [
21,
28]. This model can also be used to describe the phase transitions in CDP [
29] and ADP [
30].
According to references [
21,
28],
where J
0 is the energy of interaction of the i-th dipole, while with all other dipoles
where J(R,R’) is the energy of interaction of the dipole, located in the R-position, with another one, located in the R’-position.
Mechanical tensile stresses, acting on the embedded particles, lead to an elongation of the hydrogen bonds and, consequently, to a decrease in the tunneling frequency Ω, that leads to an increase in the TC.
Alongside this, the relation (3) makes it possible to explain the decrease in TC with a decrease in the size of embedded particles. For ultrafine particles, the sum (4) will decrease because of the limited number of interacting dipoles, which will lead to a decrease in J0 and TC.
4.3. Influence of the Piezoelectric Effect on the TC
Let us discuss the possible impact of the piezoelectric effect on the T
C. Potassium dihydrogen phosphate and ammonium dihydrogen phosphate crystals are piezoelectrics in paraelectric phase with non-zero g
36 and g
14 piezoelectric coefficients, where g
36 >> g
14 [
24,
31,
32].
The piezoelectric effect can lead to the appearance of significant electric bias fields E3 at the temperatures T ≥ TC due to the appearance of the shear component of mechanical stresses σ6 for inhomogeneous deformations of embedded particles.
We estimated the field strength E
3 = g
36σ
6 [
31,
32], and the possible increase in the T
C as a result of its action [
21]:
Taking into account that coefficient A ≈ 5 × 10
−4 K⋅cm/kV for the KDP [
33] and assuming that σ
6 ≈ 0.1σ
ef, we can find ΔT
EC ≈ 13 K.
Thus, the calculated shifting of the TC as a result of the electric field action for KDP, caused by the direct piezoelectric effect, is of the same order of magnitude as the experimentally observed value of ΔTC.
In the estimation of electric field E
3, we did not consider its decreasing due to the polarization of the glass matrix, as well as a charge leak owing to electrical conductivity of the material, although it is quite small at the temperatures near T
C [
34]. Therefore, it is obvious that the real values of the field E
3 can be less than the above estimates.
Phenomenological theory predicts the decreasing of antiferroelectric phase transition temperature under a bias electric field E [
21]. However the experimental data regarding the electric field E effect on the T
C in the ADP crystal are absent in references, probably due to the destruction of the bulk ADP samples below T
C caused by strong elastic stresses [
21]. At the same time, the experiments with antiferroelectric ADP - KDP solid solutions did not reveal affects of electric field on T
C, at least at E < 10 kV/cm [
35].
The CDP crystal has a piezoelectric effect in the ferroelectric phase only [
36]. Therefore, one can expect that the electric field acting in the CDP particles due to the piezoelectric effect will be significantly less than in KDP or ADP.
5. Conclusions
The experimental results have shown that the phase transition temperatures in ADP, KDP, and CDP particles embedded into porous glass matrices with a nanometer pore diameter (D) are higher than in bulk samples. However, the observed relationship between Tc and D is ambiguous. In our opinion, this is due D not being the only factor affecting Tc.
An analysis of the factors that affect the TC in materials under study speaks in favor of two mechanisms providing the largest contributions.
The first one is caused by dependence of the TC on tensile strains in ferroelectric particles (“pressure effect”) due to different TKL of the embedded material and the matrix. Its microscopic mechanism is the decreasing of proton tunneling due to the rise in the distances between the potential minima in double-well O-H … O hydrogen bonds.
The greatest increase in the TC due to the “pressure effect” should be expected for CDP particles, but the experimental value ΔTC for CDP particles is smaller than the same for KDP, under equal experimental conditions.
Apparently, this is due to the bias electric fields induced by elastic stresses owing to the piezoelectric effect in the paraelectric phase of potassium dihydrogen phosphate. The CDP crystal has no piezoelectric effect above the TC, so the phase transition temperature shifting in CDP particles is smaller than in the KDP.
The electric field induced by the piezoelectric effect in ADP particles, according to the above estimations, is almost double that in KDP. However it has no effect on the temperature of antiferroelectric phase transition.
Both of the above mechanisms are due to mechanical stresses, depending on both the difference between the TKL of embedded particles and the matrix, and the ratio of the layer thicknesses of the embedded substance and the interporous space. So, they depend on the filling of the matrix pores.
Together with the piezoelectric effect, another source of internal bias fields could be the chemical interactions of protons of the embedded salts and the oxygen atoms of the SiO
2 matrix. In particular, the existence of chemical interaction was revealed by analysis of IR spectra of nanocomposites based on particles of TGS and SiO
2 [
37].