Voigt effect

Last updated
Schematic of the polar Kerr effect, longitudinal Kerr effect and the Voigt effect EffetMOenanglais.png
Schematic of the polar Kerr effect, longitudinal Kerr effect and the Voigt effect

The Voigt effect is a magneto-optical phenomenon which rotates and elliptizes linearly polarised light sent into an optically active medium. [1] The effect is named after the German scientist Woldemar Voigt who discovered it in vapors. Unlike many other magneto-optical effects such as the Kerr or Faraday effect which are linearly proportional to the magnetization (or to the applied magnetic field for a non magnetized material), the Voigt effect is proportional to the square of the magnetization (or square of the magnetic field) and can be seen experimentally at normal incidence. There are also other denominations for this effect, used interchangeably in the modern scientific literature: the Cotton–Mouton effect (in reference to French scientists Aimé Cotton and Henri Mouton who discovered the same effect in liquids a few years later) and magnetic-linear birefringence , with the latter reflecting the physical meaning of the effect.

Contents

For an electromagnetic incident wave linearly polarized and an in-plane polarized sample , the expression of the rotation in reflection geometry is is:

and in the transmission geometry:

where is the difference of refraction indices depending on the Voigt parameter (same as for the Kerr effect), the material refraction indices and the parameter responsible of the Voigt effect and so proportional to the or in the case of a paramagnetic material.

Detailed calculation and an illustration are given in sections below.

Theory

Framework and coordinate system for the derivation of Voigt effect.
E
-
i
{\displaystyle {\vec {E}}_{i}}
,
E
-
r
{\displaystyle {\vec {E}}_{r}}
and
E
-
t
{\displaystyle {\vec {E}}_{t}}
are referring to the incident, reflected and transmitted electromagnetic field. Reperefig2.png
Framework and coordinate system for the derivation of Voigt effect. , and are referring to the incident, reflected and transmitted electromagnetic field.

As with the other magneto-optical effects, the theory is developed in a standard way with the use of an effective dielectric tensor from which one calculates systems eigenvalues and eigenvectors. As usual, from this tensor, magneto-optical phenomena are described mainly by the off-diagonal elements.

Here, one considers an incident polarisation propagating in the z direction: the electric field and a homogenously in-plane magnetized sample where is counted from the [100] crystallographic direction. The aim is to calculate where is the rotation of polarization due to the coupling of the light with the magnetization. Let us notice that is experimentally a small quantity of the order of mrad. is the reduced magnetization vector defined by , the magnetization at saturation. We emphasized with the fact that it is because the light propagation vector is perpendicular to the magnetization plane that it is possible to see the Voigt effect.

Dielectric tensor

Following the notation of Hubert, [2] the generalized dielectric cubic tensor take the following form:

where is the material dielectric constant, the Voigt parameter, and two cubic constants describing magneto-optical effect depending on . is the reduce . Calculation is made in the spherical approximation with . At the present moment,[ when? ] there is no evidence that this approximation is not valid, as the observation of Voigt effect is rare because it is extremely small with respect to the Kerr effect.

Eigenvalues and eigenvectors

To calculate the eigenvalues and eigenvectors, we consider the propagation equation derived from the Maxwell equations, with the convention :

When the magnetization is perpendicular to the propagation wavevector, on the contrary to the Kerr effect, may have all his three components equals to zero making calculations rather more complicated and making Fresnels equations no longer valid. A way to simplify the problem consists to use the electric field displacement vector . Since and we have . The inconvenient is to deal with the inverse dielectric tensor which may be complicated to handle. Here calculation is made in the general case which is complicated mathematically to handle, however one can follow easily the demonstration by considering .

Eigenvalues and eigenvectors are found by solving the propagation equation on which gives the following system of equation:

where represents the inverse element of the dielectric tensor , and . After a straightforward calculation of the system's determinant, one has to make a developement on 2nd order in and first order of . This led to the two eigenvalues corresponding the two refraction indices:

The corresponding eigenvectors for and for are:

Reflection geometry

Continuity relation

Knowing the eigenvectors and eigenvalues inside the material, one have to calculate the reflected electromagnetic vector usually detected in experiments. We use the continuity equations for and where is the induction defined from Maxwell's equations by . Inside the medium, the electromagnetic field is decomposed on the derived eigenvectors . The system of equation to solve is:

The solution of this system of equation are:

Calculation of rotation angle

The rotation angle and the ellipticity angle are defined from the ratio with the two following formulae:

where and represent the real and imaginary part of . Using the two previously calculated components, one obtains:

This gives for the Voigt rotation:

which can also be rewritten in the case of , , and real:

where is the difference of refraction indices. Consequently, one obtains something proportional to and which depends on the incident linear polarisation. For proper no Voigt rotation can be observed. is proportional to the square of the magnetization since and .

Transmission geometry

The calculation of the rotation of the Voigt effect in transmission is in principle equivalent to the one of the Faraday effect. In practice, this configuration is not used in general for ferromagnetic samples since the absorption length is weak in this kind of material. However, the use of transmission geometry is more common for paramagnetic liquid or cristal where the light can travel easily inside the material.

The calculation for a paramagnetic material is exactly the same with respect to a ferromagnetic one, except that the magnetization is replaced by a field ( in or ). For convenience, the field will be added at the end of calculation in the magneto-optical parameters.

Consider the transmitted electromagnetic waves propagating in a medium of length L. From equation (5), one obtains for and :

At the position z = L, the expression of is:

where and are the eigenvectors previously calculated, and is the difference for the two refraction indices. The rotation is then calculated from the ratio , with development in first order in and second order in . This gives:

Again we obtain something proportional to and , the light propagation length. Let us notice that is proportional to in the same way with respect to the geometry in reflexion for the magnetization. In order to extract the Voigt rotation, we consider , and real. Then we have to calculate the real part of (14). The resulting expression is then inserted in (8). In the approximation of no absorption, one obtains for the Voigt rotation in transmission geometry:

Illustration of Voigt effect in GaMnAs

Fig 1 : a) Experimental hysteresis cycle on a planar (Ga,Mn)As sample b) Voigt hysteresis cycle obtained by extracting the symmetric part of (a). c) Longitudinal Kerr obtained by extracting the asymmetric part of (a) Articlewikipediakerrvoigtethysteresis.png
Fig 1 : a) Experimental hysteresis cycle on a planar (Ga,Mn)As sample b) Voigt hysteresis cycle obtained by extracting the symmetric part of (a). c) Longitudinal Kerr obtained by extracting the asymmetric part of (a)
Fig 2 : a) Switching mechanism of an in-plane (Ga,Mn)As sample for a magnetic field applied along the [1-10] axis at 12 K. b) Voigt signal simulated from the mechanism showed in a) Articlewikipediasimulationhysteresis.png
Fig 2 : a) Switching mechanism of an in-plane (Ga,Mn)As sample for a magnetic field applied along the [1-10] axis at 12 K. b) Voigt signal simulated from the mechanism showed in a)

As an illustration of the application of the Voigt effect, we give an example in the magnetic semiconductor (Ga,Mn)As where a large Voigt effect was observed. [3] At low temperatures (in general for ) for a material with an in-plane magnetization, (Ga,Mn)As exhibits a biaxial anisotropy with the magnetization aligned along (or close to) <100> directions.

A typical hysteresis cycle containing the Voigt effect is shown in figure 1. This cycle was obtained by sending a linearly polarized light along the [110] direction with an incident angle of approximately 3° (more details can be found in [4] ), and measuring the rotation due to magneto-optical effects of the reflected light beam. In contrast to the common longitudinal/polar Kerr effect, the hysteresis cycle is even with respect to the magnetization, which is a signature of the Voigt effect. This cycle was obtained with a light incidence very close to normal, and it also exhibits a small odd part; a correct treatment has to be carried out in order to extract the symmetric part of the hysteresis corresponding to the Voigt effect, and the asymmetric part corresponding to the longitudinal Kerr effect.

In the case of the hysteresis presented here, the field was applied along the [1-10] direction. The switching mechanism is as follows:

  1. We start with a high negative field and the magnetization is close to the [-1-10] direction at position 1.
  2. The magnetic field is decreasing leading to a coherent magnetization rotation from 1 to 2
  3. At positive field, the magnetization switch brutally from 2 to 3 by nucleation and propagation of magnetic domains giving a first coercive field named here
  4. The magnetization stay close to the state 3 while rotating coherently to the state 4, closer from the applied field direction.
  5. Again the magnetization switches abruptly from 4 to 5 by nucleation and propagation of magnetic domains. This switching is due to the fact that the final equilibrium position is closer from the state 5 with respect to the state 4 (and so his magnetic energy is lower). This gives another coercive field named
  6. Finally the magnetization rotates coherently from the state 5 to the state 6.

The simulation of this scenario is given in the figure 2, with

As one can see, the simulated hysteresis is qualitatively the same with respect to the experimental one. Notice that the amplitude at or are approximately twice of

See also

Related Research Articles

<span class="mw-page-title-main">Pauli matrices</span> Matrices important in quantum mechanics and the study of spin

In mathematical physics and mathematics, the Pauli matrices are a set of three 2 × 2 complex matrices which are Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

<span class="mw-page-title-main">Astronomical coordinate systems</span> System for specifying positions of celestial objects

In astronomy, coordinate systems are used for specifying positions of celestial objects relative to a given reference frame, based on physical reference points available to a situated observer. Coordinate systems in astronomy can specify an object's position in three-dimensional space or plot merely its direction on a celestial sphere, if the object's distance is unknown or trivial.

<span class="mw-page-title-main">Bremsstrahlung</span> Electromagnetic radiation due to deceleration of charged particles

In particle physics, bremsstrahlung is electromagnetic radiation produced by the deceleration of a charged particle when deflected by another charged particle, typically an electron by an atomic nucleus. The moving particle loses kinetic energy, which is converted into radiation, thus satisfying the law of conservation of energy. The term is also used to refer to the process of producing the radiation. Bremsstrahlung has a continuous spectrum, which becomes more intense and whose peak intensity shifts toward higher frequencies as the change of the energy of the decelerated particles increases.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space under the operation of composition.

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.

<span class="mw-page-title-main">Propagator</span> Function in quantum field theory showing probability amplitudes of moving particles

In quantum mechanics and quantum field theory, the propagator is a function that specifies the probability amplitude for a particle to travel from one place to another in a given period of time, or to travel with a certain energy and momentum. In Feynman diagrams, which serve to calculate the rate of collisions in quantum field theory, virtual particles contribute their propagator to the rate of the scattering event described by the respective diagram. These may also be viewed as the inverse of the wave operator appropriate to the particle, and are, therefore, often called (causal) Green's functions.

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.

The Havriliak–Negami relaxation is an empirical modification of the Debye relaxation model in electromagnetism. Unlike the Debye model, the Havriliak–Negami relaxation accounts for the asymmetry and broadness of the dielectric dispersion curve. The model was first used to describe the dielectric relaxation of some polymers, by adding two exponential parameters to the Debye equation:

In general relativity, the Gibbons–Hawking–York boundary term is a term that needs to be added to the Einstein–Hilbert action when the underlying spacetime manifold has a boundary.

The Newman–Penrose (NP) formalism is a set of notation developed by Ezra T. Newman and Roger Penrose for general relativity (GR). Their notation is an effort to treat general relativity in terms of spinor notation, which introduces complex forms of the usual variables used in GR. The NP formalism is itself a special case of the tetrad formalism, where the tensors of the theory are projected onto a complete vector basis at each point in spacetime. Usually this vector basis is chosen to reflect some symmetry of the spacetime, leading to simplified expressions for physical observables. In the case of the NP formalism, the vector basis chosen is a null tetrad: a set of four null vectors—two real, and a complex-conjugate pair. The two real members often asymptotically point radially inward and radially outward, and the formalism is well adapted to treatment of the propagation of radiation in curved spacetime. The Weyl scalars, derived from the Weyl tensor, are often used. In particular, it can be shown that one of these scalars— in the appropriate frame—encodes the outgoing gravitational radiation of an asymptotically flat system.

The history of Lorentz transformations comprises the development of linear transformations forming the Lorentz group or Poincaré group preserving the Lorentz interval and the Minkowski inner product .

The Frank–Tamm formula yields the amount of Cherenkov radiation emitted on a given frequency as a charged particle moves through a medium at superluminal velocity. It is named for Russian physicists Ilya Frank and Igor Tamm who developed the theory of the Cherenkov effect in 1937, for which they were awarded a Nobel Prize in Physics in 1958.

<span class="mw-page-title-main">Dual quaternion</span> Eight-dimensional algebra over the real numbers

In mathematics, the dual quaternions are an 8-dimensional real algebra isomorphic to the tensor product of the quaternions and the dual numbers. Thus, they may be constructed in the same way as the quaternions, except using dual numbers instead of real numbers as coefficients. A dual quaternion can be represented in the form A + εB, where A and B are ordinary quaternions and ε is the dual unit, which satisfies ε2 = 0 and commutes with every element of the algebra. Unlike quaternions, the dual quaternions do not form a division algebra.

<span class="mw-page-title-main">Spinodal decomposition</span> Mechanism of spontaneous phase separation

Spinodal decomposition is a mechanism by which a single thermodynamic phase spontaneously separates into two phases. Decomposition occurs when there is no thermodynamic barrier to phase separation. As a result, phase separation via decomposition does not require the nucleation events resulting from thermodynamic fluctuations, which normally trigger phase separation.

<span class="mw-page-title-main">Gravitational lensing formalism</span>

In general relativity, a point mass deflects a light ray with impact parameter by an angle approximately equal to

<span class="mw-page-title-main">Geographical distance</span> Distance measured along the surface of the Earth

Geographical distance or geodetic distance is the distance measured along the surface of the Earth, or the shortest arch length.

<span class="mw-page-title-main">Two-body Dirac equations</span> Quantum field theory equations

In quantum field theory, and in the significant subfields of quantum electrodynamics (QED) and quantum chromodynamics (QCD), the two-body Dirac equations (TBDE) of constraint dynamics provide a three-dimensional yet manifestly covariant reformulation of the Bethe–Salpeter equation for two spin-1/2 particles. Such a reformulation is necessary since without it, as shown by Nakanishi, the Bethe–Salpeter equation possesses negative-norm solutions arising from the presence of an essentially relativistic degree of freedom, the relative time. These "ghost" states have spoiled the naive interpretation of the Bethe–Salpeter equation as a quantum mechanical wave equation. The two-body Dirac equations of constraint dynamics rectify this flaw. The forms of these equations can not only be derived from quantum field theory they can also be derived purely in the context of Dirac's constraint dynamics and relativistic mechanics and quantum mechanics. Their structures, unlike the more familiar two-body Dirac equation of Breit, which is a single equation, are that of two simultaneous quantum relativistic wave equations. A single two-body Dirac equation similar to the Breit equation can be derived from the TBDE. Unlike the Breit equation, it is manifestly covariant and free from the types of singularities that prevent a strictly nonperturbative treatment of the Breit equation. In applications of the TBDE to QED, the two particles interact by way of four-vector potentials derived from the field theoretic electromagnetic interactions between the two particles. In applications to QCD, the two particles interact by way of four-vector potentials and Lorentz invariant scalar interactions, derived in part from the field theoretic chromomagnetic interactions between the quarks and in part by phenomenological considerations. As with the Breit equation a sixteen-component spinor Ψ is used.

In the Newman–Penrose (NP) formalism of general relativity, independent components of the Ricci tensors of a four-dimensional spacetime are encoded into seven Ricci scalars which consist of three real scalars , three complex scalars and the NP curvature scalar . Physically, Ricci-NP scalars are related with the energy–momentum distribution of the spacetime due to Einstein's field equation.

In statistical mechanics, Lee–Yang theory, sometimes also known as Yang–Lee theory, is a scientific theory which seeks to describe phase transitions in large physical systems in the thermodynamic limit based on the properties of small, finite-size systems. The theory revolves around the complex zeros of partition functions of finite-size systems and how these may reveal the existence of phase transitions in the thermodynamic limit.

References

  1. Zvezdin, Anatoly Konstantinovich (1997), Taylor & Francis Group (ed.), Modern magneto-optics and magneto-optical materials : Studies in Condensed Matter, ISBN   978-0-7503-03620 .
  2. Hubert, Alex (1998), Springer (ed.), Magnetic domains, ISBN   978-3-540-85054-0
  3. Kimel (2005). "Observation of Giant Magnetic Linear Dichroism in (Ga,Mn)As". Physical Review Letters. 94 (22): 227203. Bibcode:2005PhRvL..94v7203K. doi:10.1103/physrevlett.94.227203. hdl: 2066/32798 . PMID   16090433..
  4. Shihab (2015). "Systematic study of the spin stiffness dependence on Phosphorus alloying in (Ga,Mn)As ferromagnetic semiconductor" (PDF). Applied Physics Letters. 106 (14): 142408. Bibcode:2015ApPhL.106n2408S. doi:10.1063/1.4917423.

Further reading