اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدN-body simulations in reconstruction of the kinematics of young stars in the Galaxy
366
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We try to determine the Galactic structure by comparing the observed and modeled velocities of OB-associations in the 3 kpc solar neighborhood. We made N-body simulations with a rotating stellar bar. The galactic disk in our model includes gas and stellar subsystems. The velocities of gas particles averaged over large time intervals ($sim 8$ bar rotation periods) are compared with the observed velocities of the OB-associations. Our models reproduce the directions of the radial and azimuthal components of the observed residual velocities in the Perseus and Sagittarius regions and in the Local system. The mean difference between the model and observed velocities is $Delta V=3.3$ km s$^{-1}$. The optimal value of the solar position angle $theta_b$ providing the best agreement between the model and observed velocities is $theta_b=45pm5^circ$, in good accordance with several recent estimates. The self-gravitating stellar subsystem forms a bar, an outer ring of subclass $R_1$, and slower spiral modes. Their combined gravitational perturbation leads to time-dependent morphology in the gas subsystem, which forms outer rings with elements of the $R_1$- and $R_2$-morphology. The success of N-body simulations in the Local System is likely due to the gravity of the stellar $R_1$-ring, which is omitted in models with analytical bars.
قيم البحث
اقرأ أيضاً
We use N-body simulations to examine whether a characteristic turnaround radius, as predicted from the spherical collapse model in a $rm {Lambda CDM}$ Universe, can be meaningfully identified for galaxy clusters, in the presence of full three-dimensi
onal effects. We use The Dark Sky Simulations and Illustris-TNG dark-matter--only cosmological runs to calculate radial velocity profiles around collapsed structures, extending out to many times the virial radius $R_{200}$. There, the turnaround radius can be unambiguously identified as the largest non-expanding scale around a center of gravity. We find that: (a) Indeed, a single turnaround scale can meaningfully describe strongly non-spherical structures. (b) For halos of masses $M_{200}>10^{13}M_odot$, the turnaround radius $R_{ta}$ scales with the enclosed mass $M_{ta}$ as $M_{ta}^{1/3}$, as predicted by the spherical collapse model. (c) The deviation of $R_{ta}$ in simulated halos from the spherical collapse model prediction is insensitive to halo asphericity. Rather, it is sensitive to the tidal forces due to massive neighbors when such are present. (d) Halos exhibit a characteristic average density within the turnaround scale. This characteristic density is dependent on cosmology and redshift. For the present cosmic epoch and for concordance cosmological parameters ($Omega_m sim 0.7$; $Omega_Lambda sim 0.3$) turnaround structures exhibit an average matter density contrast with the background Universe of $delta sim 11$. Thus $R_{ta}$ is equivalent to $R_{11}$ -- in a way analogous to defining the virial radius as $R_{200}$ -- with the advantage that $R_{11}$ is shown in this work to correspond to a kinematically relevant scale in N-body simulations.
We study the dynamics of faint stellar substructures around the Umbrella Galaxy, NGC 4651, which hosts a dramatic system of streams and shells formed through the tidal disruption of a nucleated dwarf elliptical galaxy. We elucidate the basic characte
ristics of the system (colours, luminosities, stellar masses) using multi-band Subaru/Suprime-Cam images. The implied stellar mass-ratio of the ongoing merger event is about 1:50. We identify candidate kinematic tracers (globular clusters, planetary nebulae, H ii regions), and follow up a subset with Keck/DEIMOS spectroscopy to obtain velocities. We find that 15 of the tracers are likely associated with halo substructures, including the probable stream progenitor nucleus. These objects delineate a kinematically cold feature in position-velocity phase space. We model the stream using single test-particle orbits, plus a rescaled pre-existing N-body simulation. We infer a very eccentric orbit with a period of roughly 0.35 Gyr and turning points at approximately 2-4 and 40 kpc, implying a recent passage of the satellite through the disc, which may have provoked the visible disturbances in the host galaxy. This work confirms that the kinematics of low surface brightness substructures can be recovered and modeled using discrete tracers - a breakthrough that opens up a fresh avenue for unraveling the detailed physics of minor merging.
In the next decade, cosmological surveys will have the statistical power to detect the absolute neutrino mass scale. N-body simulations of large-scale structure formation play a central role in interpreting data from such surveys. Yet these simulatio
ns are Newtonian in nature. We provide a quantitative study of the limitations to treating neutrinos, implemented as N-body particles, in N-body codes, focusing on the error introduced by neglecting special relativistic effects. Special relativistic effects are potentially important due to the large thermal velocities of neutrino particles in the simulation box. We derive a self-consistent theory of linear perturbations in Newtonian and non-relativistic neutrinos and use this to demonstrate that N-body simulations overestimate the neutrino free-streaming scale, and cause errors in the matter power spectrum that depend on the initial redshift of the simulations. For $z_{i} lesssim 100$, and neutrino masses within the currently allowed range, this error is $lesssim 0.5%$, though represents an up to $sim 10%$ correction to the shape of the neutrino-induced suppression to the cold dark matter power spectrum. We argue that the simulations accurately model non-linear clustering of neutrinos so that the error is confined to linear scales.
The general consensus in the N-body community is that statistical results of an ensemble of collisional N-body simulations are accurate, even though individual simulations are not. A way to test this hypothesis is to make a direct comparison of an en
semble of solutions obtained by conventional methods with an ensemble of true solutions. In order to make this possible, we wrote an N-body code called Brutus, that uses arbitrary-precision arithmetic. In combination with the Bulirsch--Stoer method, Brutus is able to obtain converged solutions, which are true up to a specified number of digits. We perform simulations of democratic 3-body systems, where after a sequence of resonances and ejections, a final configuration is reached consisting of a permanent binary and an escaping star. We do this with conventional double-precision methods, and with Brutus; both have the same set of initial conditions and initial realisations. The ensemble of solutions from the conventional simulations is compared directly to that of the converged simulations, both as an ensemble and on an individual basis to determine the distribution of the errors. We find that on average at least half of the conventional simulations diverge from the converged solution, such that the two solutions are microscopically incomparable. For the solutions which have not diverged significantly, we observe that if the integrator has a bias in energy and angular momentum, this propagates to a bias in the statistical properties of the binaries. In the case when the conventional solution has diverged onto an entirely different trajectory in phase-space, we find that the errors are centred around zero and symmetric; the error due to divergence is unbiased, as long as the time-step parameter, eta <= 2^(-5) and when simulations which violate energy conservation by more than 10% are excluded.
Stellar tidal streams are sensitive tracers of the properties of the gravitational potential in which they orbit and detailed observations of their density structure can be used to place stringent constraints on fluctuations in the potential caused b
y, e.g., the expected populations of dark matter subhalos in the standard cold dark matter paradigm (CDM). Simulations of the evolution of stellar streams in live $N$-body halos without low-mass dark-matter subhalos, however, indicate that streams exhibit significant perturbations on small scales even in the absence of substructure. Here we demonstrate, using high-resolution $N$-body simulations combined with sophisticated semi-analytic and simple analytic models, that the mass resolutions of $10^4$--$10^5,rm{M}_{odot}$ commonly used to perform such simulations cause spurious stream density variations with a similar magnitude on large scales as those expected from a CDM-like subhalo population and an order of magnitude larger on small, yet observable, scales. We estimate that mass resolutions of $approx100,rm{M}_{odot}$ ($approx1,rm{M}_{odot}$) are necessary for spurious, numerical density variations to be well below the CDM subhalo expectation on large (small) scales. That streams are sensitive to a simulations particle mass down to such small masses indicates that streams are sensitive to dark matter clustering down to these low masses if a significant fraction of the dark matter is clustered or concentrated in this way, for example, in MACHO models with masses of $10$--$100,rm{M}_{odot}$.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
