Search

Recently, Leitinger et al. (2023) and Cadelano et al (2024}) have shown that, in some Milky Way globular clusters (GCs), pristine red giant branch (RGB) stars are more centrally concentrated than enriched RGB stars. This contradicts most multiple-population formation scenarios, which predict that the enriched population 2p should initially be more concentrated than the pristine population 1P. We analyze a MOCCA GC model that exhibits a higher spatial concentration of 1P RGB stars than 2P RGB stars at 13 Gyr. The MOCCA models assume the asymptotic giant branch (AGB) pollution scenario, where 2P stars are initially more concentrated than 1P stars. Our results indicate that the observed spatial distributions of multiple populations, and possibly their kinematics, are significantly shaped by dynamical interactions. These interactions preferentially eject 2P RGB progenitors from the central regions, leading to a transient over-concentration of 1P RGB stars at late times. This effect is particularly relevant for GCs with present-day of a few $10^5 M_{\odot}$, which have retained only about 10 - 20 percent of their initial mass. Such clusters may appear dynamically young due to heating from a black hole subsystem, even if they have undergone significant mass loss and dynamical evolution. Additionally, the relatively small number of RGB stars in these clusters suggests that interpreting the spatial distributions of multiple populations solely from RGB stars may lead to biased conclusions about the overall distribution of 2P and 1P. The apparent over-concentration of the 1P relative to the 2P is likely a transient effect driven by the preferential removal of 2P RGB progenitors via strong dynamical encounters. MOCCA models of multiple stellar populations based on the AGB scenario may explain anomalous features observed in some Galactic GCs, such as NGC 3201 and NGC 6101., Comment: 5 pages, 3 figures, submitted to A&A letters

This paper investigates the formation, populations, and evolutionary paths of UltraLuminous X-ray Sources (ULXs) within Globular Clusters (GCs). ULXs, characterised by their extreme X-ray luminosities, present a challenge to our understanding of accretion physics and compact object formation. While previous studies have largely focused on field populations, this research examines the unique environment of GCs, where dynamical interactions play a significant role. Using the MOCCA Monte Carlo code, we explore how dynamics influences ULX populations within these dense stellar clusters. Our findings reveal that dynamical processes, such as binary hardening and exchanges, can both facilitate and impede ULX formation in GCs. The study explores the impact of parameters including the initial binary fraction, tidal filling, and multiple stellar populations on the evolution of ULXs. We find that non-tidally filling clusters exhibit significantly larger ULX populations compared to tidally filling ones. The results indicate that the apparent scarcity of ULXs in GCs may be related to the older stellar populations of GCs relative to the field. Furthermore, the study identifies a population of "escaper" ULXs, which originate in GCs but are ejected and emit X-rays outside the cluster. These escapers may significantly contribute to the observed field ULX population., Comment: 13 pages, 3 figures, submitted to A&A

Using the MOCCA code, we study the evolution of globular clusters with multiple stellar populations. For this purpose, the MOCCA code has been significantly extended to take into account the formation of an enriched population of stars from re-accreted gas with a time delay after the formation of the pristine population of stars. The possibility of cluster migration in the host galaxy and the fact that the pristine population can be described by a model, not in virial equilibrium are also taken into account. Gas re-accretion and cluster migration have a decisive impact on the observational parameters of clusters and the ratio of the number of objects between the pristine and enriched populations. The obtained results, together with observational data, suggest a speculative scenario that makes it possible to explain observational data, the correlation between the mass of the cluster and the ratio of the pristine to the enriched populations, and the observational fact that for some globular clusters, the pristine population is more concentrated than the enriched one. In this scenario, it is important to take into account the environment in which the cluster lives, the conditions in the galaxy when it formed, and the fact that a significant part of the globular clusters associated with the Galaxy come from dwarf galaxies that merged with the Milky Way. The initial conditions describing GCs in the simulations discussed in the paper are different from typical initial GC models that are widely used. Instead of GCs being highly concentrated and lying deep inside the Roche lobe, models that fill the Roche lobe are required. This carries strong constraints on where in the galaxy GCs are formed., Comment: 17 pages, 13 figures, submitted to A&A

We explore the evolution of various properties of multiple-population globular clusters (GCs) for a broad range of initial conditions. We simulated over 200 GC models using the MOCCA Monte Carlo code and find that present-day properties (core and half-light radii, ratio of the number of second-generation (SG) stars to the total number of stars, NSG/NTOT) of these models cover the observed values of these quantities for Milky Way GCs. Starting with a relatively small value of the SG fraction (NSG/NTOT ~ 0.25) and a SG system concentrated in the inner regions of the cluster, we find, in agreement with previous studies, that systems in which the first-generation (FG) is initially tidally filling or slightly tidally underfilling best reproduce the observed ratios of NSG/NTOT and have values of the core and half-light radii typical of those of many Galactic globular clusters. Models in which the FG is initially tidally underfilling retain values of NSG/NTOT close to their initial values. These simulations expand previous investigations and serve to further constrain the viable range of initial parameters and better understand their influence on present-day GC properties. The results of this investigation also provide the basis for our future survey aimed at building specific models to reproduce the observed trends (or lack thereof) between the properties of multiple stellar populations and other clusters properties., Comment: Accepted to Astronomy and Astrophysics (A&A)

There could be a significant population of double white dwarf binaries (DWDs) inside globular clusters (GCs), however, these are often too faint to be individually observed. We have utilized a large number GC models evolved with the Monte Carlo Cluster Simulator (MOCCA) code, to create a large statistical dataset of DWDs. These models include multiple-stellar populations, resulting in two distinct initial populations: one dense and another less dense. Due to the lower density of one population, a large number of objects escape during the early GC evolution, leading to a high mass-loss rate. In this dataset we have analysed three main groups of DWDs, namely in-cluster binaries, escaped binaries, and isolated evolution of primordial binaries. We compared the properties of these groups to observations of close and wide binaries. We find that the number of escaping DWDs is significantly larger than the number of in-cluster binaries and those that form via the isolated evolution of all promiridial binaries in our GC models. This suggests that dynamics play an important role in the formation of DWDs. For close binaries, we found a good agreement in the separations of escaped binaries and isolated binaries, but in-cluster binaries showed slight differences. We could not reproduce the observed extremely low mass WDs due to the limitations of our stellar and binary evolution prescriptions. For wide binaries, we also found a good agreement in the separations and masses, after accounting for observational selection effects. We conclude that, even though the current observational samples of DWDs are extremely biased and incomplete, our results compare reasonably well with observations., Comment: 13 pages, 8 figures

Because of their period-luminosity relation (PLR), classical Cepheids play a key role in the calibration of the extragalactic distance scale and the determination of the Hubble-Lema\^{i}tre constant $H_0$. Recent findings show that the majority of classical Cepheids should be in binary or multiple systems, which might undermine their accuracy, as the extra -- and unaccounted for -- light from the companions of Cepheids causes a shift in the PLR. We quantify this shift using synthetic populations of binary Cepheids that we developed for this purpose, as described in Paper I of this series. We find that while all PLRs are shifted toward brighter values due to the excess light from the companions, the bias in the relative distance modulus between two galaxies hosting binary Cepheids can be either positive or negative, depending on the percentage of binary Cepheids in them. If the binarity percentage in the two galaxies is similar, the effect of binarity is canceled. Otherwise, it introduces a shift in the distance modulus of the order of millimags in the near-infrared passbands and Wesenheit indices, and tens of millimags in the visual domain; its exact value depends on the variant of the synthetic population (a unique combination of metallicity, star formation history, shape and location of the instability strip, and initial parameter distributions). Such shifts in distance moduli to type Ia supernova host galaxies introduce an additional statistical error on $H_0$, which however does not prevent measuring $H_0$ with a precision of 1%., Comment: 16 pages, 11 figures, accepted for publication in The Astrophysical Journal

The HR 2562 system is a rare case where a brown dwarf companion resides in a cleared inner hole of a debris disk, offering invaluable opportunities to study the dynamical interaction between a substellar companion and a dusty disk. We present the first ALMA observation of the system as well as the continued GPI monitoring of the companion's orbit with 6 new epochs from 2016 to 2018. We update the orbital fit and, in combination with absolute astrometry from GAIA, place a 3$\sigma$ upper limit of 18.5 $M_J$ on the companion's mass. To interpret the ALMA observations, we used radiative transfer modeling to determine the disk properties. We find that the disk is well resolved and nearly edge on. While the misalignment angle between the disk and the orbit is weakly constrained due to the short orbital arc available, the data strongly support a (near) coplanar geometry for the system. Furthermore, we find that the models that describe the ALMA data best have an inner radius that is close to the companion's semi-major axis. Including a posteriori knowledge of the system's SED further narrows the constraints on the disk's inner radius and place it at a location that is in reasonable agreement with, possibly interior to, predictions from existing dynamical models of disk truncation by an interior substellar companion. HR\,2562 has the potential over the next few years to become a new testbed for dynamical interaction between a debris disk and a substellar companion., Comment: 22 pages, 17 figures

The POLISH2 optical polarimeter has been in operation at the Lick Observatory 3-m Shane telescope since 2011, and it was commissioned at the Gemini North 8-m in 2016. This instrument primarily targets exoplanets, asteroids, and the Crab pulsar, but it has also been used for a wide variety of planetary, galactic, and supernova science. POLISH2's photoelastic modulators, employed instead of rotating waveplates or ferro-electric liquid crystal modulators, offer the unprecedented ability to achieve sensitivity and accuracy of order 1 ppm (0.0001%), which are difficult to obtain with conventional polarimeters. Additionally, POLISH2 simultaneously measures intensity (Stokes I), linear polarization (Stokes Q and U), and circular polarization (Stokes V), which fully describe the polarization state of incident light. We document our laboratory and on-sky calibration methodology, our archival on-sky database, and we demonstrate conclusive detection of circular polarization of certain objects., Comment: 38 pages, 16 figures, accepted by ApJS

The majority of classical Cepheids are binary stars, yet the contribution of companions' light to the total brightness of the system has been assumed negligible and lacked a thorough, quantitative evaluation. We present an extensive study of synthetic populations of binary Cepheids, that aims to characterize Cepheids' companions (e.g. masses, evolutionary and spectral types), quantify their contribution to the brightness and color of Cepheid binaries, and assess the relevance of input parameters on the results. We introduce a collection of synthetic populations, which vary in metal content, initial parameter distribution, location of the instability strip edges, and star formation history. Our synthetic populations are free from the selection bias, while the percentage of Cepheid binaries is controlled by the binarity parameter. We successfully reproduce recent theoretical and empirical results: the percentage of binary Cepheids with main sequence (MS) companions, the contrast-mass ratio relation for binary Cepheids with MS companions, the appearance of binary Cepheids with giant evolved companions as outlier data points above the period-luminosity relation. Moreover, we present the first estimation of the percentage of binary Cepheids in the Large Magellanic Cloud and announce the quantification of the effect of binarity on the slope and zero-point of multiband period-luminosity relations, which will be reported in the next paper of this series., Comment: 29 pages, 19 figures (5 in the Appendix), accepted for ApJ, replaced with a version after editor's corrections

The discovery of neutron stars powering several ultraluminous X-ray sources (ULXs) raises important questions about the nature of the underlying population. In this paper we build on previous work studying simulated populations by incorporating a model where the emission originates from a precessing, geometrically beamed wind-cone, created by a super-critical inflow. We obtain estimates -- independent of the prescription for the precession period of the wind -- for the relative number of ULXs that are potentially visible (persistent or transient) for a range of underlying factors such as the relative abundance of black holes or neutron stars within the population, maximum precessional angle, and LMXB duty cycle. We make initial comparisons to existing data using a catalogue compiled from XMM-Newton. Finally, based on estimates for the precession period, we determine how the eROSITA all-sky survey (eRASS) will be able to constrain the underlying demographic., Comment: 14 pages, 8 Figures, Accepted submission MNRAS

The presence of radiatively driven outflows is well established in ultraluminous X-ray sources (ULXs). These outflows are optically thick and can reprocess a significant fraction of the accretion luminosity. Assuming isotropic emission, escaping radiation from the outflow's photosphere has the potential to irradiate the outer disc. Here, we explore how the atmosphere of the outer disc would respond to such irradiation, and specifically whether unstable heating may lead to significant mass loss via thermally-driven winds. We find that, for a range of physically relevant system parameters, this mass loss may actually switch off the inflow entirely and potentially drive limit-cycle behaviour (likely modulated on the timescale of the outer disc). In ULXs harbouring neutron stars, magnetic fields tend to have a slight destabilizing effect; for the strongest magnetic fields and highest accretion rates, this can push otherwise stable systems into the unstable regime. We explore the prevalence of the instability in a simulated sample of ULXs obtained from a binary population synthesis calculation. We find that almost all neutron star and black hole ULXs with Eddington-scaled accretion rates of $\dot{m}_0 < 100$ should be able to drive powerful outflows from their outer discs. Several known ULXs are expected to lie in this regime; the persistence of accretion in these sources implies the irradiation may be anisotropic which can be reconciled with the inferred reprocessed (optical) emission if some of this originates in the wind photosphere or irradiation of the secondary star., Comment: 8 pages, 9 figures, accepted to MNRAS

HR 6819 was recently reported to be a triple system with a dormant black hole (BH). The inner binary system was defined as a star estimated to be at the end of its main sequence and a dormant BH. As the inner binary is not resolved, the third component may actually just be spatially coinciding with the inner binary. In this study we test whether the system inner binary can be reconstructed using the isolated binary evolution in the Galactic field or through the dynamical evolution within globular star clusters. Our goal is to understand the formation of the HR 6819 inner binary. To simulate the inner binary evolution we assumed that the influence of the third body on the inner binary is negligible. We created synthetic populations of BH main-sequence binaries for the Galactic disc and the Galactic globular clusters to compare to the reported parameters of the HR 6819 inner binary. We have adopted very optimistic input physics, in terms of common envelope evolution and BH formation, for the formation of binaries similar to the reported inner HR 6819 binary. Despite our optimistic assumptions we cannot form systems like the inner HR 6819 binary in globular clusters. Even with our extreme assumptions, the formation of an HR 6819-like binary in the Galactic field population is not expected. We argue that if a dormant BH actually exists in the reported configuration inside HR 6819, its presence cannot easily be explained by our models based on isolated and dynamical binary evolution., Comment: 11 pages

The vast majority of binaries containing a compact object and a regular star spend most of their time in a quiescent state where no strong interactions occur between components. Detection of these binaries is extremely challenging and only few candidates have been detected through optical spectroscopy. Self-lensing represents a new means of detecting compact objects in binaries, where gravitational lensing of the light from the visible component by the compact object produces periodic optical flares. Here we show that current and planned large-area optical surveys can detect a significant number ($\sim 100$-$10,000$s) of these self-lensing binaries and provide insights into the properties of the compact lenses. We show that many of the predicted population of observable self-lensing binaries will be observed with multiple self-lensing flares; this both improves the chances of detection and also immediately distinguishes them from chance-alignment micro-lensing events. Through self-lensing we can investigate long - but previously hidden - stages of binary evolution and consequently provide new constraints on evolutionary models which impact on the number and nature of double compact object mergers., Comment: 11 pages, 6 figures, 3 tables

Although ultra-luminous X-ray sources (ULX) are important for astrophysics due to their extreme apparent super-Eddington luminosities, their nature is still poorly known. Theoretical and observational studies suggest that ULXs could be a diversified group of objects composed of low-mass X-ray binaries, high-mass X-ray binaries and marginally also systems containing intermediate-mass black holes, which is supported by their presence in a variety of environments. Observational data on the ULX donors could significantly boost our understanding of these systems, but only a few were detected. There are several candidates, mostly red supergiants (RSGs), but surveys are typically biased toward luminous near-infrared objects. Nevertheless, it is worth exploring if RSGs can be members of ULX binaries. In such systems matter accreted onto the compact body would have to be provided by the stellar wind of the companion, since a Roche-lobe overflow could be unstable for relevant mass-ratios. Here we present a comprehensive study of the evolution and population of wind-fed ULXs and provide a theoretical support for the link between RSGs and ULXs. Our estimated upper limit on contribution of wind-fed ULX to the overall ULX population is $\sim75$--$96\%$ for young ($<100$ Myr) star forming environments, $\sim 49$--$87\%$ for prolonged constant star formation (e.g., disk of Milky Way), and $\lesssim1\%$ for environments in which star formation ceased long time ($>2$ Gyr) ago. We show also that some wind-fed ULXs (up to $6\%$) may evolve into merging double compact objects (DCOs), but typical systems are not viable progenitors of such binaries because of their large separations. We demonstrate that, the exclusion of wind-fed ULXs from population studies of ULXs, might have lead to systematical errors in their conclusions., Comment: 14 pages, 9 figures, 2 tables, accepted to ApJ

We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The $J$, $H$, $K1$, and $K2$ spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challenging atmosphere to model with existing grids. Best-fit values vary from 1200 K to 1800 K for effective temperature and from 3.0 to 5.0 for log($g$), depending on which individual wavelength band is fit and which model suite is applied. The extreme redness of the companion can be partially reconciled by invoking a high-altitude layer of sub-micron dust particles, similar to dereddening approaches applied to the peculiar red field L-dwarf population. However, reconciling the HD 206893 B spectra with even those of the reddest low-gravity L-dwarf spectra still requires the contribution of additional atmospheric dust, potentially due to the debris disk environment in which the companion resides. Orbit fitting from four years of astrometric monitoring is consistent with a $\sim$30-year period, orbital inclination of 147$^{\circ}$, and semimajor axis of 10 au, well within the estimated disk inner radius of $\sim$50 au. As one of very few substellar companions imaged interior to a circumstellar disk, the properties of this system offer important dynamical constraints on companion-disk interaction and provide a benchmark for substellar and planetary atmospheric study., Comment: Accepted for publication in AJ (October 15, 2020). 35 pages, 20 figures, with tables and appendices presented in their entirety

Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the Startrack binary-star population modelling code and a detailed model of the Milky Way (MW) galaxy we calculate the expected number of detections for Gaia and LAMOST surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star-formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty which is introduced by methods which Monte Carlo sample from binary population synthesis output to produce one potential specific realisation of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realisations. Specifically, we predict $\sim41$-$340$ nBHBs to be observed by Gaia, although the numbers may drop to $\sim10$-$70$ if the recent ($\lesssim100\;$ Myr) star formation is low ($\sim1\;M_\odot$/yr ). For LAMOST we predict $\lesssim14$ detectable nBHBs, which is lower partially because its field-of-view covers just $\sim6\%$ of the Galaxy., Comment: 23 pages, 15 figures, 8 tables

We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $\sim 0.3$ mas/yr, and constrained the parallax of each star to be $\lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers., Comment: 25 pages, 12 figures

We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east side of the disk, peaking in intensity at 0.40" (44 AU) and extending out to 0.9" (100 AU). There is a spiral arm-like enhancement to the south, reminiscent of the known spiral structures on the outer rings of the disk. The location of the spiral arm is coincident with 12CO J=3-2 emission detected by ALMA, and hints at a dynamically active inner circumstellar region. Our observations also show a portion of the middle dusty ring at ~220 AU known from previous observations of this system. We fit the polarized H-band emission with a continuum radiative transfer Mie model. Our best-fit model favors an optically thin disk with a minimum dust grain size close to the blow-out size for this system: evidence of on-going dust production in the inner reaches of the disk. The thermal emission from this model accounts for virtually all of the far-infrared and millimeter flux from the entire HD 141569A disk, in agreement with the lack of ALMA continuum and CO emission beyond ~100 AU. A remaining 8-30 micron thermal excess a factor of ~2 above our model argues for a yet-unresolved warm innermost 5-15 AU component of the disk., Comment: 16 pages, 8 figures, accepted at AJ

We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: 12.8 [+5.3, -3.2] M_Jup. Our eccentricity determination of 0.12 [+0.04, -0.03] disfavors circular orbits. We consider orbit fits to several different imaging datasets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4., Comment: 29 Pages, 23 Figures. Accepted to AJ

We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems., Comment: 20 pages, 9 figures

The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions., Comment: 14 pages, 9 figures. Accepted for publication in the Astronomical Journal

We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$\sigma$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$)., Comment: 17 pages, 11 figures. Accepted for publication in the Astronomical Journal

We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in reduced GPI images. As well as introducing a small error in position angle measurements for targets observed at small zenith distances, this error led to a significant error in the previous astrometric calibration that has affected all subsequent astrometric measurements. We present a detailed description of these issues, and how they were corrected. We reduced GPI observations of calibration binaries taken periodically since the instrument was commissioned in 2014 using an updated version of the DRP. These measurements were compared to observations obtained with the NIRC2 instrument on Keck II, an instrument with an excellent astrometric calibration, allowing us to derive an updated plate scale and north offset angle for GPI. This revised astrometric calibration should be used to calibrate all measurements obtained with GPI for the purposes of precision astrometry., Comment: 36 pages, 17 figures, accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems. Updated version includes revisions made during the referee process

We explore the different formation channels of merging double compact objects (DCOs: BH-BH/BH-NS/NS-NS) that went through a ultraluminous X-ray phase (ULX: X-ray sources with apparent luminosity exceeding $10^{39}$ erg s$^{-1}$). There are many evolutionary scenarios which can naturally explain the formation of merging DCO systems: isolated binary evolution, dynamical evolution inside dense clusters and chemically homogeneous evolution of field binaries. It is not clear which scenario is responsible for the majority of LIGO/Virgo sources. Finding connections between ULXs and DCOs can potentially point to the origin of merging DCOs as more and more ULXs are discovered. We use the StarTrack population synthesis code to show how many ULXs will form merging DCOs in the framework of isolated binary evolution. Our merger rate calculation shows that in the local Universe typically 50% of merging BH-BH progenitor binaries have evolved through a ULX phase. This indicates that ULXs can be used to study the origin of LIGO/Virgo sources. We have also estimated that the fraction of observed ULXs that will form merging DCOs in future varies between 5% to 40% depending on common envelope model and metallicity., Comment: 15 pages, 7 figures, Accepted in MNRAS

We have obtained Hubble Space Telescope STIS and NICMOS, and Gemini/GPI scattered light images of the HD 191089 debris disk. We identify two spatial components: a ring resembling Kuiper Belt in radial extent (FWHM: ${\sim}$25 au, centered at ${\sim}$46 au), and a halo extending to ${\sim}$640 au. We find that the halo is significantly bluer than the ring, consistent with the scenario that the ring serves as the "birth ring" for the smaller dust in the halo. We measure the scattering phase functions in the 30{\deg}-150{\deg} scattering angle range and find the halo dust is both more forward- and backward-scattering than the ring dust. We measure a surface density power law index of -0.68${\pm}$0.04 for the halo, which indicates the slow-down of the radial outward motion of the dust. Using radiative transfer modeling, we attempt to simultaneously reproduce the (visible) total and (near-infrared) polarized intensity images of the birth ring. Our modeling leads to mutually inconsistent results, indicating that more complex models, such as the inclusion of more realistic aggregate particles, are needed., Comment: 29 pages, 18 figures, ApJ accepted

In large and complicated stellar systems like galaxies it is difficult to predict the number and characteristics of a black hole population. Such populations may be modelled as an aggregation of homogeneous (i.e. having uniform star formation history and the same initial chemical composition) stellar populations. Using realistic evolutionary models we predict the abundances and properties of black holes formed from binaries in these environments. We show that the black hole population will be dominated by single black holes originating from binary disruptions and stellar mergers. Furthermore, we discuss how black hole populations are influenced by such factors as initial parameters, metallicity, initial mass function, and natal kick models. As an example application of our results, we estimate that about 26 microlensing events to happen every year in the direction of the Galactic Bulge due to black holes in a survey like OGLE-IV. Our results may be used to perform in-depth studies related to realistic black hole populations, e.g. observational predictions for space survey missions like Gaia, or Einstein Probe. We prepared a publicly available database with the raw data from our simulations to be used for more in-depth studies., Comment: 20 pages, 16 figures, revised version

The Gemini Planet Imager (GPI) contains a 10-hole non-redundant mask (NRM), enabling interferometric resolution in complement to its coronagraphic capabilities. The NRM operates both in spectroscopic (integral field spectrograph, henceforth IFS) and polarimetric configurations. NRM observations were taken between 2013 and 2016 to characterize its performance. Most observations were taken in spectroscopic mode with the goal of obtaining precise astrometry and spectroscopy of faint companions to bright stars. We find a clear correlation between residual wavefront error measured by the AO system and the contrast sensitivity by comparing phase errors in observations of the same source, taken on different dates. We find a typical 5-$\sigma$ contrast sensitivity of $2-3~\times~10^{-3}$ at $\sim\lambda/D$. We explore the accuracy of spectral extraction of secondary components of binary systems by recovering the signal from a simulated source injected into several datasets. We outline data reduction procedures unique to GPI's IFS and describe a newly public data pipeline used for the presented analyses. We demonstrate recovery of astrometry and spectroscopy of two known companions to HR 2690 and HD 142527. NRM+polarimetry observations achieve differential visibility precision of $\sigma\sim0.4\%$ in the best case. We discuss its limitations on Gemini-S/GPI for resolving inner regions of protoplanetary disks and prospects for future upgrades. We summarize lessons learned in observing with NRM in spectroscopic and polarimetric modes., Comment: Accepted to AJ, 22 pages, 14 figures

We analyse the phenomenological implications of the two-families scenario on the merger of compact stars. That scenario is based on the coexistence of both hadronic stars and strange quark stars. After discussing the classification of the possible mergers, we turn to detailed numerical simulations of the merger of two hadronic stars, i.e., "first family" stars in which delta resonances and hyperons are present, and we show results for the threshold mass of such binaries, for the mass dynamically ejected and the mass of the disk surrounding the post-merger object. We compare these results with those obtained within the one-family scenario and we conclude that relevant signatures of the two-families scenario can be suggested, in particular: the possibility of a rapid collapse to a black hole for masses even smaller than the ones associated to GW170817; during the first milliseconds, oscillations of the postmerger remnant at frequencies higher than the ones obtained in the one-family scenario; a large value of the mass dynamically ejected and a small mass of the disk, for binaries of low total mass. Finally, based on a population synthesis analysis, we present estimates of the number of mergers for: two hadronic stars; hadronic star - strange quark star; two strange quark stars. We show that for unequal mass systems and intermediate values of the total mass, the merger of a hadronic star and a strange quark star is very likely (GW170817 has a possible interpretation into this category of mergers). On the other hand, mergers of two strange quark stars are strongly suppressed., Comment: 18 pages, 16 figures

We investigated the time-variability and spectral properties of the eclipsing X-ray source Circinus Galaxy X-1 (GG X-1), using Chandra, XMM-Newton and ROSAT. We phase-connected the lightcurves observed over 20 years, and obtained a best-fitting period $P = (25,970.0 \pm 0.1)$ s $\approx$7.2 hr, and a period derivative $\dot{P}/P = ( 10.2\pm4.6) \times 10^{-7}$ yr$^{-1}$. The X-ray lightcurve shows asymmetric eclipses, with sharp ingresses and slow, irregular egresses. The eclipse profile and duration vary substantially from cycle to cycle. We show that the X-ray spectra are consistent with a power-law-like component, absorbed by neutral and ionized Compton-thin material, and by a Compton-thick, partial-covering medium, responsible for the irregular dips. The high X-ray/optical flux ratio rules out the possibility that CG X-1 is a foreground Cataclysmic Variable; in agreement with previous studies, we conclude that it is the first example of a compact ultraluminous X-ray source fed by a Wolf-Rayet star or stripped Helium star. Its unocculted luminosity varies between $\approx$4 $\times 10^{39}$ erg s$^{-1}$ and $\approx$3 $\times 10^{40}$ erg s$^{-1}$. Both the donor star and the super-Eddington compact object drive powerful outflows: we suggest that the occulting clouds are produced in the wind-wind collision region and in the bow shock in front of the compact object. Among the rare sample of Wolf-Rayet X-ray binaries, CG X-1 is an exceptional target for studies of super-critical accretion and close binary evolution; it is also a likely progenitor of gravitational wave events., Comment: 27 pages, 11 figures, accepted by ApJ

Symbiotic stars are interacting binaries with one of the longest orbital periods. Since they can contain a massive white dwarf with a high accretion rate they are considered a possible type Ia supernovae (SNe Ia) progenitors. Among symbiotic binaries there are systems containing a Mira donor, which can have orbital periods of a few tens of years and more. This subclass of symbiotic stars due to their very large separation usually was not considered promising SNe Ia progenitors. We analysed evolution of one of the well studied symbiotic star with a Mira donor, V407 Cyg. We showed that the standard evolution model predicts that the system will not become a SN Ia. However, by simply adding a Wind Roche Lobe Overflow as one of the mass transfer modes we predict that the white dwarf in V407 Cyg will reach the Chandrasekhar limit in 40-200 Myr., Comment: 4 pages, 3 figures, submitted to MNRAS, comments welcome

We have analyzed how anisotropic emission of radiation affects the observed sample of ultraluminous X-ray sources (ULXs) by performing simulations of the evolution of stellar populations, employing recent developments in stellar and binary physics, and by utilizing a geometrical beaming model motivated by theory and observation. Whilst ULXs harboring black hole accretors (BH ULXs) are typically emitting isotropically, the majority of ULXs with neutron star accretors (NS ULXs) are found to be beamed. These findings confirm previous assertions that a significant fraction of ULXs are hidden from view due to a substantial misalignment of the emission beam and the line-of-sight. We find the total number of NS ULXs in regions with constant star formation, solar metallicity, and ages above ~1 Gyr to be higher than the BH ULXs, although observationally both populations are comparable. For lower metallicities BH ULX dominate both the total and observed ULX populations. As far as burst star-formation is concerned, young ULX populations are dominated by BH ULXs, but this changes as the population ages and, post star-formation, NS ULXs dominate both the observed and total population of ULXs. We also compare our simulation output to a previous analytical prediction for the relative ratio of BH to NS ULXs in idealized flux-limited observations and find broad agreement for all but the lowest metallicities. In so doing we find that in such surveys the observed ULX population should be heavily dominated by black-hole systems rather than by systems containing neutron stars., Comment: 11 pages, 8 figures, minor corrections, accepted to ApJ