Earthshine observations offer a unique opportunity to study Earth as an exoplanet seen from the Moon. As the Sun-Earth-Moon geometry changes, Earth can be observed as a spatially unresolved exoplanet at different phase angles, providing important context for future observations of Earth-like exoplanets. Here, we present a catalog of Earthshine polarization spectra obtained with FORS2 on the VLT, covering diverse scenes, surface conditions, cloud properties, and weather patterns for over a decade. For the first time, we model this extensive dataset in detail using a homogeneous modeling framework. Previous efforts to model some of these spectra struggled to reproduce the observed polarization continuum, even with advanced 3D radiative transfer models incorporating satellite-derived surface and atmospheric data. We improve upon this with a 3D model that includes subgrid cloud variability, wavelength-dependent surface albedo maps, and an accurate treatment of ocean glint. Our simulations successfully reproduce most observed spectra to a much higher precision than previously possible. Our statistical analysis reveals that the spectral slope in the visible can distinguish between ocean and mixed surfaces in both reflected and polarized light, which is not possible using broadband filters alone. Polarized light at large phase angles, beyond the Rayleigh scattering regime, is particularly effective in differentiating oceans from land, unlike reflected light. We also identify correlations between cloud optical thickness and the polarized spectral slope, and between cloud cover and broadband B-R differences in reflected light, demonstrating the diagnostic power of these observations. This work highlight the potential of polarization for characterizing Earth-like exoplanets. From polarization alone, we can infer oceans, vegetation, and an active water cycle, key indicators of a habitable planet.

The discovery of a significant number of rapidly rotating low mass stars showing no or few flares in TESS observations was a surprise as rapid rotation has previously been taken as implying high stellar activity. Here we present TESS and HiPERCAM $u_{s}g_{s}r_{s}i_{s}z_{s}$ observations of one of these stars LP 89--187 which has a rotation period of 0.117 d. TESS data covering three sectors (64.6 d) only show three flares which have energies a few $\times10^{33}$ erg, whilst HiPERCAM observations, which cover 0.78 of the rotation period, show no evidence for flares more energetic than $\sim10^{31}$ erg. Intriguingly, other surveys show LP 89--187 has shown weak H$\alpha$ in emission. We compare the flare energy distribution of LP 89--187 with low mass stars in the $\beta$ Pic moving group, which have an age of $\sim$24 Myr. We find LP 89--187 has a lower flare rate than the $\beta$ Pic stars. In addition, we find that TRAPPIST-1 analogue stars, which are likely significantly older than the $\beta$ Pic stars, show fewer flares with energies $>10^{33}$ erg in TESS data. We examine the relationship between amplitude and period for a sample of low mass stars and find that more rapid rotators have a higher amplitude.

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos.

In the last few decades planet search surveys have been focusing on solar type stars, and only recently the high-mass regimes. This is mostly due to challenges arising from the lack of instrumental precision, and more importantly, the inherent active nature of fast rotating massive stars. Here we report NGTS-33b (TOI-6442b), a super-Jupiter planet with mass, radius and orbital period of 3.6 $\pm$ 0.3 M$_{\rm jup}$, 1.64 $\pm$ 0.07 R$_{\rm jup}$ and $2.827972 \pm 0.000001$ days, respectively. The host is a fast rotating ($0.6654 \pm 0.0006$ day) and hot (T$_{\rm eff}$ = 7437 $\pm$ 72 K) A9V type star, with a mass and radius of 1.60 $\pm$ 0.11 M$_{\odot}$ and 1.47 $\pm$ 0.06 R$_{\odot}$, respectively. Planet structure and Gyrochronology models shows that NGTS-33 is also very young with age limits of 10-50 Myr. In addition, membership analysis points towards the star being part of the Vela OB2 association, which has an age of $\sim$ 20-35 Myr, thus providing further evidences about the young nature of NGTS-33. Its low bulk density of 0.19$\pm$0.03 g cm$^{-3}$ is 13$\%$ smaller than expected when compared to transiting hot Jupiters with similar masses. Such cannot be solely explained by its age, where an up to 15$\%$ inflated atmosphere is expected from planet structure models. Finally, we found that its emission spectroscopy metric is similar to JWST community targets, making the planet an interesting target for atmospheric follow-up. Therefore, NGTS-33b's discovery will not only add to the scarce population of young, massive and hot Jupiters, but will also help place further strong constraints on current formation and evolution models for such planetary systems.

Supermassive black hole binaries may be detectable by an upcoming suite of gravitational wave experiments. Their binary nature can also be revealed by radio jets via a short-period precession driven by the orbital motion as well as the geodetic precession at typically longer periods. We have investigated Karl G. Jansky Very Large Array (VLA) and MERLIN radio maps of powerful jet sources for morphological evidence of geodetic precession. For perhaps the best studied source, Cygnus A, we find strong evidence for geodetic precession. Projection effects can enhance precession features, for which we find indications in strongly projected sources. For a complete sample of 33 3CR radio sources we find strong evidence for jet precession in 24 cases (73 per cent). The morphology of the radio maps suggests that the precession periods are of the order of 10^6 - 10^7 yr. We consider different explanations for the morphological features and conclude that geodetic precession is the best explanation. The frequently observed gradual jet angle changes in samples of powerful blazars can be explained by orbital motion. Both observations can be explained simultaneously by postulating that a high fraction of powerful radio sources have sub-parsec supermassive black hole binaries. We consider complementary evidence and discuss if any jetted supermassive black hole with some indication of precession could be detected as individual gravitational wave source in the near future. This appears unlikely, with the possible exception of M87.

We observe the magnetic field morphology towards a nearby star-forming filamentary cloud, G202.3+2.5, by the JCMT/POL-2 850 {\mu}m thermal dust polarization observation with an angular resolution of 14.4" (~0.053 pc). The average magnetic field orientation is found to be perpendicular to the filaments while showing different behaviors in the four subregions, suggesting various effects from filaments' collision in these subregions. With the kinematics obtained by N2H+ observation by IRAM, we estimate the plane-of-sky (POS) magnetic field strength by two methods, the classical Davis-Chandrasekhar-Fermi (DCF) method and the angular dispersion function (ADF) method, B_{pos,dcf} and B_{pos,adf} are ~90 {\mu}G and ~53 {\mu}G. We study the relative importance between the gravity (G), magnetic field (B) and turbulence (T) in the four subregions, find G > T > B, G >= T > B, G ~ T > B and T > G > B in the north tail, west trunk, south root and east wing, respectively. In addition, we investigate the projection effect on the DCF and ADF methods based on a similar simulation case and find the 3D magnetic field strength may be underestimated by a factor of ~3 if applying the widely-used statistical B_{pos}-to-B_{3D} factor when using DCF or ADF method, which may further underestimate/overestimate related parameters.

We present new K-band spectroscopy for the giant elliptical galaxy M87 in the Virgo cluster, taken with the LUCI spectrograph at the Large Binocular Telescope (LBT). The new data are used to study line-strengths of K-band absorption features from different chemical species, namely Fe, Mg, Ca, Na, and CO, as a function of galactocentric distance, out to 40arcsec from the center (about half of the galaxy effective radius). The radial trends of spectral indices are compared to those for the bulge of M31, observed with the same instrument. For M87, most K-band indices exhibit flat radial profiles, with the exception of NaI2.21, that decreases outwards, with a negative radial gradient. Significant offsets are found between indices for M87 and those for the bulge of M31, the latter having weaker line-strengths for almost all features, but Fe and Ca, for which we find similar trends in both systems. We find that the behavior of CO features - most prominent in giant stars - is difficult to explain, consistent with previous results for the central regions of massive galaxies. In particular, the CO indices are stronger in M87 than M31, and do not exhibit significant radial gradients in M87, despite its IMF being bottom heavier than M31 especially in its central region. Predictions of state-of-the-art stellar population models, based on results from the optical spectral range, are able to match only the Na and Ca indices of M87, while a significant mismatch is found for all other indices. This shows that state-of-the-art stellar population models should be improved significantly in order to provide reliable constraints on the stellar population content of galaxies in the NIR spectral range.

One of the most poorly understood aspects of low-mass star formation is how multiple-star systems are formed. Here we present the results of Atacama Large Millimeter/submillimeter Array (ALMA) Band-6 observations towards a forming quadruple protostellar system, G206.93-16.61E2, in the Orion B molecular cloud. ALMA 1.3 mm continuum emission reveals four compact objects, of which two are Class I young stellar objects (YSOs), and the other two are likely in prestellar phase. The 1.3 mm continuum emission also shows three asymmetric ribbon-like structures that are connected to the four objects, with lengths ranging from $\sim$500 au to $\sim$2200 au. By comparing our data with magneto-hydrodynamic (MHD) simulations, we suggest that these ribbons trace accretion flows and also function as gas bridges connecting the member protostars. Additionally, ALMA CO J=2-1 line emission reveals a complicated molecular outflow associated with G206.93-16.61E2 with arc-like structures suggestive of an outflow cavity viewed pole-on.

The star formation efficiency (SFE) has been shown to vary across different environments, particularly within galactic starbursts and deep within the bulges of galaxies. Various quenching mechanisms may be responsible, ranging from galactic dynamics to feedback from active galactic nuclei (AGN). Here, we use spatially-resolved observations of warm ionised gas emission lines from the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope (CFHT) and cold molecular gas (CO(2-1)) from the Atacama Large Millimeter/sub-millimeter Array (ALMA) to study the SFE in the bulge of the AGN-host galaxy NGC 3169. After distinguishing star-forming regions from AGN-ionised regions using emission-line ratio diagnostics, we measure spatially-resolved molecular gas depletion times (\tau_dep = 1/SFE) with a spatial resolution of \approx 100 pc within a galactocentric radius of 1.8 kpc. We identify a star-forming ring located at radii 1.25 \pm 0.6 kpc with an average \tau_dep of 0.3 Gyr. At radii < 0.9 kpc, however, the molecular gas surface densities and depletion times increase with decreasing radius, the latter reaching approximately 2.3 Gyr at a radius \approx 500 pc. Based on analyses of the gas kinematics and comparisons with simulations, we identify AGN feedback, bulge morphology and dynamics as the possible causes of the radial profile of SFE observed in the central region of NGC 3169.

We report the discovery of ancient massive merger events in the early-type galaxies NGC 1380 and NGC 1427, members of the Fornax galaxy cluster. Both galaxies have been observed by the MUSE IFU instrument on the VLT, as part of the Fornax3D project. By fitting recently-developed population-orbital superposition models to the observed surface brightness as well as stellar kinematic, age, and metallicity maps, we obtain the stellar orbits, age and metallicity distributions of each galaxy. We then decompose each galaxy into multiple orbital-based components, including a dynamically hot inner stellar halo component which is identified as the relic of past massive mergers. By comparing to analogues from cosmological galaxy simulations, chiefly TNG50, we find that the formation of such a hot inner stellar halo requires the merger with a now-destroyed massive satellite galaxy of $3.7_{-1.5}^{+2.7} \times 10^{10}$ Msun (about $1/5$ of its current stellar mass) in the case of NGC 1380 and of $1.5_{-0.7}^{+1.6} \times10^{10}$ Msun (about $1/4$ of its current stellar mass) in the case of NGC 1427. Moreover, we infer that the last massive merger in NGC 1380 happened $\sim10$ Gyr ago based on the stellar age distribution of the re-grown dynamically cold disk, whereas the merger in NGC 1427 ended $t\lesssim 8$ Gyr ago based on the stellar populations in its hot inner stellar halo. The major merger event in NGC 1380 is the first one with both merger mass and merger time quantitatively inferred in a galaxy beyond the Local Volume. Moreover, it is the oldest and most massive merger uncovered in nearby galaxies so far.

Recent studies of high-redshift galaxies using JWST, such as GN-z11 revealed highly elevated levels of nitrogen (N). This phenomenon extends to gravitationally-lensed galaxies like the Sunburst Arc at z = 2.37, as well as to globular clusters (GCs). We propose that this originates from the presence of very massive stars (VMSs) with masses ranging from 100 to 1000\,\Msun. The He {\sc ii} observed in the Sunburst Arc could also stem from the disproportionately large contribution of VMSs. We build an entirely new Framework for massive star evolution which is no longer set by Dutch or other mass-loss "recipes" but which take the physics of $\Gamma$ or $L/M$-dependent winds into account. We discuss the mass-loss kink and the transition mass-loss rate between optically thin and thick winds, before we study the evaporative mass-loss history of VMSs. Our novel evolution models exhibit vertical evolution in the HR-diagram from the zero-age main sequence due to a self-regulatory effect driven by their wind-dominated nature, and we discuss what wind physics sets the stellar upper-mass limit. Our estimate for the Sunburst Arc in Vink (2023) suggests that the significant amounts of N found in star-forming galaxies likely arise from VMSs. We evaluate the strengths and weaknesses of previous hypotheses, including fast rotating massive stars and supermassive stars (SMSs), and we conclude that only our VMS model satisfies the relevant criteria. Finally, we advocate for the inclusion of VMSs in population synthesis and chemical evolution models, emphasizing the need for a self-consistent wind approach, which currently does not exist. Even minor inaccuracies in mass-loss rates dramatically impact the stellar evolution of VMS, as well as their ionizing and chemical feedback.

We predict quantitative mass-loss rates and terminal wind velocities for early-type supergiants and luminous blue variables (LBVs) using a dynamical version of the Monte Carlo radiative transfer method. First, the observed drop in terminal wind velocity around spectral type B1 is confirmed by the Monte Carlo method -- at the correct effective temperature of about 21 000 K. This drop in wind velocity is much steeper than would be expected from the drop in escape speed for cooler stars. The results may be particularly relevant for slow winds inferred for some High-Mass X-ray binaries. Second, the strength of the mass-loss bi-stability jump is found to be significantly larger than previously assumed. Not only could this make bi-stability braking more efficient in massive star evolution, but a rotationally-induced version of the bi-stability mechanism may now be capable of producing the correct density of outflowing disks around B[e] supergiants, although multi-dimensional modelling including the disk velocity structure is still needed. For LBVs, we find the bi-stability jump to become larger at higher metallicities, but perhaps surprisingly also larger at lower Eddington parameters. This may have consequences for the role of LBVs in the evolution of massive stars at different metallicities and Cosmic Epochs. Finally, our predicted low wind velocities may be important for explaining the slow outflow speeds of supernova type IIb/IIn progenitors, for which the direct LBV-SN link was first introduced.

Comprehensive spectral analyses of the Galactic Wolf-Rayet stars of the nitrogen sequence (i.e.\ the WN subclass) have been performed in a previous paper. However, the distances of these objects were poorly known. Distances have a direct impact on the "absolute" parameters, such as luminosities and mass-loss rates. The recent Gaia Data Release (DR2) of trigonometric parallaxes includes nearly all WN stars of our Galactic sample. In the present paper, we apply the new distances to the previously analyzed Galactic WN stars and rescale the results accordingly. On this basis, we present a revised catalog of 55 Galactic WN stars with their stellar and wind parameters. The correlations between mass-loss rate and luminosity show a large scatter, for the hydrogen-free WN stars as well as for those with detectable hydrogen. The slopes of the $\log L - \log \dot{M}$ correlations are shallower than found previously. The empirical Hertzsprung-Russell diagram (HRD) still shows the previously established dichotomy between the hydrogen-free early WN subtypes that are located on the hot side of the zero-age main sequence (ZAMS), and the late WN subtypes, which show hydrogen and reside mostly at cooler temperatures than the ZAMS (with few exceptions). However, with the new distances, the distribution of stellar luminosities became more continuous than obtained previously. The hydrogen-showing stars of late WN subtype are still found to be typically more luminous than the hydrogen-free early subtypes, but there is a range of luminosities where both subclasses overlap. The empirical HRD of the Galactic single WN stars is compared with recent evolutionary tracks. Neither these single-star evolutionary models nor binary scenarios can provide a fully satisfactory explanation for the parameters of these objects and their location in the HRD.

Interstellar complex organic molecules (COMs) in solar-like young stellar objects (YSOs), particularly within protostellar disks, are of significant interest due to their potential connection to prebiotic chemistry in emerging planetary systems. We report the discovery of a rotating feature enriched in COMs, including CH3OH, CH3CHO, and NH2CHO, in the protostellar core G192.12-11.10. By constructing a YSO model, we find that the COM-rich feature is likely located within or near the boundary of the Keplerian disk. The image synthesis results suggest that additional heating mechanisms leading to a warm ring or a warm inner disk are required to reproduce the observed emission. We discuss possible origins of the COM-rich feature, particularly accretion shocks as a plausible cause for a warm ring. Additionally, molecules such as C18O, H2CO, DCS, H2S, and OCS exhibit distinct behavior compared to CH3OH, indicating a range of physical and chemical conditions within the region. The observed kinematics of H2S and OCS suggest that OCS resides in regions closer to the central protostar than H2S, consistent with previous experimental studies.

The study of resonant oscillation modes in low-mass red giant branch stars enables their ages to be inferred with exceptional ($\sim$10%) precision, unlocking the possibility to reconstruct the temporal evolution of the Milky Way at early cosmic times. Ensuring the accuracy of such a precise age scale is a fundamental yet difficult challenge. Since the age of red giant branch stars primarily hinges on their mass, an independent mass determination for an oscillating red giant star provides the means for such assessment. We analyze the old eclipsing binary KIC10001167, which hosts an oscillating red giant branch star and is a member of the thick disk of the Milky Way. Of the known red giants in eclipsing binaries, this is the only member of the thick disk that has asteroseismic signal of high enough quality to test the seismic mass inference at the 2% level. We measure the binary orbit and obtain fundamental stellar parameters through combined analysis of light curve eclipses and radial velocities, and perform a detailed asteroseismic, photospheric, and Galactic kinematic characterization of the red giant and binary system. We show that the dynamically determined mass $0.9337\pm0.0077 \rm\ M_{\odot}$ (0.8%) of this 10 Gyr-old star agrees within 1.4% with the mass inferred from detailed modelling of individual pulsation mode frequencies (1.6%). This is now the only thick disk stellar system, hosting a red giant, where the mass has been determined both asteroseismically with better than 2% precision, and through a model-independent method at 1% precision, and we hereby affirm the potential of asteroseismology to define an accurate age scale for ancient stars to trace the Milky Way assembly history.

The proximity & low dust extinction of the Magellanic Clouds provides an ideal environment for metal-poor massive stars to be studied in detail. The HST ULLYSES initiative has provided exquisite ultraviolet spectroscopy of a large sample of OB stars in the Magellanic Clouds, and its legacy value has been enhanced through the acquisition of optical VLT/XShooter spectroscopy (XShootU). We aim to determine the spectral types and physical properties of 122 LMC & 103 SMC OB stars observed via XShootU. Physical parameters are required for these to serve as templates in spectral libraries of metal poor massive stars. We also aim to identify double-lined binaries and OeBe stars for which analysis requires non-standard treatment. We apply a pipeline designed to analyse large spectroscopic samples of hot luminous stars to XShootU spectroscopic datasets, together with grids of synthetic model spectra computed with the non-LTE atmospheric code fastwind at LMC & SMC metallicities. We determine physical and wind properties of 97 LMC & 77 SMC massive stars, ranging from O2 to B9 subtypes, representing the majority of the XShootU OB sample (OeBe & candidate SB2 systems are excluded). Results are broadly in agreement with previous optical spectroscopic studies, with evolutionary masses spanning 12-117 Msun in the LMC & 11-74 Msun in the SMC. We determine a revised Teff-spectral type calibration for Magellanic Cloud stars, identify stars with peculiar radial velocities, and compare wind properties of high luminosity O stars with dense winds, revealing ~0.27 dex higher wind momenta of LMC stars with respect to SMC counterparts. Incorporating the recent empirical metallicity dependence of Z^0.22 for wind velocities, this suggests a mass-loss dependence of Z^0.5 for luminous O stars. Studies incorporating ultraviolet mass-loss diagnostics are required for OB stars with weak winds and/or low luminosities.

The sharpest optical images of the R136 cluster in the Large Magellanic Cloud are presented, allowing for the first time to resolve members of the central core, including R136a1, the most massive star known. These data were taken using the Gemini speckle imager Zorro in medium-band filters with effective wavelengths similar to BVRI achieving angular resolutions between 30-40 mas. All stars previously known in the literature, having V&lt;16 mag within the central $2''\times2''$ were recovered. Visual companions ($\geq$40 mas; 2000 au) were detected for the WN5h stars R136 a1, and a3. Photometry of the visual companion of a1 suggests it is of mid O spectral type. Based on new photometric luminosities using the resolved Zorro imaging, the masses of the individual WN5h stars are estimated to be between 150-200 $M_{\odot}$, lowering significantly the present-day masses of some of the most massive stars known. These mass estimates are critical anchor points for establishing the stellar upper-mass function.

Using Gaia DR3 we derive new distances and luminosities for a sample of Galactic B supergiants which were thought to be post main-sequence (MS) objects from their HR diagram location beyond the terminal-age MS (TAMS). When applying the newer Gaia distances in addition to enhanced amounts of core-boundary mixing, aka convective overshooting, we show that these Galactic B supergiants are likely enclosed within the MS band, indicating an evolutionary stage of steady core hydrogen burning. We discuss the importance of considering enhanced overshooting and how vectors in the mass-luminosity plane (ML-plane) can be used to disentangle the effects of wind mass loss from interior mixing. We finish with the key message that any proposed solution to the BSG problem should consider not only an explanation for the sheer number of B supergiants inside the Hertzsprung gap, but should at the same time also account for the steep drop in rotation rates identified at spectral type B1 -- corresponding to an effective temperature of $\sim$21 kK, and for which two distinct families of solutions have been proposed.

We report the discovery of the ``mm fundamental plane of black-hole accretion'', which is a tight correlation between the nuclear 1 mm luminosity ($L_{\rm \nu, mm}$), the intrinsic $2$ -- $10$~keV X-ray luminosity ($L_{\rm X,2-10}$) and the supermassive black hole (SMBH) mass ($M_{\rm BH}$) with an intrinsic scatter ($\sigma_{\rm int}$) of $0.40$ dex. The plane is found for a sample of 48 nearby galaxies, most of which are low-luminosity active galactic nuclei (LLAGN). Combining these sources with a sample of high-luminosity (quasar-like) nearby AGN, we find that the plane still holds. We also find that $M_{\rm BH}$ correlates with $L_{\rm \nu, mm}$ at a highly significant level, although such correlation is less tight than the mm fundamental plane ($\sigma_{\rm int}=0.51$ dex). Crucially, we show that spectral energy distribution (SED) models for both advection-dominated accretion flows (ADAFs) and compact jets can explain the existence of these relations, which are not reproduced by the standard torus-thin accretion disc models usually associated to quasar-like AGN. The ADAF models reproduces the observed relations somewhat better than those for compact jets, although neither provides a perfect prediction. Our findings thus suggest that radiatively-inefficient accretion processes such as those in ADAFs or compact (and thus possibly young) jets may play a key role in both low- and high-luminosity AGN. This mm fundamental plane also offers a new, rapid method to (indirectly) estimate SMBH masses.