Peering through the veil: near-infrared photometry and extinction for the Galactic nuclear star cluster

The aims of this work are to provide accurate photometry in multiple near-infrared broadband filters, to determine the power-law index of the extinction-law toward the central parsec of the Galaxy, to provide measurements of the absolute extinction toward the Galactic center, and finally to measure the spatial variability of extinction on arcsecond scales.We use adaptive optics observations of the central parsec of the Milky Way. Absolute values for the extinction in the H, Ks, and L'-bands as well as of the power-law indices of the H to Ks and Ks to L' extinction-laws are measured based on the well-known properties of red clump stars. Extinction maps are derived based on H-Ks and Ks-L' colors. We present Ks-band photometry for ~7700 stars (H and L' photometry for a subset). From a number of recently published values we compute a mean distance of the Galactic center of R_0=8.03+-0.15 kpc, which has an uncertainty of just 2%. Based on this R_0 and on the RC method, we derive absolute mean extinction values toward the central parsec of the Galaxy of A_H=4.48+-0.13 mag, A_Ks=2.54+-0.12$ mag, and A_L'=1.27+-0.18 mag. We estimate values of the power-law indices of the extinction-law of alpha_{H-Ks}=2.21+-0.24 and alpha_{Ks-L'}=1.34+-0.29. A Ks-band extinction map for the Galactic center is computed based on this extinction law and on stellar H-Ks colors. Mean extinction values in a circular region with 0.5" radius centered on Sagittarius A* are A_{H, SgrA*}=4.35+-0.12, A_{Ks, SgrA*}=2.46+-0.03, and A_{L', SgrA*}=1.23+-0.08.Comment: accepted for publication by Astronomy & Astrophysics; please contact RS for higher quality figure

What is the Accretion Rate in Sgr A*?

The radio source Sgr A* at the center of our Galaxy is believed to be a 2.6 x 10^6 solar mass black hole which accretes gas from the winds of nearby stars. We show that limits on the X-ray and infrared emission from the Galactic Center provide an upper limit of ~ 8 x 10^{-5} solar masses per year on the mass accretion rate in Sgr A*. The advection-dominated accretion flow (ADAF) model favors a rate < 10^{-5} solar masses per year. In comparison, the Bondi accretion rate onto Sgr A*, estimated using the observed spatial distribution of mass losing stars and assuming non-interacting stellar winds, is ~ 3 x 10^{-5} solar masses per year. There is thus rough agreement between the Bondi, the ADAF, and the X-ray inferred accretion rates for Sgr A*. We discuss uncertainties in these estimates, emphasizing the importance of upcoming observations by the Chandra X-ray observatory (CXO) for tightening the X-ray derived limits.Comment: to appear in ApJ Letter

The nature of V39: an LBV candidate or LBV impostor in the very low metallicity galaxy IC 1613?

[abridged] Context: Very few examples of luminous blue variable (LBV) stars or LBV candidates (LBVc) are known, particularly at metallicities below the SMC. The LBV phase is crucial for the evolution of massive stars, and its behavior with metallicity is poorly known. V39 in IC 1613 is a well-known photometric variable, with B-band changes larger than 1mag. over its period. The star, previously proposed to be a projection of a Galactic W Virginis and an IC 1613 red supergiant, shows features that render it a possible LBVc. Method: We investigate mid-resolution blue and red VLT-VIMOS spectra of V39, covering a time span of 40 days, and perform a quantitative analysis of the combined spectrum using the model atmosphere code CMFGEN. Results: We identify strong Balmer and FeII P-Cygni profiles, and a hybrid spectrum resembling a B-A supergiant in the blue and a G-star in the red. No significant Vrad variations are detected, and the spectral changes are small over the photometric period. Our analysis places V39 in the low-luminosity part of the LBV and LBVc region, but it is also consistent with a sgB[e] star. Conclusions: The radial velocity indicates that V39 belongs to IC 1613. The lack of Vrad changes and spectroscopic variations excludes binary scenarios. The features observed are not consistent with a W Virginis star, and this possibility is also discarded. We propose that the star is a B-A LBVc or sgB[e] star surrounded by a thick disk precessing around it. If confirmed, V39 would be the lowest metallicity resolved LBV candidate known to date. Alternatively, it could represent a new transient phase of massive star evolution, an LBV impostor.Comment: In press at A&A. 10 pages, 11 figure

An ALMA 3mm continuum census of Westerlund 1

Context. Massive stars play an important role in both cluster and galactic evolution and the rate at which they lose mass is a key driver of both their own evolution and their interaction with the environment up to and including their terminal SNe explosions. Young massive clusters provide an ideal opportunity to study a co-eval population of massive stars, where both their individual properties and the interaction with their environment can be studied in detail. Aims. We aim to study the constituent stars of the Galactic cluster Westerlund 1 in order to determine mass-loss rates for the diverse post-main sequence population of massive stars. Methods. To accomplish this we made 3mm continuum observations with the Atacama Large Millimetre/submillimetre Array. Results. We detected emission from 50 stars in Westerlund 1, comprising all 21 Wolf-Rayets within the field of view, plus eight cool and 21 OB super-/hypergiants. Emission nebulae were associated with a number of the cool hypergiants while, unexpectedly, a number of hot stars also appear spatially resolved. Conclusions. We were able to measure the mass-loss rates for a unique population of massive post-main sequence stars at every stage of evolution, confirming a significant increase as stars transitioned from OB supergiant to WR states via LBV and/or cool hypergiant phases. Fortuitously, the range of spectral types exhibited by the OB supergiants provides a critical test of radiatively-driven wind theory and in particular the reality of the bi-stability jump. The extreme mass-loss rate inferred for the interacting binary Wd1-9 in comparison to other cluster members confirmed the key role binarity plays in massive stellar evolution. The presence of compact nebulae around a number of OB and WR stars is unexpected; by analogy to the cool super-/hypergiants we attribute this to confinement and sculpting of the stellar wind via interaction with the intra-cluster medium/wind. Given the morphology of core collapse SNe depend on the nature of the pre-explosion circumstellar environment, if this hypothesis is correct then the properties of the explosion depend not just on the progenitor, but also the environment in which it is located

Measuring the Black Hole Spin in Sgr A*

The polarized mm/sub-mm radiation from Sgr A* is apparently produced by a Keplerian structure whose peak emission occurs within several Schwarzschild radii (r_S=2GM/c^2) of the black hole. The Chandra X-ray counterpart, if confirmed, is presumably the self-Comptonized component from this region. In this paper, we suggest that sub-mm timing observations could yield a signal corresponding to the period P_0 of the marginally stable orbit, and therefore point directly to the black hole's spin a. Sgr A*'s mass is now known to be (2.6\pm 0.2)\times 10^6 M_\odot (an unusually accurate value for supermassive black hole candidates), for which 2.7 min<P_0<36 min, depending on the value of a and whether the Keplerian flow is prograde or retrograde. A Schwarzschild black hole (a=0) should have P_0 ~ 20 min. The identification of the orbital frequency with the innermost stable circular orbit is made feasible by the transition from optically thick to thin emission at sub-mm wavelengths. With stratification in the emitter, the peak of the sub-mm bump in Sgr A*'s spectrum is thus produced at the smallest radius. We caution, however, that theoretical uncertainties in the structure of the emission region may still produce some ambiguity in the timing signal. Given that Sgr A*'s flux at $\nu\sim 1$ mm is several Jy, these periods should lie within the temporal-resolving capability of sub-mm telescopes using bolometric detectors. A determination of P_0 should provide not only a value of a, but it should also define the angular momentum vector of the orbiting gas in relation to the black hole's spin axis. In addition, since the X-ray flux detected by Chandra appears to be the self-Comptonized mm to sub-mm component, these temporal fluctuations may also be evident in the X-ray signal.Comment: 15 pages, 1 figures. Accepted for publication in ApJ Letter

The role of the outer boundary condition in accretion disk models: theory and application

The influence of the outer boundary condition (OBC) on the dynamics and radiation of optically thin accretion flow is investigated. Bremsstrahlung and synchrotron radiations amplified by Comptonization are taken into account and two-temperature plasma assumption is adopted. The three OBCs we adopted are the temperatures of the electrons and ions and the specific angular momentum of the accretion flow at a certain outer boundary. We find that when the general parameters such as the mass accretion rate and the viscous parameter are fixed, the peak flux at various bands such as radio, IR and X-ray, can differ by as large as several orders of magnitude under different OBCs in our example. Our results indicate that OBC is both dynamically and radiatively important therefore should be regarded as a new ``parameter'' in accretion disk models. We apply the above results to the compact radio source Sgr A* and find that the discrepancy between the mass accretion rate favored by ADAF models in the literature and that favored by the three dimensional hydrodynamical simulation can be naturally resolved by seriously considering the outer boundary condition of the accretion flow.Comment: 23 pages, 9 figures,accepted by the Astrophysical Journa

Metallicity in the Galactic Center: The Arches cluster

We present a quantitative spectral analysis of five very massive stars in the Arches cluster, located near the Galactic center, to determine stellar parameters, stellar wind properties and, most importantly, metallicity content. The analysis uses a new technique, presented here for the first time, and uses line-blanketed NLTE wind/atmosphere models fit to high-resolution near-infrared spectra of late-type nitrogen-rich Wolf-Rayet stars and OfI+ stars in the cluster. It relies on the fact that massive stars reach a maximum nitrogen abundance that is related to initial metallicity when they are in the WNL phase. We determine the present-day nitrogen abundance of the WNL stars in the Arches cluster to be 1.6% (mass fraction) and constrain the stellar metallicity in the cluster to be solar. This result is invariant to assumptions about the mass-luminosity relationship, the mass-loss rates, and rotation speeds. In addition, from this analysis, we find the age of the Arches cluster to be 2-2.5Myr, assuming coeval formation

The Large-scale Bipolar Wind in the Galactic Center

During a 9-month campaign (1996--1997), the Midcourse Space Experiment (MSX) satellite mapped the Galactic Plane at mid-infrared wavelengths (4.3--21.3um). Here we report evidence for a spectacular limb- brightened, bipolar structure at the Galactic Center extending more than a degree (170 pc at 8.0 kpc) on either side of the plane. The 8.3um emission shows a tight correlation with the 3, 6 and 11 cm continuum structure over the same scales. Dense gas and dust are being entrained in a large-scale bipolar wind powered by a central starburst. The inferred energy injection at the source is ~10^54/kappa erg for which \kappa is the covering fraction of the dusty shell (kappa <= 0.1). There is observational evidence for a galactic wind on much larger scales, presumably from the same central source which produced the bipolar shell seen by MSX. Sofue has argued that the North Polar Spur -- a thermal x-ray/radio loop which extends from the Galactic Plane to b = +80 deg -- was powered by a nuclear explosion (1-30 x 10^55 erg) roughly 15 Myr ago. We demonstrate that an open-ended bipolar wind (~10^55 erg), when viewed in near-field projection, provides the most natural explanation for the observed loop structure. The ROSAT 1.5 keV diffuse x-ray map over the inner 45 deg provides compelling evidence for this interpretation. Since the faint bipolar emission would be very difficult to detect beyond the Galaxy, the phenomenon of large-scale galactic winds may be far more common than has been observed to date.Comment: 24 pages, 6 figures, aastex. High resolution figures are available at ftp://www.aao.gov.au/pub/local/jbh/astro-ph/GC/. Astrophysical Journal, accepte

Constraints on the multiplicity of the most massive stars known: R136 a1, a2, a3, and c

The most massive stars known to date are R 136 a1, a2, a3, and c within the central cluster R 136a of the Tarantula nebula in the Large Magellanic Cloud (LMC), with reported masses in excess of 150-200$M_\odot$. However, the mass estimation of these stars relies on the assumption that they are single. We collected three epochs of spectroscopy for R 136 a1, a2, a3, and c with the Space Telescope Imaging Spectrograph (STIS) of the Hubble Space Telescope (HST) in the years 2020-2021 to probe potential radial-velocity (RV) variations. We combine these epochs with an additional HST/STIS observation taken in 2012. We use cross-correlation to quantify the RVs, and establish constraints on possible companions to these stars up to periods of ~10 yr. Objects are classified as binaries when the peak-to-peak RV shifts exceed 50 km/s, and when the RV shift is significant with respect to errors. R 136 a1, a2, and a3 do not satisfy the binary criteria and are thus classified as putatively single, although formal peak-to-peak RV variability on the level 40 km/s is noted for a3. Only R 136 c is classified as binary, in agreement with literature. We can generally rule out massive companions (M2 > ~50 Msun) to R 136 a1, a2, and a3 out to orbital periods of < 1 yr (separations < 5 au) at 95% confidence, or out to tens of years (separations < ~100 au) at 50% confidence. Highly eccentric binaries (e > ~0.9) or twin companions with similar spectra could evade detection down to shorter periods (> ~10 d), though their presence is not supported by the relative X-ray faintness of R 136 a1, a2, and a3. We derive a preliminary orbital solution with a 17.2 d period for the X-ray bright binary R 136 c, though more data are needed to conclusively derive its orbit. Our study supports a lower bound of 150-200 $M_\odot$ on the upper-mass limit at LMC metallicityComment: Accepted to A&