# Scientific Papers of Josiah Willard Gibbs, Volume 2/Chapter XII

XII.
ON DOUBLE REFRACTION IN PERFECTLY TRANSPARENT MEDIA WHICH EXHIBIT THE PHENOMENA OF CIRCULAR POLARIZATION.
[American Journal of Science, ser. 3, vol. xxiii, pp. 460–476, June, 1882.]
1. In the April number of this Journal, the velocity of propagation of a system of plane waves of light, regarded as oscillating electrical fluxes, was discussed with such a degree of approximation as would account for the dispersion of colors and give Fresnel's laws of double refraction. It is the object of this paper to supplement that discussion by carrying the approximation so much further as is necessary in order to embrace the phenomena of circularly polarizing media.
2. If we imagine all the velocities in any progressive system of plane waves to be reversed at a given instant without affecting the displacements, and the system of wave-motion thus obtained to be superposed upon the original system, we obtain a system of stationary waves having the same wave-length and period of oscillation as the original progressive system. If we then reduce the magnitude of the displacements in the uniform ratio of two to one, they will be identical, at an instant of maximum displacement, with those of the original system at the same instant.
Following the same method as in the paper cited, let us especially consider the system of stationary waves, and divide the whole displacement into the regular part, represented by $\xi ,\eta ,\zeta$, and the irregular part, represented by $\xi ',\eta ',\zeta '$, in accordance with the definitions of §2 of that paper.
3. The regular part of the displacement is subject to the equations of wave-motion, which may be written (in the most general case of plane stationary waves)
| $\xi =\left(\alpha _{1}\cos {2\pi \,{\frac {u}{l}}}+\alpha _{2}\sin {2\pi \,{\frac {u}{l}}}\right)\cos {2\pi \,{\frac {t}{p}}},$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$(1) |
| $\eta =\left(\beta _{1}\cos {2\pi \,{\frac {u}{l}}}+\beta _{2}\sin {2\pi \,{\frac {u}{l}}}\right)\cos {2\pi \,{\frac {t}{p}}},$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$(1) |
| $\zeta =\left(\gamma _{1}\cos {2\pi \,{\frac {u}{l}}}+\gamma _{2}\sin {2\pi \,{\frac {u}{l}}}\right)\cos {2\pi \,{\frac {t}{p}}},$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$(1) |
$l$
$p$
$u$
$\alpha _{1},\beta _{1},\gamma _{1}$
$\xi ,\eta ,\zeta$
$\alpha _{2},\beta _{2},\gamma _{2}$
$u$
${\text{L, M, N}}$
$u$
| ${\text{L}}^{2}+{\text{M}}^{2}+{\text{N}}^{2}=1,$ | (2) |
| $u={\text{L}}x+{\text{M}}y+{\text{N}}z,$ | (3) |
| ${\text{L}}\alpha _{1}+{\text{M}}\beta _{1}+{\text{N}}\gamma _{1}=0,$${\text{L}}\alpha _{2}+{\text{M}}\beta _{2}+{\text{N}}\gamma _{2}=0.$ | (3) |
$(\xi ',\eta ',\zeta ')$
$(\xi ,\eta ,\zeta ),$
$\xi ',\eta ',\zeta '$
$\xi ,\eta ,\zeta$
$\xi ',\eta ',\zeta ',$
$\xi ,\eta ,\zeta$
* ${\frac {d\xi }{dx}},\,\,{\frac {d\eta }{dx}},\,\,{\frac {d\zeta }{dx}},\,\,{\frac {d\xi }{dy}},\,\,{\text{etc.}}$
$\xi ,\eta ,\zeta ,$
From this it follows immediately that with the same degree of approximation ${\dot {\xi '}},{\dot {\eta '}},{\dot {\zeta }},$ may be regarded, for a given point of the medium and light of a given period, as linear functions of ${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }}$ and the differential coefficients of ${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }}$ with respect to the coordinates. For these twelve quantities we shall write ${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }}$ and diff. coeff.
5. Let us now proceed to equate the statical energy of the medium at an instant of no velocity with its kinetic energy at an instant of no displacement. It will be convenient to estimate each of these quantities for a unit of volume.
6. The statical energy of an infinitesimal element of volume may be represented by $\sigma dv,$ where $\sigma$ is a quadratic function of the components of displacement ${\xi +\xi ',\eta +\eta ',\zeta +\zeta '}.$ Since for that element of volume $\xi ',\eta ',\zeta '$ may be regarded as linear functions of $\xi ,\eta ,\zeta ,$ and diff. coeff., ​we may regard $\sigma$ as a quadratic function of $\xi ,\eta ,\zeta$ and diff, coeff., or as a linear function of the seventy-eight squares and products of these quantities. But the seventy-eight coefficients by which this function is expressed will vary with the position of the element of volume with respect to the surrounding molecules.
In estimating the statical energy for any considerable space by the integral
* $\int \sigma dv,$
$\sigma$
$s$
$\xi ,\eta ,\zeta$
$\sigma ,$
* $\int sdv.$
$s$
We may divide $s$ into three parts, of which the first ($s_{'}$) contains the squares and products of $\xi ,\eta ,\zeta ,$ the second ($s_{''}$) contains the products of $\xi ,\eta ,\zeta$ with the differential coefficients, and the third ($s_{'''}$) contains the squares and products of the differential coefficients. It is evident that the average statical energy of the whole medium per unit of volume is the space-average of $s,$ and that it will consist of three parts, which are the space-averages of $s_{'},s_{''},$ and $s_{'''},$ respectively. These parts we may call ${\text{S}}_{'},{\text{S}}_{''},$ and ${\text{S}}_{'''}.$ Only the first of these was considered in the preceding paper.
Now the considerations which justify us in neglecting, for an approximate estimate, the terms of $s$ which contain the differential coefficients of $\xi ,\eta ,\zeta$ with respect to the coordinates, will apply with especial force to the terms which contain the squares and products of these differential coefficients. Therefore, to carry the approximation one step beyond that of the preceding paper, it will only be necessary to take account of $s_{'}$ and $s_{''}$ and of ${\text{S}}_{'}$ and ${\text{S}}_{''}.$
7. We may set
| $s_{'}={\text{A}}\xi ^{2}+{\text{B}}\eta ^{2}+{\text{C}}\zeta ^{2}+{\text{E}}\eta \zeta +{\text{F}}\zeta \xi +{\text{G}}\xi \eta ,$ | (5) |
${\text{A, B, C, D, E, F, G}}$
Since the average values of
* $\sin ^{2}2\pi {\frac {u}{l}},\,\cos ^{2}2\pi {\frac {u}{l}},\,\,\sin 2\pi {\frac {u}{l}}\cos 2\pi {\frac {u}{l}}$
${\tfrac {1}{2}},{\tfrac {1}{2}}$
* $\cos ^{2}2\pi {\frac {t}{p}}=1,$
| ${\begin{aligned}{\text{S}}_{'}&={\tfrac {1}{2}}({\text{A}}\alpha _{1}^{2}+{\text{B}}\beta _{1}^{2}+{\text{C}}\gamma _{1}^{2}+{\text{E}}\beta _{1}\gamma _{1}+{\text{F}}\gamma _{1}\alpha _{1}+{\text{G}}\alpha _{1}\beta _{1})\\&+{\tfrac {1}{2}}({\text{A}}\alpha _{2}^{2}+{\text{B}}\beta _{2}^{2}+{\text{C}}\gamma _{2}^{2}+{\text{E}}\beta _{2}\gamma _{2}+{\text{F}}\gamma _{2}\alpha _{2}+{\text{G}}\alpha _{2}\beta _{2}).\end{aligned}}$ | (6) |
8. The second part of the statical energy of the whole medium per unit of volume (${\text{S}}_{''}$) is the space-average of $s_{''},$ which is a linear function of the twenty-seven products of $\xi ,\eta ,\zeta$ with their differential coefficients with respect to the coordinates. Now since
* $\xi {\frac {d\xi }{dx}}={\tfrac {1}{2}}{\frac {d(\xi ^{2})}{dx}},$$\eta {\frac {d\eta }{dx}}={\tfrac {1}{2}}{\frac {d(\eta ^{2})}{dx}},$etc.,
${\text{S}}_{''}.$
* $\eta {\frac {d\zeta }{dx}}$and$\xi {\frac {d\eta }{dx}}\cdot$
* $\eta {\frac {d\zeta }{dx}}+\zeta {\frac {d\eta }{dx}}={\frac {d(\eta \zeta )}{dx}},$
$s_{''}$
${\text{S}}_{''}.$
* $\eta {\frac {d\zeta }{dx}}-\zeta {\frac {d\eta }{dx}},$$\zeta {\frac {d\xi }{dx}}-\xi {\frac {d\zeta }{dx}},$$\xi {\frac {d\eta }{dx}}-\eta {\frac {d\xi }{dx}},$etc.
$dy$
$dz$
$dx$
$s_{''}$
${\text{S}}_{''}$
* $\eta {\frac {d\zeta }{dx}}-\zeta {\frac {d\eta }{dx}}={\text{L}}\left(\eta {\frac {d\zeta }{du}}-\zeta {\frac {d\eta }{du}}\right),$
* ${\frac {2\pi {\text{L}}}{l}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2}).$
$l{\text{S}}_{''}$
${\text{L, M, N}}$
* $\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2},$$\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2},$$\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2}$
| $\Theta ={\text{L}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2})+{\text{M}}(\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2})+{\text{N}}(\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2}),$ | (7) |
| ${\text{L}}\Theta =\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2},$${\text{M}}\Theta =\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2},$${\text{N}}\Theta =\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2}.$ | (8) |
$l{\text{S}}_{''}$
${\text{L, M, N}}$
${\text{L}}\Theta ,{\text{M}}\Theta ,{\text{N}}\Theta .$
$l{\text{S}}_{''}$
${\text{L, M }}$
${\text{N.}}$
| ${\text{S}}_{''}={\frac {\Phi }{l}}\Theta ={\frac {\Phi }{l}}[{\text{L}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2})+{\text{M}}(\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2})+{\text{N}}(\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2})],$ | (9) |
$\Phi$
${\text{L, M }}$
${\text{N,}}$
9. It will be useful to consider more closely the geometrical significance of the quantity $\Theta .$ For this purpose it will be convenient to have a definite understanding with respect to the relative position of the coordinate axes.
We shall suppose that the axes of X, Y, and Z are related in the same way as lines drawn to the right, forward and upward, so that a rotation from X to Y appears clockwise to one looking in the direction of Z.
Now if from any same point, as the origin of coordinates, we lay off lines representing in direction and magnitude the displacements in all the different wave-planes, we obtain an ellipse, which we may call the displacement-ellipse. Of this, one radius vector ($\rho _{1}$) will have the components $\alpha _{1},\beta _{1},\gamma _{1}$ and another ($\rho _{2}$) the components $\alpha _{2},\beta _{2},\gamma _{2}.$ These will belong to conjugate diameters, each being parallel to the tangent at the extremity of the other. The area of the ellipse will therefore be equal to the parallelogram of which $\rho _{1}$ and $\rho _{2}$ are two sides, multiplied by $\pi .$ Now it is evident that $\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2},\,\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2},\,\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2}$ are numerically equal to the projections of this parallelogram on the planes of the coordinate axes, and are each positive or negative according as a revolution from $\rho _{1}$ to $\rho _{2}$ appears clockwise or counter-clockwise to one looking in the direction of the proper coordinate axis. Hence, $\Theta$ will be numerically equal to the parallelogram, that is, to the area of the displacement-ellipse divided by $\pi ,$ and will be positive or negative ​according as a revolution from $\rho _{1}$ to $\rho _{2}$ appears clockwise or counter-clockwise to one looking in the direction of the wave-normal. Since $\rho _{1}$ and $\rho _{2}$ are determined by displacements in planes one-quarter of a wave-length distant from each other, and the plane to which the latter relates lies on the side toward which the wave-normal is drawn, it follows that $\Theta$ is positive or negative according as the combination of displacements has the character of a right-handed or a left-handed screw.
10. The kinetic energy of the medium, which is to be estimated for an instelnt of no displacement, may be shown as in § 7 of the former paper (page 185 of this volume) to consist of two parts, of which one relates to the regular flux (${\dot {\xi }},{\dot {\eta }},{\dot {zeta}}$), and the other to the irregular flux (${\dot {\xi '}},{\dot {\eta '}},{\dot {zeta'}}$). The first, in the notation of that paper, is represented by
* ${\tfrac {1}{2}}\int ({\dot {\xi }}{\text{ Pot }}{\dot {\xi }}+{\dot {\eta }}{\text{ Pot }}{\dot {\eta }}+{\dot {\zeta }}{\text{ Pot }}{\dot {\zeta }})dv,$
* ${\frac {l^{2}}{2\pi }}\int ({\dot {\xi }}^{2}+{\dot {\eta }}^{2}+{\dot {\zeta }}^{2})dv.$
| ${\text{T}}={\frac {\pi l^{2}}{p^{2}}}(\alpha _{1}^{2}+\beta _{1}^{2}+\gamma _{1}^{2}+\alpha _{2}^{2}+\beta _{2}^{2}+\gamma _{2}^{2}).$ | (10) |
* ${\tfrac {1}{2}}\int ({\dot {\xi '}}{\text{ Pot }}{\dot {\xi '}}+{\dot {\eta '}}{\text{ Pot }}{\dot {\eta '}}+{\dot {\zeta '}}{\text{ Pot }}{\dot {\zeta '}})dv.$
${\dot {\xi '}},{\dot {\eta '}},{\dot {\zeta '}}$
${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }}$
$Pot$
${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }}$
${\text{Pot }}{\dot {\xi '}},{\text{Pot }}{\dot {\eta '}},{\text{Pot }}{\dot {\zeta '}}$
${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }}$
* ${\tfrac {1}{2}}({\dot {\xi '}}{\text{ Pot }}{\dot {\xi '}}+{\dot {\eta '}}{\text{ Pot }}{\dot {\eta '}}+{\dot {\zeta '}}{\text{ Pot }}{\dot {\zeta '}})$
${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }}$
​Yet, as in the case of the statical energy, we may substitute the average values of these coefficients for the coefficients themselves in the integral by which we obtain the energy of any considerable space. The kinetic energy due to the irregular part of the flux is thus reduced to a quadratic function of ${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }}$ and diff. coeff. which has constant coefficients for a given medium and light of a given period.
The function may be divided into three parts, of which the first contains the squares and products of ${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }},$ the second the products of ${\dot {\xi }},{\dot {\eta }},{\dot {\zeta }}$ with their differential coefficients, and the third, which may be neglected, the squares and products of the differential coefficients.
We may proceed with the reduction precisely as in the case of the statical energy, except that the differentiations with respect to the time will introduce the constant factor ${\frac {4\pi ^{2}}{p^{2}}}\cdot$ This will give for the first part of the kinetic energy of the irregular flux per unit of volume
| ${\text{T}}'={\frac {2\pi ^{2}}{p^{2}}}({\text{A}}'\alpha _{1}^{2}+{\text{B}}'\beta _{1}^{2}+{\text{C}}'\gamma _{1}^{2}+{\text{E}}'\beta _{1}\gamma _{1}+{\text{F}}'\gamma _{1}\alpha _{1}+{\text{G}}'\alpha _{1}\beta _{1})$ $+{\frac {2\pi ^{2}}{p^{2}}}({\text{A}}'\alpha _{2}^{2}+{\text{B}}'\beta _{2}^{2}+{\text{C}}'\gamma _{2}^{2}+{\text{E}}'\beta _{2}\gamma _{2}+{\text{F}}'\gamma _{2}\alpha _{2}+{\text{G}}'\alpha _{2}\beta _{2}),$ | (11) |
| ${\begin{aligned}{\text{T}}_{''}'&={\frac {4\pi ^{2}\Phi '}{p^{2}l}}\Theta \\&={\frac {4\pi ^{2}\Phi '}{p^{2}l}}[{\text{L}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2})+{\text{M}}(\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2})+{\text{N}}(\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2})],\end{aligned}}$ | (12) |
${\text{A}}',{\text{B}}',{\text{C}}',{\text{E}}',{\text{F}}',{\text{G}}'$
$\Phi '$
${\text{L, M,}}$
${\text{N}},$
12. Equating the statical and kinetic energies, we have
* ${\text{S}}_{'}+{\text{S}}_{''}={\text{T}}+{\text{T}}_{'}'+{\text{T}}_{''}',$
| ${\begin{aligned}&{\tfrac {1}{2}}({\text{A}}\alpha _{1}^{2}+{\text{B}}\beta _{1}^{2}+{\text{C}}\gamma _{1}^{2}+{\text{E}}\beta _{1}\gamma _{1}+{\text{F}}\gamma _{1}\alpha _{1}+{\text{G}}\alpha _{1}\beta _{1})\\+&{\tfrac {1}{2}}({\text{A}}\alpha _{2}^{2}+{\text{B}}\beta _{2}^{2}+{\text{C}}\gamma _{2}^{2}+{\text{E}}\beta _{2}\gamma _{2}+{\text{F}}\gamma _{2}\alpha _{2}+{\text{G}}\alpha _{2}\beta _{2})\\=&{\frac {\Phi }{l}}[{\text{L}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2})+{\text{M}}(\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2})+{\text{N}}(\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2})]\\=&{\frac {\pi l^{2}}{p^{2}}}(\alpha _{1}^{2}+\beta _{1}^{2}+\gamma _{1}^{2}+\alpha _{2}^{2}+\beta _{2}^{2}+\gamma _{2}^{2})\\+&{\frac {2\pi ^{2}}{p^{2}}}({\text{A}}'\alpha _{1}^{2}+{\text{B}}'\beta _{1}^{2}+{\text{C}}'\gamma _{1}^{2}+{\text{E}}'\beta _{1}\gamma _{1}+{\text{F}}'\gamma _{1}\alpha _{1}+{\text{G}}'\alpha _{1}\beta _{1})\\+&{\frac {2\pi ^{2}}{p^{2}}}({\text{A}}'\alpha _{2}^{2}+{\text{B}}'\beta _{2}^{2}+{\text{C}}'\gamma _{2}^{2}+{\text{E}}'\beta _{2}\gamma _{2}+{\text{F}}'\gamma _{2}\alpha _{2}+{\text{G}}'\alpha _{2}\beta _{2})\\+&{\frac {4\pi ^{2}\Phi '}{p^{2}l}}[{\text{L}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2})+{\text{M}}(\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2})+{\text{N}}(\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2})]\end{aligned}}$ | (13) |
| $a={\frac {\text{A}}{2\pi }}-{\frac {2\pi {\text{A}}'}{p^{2}}},$$b={\frac {\text{B}}{2\pi }}-{\frac {2\pi {\text{B}}'}{p^{2}}},$etc., | (14) |
| $\phi ={\frac {\Phi }{2\pi p}}-{\frac {2\pi \phi '}{p^{3}}},$ | (15) |
| ${\begin{aligned}&a\alpha _{1}^{2}+b\beta _{1}^{2}+c\gamma _{1}^{2}+e\beta _{1}\gamma _{1}+f\gamma _{1}\alpha _{1}+g\alpha _{1}\beta _{1}\\+&a\alpha _{2}^{2}+b\beta _{2}^{2}+c\gamma _{2}^{2}+e\beta _{2}\gamma _{2}+f\gamma _{2}\alpha _{2}+g\alpha _{2}\beta _{2}\\+&{\frac {2p\phi }{l}}[{\text{L}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2})+{\text{M}}(\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2})+{\text{N}}(\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2})]\\=&{\frac {l^{2}}{p^{2}}}(\alpha _{1}^{2}+\beta _{1}^{2}+\gamma _{1}^{2}+\alpha _{2}^{2}+\beta _{2}^{2}+\gamma _{2}^{2}),\end{aligned}}$ | (16) |
${a,b,c,e,f,g}$
$\phi$
${\text{L, M, N,}}$
13. Now this equation, which expresses a relation between the constants of the equations of wave-motion (1), will apply, with those equations, not only to such vibrations as actually take place, but also to such as we may imagine to take place under the influence of constraints determining the type of vibration. The free or unconstrained vibrations, with which alone we are concerned, are characterized by this, that infinitesimal variations (by constraint) of the type of vibration, that is, of the ratios of the quantities ${\alpha _{1},\beta _{1},\gamma _{1},\alpha _{2},\beta _{2},\gamma _{2},}$ will not affect the period by any quantity of the same order of magnitude. These variations must however be consistent with equations (4), which require that
| ${\text{L }}d\alpha _{1}+{\text{M }}d\beta _{1}+{\text{N }}d\gamma _{1}=0,$${\text{L }}d\alpha _{2}+{\text{M }}d\beta _{2}+{\text{N }}d\gamma _{2}=0.$ | (17) |
${\alpha _{1},\beta _{1},\gamma _{1},\alpha _{2},\beta _{2},\gamma _{2},}$
${d\alpha _{1},d\beta _{1},d\gamma _{1},d\alpha _{2},d\beta _{2},d\gamma _{2},}$
${d\alpha _{1},d\beta _{1},d\gamma _{1},}$
${d\alpha _{2},d\beta _{2},d\gamma _{2},}$
${\alpha _{1},\beta _{1},\gamma _{1},}$
${\alpha _{2},\beta _{2},\gamma _{2},}$
${\alpha _{1},\beta _{1},\gamma _{1},}$
${\alpha _{1},\beta _{1},\gamma _{1},}$
${\alpha _{2},\beta _{2},\gamma _{2},}$
${d\alpha _{1},d\beta _{1},d\gamma _{1},}$
${\alpha _{2},\beta _{2},\gamma _{2},}$
If we differentiate with respect to ${\alpha _{1},\beta _{1},\gamma _{1},}$ and write ${\alpha _{1},\beta _{1},\gamma _{1},}$ for ${d\alpha _{1},d\beta _{1},d\gamma _{1},}$ we obtain
| ${\begin{aligned}&a\alpha _{1}^{2}+b\beta _{1}^{2}+c\gamma _{1}^{2}+e\beta _{1}\gamma _{1}+f\gamma _{1}\alpha _{1}+g\alpha _{1}\beta _{1}\\+&{\frac {p\phi }{l}}[{\text{L}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2})+{\text{M}}(\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2})+{\text{N}}(\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2})]\\=&{\frac {l^{2}}{p^{2}}}(\alpha _{1}^{2}+\beta _{1}^{2}+\gamma _{1}^{2}).\end{aligned}}$ | (18) |
$\alpha _{1},\beta _{1},\gamma _{1},$
$\alpha _{2},\beta _{2},\gamma _{2}$
$d\alpha _{1},d\beta _{1},d\gamma _{1},$
| ${\begin{aligned}2a\alpha _{1}\alpha _{2}+2b\beta _{1}\beta _{2}&+2c\gamma _{1}\gamma _{2}+e(\beta _{1}\gamma _{2}+\gamma _{1}\beta _{2})+f(\gamma _{1}\alpha _{2}+\alpha _{1}\gamma _{2})\\&+g(\alpha _{1}\beta _{2}+\beta _{1}\alpha _{2})={\frac {2l^{2}}{p^{2}}}(\alpha _{1}\alpha _{2}+\beta _{1}\beta _{2}+\gamma _{1}\gamma _{2}).\end{aligned}}$ | (19) |
$\alpha _{2},\beta _{2},\gamma _{2},$
$\alpha _{2},\beta _{2},\gamma _{2}$
$d\alpha _{2},d\beta _{2},d\gamma _{2},$
| ${\begin{aligned}&a\alpha _{2}^{2}+b\beta _{2}^{2}+c\gamma _{2}^{2}+e\beta _{2}\gamma _{2}+f\gamma _{2}\alpha _{2}+g\alpha _{2}\beta _{2}\\+&{\frac {p\phi }{l}}[{\text{L}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2})+{\text{M}}(\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2})+{\text{N}}(\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2})]\\=&{\frac {l^{2}}{p^{2}}}(\alpha _{2}^{2}+\beta _{2}^{2}+\gamma _{2}^{2}).\end{aligned}}$ | (20) |
$\alpha _{2},\beta _{2},\gamma _{2}$
$\alpha _{1},\beta _{1},\gamma _{1},$
$d\alpha _{2},d\beta _{2},d\gamma _{2},$
14. The geometrical signification of our equations may now be simplified by a suitable choice of the position of the origin of coordinates, which is as yet wholly arbitrary.
We shall hereafter suppose that the origin is placed in a plane of maximum or minimum displacement, if such there are. In the case of circular polarization, in which the displacements are everywhere equal, its position is immaterial. The lines $\rho _{1}$ and $\rho _{2},$ of which $\alpha _{1},\beta _{1},\gamma _{1}$ and $\alpha _{2},\beta _{2},\gamma _{2}$ are respectively the components, will now be the semi-axes of the displacement-ellipse, and therefore at right angles. (See § 9.) The case of circular polarization will not constitute any exception. Hence,
| $\alpha _{1}\alpha _{2}+\beta _{1}\beta _{2}+\gamma _{1}\gamma _{2}=0,$ | (21) |
| $\Theta ={\text{L}}(\beta _{1}\gamma _{2}-\gamma _{1}\beta _{2})+{\text{M}}(\gamma _{1}\alpha _{2}-\alpha _{1}\gamma _{2})+{\text{N}}(\alpha _{1}\beta _{2}-\beta _{1}\alpha _{2})=\pm \rho _{1}\rho _{2},$ | (22) |
15. Equation (19) is now reduced to the form
| ${\begin{aligned}2a\alpha _{1}\alpha _{2}+2b\beta _{1}\beta _{2}&+e(\beta _{1}\gamma _{2}+\gamma _{1}\beta _{2})\\&+f(\gamma _{1}\alpha _{2}+\alpha _{1}\gamma _{2})+g(\alpha _{1}\beta _{2}+\beta _{1}\alpha _{2})=0,\end{aligned}}$ | (23) |
| $ax^{2}+by^{2}+cz^{2}+eyz+fzx+gxy,$ | (24) |
* ${\begin{aligned}2ax_{1}x_{2}+2by_{1}y_{2}+2cz_{1}z_{2}&+e(y_{1}z_{2}+z_{1}y_{2})\\&+f(z_{1}x_{2}+x_{1}z_{2})+g(x_{1}y_{2}+y_{1}x_{2})=0\end{aligned}}$
$x_{1},y_{1},z_{1}$
$x_{2},y_{2},z_{2}$
$\alpha _{1},\beta _{1},\gamma _{1}$
$\alpha _{2},\beta _{2},\gamma _{2}$
$\alpha _{1},\beta _{1},\gamma _{1}$
$\alpha _{2},\beta _{2},\gamma _{2}$
16. If we write ${\text{U}}_{1},{\text{U}}_{2}$ for the reciprocals of the semi-axes of the central section of the ellipsoid (24) by a wave-plane, ${\text{U}}_{1}$ being the reciprocal of the one to which the displacement $\alpha _{1},\beta _{1},\gamma _{1}$ is parallel, we have
| $a\alpha _{1}^{2}+b\beta _{1}^{2}+c\gamma _{1}^{2}+e\beta _{1}\gamma _{1}+f\gamma _{1}\alpha _{1}+g\alpha _{1}\beta _{1}={\text{U}}_{1}^{2}(\alpha _{1}^{2}+\beta _{1}^{2}+\gamma _{1}^{2}),$ | (25) |
$\alpha _{1},\beta _{1},\gamma _{1}.$
| $a\alpha _{2}^{2}+b\beta _{2}^{2}+c\gamma _{2}^{2}+e\beta _{2}\gamma _{2}+f\gamma _{2}\alpha _{2}+g\alpha _{2}\beta _{2}={\text{U}}_{2}^{2}(\alpha _{2}^{2}+\beta _{2}^{2}+\gamma _{2}^{2}).$ | (26) |
${\text{V}}$
| ${\text{V}}={\frac {l}{p}}\cdot$ | (27) |
| ${\text{U}}_{1}^{2}\rho _{1}^{2}\pm {\frac {\phi }{\text{V}}}\rho _{1}\rho _{2}={\text{V}}^{2}\rho _{1}^{2},$${\text{U}}_{2}^{2}\rho _{2}^{2}\pm {\frac {\phi }{\text{V}}}\rho _{1}\rho _{2}={\text{V}}^{2}\rho _{2}^{2},$ | (28) |
We may here observe that in case $\phi =0$ the solution of these ​equations is very simple. We have necessarily either $\rho _{2}=0$ and ${\text{V}}^{2}={\text{U}}_{1}^{2},$ or $\rho _{1}=0$ and ${\text{V}}^{2}={\text{U}}_{1}^{2}.$ In this case, the light is linearly polarized, and the directions of oscillation and the velocities of propagation are given by FresneFs law. Experiment has shown that this is the usual case. We wish, however, to investigate the case in which $\phi$ does not vanish. Since the term containing $\phi$ arises from the consideration of those quantities which it was allowable to neglect in the first approximation, we may assume that $\phi$ is always very small in comparison with ${\text{V}}^{3},{\text{U}}_{1}^{3},$ or ${\text{U}}_{2}^{3}.$
17. Equations (28) may be written
| ${\text{V}}^{2}-{\text{U}}_{1}^{2}=\pm {\frac {\phi }{\text{V}}}{\frac {\rho _{2}}{\rho _{1}}},$${\text{V}}^{2}-{\text{U}}_{2}^{2}=\pm {\frac {\phi }{\text{V}}}{\frac {\rho _{1}}{\rho _{2}}}\cdot$ | (29) |
| ${\text{V}}^{2}({\text{V}}^{2}-{\text{U}}_{1}^{2})({\text{V}}^{2}-{\text{U}}_{2}^{2})=\phi ^{2}.$ | (30) |
$\phi$
${\text{V}}^{2},$
${\text{U}}_{1}^{2}$
${\text{U}}_{2}^{2}$
${\text{V}}^{2}$
${\text{U}}_{1}^{2}$
${\text{U}}_{2}^{2}.$
For the numerical computation of ${\text{V}},$ when ${\text{U}}_{1},{\text{U}}_{2},$ and $\phi$ are known numerically, we may divide the equation by ${\text{V}}^{2},$ and then solve it as if the second member were known. This will give
| ${\text{V}}^{2}={\frac {{\text{U}}_{1}^{2}+{\text{U}}_{2}^{2}}{2}}\pm {\sqrt {{\frac {\phi ^{2}}{{\text{V}}^{2}}}+{\frac {({\text{U}}_{1}^{2}-{\text{U}}_{2}^{2})}{4}}}}\cdot$ | (31) |
${\text{U}}_{1}{\text{U}}_{2}$
${\text{V}}^{2}$
${\text{V}}^{2}.$
${\text{V}}^{2}$
For either value of ${\text{V}}^{2}$ we may easily find the ratio of $\rho _{1}$ to $\rho _{2},$ that is, the ratio of the axes of the displacement-ellipse, from one of equations (29), or from the equation
| ${\frac {\rho _{2}^{2}}{\rho _{1}^{2}}}={\frac {{\text{V}}^{2}-{\text{U}}_{1}^{2}}{{\text{V}}^{2}-{\text{U}}_{2}^{2}}}$ | (32) |
​In equations (29), we are to read $+$ or $-$ in the second members, according as the ray is right-handed or left-handed. (See § 16.) It follows that if the value of $\phi$ is positive, the greater velocity will belong to a right-handed ray, and the smaller to a left-handed, but if the value of $\phi$ is negative, the opposite is the case. Except when $\phi =0,$ and the polarization is linear, there will be one right-handed and one left-handed ray for any given wave-normal and period.
18. When ${\text{U}}_{1}={\text{U}}_{2},$ equations (29) give
* $\rho _{1}=\rho _{2},$${\text{V}}^{2}={\text{U}}^{2}\pm {\frac {\phi }{\text{V}}},$
${\text{U}}$
${\text{U}}_{1}$
${\text{U}}_{2}.$
${\text{U}}_{1}={\text{U}}_{2}.$
${\text{V}}_{\text{R}}$
${\text{V}}_{\text{L}},$
| ${\text{V}}_{\text{R}}^{2}={\text{U}}^{2}+{\frac {\phi }{{\text{V}}_{\text{R}}}},$${\text{V}}_{\text{L}}^{2}={\text{U}}^{2}-{\frac {\phi }{{\text{V}}_{\text{L}}}};$ | (33) |
* ${\text{V}}_{\text{R}}^{2}-{\text{V}}_{\text{L}}^{2}=\phi \left({\frac {1}{{\text{V}}_{\text{R}}}}+{\frac {1}{{\text{V}}_{\text{L}}}}\right)=\phi {\frac {{\text{V}}_{\text{R}}+{\text{V}}_{\text{L}}}{{\text{V}}_{\text{R}}{\text{V}}_{\text{L}}}},$
| ${\text{V}}_{\text{R}}-{\text{V}}_{\text{L}}={\frac {\phi }{{\text{V}}_{\text{R}}{\text{V}}_{\text{L}}}}\cdot$ | (34) |
$\theta$
| $\theta ={\frac {\pi }{p}}\left({\frac {1}{{\text{V}}_{\text{L}}}}-{\frac {1}{{\text{V}}_{\text{R}}}}\right)\cdot$ | (35) |
| $\theta ={\frac {\pi \phi }{p{\text{V}}_{\text{R}}^{2}{\text{V}}_{\text{L}}^{2}}}\cdot$ | (36) |
${\text{U}}^{4}$
${\text{V}}_{\text{R}}^{2}{\text{V}}_{\text{L}}^{2},$
| $\theta ={\frac {\pi \phi }{p{\text{U}}^{4}}}\cdot$ | (37) |
$\xi ,\eta ,\zeta ,$
$\Phi$
$\Phi '$
| $\theta ={\frac {\Phi }{2p^{2}{\text{V}}_{\text{R}}^{2}{\text{V}}_{\text{L}}^{2}}}-{\frac {2\pi ^{2}\Phi '}{p^{4}{\text{V}}_{\text{R}}^{2}{\text{V}}_{\text{L}}^{2}}}\cdot$ | (38) |
* $p={\frac {\lambda }{k}},$${\text{V}}_{\text{R}}={\frac {k}{n_{\text{R}}}},$${\text{V}}_{\text{L}}={\frac {k}{n_{\text{L}}}},$${\text{U}}={\frac {k}{n}},$
$k$
$\lambda$
$n_{\text{R}},n_{\text{L}}$
$n$
| $\theta ={\frac {\Phi n_{\text{R}}^{2}n_{\text{L}}^{2}}{2k^{2}\lambda ^{2}}}-{\frac {2\pi ^{2}\Phi 'n_{\text{R}}^{2}n_{\text{L}}^{2}}{\lambda ^{4}}}\cdot$ | (39) |
| $\theta =n_{\text{R}}^{2}n_{\text{L}}^{2}\left({\frac {\text{K}}{\lambda ^{2}}}+{\frac {{\text{K}}'}{\lambda ^{4}}}\right),$ | (40) |
${\text{K}}$
${\text{K}}'$
$n^{4}$
$n_{\text{R}}^{2}n_{\text{L}}^{2}.$
$\lambda ,$
${\text{K}}$
${\text{K}}'.$
20. But without any such assumption as that contained in the last paragraph, we may easily obtain formulæ for the experimental determination of $\Phi$ and $\Phi '$ for the optic axis of a uniaxial crystal. Considerations analogous to those of § 13 of the former paper (page 190 of this volume), show that in differentiating equation (39) we may regard $\Phi$ and $\Phi '$ as constant, although they may actually vary with $\lambda .$ This equation may be written
| ${\frac {\theta \lambda ^{2}}{n^{4}}}={\frac {\Phi }{2k^{2}}}-{\frac {2\pi ^{2}\Phi '}{\lambda ^{2}}}\cdot$ | (41) |
| ${\frac {d\left({\frac {\theta \lambda ^{2}}{n^{4}}}\right)}{d\left({\frac {1}{\lambda ^{2}}}\right)}}=-2\pi ^{2}\theta '.$ | (42) |
$\Phi '$
$\Phi$
21. If we wish to represent $\phi$ geometrically, like ${\text{U}}_{1}$ and ${\text{U}}_{2},$ we may construct the surfaces
| $\scriptstyle {\text{A}}\displaystyle x^{2}+\scriptstyle {\text{B}}\displaystyle y^{2}+\scriptstyle {\text{C}}\displaystyle z^{2}+\scriptstyle {\text{E}}\displaystyle yz+\scriptstyle {\text{F}}\displaystyle zx++\scriptstyle {\text{G}}\displaystyle xy=\pm 1,$ | (43) |
$\scriptstyle {\text{A, B,}}$
$\phi$
${\text{L}}^{2},{\text{M}}^{2},$
$\phi ,$
$\phi$
$\phi$
$\phi$
​22. The manner in which the ellipsoid (24) may be partially determined by the relations of symmetry which the medium may possess, has been sufficiently discussed in the former paper.
With respect to the quantity $\phi ,$ and the surfaces which determine it, the following principle is of fundamental importance. If one body is identical in its internal structure with the image by reflection of another, the values of $\phi$ in corresponding lines in the two bodies will be numerically equal but have opposite signs.
It follows that if a body is identical in internal structiure with its own image by reflection, the value of $\phi$ (if not zero for all directions) must be positive for some directions and negative for others. Moreover, the above described surface by which $\phi$ is represented must consist of two conjugate hyperboloids, of which one is identical in form with the image by reflection of the other. This requires that the hyperboloids shall be light cylinders with conjugate rectangular hyperbolas for bases. A crystal characterized by such properties will belong to the tetragonal system. Since $\phi =0$ for the optic axis, it would be difficult to distinguish a case of this kind from an ordinary uniaxial crystal, unless the ellipsoid (24) should approach very closely to a sphere.
It is only in the very limited case described in the last paragraph that a medium which is identical in its internal structure with its image by reflection can have the property of circular or elliptic polarization. To media which are unlike their images by reflection, and have the property of circular polarization, we may apply the following general principles.
If the medium has any axis of symmetry, the ellipsoid or hyperboloids which represent the values of $\phi$ will have an axis in the same direction. If the medium after a revolution of less than 180° about any axis is equivalent to the medium in its first position, the ellipsoid or hyperboloids will have an axis of revolution in that direction.
23. The laws of the propagation of light in plane waves, which ​have thus been derived from the single hypothesis that the disturbance by which light is transmitted consists of solenoidal electrical fluxes, and which apply to light of different colors and to the most general case of perfectly transparent and sensibly homogeneous media not subject to magnetic action, are essentially those which are generally received as embodying the results of experiment. In no particular, so far as the writer is aware, do they conflict with the results of experiment, or require the aid of auxiliary and forced hypotheses to bring them into harmony therewith.
In this respect, the electromagnetic theory of light stands in marked contrast with that theory in which the properties of an elastic solid are attributed to the ether,—a contrast which was very distinct in Maxwell's derivation of Fresnel's laws from electrical principles, but becomes more striking as we follow the subject farther into its details, and take account of the want of absolute homogeneity in the medium, so as to embrace the phenomena of the dispersion of colors and circular and elliptical polarization.
* ↑ See page 182 of this volume.
* ↑ This ellipse, which represents the simultaneous displacements in different parts of the field, will also represent the suocessive displacements at any same point in the corresponding system of progressive waves.
* ↑ See § 9 of the former paper, on page 187 of this volume.
* ↑ Compare § 11 of the former paper, page 189 of this volume.
* ↑ The reader will perceive that an earlier limitation of the position of the origin by a supposition of this nature, involving a limitation of the values of $\alpha _{1},\beta _{1},\gamma _{1},\alpha _{2},\beta _{2},\gamma _{2},$ would have been embarrassing in the operations of the last paragraph.
* ↑ We should not attribute any physical significance to the third value of ${\text{V}}.$ For this value would imply a wave-length very small in comparison with the length of ordinary waves of light, and with respeet to which our fundamental assumption that the wave-length is very great in comparison with the distances of contiguous molecules would be entirely false. Our analysis, therefore, furnishes no reason for supposing that any such velocities are possible for the propagation of electrical disturbances.
* ↑ Our experimental knowledge of circularly or elliptically polarizing media is confined to such as are optically either isotropic or uniaxial. The general theory of such media, embracing the case of two optic axes, has however been discussed by Professor von Lang ("Theorie der Circularpolarization," Sitz. Ber. Wiener Akad., vol. lxxv, p. 719). The general results of the present paper, although derived from physical hypotheees of an entirely difierent nature, are quite similar to those of the memoir cited. They would become identical, the writer believes, by the substitution of a constant for ${\frac {p\phi }{l}}$ or ${\frac {\phi }{\text{V}}}$ in the equations of this paper. (See especiaUy equations (18), (20), (28).) That a complete discussion of the subject on any theory must include the case of biaxial media having the property of circular or elliptical polarization, is evident from the consideration that it must at least be possible to produce examples of such media artificially. An isotropic or uniaxial ciystal may be made biaxial by pressure. If it has the property of circular and elliptic polarization, that property cannot be wholly destroyed by the application of small pressures.
* ↑ When the rotation of the plane of polarization appears clockwise to the observer, it has the character of a left-handed screw. But the circularly polarized ray to which ${\text{V}}_{\text{R}}$ relates, the rotation of which also appears clockwise to the observer, has the character of a right-handed screw.
* ↑ The degree of accuracy of this substitution may be shown as follows. By (33) ${\text{V}}_{\text{R}}({\text{V}}_{\text{R}}^{2}-{\text{U}}^{2})={\text{V}}_{\text{L}}({\text{U}}^{2}-{\text{V}}_{\text{L}}^{2}),$ whence ${\text{V}}_{\text{R}}^{3}+{\text{V}}_{\text{L}}^{3}=({\text{V}}_{\text{R}}+{\text{V}}_{\text{L}}){\text{U}}^{3},$ ${\text{V}}_{\text{R}}^{2}-{\text{V}}_{\text{R}}{\text{V}}_{\text{L}}+{\text{V}}_{\text{L}}^{2}={\text{U}}^{2},$ ${\text{V}}_{\text{R}}{\text{V}}_{\text{L}}={\text{U}}^{2}({\text{V}}_{\text{R}}-{\text{V}}_{\text{L}})^{2}.$
* ↑ Compare § 12 of the former paper, on page 189 of this volume.
* ↑ The necessity of the opposite signs will perhaps appear most readily from the consideration that the direction of rotation of the plane of polarization must be opposite in the two bodies.
* ↑ There is no difficulty in conceiving of the constitution of a body which would have the properties described above. Thus, we may imagine a body with molecules of a spiral form, of which one-half are right-handed and one-half left-handed, and we may suppose that the motion of electricity is opposed by a less resistance within them than without. If the axes of the right-handed molecules are parallel to the axis of X, and those of the left-handed molecules to the axis of Y, their effects would counterbalance one another when the wave-normal is parallel to the axis of Z. But when the wave-normal (of a beam of linearly polarized light) is parallel to the axis of X, the left-handed molecules would produce a left-handed (negative) rotation of the plane of polarization, the right-handed molecules having no effect; and when the wave-normal is parallel to the axis of Y, the reverse would be the case.
* ↑ The rotation of the plane of polarization which is produced by magnetic action has been discussed by Maxwell (Treatise on Electricity and Magnetism, vol. ii, chap, xxi), and by Rowland (Amer. Journ. Math., voL iii, p. 107).