Все выпуски
- 2024 Том 16
- 2023 Том 15
- 2022 Том 14
- 2021 Том 13
- 2020 Том 12
- 2019 Том 11
- 2018 Том 10
- 2017 Том 9
- 2016 Том 8
- 2015 Том 7
- 2014 Том 6
- 2013 Том 5
- 2012 Том 4
- 2011 Том 3
- 2010 Том 2
- 2009 Том 1
-
Представление групп автоморфизмами нормальных топологических пространств
Компьютерные исследования и моделирование, 2009, т. 1, № 3, с. 243-249Доказывается, что произвольная алгебраическая группа алгебраически изоморфна полной группе автоморфизмов некоторого топологического пространства (автобиекций, сохраняющих открытые множества) с нормальным типом отделимости (Т4 + Т1). Кроме того, любое непрерывное действие группы на нормальном топологическом пространстве может быть получено как действие полной группы автоморфизмов нормального топологического пространства на его подпространстве.
Representation of groups by automorphisms of normal topological spaces
Computer Research and Modeling, 2009, v. 1, no. 3, pp. 243-249Просмотров за год: 1.The famous fact [3, 5] of existence of an exact representation for any finite group in the form of the full automorphism group of a finite graph was generalize in [4]. For an arbitrary group exact representation exists in the form of the full automorphism group of Kolmogorov topological space (weak type of separability T0). For a finite group a finite space may be chosen, thus allowing to restore a finite graph with the same number of vertices and having the same automorphism group. Such topological spaces and graphs are called topological imprints and graph imprints of a group (T-imprints and G-imprints, respectively). The question of maximum type of separability of a topological space for which T-imprint can be obtained for any group is open. The author proves that the problem can be solved for the class of normal topology (maximal type of separability T4+T0). Special finite T-imprint for a symmetric group may be obtained as a discrete topology; for any other group minimal cardinality of normal T-imprint is countable. There is a generic procedure to construct a T-imprint for any group. For a finite group this procedure allows finite space partitioning into subspaces having G-imprint of the original group as their connectivity graphs.
-
Полиполярная координация и симметрии
Компьютерные исследования и моделирование, 2010, т. 2, № 4, с. 329-341Полиполярная система координат формируется семейством параметризованных по радиусу изофокусных kf-лемнискат. Как и классическая полярная система координат, она характеризует точку плоскости полиполярным радиусом ρ и полиполярным углом φ. Для любой связности семейство изометрических кривых ρ = const – лемнискат и семейство градиентных кривых φ = const являются взаимно ортогональными сопряженными координатными семействами. Рассмотрены особенности полиполярной координации, ее симметрии, а также криволинейные симметрии на многофокусных лемнискатах.
Ключевые слова: кривые, фокусы, многофокусные лемнискаты, овалы Кассини, полярная система координат, координатные семейства, группы симметрий, криволинейные симметрии.
Polypolar coordination and symmetries
Computer Research and Modeling, 2010, v. 2, no. 4, pp. 329-341Просмотров за год: 1.The polypolar system of coordinates is formed by a family of a parametrized on a radius isofocal of kf-lemniscates. As well as the classical polar system of coordinates, it characterizes a point of a plane by a polypolar radius ρ and polypolar angle φ. For anyone connectedness a family isometric of curve ρ = const – lemniscates and family gradient of curves φ = const – are mutually orthogonal conjugate coordinate families. The singularities of polypolar coordination, its symmetry, and also curvilinear symmetries on multifocal lemniscates are considered.
-
Континуальные трансформирующиеся оболочки из тонких пластин
Компьютерные исследования и моделирование, 2011, т. 3, № 1, с. 3-29Изучаются трансформирующиеся системы, собранные из трапециевидных пластин. При развертывании пакета пластинок образуется сетчатая оболочка с шестигранными ячейками. Доказывается, что при определенных соотношениях размеров граней в шестизвенниках появляются дополнительные внутренние степени свободы. Если же используются тонкие пластинки, то континуальная аппроксимация развернутой сети может интерпретироваться как оболочка с широким набором локальных кривизн. Строится кинематика континуальной модели методом подвижного репера Картана. Изучается механическое поведение континуальных сетей, если цилиндрические шарниры между пластинами выполнены из пластических материалов, обладающих памятью формы. Исследуются переходы оболочек из одной равновесной формы в другую. Показаны возможные практические применения континуальных сетей.
Ключевые слова: континуальные трансформирующиеся оболочки, репер Картана, пластические материалы, память формы, устойчивое равновесие, практическое применение.
Continuum deployable shells made of thin plates
Computer Research and Modeling, 2011, v. 3, no. 1, pp. 3-29Цитирований: 3 (РИНЦ).This paper covers deployable systems assembled from trapezium plates. When the plate package is unwrapped, a net shell with six loop cells is formed. It is proved that additional degrees of freedom appear in case of certain correlation between the sizes of the six loop faces. When thin plates were used, the continuum approximation of the deployed net could be interpreted as a shell with a wide variety of local curvatures. Kinematics of the continuum model is analyzed by the method of Cartan moving hedron. Mechanical behavior of continuum nets is studied when cylindrical hinges between the plates are completed of shape memory plastic materials. The paper researches into shell transformations from one stable form to the other. Various practical applications of the continuum nets are demonstrated.
-
Квазиклассическое приближение для многомерного нелокального уравнения Фишера–Колмогорова–Петровского–Пискунова
Компьютерные исследования и моделирование, 2015, т. 7, № 2, с. 205-219Для многомерного нелокального уравнения Фишера–Колмогорова–Петровского–Пискунова в классе траекторно-сосредоточенных функций построены квазиклассические асимптотики с точностью $O(D^{N/2})$, $N\geqslant3$. С помощью операторов симметрии получен счетный набор асимптотических решений исходного уравнения с точностью $O(D^{3/2})$. В явном виде построены асимптотические решения двумерного уравнения Фишера–Колмогорова–Петровского–Пискунова.
Ключевые слова: нелокальное уравнение Фишера–Колмогорова–Петровского–Пискунова, асимптотическое решение, система Эйнштейна–Эренфеста.
Semiclassical approximation for the nonlocal multidimensional Fisher–Kolmogorov–Petrovskii–Piskunov equation
Computer Research and Modeling, 2015, v. 7, no. 2, pp. 205-219Просмотров за год: 4.Semiclassical asymptotic solutions with accuracy $O(D^{N/2})$, $N\geqslant3$ are constructed for the multidimensional Fisher–Kolmogorov–Petrovskii–Piskunov equation in the class of trajectory-concentrated functions. Using the symmetry operators a countable set of asymptotic solutions with accuracy $O(D^{3/2})$ is obtained. Asymptotic solutions of two-dimensional Fisher–Kolmogorov–Petrovskii–Piskunov equation are found in explicit
form. -
Глобальный бифуркационный анализ рациональной системы Холлинга
Компьютерные исследования и моделирование, 2017, т. 9, № 4, с. 537-545В статье рассматривается квартичное семейство планарных векторных полей, соответствующее рациональной системе Холлинга, которая моделирует динамику популяций типа «хищник–жертва» в данной экологической или биомедицинской системе и которая обобщает классическую систему Лотки–Вольтерры. В простейших математических моделях изменение концентрации жертв в единицу времени в расчете на одного хищника, которое характеризуется так называемой функцией отклика, прямо пропорционально концентрации жертв, т. е. функция отклика в этих моделях линейная. Это означает, что в системе нет насыщения хищников, когда количество жертв достаточно велико. Однако было бы более реалистично рассматривать нелинейные и ограниченные функции отклика, и в литературе действительно используются различные виды таких функций для моделирования отклика хищников. После алгебраических преобразований рациональную систему Холлинга можно записать в виде квартичной динамической системы. Для исследования характера и расположения особых точек в фазовой плоскости этой системы используется разработанный нами метод, смысл которого состоит в том, чтобы получить простейшую (хорошо известную) систему путем обращения в нуль некоторых параметров (обычно параметров, поворачивающих поле) исходной системы, а затем последовательно вводить эти параметры, изучая динамику особых точек (как конечных, так и бесконечно удаленных) в фазовой плоскости. Используя полученную информацию об особых точках и применяя наш геометрический подход к качественному анализу, мы изучаем бифуркации предельных циклов квартичной системы. Чтобы контролировать все бифуркации предельных циклов, особенно бифуркации кратных предельных циклов, необходимо знать свойства и комбинировать действия всех параметров, поворачивающих векторное поле системы. Это может быть сделано с помощью принципа окончания Уинтнера–Перко, согласно которому максимальное однопараметрическое семейство кратных предельных циклов заканчивается либо в особой точке, которая, как правило, имеет ту же кратность (цикличность), либо на сепаратрисном цикле, который также, как правило, имеет ту же кратность (цикличность). Применяя этот принцип, мы доказываем, что квадричная система (и соответствующая рациональная система Холлинга) может иметь не более двух предельных циклов, окружающих одну особую точку.
Ключевые слова: рациональная динамическая система Холлинга, параметр поворота поля, бифуркация, особая точка, предельный цикл, принцип окончания Уинтнера–Перко.
Global bifurcation analysis of a rational Holling system
Computer Research and Modeling, 2017, v. 9, no. 4, pp. 537-545Просмотров за год: 11.In this paper, we consider a quartic family of planar vector fields corresponding to a rational Holling system which models the dynamics of the populations of predators and their prey in a given ecological or biomedical system and which is a variation on the classical Lotka–Volterra system. For the latter system, the change of the prey density per unit of time per predator called the response function is proportional to the prey density. This means that there is no saturation of the predator when the amount of available prey is large. However, it is more realistic to consider a nonlinear and bounded response function, and in fact different response functions have been used in the literature to model the predator response. After algebraic transformations, the rational Holling system can be written in the form of a quartic dynamical system. To investigate the character and distribution of the singular points in the phase plane of the quartic system, we use our method the sense of which is to obtain the simplest (well-known) system by vanishing some parameters (usually field rotation parameters) of the original system and then to input these parameters successively one by one studying the dynamics of the singular points (both finite and infinite) in the phase plane. Using the obtained information on singular points and applying our geometric approach to the qualitative analysis, we study the limit cycle bifurcations of the quartic system. To control all of the limit cycle bifurcations, especially, bifurcations of multiple limit cycles, it is necessary to know the properties and combine the effects of all of the rotation parameters. It can be done by means of the Wintner–Perko termination principle stating that the maximal one-parameter family of multiple limit cycles terminates either at a singular point which is typically of the same multiplicity (cyclicity) or on a separatrix cycle which is also typically of the same multiplicity (cyclicity). Applying this principle, we prove that the quartic system (and the corresponding rational Holling system) can have at most two limit cycles surrounding one singular point.
-
Научные и педагогические школы Александра Сергеевича Холодова
Компьютерные исследования и моделирование, 2018, т. 10, № 5, с. 561-579В развитии науки важную роль играют научные школы — объединения исследователей, связанные общей проблемой, идеями и методами, используемыми для решения проблемы. Научные школы формируются вокруг лидера и объединяющей идеи.
За время научной деятельности академика А. С. Холодова вокруг него сформировалось несколько научных школ. В обзоре делается попытка представить основные научные направления, вокруг которых сформировались яркие коллективы с общими системами взглядов и подходами к исследованиям. В обзоре отмечается эта общая основа. Во-первых, это развитие группы численных методов для решения систем дифференциальных уравнений в частных производных гиперболического типа — сеточно-характеристические методы. Во-вторых, описание численных методов в пространствах неопределенных коэф- фициентов. Этот подход развивался как для всех типов уравнений в частных производных, так и для обыкновенных дифференциальных уравнений.
На основе предложенных А. С. Холодовым численных подходов сложились научные коллективы, работающие в разных предметных областях. Это математическое моделирование динамики плазмы, динамики деформируемого твердого тела, некоторых задач биологии, биофизики, медицинской физики и биомеханики. Сравнительно новые направления — решение задач на графах (процессы транспортировки электроэнергии, моделирование транспортных потоков на дорожной сети и т. д.).
В обзоре делается попытка отследить деятельность научных школ от момента их зарождения до настоящего времени, проследить связь работ А. С. Холодова с работами его учеников и коллег. Полный обзор деятельности всех научных школ, сформировавшихся вокруг Александра Сергеевча, невозможен ввиду огромного количества и разнообразия научных результатов.
Делается также попытка связать деятельность научных школ с появлением научно-образовательной школы в Московском физико-техническом институте.
Ключевые слова: научная школа, сеточно-характеристические методы, пространства неопределенных коэффициентов, динамика плазмы, динамика деформируемого твердого тела, биомеханика, процессы на графах.
Scientific and pedagogical schools founded by A. S. Kholodov
Computer Research and Modeling, 2018, v. 10, no. 5, pp. 561-579Просмотров за год: 42.In the science development an important role the scientific schools are played. This schools are the associations of researchers connected by the common problem, the ideas and the methods used for problems solution. Usually Scientific schools are formed around the leader and the uniting idea.
The several sciences schools were created around academician A. S. Kholodov during his scientific and pedagogical activity.
This review tries to present the main scientific directions in which the bright science collectives with the common frames of reference and approaches to researches were created. In the review this common base is marked out. First, this is development of the group of numerical methods for hyperbolic type systems of partial derivatives differential equations solution — grid and characteristic methods. Secondly, the description of different numerical methods in the undetermined coefficients spaces. This approach developed for all types of partial equations and for ordinary differential equations.
On the basis of A. S. Kholodov’s numerical approaches the research teams working in different subject domains are formed. The fields of interests are including mathematical modeling of the plasma dynamics, deformable solid body dynamics, some problems of biology, biophysics, medical physics and biomechanics. The new field of interest includes solving problem on graphs (such as processes of the electric power transportation, modeling of the traffic flows on a road network etc).
There is the attempt in the present review analyzed the activity of scientific schools from the moment of their origin so far, to trace the connection of A. S. Kholodov’s works with his colleagues and followers works. The complete overview of all the scientific schools created around A. S. Kholodov is impossible due to the huge amount and a variety of the scientific results.
The attempt to connect scientific schools activity with the advent of scientific and educational school in Moscow Institute of Physics and Technology also becomes.
-
О проектировании нуля на линейное многообразие, многогранник и вершину многогранника. Ньютоновские методы минимизации
Компьютерные исследования и моделирование, 2019, т. 11, № 4, с. 563-591Рассматривается подход к построению методов решения задачи квадратичного программирования для расчета направления спуска в ньютоновских методах минимизации гладкой функции на множестве, заданном набором линейных равенств. Подход состоит из двух этапов.
На первом этапе задача квадратичного программирования преобразуется численно устойчивым прямым мультипликативным алгоритмом в эквивалентную задачу о проектировании начала координат на линейное многообразие, что определяет новую математическую формулировку двойственной квадратичной задачи. Для этого предложен численно устойчивый прямой мультипликативный метод решения систем линейных уравнений, учитывающий разреженность матриц, представленных в упакованном виде. Преимущество подхода состоит в расчете модифицированных факторов Холесского для построения существенно положительно определенной матрицы системы уравнений и ее решения в рамках одной процедуры, а также в возможности минимизации заполнения главных строк мультипликаторов без потери точности результатов. Причем изменения в позиции очередной обрабатываемой строки матрицы не вносятся, что позволяет использовать статические форматы хранения данных.
На втором этапе необходимые и достаточные условия оптимальности в форме Куна–Таккера определяют расчет направления спуска — решение двойственной квадратичной задачи сводится к решению системы линейных уравнений с симметричной положительно определенной матрицей коэффициентов для расчета множителей Лагранжа и к подстановке решения в формулу для расчета направления спуска.
Доказано, что предложенный подход к расчету направления спуска численно устойчивыми прямыми мультипликативными методами на одной итерации требует по кубическому закону меньше вычислений, чем одна итерация по сравнению с известным двойственным методом Гилла и Мюррея. Кроме того, предложенный метод допускает организацию вычислительного процесса с любой начальной точки, которую пользователь выберет в качестве исходного приближения решения.
Представлены варианты постановки задачи о проектировании начала координат на линейное многообразие, выпуклый многогранник и вершину выпуклого многогранника. Также описаны взаимосвязь и реализация методов решения этих задач.
Ключевые слова: ньютоновские методы, квадратичное программирование, двойственная квадратичная задача, разреженные матрицы, факторизация Холесского, прямой мультипликативный алгоритм, численная устойчивость, задача о проектировании нуля, линейное многообразие, вершина многогранника.
Designing a zero on a linear manifold, a polyhedron, and a vertex of a polyhedron. Newton methods of minimization
Computer Research and Modeling, 2019, v. 11, no. 4, pp. 563-591Просмотров за год: 6.We consider the approaches to the construction of methods for solving four-dimensional programming problems for calculating directions for multiple minimizations of smooth functions on a set of a given set of linear equalities. The approach consists of two stages.
At the first stage, the problem of quadratic programming is transformed by a numerically stable direct multiplicative algorithm into an equivalent problem of designing the origin of coordinates on a linear manifold, which defines a new mathematical formulation of the dual quadratic problem. For this, a numerically stable direct multiplicative method for solving systems of linear equations is proposed, taking into account the sparsity of matrices presented in packaged form. The advantage of this approach is to calculate the modified Cholesky factors to construct a substantially positive definite matrix of the system of equations and its solution in the framework of one procedure. And also in the possibility of minimizing the filling of the main rows of multipliers without losing the accuracy of the results, and no changes are made in the position of the next processed row of the matrix, which allows the use of static data storage formats.
At the second stage, the necessary and sufficient optimality conditions in the form of Kuhn–Tucker determine the calculation of the direction of descent — the solution of the dual quadratic problem is reduced to solving a system of linear equations with symmetric positive definite matrix for calculating of Lagrange's coefficients multipliers and to substituting the solution into the formula for calculating the direction of descent.
It is proved that the proposed approach to the calculation of the direction of descent by numerically stable direct multiplicative methods at one iteration requires a cubic law less computation than one iteration compared to the well-known dual method of Gill and Murray. Besides, the proposed method allows the organization of the computational process from any starting point that the user chooses as the initial approximation of the solution.
Variants of the problem of designing the origin of coordinates on a linear manifold, a convex polyhedron and a vertex of a convex polyhedron are presented. Also the relationship and implementation of methods for solving these problems are described.
-
Квазиклассические решения уравнения Гросса–Питаевского, локализованные в окрестности окружности
Компьютерные исследования и моделирование, 2009, т. 1, № 4, с. 359-365В квазиклассическом приближении показано, что для конденсата Бозе–Эйнштейна, моделируемого уравнением Гросса–Питаевского с притягивающей нелинейностью при специальной конфигурации внешнего поля магнитной ловушки, возможны неколлапсирующие солитоноподобные волновые функции.
Ключевые слова: уравнение Гросса–Питаевского, бозе-эйнштейновский конденсат, магнитооптическая ловушка.
Semiclassical solutions localized in a neighborhood of a circle for the Gross–Pitaevskii equation
Computer Research and Modeling, 2009, v. 1, no. 4, pp. 359-365Цитирований: 1 (РИНЦ).Non-collapsing soliton-like wave functions are shown to exist in semiclassical approximation for the Bose-Einstein condensate model based on the Gross–Pitaevskii equation with attractive nonlinearity and external field of magnetic trap of special form.
-
Интегрирование релятивистских волновых уравнений в космологической модели Бъянки IX
Компьютерные исследования и моделирование, 2016, т. 8, № 3, с. 433-443В работе рассматривается интегрирование уравнений Клейна–Гордона и Дирака в космологической модели Бъянки IX. При помощи метода некоммутативного интегрирования дифференциальных уравнений найдены новые точные решения для осесимметричной модели.
Метод некоммутативного интегрирования в данной задаче основан на использовании специального бесконечномерного голоморфного представления группы вращений, которое строится по невырожденной орбите коприсоединенного представления и комплексной поляризации невырожденного ковектора. Матричные элементы данного представления образуют полный и ортогональный набор и позволяют ввести обобщенное преобразование Фурье. Оператор Казимира группы вращений при этом преобразовании переходит в константу, а операторы симметрии, порожденные векторными полями Киллинга, — в линейные дифференциальные операторы первого порядка от одной зависимой переменной. Таким образом, релятивистские волновые уравнения на группе вращений допускают некоммутативную редукцию к обыкновенному дифференциальному уравнению. В отличие от широко известного метода разделения переменных метод некоммутативного интегрирования учитывает неабелеву алгебру операторов симметрии и дает решения, несущие информацию о некоммутативной симметрии задачи. Такие решения могут быть полезны для учета вакуумных квантовых эффектов и расчета конечных функций Грина методом раздвижки точек.
В работе для осесимметричной модели проведено сравнение полученных решений с известными, которые получаются методом разделения переменных. Показано, что некоммутативные решения выражаются через элементарные функции, тогда как известные решения определяются функцией Вигнера. Причем некоммутативно редуцированное уравнение Клейна–Гордона для осесимметричной модели совпадает с уравнением, редуцированным методом разделения переменных. А некоммутативно редуцированное уравнение Дирака эквивалентно редуцированному уравнению, полученному методом разделения переменных.
Ключевые слова: некоммутативное интегрирование, Бъянки IX.
Integration the relativistic wave equations in Bianchi IX cosmology model
Computer Research and Modeling, 2016, v. 8, no. 3, pp. 433-443We consider integration Clein–Gordon and Dirac equations in Bianchi IX cosmology model. Using the noncommutative integration method we found the new exact solutions for Taub universe.
Noncommutative integration method for Bianchi IX model is based on the use of the special infinite-dimensional holomorphic representation of the rotation group, which is based on the nondegenerate orbit adjoint representation, and complex polarization of degenerate covector. The matrix elements of the representation of form a complete and orthogonal set and allow you to use the generalized Fourier transform. Casimir operator for rotation group under this transformation becomes constant. And the symmetry operators generated by the Killing vector fields in the linear differential operators of the first order from one dependent variable. Thus, the relativistic wave equation on the rotation group allow non-commutative reduction to ordinary differential equations. In contrast to the well-known method of separation of variables, noncommutative integration method takes into account the non-Abelian algebra of symmetry operators and provides solutions that carry information about the non-commutative symmetry of the task. Such solutions can be useful for measuring the vacuum quantum effects and the calculation of the Green’s functions by the splitting-point method.
The work for the Taub model compared the solutions obtained with the known, which are obtained by separation of variables. It is shown that the non-commutative solutions are expressed in terms of elementary functions, while the known solutions are defined by the Wigner function. And commutative reduced by the Klein–Gordon equation for Taub model coincides with the equation, reduced by separation of variables. A commutative reduced by the Dirac equation is equivalent to the reduced equation obtained by separation of variables.
Keywords: noncommutative integration, Bianchi IX.Просмотров за год: 5. -
Классификация динамических режимов переключения намагниченности в трехслойной ферромагнитной структуре в зависимости от спин-поляризованного тока инжекции и внешнего магнитного поля. I. Продольная анизотропия
Компьютерные исследования и моделирование, 2016, т. 8, № 4, с. 605-620В приближении однородной намагниченности построена математическая модель ячейки памяти MRAM c осью анизотропии, расположенной в плоскости запоминающего ферромагнитного слоя ячейки и ориентированной параллельно ее краю (продольная анизотропия). Модель базируется на уравнении Ландау–Лифшица–Гильберта с токовым членом в форме Слончевского–Берже. Выведена система обыкновенных дифференциальных уравнений в нормальном виде, описывающая динамику намагниченности в трехслойной вентильной структуре Co/Cu/Co в зависимости от величины тока инжекции и внешнего магнитного поля, параллельного оси анизотропии магнитных слоев. Показано, что при любых токах и полях система имеет два основных состояния равновесия, расположенных на оси, совпадающей с осью анизотропии. Проведен анализ устойчивости этих состояний равновесия. Выписаны уравнения для определения дополнительных состояний равновесия. Показано, что в зависимости от величины внешнего магнитного поля и тока инжекции система может иметь всего два, четыре и шесть симметричных относительно оси анизотропии положений равновесия. Построены бифуркационные диаграммы, характеризующие основные типы динамики вектора намагниченности свободного слоя. Проведена классификация фазовых портретов на единичной сфере в зависимости от управляющих параметров (тока и поля). Изучены особенности динамики вектора намагниченности в каждой из характерных областей бифуркационной диаграммы и численно построены траектории переключения. Для построения траекторий использовался метод Рунге–Кутты. Найдены параметры, при которых существуют неустойчивые и устойчивые предельные циклы. Установлено, что неустойчивые предельные циклы существуют вокруг основного устойчивого равновесия на оси, совпадающей с осью анизотропии, а устойчивые циклы — вокруг неустойчивых дополнительных равновесий. Граница области существования устойчивых предельных циклов рассчитана численно. Обнаружены новые типы динамики под влиянием внешнего магнитного поля и спин-поляризованного тока инжекции: случайное и неполное переключение намагниченности. Аналитически определены значения пороговых токов переключения в зависимости от внешнего магнитного поля. Численно выполнены оценки времени переключения в зависимости от величин управляющих параметров.
Ключевые слова: память MRAM, одноосная анизотропия, намагниченность, свободный слой, закрепленный слой, уравнение Ландау–Лифшица–Гильберта, переключение намагниченности.
Classification of dynamical switching regimes in a three-layered ferromagnetic nanopillar governed by spin-polarized injection current and external magnetic field. I. Longitudinal anisotropy
Computer Research and Modeling, 2016, v. 8, no. 4, pp. 605-620Просмотров за год: 2. Цитирований: 6 (РИНЦ).The mathematical model of the magnetic memory cell MRAM with the in-plane anisotropy axis parallel to the edge of a free ferromagnetic layer (longitudinal anisotropy) has been constructed using approximation of uniform magnetization. The model is based on the Landau–Lifshits–Gilbert equation with the injection-current term in the Sloncžewski–Berger form. The set of ordinary differential equations for magnetization dynamics in a three-layered Co/Cu/Cu valve under the control of external magnetic field and spin-polarized current has been derived in the normal coordinate form. It was shown that the set of equations has two main stationary points on the anisotropy axis at any values of field and current. The stationary analysis of them has been performed. The algebraic equations for determination of additional stationary points have been derived. It has been shown that, depending on the field and current magnitude, the set of equations can have altogether two, four, or six stationary points symmetric in pairs relatively the anisotropy axis. The bifurcation diagrams for all the points have been constructed. The classification of the corresponding phase portraits has been performed. The typical trajectories were calculated numerically using Runge–Kutta method. The regions, where stable and unstable limit cycles exist, have been determined. It was found that the unstable limit cycles exist around the main stable equilibrium point on the axis that coincides with the anisotropy one, whereas the stable cycles surround the unstable additional points of equilibrium. The area of their existence was determined numerically. The new types of dynamics, such as accidental switching and non-complete switching, have been found. The threshold values of switching current and field have been obtained analytically. The estimations of switching times have been performed numerically.
Журнал индексируется в Scopus
Полнотекстовая версия журнала доступна также на сайте научной электронной библиотеки eLIBRARY.RU
Журнал входит в систему Российского индекса научного цитирования.
Журнал включен в базу данных Russian Science Citation Index (RSCI) на платформе Web of Science
Международная Междисциплинарная Конференция "Математика. Компьютер. Образование"