Просмотр кода страницы Двойное лучепреломление

[[Файл:Halite and Calcite in Museum.jpg|thumb|Оптические свойства [[галит]]а и [[кальцит]]а]]
[[Файл:Calcite.jpg|thumb|Двойное лучепреломление кристаллом [[кальцит]]а, положенным на бумагу с текстом]]
'''Двойно́е лучепреломле́ние''' или '''двулучепреломле́ние''' — [[Оптическое явление|оптическое свойство]] [[анизотропность|анизотропных]] материалов, в которых показатель преломления зависит от направления распространения света. В таких материалах может наблюдаться эффект расщепления [[Луч света|луча света]] на две составляющие, когда при попадании в материал образуется не один, а два преломленных луча с разным направлением и поляризацией.

Это явление описано датским учёным [[Бартолин, Расмус|Расмусом Бартолином]] на [[кристалл]]е [[кальцит|исландского шпата]] в [[1669 год]]у.

[[Файл:Birefringence.jpg|thumb|Прозрачный пластик на фоне поляризованного экрана при наблюдении через поляризатор. Цвета появляются из-за двулучепреломления материалом пластика. Объяснение эффекта в описании к файлу.]]

== Описание ==

=== Одноосные материалы ===
Простейший тип двулучепреломления наблюдается в [[Одноосные кристаллы|одноосных материалах]]. Чаще всего, это кристаллы, [[Кристаллическая решётка|решетка]] которых асимметрична, а именно вытянута или сжата в каком-либо [[Оптическая ось кристалла|направлении]]. При этом вращение вокруг этого направления (оптической оси) не меняет оптические свойства кристалла.
Поведение световой волны в такой среде зависит от направления распространения и поляризации света. Обыкновенной волной называется та, которая поляризована перпендикулярно оптической оси и направлению распространения, а поляризация необыкновенной волны перпендикулярна поляризации обыкновенной. Можно выделить три основных случая:

1) Свет распространяется вдоль оптической оси (при этом поляризация будет перпендикулярна оптической оси), тогда показатель преломления будет одинаковый для всех поляризаций, и кристалл в этом случае не отличается от изотропной среды, а между обыкновенной и необыкновенной волнами нет разницы.

 [[Файл:Линейная и круговая поляризация-ru.svg|thumb|После прохождения через четвертьволновую пластинку плоскополяризованное излучение превращается в излучение с круговой поляризацией]]
2) Свет распространяется перпендикулярно оптической оси. Тогда поляризацию можно разложить на две проекции — параллельную оптической оси и перпендикулярную. Эффективный показатель преломления будет разным для света двух ортогональных поляризаций, и при прохождении через слой (пластинку) материала может наблюдаться сдвиг по фазе между двумя компонентами. Если исходная поляризация линейная и ориентирована либо полностью вдоль, либо полностью перпендикулярно оптической оси, то на выходе из пластинки она не изменится. Однако, если исходно свет поляризован под углом к оптической оси, либо поляризация эллиптическая или циркулярная, то при прохождении через пластинку из одноосного кристалла поляризация может измениться из-за сдвига по фазе между компонентами. Сдвиг зависит от толщины пластинки, разницы между показателями преломления и длины волны света.

Пусть угол между поляризацией и оптической осью равен <math>45^o</math>. Если толщина пластинки такова, что на выходе из неё одна поляризация на четверть волны (четверть периода) отстаёт от другой, то исходная линейная поляризация превратится в круговую (такая пластинка называется четвертьволновой), если фаза одного луча отстанет от фазы другого луча на половину длины волны, то свет останется линейно поляризованным, но плоскость поляризации повернётся на некоторый угол, значение которого зависит от угла между плоскостью поляризации падающего луча и плоскостью главной оптической оси (такая пластинка называется полуволновой).

3) Свет распространяется в произвольном направлении относительно оптической оси. Тогда будет наблюдаться не один преломленный луч, а два, обладающих различной поляризацией. Направления преломленных лучей можно найти графически.

[[Файл:Birewiki1.svg|thumb|370px|Иллюстрация нахождения направления распространения обыкновенной и необыкновенной волн в одноосном кристалле]]

Математическое описание процесса является достаточно громоздким, но результат может быть наглядно проиллюстрирован с помощью построений, напоминающих иллюстрацию дифракции в кристалле с помощью [[сфера Эвальда|построения Эвальда]].

Пусть волна падает из воздуха на поверхность одноосного кристалла. Инструкция для нахождения направлений волнового и лучевого векторов для обыкновенной и необыкновенной волн для одноосного кристалла (см. рисунок, для простоты оптическая ось находится в плоскости падения). :

1. Нарисуйте поверхность кристалла горизонтально.

2. Нарисуйте в воздухе полусферу с радиусом, равным единице, и с центром, лежащим на поверхности кристалла.

2. Нарисуйте в среде полусферу с тем же центром, и радиусом, равным показателю преломления <math>n_o</math>.

3. Нарисуйте в среде эллипсоид с тем же центром, большая полуось которого ориентирована вдоль оптической оси кристалла и равна <math>n_o</math>, а малая - <math>n_e</math>.

4. Постройте падающий и отраженный лучи так, что конец падающего и начало отраженного находятся в центре сфер.

5. Проведите вертикальную линию, проходящую через пересечение отраженного луча со сферой.

6. Найдите точки пересечения линии со сферой и эллипсоидом в веществе.

7. Проведите из центра в точки пересечения направления волновых векторов обыкновенной и необыкновенной волн. Показатели преломления будут соответствовать длине этих векторов.

8. Для обыкновенной волны: вектор '''E''' должен быть перпендикулярен оптической оси и вектору '''k''', '''k'''||'''s'''.

9. Для необыкновенной волны: Лучевой вектор '''s''' должен быть перпендикулярен эллипсоиду в точке пересечения. Необыкновенный луч может не лежать в плоскости падения. Поляризация необыкновенной волны '''E''' перпендикулярна лучевому вектору '''s''' и поляризации обыкновенной волны. Вектор '''D''' перпендикулярен волновому вектору '''k'''. Векторы '''D''', '''E''', '''s''' и '''k''' необыкновенной волны должны лежать в одной плоскости<ref name="СПАУ_л11">{{cite web
 |url         = https://lms.physics.spbstu.ru/pluginfile.php/542/mod_resource/content/1/lect_18_m.pdf
 |title         = Электромагнитные волны. Волновое уравнение. Плоские волны. Поток энергии в плоской волне. Вектор Пойнтинга. Плотность потока импульса. Тензор напряжений. Световое давление. Опыты Лебедева.
 |author         = Д. А. Паршин, Г. Г. Зегря
 |work         = Электромагнитные волны. Лекция 18
 |access-date         = 2020-08-21
 |archive-date         = 2019-07-11
 |archive-url         = https://web.archive.org/web/20190711012455/http://lms.physics.spbstu.ru/pluginfile.php/542/mod_resource/content/1/lect_18_m.pdf
 |deadlink         = no
 }}</ref>.

<!--
Направление колебания [[Вектор (математика)|вектора]] электрического поля необыкновенного луча лежит в плоскости главного сечения (плоскости, проходящей через луч и оптическую ось кристалла). Оптическая ось кристалла — направление в оптически анизотропном кристалле, по которому луч света распространяется, не испытывая двойного лучепреломления.

Нарушение закона преломления света необыкновенным лучом связано с тем, что скорость распространения света (а значит и показатель преломления) волн с такой [[Поляризация электромагнитных волн|поляризацией]], как у необыкновенного луча, зависит от направления. Для обыкновенной волны скорость распространения одинакова во всех направлениях.

Можно подобрать условия, при которых обыкновенный и необыкновенный лучи распространяются по одной траектории, но с разными скоростями. Тогда наблюдается эффект изменения поляризации. Например, линейно поляризованный свет, падающий на пластинку, можно представить в виде двух составляющих (обыкновенной и необыкновенной волн), двигающихся с разными скоростями. Из-за разности скоростей этих двух составляющих на выходе из кристалла между ними будет некоторая разность фаз, и в зависимости от этой разности свет на выходе будет иметь разные поляризации. многострочный комментарий -->

=== Двуосные материалы ===
В таких кристаллах показатели преломления различны вдоль всех трех осей декартовой системы координат. Поверхность волновых векторов обладает сложной формой, но все еще существуют два выделенных направления, которые можно назвать оптическими осями, поскольку при распространении вдоль оптических осей существует только одно направление '''k'''- вектора. При этом этому направлению соответствует бесконечное множество лучевых векторов, заполняющих собой коническую поверхность, и наблюдается ''коническая рефракция''. При распространении вдоль направлений, не совпадающих с оптическими осями, наблюдается двойное лучепреломление, но в этом случае чаще всего оба луча являются необыкновенными (направление волнового и лучевого вектора не совпадает).

Двулучепреломление может наблюдаться не только в кристаллах, но и в любом материале с асимметричной структурой, например, в пластике.

== Природа явления ==
Качественно явление можно объяснить следующим образом. Из [[Уравнения Максвелла|уравнений Максвелла]] для материальной среды следует, что [[фазовая скорость]] света в среде обратно пропорциональна величине [[Диэлектрическая проницаемость|диэлектрической проницаемости]] ε среды. В некоторых кристаллах диэлектрическая проницаемость — [[тензор]]ная величина — зависит от направления электрического вектора, то есть от состояния [[Поляризация диэлектриков|поляризации]] волны, поэтому и фазовая скорость волны будет зависеть от её поляризации.

Согласно классической теории света, возникновение эффекта связано с тем, что переменное электромагнитное поле света заставляет колебаться [[электрон]]ы вещества, и эти колебания влияют на распространение света в среде, а в некоторых веществах заставить электроны колебаться проще в некоторых определённых направлениях.

== Вывод формул ==
В изотропной среде (включая свободное пространство) [[электрическая индукция]] ({{math|'''D'''}}) просто пропорциональна электрическому полю ({{math|'''E'''}}) в соответствии с {{math|'''D''' {{=}} ''ɛ'''''E'''}} где [[диэлектрическая проницаемость]] {{math|''ε''}} является просто [[скаляр]]ом (и равна {{math|''n''<sup>2</sup>''ε''<sub>0</sub>}} где {{math|''n''}} — [[показатель преломления]]). Однако, в анизотропных материалах соотношение между {{math|'''D'''}} и {{math|'''E'''}} должно описываться [[тензор]]ным уравнением:
{{EF|:|<math>\mathbf D = \boldsymbol\varepsilon \mathbf E </math>|ref=1}}

где {{math|'''ε'''}} теперь является матрицей 3&nbsp;×&nbsp;3.

Предположим что среда линейная и [[магнитная проницаемость]]: {{math|''μ'' {{=}} ''μ''<sub>0</sub>}}.

Запишем электрическое поле плоской волны с частотой {{math|''ω''}} в следующей форме:
{{EF|:|<math>\mathbf{E}=\mathbf{E}_0 e^{i(\mathbf{k} \cdot \mathbf{r}-\omega t)}</math>|ref=2}}

где {{math|'''r'''}} — радиус вектор, {{math|''t''}} — время, {{math|'''E'''<sub>0</sub>}} вектор, описывающий электрическое поле в {{math|'''r''' {{=}} 0}}, {{math|''t'' {{=}} 0}}. Найдем все возможные [[Волновой вектор|волновые векторы]] {{math|'''k'''}}. Комбинируя [[уравнения Максвелла]] для {{math|∇ × '''E'''}} и {{math|∇ × '''H'''}}, и исключая {{math|'''H''' {{=}} {{sfrac|1|''μ''<sub>0</sub>}}'''B'''}}, получаем:
{{EF|:|<math>-\nabla \times \nabla \times \mathbf{E}= \mu_0\frac{\partial^2 }{\partial t^2} \mathbf{D}</math> |ref=3a}}

Вспомним также, что в отсутствие свободных зарядов, дивергенция {{math|'''D'''}} исчезает:
{{EF|:|<math>\nabla \cdot \mathbf{D}=0 </math>|ref=3b}}

Применим [[Оператор набла#Операторы второго порядка|соотношение]] {{math|∇ × (∇ × '''A''') {{=}} ∇(∇ ⋅ '''A''') − ∇<sup>2</sup>'''A'''}} к левой части {{eqref|3a|3a}}, и воспользуемся тем, что поле представляет собой плоскую волну, а значит производная по {{math|''x''}} (например) приводит к умножению на {{math|''ik<sub>x</sub>''}}:
{{EF|:|<math> - \nabla \times \nabla \times \mathbf{E} =(\mathbf{k} \cdot \mathbf{E}) \mathbf{k} - (\mathbf{k} \cdot \mathbf{k}) \mathbf{E}</math>|{{EquationRef|3c}}}}

Правая часть {{eqref|3a|3a}} может быть выражена через {{math|'''E'''}} с помощью тензора {{math|'''ε'''}}, а производные по времени просто приводят к умножению на {{math|−''iω''}}, и тогда {{eqref|3a|3a}}:
{{EF|:|<math>(-\mathbf{k}\cdot\mathbf{k})\mathbf{E} + (\mathbf{k} \cdot \mathbf{E}) \mathbf{k}= -\mu_0 \omega^2 (\boldsymbol{\varepsilon} \mathbf{E})</math>|ref=4a}}

Применяя дифференцирование к {{eqref|3b|3b}} находим:
{{EF|:|<math>\mathbf{k} \cdot \mathbf{D} =0</math>|ref=4b}}

{{eqref|4b|Уравнение 4b}} означает, что {{math|'''D'''}} перпендикулярно направлению волнового вектора {{math|'''k'''}}, при этом это больше не верно для вектора {{math|'''E'''}} как это было бы в изотропной среде. {{eqref|4b|Уравнение 4b}} не будет использовано в дальнейшем.

Найти допустимые значения вектора {{math|'''k'''}} для данной {{math|''ω''}} проще всего в [[Прямоугольная система координат|декартовой системе координат]], в которой оси {{math|''x''}}, {{math|''y''}} и {{math|''z''}} параллельны осям симметрии кристалла (или просто выбирая ось {{math|''z''}} вдоль оптической оси одноосного кристалла). Тогда матрица для тензора {{math|'''ε'''}} будет диагональной:
{{EF|:|<math>\mathbf{\varepsilon}= \varepsilon_0 \begin{bmatrix} n_x^2 & 0 & 0 \\ 0& n_y^2 & 0 \\ 0& 0& n_z^2 \end{bmatrix} </math>|ref=4c}}

на диагонали стоят квадраты показателя преломления для поляризаций вдоль осей {{math|''x''}}, {{math|''y''}} и {{math|''z''}}. Подставляя {{math|'''ε'''}} в этой форме, и скорость света {{math|''c''}} в виде {{math|''c''<sup>2</sup> {{=}} {{sfrac|1|''μ''<sub>0</sub>''ε''<sub>0</sub>}}}}, Проекция векторного {{eqref|4a|уравнения 4a}} на ось ''x'' записывается как
{{EF|:|<math>\left(-k_x^2-k_y^2-k_z^2 \right)E_x + k_x^2 E_x + k_xk_yE_y + k_xk_zE_z = -\frac{\omega^2n_x^2}{c^2} E_x</math>|ref=5a}}

где {{math|''E<sub>x</sub>''}}, {{math|''E<sub>y</sub>''}}, {{math|''E<sub>z</sub>''}} компоненты вектора {{math|'''E'''}} и {{math|''k<sub>x</sub>''}}, {{math|''k<sub>y</sub>''}}, {{math|''k<sub>z</sub>''}} компоненты волнового вектора {{math|'''k'''}}. Запишем уравнения для всех трех проекций {{eqref|4a|ур. 4a}}:
{{EF|:|<math>\left(-k_y^2-k_z^2+\frac{\omega^2n_x^2}{c^2}\right)E_x + k_xk_yE_y + k_xk_zE_z =0</math>|ref=5b}}
{{EF|:|<math>k_xk_yE_x + \left(-k_x^2-k_z^2+\frac{\omega^2n_y^2}{c^2}\right)E_y + k_yk_zE_z =0</math>|ref=5c}}
{{EF|:|<math>k_xk_zE_x + k_yk_zE_y + \left(-k_x^2-k_y^2+\frac{\omega^2n_z^2}{c^2}\right)E_z =0</math>|ref=5d}}

Это система линейных уравнений на {{math|''E<sub>x</sub>''}}, {{math|''E<sub>y</sub>''}}, {{math|''E<sub>z</sub>''}}, которая имеет нетривиальное решение (т.е. {{math|'''E''' {{=}} 0}}) только если [[определитель]] следующей матрицы равен нулю:
{{EF|:|<math>\begin{vmatrix}
\left(-k_y^2-k_z^2+\frac{\omega^2n_x^2}{c^2}\right) & k_xk_y & k_xk_z \\
k_xk_y & \left(-k_x^2-k_z^2+\frac{\omega^2n_y^2}{c^2}\right) & k_yk_z \\
k_xk_z & k_yk_z & \left(-k_x^2-k_y^2+\frac{\omega^2n_z^2}{c^2}\right) \end{vmatrix} =0</math>|ref=6}}

Вычисляя определитель {{eqref|6|6}}, получаем
{{EF|:|<math> \frac{\omega^4}{c^4} - \frac{\omega^2}{c^2}\left(\frac{k_x^2+k_y^2}{n_z^2}+\frac{k_x^2+k_z^2}{n_y^2}+\frac{k_y^2+k_z^2}{n_x^2}\right) + \left(\frac{k_x^2}{n_y^2n_z^2}+\frac{k_y^2}{n_x^2n_z^2}+\frac{k_z^2}{n_x^2n_y^2}\right)\left(k_x^2+k_y^2+k_z^2\right)=0 </math>|ref=7}}

Уравнение {{eqref|7|7}} также называется уравнением Френеля.

=== Одноосный кристалл ===
При этом в случае одноосного материала (два диагональных элемента матрицы {{math|'''ε'''}} равны друг другу), и выбирая систему координат так, что оптическая ось направлена вдоль {{math|''z''}}, обозначим {{math|''n<sub>x</sub>'' {{=}} ''n<sub>y</sub>'' {{=}} ''n''<sub>o</sub>}} и {{math|''n<sub>z</sub>'' {{=}} ''n''<sub>e</sub>}}, выражение сводится к
{{EF|:|<math>\left(\frac{k_x^2}{n_\mathrm{o}^2}+\frac{k_y^2}{n_\mathrm{o}^2}+\frac{k_z^2}{n_\mathrm{o}^2} -\frac{\omega^2}{c^2}\right)\left(\frac{k_x^2}{n_\mathrm{e}^2}+\frac{k_y^2}{n_\mathrm{e}^2}+\frac{k_z^2}{n_\mathrm{o}^2} -\frac{\omega^2}{c^2}\right)=0</math>|ref=8}}
[[Файл:Uniaxial birefringence.svg|thumb|Поверхности волновых векторов в одноосном кристалле. Слева — положительный, справа — отрицательный.]]
Чтобы уравнение {{eqref|8|8}} выполнялось, один из множителей должен быть равен нулю. Заметим, что первый соответствует уравнению сферы, а второй — поверхности [[эллипсоид]]а в пространстве волновых векторов {{math|'''k'''}} для заданной {{math|''ω''}}. Первый множитель соответствует решению для обыкновенной волны, где показатель преломления равен {{math|''n''<sub>o</sub>}} независимо от направления, а второй — для необыкновенного. Второй множитель соответствует решению для необыкновенной волны, где эффективный показатель преломления меняется от {{math|''n''<sub>o</sub>}} до {{math|''n''<sub>e</sub>}} в зависимости от направления {{math|'''k'''}}. Для произвольного направления распространения волны возможны два вектора {{math|'''k'''}}, соответствующие двум разным поляризациям.

Для обыкновенной волны вектора {{math|'''D'''}} и {{math|'''E'''}} совпадают, а также совпадают направления волнового вектора {{math|'''k'''}} и направление [[Световой луч|лучевого вектора]] {{math|'''s'''}} в геометрической оптике (направление которого совпадает с вектором [[Групповая скорость|групповой скорости]] <math> \mathbf{v}_{gr} = \nabla_\mathbf{k} \omega</math>). Для необыкновенной волны это в общем случае не так. Рассмотрим уравнение для одноосного кристалла
{{EF|:|<math>F(k_x,k_y,k_z)=\frac{k_x^2}{n_\mathrm{e}^2}+\frac{k_y^2}{n_\mathrm{e}^2}+\frac{k_z^2}{n_\mathrm{o}^2} -\frac{\omega^2}{c^2}=0</math>|ref=9}}.
Сравним уравнение для групповой скорости <math>\mathbf{v}_{gr} = \nabla_\mathbf{k} \omega </math> с [[Нормаль#Нормаль к плоской кривой|уравнением нормали к поверхности]], заданной неявно. Поскольку уравнения совпадают с точностью до константы, лучевой вектор перпендикулярен рассматриваемому эллипсоиду.

=== Двуосный кристалл ===
Чтобы понять, как выглядит поверхность в случае, когда все диагональные элементы матрицы матрицы {{math|'''ε'''}} разные (пусть <math>n_z>n_y>n_x</math>), положим одну из компонент вектора '''k''' равной нулю (<math>k_z=0</math>) и перепишем уравнение {{eqref|7|7}}.
{{EF|:|<math> \frac{\omega^4}{c^4} - \frac{\omega^2}{c^2}\left(\frac{k_x^2+k_y^2}{n_z^2}+\frac{k_x^2}{n_y^2}+\frac{k_y^2}{n_x^2}\right) + \left(\frac{k_x^2}{n_y^2}+\frac{k_y^2}{n_x^2}\right)\frac{k_x^2+k_y^2}{n_z^2}=0 </math>|ref=10}}
Его можно разложить на множители:
{{EF|:|<math>\left( \frac{\omega^2}{c^2} - \frac{k_x^2}{n_y^2}-\frac{k_y^2}{n_x^2} \right)\left(\frac{\omega^2}{c^2}-\frac{k_x^2+k_y^2}{n_z^2}\right)=0 </math>|ref=11}} [[Файл:Conicalwiki.png|thumb|Поверхность волновых векторов в двуосном кристалле. Коническая рефракция.]]
Первый множитель представляет собой эллипс, а второй — окружность. Аналогичное разложение можно проделать для всех трех плоскостей <math>k_i=0</math>. На рисунке изображены сечения поверхностью трех координатных плоскостей в одном октанте, в остальных картина симметрична. Поверхность имеет 4 особые точки (точки самопересечения), в нашем случае лежащие в плоскости ''xz''. Через эти точки проходят две оси <math>\beta</math>, которые называются ''оптическими осями'' (либо ''бинормалями'') двуосного кристалла. Только в этих направлениях волновой вектор может иметь единственное значение. Однако, в особой точке поверхности направление нормали является неопределенным, а лучевой вектор может заполнять коническую поверхность (конус внутренней ''конической рефракции'')

== Искусственное двойное лучепреломление ==
Помимо двулучепреломляющих кристаллов двойное лучепреломление наблюдается и в [[изотропность|изотропных]] средах, помещённых в электрическое поле ([[эффект Керра]]), в магнитное поле ([[эффект Фарадея]] и [[эффект Коттона — Мутона]]), под действием механических напряжений ([[фотоупругость]]). Под действием этих факторов изначально изотропная среда меняет свои свойства и становится анизотропной. В этих случаях оптическая ось среды совпадает с направлением электрического поля, магнитного поля, направлением приложения силы.

== Положительные и отрицательные кристаллы ==
* '''Отрицательные кристаллы''' — одноосные кристаллы, в которых скорость распространения обыкновенного луча света меньше, чем скорость распространения необыкновенного луча. В кристаллографии отрицательными кристаллами называют также жидкие включения в кристаллах, имеющие ту же форму, что и сам кристалл.
* '''Положительные кристаллы''' — одноосные кристаллы, в которых скорость распространения обыкновенного луча света больше, чем скорость распространения необыкновенного луча.

== См. также ==
{{commons-inline|Birefringence}}
* [[Поляризация диэлектриков]]
* [[Эффект Коттона — Мутона]]
* [[Эффект Керра]]
* [[Эффект Поккельса]]
* [[Эффект Фарадея]]

== Литература ==
* {{Книга:Сивухин Д.В.:Оптика:2002}}
* ''Ландсберг Г. С.'' Оптика М., 2004 г.
* {{Книга:Ландау Л. Д., Лифшиц Е. М.: Электродинамика сплошных сред|2003}}
* ''Салех Б., Тейх М.'', Оптика и фотоника. Принципы и применения, пер. с англ. В 2 т.

== Примечания ==
{{примечания}}

== Ссылки ==
* Erasmus Bartholin, [https://books.google.com/books?id=F7RAAAAAcAAJ&pg=PP11#v=onepage&q&f=false Experimenta crystalli islandici disdiaclastici quibus mira & infolita refractio detegitur] (Copenhagen, Denmark: Daniel Paulli, 1669).
* Erasmus Bartholin (January 1, 1670) [https://books.google.com/books?id=Zl5FAAAAcAAJ&pg=PA2141#v=onepage&q&f=false An account of sundry experiments made and communicated by that learn'd mathematician, Dr. Erasmus Bartholin, upon a chrystal-like body, sent to him out of Island,] ''Philosophical Transactions of the Royal Society of London'', '''5''' : 2041—2048.
[[Категория:Оптические явления]]