1. Introduction
SrTiO
3 (STO), a cubic perovskite oxide, was synthesized in the early 50’th [
1] and since then was intensively investigated. A breakthrough discovery was made by K. A. Müller in 1958 [
2], when he unambiguously established the structural phase transition from cubic to tetragonal at around 105 K. Later he identified the tetragonal rotation angle of the oxygen octahedra as the order parameter [
3]. Shortly after the discovery of the 105 K structural phase transition, Cochran predicted the occurrence of zone center transverse soft optic modes in perovskite oxides like SrTiO
3 and BaTiO
3 at the Γ-point [
4]. Their existence was confirmed experimentally in the following years [
5]. Related to the rotational instability in STO is a zone boundary transverse acoustic mode at the R-point which softens with decreasing temperature to become unstable at
TS=105 K [
6]. Below
TS the mode folds back to the zone center and splits into three optic modes, two doubly degenerate and a single mode, the frequencies of which all follow the Curie Weiss law in temperature. We will show here that this point turns out to be important for the understanding of thermal conductivity TC and the birefringence Δn. Two important experiments were carried out already in the early years of research on STO, namely inelastic neutron scattering experiments where an anomalous acoustic mode dispersion was seen at small momenta
q which was interpreted as originating from mode-mode coupling [
7]. The other observation was that the dielectric permittivity did not peak at low temperatures but saturates to a plateau. This was interpreted as suppression of ferroelectricity by quantum fluctuations and was termed
quantum paraelectricity [
8].
While quantum paraelectricity attained enormous attention after its introduction, the acoustic phonon mode anomaly at small momentum was rarely addressed. It has, however, important consequences for the elastic properties and the origin of precursor phenomena as well as for dynamically driven elastically distorted domains appearing on specific length scales. We focus on this aspect by using results from theory obtained within the polarizability model and experiments from birefringence and elastic measurements. In addition, new measurements of the thermal conductivity to ultra-low temperatures are interpreted in terms of quantum effects and mode-mode coupling theory. Besides the importance of elastically distorted local regions newly obtained data from the formation of polar clusters are emphasized. Our results imply that polar and elastically distorted dynamical domains coexist at temperatures below the structural phase transition resulting from anharmonic mode-mode coupling [
9]. A comparison with experimental data is made.
2. Theory
The dispersion relations for the lowest optic and acoustic modes
along (100) have been obtained within the polarizability model [
10,
11,
12] as a function of temperature
T and the momentum
q. At high temperatures the optic and acoustic modes are well separated, and only negligible coupling or crossing occurs. Around and below the structural phase transition a tiny dip at small momentum
q appears in the acoustic mode which becomes increasingly pronounced with decreasing temperature. Simultaneously the dip shifts to smaller momentum. A special feature, which has been rarely addressed so far, is the fact that
softens at the zone boundary caused by optic-acoustic mode-mode coupling. This softening remains incomplete down to the lowest temperature of 4 K [
13]. The anomalies at small momentum can be easily identified by taking the momentum derivatives of the two modes (
Figure 1).
The
q-dependent mode-mode coupling anomaly can be associated with the deviation from linearity which appears as the inversion point in
and the maximum, i.e. the saddle point in
. This choice is justified from the fact that the acoustic mode derivative smoothly decreases with
q whereas the optic mode derivative exhibits no saddle point in the harmonic case. The corresponding momenta
qc define the length scales of elastic and polar clusters which appear simultaneously at the same temperatures. Since the coupling stems from the off-diagonal matrix elements in the dynamical matrix, an analytic expression for both derivatives are rather difficult to obtain but they can be derived numerically (see
Figure 1). Whereas elastic clusters appear on a rather large length scale as compared to polar clusters, both length scales start to diverge with decreasing temperature. A true divergence takes place in the elastically distorted areas whereas the polar clusters adopt a similar behavior as the dielectric permittivity, namely they saturate in the quantum fluctuation dominated region below ≈30 K. The calculated size of the elastically distorted areas is similar to those reported recently by inelastic neutron scattering, although by a factor of 2 smaller. This difference might be related to the too coarse
q sampling in the calculations. Above
TS precursor dynamics of the structural phase transition occurs as already reported previously for STO and several other perovskite oxides. If the linear regime for
T >
TS is subtracted from the data (inset to
Figure 2), the precursors become more apparent as indicated by the shaded area in the inset to
Figure 2. Experimentally, strong elastic precursor softening was observed at temperatures as high as ~125 K i.e. ~20 K above the actual transition point while analytical data fits show even higher precursor temperatures [
14]. Similar conclusions were reached by birefringence measurements [
15]. Theoretically, precursors have been predicted to be universal in perovskite oxides [
16]. Evidence for polar nano-domain formation has been obtained by resonant piezoelectric spectroscopy [
17], which has been interpreted in terms of ferroelastic twin walls which become polar at low temperatures.
The inset to
Figure 2 also highlights what can be expected for the birefringence data and elastic effects, namely in a rather small temperature region below
TS deviations from simple, one order parameter mean-field behavior should occur, whereas approaching the quantum fluctuation dominated regime, nonlinear elastic anomalies with the simultaneous appearance of polarization clusters are observed [
18] which is consistent with the present results (
Figure 2). It is important to emphasize that deviations from mean-field results have already been postulated by Müller and Berlinger [
19] who argued that for temperatures very close to
TS (
t =( TS –
T)/
TS > 0.9) the order parameter can be described by a critical exponent β≈0.33(2). Later, Salje et al. [
20] showed that criticality occurs only in a very small interval near
Ts while the temperature evolution of the excess entropy is fully compatible with a mean-field, near tricritical Landau potential. A weak singularity was seen for a maximum temperature interval of 2K. Below 103 K (0.981
TS) the order parameter follows a mean field behavior which is very well described by a 2-4-6 Landau behavior. Birefringence data by Geday and Glazer [
21] concurred and did not find deviations from a mean field behavior
A closer look into the cluster formation below
TS and especially in the quantum regime reveals that three novel regions in the cluster size can be identified (
Figure 3) which can be related to the twin boundary formation below
TS. For temperatures 60 K <
T <
TS a small but smooth increase in both cluster states takes place. Below 60 K and more pronounced below 40 K especially the elastic properties are strongly affected, and a rapid increase of the elastically deformed clusters takes place. The polar regions exhibit a similar but much less pronounced behavior. Below some 25 K to 20 K a steady state is reached characterized by scattering in the size and distribution around an average value of the polar elastic domains.
In the quantum paraelectric state a “quantum domain glass” at
T < 40 K shows intense relaxation and temperature hysteresis of its nanostructure. This includes a high correlation between domains so that domains float collectively in a complex, smooth landscape with long relaxation times. In the “quantum domain solid” state below 25 K this correlation becomes very strong. Approaching zero Kelvin leads to a predominance of large coherently moving clusters [
13].
Motivated by recent inelastic neutron scattering data [
22], the soft modes observed in STO have been reinvestigated. Special interest [
22] was devoted to the transverse optic/acoustic mode-mode coupling region which followed to the quantum fluctuation dominated regime. The conclusions were that the ground state of STO is formed by a hybridized optical-acoustic phonon mode where quantum fluctuations are accompanied by fluctuating domains of mesoscopic length. In addition, the softening of the M and Γ point modes was followed over a large temperature range with the observation that their temperature dependence perfectly overlap on the whole temperature range when shifted by 9 meV. The origin of this observation has not been addressed yet. As is shown in
Figure 4, this fact finds a natural explanation since both modes are dominated by the polarizability coordinate and consequently adopt analogous temperature behavior. In addition, also the soft acoustic zone boundary mode related to the TO Γ-point mode follows this temperature dependence, such that all three modes of STO overlap if shifted by a mode-specific energy. At very low energies the latter two modes saturate as a consequence of the onset of quantum fluctuations. A distinct difference between the three modes lies in their high temperature properties where deviations from mean-field behavior and saturation start near
TS for the transverse optic mode. The R-point acoustic mode shows first deviations from mean field behavior below ca. 300K. This temperature evolution has been predicted in Refs. 9 – 11 and identified for the first time in the antimony-sulfo-iodide SbSI [
23] and later also in STO [
24].
3. Thermal Conductivity
We now turn to the discussion of the thermal conductivity (TC). Experimentally the TC was measured with a conventional two-thermometer one heater arrangement using a commercial system (Quantum Design, PPMS) Our data (see
Figure 5) which cover the temperature range from room temperature down to 2 K are similar to the results reported by Jaoui,
et al. [
25], Martelli, et al. [
26] and Steigmeier [
27]
The TC is a function of the specific heat c
v, the phonon mean free path and the phonon group velocity. All paramters are specifically temperature dependent. The transverse optic soft mode and the related acoustic mode dominate TC at low temperatures since the crossing wave vector avoidance rule produces the dip in the acoustic mode dispersion and defines the mean-free path. The group velocities for these two modes are readily obtained from their dispersion where the one related to
is much larger than for the optic mode. Acoordingly, the contribution from the optic mode can be neglected for this discussion. A strong increase of TC [
28,
29,
30] is seen near 40 K. Around
TS a rather broad anomaly (see inset to
Figure 5) in TC appears which signals the phase transition and possibly the nucleation of domain structures in the ferroelastic phase. In
Figure 6 the frequencies of the TA and TO modes are plotted as function of temperature where the shaded regions labeled I, II, III denote the quantum domain solid, the quantum domain glass regime, and the region around the phase transition temperature where precursors and deviations from mean-field behavior are highlighted. Note, that above
TS ≈ 105 K the calculated TC is in good agreement with the measured data [
12]. For
T >
TS substantial corrections are needed which relate to the R-point acoustic soft mode.
Above
TS the R-point acoustic mode softens substantially and completely dominates the TC. This can be shown by replotting the inverse of TC (
Figure 5) as a function of temperature (
Figure 7). Above
TS TC follows a Curie-Weiss law like the soft mode and shows deviations from it in the precursor region. Here the effect of other modes is apparently less relevant for the heat transport which is dominated by the acoustic R-point mode. Below
TS this mode converts to optic modes and a more complex temperature evolution of TC is seen.
From these results we conclude that novel temperature scales have been observed which are needed to explain experimental observations Especially, we have shown that the quantum fluctuation dominated region has profound effects on TC followed by a crossover state where the lowest transverse optic and the related acoustic modes are the most important ingredients for the understanding of the dynamics and the TC. Above TS the zone boundary soft acoustic mode at the R-point gains importance and is the only relevant contributing to TC.
4. Birefringence
We now discuss birefringence measurements at low temperatures. These experiments have been performed analogously to the procedure reported in [
20]. They are now analyzed in somewhat deeper detail and with emphasis on the critical regions discussed above. The results are displayed in
Figure 8 for temperatures 80 K <
T < 180 K. Below
TS the data can be fitted linearly in accordance with simple Landau theory. Close to
TS weak deviations are observable which cover a small interval around
TS, small compared with results from measurements of the excess entropy and the structural order parameter
Q [
20]. Note that the birefringence data reveal small variations in
TS with variable sample thickness.
Nominally, Δn should be zero above the phase transition temperature small but finite signals are observed at
T >>
TS defining a wide precursor regime. This behavior is analogous to elastic anomalies observed in [
17]. The data above
TS can be described in a related manner as the TC, namely as a function of
t = (1 -
T/
TS) which is shown in
Figure 8.
Figure 9.
Double logarithmic plot of the birefringence Δn as function of the reduced temperature t. If it is assumed that below the continuous transition the critical exponent can be computed from the relationship Δn(T) ~ (TS – T)2β we find β to be almost 0.5 below 99 K.
Figure 9.
Double logarithmic plot of the birefringence Δn as function of the reduced temperature t. If it is assumed that below the continuous transition the critical exponent can be computed from the relationship Δn(T) ~ (TS – T)2β we find β to be almost 0.5 below 99 K.
An inversion point in the several data around identical temperatures. TC in
Figure 7 follows a linear increase at higher temperatures which can be fitted by a Curie-Weiss law with a critical exponent
β≈0.498, which is very close to the expected mean-field behavior. The question arises as to why the birefringence should exhibit the same temperature dependence as the zone boundary transverse acoustic mode? Whereas the TC dependence can be well explained by being dominated by this mode, this is not necessarily true for Δn. Here the softening of
might lead to elastically distorted regions which adopt lower symmetry than the bulk and correspondingly allow for the observed weak birefringence. On the other hand, defects might also be a possible source of symmetry lowering, at present it is yet not clear which effect dominates.