Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTheory of Stellar Population Synthesis with an application to N-Body simulations
107
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Aims. We present here a new theoretical approach to population synthesis. The aim is to predict colour magnitude diagrams (CMDs) for huge numbers of stars. With this method we generate synthetic CMDs for N-body simulations of galaxies. Sophisticated hydrodynamic N-body models of galaxies require equal quality simulations of the photometric properties of their stellar content. The only prerequisite for the method to work is very little information on the star formation and chemical enrichment histories, i.e. the age and metallicity of all star-particles as a function of time. The method takes into account the gap between the mass of real stars and that of the star-particles in N-body simulations, which best correspond to the mass of star clusters with different age and metallicity, i.e. a manifold of single stellar sopulations (SSP). Methods. The theory extends the concept of SSP to include the phase-space (position and velocity) of each star. Furthermore, it accelerates the building up of simulated CMD by using a database of theoretical SSPs that extends to all ages and metallicities of interest. Finally, it uses the concept of distribution functions to build up the CMD. The technique is independent of the mass resolution and the way the N-body simulation has been calculated. This allows us to generate CMDs for simulated stellar systems of any kind: from open clusters to globular clusters, dwarf galaxies, or spiral and elliptical galaxies. Results. The new theory is applied to an N-body simulation of a disc galaxy to test its performance and highlight its flexibility.
rate research
Read More
Self-interacting dark matter (SIDM) models have the potential to solve the small-scale problems that arise in the cold dark matter paradigm. Simulations are a powerful tool for studying SIDM in the context of astrophysics, but it is numerically challenging to study differential cross-sections that favour small-angle scattering, as in light-mediator models. Here, we present a novel approach to model frequent scattering based on an effective drag force, which we have implemented into the N-body code gadget-3. In a range of test problems, we demonstrate that our implementation accurately models frequent scattering. Our implementation can be used to study differences between SIDM models that predict rare and frequent scattering. We simulate core formation in isolated dark matter haloes, as well as major mergers of galaxy clusters and find that SIDM models with rare and frequent interactions make different predictions. In particular, frequent interactions are able to produce larger offsets between the distribution of galaxies and dark matter in equal-mass mergers.
We present SPISEA (Stellar Population Interface for Stellar Evolution and Atmospheres), an open-source Python package that simulates simple stellar populations. The strength of SPISEA is its modular interface which offers the user control of 13 input properties including (but not limited to) the Initial Mass Function, stellar multiplicity, extinction law, and the metallicity-dependent stellar evolution and atmosphere model grids used. The user also has control over the Initial-Final Mass Relation in order to produce compact stellar remnants (black holes, neutron stars, and white dwarfs). We demonstrate several outputs produced by the code, including color-magnitude diagrams, HR-diagrams, luminosity functions, and mass functions. SPISEA is object-oriented and extensible, and we welcome contributions from the community. The code and documentation are available on GitHub and ReadtheDocs, respectively.
The gravitational properties of a torus are investigated. It is shown that a torus can be formed from test particles orbiting in the gravitational field of a central mass. In this case, a toroidal distribution is achieved because of the significant spread of inclinations and eccentricities of the orbits. To investigate the self-gravity of the torus we consider the $N$-body problem for a torus located in the gravitational field of a central mass. It is shown that in the equilibrium state the cross-section of the torus is oval with a Gaussian density distribution. The dependence of the obscuring efficiency on torus inclination is found.
We present results of an infrared study of the molecular cloud Lynds 1340, forming three groups of low and intermediate-mass stars. Our goals are to identify and characterise the young stellar population of the cloud, study the relationships between the properties of the cloud and the emergent stellar groups, and integrate L1340 into the picture of the star-forming activity of our Galactic environment. We selected candidate young stellar objects from the Spitzer and WISE data bases using various published color criteria, and classified them based on the slope of the spectral energy distribution. We identified 170 Class II, 27 Flat SED, and Class 0/I sources. High angular resolution near-infrared observations of the RNO 7 cluster, embedded in L1340, revealed eight new young stars of near-infrared excess. The surface density distribution of young stellar objects shows three groups, associated with the three major molecular clumps of L1340, each consisting of less than 100 members, including both pre-main sequence stars and embedded protostars. New Herbig--Haro objects were identified in the Spitzer images. Our results demonstrate that L1340 is a prolific star-forming region of our Galactic environment in which several specific properties of the intermediate-mass mode of star formation can be studied in detail.
$N$-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be established. We address this reliability by running a three-body planetary system over about $200$ e-folding times. Using nearby initial conditions, we can construct statistics of the long-term phase-space structure and compare to rough estimates of resonant widths of the system. We compared statistics for a wide range of numerical methods, including a Runge--Kutta method, Wisdom--Holman method, symplectic corrector methods, and a method by Laskar and Robutel. Improving an integrator did not increase the phase space accuracy, but simply increasing the number of initial conditions did. In fact, the statistics of a higher order symplectic corrector method were inconsistent with the other methods in one test.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
