Review Article  Open Access
B. K. Ponomarev, A. Zhukov, "Magnetic and Magnetoelectric Properties of Rare Earth Molybdates", Physics Research International, vol. 2012, Article ID 276348, 22 pages, 2012. https://doi.org/10.1155/2012/276348
Magnetic and Magnetoelectric Properties of Rare Earth Molybdates
Abstract
We present results on ferroelectric, magnetic, magnetooptical 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 byproduct 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 K (GMO) and 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.
1. Introduction
At the end of 1960s, there was a great interest in the rare earth molybdates family (RMO) ( and ) because these compounds exhibited the phenomena of ferroelectricity and ferroelasticity [1].
RMO from (PMO) to TMO crystallizes to the tetragonal structure with a space group P (point symmetry group is ). Their melting points are 1045°C for PMO and 1172°C for TMO. (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 hightemperature 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 ferroelasticferroelectric 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/cm^{2}. The dielectric permeability peak at the Curie temperature is of ≈10 for GMO.
The ferroelastic properties are well exhibited. In a singledomain 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(MoO_{4})_{3} (DGMO), TbGd(MoO_{4})_{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 field along the [001] and [110] axes at and . The experimental field dependences of the wave numbers of the absorption peaks were obtained. The analysis of the dependences shows that the excited multiplet ^{5} of 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 were obtained from the experimental data in the magnetic field along the [001] axis for the ^{5} multiplet. The experimental dependences were described well by the theory of the magnetism of singlets.
The dependences of the photoinduced voltage in TMO upon the time of the illumination by a laser beam were measured at room temperature at . The power of the laser radiation was . 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 (). The lowfrequency 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 C_{2} for the nearest surrounding of ions in GMO was established experimentally. The calculated value of the spontaneous polarization, assuming point charges and , was 260 · 10^{−9} Coul/cm^{2}. This result can be considered as being in a satisfactory agreement with the experimental value . 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 sites is assumed to be C_{2}. The experimental data discussed below relates mostly to , , and . Table 1 contains the quantum numbers for corresponding free ions in the ground state.

It is worth mentioning that starting from 90th magnetic shapememory alloys (MSMAs) attract special interest owing to significant magneticfieldinduced strain (MFIS), also referred as the “magnetic shapememory effect” originated from coupling between magnetic and structural ordering. Such strain arises through the magneticfieldinduced motion of twin boundaries [3]. Since the magnetic shapememory 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 shapememory effects cause magnetic fieldinduced superelasticity, which is the magnetically induced recovery of a large mechanically induced deformation [4]. It is worth mentioning that, in fact, magneticfieldinduced strain (MFIS) effect has been previously observed and described in the same terms at the beginning of 90th in ferroelectrics (rareearth molybdates) [5]. Particularly, terbium molybdate shows so strong fieldinduced strain, that the spontaneous destruction in a constant magnetic field of along the [100] axis at has been observed [5], but at that time, most attention has been paid to magneticfieldinduced 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 magneticfieldinduced 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 family (RMO) (, and ).
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 ferroelasticferroelectric GMOtype structure. The  transition temperature is also the Curie temperature of the orthorhombic RMO. The spontaneous polarization appears at the transition as a byproduct 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 farinfrared reflectivity and transmission measurements in GMO.
Analogous result was obtained for TMO from neutronscattering experiments [8, 9]. In [6], from the examination of the soft mode, it was shown that the primaryorder parameter, which necessarily shows large fluctuations, near , 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 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 and point group . Below the Curie temperature, it is orthorhombic with space group and point group . The reduction in point group symmetry from to 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 and axes of the tetragonal phase become the [100] and [010] axes of the orthorhombic phase. These orthorhombic axes are unequal. The parameter along the [010] axis is larger than the parameter along the [100] axis: . So, the crystal changes its shape during the transition. The axes [100] and [010] of the tetragonal phase become the axes 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 . The relative change of the [001] axis parameter during the transition is smaller than , 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 between and 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 polarization in GMO was measured in [13] over the temperature range from up to the Curie point . The electric polarization takes values of at and at . The abrupt jump in the curve was observed at with a hysteresis loop of width. This is a direct indication on the phase transition of the first order. Room temperature values of in SMO, EMO, GMO, and TMO are in the range of for EMO and for SMO [14].
The temperature dependences of electric and elastic properties of GMO were measured in [15]. The dielectric permittivity of the clamped crystal does not depend on the temperature and shows no anomaly at the Curie point, while the elastic constant 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 , hardtocontrol 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 .
The presence of the domain walls significantly influences the resonant behavior of a crystal.
3. Magnetic Properties
Magnetic properties of a singlecrystal and singledomain spherical sample of GMO were investigated in [18–22]. In [18–20], the magnetization along the and [100] axes, respectively, was measured at temperatures from 0.3 K to 4.2 K in magnetic field up to . The saturation magnetization values at low temperatures and in high magnetic field along all crystal directions correspond to the magnetic moment which is typical for ion in the ground state ^{8}. At temperatures above ≈1 K, the field dependences of the magnetization show a typically paramagnetic behavior. At , the field dependence of the magnetization exhibits a weak anomaly at H = (23) kOe that indicates on the field induced transformation in the antiferromagnetic structure from the antiparallel to parallel orientation of ionic magnetic moments [18]. At , 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.
Figure 1 shows the field dependences of the magnetization of GMO at 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 that corresponds to a free ion in the ground state ^{8} 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 is saturated at the magnetic field . The saturation magnetization is . This value agrees well with the magnetic moment of ion in the ground state ^{8}. 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 2 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: In the state of the saturation (), the anisotropy energy is ion. This value is higher than the anisotropy constant of ferromagnetic iron at room temperature (). This is the energy needed to align all magnetic moments of ions in the paramagnetic GMO along the hard direction [100] at .
Figure 3 shows the dependence of the anisotropy energy in GMO upon the reduced magnetization along [100]. is derived from by introducing a new variable .
It is seen that this dependence is described well by the relation: Here, 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 . The temperature of the jump was . The magnetic field is the critical field of the destruction of the antiferromagnetic order in GMO at . The magnetization was at the beginning of the jump and ≈30 kG·cm^{3}/mole Gd^{3+} 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 spinflop.
The heat capacity of GMO was measured in the magnetic field up to 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 singlecrystal and singledomain spherical sample of TMO were measured over the range (0.4–4.2) K with fields up to 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 because of the mechanical destruction of the sample.
According to [23], the temperature dependences of the heat capacity at zero field and at show type anomalies with pronounced maxima at and , 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 . All of the above facts indicate on the antiferromagnetic ordering of the magnetic moments along the easy [001] axis with the Neel temperature 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 saturation differs noticeably from for a usual paramagnetic substance with ion in the ground state ^{7}F_{6}. On the basis of this difference, the authors in [23] conclude that the total mechanical moment of the ground state doublet is . 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 of 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 selfdestruction of TMO single crystal. Our later measurements showed that a sample of TMO at certain conditions can survive at T = 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 Figures 2 and 3. The estimate of the anisotropy energy of TMO at and K gives erg/g. From Figure 2, the corresponding value for GMO is ≈10^{5} 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 ^{159}Tb^{3+} becomes very long at high fields and low temperatures. For example, in the magnetic field at the temperature , 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 singlecrystal singledomain samples of TMO was measured in pulsed magnetic fields up to with pulse duration of 0.01 s. The magnetization curves were measured at and 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 , , the magnetization reached and at , , it was . Along the [010] axis, the magnetization was at , , and at , . 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 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 at .
Figure 4 shows the magnetization curves of TMO in the plane (001) at . At the beginning of the pulse of the magnetic field, the [100] axis of the sample was oriented along the field. At , 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 , 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 singlecrystal singledomain sample of TMO was measured in the direct current magnetic field up to in the temperature range T = (4.2–300) K. The magnetization curves at along the directions 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 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.
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 . This can be explained by the temperature dependence of the critical field of the [100][010] interchange and the magnetocaloric effect in pulsed fields. Our measurements of the magnetoelectric effect showed that the critical field decreases from approximately at down to at . The [100][010] interchange of Figure 4 was observed because of increasing the temperature due to the magnetocaloric effect and the corresponding decreasing of .
It is seen that, along the easy direction [001], the saturation magnetization per one terbium ion is not equal to the magnetic moment of a free ion . This is the result of the interaction of the angular momentum of 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: The Hamiltonian of the problem is is the Hamiltonian of the crystal field; , , and are operators of projections of the total mechanical momentum on the , , and 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 free ion in the magnetic field along axis. This is the case of the full degeneracy of the orbital multiplet: From (5) and (6), the Hamiltonian of a free ion is obtained: Here, is the angle in the , plane. The dependent parts of normalized wave functions of for degenerate levels of the multiplet of a free ion are [27] 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 Solutions of (10) give the eigenvalues of the energy of a free ion: The components of the degenerate multiplet (11) possess the field independent magnetic moments: Using (11) and (12) and the statistics of MaxwellBoltzman, one obtains the Brillouin function to describe the field and temperature dependence of the magnetization in a system of free ions: The total mechanical momentum of ion in ^{7}F_{6} ground state is and are the angular and spin moments, respectively. The factor is [28] The corresponding saturation magnetic moment is .
The singlet magnetism theory describes the magnetic properties of ions in the crystal field. It is based on the JahnTeller 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 electrons unavoidably have the singlet ground state [31]. If two lowest singlets are separated by a narrow gap of the order of 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 given by (12). The forms of the combinations are determined by the symmetry of the crystal field and of the operators: , , and .
The point symmetry group of the surrounding of ion in TMO is C_{2} [2]. For this symmetry, the Hamiltonian matrix for the magnetic field along z axis is Here, is the splitting of energy levels of the quasidoublet under consideration.
is a nondiagonal matrix element of the total mechanical momentum operator. It is given by the relation: and are the wave functions of the lower and upper levels of the quasidoublet, respectively. The difference between and is not too large. This result indicates that the main harmonics that form the quasidoublet are and . From here, it follows 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 From (20), one obtains the eigenvalues of the energy of the quasidoublet in the magnetic field : Here, —the energy of the lower level of the quasidoublet; —the energy of the upper level of the quasidoublet. Magnetic moments of the components of the quasidoublet are given by the following relation: 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 .
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 ion in the ground state ^{8}S_{7/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 MaxwellBoltzman 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: This is a twolevel 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 of ion in TMO is an octahedron formed by seven oxygen ions [2]. The point symmetry group of the octahedron can be described as C_{2}. 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 and are zero in the two singlet model. Nonzero values for and appear when taking account of the nearest excited singlet. This is a threelevel model of singlet magnetism theory.
The sites of 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 for ions, one obtains in the threelevel modification of the theory the following expression for the field and the temperature dependence of the magnetization: Here, indicates the direction in the crystal and takes values and [100].
indicates the type of the site. . is a projection of the field on the direction . is a projection of the magnetic moment per one ion on the direction of the magnetic field. are the saturation magnetic moments of ions at sites of two types for main directions in the crystal. —Van Vleck susceptibility along the direction . —Boltzmann constant. and are coefficients of the theory that can be found from energy gaps and between the singlets. is the energy gap between the th excited singlet and the ground singlet for the sites of the first type, is the same for the sites of the second type. and can be taken from the photoluminescence spectrum of ion in TMO [5]: ; . Using these values of and , the coefficients and were found. Then, the magnetization curves were calculated at the temperature for the directions and [001] by fitting expression (24) to the corresponding experimental curves. In these calculations, the coefficients and were used as the fitting parameters. It was found that the values along the [001] axis are identical for the nonequivalent sites. Along the [100] axis, the saturation magnetic moments for nonequivalent sites are different: . For the [010] direction, the opposite nonequality occurs: . This result reflects the difference in the symmetry of the nonequivalent sites of ions. A good agreement with the experiment takes place at for main directions. One can see that the difference between the experimental saturation magnetic moment of TMO and the magnetic moment of free ion in the ground state ^{7}F_{6} is explained completely by the singlet magnetism theory. So, the additional assumption about the ground state of ion in TMO with [23] is not necessary.
The comparison of the magnetization curves of GMO and TMO indicates definitely the role of the angular momentum of the trivalent rare earth ion in the formation of the magnetic properties in rare earth molybdates.
The value determines the main contribution to the energy of the natural magnetic crystalline anisotropy. In GMO, ion in the ground state ^{8}S_{7/2} has , and the energy of anisotropy in the magnetically saturated state is . The magnetization reaches the saturation along hard directions [100] and [010] at (see Figures 1 and 2). In TMO, ion in the ground state ^{7}F_{6} has , and the magnetization along hard directions [100] and [010] is far from saturation even in the magnetic field (see Figure 5). The energies of anisotropy at kOe found using relation (1) are erg⁄g and 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 electrons. Magnetizing the crystal containing nonspherical ions like 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 . 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 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 4felectrons. 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 . The first experimental observations of linear MEE in 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 tetragonal was published in [36].
In [37], nonlinear MEE was observed in paramagnetic TMO at . 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.
Figure 6 shows the field dependence of the magnetically induced electric polarization (MEP) in TMO at [38]. The field was applied along the orthorhombic axis . MEP was measured along the axis. As the spontaneous polarization is directed along the axis, this result can be considered as the rotation of the in the magnetic field.
The dependence of MEP upon the magnetic field is quadratic with a good accuracy. MEP is at . It is small compared to the spontaneous electric polarization [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: For a usual paramagnetic crystal, Here correspond to the , , axes, respectively; 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 ions surrounding. As it was noted in Section 3, the nearest surrounding of ion in TMO is an octahedron formed by seven oxygen ions [2]. The octahedron has the C_{2} 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 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 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 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 individual 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 4felectrons. This charge cloud is coupled strongly with an ionic magnetic moment due to a strong spinorbit 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 piezoeffect giving rise to MEE. The effect is large in those molybdates, in which rare earth ions have large magnetic moments and nonzero angular moments .
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 Tb^{3+} 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 piezoeffect. 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 . This prediction was confirmed experimentally in [38].
Figure 7 shows the experimental dependence in TMO. MEP was measured along the [110] axis in pulsed magnetic field of 104.5 kOe directed along the [] axis. Solid line presents the leastsquare fitting. One can see that the prediction is fulfilled well.
The movement of the ferroelectric domain walls initiated in (TGMO) by the pulsed magnetic field of at was observed in [40, 41]. It was the first observation of the reversal of the spontaneous electric polarization due to the magnetic field.
Figure 8 shows the field dependence of the electric polarization in a singledomain sample of TMO at [42]. The polarization was measured along the 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 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 and by the mutual interchange of and axes. In orthorhombic rare earth molybdates, the ferroelectric domains are identical to the structural twin domains.
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 () 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 .
The process of the displacement of the domain wall starts when the new domain reaches the critical size and begins to grow. In Figure 8, the start of the displacement takes place at . 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 [] 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 , the opposite domains are equal in size (). 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 , 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 Figure 8 showed that after the jump the axes and 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 dependence was measured in [42]. While decreasing the temperature from down to , the critical field decreases from to . Further decrease of the temperature results in an increase of the critical field up to at . At lower temperatures down to , the reversal of the electric polarization was not observed up to .
During the growth of the new domain in TMO, small rotations of the tetrahedrons 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.
Figure 9 shows the field dependence of MEP with a sharp break down of MEP at . 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 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 mm^{3}. The large face of the sample was the (001) plane and the small faces were the (110) and () 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 along the [100] axis without the reversal.
The anisotropy of MEE in TMO was investigated experimentally in [43] in the magnetic field up to at temperatures ; and at up to . 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: Here, is the angle between the magnetic field and the tetragonal axis; is the constant of the theory.
Figure 10 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 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 , curves 1 and 6). MEP is negative in the magnetic field directed along the [100] axis (a negative paraprocess at , curves 5 and 10).
The field dependences do not agree with (28) and what is more unusual is that the itself depends upon the direction of the field.
The angular dependences of MEP in TMO do not agree with (28) neither.
Figure 11 displays some of the experimental angular dependences of MEP in TMO at temperatures (curves 1−3) and (curve 4). Curves 1, 2, 3 were measured at , , and , respectively. Curve 4 was measured at . 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 . Their arrangement is not symmetric relative to contrary to (28). The absolute values of and are different. The values of . 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].
The magneticfieldinduced electric polarization in crystals (MEP) is, like the magnetostriction and magnetic linear birefringence, an evenparity 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 dions. The situation with fions 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 4fshell. Their calculation requires knowledge of the energy levels and wave functions of these ions. It was established in [26] that the ground state of ions in TMO is apparently a combination of three closelying singlets, which belongs to different irreducible representations of the symmetry group of the ion environment and are separated by in energy.
In this case, at the magnetic field in the (001) plane, the following relation takes place [43]: Here, is the temperature. It should not be very high, but the condition should be fulfilled. is the temperature of the antiferromagnetic ordering. In TMO, [23].
, and are constants of the theory which do not depend on the temperature; is the magnetic field component along the [010] axis: is the magnetization component along the [010] axis.
is the magnetic field component along the [100] axis: is the magnetization component along the [100] axis.
and are the Van Vleck susceptibilities along the [010] and [100] axes, respectively. and are the magnetization components of the ion quasitriplet in the (001) plane. The values of the constants , and were determined in [43] using corresponding experimental dependences from [43] and from [26]. The solid lines 1, 2, and 3 on Figure 11 were calculated using (29), the values of , , 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 . At , there is no agreement between the theory and the experiment because the theory does not take account of the antiferromagnetic order that arises at .
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 in (GMO) at 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 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 the value of ≈10^{5} 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 .
At and , MEP was . 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 and . The experimental angular dependences of MEP at and various field values are shown on Figure 12. The solid lines are calculated by (28) shifted along the axis by . One can see that the experimental dependences agree with the relation (28) through the whole field range investigated. The experimental dependences at 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 and satisfactorily.
The experimental field dependences in GMO at and are shown in Figure 13 for two orientations of the field along the [010] axis (filled circles) and along the [100] axis (open circles).
It is seen that, for both directions of the field, the experimental isotherms agree with the relation (28) only at , that is, . At higher fields, the disagreement increases. At , 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 in GMO can be described by the relation analogous to that obtained in [49] for the single ion magnetostriction: Here, is a normalized hyperbolic Bessel function of th order: is the Legendre polynomial. is the inverse Langevin function of the reduced magnetization . For the case considered, has been defined by the secondorder symmetry of the angular dependence obtained from Figure 12. For , the reduced hyperbolic Bessel function is given by the relation: 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 . 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 , open symbols refer to . The experimental values of MEP do not exceed .
Calculated dependences are shown by solid lines in Figure 14 for T = 4.2 K (curves 1 and 2). Coefficients are and for the [100] and [010] axes, respectively. For , the divergence between the experiment and the theory is . So, the relation (32) obtained from the theory of the single ion magnetostriction agrees to the experiment well at . For m = 0.95, . The value in the case of GMO corresponds to the field . 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 near indicates some additional effects that should be considered. The dependences 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 paramagnetismantiferromagnetism 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 in (SMO) was investigated in direct current magnetic field up to at temperatures from to [50]. MEP was measured along the [001] axis in the magnetic field along the [010] axis and [100] axis.
Figure 15 represents a cycle of the MEP variation (abcde) observed in a single domain sample of SMO at 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 . The branch ab was passed for 265 s. It represents the experimental dependence of MEP in the increasing field . At the b point (), the MEP was . At this point, the magnetic field was fixed and remained unvaried for 232 s. During this time, the MEP relaxed to (point c). At point c, the magnetic field started to decrease at a rate . The descending branch cd represents the experimental dependence of MEP in the decreasing field . At point d, the magnetic field reached zero and was fixed. The MEP was found to be at this point. After fixing the magnetic field at zero, the MEP was measured for . During this time, it relaxed to .
This variation of the MEP is plotted as de on Figure 15. 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 . 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: The relaxation time increased with decreasing temperature. It depended upon the magnetic field only weakly.
Table 2 presents the values of obtained at various temperatures in the field .

The 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 and can be approximately described by the expression where is a constant. In this case, a timedependent contribution could be eliminated from the experimental field dependences using the relations: Here, does not depend on the time. The and arrows denote the values obtained under increasing and decreasing field respectively; is the maximum field reached in the measurements; is the rate of the magnetic field variation with the time; is the true MEP at the maximum field calculated from the experimental curve using (36); 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 timedependent contribution from the experimental dependences yields a hysteresisfree 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 Figure 16. The dependences 1 and 6 obtained at 4.4 K are fitted well by the relation (28) at and . At higher fields and lower temperatures, this relation fails.
The data described in Section 4 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.
Table 3 shows the values of MEP in the three rare earth molybdates investigated here at fixed temperature and field values. The values of and were taken from [51]. There is a clear qualitative correlation between the values of the angular momentum and the values of . GMO with zero angular momentum of ion possesses the least value. The angular momentum characterizes the spatial anisotropy of an ion. It is obvious that the rotation of an anisotropic ion with produces more strains in a crystal lattice than the rotation of a spherically symmetrical ion with . So, the change of the electric polarization (MEP = ) is larger for than for . There is some incongruity in this conclusion: the angular momentum of ion is larger than that of ion, but the 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 of ions under consideration. Here, is factor; is the total mechanical moment of ion; is the spin mechanical moment of ion. The sign “−” corresponds to light rare earth ions from to , and the sign “+” corresponds to heavy rare earth ions from to . is the Bohr magneton. One can see that the magnetic moment of ion is larger than that one of ion. So, one should conclude that the total effect of the space anisotropy and the magnetic moment of ion is larger than that one of 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.
5. MagnetoOptical Properties
In [26], the experimental dependences of the magnetization in terbium molybdate  (TMO) upon the magnetic field and the temperature were described satisfactorily using the theory of the magnetism of singlets [31] at low temperatures . At higher temperatures, the agreement was not good. It was due to increasing population of highest levels of ^{7}F_{6} multiplet of ions. For the calculation of the magnetization at high temperatures, the highest levels of ^{7}F_{6} multiplet should be known. To calculate the spectrum of ^{7}F_{6} multiplet, one should know the crystal field acting on 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 ^{5}D_{4} can simplify the calculations of the crystal field parameters in the case under investigation. The spectrum of ^{5}D_{4} 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 wellresolved absorption peaks were observed. It was established that those peaks correspond to the transitions ^{7}F_{6}^{5}D_{4} in ions.
In [53], the absorption spectra were measured in TMO in the nonpolarized light in the magnetic field along the [001] and [110] axes at and . A wave vector of the light was parallel to the magnetic field. At zero field, 14 wellresolved peaks with wave numbers () 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 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 field upon the direction of the field. The analysis of these dependences showed that the excited multiplet ^{5}D_{4} of 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 of the excited multiplet ^{5}D_{4} of ion in the crystal field with the point group symmetry C_{2} was found for the magnetic field up to 200 kOe along the [001] axis. The experimental dependences 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 for 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 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 . 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/cm^{2}. 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 length induced in TMO the luminescence in the green part of the visible spectrum with the wave length ≈. The quantum output constituted 20%.
In [55], the dependences of the photoinduced voltage in TMO upon the time of the illumination by a laser beam were measured at room temperature at . The power of the laser radiation was . 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 that was directed opposite to the spontaneous electric polarization . This effect was observed at a small duration of illumination . It was called APV1. The effect that was observed in [54] was called APV2. The sign of APV1 was opposite to the sign of APV2. The voltage and the relaxation time of APV1 were essentially smaller than those of APV2. That is why the APV1 effect was not observed on the background of the APV2 effect in [54].
The first observation of APVeffect was reported for a single crystal sample of ferroelectric [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: The effect APV1 is well described by the expressions (38) and (39). The effect APV2 does not follow these expressions.
Figure 17 shows experimental time dependences of the photoinduced current (curve 1) and the photoinduced voltage (curve 2).
These curves were obtained on the sample number 1 with size . 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 the face of the sample by a cylinder lens to the spot with dimensions . 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/cm^{2}.
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 ; . The stepped form of the curve 1 is due to the resolution of the equipment. When the light is switched on, the negative voltage appears. This voltage is due to the photoinduced electric polarization that is directed opposite to the spontaneous electric polarization . The value corresponds to the decreasing of by the value . When the light is switched off, the photoinduced voltage decreases exponentially with the time constant . The photoinduced voltage and the photoinduced current on Figure 17 are due to the effect APV1.
Figure 18 shows the dependences 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 , the value reaches a minimum of and begins to increase with increasing .
At , the sign of is changed. The positive sign of at corresponds to the coincidence of the directions of and . The tendency to the saturation of is seen at large values of . 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 aftereffect takes place. In other words, the long negative pulse of the polarization due to the effect APV1 exists on the background of slow positive drift of the polarization. This drift is due to the effect APV2.
Curve 2 in Figure 18 was measured on the sample number 1 with a homogeneous distribution of the intensity across the face of the sample. The crosssection of the laser beam had a circular form with . 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 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 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 . The relative directions of the crystal axes and of the laser beam were the same as for sample number 1. The crosssection 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 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 APV2 is so small that one can neglect it.
Figure 19 shows curve 1 from Figure 18 contracted along the axis and its prolongation for large 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 at and then decreases.
The value of corresponds to the photopolarization . The relaxation time is . It is more than the time constant of the input circuit of the electrometer . Due to this, the dependence 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 APV2 effect for the present experiment was achieved at . Due to the relation between and , this value of should be considered as a bottom estimate.
The two effects described above show quite different dependences upon the conditions of the experiment.
The APV1 effect exists both at homogeneous and nonhomogeneous illuminations (Figure 18). It reaches the saturation at (Figure 17). The APV2 effect exists only at nonhomogeneous illumination.
It was not observed at a homogeneous illumination (Figure 18, curve 3). It depends on the time and the geometry of the beam crosssection (Figure 18).
The time constant of the APV1 effect is . From Figures 17 and 18, one can see that APV1 effect reaches the saturation in a relatively short time and after this the voltage due to APV1 effect is constant. Contrary to APV1 effect, the voltage due to APV2 effect grows for a very long time of the order of several thousand seconds showing only very weak tendency to the saturation.
The APV2 effect possesses the aftereffect: the value increases during several thousand seconds after the light is switched off (Figure 19). No aftereffect was observed while measuring the APV1 effect.
The APV1 effect has a feature that is common to photoinduced 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 APV2 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: APV1 and APV2. They have quite different mechanisms.
In [55], it was concluded that the APV1 effect originates due to the transitions of ions from the ground state ^{7}F_{6} to the excited state ^{5}D_{4} induced by the laser irradiation. The saturation of the APV1 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 of the APV1 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 ion in the excited state ^{5}D_{4}.
The APV2 effect originates from the autoionization of the excited 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 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 pair appears instead of a pair . 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 (GMO) near the Curie temperature. The temperature dependence of the intensity of the acoustic main resonant peak was measured in the temperature range including the Curie temperature, . The lowfrequency 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 mm^{3}. The diagonals of the (001) face were parallel to the orthorhombic [100] and [010] axis.
The piezoacoustic resonance was investigated. The scheme recommended in the list of standards on measurements of piezoelectric crystals [64] was used. Acoustic oscillations in the sample were generated using selfexciting 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 were measured at fixed amplitude of the output voltage of the generator in the frequency range (20–300) kHz at . Under the influence of an alternative voltage, the alternative strains arise in a piezoelectric crystal as a consequence of the inverse piezoeffect. Under the influence of these strains, the alternative electric polarization arises in the crystal as a consequence of the direct piezoeffect. 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 , an inductivity , and a capacity that are connected in series [65, 66]. The capacity depends upon the elastic and piezoelectric properties of the crystal, its dielectric permittivity, the dimensions of the sample, and the oscillation mode. The inductivity depends upon the crystal density, its piezoelectric constants and the dimensions of the sample. The resistance 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 sample (the resistance in the absence of the oscillations), and the static capacity of the sample (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 . 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 , and the resonant maximum appears on the frequency dependence at The value of the maximum is given by the following equation: Here, is the temperature; the static capacity of the sample; the input capacity of the oscilloscope;, , are the resistance, the inductivity, and the capacity of the equivalent contour, respectively. Equation (41) is derived for the case ; . is the active input resistance of the oscilloscope. It was assumed that the values , , and do not shift noticeably the resonance term .
Figure 20 shows the experimental temperature dependence of the resonant peak intensity .
At temperatures higher than ≈130°C, the intensity 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. Figure 20 shows maximum (open circles) and minimum (filled circles) values of . The maximum and the minimum values of were found from the observation during two minutes at fixed values of the temperature and the frequency. The variations of due to the fluctuations reached . The fluctuations could be observed only in the limit of the resonant peak width . Out of the resonance, they were zero with the accuracy of . Curves are drawn only for the eyes guide.
Figure 21 shows the time dependence of the resonant peak intensity at and . The oscillating line is a spline. The horizontal line is a mean value of over the time of observations. The mean square fluctuation of is ≈33 mV.
It constitutes ≈15% of the mean value . 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 can be seen but it is impossible to fix them visually. Observed fluctuations of 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 at the resonance is related to the parameters of the equivalent contour by (41). It is seen from (41) that the value of and consequently the values of its fluctuations can be presented as a function of the only variable . 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.
Acknowledgments
A. Zhukov wishes to underline that most scientific contribution in the field of RMO belongs to Prof. B. K. Ponomarev, who recently passed away. Consequently, the role of A. Zhukov, apart of magnetic measurements performed in 1990s, was mostly the preparation of this publication. This work was supported by EU ERANET Program under Project “SoMaMicSens” (MANUNET2010Basque3), by Spanish Ministry of Science and Innovation, MICINN under Project MAT201018914 and by the Basque Government under Saiotek 11 MIMAGURA Project (SPE11UN087). A. Zhukov wishes to acknowledge the support of the Basque Government under Program of Mobility of the Investigating Personnel of the Department of Education, Universities and Investigation, Grant nos. MV2009221 and MV2010231.
References
 L. H. Brixner, J. R. Barkley, and W. Jeitschko, “Rareearth molibdates (VI),” in Handbook on the Physics and Chemistry of Rare Earths, K. A. Gschneidner Jr. and L. Eyring, Eds., vol. chapter 30, pp. 610–655, NorthHolland Publishing Company, 1979. View at: Google Scholar
 E. T. Keve, S. C. Abrahams, and J. L. Bernstein, “Ferroelectric ferroelastic paramagnetic betaGd_{2}(MoO_{4})_{3} crystal structure of the transitionmetal molybdates and tungstates. VI,” The Journal of Chemical Physics, vol. 54, no. 7, pp. 3195–3206, 1971. View at: Google Scholar
 K. Ullakko, J. K. Huang, C. Kantner, R. C. O'Handley, and V. V. Kokorin, “Large magneticfieldinduced strains in Ni2MnGa single crystals,” Applied Physics Letters, vol. 69, no. 13, pp. 1966–1968, 1996. View at: Google Scholar
 D. C. Dunand and P. Müllner, “Size effects on magnetic actuation in NiMnGa shapememory alloys,” Advanced Materials, vol. 23, no. 2, pp. 216–232, 2011. View at: Publisher Site  Google Scholar
 Boris K. Ponomarev, Sergey A. Ivanov, Valery D. Negrii et al., “Giant magnetic anisotropy in paramagnetic Tb_{2}(MoO_{4})_{3},” Ferroelectrics, vol. 151, no. 1–4, pp. 103–108, 1994. View at: Google Scholar
 A. P. Levanyuk and D. G. Sannikov, “On closed by the temperature phase transitions of the second type,” JETP Letters, vol. 55, pp. 256–265, 1968. View at: Google Scholar
 J. Petzelt and V. Dvořak, “New type of ferroelectric soft mode in gadolinium molybdate,” Physica Status Solidi (b), vol. 46, pp. 413–423, 1971. View at: Google Scholar
 J. D. Axe, B. Dorner, and G. Shirane, “Mechanism of the ferroelectric phase transformation in rareearth molybdates,” Physical Review Letters, vol. 26, no. 9, pp. 519–523, 1971. View at: Publisher Site  Google Scholar
 B. Dorner, J. D. Axe, and G. Shirane, “Neutronscattering study of the ferroelectric phase transformation in Tb_{2}(MoO_{4})_{3},” Physical Review B, vol. 6, no. 5, pp. 1950–1963, 1972. View at: Publisher Site  Google Scholar
 W. Jeitschko, “The crystal structure of ferroelectric gadolinium molybdate, Gd_{2}(MoO_{4})_{3},” Die Naturwissenschaften, vol. 57, no. 11, p. 544, 1970. View at: Publisher Site  Google Scholar
 W. Jeitschko, “A comprehensive Xray study of the ferroelectricferroelastic and paraelectricparaelastic phases of Gd_{2}(MoO_{4})_{3},” Acta Crystallographica B, vol. 28, pp. 60–76, 1972. View at: Google Scholar
 K. Aizu, “Possible species of "Ferroelastic" crystals and of simultaneously ferroelectric and ferroelastic Crystals,” Journal of the Physical Society of Japan, vol. 27, no. 2, pp. 387–396, 1969. View at: Google Scholar
 E. Sawaguchi and L. E. Cross, “Spontaneous polarization of Gd_{2}(MoO_{4})_{3},” Journal of Applied Physics, vol. 44, no. 6, pp. 2541–2544, 1973. View at: Publisher Site  Google Scholar
 H. J. Borhardt and P. E. Bierstedt, “Ferroelectric rareearth molybdates,” Journal of Applied Physics, vol. 38, no. 5, pp. 2057–2061, 1967. View at: Google Scholar
 L. E. Cross, A. Fouskova, and S. E. Cummins, “Gadolinium molybdate, a new type of ferroelectric crystal,” Physical Review Letters, vol. 21, no. 12, pp. 812–814, 1968. View at: Publisher Site  Google Scholar
 F. G. Ullman, B. N. Ganguly, and J. R. Zeidler, “Pyroelectric detection properties of gadolinium molybdate (gmo),” Journal of Electronic Materials, vol. 1, no. 3, pp. 425–434, 1972. View at: Publisher Site  Google Scholar
 F. G. Ullman, K. M. Cheung, G. A. Rakes, and B. N. Ganguly, “Reexamination of the pyroelectric current anomaly in gadolinium molybdatet,” Ferroelectrics, vol. 9, pp. 63–64, 1975. View at: Google Scholar
 R. A. Fisher, E. W. Hornung, G. E. Brodale, and W. F. Giauque, “Magnetothermodynamics of antiferromagnetic, ferroelectric βGd_{2}(M.oO_{4})_{3}. I. Heat capacity, entropy, magnetic moment of the electrically polarized form from 0.4 to 4.2°K with Fields to 90 kG along the c crystal axis,” The Journal of Chemical Physics, vol. 56, no. 1, pp. 193–212, 1972. View at: Google Scholar
 E. W. Hornung, G. E. Brodale, R. A. Fisher, and W. F. Giauque, “Magnetothermodynamics of antiferromagnetic, ferroelectric βGd_{2}(MoO_{4})_{3}. II. Heat capacity, entropy, magnetic moment of the electrically polarized form from 0.4 to 4.2°K with fields to 90 kG along the a crystal axis,” The Journal of Chemical Physics, vol. 56, no. 10, pp. 5007–5018, 1972. View at: Google Scholar
 G. E. Brodale, R. A. Fisher, E. W. Hornung, and W. F. Giauque, “Magnetothermodynamics of antiferromagnetic, ferroelectric βGd_{2}(MoO_{4})_{3}. III. Heat capacity, entropy, magnetic moment of the electrically polarized form from 0.4 to 4.2°K with fields to 90 kG along the b crystal axis,” The Journal of Chemical Physics, vol. 56, no. 12, pp. 6118–6125, 1972. View at: Google Scholar
 R. A. Fisher, E. W. Hornung, G. E. Brodale, and W. F. Glauque, “Magnetothermodynamics of antiferromagnetic, ferroelectric, ferroelastic βGd_{2}(MoO_{4})_{3}, IV Thermodynamic temperature and other properties without heat introduction below 0.5°K. Fields to 10 kG along the c^{+} crystal axis,” The Journal of Chemical Physics, vol. 59, no. 11, pp. 5796–5809, 1973. View at: Google Scholar
 R. A. Fisher, E. W. Hornung, G. E. Brodale, and W. F. Giauque, “Magnetothermodynamics of antiferromagnetic, polarized ferroelectric, ferroelastic βGd_{2}(MoO_{4})_{3}. V. Thermodynamic temperature and other properties with heat introduction below 0.5°K. Fields to 5 kG along the b crystal axis,” The Journal of Chemical Physics, vol. 69, no. 6, pp. 2892–2900, 1978. View at: Google Scholar
 R. A. Fisher, E. W. Hornung, G. E. Brodale, and W. F. Giauque, “Magnetothermodynamics of ferroelectric, ferroelastic, antiferromagnetic βterbium molybdate. I. Heat capacity, entropy, magnetic moment of the electrically polarized form from 0.4 to 4.2 °K with fields to 90 kG along the c crystal axis,” The Journal of Chemical Physics, vol. 63, no. 3, pp. 1295–1308, 1975. View at: Google Scholar
 B. K. Ponomarev, Yu. F. Popov, and D. S. Red'kin, “Magnetostriction of paramagnetic terbium molybdate in fields up to 150 kOe,” JETP Letters, vol. 57, pp. 483–486, 1993. View at: Google Scholar
 B. K. Ponomarev, S. A. Ivanov, B. S. Red'Kin, and V. N. Kurlov, “Irreversible alterations of ferroelectric domain structure in paramagnetic rare earth molybdates induced by a magnetic field,” Journal of Applied Physics, vol. 75, no. 12, pp. 8004–8007, 1994. View at: Publisher Site  Google Scholar
 B. K. Ponomarev, A. I. Popov, J. Van Tol et al., “Magnetism of singlets in terbium molybdate,” Journal of Magnetism and Magnetic Materials, vol. 258259, pp. 510–512, 2003. View at: Publisher Site  Google Scholar
 L. D. Landau and E. M. Lifshits, Kvantovaya Mekhanika, Nauka, Moscow, Russia, 1974.
 J. Samuel Smart, Effective Field Theories of Magnetism, IBM Watson Research Center, New York, NY, USA, 1966.
 H. A. Jahn and E. Teller, “Stability of polyatomic molecules in degenerate electronic states. I. Orbital degeneracy,” Proceedings of the Royal Society A, vol. 161, pp. 220–235, 1937. View at: Google Scholar
 H. A. Kramers, “Theorie generale de la rotation paramagnetique dans les cristaux,” Proceedings Academy of Science, Amsterdam, vol. 33, pp. 959–972, 1930. View at: Google Scholar
 A. K. Zvezdin, V. M. Matveev, A. A. Mukhin, and A. I. Popov, Redkozemel'nye Iony v MagnitoUporyadochennykh Kistallakh, Nauka, Moscow, Russia, 1985.
 L. D. Landau and E. M. Lifshits, Electrodynamics of Continuum, Gostekhizdat, Moscow, Russia, 1957.
 I. E. Dzyaloshinskii, “Topological model of a heisenberg spin glass,” JETP Letters, vol. 37, pp. 881–882, 1959. View at: Google Scholar
 D. N. Astrov, “The magnetoelectric effect in antiferromagnetics,” JETP Letters, vol. 38, pp. 984–985, 1960. View at: Google Scholar
 D. N. Astrov, “Magnetoelectric effect in chromium oxide,” JETP Letters, vol. 40, pp. 1035–1041, 1961. View at: Google Scholar
 S. L. Hou and N. Bloembergen, “Paramagnetoelectric effects in NiSO_{4} 6H_{2}O,” Physical Review, vol. 138, no. 4A, pp. A1218–A1226, 1965. View at: Publisher Site  Google Scholar
 S. A. Ivanov, V. N. Kurlov, B. K. Ponomarev, and B. S. Red'kin, “Magnetoelectric effect in terbium molibdate,” JETP Letters, vol. 52, no. 7, pp. 1003–1005, 1990. View at: Google Scholar
 B. K. Ponomarev, S. A. Ivanov, B. S. Red'kin, and V. N. Kurlov, “Magnetoelectrical effect in paramagnetic rareearth molybdates,” Physica B, vol. 177, no. 1–4, pp. 327–329, 1992. View at: Google Scholar
 I. E. Chupis, “Effect of magnetic field on ferroelectric properties of rareearthe molobdates,” Fizika Nizkikh Temperatur, vol. 21, no. 9, pp. 941–4695, 1995. View at: Google Scholar
 B. K. Ponomarev, S. A. Ivanov, B. S. Red'kin, and V. N. Kurlov, “Effect of a magnetic field on the ferroelectric domain structure of TbGd(MoO_{4})_{3},” JETP Letters, vol. 55, pp. 353–361, 1992. View at: Google Scholar
 B. K. Ponomarev, V. D. Negrii, B. S. Red'kin, and Yu F. Popov, “Magnetoelectroelastic effects in some rare earth molybdates and related properties,” Journal of Physics D: Applied Physics, vol. 27, no. 10, pp. 1995–2001, 1994. View at: Publisher Site  Google Scholar
 H. Wiegelmann, B. K. Ponomarev, J. van Tol, A. G. M. Jansen, P. Wyder, and B. S. Red'kin, “Magnetoelectric properties of ferroelectric rare earth molybdates,” Ferroelectrics, vol. 183, pp. 195–204, 1996. View at: Google Scholar
 B. K. Ponomarev, A. I. Popov, E. Steep et al., Fizika Tverdogo Tela, vol. 47, no. 7, p. 1326, 2005.
 N. F. Vedernikov, A. K Zvezdin, R. Z. Levitin, and A. I. Popov, “Magnetic linear birefringence of rareearth garnets,” JETP Letters, vol. 93, p. 2161, 1987, Soviet Physics—JETP, vol. 66, no. 6, p. 1233, 1987. View at: Google Scholar
 N. S. Akulov, Ferromagnetics, Gostekhizdat, Moscow, Russia, 1939.
 E. Callen and H. B. Callen, “Magnetostriction, forced magnetostriction, and anomalous thermal expansion in ferromagnets,” Physical Review, vol. 139, no. 2A, pp. A455–A471, 1965. View at: Publisher Site  Google Scholar
 B. K. Ponomarev, B. S. Red'kin, H. Wiegelmann, A. G. M. Jansen, P. Wyder, and J. van Tol., “Magnetoelectric effect in orthorhombic Gd_{2}(MoO_{4})_{3},” Ferroelectric letters, vol. 18, no. 34, pp. 133–140, 1994. View at: Google Scholar
 B. K. Ponomarev, E. Steep, H. Wiegelmann, A. G. M. Jansen, P. Wyder, and B .S. Red'kin, “Anisotropy of the magnetoelectric effect in ${\beta}^{\prime}$Gd_{2}(MoO_{4})_{3},” Physics of the Solid State, vol. 42, no. 4, pp. 716–738, 2000. View at: Google Scholar
 H. B. Callen and E. Callen, “The present status of the temperature dependence of magnetocrystalline anisotropy, and the l(l+1) 2 power law,” Journal of Physics and Chemistry of Solids, vol. 27, no. 8, pp. 1271–1285, 1966. View at: G