Stellar Number Density Research Articles

An Unexpected Detection of Bifurcated Blue Straggler Sequences in the Young Globular Cluster NGC 2173∗

Abstract The bifurcated patterns in the color–magnitude diagrams of blue straggler stars (BSSs) have attracted significant attention. This type of special (but rare) pattern of two distinct blue straggler sequences is commonly interpreted as evidence that cluster core-collapse-driven stellar collisions are an efficient formation mechanism. Here, we report the detection of a bifurcated blue straggler distribution in a young Large Magellanic Cloud cluster, NGC 2173. Because of the cluster’s low central stellar number density and its young age, dynamical analysis shows that stellar collisions alone cannot explain the observed BSSs. Therefore, binary evolution is instead the most viable explanation of the origin of these BSSs. However, the reason why binary evolution would render the color–magnitude distribution of BSSs bifurcated remains unclear.

Hot Jupiters Driven by High-eccentricity Migration in Globular Clusters

Abstract Hot Jupiters (HJs) are short-period giant planets that are observed around of solar-type field stars. One possible formation scenario for HJs is high-eccentricity (high-e) migration, in which the planet forms at much larger radii, is excited to high eccentricity by some mechanism, and migrates to its current orbit due to tidal dissipation occurring near periapsis. We consider high-e migration in dense stellar systems such as the cores of globular clusters (GCs), in which encounters with passing stars can excite planets to the high eccentricities needed to initiate migration. We study this process via Monte Carlo simulations of encounters with a star+planet system including the effects of tidal dissipation, using an efficient regularized restricted three-body code. HJs are produced in our simulations over a significant range of the stellar number density . Assuming the planet is initially on a low-eccentricity orbit with semimajor axis 1 au, for the encounter rate is too low to induce orbital migration, whereas for HJ formation is suppressed because the planet is more likely ejected from its host star, tidally disrupted, or transferred to a perturbing star. The fraction of planets that are converted to HJs peaks at for intermediate number densities of . Warm Jupiters, giant planets with periods between 10 and 100 days, are produced in our simulations with an efficiency of up to . Our results suggest that HJs can form through high-e migration induced by stellar encounters in the centers of of dense GCs, but not in their outskirts where the densities are lower.

Galactic disk structure as revealed by LAMOST A stars

AbstractUsing the spectroscopic distances of over 0.12 million A-type stars selected from the LAMOST Spectroscopic Survey of the Galactic Anti-center (LSS-GAC), we map their three-dimensional number density distributions in the Galaxy. These stellar number density maps allow an investigation of the Galactic young age thin disk structure with no a priori assumptions about the functional form of its components. The data show strong evidence for a significant flaring young disk. A more detail analysis show that the stellar flaring have different behaviours between the Northern and the Southern Galactic disks. The maps also reveal spatially coherent, kpc-scale stellar substructure in the thin disk. Finally, we detect the Perseus arm stellar overdensity at R ~ 10 kpc.

ASTROBIOLOGICAL EFFECTS OF GAMMA-RAY BURSTS IN THE MILKY WAY GALAXY

A planet having protective ozone within the collimated beam of a Gamma Ray Burst (GRB) may suffer ozone depletion, potentially causing a mass extinction event to existing life on a planet's surface and oceans. We model the dangers of long GRBs to planets in the Milky Way and utilize a static statistical model of the Galaxy that matches major observable properties, such as the inside-out star formation history, metallicity evolution, and 3-dimensional stellar number density distribution. The GRB formation rate is a function of both the star formation history and metallicity; however, the extent to which chemical evolution reduces the GRB rate over time in the Milky Way is still an open question. Therefore, we compare the damaging effects of GRBs to biospheres in the Milky Way using two models. One model generates GRBs as a function of the inside-out star formation history. The other model follows the star formation history, but generates GRB progenitors as a function of metallicity, thereby favoring metal-poor host regions of the Galaxy over time. If the GRB rate only follows the star formation history, the majority of the GRBs occur in the inner Galaxy. However, if GRB progenitors are constrained to low metallicity environments, then GRBs only form in the metal-poor outskirts at recent epochs. Interestingly, over the past 1 Gyr, the surface density of stars (and their corresponding planets) that survive a GRB is still greatest in the inner galaxy in both models. The present day danger of long GRBs to life at the solar radius ($R_\odot=8$ kpc) is low. We find that at least $\sim$65% of stars survive a GRB over the past 1 Gyr. Furthermore, when the GRB rate was expected to have been enhanced at higher redshifts, such as $z\gtrsim0.5$, our results suggest that a large fraction of planets would have survived these lethal GRB events.

Dust in three dimensions in the Galactic plane

We present three dimensional maps in monochromatic extinction $A_{\rm 0}$ and the extinction parameter $R_0$ within a few degrees of the Galactic plane. These are inferred using photometry from the Pan-STARRS1 and Spitzer Glimpse surveys of nearly $20$ million stars located in the region $l = 0-250^{\circ}$ and from $b = -4.5^{\circ}$ to $b=4.5^{\circ}$. Given the available stellar number density, we use an angular resolution of $7' \times 7'$ and steps of $1{\rm mag}$ in distance modulus. We simultaneously estimate distance modulus and effective temperature $T_{\rm eff}$ alongside the other parameters for stars individually using the method of \citet{Hanson2014} before combining these estimates to a complete map. The full maps are available via the MNRAS website.

Low-end mass function of the Quintuplet cluster

The Quintuplet and Arches clusters, which were formed in the harsh environment of the Galactic Center (GC) a few million years ago, have been excellent targets for studying the effects of a star-forming environment on the initial mass function (IMF). In order to estimate the shape of the low-end IMF of the Arches cluster, Shin & Kim devised a novel photometric method that utilizes pixel intensity histograms (PIHs) of the observed images. Here, we apply the PIH method to the Quintuplet cluster and estimate the shape of its low-end IMF below the magnitude of completeness limit as set by conventional photometry. We found that the low-end IMF of the Quintuplet is consistent with that found for the Arches cluster--Kroupa MF, with a significant number of low-mass stars below 1 $M_{\odot}$. We conclude that the most likely IMFs of the Quintuplet and the Arches clusters are not too different from the IMFs found in the Galactic disc. We also find that the observed PIHs and stellar number density profiles of both clusters are best reproduced when the clusters are assumed to be at three-dimensional distances of approximately 100 pc from the GC.

Number density distribution of near-infrared sources on a sub-degree scale in the Galactic center: Comparison with the Fe xxv Kα line at 6.7 keV

Abstract The stellar distribution derived from an H- and KS-band survey of the central region of our Galaxy is compared with the Fe xxv Kα (6.7 keV) line intensity observed with the Suzaku satellite. The survey is for the galactic coordinates |l| ≲ 3$_{.}^{\circ}$0 and |b | ≲ 1$_{.}^{\circ}$0 (equivalent to 0.8 kpc × 0.3 kpc for R⊙ = 8 kpc), and the number-density distribution N(KS,0; l, b) of stars is derived by using the extinction-corrected magnitude KS,0 = 10.5. This is deep enough to probe the old red-giant population and in turn to estimate the (l, b) distribution of faint X-ray point sources such as coronally active binaries and cataclysmic variables. In the Galactic plane (b = 0°), N(10.5; l, b) increases in the direction of the Galactic center as |l|−0.30±0.03 in the range of − 0$_{.}^{\circ}$1 ≥ l ≥ − 0$_{.}^{\circ}$7, but this increase is significantly slower than the increase (|l|−0.44±0.02) of the Fe xxv Kα line intensity. If normalized with the ratios in the outer region 1$_{.}^{\circ}$5 ≤ |l| ≤ 2$_{.}^{\circ}$8, where faint X-ray point sources are argued to dominate the diffuse Galactic X-ray ridge emission, the excess of the Fe xxv Kα line intensity over the stellar number density is at least a factor of two at |l| = 0$_{.}^{\circ}$1. This indicates that a significant part of the Galactic-center diffuse emission arises from a truly diffuse optically thin thermal plasma, and not from an unresolved collection of faint X-ray point sources related to the old stellar population.

An analysis of the evolving comoving number density of galaxies in hydrodynamical simulations

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.

Probing massive stars around gamma-ray burst progenitors

Long Gamma-Ray Bursts (GRBs) are produced by ultra-relativistic jets launched from core collapse of massive stars. Most massive stars form in binaries and/or in star clusters, which means that there may be a significant external photon field (EPF) around the GRB progenitor. We calculate the inverse-Compton scattering of EPF by the hot electrons in the GRB jet. Three possible cases of EPFs are considered: the progenitor is (I) in a massive binary system, (II) surrounded by a Wolf-Rayet-star wind, and (III) in a dense star cluster. Typical luminosities of 10^46 - 10^50 erg/s in the 1 - 100 GeV band are expected, depending on the stellar luminosity, binary separation (I), wind mass loss rate (II), stellar number density (III), etc. We calculate the lightcurve and spectrum in each case, taking fully into account the equal-arrival time surfaces and possible pair-production absorption with the prompt gamma-rays. Observations can put constraints on the existence of such EPFs (and hence on the nature of GRB progenitors) and on the radius where the jet internal dissipation process accelerates electrons.

Galactic disk bulk motions as revealed by the LSS-GAC DR2

We report a detailed investigation of the bulk motions of the nearby Galactic stellar disk, based on three samples selected from the LSS-GAC DR2: a global sample containing 0.57 million FGK dwarfs out to ∼2 kpc, a local subset of the global sample consisting of ∼5400 stars within 150 pc, and an anti-center sample containing ∼4400 AFGK dwarfs and red clump stars within windows a few degrees wide centered on the Galactic Anti-center. The global sample is used to construct a three-dimensional map of bulk motions of the Galactic disk from the solar vicinity out to ∼2 kpc with a spatial resolution of ∼250 pc. Typical values of the radial and vertical components of bulk motion range from −15 km s−1 to 15 km s−1; in contrast, the lag behind the circular motion dominates the azimuthal component by up to ∼15 km s−1. The map reveals spatially coherent, kpc-scale stellar flows in the disk, with typical velocities of a few tens of km s−1. Bending- and breathing-mode perturbations are clearly visible, and vary smoothly across the disk plane. Our data also reveal higher-order perturbations, such as breaks and ripples, in the profiles of vertical motion versus height. From the local sample, we find that stars from different populations exhibit very different patterns of bulk motion. Finally, the anti-center sample reveals a number of peaks in stellar number density in the line-of-sight velocity versus distance distribution, with the nearer ones apparently related to the known moving groups. The “velocity bifurcation” reported by Liu et al. at Galactocentric radii 10–11 kpc is confirmed. However, just beyond this distance, our data also reveal a new triple-peaked structure.

Complex Stellar System ESO65SC03: Open Cluster or Remnant?

AbstractWe present a complete spatial and dynamical study of the poorly populated stellar system ESO65SC03. The radial distribution of the system gives a core and cluster radii of 1.10±0.63 and 5.36±0.24 arcmin, respectively. The SNDP does not show any clear enhancement of the surface stellar number density between the stars of the system and the field regions. We derive the optimum isochrone solution for a particular grid size in the CMD using the statistical cleaning procedure. Using the statistically cleaned CMDs, we find the distance modulus, (m − M)0, and reddening, E(B − V), of the system to be 11.8±0.2 and 0.45 mag, respectively. The mean proper motion of this system is − 5.37±0.81 mas yr−1 and 0.31±0.40 in RA and DEC directions, respectively. The mean proper motion of this system is found to be almost similar to the field region. The mass function for the brighter stars is found to be too high for the system to be an open cluster. These combined results place constraints on whether stellar system ESO65SC03 is a POCR or an Asterism. Our understanding is that the ESO65SC03 is in a stage of POCR by losing their main-sequence stars in the dynamic evolution processes.

A census of stellar mass in ten massive haloes atz~ 1 from the GCLASS Survey

We study the stellar mass content of massive haloes in the redshift range 0.86<z<1.34, by measuring: (1) The stellar mass in the central galaxy versus total dynamical halo mass. (2) The total stellar mass (including satellites) versus total halo mass. (3) The radial stellar mass and number density profiles for the ensemble halo. We use a Ks-band selected catalogue for the 10 clusters in the Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS), with photometric redshifts and stellar masses measured from 11-band SED fitting. Combining the photometric catalogues with the deep spectroscopic component of GCLASS, we correct the cluster galaxy sample for interlopers. We also perform a dynamical analysis of the cluster galaxies to estimate the halo mass M200 for each cluster based on a measurement of its velocity dispersion. (1) We find that the central galaxy stellar mass fraction decreases with total halo mass, and that this is in reasonable quantitative agreement with measurements from abundance matching studies at z~1. (2) The total stellar mass fractions of these systems decrease with halo mass, indicating that lower mass systems are more efficient at transforming baryons into stars. We find the total stellar mass to be a good proxy for total halo mass, with a small intrinsic scatter. When we compare these results with literature measurements, we find that the stellar mass fraction at fixed halo mass shows no significant evolution in the range 0<z<1. (3) We measure a relatively high NFW concentration parameter c ~ 7 for the stellar mass distribution in these clusters. A simple model shows that the stellar mass content of GCLASS can evolve in typical distributions observed at lower redshifts if the clusters primarily accrete stellar mass onto the outskirts. [Abridged]

THE STELLAR NUMBER DENSITY DISTRIBUTION IN THE LOCAL SOLAR NEIGHBORHOOD IS NORTH-SOUTH ASYMMETRIC

We study the number density distribution of a sample of K and M dwarf stars, matched North and South of the Galactic plane within a distance of 2 kpc from the sun, using observations from the Ninth Data Release of the Sloan Digital Sky Survey. We determine distances using the photometric parallax method, and in this context systematic effects exist which could potentially impact the determination of the number density profile with height from the Galactic plane --- and ultimately affect a number density North-South asymmetry. They include: (i) the calibration of the various photometric parallax relations, (ii) the ability to separate dwarfs from giants in our sample, (iii) the role of stellar population differences such as age and metallicity, (iv) the ability to determine the offset of the sun from the Galactic plane, and (v) the correction for reddening from dust in the Galactic plane, though our stars are at high Galactic latitudes. We find the various analyzed systematic effects to have a negligible impact on our observed asymmetry, and using a new and larger sample of stars we confirm and refine the earlier discovery of Widrow et al. of a significant Galactic North-South asymmetry in the stellar number density distribution.

The origin of the Galactic center diffuse X-ray emission investigated by near-infrared imaging and polarimetric observations

AbstractThe origin of the Galactic center diffuse X-ray emission (GCDX) is still under intense investigation. We have found a clear excess in a longitudinal GCDX profile over a stellar number density profile in the nuclear bulge region, suggesting a significant contribution of diffuse, interstellar hot plasma to the GCDX. We have estimated that contributions of an old stellar population to the GCDX are ∼50% and ∼20% in the nuclear stellar disk and nuclear star cluster, respectively. Our near-infrared polarimetric observations show that the GCDX region is permeated by a large scale, toroidal magnetic field. Together with observed magnetic field strengths in nearly energy equipartition, the interstellar hot plasma could be confined by the toroidal magnetic field.

MAGNETICALLY CONFINED INTERSTELLAR HOT PLASMA IN THE NUCLEAR BULGE OF OUR GALAXY

The origin of the Galactic center diffuse X-ray emission (GCDX) is still under intense investigation. In particular, the interpretation of the hot (kT ~ 7 keV) component of the GCDX, characterised by the strong Fe 6.7 keV line emission, has been contentious. If the hot component originates from a truly diffuse interstellar plasma, not a collection of unresolved point sources, such plasma cannot be gravitationally bound, and its regeneration would require a huge amount of energy. Here we show that the spatial distribution of the GCDX does NOT correlate with the number density distribution of an old stellar population traced by near-infrared light, strongly suggesting a significant contribution of the diffuse interstellar plasma. Contributions of the old stellar population to the GCDX are implied to be about 50 % and 20 % in the Nuclear stellar disk and Nuclear star cluster, respectively. For the Nuclear stellar disk, a scale height of 0.32 +- 0.02 deg is obtained for the first time from the stellar number density profiles. We also show the results of the extended near-infrared polarimetric observations in the central 3 deg * 2 deg region of our Galaxy, and confirm that the GCDX region is permeated by a large scale, toroidal magnetic field as previously claimed. Together with observed magnetic field strengths close to energy equipartition, the hot plasma could be magnetically confined, reducing the amount of energy required to sustain it.

The SLoWPoKES catalog of low‐mass ultra–wide binaries: A cool stars resource for testing fundamental properties and for constraining binary formation theory

AbstractWe present results from the Sloan Low‐mass Wide Pairs of Kinematically Equivalent Stars (SLoWPoKES) catalog of ultra‐wide (103–105.5 AU), low‐mass (K5–M7) common proper motion binaries. We constructed a Galactic model, based on empirical stellar number density and 3D velocity distributions, to select bona fide pairs with probability of chance alignment &lt;5 %, making SLoWPoKES an efficient sample for followup observations. Our initial catalog contains 1342 disk dwarf, subdwarf, and white dwarf‐red dwarf systems and is the largest collection of low‐mass, wide binaries ever assembled. The diversity – in mass, metallicity, age, and evolutionary states – of SLoWPoKES pairs makes it a valuable resource of coeval laboratories to examine and constrain the physical properties of low‐mass stars. SLoWPoKES pairs show signatures of two (or more) formation modes in the distribution of the physical separation and higher‐order multiplicity. Neither dynamical dissipation of primordial triples/quadruples or dynamical capture of ejected stars can explain the observed populations by itself. We use follow‐up spectroscopic observations to recalibrate the metallicity‐sensitive ζTiO/CaH index by assuming that both members of the binary system have the same composition. Our new formulation is a significantly better tracer of absolute metallicity, particularly for the early‐type M dwarfs. The catalogs are publicly available on a custom data visualization portal. (© 2013 WILEY‐VCH Verlag GmbH &amp; Co. KGaA, Weinheim)

CONSTRAINTS ON THE SHAPE OF THE MILKY WAY DARK MATTER HALO FROM JEANS EQUATIONS APPLIED TO SLOAN DIGITAL SKY SURVEY DATA

We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body + SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of ~20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q_DM = 0.47 +/- 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST.

The S-star cluster at the center of the Milky Way

Sagittarius A*, the super-massive black hole at the center of the Milky Way, is surrounded by a small cluster of high velocity stars, known as the S-stars. We aim to constrain the amount and nature of stellar and dark mass associated with the cluster in the immediate vicinity of Sagittarius A*. We use near-infrared imaging to determine the $K_\mathrm{s}$-band luminosity function of the S-star cluster members, and the distribution of the diffuse background emission and the stellar number density counts around the central black hole. This allows us to determine the stellar light and mass contribution expected from the faint members of the cluster. We then use post-Newtonian N-body techniques to investigate the effect of stellar perturbations on the motion of S2, as a means of detecting the number and masses of the perturbers. We find that the stellar mass derived from the $K_\mathrm{s}$-band luminosity extrapolation is much smaller than the amount of mass that might be present considering the uncertainties in the orbital motion of the star S2. Also the amount of light from the fainter S-cluster members is below the amount of residual light at the position of the S-star cluster after removing the bright cluster members. If the distribution of stars and stellar remnants is strongly enough peaked near Sagittarius A*, observed changes in the orbital elements of S2 can be used to constrain both their masses and numbers. Based on simulations of the cluster of high velocity stars we find that at a wavelength of 2.2 $\mu$m close to the confusion level for 8 m class telescopes blend stars will occur (preferentially near the position of Sagittarius A*) that last for typically 3 years before they dissolve due to proper motions.

The impact of high spatial frequency atmospheric distortions on weak-lensing measurements

High precision cosmology with weak gravitational lensing requires a precise measure of the Point Spread Function across the imaging data where the accuracy to which high spatial frequency variation can be modelled is limited by the stellar number density across the field. We analyse dense stellar fields imaged at the Canada-France-Hawaii Telescope to quantify the degree of high spatial frequency variation in ground-based imaging Point Spread Functions and compare our results to models of atmospheric turbulence. The data shows an anisotropic turbulence pattern with an orientation independent of the wind direction and wind speed. We find the amplitude of the high spatial frequencies to decrease with increasing exposure time as $t^{-1/2}$, and find a negligibly small atmospheric contribution to the Point Spread Function ellipticity variation for exposure times $t>180$ seconds. For future surveys analysing shorter exposure data, this anisotropic turbulence will need to be taken into account as the amplitude of the correlated atmospheric distortions becomes comparable to a cosmological lensing signal on scales less than $\sim 10$ arcminutes. This effect could be mitigated, however, by correlating galaxy shear measured on exposures imaged with a time separation greater than 50 seconds, for which we find the spatial turbulence patterns to be uncorrelated.

A Model of Habitability Within the Milky Way Galaxy

We present a model of the galactic habitable zone (GHZ), described in terms of the spatial and temporal dimensions of the Galaxy that may favor the development of complex life. The Milky Way galaxy was modeled using a computational approach by populating stars and their planetary systems on an individual basis by employing Monte Carlo methods. We began with well-established properties of the disk of the Milky Way, such as the stellar number density distribution, the initial mass function, the star formation history, and the metallicity gradient as a function of radial position and time. We varied some of these properties and created four models to test the sensitivity of our assumptions. To assess habitability on the galactic scale, we modeled supernova rates, planet formation, and the time required for complex life to evolve. Our study has improved on other literature on the GHZ by populating stars on an individual basis and modeling Type II supernova (SNII) and Type Ia supernova (SNIa) sterilizations by selecting their progenitors from within this preexisting stellar population. Furthermore, we considered habitability on tidally locked and non-tidally locked planets separately and studied habitability as a function of height above and below the galactic midplane. In the model that most accurately reproduces the properties of the Galaxy, the results indicate that an individual SNIa is ∼5.6× more lethal than an individual SNII on average. In addition, we predict that ∼1.2% of all stars host a planet that may have been capable of supporting complex life at some point in the history of the Galaxy. Of those stars with a habitable planet, ∼75% of planets are predicted to be in a tidally locked configuration with their host star. The majority of these planets that may support complex life are found toward the inner Galaxy, distributed within, and significantly above and below, the galactic midplane.