Расчёт диэлектрической проницаемости исп. ур Клаузиуса—Мосотти. Всегда получается единица О_о

@dan24 · Регистрация: 19.01.2013

Author24 — интернет-сервис помощи студентам

Что есть:
1) Есть программа Gaussian. В ней есть изучаемый кластер в 100 атомов ( плоская структура типа графена в узлах которой иногда встречается азот. В середине кластера помещено два атома ферромагнитных металлов ( Ni и Co, например).
2) Для данного кластера получен в программе Gaussian дипольный момент $https://www.cyberforum.ru/cgi-bin/latex.cgi?\mu$ , поляризуемость $https://www.cyberforum.ru/cgi-bin/latex.cgi?\alpha$ и гиперполяризуемость $https://www.cyberforum.ru/cgi-bin/latex.cgi?\beta$ . Из этого всего нам понадобится только поляризуемость $https://www.cyberforum.ru/cgi-bin/latex.cgi?\alpha$ .

Я прочитал статью: Ссылка на DOI статьи Также см. статью в приложении 1.
Там делается мега финт ушами и применяется к смолам т.е. к твёрдым веществам ур Клаузиуса—Моссотти, предназначенное для газов.
Мне очень надо посчитать диэлектрическую проницаемость любым способом с любой погрешностью ( потом разберём как уменьшить погрешность) и я решил сделать также как с статье.

Вводная к вопросу ( можно и не читать):
Ур Клаузиуса—Моссотти:
$https://www.cyberforum.ru/cgi-bin/latex.cgi?(\epsilon-1)/(\epsilon+2)=(N*\alpha )/(3*{\epsilon }_{0})$
Подробнее: ОСНОВЫ ДИЭЛЕКТРИЧЕСКОЙ СПЕКТРОСКОПИИ Приложение 2.
Ещё более подробно: Физика диэлектриков, Поплавский Ю.М. 1980г. Приложение 3.
При этом в зарубежной статье формула имеет вид:
$https://www.cyberforum.ru/cgi-bin/latex.cgi?(\epsilon-1)/(\epsilon-2)=(N*\alpha )/(3*{\epsilon }_{0})$ ( эпсилон не +2, а -2 ). И следующие формулы 2,3 на стр 525 также содержат ошибки ( как мне показалось).

Вопрос.
Пусть

N*\alpha )/(3*{\epsilon }_{0})=l[/LATEX] , чтобы было удобно смотреть, тогда получим то, что на картинке 1

Что получается:
Какой бы ни был параметр l ( он там немного по другому считается, это не суть. в любом случае не получается) всегда выходит значение
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\epsilon=(2*l+1)/(l-2) =1$
Всегда!!!!!
Если брать в инзачальной форме $https://www.cyberforum.ru/cgi-bin/latex.cgi?\epsilon +2$ , то всегда получается 2. Короче как ни крути всегда выходит константа. Что за дела такие?!

( Возможно вопрос больше по самому мэплу. Пока не понятно то ли я как-то не правильно по формуле рассчитываю, то ли мэпл касячит).

Мэпловский файл в архиве.

Я искал ответ на этот вопрос где-то неделю. Я написал сюда от безысходности.

Очевидно, что что-бы мы не вычислили в правой части исходного уравнения, в СИ эта величина будет порядка 10^-23. Вычитая из этой мелочи единицу MAPLE получает единицу.
В общем я не знаю можно ли так делать, но если перевести все величины в сгс, а не в СИ, то величина справа в уравнении Клаузиуса—Мосотти становится порядка 10^0 ( Не знаю какой это порядок. Нормальный короче ).
Например у меня получилось число справа 1.24.
Далее вычисляя диэлектрическую проницаемость мы получаем число отличное от 1 или 2. Я получил 6.
В чём прелесть?! В том, что диэлектрическая проницаемость - безразмерная величина. И в этом я считаю есть читерство и петушковая магия.
В СИ мы вычитаем из единицы 1.380649000*10^(-23), а в СГС 1.24. Я думаю разница огромна.

Не по теме:

Короче не знаю как работает эта физика.... Мне всё равно деньги не дадут на реальные эксперименты. Так что вообще всё равно как там в реальности)))). Но очень интересно как так выходит то....
В общем вопрос типа решён(.

	23.11.2019, 16:03