Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

MaltsevSciD

.pdf
Скачиваний:
19
Добавлен:
28.03.2016
Размер:
4.49 Mб
Скачать

Интенсивность, отн. ед.

2.0

1.5

1.0

0.5

0.010

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

20

 

 

 

 

 

 

 

 

 

25

 

 

 

 

 

 

30

35

40

45

 

 

 

Угол рассеяния, градусы

 

 

 

 

Рис. 3.14. Индикатриса модифицированного цилиндра.

91

 

1

 

 

 

 

 

 

 

 

 

 

0.1

 

β = 900

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.01

 

 

 

 

 

 

 

 

 

units

1E-3

 

 

 

β = 600

 

 

 

 

 

 

 

 

 

 

 

 

 

arb.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Intensity,

1E-4

 

 

 

 

 

 

 

 

 

1E-5

 

 

 

 

 

 

 

 

 

 

1E-6

 

β = 00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1E-7

 

 

 

 

 

 

 

 

 

 

16

18

20

22

24

26

28

30

32

34

Scattering angle, degrees

Рис. 3.15. Индикатрисы эритроцита рассчитанные по ВКБ приближению для разной ориентации относительно направления падающего излучения.

92

Рассмотрим еще один показательный пример. Рассчитаем индикатрисы одиночного эритроцита в зависимости от его ориентации относительно падающего излучения. Форма эритроцита задается формулой [101]:

z2=(0.86)2(1-x2)(0.01384083+0.2842917x2+0.01306932x4). (3.16)

Это уравнение описывает профиль эритроцита в плоскости перпендикулярной оси симметрии эритроцита. Типичный объем эритроцита равен 85 мкм3, что соответствует длине 6.5 мкм вдоль длинной оси эритроцита и размеру перетяжки 1.3 мкм. Реальная часть показателя преломления эритроцита лежит между 1.40 и 1.42 для длины волны λ =0.6328 мкм [14, 102].

Расчет проводился по формуле ВКБ приближения (3.11). Результат показан на Рис. 3.15. Как видно из представленных данных, контраст индикатрисы очень чувствителен к ориентации частицы.

3.3.3.2. Функция распределения плотности набега фазы частицы.

93

Рис. 3.16. Геометрия рассеяния света на частице произвольной формы.

94

На основе результатов предыдущего раздела попытаемся найти зависящую от характеристик частицы функцию, параметры которой влияют на контраст индикатрисы. В разделе 3.3.2.2 было показано, что контраст индикатрисы сферической частицы зависит в основном от параметра набега фазы ρ=(2π/λ) dn0(m-1), где m=n/n0, n - показатель преломления частицы, n0 - показатель преломления среды, d - диаметр шара. Рассмотрим функцию распределения набега фазы по сечению частицы. Для этого поместим частицу в систему координат (Рис. 3.16). Выберем направление по оси Z, вдоль которого будем считать величину ρ=(2π/λ) (z2-z1)n0(m-1) (где z1 и z2 - координаты точек пересечения прямой, параллельной оси Z, с поверхностью частицы). Таким образом, зная форму частицы и зависимости z1 и z2 от координат x и y, можно найти функцию ρ = ρ(x, y). Из её свойств выделим, что она однозначная и определена только для точек проекции частицы на плоскость XY. Набег фазы не может быть больше 2π . Следовательно, окончательный вид ρ(x, y) будет такой:

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

(3.17)

 

 

 

 

 

 

 

 

 

 

 

 

ρ(x,y) = mod

 

n0 (z2

(x,y) z1

(x,y))(m − 1), 2

π

,

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где функция mod(a, b) означает остаток от деления а на b.

Теперь, зная функцию ρ(x, y),

можно найти её градиент на плоскости

 

 

 

 

 

 

 

 

 

 

!

 

∂ρ ∂ρ

(конечно, в тех точках, где он существует) : ρ(x,y) =

 

 

,

.

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

Тогда предлагается ввести следующую функцию распределения плотности

набега фазы частицы:

 

 

 

 

 

 

 

 

 

 

 

 

 

const

 

 

 

 

 

 

 

 

(3.18)

ν (ρ) =

 

 

 

 

 

 

 

 

 

 

 

min{

 

! ρ'(x,y)}

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ'= ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где в числителе стоит нормировочная константа, а выражение в

знаменателе означает, что если существуют точки (x, y), в которых набег фазы

!

одинаковый, то берётся та точка, в которой ρ минимален.

Можно определить характеристики данного распределения. Так средняя величина распределения определяется выражением

95

ρmax (3.19)

ρ= ρν(ρ)dρ,

ρmin

а ширина функции ν (ρ) вычисляется по формуле

ρmax (3.20)

∆ ν = ρ− ρν (ρ)dρ,

ρmin

3.3.3.3. Функция распределения плотности набега фазы частиц разной формы.

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

всех точек проекции цилиндра на плоскость XY набег фазы одинаков.

!

Следовательно, ρ для любой пары (x, y) равен нулю. Функция ν (ρ) определена только в точке ρ = (2π/λ) hn0(m-1) (где h- высота цилиндра) и равна в ней бесконечности. Ширина этой функции равна нулю (Рис. 3.17 a)). Для модифицированного цилиндра и куба (Рис. 3.11 c)) и (Рис. 3.11 d)) аналогичная ситуация: набег фазы для любой точки проекции фигуры на плоскость XY имеет одну величину. Следовательно, ширина функции распределения равна нулю.

Теперь рассмотрим

сдвоенный конус, для которого величина ρ = 2(a-

(x2+y2)1/2)(2π/λ) n0(m-1), где а- радиус и высота одного конуса. Тогда

 

 

2x

 

2y

 

 

 

(3.21)

! ρ(x,y) =

 

 

 

 

 

,

2

kn

0

(m − 1) ,

 

 

x

2

+ y

2

x

2

+

y

 

 

 

 

 

 

 

 

 

 

 

!

 

 

4x2

 

 

 

4y2

 

 

(3.22)

=

 

 

+

 

kn0 (m − 1) = 2kn0 (m − 1) ,

ρ(x,y)

x2 +

y2

 

x

2 +

 

 

 

 

 

 

y2

 

 

И функция ν (ρ), нормированная на единицу, равна 1/(2kn0(m-1)). По формуле (3.20) находим ширину ∆ν = akn0(m-1). На Рис. 3.17 b) приведены графики функций для сдвоенных конусов с а=0.8 мкм, но с разными показателями преломления. Сравнивая графики на Рис. 3.13 и Рис. 3.17 b),

96

можно сделать вывод, что при увеличении ширины функции распределения ν (ρ) контраст индикатрисы уменьшается, доходя до нуля.

Для сферической частицы ρ(x,y) =

2

a2

 

x2

 

y2 kn0 (m − 1) , и

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

(3.23)

! ρ(x,y) =

 

 

2x

 

 

 

 

 

 

2y

 

 

 

 

 

 

 

2

 

 

2

 

+

 

 

2

 

 

2

 

kn

0

(m − 1) =

 

 

a

x

2

y

 

 

 

a

x

2

y

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ2

− ρ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2kn0 (m − 1)

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где ρ0 =

(2π/λ)

2an0(m-1)

 

-

максимальный набег фазы луча, проходящего

через центр сферической частицы. Функция распределения, нормированная на единицу:

 

1/ρ0

(3.24)

ν (ρ) =

,

(ρ /ρ)2

 

− 1

 

0

 

На Рис. 3.18 показаны графики функций распределения ν (ρ) для сферических частиц, у которых максимальные набеги фаз ρ0 больше и меньше

2π . Ширина этой функции в зависимости от ρ0, посчитанная по формуле (3.20), показана на Рис. 3.19.

97

ν(ρ)

8

6

4

2

0

0

1

2

3

4

5

6

7

 

 

 

Набег фазы

 

 

 

 

 

 

 

 

 

 

 

 

a)

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1.34

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ν(ρ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,5

1,0

1,5

 

2,0

2,5

 

3,0

 

0,0

 

 

 

Набег фазы

b)

Рис. 3.17. Функция распределения плотности набега фазы для цилиндра a) и сдвоенного конуса b).

98

Частица, имеющая форму эритроцита (формула (3.16), описывающая его поверхность, немного изменена для 0 r 3.15, где 3.15 мкм - настоящий радиус клетки эритроцита, m = 1.41 - его показатель преломления, показатель преломления среды - 1.333), обладает функцией распределения плотности набега фазы, которая представлена на Рис. 3.20. Рассчитанная по формуле (3.20) ширина функции плотности набега фазы и средняя величина набега фазы равны:

∆ν

= 0.221, ρ = 0.813

(3.25)

 

99

 

1,4

 

 

 

 

 

 

 

 

1,2

 

Шар

 

 

 

 

 

 

 

ρ0<2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,0

 

ρ0>2π

 

 

 

 

 

ν(ρ)

0,8

 

 

 

 

 

 

 

0,6

 

 

 

 

 

 

 

 

0,4

 

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

 

 

0,00

1

2

3

4

5

6

7

 

 

 

Набег фазы

 

 

 

Рис. 3.18. Функция распределения плотности набега фазы для сферических частиц.

100

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]