Search

We aim to trace the evolution of eight edge-on star-forming disc galaxies through the analysis of stellar population properties of their thin and thick discs. These galaxies have relatively low stellar masses (4 $\times$ 10$^9$ to 6 $\times$ 10$^{10}$ $M_{\odot}$). We use Multi-Unit Spectroscopic Explorer (MUSE) observations and full-spectrum fitting to produce spatially resolved maps of ages, metallicities and [Mg/Fe]-abundances and extract the star formation histories of stellar discs. Our maps show thick discs that are on average older, more metal-poor and more [Mg/Fe]-enhanced than thin discs. However, age differences between thin and thick discs are small (around 2 Gyr) and the thick discs are younger than previously observed in more massive and more quiescent galaxies. Both thin and thick discs show mostly sub-solar metallicities, and the vertical metallicity gradient is milder than previously observed in similar studies. [Mg/Fe] differences between thick and thin discs are not sharp. The star formation histories of thick discs extend down to recent times, although most of the mass in young stars was formed in thin discs. Our findings show thick discs that are different from old quiescent thick discs previously observed in galaxies of different morphologies and/or different masses. We propose that thick discs in these galaxies did not form quickly at high redshift, but slowly in an extended time. Also, the thin discs formed slowly, however, a larger mass fraction was created at very recent times., Comment: 20 pages, 9 figures, submitted to A&A

Understanding how gas flows into galactic centres, fuels the AGN, and is in turn expelled back through feedback processes is of great importance to appreciate the role AGN play in the growth and evolution of galaxies. We use the MUSE-AO optical spectra of the inner 7".5x 7".5 of the nearby Seyfert 1 galaxy NGC 4593 to characterise its ionised gas kinematics. We fit single-Gaussian components to the [O iii]$\lambda$ 5007 and [N ii]$\lambda$ 6583 emission lines, and double-Gaussian components to H$\alpha$ and H$\beta$ to determine the main ionisation mechanism of the gas. To determine the kinematics of the ionised gas, we fit double-Gaussian components to the [O iii]$\lambda$ 5007 line. Based on the stellar kinematic maps, we confirm the presence of a rotating disc. For the ionised gas, we find high-velocity dispersion values of up to $\mathrm{200-250\,km\,s^{-1}}$ that show that part of the gas is highly perturbed. The dominant ionisation mechanism of the gas is AGN photoionisation, which reaches the highest values within the innermost 4arcsec (680pc) diameter of the galaxy. At larger radii, the emission line ratios correspond to values in the composite region of the BPT diagram. The broad-component of [O iii]$\lambda$ 5007 shows blue-shifted velocities on the east side of the central 2arcsec (340pc), which spatially coincide with a region of high velocity-dispersion. This confirms the presence of outflowing gas. We estimate a mass outflow rate and kinetic power of ${\dot{M}}\geq {0.048 \,M_\odot\,{{\rm yr}^{-1}} }$ and ${\dot{E}_{\text{kin}}} \geq 4.09 \times 10^{39} \, \text{erg} \, \text{s}^{-1}$. The derived mass outflow rate is consistent with that expected from empirical relations between mass outflow rate and AGN luminosity for a low-luminosity AGN such as NGC 4593., Comment: Accepted for publication in A&A

The ionization of diffuse gas located far above the energetic midplane OB stars poses a challenge to the commonly accepted notion that radiation from OB stars is the primary ionization source for gas in galaxies. We investigated the sources of ionizing radiation, specifically leaking midplane HII regions and/or in situ hot low-mass evolved stars (HOLMES), in extraplanar diffuse ionized gas (eDIG) in a sample of eight nearby (17-52 Mpc) edge-on disk galaxies observed with the Multi Unit Spectroscopic Explorer (MUSE). We constructed a model for the photoionization of eDIG clouds and the propagation of ionizing radiation through the eDIG using subsequent runs of Cloudy photoionization code. Our model includes radiation originating both from midplane OB stars and in situ evolved stars and its dilution and processing as it propagates in the eDIG. We fit the model to the data using the vertical line ratio profiles of our sample galaxies, and find that while the ionization by in situ evolved stars is insignificant for most of the galaxies in our sample, it may be able to explain the enhanced high-ionization lines in the eDIG of the green valley galaxy ESO 544-27. Our results show that while leaking radiation from midplane HII regions is the primary ionization source for eDIG, in situ evolved stars can play a significant part in ionizing extraplanar gas in galaxies with low star forming rates., Comment: 19 pages, 22 figures, accepted to Astronomy & Astrophysics

Highly inclined (edge-on) disk galaxies offer the unique perspective to constrain their intrinsic flattening, $c/a$, where $c$ and $a$ are respectively the vertical and long radial axes of the disk measured at suitable stellar densities. The ratio $c/a$ is a necessary quantity in the assessment of galaxy inclinations, three-dimensional structural reconstructions, total masses, as well as a constraint to galaxy formation models. 3.6 micron maps of 133 edge-on spiral galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G) and its early-type galaxy extension are used to revisit the assessment of $c/a$ free from dust extinction and away from the influence of a stellar bulge. We present a simple definition of $c/a$ and explore trends with other galactic physical parameters: total stellar mass, concentration index, total HI mass, mass of the central mass concentration, circular velocity, model-dependent scales, as well as Hubble type. Other than a dependence on early/late Hubble types, and a related trend with light concentration, no other parameters were found to correlate with the intrinsic flattening of spiral galaxies. The latter is mostly constant with $\langle c/a \rangle$ = 0.124 $\pm$ 0.001 (stat) $\pm$ 0.033 (intrinsic/systematic) and greater for earlier types., Comment: 27 full pages, 14 figures, 6 tables, 1 appendix. Submitted to AAS Journals. Comments welcome

Metallicity measurements in galaxies can give valuable clues about galaxy evolution. One of the mechanisms postulated for metallicity redistribution in galaxies is gas flows induced by AGN, but the details of this process remain elusive. We report the discovery of a positive radial gradient in the gas-phase metallicity of the narrow line region of the Seyfert 2 galaxy NGC 7130, which is not found when considering the star-forming components in the galaxy disk. To determine gas-phase metallicities for each kinematic component, we use both active galactic nuclei (AGN) and star-forming (SF) strong-line abundance relations, as well as BPT diagnostic diagrams. These relations involve sensitive strong emission lines, namely [OIII]5007, [NII]6584, H$\alpha$, H$\beta$, [SII]6716, and [SII]6731, observed with the adaptive-optics-assisted mode of the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT). The presence of a positive radial metallicity gradient only in the AGN ionized component suggests that metals may be transported from central areas to its purlieus by AGN activity., Comment: Accepted for publication in Astronomy & Astrophysics (A&A)

We investigate the star formation occurring in the Planck Galactic cold clump PGCC 120.69+2.66. Near-infrared JHKs images and K-band spectroscopy obtained with NOTCam at the Nordic Optical Telescope complemented with archive data are used to study the stellar content. In addition, millimetre line CO and CS spectra were obtained with the Onsala 20 m telescope, and sub-millimetre continuum SCUBA archive data are used to characterise the host molecular cloud. We identify a molecular cloud core traced by CO and CS emission at a distance of 1.1 kpc. In the region studied, we identify 5 submm continuum cores. Embedded in and around these dense submm cores, we find 38 young stellar objects, classified as 9 Class I, 8 Class II, and 21 near-IR excess or variability sources, accompanied by bipolar nebulosities and signs of protostellar jets. Furthermore, a very bright and reddened source is found towards this molecular cloud core. Even though its location appears to suggest its association to the star formation region, its infrared spectral type is compatible with a red supergiant, hidden behind 36 mag of visual extinction., Comment: 12 pages, 9 figures

Star formation estimates based on the counting of YSOs is commonly applied to nearby star-forming regions in the Galaxy. With this method, the SFRs are measured using the counts of YSOs in a particular protostellar Class, a typical protostellar mass, and the lifetime associated with this Class. However, the assumptions underlying the validity of the method such as that of a constant star formation history (SFH) and whether the method is valid for all protostellar Classes has never been fully tested. In this work, we use Monte Carlo models to test the validity of the method. We build synthetic clusters in which stars form at times that are randomly drawn from a specified SFH. The latter is either constant or time-dependent with a burst like behavior. The masses of the protostars are randomly drawn from an IMF which can be either similar to that of the Milky Way field or be variable . For each star in every cluster, the lifetimes associated with the different protostellar classes are also randomly drawn from Gaussian distribution functions centered around their most likely value as suggested by the observations. We find that only the SFR derived using the Class 0 population can reproduce the true SFR at all epochs, and this is true irrespective of the shape of the SFH. For a constant SFH, the SFR derived using the more evolved populations of protostars (Classes I, F, II, and III) reproduce the real SFR only at later epochs which correspond to epochs at which their numbers have reached a steady state. For a time-dependent burst-like SFH, all SFR estimates based on the number counts of the evolved populations fail to reproduce the true SFR. We also show how the offsets between Class I and Class II based SFRs and the true SFR plotted as a function of the number ratios of Class I and Class II versus Class III YSOs can be used in order to constrain the SFH of observed molecular clouds., Comment: Submitted. Comments are welcome

High-latitude intermediate-velocity clouds (IVCs) are part of the Milky Way's HI halo and originate from either a galactic fountain process or extragalactic gas infall. They are partly molecular and can most of the time be identified in CO. Some of these regions also exhibit high-velocity cloud (HVC) gas, which is mostly atomic, and gas at local velocities (LVCs), which is partly atomic and partly molecular. We conducted a study on the IVCs Draco and Spider, both were exposed to a very weak UV field, using the receiver upGREAT on SOFIA. The 158 micron line of ionized carbon (CII) was observed, and the results are as follows: In Draco, the CII line was detected at intermediate velocities (but not at local or high velocities) in four out of five positions. No CII emission was found at any velocity in the two observed positions in Spider. To understand the excitation conditions of the gas in Draco, we analyzed complementary CO and HI data as well as dust column density and temperature maps from Herschel. The observed CII intensities suggest the presence of shocks in Draco that heat the gas and subsequently emit in the CII cooling line. These shocks are likely caused by the fast cloud's motion toward the Galactic plane that is accompanied by collisions between HI clouds. The nondetection of CII in the Spider IVC and LVC as well as in other low-density clouds at local velocities that we present in this paper (Polaris and Musca) supports the idea that highly dynamic processes are necessary for CII excitation in UV-faint low-density regions., Comment: in press Astronomy and Astrophysics

The universality or non-universality of the initial mass function (IMF) has significant implications for determining star formation rates and star formation histories from photometric properties of stellar populations. We reexamine whether the IMF is deficient in high-mass stars (top-light) in the low-density environment of the outer disk of M83 and constrain the shape of the IMF therein. Using archival Galaxy Evolution Explorer (GALEX) far ultraviolet (FUV) and near ultraviolet (NUV) data and new deep OmegaCAM narrowband H$\alpha$ imaging, we constructed a catalog of FUV-selected objects in the outer disk of M83. We counted H$\alpha$-bright clusters and clusters that are blue in FUV$-$NUV in the catalog, measured the maximum flux ratio $F_{\mathrm{H}\alpha}/f_{\lambda \mathrm{FUV}}$ among the clusters, and measured the total flux ratio $\Sigma F_{\mathrm{H}\alpha}/\Sigma f_{\lambda \mathrm{FUV}}$ over the catalog. We then compared these measurements to predictions from stellar population synthesis models made with a standard Salpeter IMF, truncated IMFs, and steep IMFs. We also investigated the effect of varying the assumed internal extinction on our results. We are not able to reproduce our observations with models using the standard Salpeter IMF or the truncated IMFs. It is only when assuming an average internal extinction of $0.10 < A_{\mathrm{V}} < 0.15$ in the outer disk stellar clusters that models with steep IMFs ($\alpha > 3.1$) simultaneously reproduce the observed cluster counts, the maximum observed $F_{\mathrm{H}\alpha}/f_{\lambda \mathrm{FUV}}$, and the observed $\Sigma F_{\mathrm{H}\alpha}/\Sigma f_{\lambda \mathrm{FUV}}$. Our results support a non-universal IMF that is deficient in high-mass stars in low-density environments., Comment: 18 pages, 15 figures, accepted to Astronomy & Astrophysics

We show that the rotation curves of 16 nearby disc galaxies in the THINGS sample and the Milky Way can be described by the NFW halo model and by the Bosma effect at approximately the same level of accuracy. The latter effect suggests that the behavior of the rotation curve at large radii is determined by the rescaled gas component and thus that dark matter and gas distributions are tightly correlated. By focusing on galaxies with exponential decay in their gas surface density, we can normalize their rotation curves to match the exponential thin disc model at large enough radii. This normalization assumes that the galaxy mass is estimated consistently within this model, assuming a thin disc structure. We show that this rescaling allows us to derive a new version of the Tully-Fisher (TF) relation, the Bosma TF relation that nicely fit the data. In the framework of this model, the connection between the Bosma Tully-Fisher (TF) relation and the baryonic TF relation can be established by considering an additional empirical relation between the baryonic mass and the total mass of the disc, as measured in the data., Comment: 21 pages, 12 figures. Monthly Notices of the Royal Astronomical Society, in the press

A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yield an impressive dataset of ~120 billion sources and map the structure of our Galaxy. The effort would foster subsequent expansions in numerous dimensions (spatial, depth, wavelengths, epochs). Importantly, the survey would benefit from early defintion by the community, namely because the Galactic disk is a complex environment, and different science goals will require trade offs., Comment: Submitted to the Roman Project on October 22 2021 in response to a call for white papers on early-definition Astrophysics opportunity

In the standard cosmological model of galaxy evolution, mergers and interactions play a fundamental role in shaping galaxies. Galaxies that are currently isolated are thus interesting, allowing us to distinguish between internal and external processes affecting the galactic structure. However, current observational limits may obscure crucial information in the low-mass or low-brightness regime. We use optical imaging of a subsample of the AMIGA catalogue of isolated galaxies to explore the impact of different factors on the structure of these galaxies. We study the type of disc break as a function of the degree of isolation and the presence of interaction indicators like tidal streams or plumes only detectable in the low surface brightness regime. We present deep optical imaging of a sample of 25 isolated galaxies. Through careful data processing and analysis techniques, the surface brightness limits achieved are comparable to those to be obtained on the 10-year LSST coadds. The extreme depth of our imaging allows us to study the interaction signatures of 20 galaxies, given that the presence of Galactic cirrus is a strong limiting factor in the characterisation of interactions for the remaining 5 of them. We detect previously unreported interaction features in 8 (40%) galaxies in our sample. We identify 9 galaxies (36%) showing an exponential disc (Type I), 14 galaxies (56%) with down-bending (Type II) profile and only 2 galaxies (8%) with up-bending (Type III) profiles. Isolated galaxies have considerably more purely exponential discs and fewer up-bending surface brightness profiles than field or cluster galaxies. We suggest that major mergers produce up-bending profiles while a threshold in star formation probably forms down-bending profiles. Unperturbed galaxies, evolving slowly with a low star formation rate could cause the high rate of Type I discs in isolated galaxies observed., Comment: 20 pages, 7 figures, and 2 tables (incl. 2 appendix). Accepted for publication in A&A. Images and profiles will be made publicly available upon publication

Context: The hierarchical model of galaxy formation suggests that galaxies are continuously growing. However, our position inside the Milky Way prevents us from studying the disk edge. Truncations are low surface brightness features located in the disk outskirts of external galaxies. They indicate where the disk brightness abruptly drops and their location is thought to change dynamically. In previous analyses of Milky Way-like galaxies, truncations were detected up to 3 kpc above the mid-plane but whether they remain present beyond that height remains unclear. Aims: Our goal is to determine whether truncations can be detected above 3 kpc height in the Milky Way-like galaxy NGC 4565, thus establishing the actual disk thickness. We also aim to study how the truncation relates to disk properties such as star formation activity or the warp. Methods: We perform a vertical study of the disk of NGC 4565 edge in unprecedented detail. We explore the truncation radius at different heights above/below the disk mid-plane (0

Disk galaxies are typically in a stable configuration where matter moves in almost closed circular orbits. However, non-circular motions caused by distortions, warps, lopsidedness, or satellite interactions are common and leave distinct signatures on galaxy velocity maps. We develop an algorithm that uses an ordinary least square method for fitting a non-axisymmetric model to the observed two-dimensional line-of-sight velocity map of an external galaxy, which allows for anisotropic non-circular motions. The method approximates a galaxy as a flat disk, which is an appropriate assumption for spiral galaxies within the optical radius where warps are rare. In the outer parts of HI distributions, which may extend well into the warp region, we use this method in combination with a standard rotating tilted ring model to constrain the range of radii where the flat disk assumption can be conservatively considered valid. Within this range, the transversal and radial velocity profiles, averaged in rings, can be directly reconstructed from the velocity map. The novelty of the algorithm consists in using arc segments in addition to rings: in this way spatial velocity anisotropies can be measured in both components, allowing for the reconstruction of angularly resolved coarse-grained two-dimensional velocity maps. We applied this algorithm to 25 disk galaxies from the THINGS sample for which we can provide 2D maps of both velocity components., Comment: 42 pages, 37 figures, Monthly Notices of the Royal Astronomical Society Main Journal in the press (2023)

According to the $\Lambda$CDM cosmology, present-day galaxies with stellar masses $M_\star>10^{11} {\rm M}_\odot$ should contain a sizable fraction of dark matter within their stellar body. Models indicate that in massive early-type galaxies (ETGs) dark matter should account for $\sim60\%$ of the dynamical mass within five effective radii ($5 R_{\rm e}$). Most massive ETGs have been shaped through a two-phase process: the rapid growth of a compact core was followed by the accretion of an extended envelope through mergers. The exceedingly rare galaxies that have avoided the second phase, the so-called relic galaxies, are thought to be the frozen remains of the massive ETG population at $z\gtrsim2$. The best relic galaxy candidate discovered to date is NGC 1277, in the Perseus cluster. We used deep integral field GCMS data to revisit NGC 1277 out to an unprecedented radius of 6 kpc (corresponding to $5 R_{\rm e}$). By using Jeans anisotropic modelling we find a negligible dark matter fraction within $5 R_{\rm e}$ ($f_{\rm DM}(5 R_{\rm e})<0.05$; two-sigma confidence level), which is in tension with the expectation. Since the lack of an extended envelope would reduce dynamical friction and prevent the accretion of an envelope, we propose that NGC 1277 lost its dark matter very early or that it was dark matter deficient ab initio. We discuss our discovery in the framework of recent proposals suggesting that some relic galaxies may result from dark matter stripping as they fell in and interacted within galaxy clusters. Alternatively, NGC 1277 might have been born in a high-velocity collision of gas-rich proto-galactic fragments, where dark matter left behind a disc of dissipative baryons. We speculate that the relative velocities of $\approx2000 {\rm km/s}$ required for the latter process to happen were possible in the progenitors of the present-day rich galaxy clusters., Comment: Accepted for publication in A&A

The Gaia-ESO survey sample of massive OB stars in the Carina Nebula consists of 234 stars. The addition of brighter sources from the Galactic O-Star Spectroscopic Survey and additional sources from the literature allows us to create the most complete census of massive OB stars done so far in the region. It contains a total of 316 stars, being 18 of them in the background and four in the foreground. Of the 294 stellar systems in Car OB1, 74 are of O type, 214 are of non-supergiant B type and 6 are of WR or non-O supergiant (II to Ia) spectral class. We identify 20 spectroscopic binary systems with an O-star primary, of which 6 are reported for the first time, and another 18 with a B-star primary, of which 13 are new detections. The average observed double-lined binary fraction of O-type stars in the surveyed region is 0.35, which represents a lower limit. We find a good correlation between the spectroscopic n-qualifier and the projected rotational velocity of the stars. The fraction of candidate runaways among the stars with and without the n-qualifier is 4.4% and 2.4%, respectively, although non resolved double-lined binaries can be contaminating the fast rotators sample., Comment: Accepted for publication in A&A, 62 pages, 8 figures

The HH 24 complex harbors five collimated jets emanating from a small protostellar multiple system. We have carried out a multi-wavelength study of the jets, their driving sources, and the cloud core hosting the embedded stellar system, based on data from the HST, Gemini, Subaru, APO 3.5m, VLA, and ALMA telescopes. The data show that the multiple system, SSV 63, contains at least 7 sources, ranging in mass from the hydrogen-burning limit to proto-Herbig Ae stars. The stars are in an unstable non-hierarchical configuration, and one member, a borderline brown dwarf, is moving away from the protostellar system with 25 km/s, after being ejected about 5,800 yr ago as an orphaned protostar. Five of the embedded sources are surrounded by small, possibly truncated, disks resolved at 1.3 mm with ALMA. Proper motions and radial velocities imply jet speeds of 200-300 km/s. The two main HH 24 jets, E and C, form a bipolar jet system which traces the innermost portions of parsec-scale chains of Herbig-Haro and H2 shocks with a total extent of at least 3 parsec. H2CO and C18O observations show that the core has been churned and continuously fed by an infalling streamer. 13CO and 12CO trace compact, low-velocity, cavity walls carved by the jets and an ultra-compact molecular outflow from the most embedded object. Chaotic N-body dynamics likely will eject several more of these objects. The ejection of stars from their feeding zones sets their masses. Dynamical decay of non-hierarchical systems can thus be a major contributor to establishing the initial mass function., Comment: 57 pages, 61 figures, 12 tables; accepted to Astron.J

Thin and thick disks are found in most spiral galaxies, yet their formation scenarios remain uncertain. Whether thick disks form through slow or fast, internal or environmental, processes is unclear. The physical origin of outer truncations in thin and thick disks, observed as a drop in optical and near-infrared (NIR) surface brightness profiles, is also a much debated topic. These truncations have been linked to star formation (SF) thresholds in Milky-Way type galaxies, but no such connection has been made for their low-mass counterparts or in thick disks. Our photometric analysis of the edge-on galaxy UGC 7321 offers a possible breakthrough. This well-studied diffuse, isolated, bulgeless, ultra-thin galaxy is thought to be under-evolved both dynamically and in SF. It is an ideal target to disentangle internal effects in the formation of thick disks and truncations. Our axial light profiles from deep far- and near-ultraviolet (UV; GALEX) images, tracing recent SF, and optical (DESI grz) and NIR (Spitzer 3.6 microns) images, tracing old stellar populations, enable a detailed identification of an outer truncation in all probed wavelengths in both the thin and thick disks. After deprojecting to a face-on view, a sharp truncation signature is found at a stellar density of roughly 1.5 solar masses per square parsec, in agreement with theoretical expectations of gas density SF thresholds. The redder colours beyond the truncation radius are indicative of stellar migration towards the outer regions. We thus show that thick disks and truncations can form via internal mechanisms alone, given the pristine nature of UGC 7321. We report the discovery of a truncation at and above the mid-plane of a diffuse galaxy that is linked to a SF threshold; this poses a constraint on physically-motivated disk size measurements among low-mass galaxies., Comment: Accepted for publication in A&A (August 29, 2022). 7 pages, 3 figures

We investigate the young stellar populations associated with DR 18 and ECX 6-21, which are two isolated globules irradiated by the O-type stars of the Cygnus OB2 association. Both are HII regions containing obvious tracers of recent and ongoing star formation. We also study smaller isolated molecular structures in their surroundings. Both globules contain their own embedded populations, with a higher fraction of the less-evolved classes. Masses and temperatures are estimated under the assumption of a common age of 1 Myr, which has been found to appropriately represent the general Cygnus OB2 YSO population but is most probably an overestimate for both globules, especially ECX 6-21. The early-B star responsible for the erosion of DR 18 is found to be part of a small aggregate of intermediate-mass stars still embedded in the cloud, which probably contains a second site of recent star formation, also with intermediate-mass stars. We confirm the two main star forming sites embedded in ECX 6-21 described in previous works, with the southern site being more evolved than the northern site. We also discuss the small globule ECX 6-21-W ($= G79.8+1.2$), and propose that its non thermal radio spectrum is due to synchrotron emission from an embedded jet, whose existence is suggested by our observations. The extreme youth of some of the YSOs suggests that star formation in both globules started after they became externally irradiated. The populations of both globules are not found to be particularly rich, but they contain stars with estimated masses similar or above that of the Sun in numbers that hint at some differences with respect to the star formation process taking place in more quiescent regions where low-mass stars dominate, which deeper observations may confirm., Comment: Accepted for publication by Astronomy and Astrophysics

W40 is a heavily obscured bipolar HII region projected in the direction of the Aquila Rift and ionized by hot stars in a central, partly embedded cluster. The study of the cluster and its surroundings has been greatly hampered thus far by the strong extinction in the region. We use the Gaia eDR3 catalog to establish astrometric membership criteria based on the population of the W40 central cluster, reassess the distance of the region, and identify in this way new members, both inside and outside the cluster. We obtain visible spectroscopy in the red spectral region to classify both known and new members, complemented with Gaia and Spitzer photometry to assess the evolutionary status of the stellar population. We derive a high-confidence geometric distance to the W40 region of 502 pc $\pm$ 4 pc and confirm the presence of a comoving extended population of stars at the same distance, spreading over the whole projected area of the HII region and beyond. Spectral classifications are presented for 21 members of the W40 region, 10 of them belonging to the central cluster. One of the newly identified B stars in the extended population is clearly interacting with the shell surrounding the HII region, giving rise to a small arc-shaped nebula that traces a bow shock. The infrared excess properties suggest that the extended population is significantly older ($\sim 3$ Myr) than the W40 central cluster ($< 1$ Myr). The area currently occupied by the W40 HII region and its surroundings has a history of star formation extending at least several million years in the past, of which the formation of the W40 central cluster and the subsequent HII region is one of the latest episodes. The newly determined distance suggests that W40 is behind, and physically detached from, a pervasive large dust layer which is some 60 pc foreground to it as determined by previous studies., Comment: Submitted to Astronomy and Astrophysics

2MASS J20395358+4222505 is an obscured early B supergiant near the massive OB star association Cyg OB2. Despite its bright infrared magnitude (K$_{s}$=5.82) it has remained largely ignored because of its dim optical magnitude (B=16.63, V=13.68). In a previous paper we classified it as a highly reddened, potentially extremely luminous, early B-type supergiant. We obtained its spectrum in the U, B and R spectral bands during commissioning observations with the instrument MEGARA@GTC. It displays a particularly strong H$\alpha$ emission for its spectral type, B1 Ia. The star seems to be in an intermediate phase between super- and hypergiant, a group that it will probably join in the near (astronomical) future. We observe a radial velocity difference between individual observations and determine the stellar parameters, obtaining T$_{eff}$ = 24000 K, logg$_{c}$= 2.88 $\pm$ 0.15. The rotational velocity found is large for a B-supergiant, vsini= 110 $\pm$ 25 km s$^{-1}$. The abundance pattern is consistent with solar, with a mild C underabundance (based on a single line). Assuming that J20395358+4222505 is at the distance of Cyg OB2 we derive the radius from infrared photometry, finding R= 41.2 $\pm$ 4.0 R$_{\odot}$, log(L/L$_{\odot}$)= 5.71 $\pm$ 0.04 and a spectroscopic mass of 46.5 $\pm$ 15.0 M$_{\odot}$. The clumped mass-loss rate (clumping factor 10) is very high for the spectral type, $\dot{M}$ = 2.4x10$^{-6}$ M$_{\odot}$ a$^{-1}$. The high rotational velocity and mass-loss rate place the star at the hot side of the bi-stability jump. Together with the nearly solar CNO abundance pattern, they may also point to evolution in a binary system, J20395358+4222505 being the initial secondary., Comment: Accepted for publication in MNRAS, 13 pages, 9 figures

The Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G) is a detailed study of over 2300 nearby galaxies in the near-infrared (NIR), which has been critical to our understanding of the detailed structures of nearby galaxies. Because the sample galaxies were selected only using radio-derived velocities, however, the survey favored late-type disk galaxies over lenticulars and ellipticals. A follow-up Spitzer survey was conducted to rectify this bias, adding 465 early-type galaxies (ETGs) to the original sample, to be analyzed in a manner consistent with the initial survey. We present the data release of this ETG extension, up to the third data processing pipeline (P3): surface photometry. We produce curves of growth and radial surface brightness profiles (with and without inclination corrections) using reduced and masked Spitzer IRAC 3.6$\mu$m and 4.5$\mu$m images produced through Pipelines 1 and 2, respectively. From these profiles, we derive the following integrated quantities: total magnitudes, stellar masses, concentration parameters, and galaxy size metrics. We showcase NIR scaling relations for ETGs among these quantities. We examine general trends across the whole S$^{4}$G and ETG extension among our derived parameters, highlighting differences between ETGs and late-type galaxies (LTGs). ETGs are, on average, more massive and more concentrated than LTGs, and also show subtle distinctions among ETG morphological sub-types. We also derive the following scaling relations and compare with previous results in visible light: mass--size (both half-light and isophotal), mass--concentration, mass--surface brightness (central, effective, and within 1 kpc), and mass--color. We find good agreement with previous works, though some relations (e.g., mass--central surface brightness) will require more careful multi-component decompositions to be fully understood., Comment: 25 pages, 17 figures, accepted for publication in A&A

We investigate the origin of the extraplanar diffuse ionized gas (eDIG) and its predominant ionization mechanisms in five nearby (17-46 Mpc) low-mass ($10^9\text{-}10^{10}$ $M_{\odot}$) edge-on disk galaxies: ESO 157-49, ESO 469-15, ESO 544-27, IC 217, and IC 1553. We acquired Multi Unit Spectroscopic Explorer (MUSE) integral field spectroscopy and deep narrowband H$\alpha$ imaging of our sample galaxies. To investigate the connection between in-plane star formation and eDIG, we perform a photometric analysis of our narrowband H$\alpha$ imaging. We measure eDIG scale heights of $h_{z\text{eDIG}} = 0.59 \text{-} 1.39$ kpc and find a positive correlation between them and specific star formation rates. In all galaxies, we also find a strong correlation between extraplanar and midplane radial H$\alpha$ profiles. Using our MUSE data, we investigate the origin of eDIG via kinematics. We find ionized gas rotation velocity lags above the midplane with values between 10 and 27 km s$^{-1}$ kpc$^{-1}$. While we do find hints of an accretion origin for the ionized gas in ESO 157-49, IC 217, and IC 1553, overall the ionized gas kinematics of our galaxies do not match a steady galaxy model or any simplistic model of accretion or internal origin for the gas. We also construct standard diagnostic diagrams and emission-line maps (EW(H$\alpha$), [NII]/H$\alpha$, [SII]//H$\alpha$, [OIII]/H$\beta$) and find regions consistent with mixed OB star and hot low-mass evolved stars (HOLMES) ionization, and mixed OB-shock ionization. Our results suggest that OB stars are the primary driver of eDIG ionization, while both HOLMES and shocks may locally contribute to the ionization of eDIG to a significant degree. Despite our galaxies' similar structures and masses, we find a surprisingly composite image of ionization mechanisms and a multifarious origin for the eDIG., Comment: 21 pages, 14 figures, accepted to Astronomy & Astrophysics

IRAS 20319+3958 in Cygnus X South is a rare example of a free-floating globule (mass ~240 Msun, length ~1.5 pc) with an internal HII region created by the stellar feedback of embedded intermediate-mass stars, in particular, one Herbig Be star. Here, we present a Herschel/HIFI CII 158 mu map of the whole globule and a large set of other FIR lines (mid-to high-J CO lines observed with Herschel/PACS and SPIRE, the OI 63 mu line and the CO 16-15 line observed with upGREAT on SOFIA), covering the globule head and partly a position in the tail. The CII map revealed that the whole globule is probably rotating. Highly collimated, high-velocity CII emission is detected close to the Herbig Be star. We performed a PDR analysis using the KOSMA-tau PDR code for one position in the head and one in the tail. The observed FIR lines in the head can be reproduced with a two-component model: an extended, non-clumpy outer PDR shell and a clumpy, dense, and thin inner PDR layer, representing the interface between the HII region cavity and the external PDR. The modelled internal UV field of ~2500 Go is similar to what we obtained from the Herschel FIR fluxes, but lower than what we estimated from the census of the embedded stars. External illumination from the ~30 pc distant Cyg OB2 cluster, producing an UV field of ~150-600 G0 as an upper limit, is responsible for most of the CII emission. For the tail, we modelled the emission with a non-clumpy component, exposed to a UV-field of around 140 Go., Comment: accepted by Astronomy & Astrophysics

Context: Th 28 is a Classical T Tauri star in the Lupus 3 cloud which drives an extended bipolar jet. Previous studies of the inner jet identified signatures of rotation around the outflow axis, a key result for theories of jet launching. Thus this is an important source in which to investigate the poorly understood jet launching mechanism. We investigate the morphology and kinematics of the Th 28 micro-jets with the aim of characterizing their structure and outflow activity, using optical integral-field spectroscopy observations obtained with VLT/MUSE. We use spectro-imaging and position-velocity maps to investigate the kinematic and morphological features of the jet, and obtain a catalogue of emission lines in which the jet is visible. A Lucy-Richardson deconvolution procedure is used to differentiate the structure of the inner micro-jet region. Spatial profiles extracted perpendicular to the jet axis are fitted to investigate the jet width, opening angle and the evolution of the jet axis. We confirm the previously identified knot HHW$_{2}$ within the red-shifted jet and identify three additional knots in each lobe for the first time. We also find [O III]$\lambda$5007 emission from the blue-shifted micro-jet including the knot closest to the star. Proper motions for the innermost knots on each side are estimated and we show that new knots are ejected on an approximate timescale of 10-15 years. The jet axis centroids show a point-symmetric wiggle within the inner portion of both micro-jets indicating precession. We use the jet shape to measure a precession period of 8 years, with a half-opening angle < 0.6$^{\circ}$. This may provide an alternative explanation for the rotation signatures previously reported. We find the jet shape to be compatible with precession due to a brown dwarf companion orbiting at a separation $\leq$ 0.3 au., Comment: Accepted by Astronomy & Astrophysics; 21 pages, 23 figures

One of the suggested thick disc formation mechanisms is that they were born quickly and in situ from a turbulent clumpy disc. Subsequently, thin discs formed slowly within them from leftovers of the turbulent phase and from material accreted through cold flows and minor mergers. In this letter, I propose an observational test to verify this hypothesis. By combining thick disc and total stellar masses of edge-on galaxies with galaxy stellar mass functions calculated in the redshift range of $z\leq3.0$, I derived a positive correlation between the age of the youngest stars in thick discs and the stellar mass of the host galaxy; galaxies with a present-day stellar mass of $\mathcal{M}_\star(z=0)<10^{10}\,\mathcal{M}_\odot$ have thick disc stars as young as $4-6\,{\rm Gyr}$, whereas the youngest stars in the thick discs of Milky-Way-like galaxies are $\sim10\,{\rm Gyr}$ old. I tested this prediction against the scarcely available thick disc age estimates, all of them are from galaxies with $\mathcal{M}_\star(z=0)\gtrsim10^{10}\,\mathcal{M}_\odot$, and I find that field spiral galaxies seem to follow the expectation. On the other hand, my derivation predicts ages that are too low for the thick discs in lenticular galaxies, indicating a fast early evolution for S0 galaxies. I propose the idea of conclusively testing whether thick discs formed quickly and in situ by obtaining the ages of thick discs in field galaxies with masses of $\mathcal{M}_\star(z=0)\sim10^{9.5}\,\mathcal{M}_\odot$ and by checking whether they contain $\sim5\,{\rm Gyr}$-old stars., Comment: Letter accepted for publication in A&A

The dominant mechanism leading to the formation of brown dwarfs (BDs) remains uncertain. The most direct keys to formation, which are obtained from younger objects (pre-BD cores and proto-BDs), are limited by the very low number statistics available. We aim to identify and characterize a set of pre- and proto-BDs as well as Class II BDs in the Lupus 1 and 3 molecular clouds to test their formation mechanism. We performed ALMA band 6 (1.3 mm) continuum observations of a selection of 64 cores previously identified from AzTEC/ASTE data (1.1 mm), along with previously known Class II BDs in the Lupus 1 and 3 molecular clouds. Surveyed archival data in the optical were used to complement these observations. We expect these ALMA observations prove efficient in detecting the youngest sources in these regions, since they probe the frequency domain at which these sources emit most of their radiation. We detected 19 sources from 15 ALMA fields. Considering all the pointings in our observing setup, the ALMA detection rate was $\sim$23% and the derived masses of the detected sources were between $\sim$0.18 and 124 $\mathrm{M_{Jup}}$. We classified these sources according to their spectral energy distribution as 5 Class II sources, 2 new Class I/0 candidats, and 12 new possible pre-BD or deeply embedded protostellar candidates. We detected a promising candidate for a Class 0/I proto-BD source and inferred the disk dust mass of a bona fide Class II BD. The pre-BD cores might be the byproduct of an ongoing process of large-scale collapse. The Class II BD disks follow the correlation between disk mass and the mass of the central object that is observed at the low-mass stellar regime. We conclude that it is highly probable that the sources in the sample are formed as a scaled-down version of low-mass star formation, although disk fragmentation may be responsible for a considerable fraction of BDs., Comment: 21 pages, 10 figures, 13 tables, 3 Appendixes. Accepted in A&A

AGN are a key ingredient for understanding galactic evolution. AGN-driven outflows are one of the manifestations of feedback. The AO mode for MUSE at the VLT permits to study the innermost tens of parsecs of nearby AGN in the optical. We present a detailed analysis of the ionised gas in the central regions of NGC 7130, an archetypical composite Seyfert and nuclear starburst galaxy. We achieve an angular resolution of 0.17$^{\prime\prime}$ (50 pc). We performed a multi-component analysis of the main ISM lines and identified nine kinematic components, six of which correspond to the outflow. The outflow is biconic and has velocities of a few $100\,{\rm km\,s^{-1}}$ with respect to the disc. We decompose the approaching side of the outflow into a broad and a narrow component with typical velocity dispersions below and above $\sim200\,{\rm km\,s^{-1}}$, respectively. The blueshifted narrow component has substructure, in particular a collimated plume aligned with the radio jet, indicating that it may be jet-powered. The redshifted lobe is composed of two Narrow Components and a Broad Component. An additional redshifted component is seen outside the main outflow axis. Line ratio diagnostics indicate that the outflow gas in the main axis is AGN-powered whereas the off-axis component has LINER properties. The ionised gas mass outflow rate is $\dot{M}=1.2\pm0.7\,M_{\odot}\,{\rm yr^{-1}}$ and the kinetic power is $\dot{E}_{\rm kin}=(2.7\pm2.0)\times10^{41}\,{\rm erg\,s^{-1}}$, which corresponds to $F_{\rm kin}=0.12\pm0.09\%$ of the bolometric AGN power. The combination of high angular resolution integral field spectroscopy and a careful multi-component decomposition allows a uniquely detailed view of the outflow in NGC 7130, illustrating that AGN kinematics are more complex than traditionally derived from less sophisticated data and analyses. (abridged), Comment: Published in A&A. This version of the preprint has been corrected with the modifications appearing in the A&A corrigendum

The Cygnus region, which dominates the local spiral arm of the Galaxy, is one of the nearest complexes of massive star formation. Its massive stellar content, regions of ongoing star formation, and molecular gas have been studied in detail. However, little is known of the history of the region beyond the past 10 Myr. The brightness and spectroscopic characteristics of red supergiants make it easy to identify them and build up a virtually complete sample of such stars at the distance of the Cygnus region, thus providing a record of massive star formation extending several tens of Myr into the past, a period inaccessible through the O and early B stars observable at present. We have made a selection of a sample of bright, red stars in an area of 84 square degrees covering the whole present extension of the Cygnus region. We have obtained spectroscopy in the red visible range allowing an accurate, homogeneous spectral classification as well as a reliable separation between supergiants and other cool stars. Our data are complemented with Gaia Data Release 2 astrometric data. We have identified 29 red supergiants in the area, 17 of which had not been previously classified as supergiants. Twenty-four of the 29 most likely belong to the Cygnus region and four of the remaining to the Perseus arm. We have used their derived luminosities and masses to infer the star formation history of the region. Intense massive star formation activity is found to have started approximately 15 Myr ago, and we find evidence for two other episodes, one taking place between 20 and 30 Myr ago and another one having ended approximately 40 Myr ago. There are small but significant differences between the kinematic properties of red supergiants younger or older then 20 Myr, hinting that stars of the older group were formed outside the precursor of the present Cygnus complex, possibly in the Sagittarius-Carina arm., Comment: Accepted for publication in Astronomy and Astrophysics

The Vanishing & Appearing Sources during a Century of Observations (VASCO) project investigates astronomical surveys spanning a time interval of 70 years, searching for unusual and exotic transients. We present herein the VASCO Citizen Science Project, which can identify unusual candidates driven by three different approaches: hypothesis, exploratory, and machine learning, which is particularly useful for SETI searches. To address the big data challenge, VASCO combines three methods: the Virtual Observatory, user-aided machine learning, and visual inspection through citizen science. Here we demonstrate the citizen science project and its improved candidate selection process, and we give a progress report. We also present the VASCO citizen science network led by amateur astronomy associations mainly located in Algeria, Cameroon, and Nigeria. At the moment of writing, the citizen science project has carefully examined 15,593 candidate image pairs in the data (ca. 10% of the candidates), and has so far identified 798 objects classified as "vanished". The most interesting candidates will be followed up with optical and infrared imaging, together with the observations by the most potent radio telescopes., Comment: 16 pages, 7 figures, contribution to IAC 2020 conference

We investigate the spatial distribution of star formation (SF) within bars of nearby disk galaxies (inclination $< 65^{\circ}$) from the S$^4$G survey. We use archival GALEX far- and near-UV imaging for 772 barred galaxies. We also assemble a compilation of continuum-subtracted H$\alpha$ images for 433 barred galaxies, of which 70 are produced by ourselves from ancillary photometry and MUSE/CALIFA IFU data cubes. We employ two complementary approaches: i) the analysis of bar/disk stacks built from co-added UV images of hundreds of galaxies; and ii) the classification of the morphology of ionised regions in galaxies into three main SF classes: A) only circumnuclear SF, B) SF at the bar ends, but not along the bar, and C) SF along the bar. Lenticular galaxies typically belong to SF class A: this is probably related to bar-induced SF quenching. The distribution of SF class B peaks for early- and intermediate-type spirals: this most likely results from the interplay of gas flow, shocks, and enhanced shear in centrally concentrated galaxies with large bar amplitudes. Late-type galaxies are mainly assigned to SF class C: we argue that this is a consequence of low shear. In bar stacks of spirals, the UV emission traces the stellar bars and dominates on their leading side, as witnessed in simulations. For early-types, the central UV emission is $\sim$0.5 mag brighter in strongly barred galaxies, relative to their weakly barred counterparts: this is related to the efficiency of strong bars sweeping the disk gas and triggering central starbursts. We also show that the distributions of SF in inner ringed galaxies are broadly the same in barred and non-barred galaxies, including a UV/H$\alpha$ deficit in the middle part of the bar: this hints at the effect of resonance rings trapping gas. Distinct distributions of SF within bars are reported in galaxies of different morphological types (Abridged)., Comment: Accepted for publication in A&A (September 2, 2020). 18 pages, 14 figures, 2 tables, plus appendices (16 pages, 8 figures, 1 table)

Cygnus OB2 provides a unique insight into the high-mass stellar content in one of the largest groups of young massive stars in our Galaxy. Although several studies of its massive population have been carried out over the last decades, an extensive spectroscopic study of the whole known O-star population in the association is still lacking. In this work, we created the most complete spectroscopic census of O stars carried out so far in Cygnus OB2 using already existing and new spectroscopy. We present the spectra for 78 O-type stars, from which we identify new binary systems, obtain the distribution of rotational velocities, and determine the main stellar parameters for all the stars in the region that have not been detected as double-line spectroscopic binaries. We also derive radii, luminosities, and masses for those stars with reliable Gaia astrometry, in addition to creating the Hertzsprung-Russell Diagram to interpret the evolutionary status of the association. This work has shown the improvement reached when using accurate spectroscopic parameters and astrometry for the interpretation of the evolutionary status of a population, revealing, in the case of Cygnus OB2, at least two star-forming bursts at $\sim$3 and $\sim$5 Myr. We find an apparent deficit of very fast rotators in the distribution of rotational velocities. The inspection of the dynamical distribution of the sample has allowed us to identify nine O stars with peculiar proper motions and discuss a possible dynamical ejection scenario or past supernova explosions in the region., Comment: Accepted for publication in A&A, 25 pages, 11 figures

Very low-mass stars are known to have jets and outflows, which is indicative of a scaled-down version of low-mass star formation. However, only very few outflows in very low-mass sources are well characterized. We characterize the bipolar molecular outflow of the very low-mass star Par-Lup3-4, a 0.12 M$_{\odot}$ object known to power an optical jet. We observed Par-Lup3-4 with ALMA in Bands 6 and 7, detecting both the continuum and CO molecular gas. In particular, we studied three main emission lines: CO(2-1), CO(3-2), and $^{13}$CO(3-2). Our observations reveal for the first time the base of a bipolar molecular outflow in a very low-mass star, as well as a stream of material moving perpendicular to the primary outflow of this source. The primary outflow morphology is consistent with the previously determined jet orientation and disk inclination. The outflow mass is $9.5\times10^{-7}\mathrm{M}_{\odot}$ , with an outflow rate of $4.3\times10^{-9}\mathrm{M}_{\odot}\mathrm{yr}^{-1}$ A new fitting to the spectral energy distribution suggests that Par-Lup3-4 may be a binary system. We have characterized Par-Lup3-4 in detail, and its properties are consistent with those reported in other very low-mass sources. This source provides further evidence that very low-mass sources form as a scaled-down version of low-mass stars., Comment: 20 pages, 11 figures, 5 tables. Accepted in A&A

A statistical study of intermediate Palomar Transient Factory supernovae (SNe) in Type 1 AGN has shown a major deficit of supernovae around Type 1 AGN host galaxies, with respect to Type 2 AGN hosts. The aim of this work is to test whether there is any preference for Type 1 AGN to host SN of a specific kind. Through the analysis of SN occurrence and their type (thermonuclear vs core-collapse), we can directly link the type of stars producing the SN events, thus this is an indirect way to study host galaxies in Type 1 AGN. We examine the detection fractions of SNe, the host galaxies and compare the sample properties to typical host galaxies in the Open Supernova Catalog (OSC; Guillochon et al. 2017). The majority of the host galaxies in the AGN sample are late-type, similar to typical galaxies hosting SN within the OSC. The findings are supportive of a deficiency of SNe near Type 1 AGN, although we cannot with certainty assess the overall detection fractions of SNe in Type 1 AGN relative to other SN host galaxies. We can state that Type 1 AGN has equal detection fractions of thermonuclear vs core-collapse SNe. However, we note the possibility of a higher detection rate of core-collapse supernovae in Type-1 AGN with insecure AGN classifications., Comment: 12 pages, 5 figures, 3 tables. Accepted for publication in MNRAS

Very few examples are known of red supergiant runaways, all of them descending from the more massive O-type precursors, but none from the lower mass B-type precursors, although runaway statistics among B-type stars suggest that K-type runaways must be relatively numerous. We study HD 137071, a star that has been considered so far as a normal K-type red giant. Its parallax measured by Gaia and the derived luminosity suggest that it is actually a supergiant, whereas its derived distance to the galactic plane and its spatial velocity of 54.1 km s$^{-1}$ with respect to the local standard of rest suggest that it is also a runaway star. However, intrinsic limitations in determining the trigonometric parallaxes of cool supergiants, even in the Gaia era, require accurate spectral classifications for confirmation. We reliably classify HD 137071 as a K4II star establishing its membership to the extreme Population I, which is in agreement with the luminosity derived using the Gaia DR2 parallax measurement. Kinematical data from the Gaia DR2 catalog confirm its high spatial velocity and its runaway nature. Combining the spectral classification with astrometric information, a state-of-the-art galactic potential model, and evolutionary models for high-mass stars we trace the motion of HD 137071 back to the proximities of the galactic plane and speculate on which of the two proposed mechanisms for the production of runaway stars may be responsible for the high velocity of HD 137071. The available data favor the formation of HD 137071 in a massive binary system where the more massive companion underwent a supernova explosion about 32 Myr ago., Comment: 6 pages, 4 figures. To be published in Astronomy and Astrophysics

Using SPITZER 3.6$\mu$m imaging, we investigate the physical and data-driven origins of up-bending (Type III) disk breaks. We apply a robust new break-finding algorithm to 175 low-inclination disk galaxies previously identified as containing Type III breaks, classify each galaxy by its outermost re-classified (via our new algorithm) break type, and compare the local environments of each resulting subgroup. Using three different measures of the local density of galaxies, we find that galaxies with extended outer spheroids (Type IIIs) occupy the highest density environments in our sample, while those with extended down-bending (Type II) disks and symmetric outskirts occupy the lowest density environments. Among outermost breaks, the most common origin of Type III breaks in our sample is methodological; the use of elliptical apertures to measure the radial profiles of asymmetric galaxies usually results in features akin to Type III breaks., Comment: 4 pages, 1 figure, for IAU Symposium 355

In this paper we report the current status of a new research program. The primary goal of the "Vanishing & Appearing Sources during a Century of Observations" (VASCO) project is to search for vanishing and appearing sources using existing survey data to find examples of exceptional astrophysical transients. The implications of finding such objects extend from traditional astrophysics fields to the more exotic searches for evidence of technologically advanced civilizations. In this first paper we present new, deeper observations of the tentative candidate discovered by Villarroel et al. (2016). We then perform the first searches for vanishing objects throughout the sky by comparing 600 million objects from the US Naval Observatory Catalogue (USNO) B1.0 down to a limiting magnitude of $\sim 20 - 21$ with the recent Pan-STARRS Data Release-1 (DR1) with a limiting magnitude of $\sim$ 23.4. We find about 150,000 preliminary candidates that do not have any Pan-STARRS counterpart within a 30 arcsec radius. We show that these objects are redder and have larger proper motions than typical USNO objects. We visually examine the images for a subset of about 24,000 candidates, superseding the 2016 study with a sample ten times larger. We find about $\sim$ 100 point sources visible in only one epoch in the red band of the USNO which may be of interest in searches for strong M dwarf flares, high-redshift supernovae or other catagories of unidentified red transients., Comment: SETI meets time domain astronomy. Accepted into the Astronomical Journal

During the 2017 record-breaking burning season in Canada / United States, intense wild fires raged during the first week of September in the Pacific northwestern region (British Columbia, Alberta, Washington, Oregon, Idaho, Montana and northern California) burning mostly temperate coniferous forests. The heavy loads of smoke particles emitted in the atmosphere reached the Iberian Peninsula (IP) a few days later on 7 and 8 September. Satellite imagery allows to identify two main smoke clouds emitted during two different periods that were injected and transported in the atmosphere at several altitude levels. Columnar properties on 7 and 8 September at two Aerosol Robotic Network (AERONET) mid-altitude, background sites in northern and southern Spain are: aerosol optical depth (AOD) at 440 nm up to 0.62, Angstrom exponent of 1.6-1.7, large dominance of small particles (fine mode fraction > 0.88), low absorption AOD at 440 nm (<0.008) and large single scattering albedo at 440 nm (>0.98). Profiles from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) show the presence of smoke particles in the stratosphere during the transport, whereas the smoke is only observed in the troposphere at its arrival over the IP. Portuguese and Spanish ground lidar stations from the European Aerosol Research Lidar Network / Aerosols, Clouds, and Trace gases Research InfraStructure Network (EARLINET/ACTRIS) and the Micro-Pulse Lidar NETwork (MPLNET) reveal smoke plumes with different properties: particle depolarization ratio and color ratio, respectively, of 0.05 and 2.5 in the mid troposphere (5-9 km) and of 0.10 and 3.0 in the upper troposphere (10-13 km). In the mid troposphere the particle depolarization ratio does not seem time-dependent during the transport whereas the color ratio seems to increase (larger particles sediment first)., Comment: 19 pages, 10 figures

We investigate the silicate feature of the two Class I components of HH250-IRS, a resolved binary system with a separation of $0''53$ driving a Herbig-Haro flow. Each component has its own circumstellar envelope, and the system is surrounded by a circumbinary disk. We have carried out low resolution spectroscopy in the 8-13$\mu$m range using VISIR at ESO's Very Large Telescope. The silicate features of both sources are clearly different. The NW component has a broad, smooth absorption profile lacking structure. The SE component shows the silicate feature in emission, with structure longwards of 9.5$\mu$m indicating the presence of crystalline dust in the dominant form of forsterite. The apparent lack of an absorption feature caused by foreground dust is probably due to the filling of the band with emission by amorphous silicates in the envelope of the object. Despite their virtually certain coevality, the differences in the components of the HH250-IRS binary are most likely due to markedly different circumstellar environments. The NW component displays an unevolved envelope, whereas dust growth and crystallization has taken place in the SE component. The weak or absent signatures of enstatite in the latter are fairly unusual among envelopes with crystalline dust, and we tentatively relate it to a possible wide gap or an inner truncation of the disk already hinted in previous observations by a drop in the $L'$-band flux, which might indicate that the SE component could actually be a very close binary. We speculate that the clear differences between the silicate feature spectra of both components of HH250-IRS may be due either to disk evolution sped up by multiplicity, or by accretion variability leading to episodes of crystal formation., Comment: Accepted by Astronomy and Astrophysics

Aims: Using a sample of 175 low-inclination galaxies from the S$^{4}$G, we investigate the origins of up-bending (Type III) breaks in the 3.6 $\mu$m surface brightness profiles of disk galaxies. Methods: We re-analyze a sample of previously identified Type III disk break-hosting galaxies using a new, unbiased break-finding algorithm, which uncovered many new, sometimes subtle disk breaks across the whole sample. We classify each break by its likely origin through close examination of the galaxy images across wavelengths, and compare samples of galaxies separated by their outermost identified break types in terms of their stellar populations and local environments. Results: We find that more than half of the confirmed Type III breaks in our sample can be attributed to morphological asymmetry in the host galaxies. As these breaks are mostly an artifact of the azimuthal averaging process, their status as physical "breaks" is questionable. Such galaxies occupy some of the highest density environments in our sample, implying that much of this asymmetry is the result of tidal disturbance. Additionally, we find that Type III breaks related to extended spiral arms or star formation often host down-bending (Type II) breaks at larger radius which were previously unidentified. Such galaxies reside in the lowest density environments in our sample, in line with previous studies that found a lack of Type II breaks in clusters. Galaxies occupying the highest density environments most often show Type III breaks associated with outer spheroidal components. Conclusions: We find that Type III breaks in the outer disks of galaxies arise most often through environmental influence: either tidal disturbance (resulting in disk asymmetry) or heating through, e.g., galaxy harrassment (leading to spheroidal components). Galaxies hosting the latter break types also show... (abstract continues), Comment: 21 pages, 10 figures, accepted for publication in A&A

Thick discs are nearly ubiquitous components of the discs of present-day galaxies. It has been proposed that a fraction of their stars has been accreted. Here, we aim to find whether accretion of satellites is the main thick disc formation mechanism. To do so, we observed a sample of eight nearby edge-on galaxies with the MUSE integral field unit at the VLT. Six of the galaxies have a distinct thick disc. We derived thick disc velocities and velocity dispersions for the galaxies in our sample. We devise a formalism to estimate the fractions of retrograde material in the thick discs by using kinematical maps and thin/thick dis decompositions. None of the galaxies in our sample shows strong evidence for retrograde material at large distances from the centre. Including those found in the literature, there are seventeen thick discs with studied kinematics, with only one showing unambiguous signatures of retrograde material. Literature numerical studies of dynamical friction allow us to estimate that at the current cosmic time about one in six mergers for which the stars of the accreted galaxy ended in a thick disc were retrograde. This is in tension with the observed fraction of 1/17 of galaxies with a partly retrograde thick disc. We conclude that satellite accretion is not favoured by observations to be the main thick disk formation mechanism., Comment: Accepted for publication in A&A

Observations of the galactic disk at mid-infrared and longer wavelengths reveal a wealth of structures indicating the existence of complexes of recent massive star formation. However, little or nothing is known about the stellar component of those complexes. We have carried out observations aiming at the identification of early-type stars in the direction of the bright infrared source RAFGL~5475, around which several interstellar medium structures usually associated with the presence of massive stars have been identified. Our observations have the potential of revealing the suspected but thus far unknown stellar component of the region around RAFGL~5475. We have carried out near-infrared imaging observations ($JHK_S$ bands) designed to reveal the presence of early-type stars based on their positions in color-color and color-magnitude diagrams centered on the location of RAFGL~5475. We took into account the possibility that candidates found might belong to a foreground population physically related either to M16 or M17, two giant HII regions lying midway between the Sun and RAFGL~5475. The near-infrared color-color diagram shows clear evidence for the presence of a moderately obscured population of early-type stars in the region imaged. By studying the distribution of extinction in their direction and basic characteristics of the interstellar medium we show that these new early-type stars are most likely associated with RAFGL~5475. By investigating the possible existence of massive early-type stars in the direction of RAFGL~5475 we have discovered the existence of a new OB association. A very preliminary assessment of its contents suggests the presence of several O-type stars, some of them likely to be associated with structures in the interstellar medium. The new association is located at 4 kpc from the Sun in the Scutum-Centaurus arm., Comment: 8 pages, 7 figures; accepted by Astronomy and Astrophysics

Galaxy velocity maps often show the typical pattern of a rotating disk, consistent with the dynamical model where emitters rotate in circular orbits around the galactic center. The simplest template used to fit these maps consists in the rotating disk model (RDM) where the amplitude of circular velocities is fixed by the observed velocity profile along the kinematic axis. A more sophisticated template is the rotating tilted-ring model (RTRM) that takes into account the presence of warps and allows a radius-dependent orientation of the kinematic axis. In both cases, axisymmetry is assumed and residuals between the observed and the model velocity fields are interpreted as noncircular motions. We show that if a galaxy is not axisymmetric, there is an intrinsic degeneracy between a rotational and a radial velocity field. We then introduce a new galaxy template, the radial ellipse model (REM), that is not axisymmetric and has a purely radial velocity field with an amplitude that is correlated with the major axis of the ellipse. We show that best fits to the observed two-dimensional velocity fields of 28 galaxies extracted from the THINGS sample with both the REM and the RDM give residuals with similar amplitudes, where the REM residuals trace nonradial motions. Best fits obtained with the RTRM, because of its larger number of free parameters, give the smallest residuals: however, we argue that this does not necessarily imply that the RTRM gives the most accurate representation of a galaxy velocity field. Instead, we show that this method is not able to disentangle between circular and radial motions for the case of nonaxisymmetric systems. Finally, we consider the physical motivation of the REM, and discuss how the interpretation of galactic dynamics changes if one assumes that the main component of a galaxy velocity field is modeled as a RDM/RTRM or as a REM., Comment: 26 pages, 32 figures. Astronomy and Astrophysics in the press

We present the discovery of a small kinematically decoupled core of 0.2$^{\prime\prime}$ (60 pc) in radius as well as an outflow jet in the archetypical AGN-starburst "composite" galaxy NGC 7130 from integral field data obtained with the adaptive optics-assisted MUSE-NFM instrument on the VLT. Correcting the already good natural seeing at the time of our science verification observations with the four-laser GALACSI AO system, we reach an unprecedented spatial resolution at optical wavelengths of around 0.15$^{\prime\prime}$. We confirm the existence of star-forming knots arranged in a ring of 0.58$^{\prime\prime}$ (185 pc) in radius around the nucleus, previously observed from UV and optical Hubble Space Telescope and CO(6-5) ALMA imaging. We determine the position of the nucleus as the location of a peak in gas velocity dispersion. A plume of material extends towards the NE from the nucleus until at least the edge of our field of view at 2$^{\prime\prime}$ (640 pc) radius which we interpret as an outflow jet originating in the AGN. The plume is not visible morphologically, but is clearly characterised in our data by emission-line ratios characteristic of AGN emission, enhanced gas velocity dispersion, and distinct non-circular gas velocities. Its orientation is roughly perpendicular to the line of nodes of the rotating host galaxy disc. A circumnuclear area of positive and negative velocities of 0.2$^{\prime\prime}$ in radius indicates a tiny inner disc, which can only be seen after combining the integral field spectroscopic capabilities of MUSE with adaptive optics., Comment: Accepted for publication in A&A letters

Cygnus OB2 is a rich OB association in the Galaxy which has experienced intense star formation in the last 20-25 Myr. Its stellar population shows a correlation between age and Galactic longitude. Exploring the chemical composition of its stellar content we will be able to check the degree of homogeneity of the natal molecular cloud and possible effects of self-enrichment processes. Our aim is to determine silicon and oxygen abundances for a sample of eight early-type slow rotators in Cygnus OB2 in order to check possible inhomogeneities across the whole association and whether there exists a correlation of chemical composition with Galactic longitude. We have performed a spectroscopic analysis of a sample of late O and early B stars with low rotational velocity, which have been chosen so as to cover the whole association area. We have carried out an analysis based on equivalent widths of metal lines, the wings of the H Balmer lines and FASTWIND stellar atmosphere models to determine their stellar fundamental parameters as well as the silicon and oxygen surface abundances. We derive a rather homogeneous distribution of silicon and oxygen abundances across the region, with average values of 12+log(Si/H)=7.53$\pm$0.08 dex and 12+log(O/H)=8.65$\pm$0.12 dex. We find a homogeneous chemical composition in Cygnus OB2 with no clear evidence for significant chemical self-enrichment, despite indications of strong stellar winds and possible supernovae during the history of the region. Comparison with different scenarios of chemical enrichment by stellar winds and supernovae point to star forming efficiencies not significantly above 10%. The degree of homogeneity that we find is consistent with the observed Milky Way oxygen gradient based on HII regions. We also find that the oxygen scatter within Cygnus OB2 is at least of the same order than among HII regions at similar Galactocentric distance., Comment: 15 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

It has been suggested that the well-studied giant HII regions M16 and M17 may have had a common origin, being an example of large-scale triggered star formation. While some features of the distribution of the interstellar medium in the region support this interpretation, no definitive detection of an earlier population of massive stars responsible for the triggering has been made thus far. We have carried out observations looking for red supergiants in the area covered by a giant shell seen in HI and CO centered on galactic coordinates $l \sim 14^\circ 5$, $b\sim +1^\circ$ that peaks near the same radial velocity as the bulk of the emission from both giant HII regions, which are located along the shell. Red supergiants have ages in the range expected for the parent association whose most massive members could have triggered the formation of the shell and of the giant HII regions along its rim. Out of a sample of 37 bright red stars, we identify four red supergiants that confirm the existence of massive stars in the age range between $\sim 10$ and $\sim 30$~Myr in the area. At least three of them have Gaia DR2 parallaxes consistent with them being at the same distance as M16 and M17. The evidence of past massive star formation within the area of the gaseous shell lends support to the idea that it was formed by the combined action of stellar winds and ionizing radiation of the precursors of the current red supergiants. These could be the remnants of a richer population, whose most massive members have exploded already as core-collapse supernovae. The expansion of the shell against the surrounding medium, perhaps combined with the overrun of preexisting clouds, is thus a plausible trigger of the formation of a second generation of stars currently responsible for the ionization of M16 and M17., Comment: 7 pages, 3 figures; accepted for publication in Astronomy and Astrophysics

Extensive surveys of star-forming regions with Spitzer have revealed populations of disk-bearing young stellar objects. These have provided crucial constraints, such as the timescale of dispersal of protoplanetary disks, obtained by carefully combining infrared data with spectroscopic or X-ray data. While observations in various regions agree with the general trend of decreasing disk fraction with age, the Lupus V and VI regions appeared to have been at odds, having an extremely low disk fraction. Here we show, using the recent Gaia data release 2 (DR2), that these extremely low disk fractions are actually due to a very high contamination by background giants. Out of the 83 candidate young stellar objects (YSOs) in these clouds observed by Gaia, only five have distances of 150 pc, similar to YSOs in the other Lupus clouds, and have similar proper motions to other members in this star-forming complex. Of these five targets, four have optically thick (Class II) disks. On the one hand, this result resolves the conundrum of the puzzling low disk fraction in these clouds, while, on the other hand, it further clarifies the need to confirm the Spitzer selected diskless population with other tracers, especially in regions at low galactic latitude like Lupus V and VI. The use of Gaia astrometry is now an independent and reliable way to further assess the membership of candidate YSOs in these, and potentially other, star-forming regions., Comment: Accepted for publication on Astronomy&Astrophysics Letters

The central area (40"x40") of the bipolar nebula S106 was mapped in the OI line at 63.2 micron with high angular (6") and spectral resolution, using GREAT on board SOFIA. The OI emission distribution is compared to the CO 16-15, CII 158 micron, and CO 11-10 lines, mm-molecular lines, and continuum. It is composed of several velocity components in the range from -30 km/s to 25 km/s. The high-velocity blue- and redshifted emission can be explained as arising from accelerated photodissociated (PDR) gas associated with a dark lane close to the massive binary system S106 IR, and from shocks caused by the stellar wind and/or a disk--envelope interaction. At velocities from -9 to -4 km/s and 0.5 to 8 km/s line wings are observed that we attribute to cooling in PDRs created by the ionizing radiation impinging on the cavity walls. The bulk velocity range is dominated by PDR emission from the clumpy molecular cloud. Modelling the emission in the different velocity ranges with the KOSMA-tau code constrains a radiation field chi of a few times 10^4 and densities n of a few times 10^4 cm^-3. Considering self-absorption of the OI line results in higher densities (up to 10^6 cm^-3) only for the gas component seen at high blue- and red velocities. The dark lane has a mass of 275 Msun and shows a velocity difference of 1.4 km/s along its projected length of 1 pc, determined from H13CO+ 1-0 mapping. It can be interpreted as a massive accretion flow, or the remains of it, linked to S106 IR/FIR. The most likely explanation is that the binary system is at a stage of its evolution where gas accretion is counteracted by the stellar winds and radiation, leading to the very complex observed spatial and kinematic emission distribution of the various tracers., Comment: Final version as it will appear in Asronomy&Astrophysics, accepted 22.5.2018

S106 is one of the best known bipolar HII regions, thoroughly studied and modelled at infrared, submillimeter and millimeter wavelengths, and it is one of the nearest examples of the late stages of massive star formation in which the newly formed star that ionizes it is still surrounded by vast amounts of gas and dust. However, little is known about its heavily obscured central source, S106IR. The possible binarity of the central source is investigated, which is considered to be likely given the high binarity fraction among massive stars. We have carried out visible and near-infrared photometric monitoring looking for short-term variability, with special interest in that related to the presence of a close binary companion to S106IR that may produce periodic eclipses or tidal distortion of the shape of the members of the system. A periodic variability of S106IR in the J band is found with a period of 5.0 days and an amplitude of about 0.1 mag. The light curve displays a slow rise from minimum to maximum followed by a steep decrease, and can be well reproduced by a close binary system composed of two stars with different luminosity orbiting each other in an elliptical orbit of moderate eccentricity. S106IR also shows hints of short-term variability possibly related to accretion. We also report variability of four other stars previously classified as members of the S106 cluster, all of which are strong X-ray emitters. The newly discovered close binarity of S106IR adds a new element to the modeling of the nebula and to the understanding of the dynamics of the gas around the ionizing source, which suggests that the components of the binary are accreting via a circumbinary disk. Binarity also helps to explain the apparent mismatch between the spectral type of the ionizing source inferred from the nebular spectrum and its high brightness at near-infrared wavelengths., Comment: Accepted for publication in Astronomy and Astrophysics

Herbig-Haro flows are signposts of recent major accretion and outflow episodes. We aim to determine the nature and properties of the little-known outflow source HH 250-IRS, which is embedded in the Aquila clouds. We have obtained adaptive optics-assisted L-band images with the NACO instrument on the Very Large Telescope (VLT), together with N- and Q-band imaging with VISIR also on the VLT. Using the SINFONI instrument on the VLT we carried out H- and K-band integral field spectroscopy of HH 250-IRS, complemented with spectra obtained with the SpeX instrument at the InfraRed Telescope Facility (IRTF) in the JHKL bands. Finally, the SubMillimeter Array (SMA) interferometer was used to study the circumstellar environment of HH 250-IRS at 225 and 351 GHz with CO (2-1) and CO (3-2) maps and 0.9 mm and 1.3 mm continuum images. The HH 250-IRS source is resolved into a binary with 0''53 separation, corresponding to 120 AU at the adopted distance of 225 pc. The individual components show heavily veiled spectra with weak CO absorption indicative of late-type stars. Both are Class I sources, but their spectral energy distributions between 1.5 $\mu$m and 19 $\mu$m differ markedly and suggest the existence of a large cavity around one of the components. The millimeter interferometric observations indicate that the gas mainly traces a circumbinary envelope or disk, while the dust emission is dominated by one of the circumstellar envelopes. HH 250-IRS is a new addition to the handful of multiple systems where the individual stellar components, the circumstellar disks and a circumbinary disk can be studied in detail, and a rare case among those systems in which a Herbig-Haro flow is present., Comment: Accepted for publication by Astronomy and Astrophysics