Все выпуски
- 2024 Том 16
- Номер 1 (специальный выпуск)
- 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
-
Смешанный алгоритм расчета динамики переноса заряда в ДНК на больших временных интервалах
Компьютерные исследования и моделирование, 2010, т. 2, № 1, с. 63-72Перенос заряда в ДНК моделируется с помощью дискретной модели Холстейна «квантовая частица + классическая цепочка сайтов + взаимодействие». Влияние температуры термостата учитывается с помощью случайной силы, действующей на классические сайты (уравнение Ланжевена). Таким образом, динамика распространения заряда вдоль цепочки описывается системой ОДУ со случайной правой частью. Для интегрирования таких систем обычно применяют алгоритмы 1 или 2 порядка. Мы разработали смешанный алгоритм, имеющий 4 порядок точности по быстрым «квантовым» переменным (заметим, что в «квантовой» подсистеме должно соблюдаться условие: «сумма вероятностей нахождения заряда на сайте постоянна по времени») и 2 порядок по медленным «классическим» переменным, на которые действует случайная сила. Алгоритм позволяет считать на бóльших временах, чем стандартные. В качестве примера приведен модельный расчет развала полярона в однородной цепочке под действием температурных флуктуаций.
Ключевые слова: ДНК, модель Холстейна, уравнение Ланжевена, алгоритм интегрирования ОДУ со случайной правой частью.
Mixed algorithm for modeling of charge transfer in DNA on long time intervals
Computer Research and Modeling, 2010, v. 2, no. 1, pp. 63-72Просмотров за год: 2. Цитирований: 2 (РИНЦ).Charge transfer in DNA is simulated by a discrete Holstein model «quantum particle + classical site chain + interaction». Thermostat temperature is taken into account as stochastic force, which acts on classical sites (Langevin equation). Thus dynamics of charge migration along the chain is described by ODE system with stochastic right-hand side. To integrate the system numerically, algorithms of order 1 or 2 are usually applied. We developed «mixed» algorithm having 4th order of accuracy for fast «quantum» variables (note that in quantum subsystem the condition «sum of probabilities of charge being on site is time-constant» must be held), and 2nd order for slow classical variables, which are affecting by stochastic force. The algorithm allows us to calculate trajectories on longer time intervals as compared to standard algorithms. Model calculations of polaron disruption in homogeneous chain caused by temperature fluctuations are given as an example.
-
Эффективные генераторы псевдослучайных чисел при молекулярном моделировании на видеокартах
Компьютерные исследования и моделирование, 2011, т. 3, № 3, с. 287-308Динамика Ланжевена, метод Монте-Карло и моделирование молекулярной динамики в неявном растворителе требуют больших массивов случайных чисел на каждом шаге расчета. Мы исследовали два подхода в реализации генераторов на графических процессорах. Первый реализует последовательный алгоритм генератора на каждом потоке в отдельности. Второй основан на возможности взаимодействия между потоками и реализует общий алгоритм на всех потоках в целом. Мы покажем использование этих подходов на примере алгоритмов Ran 2, Hybrid Taus и Lagged Fibonacci. Для проверки случайности полученных чисел мы использовали разработанные генераторы при моделировании динамики Ланжевена N независимых гармонических осцилляторов в термостате. Это позволило нам оценить статистические характеристики генераторов. Мы также исследовали производительность, использование памяти и ускорение, получаемое при переносе алгоритма с центрального на графический процессор.
Ключевые слова: псевдослучайные числа, графический процессор, генератор, молекулярное моделирование.
Efficient Pseudorandom number generators for biomolecular simulations on graphics processors
Computer Research and Modeling, 2011, v. 3, no. 3, pp. 287-308Просмотров за год: 11. Цитирований: 2 (РИНЦ).Langevin Dynamics, Monte Carlo, and all-atom Molecular Dynamics simulations in implicit solvent require a reliable source of pseudorandom numbers generated at each step of calculation. We present the two main approaches for implementation of pseudorandom number generators on a GPU. In the first approach, inherent in CPU-based calculations, one PRNG produces a stream of pseudorandom numbers in each thread of execution, whereas the second approach builds on the ability of different threads to communicate, thus, sharing random seeds across the entire device. We exemplify the use of these approaches through the development of Ran2, Hybrid Taus, and Lagged Fibonacci algorithms. As an application-based test of randomness, we carry out LD simulations of N independent harmonic oscillators coupled to a stochastic thermostat. This model allows us to assess statistical quality of pseudorandom numbers. We also profile performance of these generators in terms of the computational time, memory usage, and the speedup factor (CPU/GPU time).
-
Движение открытых состояний ДНК под действием случайной силы
Компьютерные исследования и моделирование, 2015, т. 7, № 6, с. 1295-1307Известно, что в нативном состоянии молекула ДНК всегда содержит некоторое количество локально расплетенных участков, часто называемых открытыми состояниями ДНК. Считается, что эти состояния играют важную роль в ДНК-белковом узнавании, и изучение их динамики может пролить дополнительный свет на механизмы регуляции транскрипции и репликации. В этой статье мы рассматриваем влияние термостата на движение открытых состояний в искусственной последовательности, состоящей из четырех однородных областей. Мы построим энергетический профиль этой последовательности и исследуем траектории движения открытых состояний в этом профиле под действием случайной силы.
Motion of DNA open states influenced by random force
Computer Research and Modeling, 2015, v. 7, no. 6, pp. 1295-1307Просмотров за год: 3.It is known that in the native state the DNA molecule always contains some amount of locally unwound regions, often called the open states of DNA. It is believed that these states play an important role in DNA-protein recognition and that the study of the open states dynamics may shed further light on the mechanisms of regulation of transcription and replication. In this paper we consider the effect of the thermostat on the movement of the open states in the artificial sequence consisting of four homogeneous regions. We construct the energetic profile of the sequence and investigate the trajectories of the movement of the open states under the action of a random force.
-
Трехмерное молекулярно-динамическое моделирование термодинамического равновесия нагретого никеля
Компьютерные исследования и моделирование, 2015, т. 7, № 3, с. 573-579Представленная работа посвящена молекулярно-динамическому моделированию процессов термического воздействия на металлический образец, который состоит из атомов никеля. Для решения этой задачи используется континуальная математическая модель, основанная на уравнениях классической механики Ньютона, выбран численный метод, использующий в основе схему Верле, предложен параллельный алго- ритм и осуществлена его реализация в рамках MPIи OpenMP. С помощью разработанной параллельной программы было проведено исследование термодинамического равновесия атомов никеля при условии нагрева образца до желаемой температуры. В численных экспериментах определены оптимальные параметры методики расчета и физические параметры исследуемого процесса. Полученные численные результаты хорошо согласуются с известными теоретическими и экспериментальными данными.
Ключевые слова: молекулярно-динамическое моделирование, никель, ППА, температура, термостат, уравнения Ньютона, параллельные алгоритмы и программы, MPI, OpenMP.
3D molecular dynamic simulation of thermodynamic equilibrium problem for heated nickel
Computer Research and Modeling, 2015, v. 7, no. 3, pp. 573-579Просмотров за год: 2.This work is devoted to molecular dynamic modeling of the thermal impact processes on the metal sample consisting of nickel atoms. For the solution of this problem, a continuous mathematical model on the basis of the classical Newton mechanics equations has been used; a numerical method based on the Verlet scheme has been chosen; a parallel algorithm has been offered, and its realization within the MPI and OpenMP technologies has been executed. By means of the developed parallel program, the investigation of thermodynamic equilibrium of nickel atoms’ system under the conditions of heating a sample to desired temperature has been executed. In numerical experiments both optimum parameters of calculation procedure and physical parameters of analyzed process have been defined. The obtained numerical results are well corresponding to known theoretical and experimental data.
Журнал индексируется в Scopus
Полнотекстовая версия журнала доступна также на сайте научной электронной библиотеки eLIBRARY.RU
Журнал входит в систему Российского индекса научного цитирования.
Журнал включен в базу данных Russian Science Citation Index (RSCI) на платформе Web of Science
Международная Междисциплинарная Конференция "Математика. Компьютер. Образование"