اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدTowards investigation of evolution of dynamical systems with independence of time accuracy: more classes of systems
30
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
The recently developed method (Paper 1) enabling one to investigate the evolution of dynamical systems with an accuracy not dependent on time is developed further. The classes of dynamical systems which can be studied by that method are much extended, now including systems that are; (1) non-Hamiltonian, conservative; (2) Hamiltonian with time-dependent perturbation; (3) non-conservative (with dissipation). These systems cover various types of N-body gravitating systems of astrophysical and cosmological interest, such as the orbital evolution of planets, minor planets, artificial satellites due to tidal, non-tidal perturbations and thermal thrust, evolving close binary stellar systems, and the dynamics of accretion disks.
قيم البحث
اقرأ أيضاً
Extensive N-body simulations are among the key means for the study of numerous astrophysical and cosmological phenomena, so various schemes are developed for possibly higher accuracy computations. We demonstrate the principal possibility for revealin
g the evolution of a perturbed Hamiltonian system with an accuracy independent on time. The method is based on the Laplace transform and the derivation and analytical solution of an evolution equation in the phase space for the resolvent and using computer algebra.
We describe AMUSE, the Astrophysical Multipurpose Software Environment, a programming framework designed to manage multi-scale, multi-physics simulations in a hierarchical, extensible, and internally consistent way. Constructed as a collection of ind
ividual modules, AMUSE allows computational tools for different physical domains to be easily combined into a single task. It facilitates the coupling of modules written in different languages by providing inter-language tools and a standard programming interface that represents a balance between generality and computational efficiency. The framework currently incorporates the domains of stellar dynamics, stellar evolution, gas dynamics, and radiative transfer. We present some applications of the framework and outline plans for future development of the package.
Asteroid material is detected in white dwarfs (WDs) as atmospheric pollution by metals, in the form of gas/dust discs, or in photometric transits. Within the current paradigm, minor bodies need to be scattered, most likely by planets, into highly ecc
entric orbits where the material gets disrupted by tidal forces and then accreted onto the star. This can occur through a planet-planet scattering process triggered by the stellar mass loss during the post main-sequence evolution of planetary systems. So far, studies of the $N$-body dynamics of this process have used artificial planetary system architectures built ad hoc. In this work, we attempt to go a step further and study the dynamical instability provided by more restrictive systems, that, at the same time allow us an exploration of a wider parameter space: the hundreds of multiple planetary systems found around main-sequence (MS) stars. We find that most of our simulated systems remain stable during the MS, Red and Asymptotic Giant Branch and for several Gyr into the WD phases of the host star. Overall, only $approx$ 2.3$%$ of the simulated systems lose a planet on the WD as a result of dynamical instability. If the instabilities take place during the WD phase most of them result in planet ejections with just 5 planetary configurations ending as a collision of a planet with the WD. Finally 3.2$%$ of the simulated systems experience some form of orbital scattering or orbit crossing that could contribute to the pollution at a sustained rate if planetesimals are present in the same system.
We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of sph
erical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures, and the introduction of a parallel random number generation scheme, as well as a parallel sorting algorithm, required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. The implementation uses the Message Passing Interface (MPI) library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude, from 10^5 to 10^7. We find that our results are in good agreement with self-similar core-collapse solutions, and the core collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within less than 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N=10^5, 128 for N=10^6 and 256 for N=10^7. The runtime reaches a saturation with the addition of more processors beyond these limits which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60x, 100x, and 220x, respectively.
In this thesis we study the strongly-correlated-electron physics of the longstanding H-Tc-superconductivity problem using a non-perturbative method, the Dynamical Mean Field Theory (DMFT), capable to go beyond standard perturbation-theory techniques.
DMFT is by construction a local theory which neglects spatial correlation. Experiments have however shown that the latter is a fundamental property of cuprate materials. In a first step, we approach the problem of spatial correlation in the normal state of cuprate materials using a phenomenological Fermi-Liquid-Boltzmann model. We then introduce and develop in detail an extension to DMFT, the Cellular Dynamical Mean Field Theory (CDMFT), capable of considering short-ranged spatial correlation in a system, and we implement it using the exact diagonalization algorithm . After benchmarking CDMFT with the exact one-dimensional solution of the Hubbard Model, we employ it to study the density-driven Mott metal-insulator transition in the two-dimensional Hubbard Model, focusing in particular on the anomalous properties of the doped normal state close to the Mott insulator. We finally study the superconducting state. We show that within CDMFT the one-band Hubbard Model supports a d-wave superconductive state, which strongly departs from the standard BCS theory. We conjecture a link between the instabilities found in the normal state and the onset of superconductivity.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
