к библиотеке   к оглавлению   FAQ по эфирной физике   ТОЭЭ   ТЭЦ   ТПОИ   ТИ  

РЕАЛЬНАЯ ФИЗИКА

Глоссарий по физике

А   Б   В   Г   Д   Е   Ж   З   И   К   Л   М   Н   О   П   Р   С   Т   У   Ф   Х   Ц   Ч   Ш   Э   Ю   Я  

Спиновые волны

Спиновые волны - волны нарушений магнитной упорядоченности в ферро- , антиферро- и ферримагнетиках. Спины атомов в этих веществах и связанные с ними магн. моменты в осн. состоянии упорядочены. Отклонение магн. момента от преимущественного направления не локализуется на атоме, а в виде волны распространяется в среде.

Спиновые волны - элементарное возбуждение магн. системы в магнитоупорядоченной среде; квазичастицы, соответствующие спиновые волны, наз. магнонами. Существование спиновых волн в ферромагнетиках предсказано Ф. Блохом (F. Bloch) в 1930. Вся совокупность экспе-рим. факторов о поведении магнитоупорядоченных тел при температуре Т значительно ниже температуры Кюри Тс (или температуры Нееля TN)свидетельствует о существовании спиновых волн (в частности, Блоха закон).

Спиновые волны, как всякая волна в кристалле, характеризуется законом дисперсии - зависимостью её частоты w от квазиволнового вектора k. Энергия8060-29.jpg и квазиимпульс р магнона равны:8060-30.jpg ,8060-31.jpg. Кристалл с N магнитными под решетками имеет N типов (ветвей, мод) спиновых волн с разл. законами дисперсии: ш = wi(k); i = 1, 2, ..., N.

Классическое описание. Спиновые волны допускают наглядную классич. интерпретацию: рассмотрим цепочку атомов, расстояние между к-рыми а, в магн. поле Н. Если волновой вектор k = 0, то все спины синфазно прецессируют вокруг Н с частотой w0 (однородная прецессия). При8060-32.jpg прецессия спинов неоднородна - разные спины повёрнуты на разные углы, разность углов поворота равна ka (рис. 1). Частота неоднородной прецессии8060-34.jpg . В реальных системах малые колебания магн. моментов атомов осуществляются в виде волн неоднородной прецессии.
8060-33.jpg

Рис. 1. Спиновая волна в линейной цепочке спинов: a - вид цепочки спинов в перспективе (сбоку); б - вид цепочки спинов сверху; волна изображена линией, проходящей через концы спиновых векторов.

В случае длинных волн колебания магн. моментов можно описывать как колебания макроскопич. векторов - плотностей магн. моментов (намагниченностей)подрешёток Mi (r,t) - функций координаты r и времени t. При неоднородной прецессии длины векторов8060-35.jpg8060-36.jpg, где8060-37.jpg - магн. момент атома i-й подрешётки, v0 - объём ячейки кристалла, сохраняются;8060-38.jpg - интегралы движения. Законы дисперсии длинноволновых спиновых волн определяются из квантового уравнения:
8060-39.jpg

Здесь8060-40.jpg - магнитомеханическое отношение (без учёта спин-орбитальной связи); е, т - заряд и масса электрона,8060-41.jpg - энергия взаимодействия магн. моментов подрешёток, Нэф - эфф. магн. поле (см. ниже). Осн. состояние определяется условием коллинеарности намагниченностей Mi0 и эфф. магн. полей8060-42.jpg:
8060-43.jpg

Линеаризация ур-ния (1) с учётом (2) приводит к системе ур-ний для пере.м. составляющих магн. моментов:
8060-44.jpg

Поле Нэф, кроме пост. магн. поля Н, содержит перем. часть h - магн. поле, связанное с взаимодействием между подрешётками и с неоднородностью их намагниченностей. Т. к. частоты спиновых волн невелики, то для определения h можно воспользоваться ур-ниями магнитостатики:
8060-45.jpg

Магн. поле h осуществляет магнитодипольное взаимодействие между колеблющимися магн. моментами. Линеаризов. ур-ния (1) совместно с (4) суть ур-ния спиновых волн. Подставляя в них
8060-46.jpg

получаем алгебраич. систему ур-нии относительно амплитуд спиновых волн mi0. Равенство нулю детерминанта этой системы приводит к ур-нию N-го порядка относительно w2. Его решения определяют законы дисперсии спиновых волн при8060-47.jpg

Обычно в магнитоупорядоченных средах гл. роль во взаимодействии между магн. моментами атомов играет обменное взаимодействие ,изотропное относительно однородного поворота магн. моментов атомов. Магн. порядок появляется в результате спонтанного нарушения симметрии обменного взаимодействия. Энергия обменного взаимодействия соседних атомов8060-48.jpg порядка температуры Кюри Тс (температуры TN); знак J выбирается так, что при J > 0 обменное взаимодействие благоприятствовало бы ферромагн. упорядочению, а при J < 0 - антиферромагнитному.

Ветви спиновых волн. Число ветвей спиновых волн равно числу магн. подрешёток. Это обусловлено прецессионным характером движения магн. моментов подрешёток. Ветви спиновых волн принято делить на акустические и оптические аналогично колебаниям кристаллической решётки. Если пренебречь малыми (по сравнению с обменными), т. н. релятивистскими, взаимодействиями (зеемановским с постоянным магн. полем, спин-орбитальным - источником энергии магнитной анизотропии, магнитодипольным и др.), то акустич. типы спиновых волн представляют собой голдстоуновские моды ,т. е. в их энергетич. спектре при k = 0 щель отсутствует. Частоты акустич. спиновых волн стремятся к 0 с ростом длины волны8060-49.jpg . Их число и характер закона дисперсии w(k) при8060-50.jpg зависят от структуры осн. состояния магнетика, причём при любом кол-ве подрешёток число акустич. мод8060-51.jpg3. У одноподрешёточного ферромагнетика одна акустич. мода с8060-52.jpgпри8060-53.jpg; у двухподрешёточного антиферромагнетика 2 вырожденные акустич. моды с8060-54.jpg . В ферромагнетике магнон напоминает нерелятивистскую частицу с энергией8060-55.jpg, в антиферромагнетике - акустич. фонон с8060-56.jpg (т, и - масса частицы и скорость звука). Примеры магнетиков, имеющих 3 акустич. ветви в спектре спиновых волн, - многоподрешёточные антиферромагнетики с неколлинеарным расположением магн. моментов в упорядоченном состоянии при Н = 0 (UO2, CsNiCl3, CsMnBr3 и др.). Учёт релятивистских взаимодействий приводит к возникновению энергетич. щелей в спектре акустич. ветвей спиновых волн 8060-57.jpg (8060-58.jpg - частоты однородной прецессии). Когда в спектре спиновых волн есть оптич. моды, их частоты однородной прецессии8060-59.jpg

Дисперсия. спиновых волн являются причиной зависимости тензора магнитной восприимчивости8060-60.jpg от волнового вектора8060-61.jpg (см. Дисперсия пространственная ).Частотная дисперсия (зависимость8060-62.jpg от w) является следствием прецессии магн. моментов подрешёток. Тензор8060-63.jpg определяется в результате решения ур-ния (1), а Максвелла уравнения дают возможность найти связь между ш и k, т. е. законы дисперсии спиновых волн, учитывающие конечность скорости света. При8060-64.jpg они отличаются от законов дисперсии, полученных на основе ур-ний магнитостатики (4), малыми поправками, к-рые иногда существенны, напр. при описании взаимодействия спиновых волн с электронами проводимости в металлах и полупроводниках.

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

С ростом k (ak ~ 1) проявляется дискретная (кристаллич.) структура магнетиков. Для получения законов дисперсии, справедливых при произвольном значении ak, обычно используют приближённые представления спиновых операторов8060-65.jpg через операторы рождения8060-66.jpg и уничтожения8060-67.jpg магнонов, подчиняющиеся бозевским правилам коммутации (преобразование Хольштейна - Примакова):
8060-68.jpg

Здесь индекс I нумерует атомы, координатные оси выбраны так, чтобы ось z для каждого атома была направлена вдоль равновесного положения спина. Из правил коммутации для , al следует, что - любое целое число от 08060-69.jpg до8060-70.jpg, хотя по физич.8060-71.jpg смыслу8060-72.jpg Вблизи основного состояния ср. значение nl значительно меньше sl и приближённые ф-лы (6) пригодны для вычисления спектра тем точнее, чем больше sl (в квантовомеханич. пределе8060-73.jpg ). Однако и при sl>l ~ 1 частоты спиновых волн, как правило, лишь небольшими поправками отличаются от значений, найденных с помощью (6).

Магнонный спектр. Теоретич. рассмотрение позволяет вычислить энергию магнонов при любом k. Это приводит к периодич. зависимости
8060-74.jpg

где Ь - произвольный вектор обратной решётки. Так, гамильтониан одноподрешёточного магнетика
8060-75.jpg

Здесь J(Rlm) - обменный интеграл между l-м и m-м атомами, Rlm - вектор, соединяющий эти атомы,8060-76.jpg - магн. момент атома. С помощью (7) и (8) (пренебрегая взаимодействием между магнонами) можно получить спектр магнонов:
8060-77.jpg

Ширина магнонной энергетич. зоны8060-78.jpg -8060-79.jpg , где8060-80.jpg, равна:
8060-81.jpg

(J - обменный интеграл для ближайших соседей). Соотношение8060-82.jpg - общее свойство магнонных зон. Магнитный момент магнона. Зависимость энергии магнона от магн. поля Н означает, что магнон обладает магн. моментом:
8060-83.jpg

В простейшем случае чисто обменного одноподрешёточного ферромагнетика магн. момент магнона равен магн. моменту атома и направлен против равновесной намагниченности. Увеличение числа магнонов приводит к уменьшению величины спонтанной намагниченности магнетика. В многоподрешёточных магнетиках рост числа магнонов уменьшает намагниченность подрешёток.

В магн. металлах (Fe, Co, Ni и др.), где за магн. свойства ответственны d-электроны, в формировании спектра спиновых волн принимают участие нелокализов. электроны проводимости .В длинноволновом пределе8060-84.jpg G. в. в магн. металле - одна из ветвей колебания ферми-жидкости.

Газ магнонов. Магноны являются бозонами. При конечной температуре8060-85.jpg магнонов много. Их число NM пропорц. объёму тела V и растёт с ростом Т:
8060-86.jpg

-для ферромагнетиков,
8060-87.jpg

- для антиферромагнетиков.

Мн. свойства магнетиков при8060-88.jpg удобно описывать, считая, что спиновые волны представляют собой почти идеальный газ магнонов (см. Вырожденный газ). Химический потенциал газа магнонов равен 0, т. к. число магнонов не сохраняется; равновесная функция распределения магнонов по энергиям:
8060-89.jpg

Ф-ла (11) позволяет вычислить температурную зависимость термодинамич. характеристик магнетика (намагниченности, теплоёмкости, магн. восприимчивости и др.). Получающиеся выражения тем точнее, чем идеальнее газ магнонов. Неидеальность - результат взаимодействия магнонов друг с другом, с др. квазичастицами ,(с фононами, электронами). С ростом Т число любых квазичастиц растёт, их взаимодействие становится столь существенным, что представление об идеальном газе магнонов перестаёт быть справедливым. Кроме того, может нарушиться условие квазистационарности спиновых волн 8060-90.jpg , где8060-91.jpg- время жизни магнона. Поэтому простейшая концепция газа магнонов применима при8060-92.jpg. При этом важную роль играют низкочастотные (релятивистские) магноны; при8060-93.jpg их значительно больше, чем обменных (последних экспоненциально мало). Однако учёт изменения спектра магнонов при повышении температуры позволяет обобщить концепцию газа магнонов практически на широкий диапазон Т, включающий Тс.

Влияние спиновых волн на кинетические свойства магнетиков. Спиновые волны позволяют описать не только термодинамич. (равновесные) свойства магнетиков, но и их кинетические и резонансные свойства. В теплопроводности магнетиков наряду с фононами и электронами (для проводников) принимают участие магноны: один из механизмов затухания звука - рассеяние звуковых волн на магнонах; в магн. металлах и полупроводниках рассеяние электронов на магнонах - один из механизмов электросопротивления; ферро- и антиферромагнитный резонансы можно представить как превращение фотона в магнон, при ферроакустич. резонансе в магнон превращается фонон.

Для описания кинетических и резонансных процессов существенно время жизни магнона8060-94.jpg. Среди процессов, определяющих время жизни магнонов, выделяют собств. процессы, характерные для идеального кристалла (магнон-магнонные, магнон-фононные и др. взаимодействия), и несобственные (рассеяние магнонов на примесях, дислокациях, границах кристаллитов и поверхности образца).

Взаимодействие магнонов друг с другом и с др. квазичастицами может привести не только к их рассеянию, но и к перестройке их спектра. С возрастанием числа магнонов (NM)наблюдается нелинейный (по NМ) сдвиг частоты спиновой волны. Учёт членов ф-лы (8), «отброшенных» при получении ф-лы (9), приводит к взаимодействию магнонов, носящему характер притяжения. В результате притяжения между магнонами может образоваться своеобразный спиновый комплекс - двухчастичное связанное состояние. В частности, в ферромагнетике, состоящем из атомов со спином 1/2, возникает возбуждение, соответствующее движению по кристаллу двух спинов, связанных между собой и перевёрнутых относительно вектора намагниченности. Как правило, спиновые комплексы образуют магноны с энергией8060-95.jpg [их роль при8060-96.jpgневелика].

Резонанс между спиновой волной и волной колебания др. природы (напр., звуковой) может привести к «расталкиванию» ветвей, что проявляется в существовании гибридных колебаний, напр. магнитоупругих (см. магнитоупругие волны, магнитоупругое взаимодействие).

Экспериментальные методы. Первыми эксперим. методами исследования спиновых волн были измерения температурной зависимости термодинамич. характеристик - намагниченности, магн. части теплоёмкости (рис. 2, 3),
8060-97.jpg

Рис. 2. Температурная зависимость намагниченности ферромагнитного соединения ЕиО. Сплошная кривая - расчёт М(Т) по теории спиновых волн.
8060-98.jpg

Рис. 3. Температурная зависимость магнитной части теплоёмкости См легкоплоскостного антиферромагнетика МnСО3. При низких температурах См = aТ3, резкое отклонение от этого закона происходит при Т > 6 К, соответствующей «включению» второй ветви спектра.

Неупругое рассеяние нейтронов является наиб. информативным методом, позволяющим определить закон дисперсии спиновых волн и оценить время жизни всех типов магнонов. Использование поляризованных нейтронов, кроме того, даёт возможность получить сведения о поляризации спиновых волн. Исследованы спектры сотен магнетиков, в т. ч. сложных (рис. 4, 5).
8060-99.jpg

Рис. 4. Спектр спиновых волн ферромагнитного кобальтового сплава (92% Со, 8% Fе), полученных с помощью неупругого рассеяния нейтронов.

Неупругое рассеяние нейтронов не позволяет исследовать спектр спиновых волн при предельно малых квазиволновых векторах k, т. к. в этом случае пик неупругого рассеяния накладывается на пик упругого рассеяния (см. Магнитная нейтронография ).Ферро- и антиферромагн. резонансы дают возможность измерить значение частот однородной прецессии w0, т. е. щелей8060-100.jpg в спектре магнонов. Для исследования нач. участка спектра (8060-101.jpg см-1) используют резонанс на стоячих спиновых волн в пластинах, параметрич. возбуждение спиновых волн эл--магн. полем, а также неупругое рассеяние света (Мандельштама - Бриллюэна рассеяние). Каждый из методов не универсален, но в совокупности они позволили с большой полнотой определить спектр спиновых волн многих магнитоупорядоченных кристаллов.
8060-102.jpg

Рис. 5. Спектр спиновых волн в кубическом антиферромагнетике RbMnF3, установленный методом неупругого рассеяния нейтронов; кривые - расчёты спектров в предположении, что |J|/k = 3,4 К.

8060-105.jpg

Рис. 6. Теоретический спектр спиновых волн в железоиттриевом гранате.

Длинноволновые участки спектра спиновых волн нек-рых веществ: 1) одноподрешёточный кубич. ферромагнетик (N = 1)
8060-103.jpg

Здесь М - намагниченность насыщения, Nz - размагничивающий фактор, - угол между намагниченностью М и волновым8060-104.jpg вектором k спиновой волны. Коэф. wобм характеризует роль обменного взаимодействия магн. атомов, коэф. wм - магнитодипольного взаимодействия. Ф-ла (12) описывает также акустич. ветвь спиновых волн ферримагнетиков, в частности железоиттриевого граната (ЖИГ), у к-рого 20 подрешёток и соответственно 20 ветвей спиновых волн (рис. 6). В табл. 1 приведены константы акустич. ветви спиновых волн ЖИГ:

Табл. 1.
8060-106.jpg

Табл. 2.
8060-107.jpg

В ЖИГ наиб. исследованы процессы релаксации спиновых волн. В чистых монокристаллах теоретич. значения времён жизни релятивистских магнонов согласуются с экспериментом. При комнатной температуре (300 К)8060-108.jpg 2,6*106 с-1 при8060-109.jpg

2) Двухподрешёточные одноосные антиферромагнетики с магнитной анизотропией типа «лёгкая плоскость» имеют 2 акустич. ветви спиновых волн (Н параллельно лёгкой плоскости):
8060-110.jpg

Здесь НА, НЕ - поля анизотропии и обмена, НД - т. н. поле Дзялошинского, описывающее силу, приводящую к слабому ферромагнетизму,8060-111.jpg - слагаемое, определяемое слабыми взаимодействиями (сверхтонким, магнитоупругим),8060-112.jpg - константы неоднородного обмена (8060-113.jpg - вдоль оси симметрии кристалла,8060-114.jpg- перпендикулярно к оси; табл. 2).

3) Двухподрешёточные антиферромагнетики с магн. анизотропией типа «лёгкая ось» имеют 2 акустич. ветви спиновых волн, вырожденных при Н = 0:
8060-115.jpg

(Н параллельно «лёгкой оси»). Величина щели при8060-116.jpg для большинства исследованных легкоосных антиферромагнетиков лежит в диапазоне 100 - 1000 ГГц.

Спиновые волны в низкоразмерных системах, в кристаллах с большой энергией магнитной анизотропии, в поликристаллах. В двумерных и одномерных системах, описываемых моделью Гейзенберга, Спиновые волн нельзя трактовать как малое колебание, т. к. даже при Т = Тс магн. упорядочение не наступает (в согласии с Мёрмина - Вагнера теоремой). В подобных магнетиках при Т - Тс возникают бесщелевые возбуждения - спиновые волн, у к-рых скорость (если8060-117.jpg ) или эфф. масса (если8060-118.jpg) служит осн. характеристикой, отличающей низкотемпературную фазу (Т < Тс)от высокотемпературной (Т > Тс).

В нек-рых кристаллах (напр., CsCoCl3, FeF2) энергия магн. анизотропии не мала по сравнению с обменной энергией. При этом структура осн. состояния и спектр спиновых волн зависят от конкретного соотношения между обменной энергией и энергией анизотропии. Характерная особенность - сложная зависимость магн. характеристик от магн. поля, перестройка осн. состояния под действием магн. поля.

Длинноволновые спиновые волны 8060-119.jpg сохраняют смысл в поликристаллах .Дополнительное (по сравнению с монокристаллами) затухание спиновых волн связано с рассеянием на границах кристаллитов.

Спиновые волны в парамагнитных металлах и газах. В парамагнитных металлах спиновые волны предсказаны В. П. Силиным в 1960, обнаружены экспериментально в 1967. В немагн. металлах спиновые волны - колебания спиновой плотности электронов проводимости, обусловленные обменным взаимодействием между ними. Спиновые волны в немагн. металлах проявляются, напр., в селективной прозрачности металлич. пластин для эл--магн. волн с частотами, близкими частоте электронного парамагн. резонанса.

В классическом (невырожденном) газе частиц, обладающих спинами, наряду с упругими волнами за счёт обменного взаимодействия между атомами могут распространяться своеобразные волны, также называемые спиновыми. Они предсказаны в 1981, обнаружены в атомарном водороде из Не в 1984.

Литература по спиновые волнам

  1. Ахиезер А. И., БарьяхтарВ. Г., Пелетминский С. В., Спиновые волны, М., 1967;
  2. Вонсовский С. В., Магнетизм, М., 1971;
  3. Уайт Р., Квантовая теория магнетизма, пер. с англ., 2 изд., М., 1985;
  4. Косевич А. М., Иванов Б. А., Ковалев А. С., Нелинейные волны намагниченности. Динамические и топологические солитоны, К., 1983;
  5. Львов В. С., Нелинейные спиновые волны, М., 1987;

см. также лит. к ст. Ферромагнетизм, Антиферромагнетизм.

Л. А. Прозорова

к библиотеке   к оглавлению   FAQ по эфирной физике   ТОЭЭ   ТЭЦ   ТПОИ   ТИ  

Знаете ли Вы, что такое "усталость света"?
Усталость света, анг. tired light - это явление потери энергии квантом электромагнитного излучения при прохождении космических расстояний, то же самое, что эффект красного смещения спектра далеких галактик, обнаруженный Эдвином Хабблом в 1926 г.
На самом деле кванты света, проходя миллиарды световых лет, отдают свою энергию эфиру, "пустому пространству", так как он является реальной физической средой - носителем электромагнитных колебаний с ненулевой вязкостью или трением, и, следовательно, колебания в этой среде должны затухать с расходом энергии на трение. Трение это чрезвычайно мало, а потому эффект "старения света" или "красное смещение Хаббла" обнаруживается лишь на межгалактических расстояниях.
Таким образом, свет далеких звезд не суммируется со светом ближних. Далекие звезды становятся красными, а совсем далекие уходят в радиодиапазон и перестают быть видимыми вообще. Это реально наблюдаемое явление астрономии глубокого космоса. Подробнее читайте в FAQ по эфирной физике.

НОВОСТИ ФОРУМАФорум Рыцари теории эфира
Рыцари теории эфира
 24.07.2019 - 05:16: НОВЫЕ ТЕХНОЛОГИИ - New Technologies -> ПРОБЛЕМА ИСКУССТВЕННОГО ИНТЕЛЛЕКТА - Карим_Хайдаров.
22.07.2019 - 18:29: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Вячеслава Осиевского - Карим_Хайдаров.
21.07.2019 - 14:44: ЭКСПЕРИМЕНТАЛЬНАЯ ФИЗИКА - Experimental Physics -> ИНЕРЦИОИДНЫЕ ГЕНЕРАТОРЫ ЭНЕРГИИ - Карим_Хайдаров.
21.07.2019 - 14:42: ГЕОФИЗИКА И ФИЗИКА ПЛАНЕТ - Geophysics and planetology -> Сезонные колебания уровня вод морей и океанов - Карим_Хайдаров.
21.07.2019 - 09:44: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> ЗА НАМИ БЛЮДЯТ - Карим_Хайдаров.
21.07.2019 - 09:43: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Галины Царёвой - Карим_Хайдаров.
20.07.2019 - 05:34: ЭКОНОМИКА И ФИНАНСЫ - Economy and Finances -> КОЛЛАПС МИРОВОЙ ФИНАНСОВОЙ СИСТЕМЫ - Карим_Хайдаров.
20.07.2019 - 05:30: ЭКОЛОГИЯ - Ecology -> Биологическая безопасность населения - Карим_Хайдаров.
16.07.2019 - 10:00: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> РАСЧЕЛОВЕЧИВАНИЕ ЧЕЛОВЕКА. КОМУ ЭТО НАДО? - Карим_Хайдаров.
16.07.2019 - 09:58: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от О.Н. Четвериковой - Карим_Хайдаров.
12.07.2019 - 17:46: ФИЗИКА ЭФИРА - Aether Physics -> Понятие времени и эфир - Владимир_Афонин.
11.07.2019 - 07:14: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> Проблема государственного терроризма - Карим_Хайдаров.
Bourabai Research Institution home page

Bourabai Research - Технологии XXI века Bourabai Research Institution