# Magnetic and Magnetoelectric Properties of Rare Earth Molybdates

Authors: B. K. Ponomarev; A. Zhukov
Journal: Physics Research International
 (2012)

---

## Abstract

We present results on ferroelectric, magnetic, magneto-optical properties and magnetoelectric effect of rare earth molybdates (gadolinium molybdate, GMO, and terbium molybdate, TMO, and samarium molybdate, SMO), belonging to a new type of ferroelectrics predicted by Levanyuk and Sannikov. While cooling the tetragonalβ-phase becomes unstable with respect to two degenerate modes of lattice vibrations. The β-β′ transition is induced by this instability. The spontaneous polarization appears as a by-product of the lattice transformation. The electric order in TMO is of antiferroelectric type. Ferroelectric and ferroelastic GMO and TMO at room temperature are paramagnets. At low temperatures GMO and TMO are antiferromagnetic with the Neel temperatures TN=0.3 K (GMO) and TN=0.45 K (TMO). TMO shows the spontaneous destruction at 40 kOe magnetic field. Temperature and field dependences of the magnetization in TMO are well described by the magnetism theory of singlets at 4.2 K ≤ T ≤ 30 K. The magnetoelectric effect in SMO, GMO and TMO, the anisotropy of magnetoelectric effect in TMO at T = (1.8–4.2) K, the Zeeman effect in TMO, the inversion of the electric polarization induced by the laser beam are discussed. The correlation between the magnetic moment of rare earth ion and the magnetoelectric effect value is predicted. The giant fluctuations of the acoustic resonance peak intensity near the Curie point are observed.


---

## Body

## 1. Introduction

At the end of 1960s, there was a great interest in the rare earth molybdates familyR2(MoO4)3 (RMO) (R=Pr,Nd,Sm,Eu,Gd,Tb and Dy) because these compounds exhibited the phenomena of ferroelectricity and ferroelasticity [1].RMO fromPr2(MoO4)3 (PMO) to TMO crystallizes to the tetragonal β-structure with a space group P4̅2m (point symmetry group is 4̅2m). Their melting points are 1045°C for PMO and 1172°C for TMO. Dy2(MoO4)3 (DMO) crystallizes to cubic γ-phase at 1222°C and transforms to β-structure at 1030°C. While cooling the family undergoes a transformation from the tetragonal β-phase to the monoclinic α-Phase. The temperatures of the β-α transformation are 987°C for PMO and 805°C for DMO. However, the transformation is sluggish and, therefore, the high-temperature phase can be quenched in. If the thermodynamically metastable β-phases of RMO are further cooled, they undergo the second transformation leading to the lower symmetry ferroelastic-ferroelectric orthorhombic Pba2  β′-structures (point symmetry group mm2). They are also thermodynamically metastable. The temperatures of β-β′ phase transitions are 235°C for PMO and 145°C for DMO.The single crystal samples of RMO are transparent in the visible light.Bothβ and β′ phases are piezoelectric. The β′ phase is ferroelectric. The spontaneous electric polarization values are of ≈2 × 10−7 Coul/cm2. The dielectric permeability peak at the Curie temperature is of ≈10 for GMO.The ferroelastic properties are well exhibited. In a single-domain sample of (001) cut, the switching of the ferroelectric domain can be observed visually in polarized light at compressing the sample along the [010] axis.Orthorhombicβ′-phases of rare earth molybdates are paramagnetic down to temperatures below l K.In RMO with R = Sm, Gd, Tb and in mixed molybdates DyGd(MoO4)3 (DGMO), TbGd(MoO4)3 (TGMO) the nonlinear magnetoelectric effect was observed experimentally. In TMO and in TGMO, the switching of ferroelectric domains by the magnetic field was observed experimentally. These are the first substances in which the possibility was found to switch the ferroelectric domains by the magnetic field.The longitudinal Zeeman effect was investigated in TMO in the magnetic fieldH≤20T along the [001] and [110] axes at T=4.2K and 1.7K. The experimental field dependences of the wave numbers of the absorption peaks νi(H) were obtained. The analysis of the dependences shows that the excited multiplet 5D4 of Tb3+ ion in the crystal field of TMO can be considered as consisting from five singlets and two quasidoublets. The field dependences of the energy levels Ei(H) were obtained from the experimental data in the magnetic field along the [001] axis for the 5D4 multiplet. The experimental dependences Ei(H) were described well by the theory of the magnetism of singlets.The dependences of the photoinduced voltage in TMO upon the timet of the illumination by a laser beam were measured at room temperature at 0.1s≤t≤4500s. The power of the laser radiation was 0.4W. The distribution of the intensity along the sample was varied from highly inhomogeneous to homogeneous. It was found that the photoinduced voltages of two types with the opposite signs appear due to the inhomogeneous illumination. The model was proposed to explain the appearance of the photoinduced voltage due to the inhomogeneous illumination.The qualitative explanation of the mechanisms of these two effects was given. To make quantitative estimates, it is necessary to measure the dependences of the effect upon the temperature of the sample and upon the intensity of the light.The temperature dependence of the intensity of the main acoustic resonant peak in GMO was measured at temperatures from 22°C to 165°C. This temperature range includes the Curie point (TC=159°C). The low-frequency fluctuations of the amplitude of the resonant current in GMO were observed near the Curie point. The corresponding variations of the measured voltage were 0.1 Volt. The typical times of the fluctuations were (10–100)  sec. The relative values of the corresponding fluctuations of the piezoacoustic impedance in GMO are two orders of magnitude larger than the values of the fluctuations of the physical parameters in the solids that were known before. Such fluctuations of the physical properties of the solids were observed visually never before.In [2], the point symmetry group C2 for the nearest surrounding of Gd3+ ions in GMO was established experimentally. The calculated value of the spontaneous polarization, assuming point charges Gd3+ and (MoO4)2-, was 260 · 10−9 Coul/cm2. This result can be considered as being in a satisfactory agreement with the experimental value 200·10-9Coul/cm2. All rare earth molybdates are isomorphic to GMO. So, the rare earth ions in RMO crystal lattice are assumed to be trivalent and the point symmetry group of the local crystal field at R3+ sites is assumed to be C2. The experimental data discussed below relates mostly to Sm2(MoO4)3, Gd2(MoO4)3, and Tb2(MoO4)3. Table 1 contains the quantum numbers for corresponding free R3+ ions in the ground state.Table 1The ionThe ground  stateThe spin,SThe angular Momentum,LThe total Momentum,JSm3+6H5/25/255/2Gd3+8S7/27/207/2Tb3+7F6336It is worth mentioning that starting from 90-th magnetic shape-memory alloys (MSMAs) attract special interest owing to significant magnetic-field-induced strain (MFIS), also referred as the “magnetic shape-memory effect” originated from coupling between magnetic and structural ordering. Such strain arises through the magnetic-field-induced motion of twin boundaries [3]. Since the magnetic shape-memory effect is useful for actuation purposes, the inverse effect may be utilized for sensing and energy harvesting applications [4]. Consequently, the direct and inverse magnetic shape-memory effects cause magnetic field-induced superelasticity, which is the magnetically induced recovery of a large mechanically induced deformation [4]. It is worth mentioning that, in fact, magnetic-field-induced strain (MFIS) effect has been previously observed and described in the same terms at the beginning of 90-th in ferroelectrics (rare-earth molybdates) [5]. Particularly, terbium molybdate shows so strong field-induced strain, that the spontaneous destruction in a constant magnetic field of 40kOe along the [100] axis at T=1.5K has been observed [5], but at that time, most attention has been paid to magnetic-field-induced electric polarization and magnetoelectric effects.Latter family of magnetic shape memory alloys (MSMA) has been introduced, where the coupling between magnetic and structural ordering in conjunction with the magnetic and structural transformations giving rise to various functional properties: the magnetocaloric effect, the magnetic-field-induced martensitic transformation, as well as its reverse transformation, giant magnetoresistance, and electric polarization [4].In this paper, we will present the review on ferroelectric, magnetic, electric, and structural properties of the rare earth molybdates familyR2(MoO4)3 (RMO) (R=Pr,Nd,Sm,Eu,Gd,Tb, and Dy).


## 2. Ferroelectric Properties

Here, only a short description of ferroelectric properties of RMO will be given. The detailed review of this question is available in [1].RMO belongs to a new type of ferroelectrics (improper ferroelectrics) theoretically predicted by Levanyuk and Sannikov [6].While cooling the tetragonalβ-phase becomes unstable with respect to two degenerate modes (soft modes) of lattice vibrations, this instability induces the displacive phase transition of the first order to the ferroelastic-ferroelectric β′-GMO-type structure. The β-β′ transition temperature is also the Curie temperature of the orthorhombic β′-RMO. The spontaneous polarization appears at the transition as a by-product of the lattice transformation and cannot be considered as the order parameter. This conclusion was obtained for GMO by Petzelt and Dvořak [7] using group analysis. They confirmed this result experimentally by far-infrared reflectivity and transmission measurements in GMO.Analogous result was obtained for TMO from neutron-scattering experiments [8, 9]. In [6], from the examination of the soft mode, it was shown that the primary-order parameter, which necessarily shows large fluctuations, near TC, have no macroscopic polarization but have instead antipolar. The freezing in of this antiferroelectric static displacement couples to the shear strain which in turn produces the electric polarization by piezoelectric coupling.The displacive mechanism ofβ-β′ transition is valid for all the members of RMO family because they are isomorphic to GMO and TMO.The crystal structure ofβ′-GMO was solved in [2, 10, 11]. Above the Curie temperature, the structure is tetragonal with space group P4̅21m and point group 4̅21m. Below the Curie temperature, it is orthorhombic with space group Pba2 and point group mm2. The reduction in point group symmetry from 4̅21m to mm2 allows the formation of twin domains. They can be distinguished in polarized light. After cooling through the Curie point, the RMO single crystals always are in a polydomain state. When cooled through the Curie temperature, the equivalent [110] and [11̅0] axes of the tetragonal phase become the [100] and [010] axes of the orthorhombic phase. These orthorhombic axes are unequal. The b parameter along the [010] axis is larger than the a parameter along the [100] axis: (b-a)/a≈3×10-3. So, the crystal changes its shape during the transition. The axes [100] and [010] of the tetragonal phase become the axes [11̅0] and [110] of the orthorhombic phase and the right angle between them changes by the value of ≈10′. This is the shear angle of the transition. It is determined by the relation (b-a)/a. The relative change of the [001] axis parameter Δc/c during the transition is smaller than 10-4, but the structure of the [001] axis is changed drastically. It is not polar in the tetragonal phase and it becomes polar in the orthorhombic phase. This is the origin of the ferroelectric order.The applied electric field along the [001] axis can reverse the [001] axis. This is accompanied by the mutual interchange of the orthorhombic [100] and [010] axes.Due to the difference betweenb and a parameters, the mechanical compression along the [010] axis switches the orientational state. This is the origin of the ferroelasticity. A material is said to be ferroelastic if it has two or more stable orientational states in the absence of the external mechanical stress and if it can be reproducibly transformed from one to another of these states by the application of mechanical stress. Ferroelasticity was established in RMO by Aizu [12].Ferroelasticity can be used to transform a multidomain sample into a single domain state. In a polydomain state due to a very small value of the shear angle, the direction of the [010] axes in domains with positive [001] axes almost coincides with the [100] axes in domains with negative [001] axes. So, a compression along [010] axes of positive domains makes them unstable. Simultaneously, this compression makes the negative domains more stable. Then, the domains with the [100] axes along the compression should grow and the domains with the [010] axes along the compression should disappear. This is the way to obtain a single domain sample.The temperature dependence of the spontaneous electric polarizationPS in GMO was measured in [13] over the temperature range from 4.2K up to the Curie point TC=432.3K. The electric polarization takes values of PS≈290×10-9Coul/cm2 at T=4.2K and PS≈40×10-9Coul/cm2 at T=TC=432.3K. The abrupt jump in the PS(T) curve was observed at T=TC=432.3K with a hysteresis loop of ΔT≈0.1K width. This is a direct indication on the phase transition of the first order. Room temperature values of PS in SMO, EMO, GMO, and TMO are in the range of PS≈140×10-9Coul/cm2 for EMO and PS≈240×10-9Coul/cm2 for SMO [14].The temperature dependences of electric and elastic properties of GMO were measured in [15]. The dielectric permittivity of the clamped crystal ε33κ does not depend on the temperature and shows no anomaly at the Curie point, while the elastic constant c66E shows strong temperature dependence with a marked anomaly at the Curie point. The dielectric permittivity at constant and zero stress shows a weak anomaly at the Curie point due to the piezoelectric coupling. This result indicates that the ferroelectric ordering in GMO is a secondary effect. It appears as a consequence of the structural β-β′ phase transition.The pyroelectric properties of GMO were studied in [16, 17]. Near TC, hard-to-control spurious domain nucleation was observed. It was explained by small thermal gradients caused by the thermal signal being detected.Bothβ and β′ phases of RMO are piezoelectric [17]. The electromechanical coupling coefficient is 4% at room temperature and 22% near TC.The presence of the domain walls significantly influences the resonant behavior of a crystal.


## 3. Magnetic Properties

Magnetic properties of a single-crystal and single-domain spherical sample of GMO were investigated in [18–22]. In [18–20], the magnetization along the [001],[010] and [100] axes, respectively, was measured at temperatures from 0.3 K to 4.2 K in magnetic field up to 90kOe. The saturation magnetization values at low temperatures and in high magnetic field along all crystal directions correspond to the magnetic moment μS=7μB which is typical for Gd3+ ion in the ground state 8S7/2. At temperatures above ≈1 K, the field dependences of the magnetization show a typically paramagnetic behavior. At T=0.5K, the field dependence of the magnetization exhibits a weak anomaly at H = (2-3) kOe that indicates on the field induced transformation in the antiferromagnetic structure from the antiparallel to parallel orientation of ionic magnetic moments [18]. At T=0.35K, this anomaly becomes more pronounced. On the basis of these data, the authors conclude that GMO becomes antiferromagnetic at temperatures lower than ≈0.3 K.Figure1 shows the field dependences of the magnetization of GMO at T=1K in the magnetic field along the [001] axis (open circles), the [100] axis (triangles), and the [010] axis (crosses). These dependences were calculated from the tables published in [18–20]. The solid line is the Brillouin function B7/2(μSH/kT) that corresponds to a free Gd3+ ion in the ground state 8S7/2 as it should be in a usual paramagnetic substance. The Brillouin function describes qualitatively the experimental magnetization curve along the [001] axis. One can see that the magnetization of the paramagnetic GMO at T=1K is saturated at the magnetic field H≈60kOe. The saturation magnetization is MS=39kG·cm3/moleGd3+=6.98μB/Gd3+ion. This value agrees well with the magnetic moment μS=7μB of Gd3+ ion in the ground state 8S7/2. The dependences of the magnetization in the magnetic field along the [100] and [010] axes are placed noticeably lower than for [001]. It is due to the natural magnetic crystalline anisotropy. The magnetization curves along the [100] and [010] axes are indistinguishable. It means that the [001] axis is the easy magnetization direction in GMO.Figure 1
The field dependences of the magnetization of GMO atT = 1 K in the magnetic field directed along the [001] axis (open circles), the [100] axis (triangles) and the [010] axis (crosses).Figure2 shows the experimental field dependence of the anisotropy energy for the magnetic field along the [100] axis. It was calculated from the data of Figure 1 using the relation:(1)Ea(H)=∫0H[M[001](H)-M[100](H)]dH.
In the state of the saturation (H=HS≈60kOe), the anisotropy energy is Ea(HS)=ES≈105erg/g≈4.5×10-5eV/Gd3+ ion. This value is higher than the anisotropy constant of ferromagnetic iron at room temperature (K1≈0.6×105erg/g). This is the energy needed to align all magnetic moments of Gd3+ ions in the paramagnetic GMO along the hard direction [100] at T=1K.Figure 2
The field depends on the natural magnetic crystalline anisotropy energy upon the magnetic field along the [100] axis.Figure3 shows the dependence of the anisotropy energy in GMO Ea(mr) upon the reduced magnetization mr=M[100]/MS along [100]. Ea(mr) is derived from(2)Ea(M[100])=∫0M[100][H[100](M)-H[001](M)]dM;
by introducing a new variable mr=M[100]/MS.Figure 3
The energy of the natural magnetic crystalline anisotropy of GMO versus the reduced magnetization. The magnetic field is directed along the [100] axis.It is seen that this dependence is described well by the relation:(3)Ea(mr)=ESmr4.
Here, ES is the anisotropy energy of a paramagnetic substance in a saturated state.In [21], the magnetization curves of GMO were measured in the magnetic field along the [001] axis in the temperature range (0.1–0.5) K using the method of adiabatic demagnetizing. The sharp jump was observed on the isentropic magnetization curve in the magnetic field H=1.615kOe. The temperature of the jump was 0.11K. The magnetic field H=1.615kOe is the critical field of the destruction of the antiferromagnetic order in GMO at T=0.11K. The magnetization was 350G·cm3/moleGd3+ at the beginning of the jump and ≈30 kG·cm3/mole Gd3+ at the end of the jump. The magnetization curve between these two points was a vertical straight line. So, the phase transition from the antiferromagnetic to the paramagnetic state in GMO goes without the stage of the spin-flop.The heat capacity of GMO was measured in the magnetic field up to90kOe at the temperatures (0.1–4.2) K in papers [18–22]. The thermodynamic functions were calculated using the measured values of the magnetization and the heat capacity.In [23], the magnetization and heat capacity of a single-crystal and single-domain spherical sample of TMO were measured over the range (0.4–4.2) K with fields up to 90kOe directed along the [001] axis. The magnetization was measured also along the [100] axis, but these measurements could not be performed in the field above 25kOe because of the mechanical destruction of the sample.According to [23], the temperature dependences of the heat capacity at zero field and at H=0.5kOe show λ-type anomalies with pronounced maxima at T=0.44K and 0.45K, respectively. As the magnetic field is raised, the λ maxima shift toward the lower temperatures. The temperature dependences of the magnetization and magnetic susceptibility show the broad maxima around T=1.6K. All of the above facts indicate on the antiferromagnetic ordering of the Tb3+ magnetic moments along the easy [001] axis with the Neel temperature TN≈(0.44–1.6) K.Above these temperatures, the magnetization curves along the [001] axis look like for a paramagnetic substance at very low temperatures and in extremely high magnetic fields; but the experimental value of the magnetic saturationMS=7.7μB/Tb3+ differs noticeably from MS=9μB/Tb3+ for a usual paramagnetic substance with Tb3+ion in the ground state 7F6. On the basis of this difference, the authors in [23] conclude that the total mechanical moment of the ground state doublet is Jz=±5. In the following, we will show on the basis of our magnetic measurements and the observations of Zeeman effect that it is not correct. The unusual value of the magnetic moment in TMO can be explained strictly by the quantum mechanical theory of magnetism of singlets with nondiagonal matrix elements in the angular momentum operator ofTb3+ ion.The magnetization curves along the [100] axis were measured only up to 25 kOe. The crystal was broken when the magnetic field was increased above 40 kOe. The authors explain the destruction of the sample by the magnetic torque. It is not a complete explanation of the self-destruction of TMO single crystal. Our later measurements showed that a sample of TMO at certain conditions can survive atT = 4.2 K even in the field of 200 kOe along the [100] axis.The absence of the magnetization data for TMO along the [100] axis at high magnetic fields makes it impossible to plot the anisotropy energy dependencies upon the magnetic field and the magnetization analogous to those for GMO on Figures2 and 3. The estimate of the anisotropy energy of TMO at H=25kOe and T=1 K gives Ea≈106 erg/g. From Figure 2, the corresponding value for GMO is ≈105 erg/g, that is, ten times lower than for TMO.It was established from the measurements of the heat capacity in [23] that the relaxation time of the nuclear spin system of 159Tb3+ becomes very long at high fields and low temperatures. For example, in the magnetic field H=90kOe at the temperature T≈1K, the order of magnitude of the relaxation time is several hundred seconds. The heat capacity values taken below those temperatures are not in equilibrium. The higher the magnetic field the higher is, the temperature below which one cannot reach the equilibrium between the nuclear system and the lattice formed by electronic shells of ions. This phenomenon complicates the magnetic measurements in TMO at low temperatures and in high magnetic fields.In [5], the magnetization of single-crystal single-domain samples of TMO was measured in pulsed magnetic fields up to 300kOe with pulse duration of 0.01 s. The magnetization curves were measured at T=4.2K and 78K along [001] axis and in the (001) plane. The magnetization process was adiabatic. The magnetization was highly anisotropic. The magnetization curves along the [001] axis were close to the paramagnetic saturation at both temperatures. At T=78K, H=300kOe, the magnetization reached 7.2μB/Tb3+ and at T=4.2K, H=240kOe, it was 7.8μB/Tb3+. Along the [010] axis, the magnetization was 3.4μB/Tb3+ at T=78K, H=300kOe, and 4.2μB/Tb3+ at T=4.2K, H=240kOe. The energy of the natural magnetic crystalline anisotropy of paramagnetic TMO was estimated qualitatively using the magnetization curves along the [001] axis and in the (001) plane. It proved to be of the order of Ea≈2·107erg/g at both temperatures. The magnetic anisotropy of the kind is unusual for a paramagnetic substance. It is the order of magnitude larger than the magnetic anisotropy constant of ferromagnetic cobalt K1≈106erg/g at 78K.Figure4 shows the magnetization curves of TMO in the plane (001) at T=4.2K. At the beginning of the pulse of the magnetic field, the [100] axis of the sample was oriented along the field. At H=190kOe, there was an anomaly on the magnetization curve. In the decreasing field, the magnetization curve was placed above that in the increasing field. Inspection of the sample in the polarized light after these measurements showed that the [100] and [010] axes were interchanged as compared to the state of the sample before switching on the pulse of the field. So, one could get the difference between the magnetization values for the field along the [100] and [010] axes. At H=180kOe, it was ≈0.6μB/Tb 3+. This interchange was observed also in [24, 25] from the measurements of the magnetostriction and the magnetoelectric effect. The interchange took place when the applied magnetic field was directed along the [100] axis and reached some critical value. In [26], the magnetization of a single-crystal single-domain sample of TMO was measured in the direct current magnetic field up to 190kOe in the temperature range T = (4.2–300) K. The magnetization curves at T=4.2K along the directions [001],[010] and [100] are shown on Figure 5. The magnetization curves along [010] and [100] on Figure 5 differ noticeably from the corresponding curves on Figure 4. The magnetization on Figure 4 is smaller than that on Figure 5. It is because the curves on Figure 5 are isotherms and those on Figure 4 are isentropes. The adiabatic process of magnetizing the sample on Figure 4 is accompanied by an increase of the temperature due to the magnetocaloric effect. According to our estimations, the temperature of the sample at H=240kOe on Figure 4 is ≈(10–20) K higher than at zero field. Respectively, the magnetization of Figure 4 is smaller than that of Figure 5.Figure 4
The magnetization isentropes in TMO at the initial temperature of the sample,T=4.2K. On the lower branch, the magnetic field is directed along the [100] axis. The anomaly at H≈190kOe shown by the arrow is due to the interchange of the [100] and [010] axes. On the upper branch, the field is parallel to the [010] axis.Figure 5
The magnetization curves of TMO atT=4.2K along main crystallographic directions. Solid lines are calculated using the theory of magnetism of singlets.There is no evidence of the interchange of the [100]↔[010] axes of Figure 5. On the contrary, Figure 4 shows the obvious transition of the sample from the state with the [100] axis parallel to the field into the state with the [010] axis parallel to the field at H=190kOe. This can be explained by the temperature dependence of the critical field of the [100]↔[010] interchange Hcr(T) and the magnetocaloric effect in pulsed fields. Our measurements of the magnetoelectric effect showed that the critical field Hcr(T) decreases from approximately 130kOe at T=10K down to 30kOe at T=30K. The [100]↔[010] interchange of Figure 4 was observed because of increasing the temperature due to the magnetocaloric effect and the corresponding decreasing of Hcr(T).It is seen that, along the easy direction [001], the saturation magnetization per one terbium ionMS≈8μB/Tb3+ is not equal to the magnetic moment of a free Tb3+ ion MTb3+=9μB/Tb3+. This is the result of the interaction of the angular momentum ofTb3+ ion with the crystal field. So, to describe the experiment, one should use the singlet magnetism theory.Before describing this theory, it is very useful to remind the main starting points of the standard theory of paramagnetism of noninteracting ions.The Zeeman energy is given by the following equation:(4)V̂Z=-μ̂LH⃗-μ̂SH⃗=-μB(L̂+2Ŝ)H⃗=-μBgJĴH⃗.
The Hamiltonian of the problem is(5)Ĥ=ĤCR+V̂Z=HCR-gJμB(HxĴx+HyĴy+HzĴz);ĤCR is the Hamiltonian of the crystal field; Ĵx, Ĵy, and Ĵz are operators of projections of the total mechanical momentum on the x, y, and z axes. This consideration gives only one of three contributions to the magnetic moment of the ion—the paramagnetism of orientation. We do not consider here the paramagnetism of polarization according to Van Vleck and the diamagnetism of precession.Consider a freeTb3+ ion in the magnetic field along z-axis. This is the case of the full degeneracy of the orbital multiplet:(6)ĤCR=0;Hx=Hy=0.
From (5) and (6), the Hamiltonian of a free ion is obtained:(7)Ĥ=-gJμBHzĴz;(8)Ĵz=-i∂∂φ.
Here, φ is the angle in the x,y plane. The φ-dependent parts of normalized wave functions of Ĵz for degenerate levels of the multiplet of a free ion are [27](9)|J,m〉=Φm(φ)=exp⁡(imφ)2π.
The fact is that a wave function is given by only one exponent plays a decisive role.Due to this fact, only diagonal elements of the Hamiltonian (7) can have non zero values. So, the secular equation is(10)∏m=J-J(gμBmHz-Em)=0.
Solutions of (10) give the eigenvalues of the energy of a free ion:(11)Em=gμBmHz;(12)m=-J,(-J+1),…,0,…,(J-1),J.
The components of the degenerate multiplet (11) possess the field independent magnetic moments:(13)Mzm=-∂Em∂Hz=-gμBm;
Using (11) and (12) and the statistics of Maxwell-Boltzman, one obtains the Brillouin function to describe the field and temperature dependence of the magnetization in a system of free ions:(14)Mz(H,T)=gJμBBJ(gJμBHzkT)=gJμB[2J+12Jcth(2J+12J×gJμBHzkT)-12Jcth(12J×gJμBHzkT)].
The total mechanical momentum of Tb3+ ion in 7F6 ground state is(15)J=L+S=3+3=6;L and S are the angular and spin moments, respectively. The g-factor is [28](16)g=1+J(J+1)+S(S+1)-L(L+1)2J(J+1)=32.
The corresponding saturation magnetic moment is μS=gJμB=9μB.The singlet magnetism theory describes the magnetic properties of ions in the crystal field. It is based on the Jahn-Teller theorem [29]. According to this theorem, the surrounding of the ion cannot be stable if the ground state of the ion in the crystal field is degenerate (excluding the case of Kramers degeneracy for ions with an odd number of outer electrons [30]). The spontaneous rearrangement of the crystal lattice goes on because the development of the instability reduces the symmetry of the crystal field and removes the degeneracy of the ground state. It means that rare earth ions with the even number of 4f electrons unavoidably have the singlet ground state [31]. If two lowest singlets are separated by a narrow gap of the order of 2Δ0≈(1–5)cm-1 and are much more away from other singlets, then the former two singlets can be considered as a quasidoublet. The magnetic moment of the ion with the quasidoublet as a ground state should be calculated by the singlet magnetism theory [31].The main distinction of this theory from the standard theory of the paramagnetism of free ions (6)–(16) is that the wave functions of the orbital singlets cannot be described by only one exponent like (9). The crystal field splits the degenerate multiplet of a free ion that is given by (11)–(12) and mixes the wave functions (9) of these levels. The wave functions of the singlets are linear combinations of the functions given by (9) with the values of m given by (12). The forms of the combinations are determined by the symmetry of the crystal field and of the operators: Ĵx, Ĵy, and Ĵz.The point symmetry group of the surrounding ofTb3+ ion in TMO is C2 [2]. For this symmetry, the Hamiltonian matrix for the magnetic field along z axis is(17)|-Δ0gJeffμBHzgJeffμBHz+Δ0|.
Here, 2Δ0 is the splitting of energy levels of the quasi-doublet under consideration.Jeff is a nondiagonal matrix element of the total mechanical momentum operator. It is given by the relation:(18)Jeff=〈ψz|Ĵx|ψS〉=6a6Sa6z+4a4Sa4z+2a2Sa2z;ψS and ψz are the wave functions of the lower and upper levels of the quasidoublet, respectively. The difference between MS and MTb3+ is not too large. This result indicates that the main harmonics that form the quasidoublet are (|6,6〉-|6,-6〉)/2 and (|6,6〉+|6,-6〉)/2. From here, it follows
(19)ψS=12[a6S(|6,6〉+|6,-6〉)+a4S(|6,4〉+|6,-4〉)+a2S(|6,2〉+|6,-2〉)+a0S];ψz=12[a6z(|6,6〉-|6,-6〉)+a4z(|6,4〉-|6,-4〉)+a2z(|6,2〉-|6,-2〉)].
One can see from (17) that for an ion with the ground state in the form of a quasidoublet the diagonal matrix elements do not depend upon the magnetic field and the nondiagonal elements are not zero and depend upon the magnetic field. This is a main distinction of the magnetic properties of the spectrum of singlets compared to the degenerate spectrum (11). The secular equation of the Hamiltonian (17) is(20)|-Δ0-EgJeffμBHzgJeffμBHz+Δ0-E|=0.
From (20), one obtains the eigenvalues of the energy of the quasi-doublet E1,2 in the magnetic field Hz:(21)E1,2=∓Δ02+(gJeffμBHz)2.
Here, E1—the energy of the lower level of the quasidoublet; E2—the energy of the upper level of the quasidoublet. Magnetic moments of the components of the quasidoublet are given by the following relation:(22)Mz1,2=-∂E1,2∂Hz=±2g2Jeff2μB2Hz2Δ02+(gJeffμBHz)2.
There is a principal distinction between ions with the nondegenerate spectrum (21) and ions with the degenerate spectrum (11). The energy levels of ions with the degenerate spectrum (11) have the field independent magnetic moments (13). These moments are not zero in the absence of the magnetic field if m≠0.As it is seen from (22), the magnetic moments of the levels of the quasidoublet essentially depend upon the magnetic field. They are zero at zero magnetic field. In other words, the ion with the nondegenerate spectrum of singlets can be magnetically polarized in the magnetic field. This is the paramagnetism of the polarization. So we came to the recognition of a very important fact that the idea of some definite and field independent magnetic moment of the ion is applicable only to a narrow class of problems in which the potential of the ion has the spherical symmetry and the spectrum of the ion is degenerate. The sample of the kind is GMO where Gd3+ ion in the ground state 8S7/2  has zero angular momentum and hence the spherical symmetry of the ionic potential. The overwhelming majority of problems, and TMO among them, insist on taking account of the field dependence of the ionic magnetic moment. This is done in the singlet magnetism theory by the relation (22).Using the Maxwell-Boltzman statistics, one can obtain the field and temperature dependence of the magnetization per ion for a singlet magnetic substance in the field directed along the [001] axis:(23)Mz(H,T)=g2Jeff2μB2HzΔ02+(gJeffμBHz)2⋅thΔ02+(gJeffμBHz)2kT.
This is a two-level version of the singlet magnetism theory. It does not describe the anisotropy of the experimental results shown in Figure 5. It gives a zero value for the magnetization along the [100] and [010] axes. To describe the anisotropy, one should take account of not less than three lowest singlets of the spectrum.The nearest surrounding ofTb3+ ion in TMO is an octahedron formed by seven oxygen ions [2]. The point symmetry group of the octahedron can be described as C2. The axis of the second order is a polar axis of the octahedron. Due to this symmetry, the [100] and [010] components of the magnetic moment Mx and My are zero in the two singlet model. Nonzero values for Mx and My appear when taking account of the nearest excited singlet. This is a three-level model of singlet magnetism theory.The sites ofTb3+ ions in TMO are of two types. The directions of the polar axes of the octahedrons at nonequivalent sites are approximately opposite to each other and parallel or antiparallel (approximately) to the [001] axis. The octahedrons at the nonequivalent sites are rotated with respect to each other by the angle of approximately 90° around the [001] axis. It means that the axes of the easy magnetization in the (001) plane should be different for ions disposed at the nonequivalent sites: the angle between them should be of 90°.Taking account of the nonequivalent sites forTb3+ ions, one obtains in the three-level modification of the theory the following expression for the field and the temperature dependence of the magnetization:(24)Mp(Hp,T)=∑i=1,2μip2Hpth(Bip/kT)/Bip2(1+exp⁡(dip/kT)/2ch(Bip/kT))+χpHp.
Here, p indicates the direction in the crystal and takes values [001],[010] and [100].i indicates the type of the site. Bip=bip2+(μipHp)2. Hp is a projection of the field on the direction p. Mp is a projection of the magnetic moment per one Tb3+ ion on the direction of the magnetic field. μip are the saturation magnetic moments of Tb3+ ions at sites of two types for main directions in the crystal. χp—Van Vleck susceptibility along the direction p. k—Boltzmann constant. bip and dip are coefficients of the theory that can be found from energy gaps Δn and δm between the singlets. Δn is the energy gap between the n-th excited singlet and the ground singlet for the sites of the first type, δm is the same for the sites of the second type. Δn and δm can be taken from the photoluminescence spectrum of Tb3+ ion in TMO [5]: Δ1≈δ1=(2.7±0.35)cm-1; Δ2≈δ2=(5.4±0.35)cm-1. Using these values of Δn and δm, the coefficients bip and dip were found. Then, the magnetization curves Mp(Hp) were calculated at the temperature T=4.2K for the directions [001],[010] and [001] by fitting expression (24) to the corresponding experimental curves. In these calculations, the coefficients μip and χp were used as the fitting parameters. It was found that the values μ1[001]=μ2[001]=7.9μB/Tb3+ along the [001] axis are identical for the nonequivalent sites. Along the [100] axis, the saturation magnetic moments for nonequivalent sites are different: μ1[100]=6.4μB/Tb3+≫μ2[100]=1μB/Tb3+. For the [010] direction, the opposite nonequality occurs: μ1[010]=0.05μB/Tb3+≪μ2[010]=6.9μB/Tb3+. This result reflects the difference in the symmetry of the nonequivalent sites of Tb3+ ions. A good agreement with the experiment takes place at T=4.2K for main directions. One can see that the difference between the experimental saturation magnetic moment of TMO and the magnetic moment of free Tb3+ ion in the ground state 7F6 is explained completely by the singlet magnetism theory. So, the additional assumption about the ground state of Tb+ ion in TMO with J=±5 [23] is not necessary.The comparison of the magnetization curves of GMO and TMO indicates definitely the role of the angular momentumL of the trivalent rare earth ion in the formation of the magnetic properties in rare earth molybdates.TheL value determines the main contribution to the energy Ea of the natural magnetic crystalline anisotropy. In GMO, Gd3+ ion in the ground state 8S7/2 has L=0, and the energy of anisotropy in the magnetically saturated state is Ea≈105erg/g. The magnetization reaches the saturation along hard directions [100] and [010] at H≈60kOe (see Figures 1 and 2). In TMO, Tb3+ ion in the ground state 7F6 has L=3, and the magnetization along hard directions [100] and [010] is far from saturation even in the magnetic field H=190kOe (see Figure 5). The energies of anisotropy at H=190 kOe found using relation (1) are Ea[100]≈1.13·107 erg⁄g and Ea[010]≈0.96·107 erg⁄g that is two orders of magnitude larger than in GMO. There is a strong correlation between the angular momentum value and the anisotropy energy. The ion with the nonzero angular momentum has a nonspherical charge cloud of the outer 4f-electrons. Magnetizing the crystal containing nonspherical ions like Tb3+ rotates these ions and hence changes the coupling between the ions and the crystal field much more than in the crystal with spherically isotropic ions like Gd3+. So, TMO is magnetically more anisotropic than GMO. The anisotropy energy values for TMO found in [26] agree by the order of magnitude with the value Ea≈2·107 erg/g found from the measurements in pulsed fields in [5]. The difference is mainly due to the adiabatic character of magnetizing in [5] in contradiction to the isothermal one in [26].The samples investigated in [5, 26] were not destroyed in contradiction to the result of [23] despite that magnetic fields in [5, 27] were three times higher than in [23]. It means that the magnetic torque acting on the sample in [5, 26] was also larger than in [23]. So, the reason of the destroying the sample in [23] was not the magnetic torque in despite of conclusion made in [23].The rotation of the nonspherical rare earth ions while magnetizing RMO crystals distorts the crystal lattice of RMO. The distortions of a ferroelectric crystal lattice unavoidably change its electric polarization due to the piezoelectric effect. This is the origin of the magnetoelectric effect to be discussed in the following section.


## 4. Magnetoelectric Effect

The presence of rare earth ions with nonzero angular momentum and nonzero total momentum in the crystal lattice of a ferroelectric substance leads to the large magnetoelectric effect (MEE). Rare earth ions with nonzero angular momentum possess the spatially anisotropic charge cloud of outer 4f-electrons. The spatial anisotropy of the magnetic ions determines large values of the natural magnetic crystalline anisotropy and the magnetostriction. These properties in the ferroelectric crystal ensure the strong dependence of the electric polarization upon the magnetic field, that is, MEE.MEE was predicted in [32], but the consideration of the problem was restricted by the linear magnetoelectric effect. It was established that this effect is possible only in the antiferromagnetic crystals. In [33], the prediction of linear MEE was given for a specific substance Cr2O3. The first experimental observations of linear MEE in Cr2O3 were done in [34, 35].The most important point of MEE in rare earth molybdates is the absence of the restrictions related to the magnetic structure. The general consideration of the problem shows that nonlinear MEE in rare earth molybdates can exist in the paramagnetic state.The first observation of nonlinear MEE in the paramagnetic tetragonalNiSO4·6H2O was published in [36].In [37], nonlinear MEE was observed in paramagnetic TMO at T=77K. That was the first observation of MEE in the rare earth molybdate family. In [38], MEE was measured in TMO and GMO in pulsed magnetic fields up to 110 kOe at T = (77–290) K.Figure6 shows the field dependence of the magnetically induced electric polarization (MEP) in TMO at T=77K [38]. The field was applied along the orthorhombic axis [11̅0]. MEP was measured along the [110] axis. As the spontaneous polarization P̅S is directed along the [001] axis, this result can be considered as the rotation of the P̅S in the magnetic field.Figure 6
The field dependence of the magnetically induced electric polarization in TMO along the [110] axis atT=77K.The pulsed magnetic field is directed along the [11̅0] axis. The solid line is the least-square fitting by P=p0H2. p0=4.12·10-19Coul·cm-2·Oe-2.The dependence of MEP upon the magnetic field is quadratic with a good accuracy. MEP isP=4.12·10-9Coul/cm2 at H=10T. It is small compared to the spontaneous electric polarization PS=180·10-9Coul/cm2 [14].MEP in GMO was equal to zero in the limit of the experimental error.The standard symmetry consideration [39] shows that, at a given geometry of the experiment, the components of the magnetically induced electrical polarization in the orthorhombic paramagnetic crystal are proportional to the products of the magnetization components:(25)Px~Mx⋅Mz;Py~My⋅Mz.
For a usual paramagnetic crystal,(26)Mz=χzzHz.
Here x,y,z correspond to the [100], [010], [001] axes, respectively; χzz is a component of the diagonal tensor of the magnetic susceptibility. It follows from (25) and (26) that the effect on Figure 6 should be zero.The contradiction between the result of calculations [39] and the experimental observation of the nonzero effect in [37, 38] can be eliminated by taking account of the local symmetry of Tb3+ ions surrounding. As it was noted in Section 3, the nearest surrounding of Tb3+ ion in TMO is an octahedron formed by seven oxygen ions [2]. The octahedron has the C2 point symmetry group. The axis of the second order is a polar axis of the octahedron. The symmetry axes of the octahedron do not coincide with the symmetry axes of the crystal. The polar axis of the octahedron makes an angle with the [001] axis, and the magnetic field applied in the (001) plane has a nonzero projection on the polar axis of the octahedron. According to the theory of magnetism of singlets, the magnetic field applied along the polar axis of the octahedron magnetically polarizes Tb3+ ion with the magnetic moment along the polar axis. So, the magnetic field applied in the (001) plane of the crystal induces a magnetic moment on Tb3+ ion along the polar axis of the octahedron. This moment, in turn, has a nonzero projection on the [001] axis of the crystal and, consequently, induces the electric polarization components along [100] and [010] axes according to (25). Hence, the bulk magnetic susceptibility tensor of TMO single crystal should have nonzero nondiagonal components and (26) should be replaced by(27)Mz=(χxz+χyz+χzz)Hz.
Then, there would be no contradiction between the measurements [37, 38] and the calculations [39]. The relation (27) should be verified experimentally.The magnetic susceptibility tensor of an individualTb3+ ion with the main directions along the symmetry axes of the oxygen octahedron has only diagonal components.At the beginning of this section, it was pointed out on the role of the angular momentum of a rare earth ion. A rare earth ion with a nonzero angular momentum has a nonspherical charge cloud of magnetoactive 4f-electrons. This charge cloud is coupled strongly with an ionic magnetic moment due to a strong spin-orbit interaction. It means that the rotation of the ionic magnetic moment by means of the magnetic field should result in the rotation of the charge cloud. The rotation of a nonspherical charge cloud induces strains in a crystal lattice—the magnetostriction. This should change the electric polarization of a crystal due to the piezo-effect giving rise to MEE. The effect is large in those molybdates, in which rare earth ions have large magnetic momentsμS=gJμB and nonzero angular moments L.In line with this mechanism, it was suggested that the MEP was due to the change of the spontaneous electric polarization of TMO induced by the single ion magnetostriction of Tb3+ ions [38]. This suggestion is the base of the magnetostriction model of MEE in rare earth molybdates. According to this model, the dependence of MEP upon the inverse temperature squared should be linear at high temperatures and at not very strong field. This can be shown as follows. MEP is caused by the piezoe-ffect. So it depends linearly on the strains. The strains are due to the magnetostriction. The magnetostriction at high temperatures and at not very strong magnetic field is proportional to the magnetization squared. The magnetization at these conditions is inversely proportional to the temperature. This leads to the linear dependence of MEP upon the inverse temperature squared P(T-2). This prediction was confirmed experimentally in [38].Figure7 shows the experimental dependence P(T-2) in TMO. MEP was measured along the [110] axis in pulsed magnetic field of 104.5 kOe directed along the [11̅0] axis. Solid line presents the least-square fitting. One can see that the prediction is fulfilled well.Figure 7
The dependence of the electrical polarization in TMO along the [110] axis in pulsed magnetic field of 104.5 kOe directed along [11̅0] axis.The movement of the ferroelectric domain walls initiated inTbGd(MoO4)3 (TGMO) by the pulsed magnetic field of 100kOe at T=78K was observed in [40, 41]. It was the first observation of the reversal of the spontaneous electric polarization due to the magnetic field.Figure8 shows the field dependence of the electric polarization in a single-domain sample of TMO at T=20K [42]. The polarization was measured along the [001] axis that is along the spontaneous electric polarization. In the initial state, the magnetic field was directed along the axis [100]. The negative sign was ascribed to the initial value of the spontaneous electric polarization on Figure 8. It is seen from Figure 8 that the magnetic field directed along the [100] axis decreases the absolute value of the spontaneous polarization. This decreasing is due to the negative paraprocess in the ferroelectric subsystem that is the decreasing of the parameter of the ferroelectric order. The negative paraprocess continues up to a critical field H=HC at which the magnetoelastic strains become large enough to nucleate a new domain. The new domain differs from the old one by the opposite direction of P̅S and by the mutual interchange of [100] and [010] axes. In orthorhombic rare earth molybdates, the ferroelectric domains are identical to the structural twin domains.Figure 8
The field dependence of the electric polarization in TMO measured along the [001] axis.T=20K.At H<50kOe, the field was directed along the [100] axis. At H>50kOe, the field was directed along the [010] axis.The samples investigated in [42] were put into the single domain state by the mechanical compression along the [010] axis. In the samples of this kind, the nucleation of the reversed domain usually goes on by arising of a flat domain wall in the (110) or (11̅0) plane. This is the beginning of the reversal of electric polarization. The reversal results in a jump of the electric polarization near a critical value of H=HC≈50kOe.The process of the displacement of the domain wall starts when the new domain reaches the critical size and begins to grow. In Figure8, the start of the displacement takes place at H≈45kOe. The new domain grows at the expense of the old one. During the growth, the domain wall move in the direction normal to its plane that is along the [110] or [11̅0] axes. At first, when the new domain is smaller than the old one, the absolute value of polarization of the sample being the sum of opposite polarizations of the two domains is decreasing. At H≈53kOe, the opposite domains are equal in size (P=0). Further, the absolute value of the polarization began to grow. If the critical field is applied for a long time, the new domain occupies all the volume of the sample. At H≈57kOe, the reversal of the electric polarization is completed and the displacement process ceases.After this, the increase of the polarization goes on by the mechanism of the positive paraprocess.The examination of the sample in polarized light after the measurements of Figure8 showed that after the jump the axes [100] and [010] were mutually interchanged. It is known that this interchange is accompanied by the reversal of the electric polarization. This supports the above interpretation of Figure 8.These experiments were the first to realize the reversal of the electric polarization in a ferroelectric crystal without using the electric field. The temperature dependenceHC(T) was measured in [42]. While decreasing the temperature from 100K down to 30K, the critical field decreases from 60kOe to 30kOe. Further decrease of the temperature results in an increase of the critical field up to 130kOe at T=15K. At lower temperatures down to 4.2K, the reversal of the electric polarization was not observed up to 200kOe.During the growth of the new domain in TMO, small rotations of the tetrahedronsMoO4 take place [2]. These rotations create defects in the crystal lattice. The more defects are in the lattice the larger is the critical field. That is why every reversal of the spontaneous electric polarization increases the critical field. When the number of the reversal cycles was large enough the sample was mechanically destroyed into small particles at the next in turn reversal.Figure9 shows the field dependence of MEP with a sharp break down of MEP at H=71.5kOe. It is due to the mechanical destruction of the sample. Before the destruction, the spontaneous polarization of the sample was reversed 5 times at temperatures from 98 K to 20 K. One can compare this number of complete reversals with the only attempt to magnetize the sample of TMO along the [100] axis in Fishers experiment [23]. That attempt failed at H=40kOe because of the destruction of the sample. To the authors’ opinion, the reasons of the destruction were the spherical shape and the pin holes of Fishers sample. The pin holes were the origins of a large number of the domain walls of different orientations. The domain walls intersected each other and very large stresses were concentrated at those intersections. These stresses destroyed the sample. The spherical shape stopped the movement of the domain walls because when a flat domain wall moves through the sphere, the area of the wall has to change its size and this also generates the stresses with the evident destruction. The samples used in [42] were rectangular parallelepipeds. They had dimensions 7×7×1 mm3. The large face of the sample was the (001) plane and the small faces were the (110) and (11̅0) faces. At this orientation, the domain walls were always parallel to the small faces and normal to the large face of the sample. This orientation ensured the easiest conditions for the domain walls to move. So, the samples sustained several reversals of the spontaneous electric polarization and several tens of magnetizing up to H=200kOe along the [100] axis without the reversal.Figure 9
The field dependence of MEP in TMO atT=5K. The field is directed along the [100] axis. The positive direction of the spontaneous polarization coincides with that one of the ΔP axis. A sharp breakdown of MEP at H=71.5kOe is due to the mechanical destruction of the sample.The anisotropy of MEE in TMO was investigated experimentally in [43] in the magnetic field up to 60kOe at temperatures 4.2K; 1.8K and at T=0.39K up to H=100kOe. MEP was measured along the [001] axis. The magnetic field was applied in the (001) plane. The angle φ between the [010] axis and the magnetic field was varied from 0 to 90°. The experimental angular dependences of MEP did not obey the relations of the standard theory of the anisotropy of MEE in a paramagnetic crystal derived in [36]. According to [36], the dependence of MEP in a tetragonal paramagnetic crystal upon the strength of the magnetic field and upon its direction can be described by the relation:(28)ΔP=ξH2⋅sin⁡2φ.
Here, φ is the angle between the magnetic field and the tetragonal axis; ξ is the constant of the theory.Figure10 displays some of the experimental isotherms of MEP measured along the [001] axis in TMO at 4.2 K (curves 1–5) and at 0.39 K (curves 6–10) for various directions of the field in the (001) plane. The sign of ΔP(H) in Figure 10 was determined by comparing it with the sign of the jump of the electric polarization when the spontaneous polarization switches in the critical magnetic field (see Figure 8). The field dependences of MEP are nonlinear. MEP is positive in the magnetic field directed along the [010] axis (a positive paraelectric process in the ferroelectric subsystem at φ=0, curves 1 and 6). MEP is negative in the magnetic field directed along the [100] axis (a negative paraprocess at φ=90°, curves 5 and 10).Figure 10
MEP in TMO induced along the[001] axis by the magnetic field applied in the (001) plane for various values of the angle φ between the field and the [010] axis at T=4.2K (curves 1–5) and T=0.39K (curves 6–10); (1) φ=0°; (2) 35.5°; (3) 57.7°; (4) 66.6°; (5) 90°; (6) 0°; (7) 35.5; (8) 57.7°; (9) 68.8°; (10) 90°.The field dependencesΔP(H) do not agree with (28) and what is more unusual is that the ΔP(H) itself depends upon the direction of the field.The angular dependences of MEP in TMO do not agree with (28) neither.Figure11 displays some of the experimental angular dependences of MEP in TMO at temperatures T=4.2K (curves 1−3) and T=0.39K (curve 4). Curves 1, 2, 3 were measured at H=20kOe, 40kOe, and 60kOe, respectively. Curve 4 was measured at H=100kOe. The line 4 is drawn for the eyes guide. The lines 1–3 are calculated following the theory of mean quadrupole moments [43, 44]. The points refer to the experimental values of ΔP(H,φ). Their arrangement is not symmetric relative to φ=45° contrary to (28). The absolute values of ΔP(φ=0) and ΔP(φ=90°) are different. The values of ΔP(φ=45°)≠0. So, (28) of the theory [36] does not describe MEE in TMO. Obviously, a new approach is needed to describe MEE in TMO. This approach was developed in [44] and applied to the problem in [43] using the experimental field dependences of the magnetization [26].Figure 11
The angular dependences of MEP in TMO. MEP was measured along the[001] axis. The field was applied in the (001) plane. φ is the angle between the [010] axis and the field. Curves 1, 2, and 3 were measured at T=4.2K in various fields: curve 1 at H=20kOe; 2—40kOe; 3—60kOe. Curve 4 was measured at T=0.39K and H=100kOe.The magnetic-field-induced electric polarization in crystals (MEP) is, like the magnetostriction and magnetic linear birefringence, an even-parity magnetic effect. The classical theory of even effects was developed in [45, 46]. This theory describes well the properties of weakly anisotropic compounds, in particular, of those containing d-ions. The situation with f-ions is more complicated. According to [44], even magnetic effects in rare earth compounds are governed by magnetic susceptibilities of the mean multipole moments of the rare earth ions 4f-shell. Their calculation requires knowledge of the energy levels and wave functions of these ions. It was established in [26] that the ground state of Tb3+ ions in TMO is apparently a combination of three close-lying singlets, which belongs to different irreducible representations of the symmetry group of the ion environment and are separated by 2.7cm-1 in energy.In this case, at the magnetic field in the (001) plane, the following relation takes place [43]:(29)ΔPz(T,H,φ)=AHxMx(T)+BHyMy(T)+CHxMy(T)+DHyMx(T)+EHx2+FHxHy+GHy2.
Here, T is the temperature. It should not be very high, but the condition T≫TN should be fulfilled. TN is the temperature of the antiferromagnetic ordering. In TMO, TN=0.45K [23].A,B,C,D,E,F, and G are constants of the theory which do not depend on the temperature; Hx=H[010] is the magnetic field component along the [010] axis:      (30)Mx=M[010]-χVV[010]H[010];M[010] is the magnetization component along the [010] axis.Hy=H[100] is the magnetic field component along the [100] axis:  (31)My=M[100]-χVV[100]H[100];M[100] is the magnetization component along the [100] axis.χVV[010] and χVV[100] are the Van Vleck susceptibilities along the [010] and [100] axes, respectively. Mx and My are the magnetization components of the Tb3+ ion quasitriplet in the (001) plane. The values of the constants A,B,C,D,E,F, and G were determined in [43] using corresponding experimental dependences ΔP(H) from [43] and Mx,y(Hx,y) from [26]. The solid lines 1, 2, and 3 on Figure 11 were calculated using (29), the values of A,B,C,D,E,F, G, and the dependences (30) and (31) that were calculated in turn for various values of φ from the experimental data of [26]. It is seen that the theory agrees well with the experiment. A satisfactory agreement was obtained in analogous calculations for T=1.8K. At T=0.39K, there is no agreement between the theory and the experiment because the theory does not take account of the antiferromagnetic order that arises at TN=0.45K.Four mechanisms of changing the electric polarization were established in TMO: the positive and negative paraelectric processes, the rotation of the electric polarization, and the displacement of the domain wall.As it was noted above, MEE inGd2(MoO4)3 (GMO) at T=77K was not observed. This fact is in a qualitative agreement with the magnetostrictive model [38]. However, MEE in GMO is not forbidden utterly by this model. In spite of a zero angular momentum of Gd3+ ion, the natural magnetic crystalline anisotropy in GMO is not zero. The calculations based on data of Figure 1 derived from results of [18–20] give for the anisotropy energy of GMO in the magnetically saturated state at T=1K the value of ≈105 erg/g. It is higher than the anisotropy constant in iron at room temperature. It is well known that the existence of the magnetic crystalline anisotropy involves the existence of the magnetostriction. From these reasoning, it follows that MEE can exist in GMO at low temperatures. MEP in GMO was measured in [47] along the [001] axis at T = (5–300) K for two directions of the magnetic field along the [100] and [010] axes and for H≤190kOe.AtT=5K and H=190kOe, MEP was ΔP≈0.3·10-9Coul/cm2. This fact confirmed the magnetostrictive model of MEE. In [48], the anisotropy of MEE in GMO was investigated in the magnetic field up to 200 kOe at temperatures T=4.2K and 0.4K. The experimental angular dependences of MEP at T=0.4K and various field values are shown on Figure 12. The solid lines are calculated by (28) shifted along the φ axis by φ0=π/2. One can see that the experimental dependences ΔP(φ) agree with the relation (28) through the whole field range investigated. The experimental dependences ΔP(φ) at T=4.2K agree with (28) as well. It is evident that at higher temperatures, the theory [36] developed for a paramagnetic substance should work also. So, (28) of the standard theory [36] describes the angular dependences of MEE in GMO satisfactorily in the ranges of variables H=(0–200)kOe and T>0.4K satisfactorily.Figure 12
Angular dependences of MEP in GMO at various values of the magnetic fields.φ° is the angle between the [010] axis and the magnetic field. T=0.4K. Curve (1) H=10kOe; (2) 100kOe; (3) 200kOe.The experimental field dependencesΔP(H2) in GMO at T=4.2K and 0<H<20kOe are shown in Figure 13 for two orientations of the field along the [010] axis (filled circles) and along the [100] axis (open circles).Figure 13
The dependences of MEP in GMO upon the magnetic field squared atT=4.2K. Filled circles were obtained in the magnetic field along the [010] axis; open ones—the field along the [100] axis.It is seen that, for both directions of the field, the experimental isotherms agree with the relation (28) only at H2<40(kOe)2, that is, H<6.3kOe. At higher fields, the disagreement increases. At T=0.4K, the theory [36] fails utterly. So, the theory [36] is practically inapplicable to the description of MEE in GMO at low temperatures. The main origin of the magnetoelectric effect in rare earth molybdates is the single ion magnetostriction of rare earth ions [37, 38]. So, we supposed that the experimental dependence P(m) in GMO can be described by the relation analogous to that obtained in [49] for the single ion magnetostriction:(32)ΔP(m(H,T))=ΔP0⋅Îl+1/2[L-1(m(H,T))].
Here, Îl+1/2(x) is a normalized hyperbolic Bessel function of l-th order:(33)Îl+1/2(x)=[Il+1/2(x)]/[I1/2(x)];Il+1/2(x)=∫-11Pl(y)exp⁡(xy)dy;Pl(y) is the Legendre polynomial. L-1(m) is the inverse Langevin function of the reduced magnetization m=M/MS. For the case considered, l=2 has been defined by the second-order symmetry of the angular dependence P(φ) obtained from Figure 12. For l=2, the reduced hyperbolic Bessel function is given by the relation:(34)Î5/2(x)=1+3x2-3xcthx.
Relation (32) was used to describe the field dependence of MEE in GMO. Results are shown on Figure 14. Points on this figure present experimental dependences of MEP in GMO upon the reduced magnetization m. The circles present MEP values measured in the field directed along the [010] axis. The squares correspond to the field along the [100] axis. Filled symbols refer to the temperature T=0.4K, open symbols refer to T=4.2K. The experimental values of MEP do not exceed ΔP≈0.3·10-9Coul/cm2.Figure 14
Dependences of MEP along the [001] axis in GMO upon the reduced magnetizationm=M/MS at various temperatures and magnetic field directions. T=4.2K: curve (1) the magnetic field along the [100] axis; curve (2) the field along the [010] axis. T=0.4K: curve (3) the field along the [100] axis; curve (4) the field along the [010] axis.Calculated dependencesΔP(m) are shown by solid lines in Figure 14 for T = 4.2 K (curves 1 and 2). Coefficients ΔP0 are -0.18·10-9C/cm2 and 0.17·10-9C/cm2 for the [100] and [010] axes, respectively. For 0<m<0.8, the divergence between the experiment and the theory is δ(ΔP)/(ΔP)<0.02. So, the relation (32) obtained from the theory of the single ion magnetostriction agrees to the experiment well at 0<m<0.8. For m = 0.95,   δ(ΔP)/(ΔP)=0.09. The value m=0.95 in the case of GMO corresponds to the field H≈56,7kOe. It is seen that the range of validity of relation (32) is wider than relation (28) [36]. Nevertheless, the most part of the field range is not described.The anomaly nearm=1 indicates some additional effects that should be considered. The dependences P(m) at T = 0.4 K cannot be explained by the theory of the single ion magnetostriction in its present form because this theory does not take account of the antiferromagnetic structure. The magnetic phase transition paramagnetism-antiferromagnetism that takes place in GMO at T = 0.3 K may affect the results of the measurements at T = 0.4 K. The solid lines 3 and 4 are drawn only for the eyes guide.MEE inSm2(MoO4)3 (SMO) was investigated in direct current magnetic field up to 200kOe at temperatures from 4.4K to 0.43K  [50]. MEP was measured along the [001] axis in the magnetic field along the [010] axis and [100] axis.Figure15 represents a cycle of the MEP variation (abcde) observed in a single domain sample of SMO at T=1.4K  in the magnetic field directed along the 010] axis. At the ascending branch of the cycle (ab), the magnetic field was increased at a rate dH/dt=0.755kOe/s. The branch ab was passed for 265 s. It represents the experimental dependence of MEP in the increasing field ΔPe↑(H). At the b point (H=200kOe), the MEP was ΔPe↑=-5.46·10-9Coul/cm2. At this point, the magnetic field was fixed and remained unvaried for 232 s. During this time, the MEP relaxed to ΔPe↓=-3.01·10-9Coul/cm2 (point c). At point c, the magnetic field started to decrease at a rate dH/dt=-0.755kOe/s. The descending branch cd represents the experimental dependence of MEP in the decreasing field ΔPe↓(H). At point d, the magnetic field reached zero and was fixed. The MEP was found to be ΔPe↓=3.12·10-9Coul/cm2 at this point. After fixing the magnetic field at zero, the MEP was measured for 338s. During this time, it relaxed to ΔP=1.4·10-9Coul/cm2.Figure 15
Field dependences of MEP inSm2(MoO4)3 at T=1.4K. The field is directed along the [010] axis. Curves ab and cd present experimental dependences ΔPe↑↓(H) in increasing and decreasing field, respectively. Straight lines bc and de present the relaxation of ΔP with the time at H=200kOe and H=0, respectively. Curves af and fa are calculated from the curves ab and cd using (36) and (37), respectively.This variation of the MEP is plotted as de on Figure15. Thus, the MEP in SMO relaxes in a fixed magnetic field. The relaxation is not due to the charge leakage through the electrometer input circuit because the times of measurements at a fixed field 232 s and 338 s are two orders of magnitude shorter than the time constant of the electrometer τem>104s. Analogous measurements of the MEP cycles were performed for other magnetic field orientations and temperatures. The time dependences of MEP on the segments bc and de were fitted well by the exponents:(35)ΔPe(H=const,t)=ΔP0exp⁡(-tτ)+PR.
The relaxation time τ increased with decreasing temperature. It depended upon the magnetic field only weakly.Table2 presents the values of τ obtained at various temperatures in the field 200kOe.Table 2
Relaxation timeτ of MEP in SMO at various temperatures in the field 200 kOe.
T (K)The field along [010],τ (s)The field along [100],τ (s)4.41271282.62791921.43374030.94904900.43534828The relaxation of the MEP noticeably affected the MEP field dependences because the measurement time 265 s was comparable with the relaxation timeτ. Our measurements showed that the relation between ΔP0 and PR can be approximately described by the expression PR≈β(ΔP0) where β≪1 is a constant. In this case, a time-dependent contribution could be eliminated from the experimental field dependences ΔPe↑↓(H) using the relations:(36)ΔP↑(H)=ΔPe↑(H)+∫0HΔPe↑(H)ατ(1+β)dH;(37)ΔP↑(H)=ΔPe↑(H)-∫HmHΔPe↓(H)ατ(1+β)dH+P↑(Hm)-Pe↓(Hm).
Here, P↑↓(H) does not depend on the time. The ↑ and ↓ arrows denote the values obtained under increasing and decreasing field respectively; Hm is the maximum field reached in the measurements; α=dH/dt is the rate of the magnetic field variation with the time; P↑(Hm) is the true MEP at the maximum field calculated from the experimental curve Pe↑(H) using (36); Pe↑(Hm) is the experimental value of MEP at the start of the decrease in the field (point c in the cycle in Figure 15). Curves ab and cd were used to obtain the af and fa curves (Figure 15) by means of (36) and (37), respectively. One can see that elimination of the time-dependent contribution from the experimental ΔPe↑↓(H) dependences yields a hysteresis-free P(H) relation. Analogous results were obtained at the temperatures 0.9; 1.4; 2.6; 4.4 K in the field along the [010] and [100] axes.Corresponding curves are shown in Figure16. The dependences 1 and 6 obtained at 4.4 K are fitted well by the relation (28) at φ=0;π/2 and H<50kOe. At higher fields and lower temperatures, this relation fails.Figure 16
Field dependences of MEP in SMO after the elimination of the time-dependent contribution.  Curves 1–5 are measured in the field along the [010] axis; curves 6–10 are measured in the field along the [100] axis. The temperatures: curves (1, 6)T = 4.4 K; (2, 7) 2.6 K; (3, 8) 1.4 K; (4, 9) 0.9 K; (5, 10) 0.43K.The data described in Section4 qualitatively confirm the magnetostriction mechanism of MEE in rare earth molybdates. This mechanism is based on the assumption that the magnetostriction affects the electric polarization of these ferroelectric compounds.Table3 shows the values of MEP in the three rare earth molybdates investigated here at fixed temperature and field values. The values of g and J were taken from [51]. There is a clear qualitative correlation between the values of the angular momentum L and the values of ΔP. GMO with zero angular momentum of R3+ ion possesses the least ΔP value. The angular momentum L characterizes the spatial anisotropy of an ion. It is obvious that the rotation of an anisotropic ion with L≠0 produces more strains in a crystal lattice than the rotation of a spherically symmetrical ion with L=0. So, the change of the electric polarization (MEP = ΔP) is larger for L≠0 than for L=0. There is some incongruity in this conclusion: the angular momentum of Sm3+ ion is larger than that of Tb3+ ion, but the ΔP value in SMO is the order of magnitude smaller than that in TMO. This fact can be explained qualitatively by relations between the ionic magnetic moments gJμB of R3+ ions under consideration. Here, g is g-factor; J=L±S is the total mechanical moment of R3+ ion; S is the spin mechanical moment of R3+ ion. The sign “−” corresponds to light rare earth ions from Ce to Eu, and the sign “+” corresponds to heavy rare earth ions from Tb to Yb. μB is the Bohr magneton. One can see that the magnetic moment of Tb3+ ion is larger than that one of Sm3+ ion. So, one should conclude that the total effect of the space anisotropy and the magnetic moment of Tb3+ ion is larger than that one of Sm3+ ion. It can be because of that the ion with the larger magnetic moments rotates by the larger angle than the ion with a smaller magnetic moment in the same magnetic field. The final result is determined by the quantitative contributions of the spatial anisotropy of ions and the angles of their rotations due to the interaction of their magnetic moments with the applied magnetic field.Table 3
A comparison of MEP absolute values in TMO, GMO, and SMO atT=4.2 K and H=100 kOe.
CompoundFigureT (K)ΔP (10−9 Coul/cm2)L of R3+ iongJ of R3+  ionTMOFigure110.39≈(42–57)39GMOFigure120.4≈0.2507SMOFigure160.43≈350.7


## 5. Magneto-Optical Properties

In [26], the experimental dependences of the magnetization in terbium molybdate β′-Tb2(MoO4)3 (TMO) upon the magnetic field and the temperature M(H,T) were described satisfactorily using the theory of the magnetism of singlets [31] at low temperatures T<30K. At higher temperatures, the agreement was not good. It was due to increasing population of highest levels of 7F6 multiplet of Tb3+ ions. For the calculation of the magnetization at high temperatures, the highest levels of 7F6 multiplet should be known. To calculate the spectrum of 7F6 multiplet, one should know the crystal field acting on Tb3+-ion. Now, there is no theory that enables one to calculate the crystal field parameters from the first principles. So, the investigations that allow obtain directly the data about the spectra of multiplets are very useful. For example, the knowledge of the spectrum of the excited multiplet 5D4 can simplify the calculations of the crystal field parameters in the case under investigation. The spectrum of 5D4 multiplet can be found from the experimental data for the optical absorption in TMO at low temperatures. The optical absorption in TMO was studied at zero magnetic field in [52]. 14 well-resolved absorption peaks were observed. It was established that those peaks correspond to the transitions 7F6⇒5D4 in Tb3+ ions.In [53], the absorption spectra were measured in TMO in the nonpolarized light in the magnetic field H≤200kOe along the [001] and [110] axes at T=4.2K and T=1.7K. A wave vector of the light was parallel to the magnetic field. At zero field, 14 well-resolved peaks with wave numbers ν=20519cm-1–20655cm-1  (λ=4873.5Å–4841.4Å) were observed. The magnetic field shifted the spectrum toward large ν values. The spectra were measured at 29 different field values. The experimental field dependences of the wave numbers of the absorption peaks νi(H) were obtained for both directions of the field.The spectra did not depend upon the temperature. There was a pronounced dependence of the peak numbers and of the rate of their shift with the fielddνi/dH upon the direction of the field. The analysis of these dependences showed that the excited multiplet 5D4 of Tb3+ ion in the crystal field of TMO could be considered as consisting of five singlets and two quasidoublets. On the base of the measurements and analysis of the absorption spectra in terbium molybdate, the field dependence of the energy spectrum Ei(H) of the excited multiplet 5D4 of Tb3+ ion in the crystal field with the point group symmetry C2 was found for the magnetic field up to 200 kOe along the [001] axis. The experimental dependences Ei(H) could be described well by the theory of the magnetism of singlets.The number of the absorption peaks was explained by the existing of two types of nonequivalent sites forTb3+ ions in the crystal lattice of TMO.


## 6. Anomalously High Photo Voltages

In [54], an anomalously high photoinduced voltage (APV effect) in TMO was observed. The illumination of a single crystal of TMO by the beam of the argon laser with a wave length of λ= 4880 Å and a power of 0.5W induced in the ferroelectric single domain sample an additional electric polarization. The wave length corresponded to a proper absorption band of TMO. The laser beam was directed along the [110] axis. The surface of polarization of the laser beam was normal to the (001) plane. The direction of the photoinduced polarization was the same as that of the spontaneous polarization. The laser beam was focused to a spot with an area 10-2cm2. It means that the distribution of the intensity of the light was highly inhomogeneous. The average value of the intensity in the illuminated spot was ≈50 W/cm2. The photoinduced electric polarization of the crystal depended upon the time of the illumination. After the illumination during 4.5 h, the photoinduced polarization was approximately three times larger than the spontaneous electric polarization and showed the tendency to the saturation.The laser beam with the wave length4880Å induced in TMO the luminescence in the green part of the visible spectrum with the wave length ≈5430Å. The quantum output constituted 20%.In [55], the dependences of the photoinduced voltage in TMO upon the time t of the illumination by a laser beam were measured at room temperature at 0.1s≤t≤4500s. The power of the laser radiation was 0.4W. The distribution of the intensity along the sample was varied from highly inhomogeneous to homogeneous. It was found that the photoinduced voltages of two types with opposite signs appear due to the inhomogeneous illumination.In [55], the APV effect in TMO was measured at heightened sensibility, respectively, to the photoinduced electric polarization of the crystal.Due to this, it was possible to observe relatively small contribution to the photoinduced polarizationΔPph that was directed opposite to the spontaneous electric polarization PS. This effect was observed at a small duration of illumination t<80s. It was called APV-1. The effect that was observed in [54] was called APV-2. The sign of APV-1 was opposite to the sign of APV-2. The voltage and the relaxation time of APV-1 were essentially smaller than those of APV-2. That is why the APV-1 effect was not observed on the background of the APV-2 effect in [54].The first observation of APV-effect was reported for a single crystal sample of ferroelectricSbSI0.35Br0.65 [56]. A phenomenological theory of the APV effect was described in [57]. In works published up to now, the diffusion contribution to the photoinduced current was ignored [57, 58]. According to [57], the time dependences of the photoinduced current and the photoinduced voltage are given by the following expressions:(38)i(t)=i0exp⁡(-tτ);(39)V(t)=V0[1-exp⁡(-tτ)].
The effect APV-1 is well described by the expressions (38) and (39). The effect APV-2 does not follow these expressions.Figure17 shows experimental time dependences of the photoinduced current i(t) (curve 1) and the photoinduced voltage V(t) (curve 2).Figure 17
The dependence upon the time of the photoinduced current (curve 1) and the  photoinduced voltage (curve 2) in TMO.  Sample number 1. The cross-section of the beam has size0.5mm×2mm.These curves were obtained on the sample number 1 with size6×7×1mm3. The large face of the sample was parallel to the plane (001).The [100] and [010] crystal axes were directed along the bisectors of the right angles of the (001) sides. The laser beam was focused on the7×1mm2 face of the sample by a cylinder lens to the spot with dimensions 2×0.5mm2. The large size of the spot was along the [001] axis of the crystal. The intensity of the light in the spot was ≈40 W/cm2.Arrows show the time when the light was switched on and off. These dependences are in a good agreement with expressions (38) and (39) respectively. The constants are i0≈0.2·10-9A; V0≈-10V;τ≈5s. The stepped form of the curve 1 is due to the resolution of the equipment. When the light is switched on, the negative voltage V(t) appears. This voltage is due to the photoinduced electric polarization ΔP⃗ph that is directed opposite to the spontaneous electric polarization P⃗S. The value V0≈-10V corresponds to the decreasing of P⃗S by the value ΔPph=4.8·10-9C/cm2. When the light is switched off, the photoinduced voltage decreases exponentially with the time constant τ≈5s. The photoinduced voltage and the photoinduced current on Figure 17 are due to the effect APV-1.Figure18 shows the dependences ΔPph(t) for various distributions of the light intensity along the sample. The switching on and off the light are shown by the arrows. Curve 1 was measured on the sample number1 at the same conditions as in Figure 17. One can see that the type of the dependence shown in Figure 17 occurs only during short time of ≈200 s after the light is switched on. At t≈(250–300)s, the ΔPph(t) value reaches a minimum of ΔPph≈-4.8·10-9C/cm2 and begins to increase with increasing t.Figure 18
The dependence of the photoinduced electric polarization in TMO upon the time. Curve 1. Sample number 1. A nonhomogeneous illumination. Curve 2. Sample number 1. A homogeneous illumination of the face of the sample. A nonhomogeneous distribution of the intensity in the sample due to the absorption. Curve 3. Sample number 2. An approximately homogeneous distribution of the intensity in the sample.Att≈1700s, the sign of ΔPph(t) is changed. The positive sign of ΔPph(t) at t>1700s corresponds to the coincidence of the directions of P⃗S and ΔP⃗ph. The tendency to the saturation of ΔPph(t) is seen at large values of t. At switching off the light, the positive jump of the photoinduced polarization occurs. It is of the same absolute value as the negative jump at switching on the light. After this jump, the photoinduced polarization continues to increase slowly: an after-effect takes place. In other words, the long negative pulse of the polarization due to the effect APV-1 exists on the background of slow positive drift of the polarization. This drift is due to the effect APV-2.Curve 2 in Figure18 was measured on the sample number 1 with a homogeneous distribution of the intensity across the 7×1mm2 face of the sample. The cross-section of the laser beam had a circular form with diameterϕ=7mm. So, the intensity of the illumination was lower than for curve 1, and the photoinduced polarization was smaller. One can see that curve 2 shows a smaller relative contribution of the slow positive drift to the resulting polarization than curve 1. It should be pointed out that in spite of the homogeneous illumination of the 7×1mm2 side the distribution of the intensity of the light was not homogeneous inside the sample along the direction of the propagation of the light due to the absorption. To obtain the ΔPph(t) dependence at still more homogeneous illumination, the measurements were performed on the sample number 2. The sample number 2 was a cube of the size of 2×2×2mm3. The relative directions of the crystal axes and of the laser beam were the same as for sample number 1. The cross-section of the beam was as for curve 2.As the illuminated face of the sample number 2 was smaller than that for sample number 1, theΔPph(t) value also was smaller. As the size of the sample number 2 along the direction of the propagation of the light was three times smaller than for the sample number 1, the relative nonhomogeneity of the illumination was also smaller. Curve 3 in Figure 18 was measured on the sample number 2. One can see that in this case the change of the sign of the voltage was not observed during all the time of the illumination. At these conditions, the positive drift due to APV-2 is so small that one can neglect it.Figure19 shows curve 1 from Figure 18 contracted along the t-axis and its prolongation for large t values that describes the relaxation of the photopolarization. The relaxation curve gives qualitative information that is important for explanation of the APV effect in TMO. The curve shows a very durable aftereffect. The increase of the photoinduced voltage continues during ≈3500 s after the light is switched off. The voltage reaches the maximum Vmax≈35V at t≈8000s and then decreases.Figure 19
The relaxation of the photoinduced voltage in terbium molybdate after  measuring curve 1 on Figure18.The value ofVmax≈35V corresponds to the photopolarization ΔPph2≈17·10-9C·cm-2. The relaxation time is τrelax≈20h=7.2×104s. It is more than the time constant of the input circuit of the electrometer τinput. Due to this, the dependence V(t) of the relaxation process is not reproduced correctly. So, only the initial part of the relaxation curve is shown in Figure 19. The complete relaxation of the APV-2 effect for the present experiment was achieved at tR≈150h=5.4×105s. Due to the relation between τinput and τrelax, this value of tR should be considered as a bottom estimate.The two effects described above show quite different dependences upon the conditions of the experiment.The APV-1 effect exists both at homogeneous and nonhomogeneous illuminations (Figure18). It reaches the saturation at t≈40s (Figure 17). The APV-2 effect exists only at nonhomogeneous illumination.It was not observed at a homogeneous illumination (Figure18, curve 3). It depends on the time and the geometry of the beam cross-section (Figure 18).The time constant of the APV-1 effect isτ=5s. From Figures 17 and 18, one can see that APV-1 effect reaches the saturation in a relatively short time and after this the voltage due to APV-1 effect is constant. Contrary to APV-1 effect, the voltage due to APV-2 effect grows for a very long time of the order of several thousand seconds showing only very weak tendency to the saturation.The APV-2 effect possesses the after-effect: thePph2 value increases during several thousand seconds after the light is switched off (Figure 19). No after-effect was observed while measuring the APV-1 effect.The APV-1 effect has a feature that is common to photo-induced currents measured in [59–62]. It is a strong correlation between the switching the light off and the change of the direction of the photoinduced current. The photoinduced voltage due to this effect begins to decrease immediately after the light is switched off. At the same time, the direction of the corresponding current is changed. The APV-2 effect is essentially different from those observed in [59–62]. The photoinduced voltage due to this effect is growing during a long time (almost an hour) after the light is switched off. The photoinduced current due to this effect changes its sign with a very long delay after the light is switched off. So, it is ascertained that there are two different effects: APV-1 and APV-2. They have quite different mechanisms.In [55], it was concluded that the APV-1 effect originates due to the transitions of Tb3+ ions from the ground state 7F6 to the excited state 5D4 induced by the laser irradiation. The saturation of the APV-1 effect is reached when the dynamic equilibrium between the numbers of straightforward and backward transitions is established. The lattice accommodation explains the large value of the time constant τ≈5s of the APV-1 effect. The process of the accommodation takes some time. That is why the time constant τ is essentially larger than the life time of a free Tb3+ ion in the excited state 5D4.The APV-2 effect originates from the autoionization of the excitedTb3+ ions and the subsequent diffusion of the knocked out electrons in the electric field of the spontaneous polarization through the crystal. The autoionization can be caused by the internal electrical field related to the ferroelectric order. The autoionization results in appearance of a Tb4+ hole and a nonequilibrium electron. If the light is nonhomogeneous, the diffusion of these carriers from illuminated part of the sample to the nonilluminated one gives rise to the photoinduced current. The sign of the effect is determined by the fact that a pairTb4+-(MoO4)2- appears instead of a pair Tb3+-(MoO4)2-. So, autoionization increases the resulting electrical dipole moment, and the polarization of the sample increases as well.


## 7. Fluctuations Near theβ-β′ Phase Transition

In [63], giant fluctuations of the amplitude of the acoustic resonant peak were observed experimentally in orthorhombic Gd2(MoO4)3 (GMO) near the Curie temperature. The temperature dependence of the intensity of the acoustic main resonant peak was measured in the temperature range (22–165)°C including the Curie temperature, TC=159°C. The low-frequency fluctuations of the amplitude of the resonant current through the sample were observed in the vicinity of the Curie temperature. The sample was cut in the (001) plane. Its dimensions were 7×7×1mm3. The diagonals of the (001) face were parallel to the orthorhombic [100] and [010] axis.The piezo-acoustic resonance was investigated. The scheme recommended in the list of standards on measurements of piezo-electric crystals [64] was used. Acoustic oscillations in the sample were generated using self-exciting method. An oscilloscope was used as a detector of oscillations in the sample. Frequency dependences of the amplitude of the sinusoidal input voltage of the oscilloscope V1 were measured at fixed amplitude of the output voltage of the generator V0 in the frequency range (20–300) kHz at T=(22–165)°C. Under the influence of an alternative voltage, the alternative strains arise in a piezoelectric crystal as a consequence of the inverse piezoe-ffect. Under the influence of these strains, the alternative electric polarization arises in the crystal as a consequence of the direct piezoe-ffect. This alternative polarization is equivalent to appearing of an additional alternative electric current. The dependence of this additional current upon the frequency and the value of the voltage of the generator is described in terms of an equivalent electrical oscillating contour. This contour consists of an active resistance R, an inductivity L, and a capacity C that are connected in series [65, 66]. The capacity C depends upon the elastic and piezoelectric properties of the crystal, its dielectric permittivity, the dimensions of the sample, and the oscillation mode. The inductivity L depends upon the crystal density, its piezoelectric constants and the dimensions of the sample. The resistance R is determined by the dissipation of the oscillation energy due to the internal friction and depends on the viscosity, the piezoelectric constants, and the dimensions of the sample.The equivalent contour, the static resistance of the sampleR1 (the resistance in the absence of the oscillations), and the static capacity of the sample C1 (the capacity in the absence of the oscillations) are connected in parallel. When the generator frequency is equal to the resonant frequency of the equivalent contour, the impedance of the contour decreases and becomes equal to its active resistance R. The decreasing of the sample impedance results in an increase of the current through the sample. This results in an increase of the input voltage of the oscilloscope V1, and the resonant maximum appears on the frequency dependence V1(f) at(40)f=fR=12πLC.
The value of the maximum is given by the following equation:(41)V1[fR(T)]=VR(T)=V0[1+ξ(T)C12]ξ(T)C22+[1+ξ(T)C1(C1+C2)]2.
Here, T is the temperature; C1 the static capacity of the sample; C2 the input capacity of the oscilloscope;(42)ξ(T)=[R(T)]2L(T)C(T)=4π2[fR(T)]2[R(T)]2;R, L, C are the resistance, the inductivity, and the capacity of the equivalent contour, respectively. Equation (41) is derived for the case R1≫1/(ωC1); R2≫1/(ωC2). R2  is the active input resistance of the oscilloscope. It was assumed that the values R, C1, and C2 do not shift noticeably the resonance term ω2LC=1.Figure20 shows the experimental temperature dependence of the resonant peak intensity VR(T).Figure 20
The temperature dependences of the maximum (open circles) and minimum (filled circles) intensities of the fluctuating resonant peakVR. Lines are drawn for the eyes guide.At temperatures higher than≈130°C, the intensity VR(T) displays large and slow fluctuations at constant values of the temperature and the frequency. These fluctuations can be fixed visually.The time of fluctuations reaches tens of seconds. Figure20 shows maximum (open circles) and minimum (filled circles) values of VR. The maximum and the minimum values of VR were found from the observation during two minutes at fixed values of the temperature and the frequency. The variations of VR due to the fluctuations reached 0.1Volt. The fluctuations could be observed only in the limit of the resonant peak width ΔfR≈1kHz. Out of the resonance, they were zero with the accuracy of 0.008VR≈0.0003Volt. Curves are drawn only for the eyes guide.Figure21 shows the time dependence of the resonant peak intensity VR(t) at T=154.25°C and fR=154.9kHz. The oscillating line is a spline. The horizontal line is a mean value of VR over the time of observations. The mean square fluctuation of VR is ≈33 mV.Figure 21
Fluctuations atT=154.25°C and  fR=154.9kHz. The oscillating line is a spline. The horizontal line is a mean value of  VR over the time of observations.It constitutes≈15% of the mean value V̅R≈0.22Volt. The real frequency spectrum of the fluctuations is wider than that of the curve on Figure 21. Rapid fluctuations occur with periods much less than one second. So fast variations of VR(t) can be seen but it is impossible to fix them visually. Observed fluctuations of VR(t) occurs due to the fluctuations of the resonant current through the sample. When the temperature of the sample is near the temperature of the phase transition, various parameters of the crystal fluctuate and lead to the fluctuations of the impedance of the equivalent contour. This results in the fluctuations of the current through the sample. The value VR at the resonance is related to the parameters of the equivalent contour by (41). It is seen from (41) that the value of VR and consequently the values of its fluctuations can be presented as a function of the only variable ξ=R2/LC. It can be shown that ξ depends on four material constants. They are the elastic constant, the dielectric permeability, the coefficient of the internal friction, and the piezoelectric constant. ξ depends also upon the density ρ, but fluctuations of the density cannot reach so large values as shown on Figure 21. The observed fluctuations of the piezoacoustic impedance in GMO exceeds fluctuations of parameters of a crystal reported up to now [67] by the two orders of magnitude.



---

*Source: 276348-2012-05-09.xml*