Search

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

The Cygnus complex is one of the most powerful star forming regions at a close distance from the Sun (~1.4 kpc). Its richest OB association Cygnus OB2 is known to harbor many tens of O-type stars and hundreds of B-type stars, providing a large homogeneous population of OB stars that can be analyzed. Many studies of its massive population have been developed in the last decades, although the total number of OB stars is still incomplete. Our aim is to increase the sample of O and B members of Cygnus OB2 and its surroundings by spectroscopically classifying 61 candidates as possible OB-type members of Cygnus OB2. We have obtained new blue intermediate-resolution spectra suitable for spectral classification of the 61 candidates in Cygnus OB2 and surroundings. We thus performed a spectral classification of the sample using He I-II and metal lines rates, as well as the Marxist Ghost Buster (MGB) software for O-type stars and the IACOB standards catalog for B-type stars. Out of the 61 candidates, we have classified 42 stars as new massive OB-type stars, earlier than B3, in Cygnus OB2 and surroundings, including 11 O-type stars. The other candidates are discarded as they display later spectral types inconsistent with membership in the association. However, the magnitude cutoff and dust extinction introduce an incompleteness. Many O and early B stars at B > 16 mag are still undiscovered in the region. Finally, we have studied the age and extinction distribution of our sample within the region, placing them in the Hertzsprung-Russell Diagram using different stellar models in order to assess age uncertainties. Massive star formation in Cygnus OB2 seems to have proceeded from lower to higher Galactic longitudes, regardless of the details of the models used. The correlation between age and Galactic longitude previously found in the region is now confirmed., Comment: 20 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

Recent studies have made the community aware of scattered light when examining low-surface-brightness galaxy features such as thick discs. In our past studies of the thick discs of edge-on galaxies in the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G) we modelled the point spread function as a Gaussian. We re-examine our results using a revised point spread function model that accounts for extended wings out to more than 2.5arcmin. We study the $3.6\mu{\rm m}$ images of 141 edge-on galaxies from the S$^4$G. We decompose the surface brightness profiles of the galaxies perpendicular to their mid-planes assuming that discs are made of two stellar discs in hydrostatic equilibrium. We decompose the axial surface brightness profiles of galaxies to model the central mass concentration - described by a S\'ersic function - and the disc - described by a broken exponential disc. Our improved treatment confirms the ubiquity of thick discs. The main difference between our current fits and those presented before is that now the scattered light from the thin disc dominates the surface brightness at levels below $\mu\sim26\,{\rm mag\,arcsec^{-2}}$. This does not affect drastically any of our previously presented results: 1) Thick discs are nearly ubiquitous. They are not an artefact caused by scattered light as has been suggested elsewhere. 2) Thick discs have masses comparable to those of thin discs in low-mass galaxies - circular velocities $v_{\rm c}<120\,{\rm km\,s^{-1}}$ - whereas they are typically less massive than the thin discs in high-mass galaxies. 3) Thick discs and central mass concentrations seem to have formed at the same epoch from a common material reservoir. 4) Roughly 60% of the up-bending breaks in face-on galaxies are caused by the superposition of a thin and a thick disc where the scale-length of the latter is the largest. (Abridged), Comment: Accepted for publication in A&A. 26 pages and 22 figures. The appendices are truncated due to size limitations. A full version of the paper can be accessed at http://research.iac.es/galeria/sebastiencomeron/psf_v11.pdf

Active Galactic Nuclei (AGN) are extremely powerful cosmic objects, driven by accretion of hot gas upon super-massive black holes. The zoo of AGN classes are divided into two major groups, with Type-1 AGN displaying broad Balmer emission lines and Type-2 narrow ones. For a long time it was believed that a Type-2 AGN is a Type-1 AGN viewed through a dusty kiloparsec-size torus, but an emerging body of observations suggests more than just the viewing angle matters. Here we report significant differences in supernova counts and classes in the first study to date of supernovae near Type-1 and Type-2 AGN host galaxies, using data from the intermediate Palomar Transient Factory, the Sloan Digital Sky Survey Data Release 7 and Galaxy Zoo. We detect many more supernovae in Type-2 AGN hosts (size of effect $\sim$ 5.1$\sigma$) compared to Type-1 hosts, which shows that the two classes of AGN are located inside host galaxies with different properties. In addition, Type-1 and Type-2 AGN that are dominated by star formation according to WISE colours $m_{W1} - m_{W2} < 0.5$ and are matched in 22 $\mu$m absolute magnitude differ by a factor of ten in $L$[OIII]$\lambda$5007 luminosity, suggesting that when residing in similar type of host galaxies Type-1 AGN are much more luminous. Our results demonstrate two more factors that play an important role in completing the current picture: the age of stellar populations and the AGN luminosity. This has immediate consequences for understanding the many AGN classes and galaxy evolution., Comment: Accepted to the Astrophysical Journal (ApJ)

Globules and pillars, impressively revealed by the Spitzer and Herschel satellites, for example, are pervasive features found in regions of massive star formation. We studied the globule IRAS 20319+3958 in Cygnus X by means of visible and near-infrared imaging and spectroscopy, complemented with mid-infrared Spitzer/IRAC imaging, in order to obtain a census of its stellar content and the nature of its embedded sources. Our observations show that the globule contains an embedded aggregate of about 30 very young ($\lesssim 1$~Myr) stellar objects, for which we estimate a total mass of ~ 90 M_sun. The most massive members are three systems containing early B-type stars. Two of them most likely produced very compact HII regions, one of them being still highly embedded and coinciding with a peak seen in emission lines characterising the photon dominated region (PDR). Two of these three systems are resolved binaries, and one of those contains a visible Herbig Be star. An approximate derivation of the mass function of the members of the aggregate gives hints of a slope at high masses shallower than the classical Salpeter slope, and a peak of the mass distribution at a mass higher than that at which the widely adopted log-normal initial mass function peaks. The emission distribution of H$_2$ and Br gamma, tracing the PDR and the ionised gas phase, respectively, suggests that molecular gas is distributed as a shell around the embedded aggregate, filled with centrally-condensed ionised gas. Both, the morphology and the low excitation of the HII region, indicate that the sources of ionisation are the B stars of the embedded aggregate, rather than the external UV field caused by the O stars of Cygnus OB2. The youth of the embedded cluster, combined with the isolation of the globule, suggests that star formation in the globule was triggered by the passage of the ionisation front., Comment: 17 pages, 17 figures, 4 tables, Accepted for publication in A&A

Although mid-to-late type M dwarfs are the most common stars in our stellar neighborhood, our knowledge of these objects is still limited. Open questions include the evolution of their angular momentum, internal structures, dust settling in their atmospheres, age dispersion within populations. In addition, at young ages, late-type Ms have masses below the hydrogen burning limit and therefore are key objects in the debate on the brown dwarf mechanism of formation. In this work we determine and study in detail the physical parameters of two samples of young, late M-type sources belonging to either the Chamaeleon I Dark Cloud or the TW Hydrae Association and compare them with the results obtained in the literature for other young clusters and also for older, field, dwarfs. We used multi-wavelength photometry to construct and analyze SEDs to determine general properties of the photosphere and disk presence. We also used low resolution optical and near-infrared spectroscopy to study activity, accretion, gravity and effective temperature sensitive indicators. We propose a VO-based spectral index that is both temperature and age sensitive. We derived physical parameters using independent techniques confirming the already common feature/problem of the age/luminosity spread. In particular, we highlight two brown dwarfs showing very similar temperatures but clearly different surface gravity (explained invoking extreme early accretion). We also show how, despite large improvement in the dust treatment in theoretical models, there is still room for further progress in the simultaneous reproduction of the optical and near-infrared features of these cold young objects., Comment: 26 pages, 22 figures, 10 tables. Accepted in MNRAS

ESO 243-49 is a high-mass (circular velocity $v_{\rm c}\approx200\,{\rm km\,s^{-1}}$) edge-on S0 galaxy in the Abell 2877 cluster at a distance of $\sim95\,{\rm Mpc}$. To elucidate the origin of its thick disc, we use MUSE science verification data to study its kinematics and stellar populations. The thick disc emits $\sim80\%$ of the light at heights in excess of $3.5^{\prime\prime}$ ($1.6\,{\rm kpc}$). The rotation velocities of its stars lag by $30-40\,{\rm km\,s^{-1}}$ compared to those in the thin disc, which is compatible with the asymmetric drift. The thick disc is found to be more metal-poor than the thin disc, but both discs have old ages. We suggest an internal origin for the thick disc stars in high-mass galaxies. We propose that the thick disc formed either ${\rm a)}$ first in a turbulent phase with a high star formation rate and that a thin disc formed shortly afterwards, or ${\rm b)}$ because of the dynamical heating of a thin pre-existing component. Either way, the star formation in ESO 243-49 was quenched just a few Gyrs after the galaxy was born and the formation of a thin and a thick disc must have occurred before the galaxy stopped forming stars. The formation of the discs was so fast that it could be described as a monolithic collapse where several generations of stars formed in a rapid succession., Comment: Accepted for publication in A&A. The reduced data-cube as well as the data necessary to build the kinematic and stellar population maps are available at https://etsin.avointiede.fi/dataset/urn-nbn-fi-csc-kata20160924142911632378

The radiative feedback of massive stars on molecular clouds creates pillars, globules and other features at the interface between the HII region and molecular cloud. We present here Herschel observations between 70 and 500 micron of the immediate environment of the Cygnus OB2 association, performed within the HOBYS program. All structures were detected based on their appearance at 70 micron, and have been classified as pillars, globules, evaporating gasous globules (EGGs), proplyd-like objects, and condensations. From the 70 and 160 micron flux maps, we derive the local FUV field on the PDR surfaces. In parallel, we use a census of the O-stars to estimate the overall FUV-field, that is 10^3-10^4 G_0 close to the central OB cluster (within 10 pc) and decreases down to a few tens G_0, in a distance of 50 pc. From a SED fit to the four longest Herschel wavelengths, we determine column density and temperature maps and derive masses, volume densities and surface densities for these structures. We find that the morphological classification corresponds to distinct physical properties. Pillars and globules have the longest estimated photoevaporation lifetimes, a few 10^6 yr, while all other features should survive less than that. These lifetimes are consistent with that found in simulations of turbulent, UV-illuminated clouds. We propose a tentative evolutionary scheme in which pillars can evolve into globules, which in turn then evolve into EGGs, condensations and proplyd-like objects., Comment: Received February 18, 2016; accepted April 5, 2016, by A&A

We present a detailed study of the shape of the innermost part of the rotation curves of a sample of 29 nearby spiral galaxies, based on high angular and spectral resolution kinematic Halpha Fabry-Perot observations. In particular, we quantify the steepness of the rotation curve by measuring its slope dRvc(0). We explore the relationship between the inner slope and several galaxy parameters, such as stellar mass, maximum rotational velocity, central surface brightness ({\mu}0), bar strength and bulge-to-total ratio. Even with our limited dynamical range, we find a trend for low-mass galaxies to exhibit shallower rotation curve inner slopes than high-mass galaxies, whereas steep inner slopes are found exclusively in high-mass galaxies. This trend may arise from the relationship between the total stellar mass and the mass of the bulge, which are correlated among them. We find a correlation between the inner slope of the rotation curve and the morphological T-type, complementary to the scaling relation between dRvc(0) and {\mu}0 previously reported in the literature. Although we find that the inner slope increases with the Fourier amplitude A2 and decreases with the bar torque Qb, this may arise from the presence of the bulge implicit in both A2 and Qb. As previously noted in the literature, the more compact the mass in the central parts of a galaxy (more concretely, the presence of a bulge), the steeper the inner slopes. We conclude that the baryonic matter dominates the dynamics in the central parts of our sample galaxies., Comment: 11 pages, 11 figures, accepted for publication in MNRAS

Using 3.6 and 4.5$\mu$m images of 73 late-type, edge-on galaxies from the S$^4$G survey, we compare the richness of the globular cluster populations of these galaxies to those of early type galaxies that we measured previously. In general, the galaxies presented here fill in the distribution for galaxies with lower stellar mass, M$_*$, specifically $\log({\rm M}_*/{\rm M}_\odot) < 10$, overlap the results for early-type galaxies of similar masses, and, by doing so, strengthen the case for a dependence of the number of globular clusters per $10^9\ {\rm M}_\odot$ of galaxy stellar mass, T$_{\rm N}$, on M$_*$. For $8.5 < \log ({\rm M}_*/{\rm M}_\odot) < 10.5$ we find the relationship can be satisfactorily described as T$_{\rm N} = ({\rm M}_*/10^{6.7})^{-0.56}$ when M$_*$ is expressed in solar masses. The functional form of the relationship is only weakly constrained and extrapolation outside this range is not advised. Our late-type galaxies, in contrast to our early-types, do not show the tendency for low mass galaxies to split into two T$_{\rm N}$ families. Using these results and a galaxy stellar mass function from the literature, we calculate that in a volume limited, local Universe sample, clusters are most likely to be found around fairly massive galaxies (M$_* \sim 10^{10.8}$ M$_\odot$) and present a fitting function for the volume number density of clusters as a function of parent galaxy stellar mass. We find no correlation between T$_{\rm N}$ and large-scale environment, but do find a tendency for galaxies of fixed M$_*$ to have larger T$_{\rm N}$ if they have converted a larger proportion of their baryons into stars., Comment: Accepted for publication in ApJ, 13 pages

The formation mechanisms of thick discs are under discussion. Thick discs might have formed either at high redshift on a short time-scale or might have been built slowly over time. They may have an internal or an external origin. Here we study in detail the kinematics and the stellar populations of the thick disc of ESO533-4. ESO533-4 is a nearby bulgeless galaxy. We present the first ever IFU study of an edge-on galaxy with enough depth to study the thick disc. We exposed ESO533-4 with VIMOS@VLT for 6.5hours. The FOV covered an axial extent 0.1-0.7r_25 (1-7kpc). We used pPXF and the MILES library to obtain velocity and stellar population maps. We compared our kinematic data with simple GADGET-2 models. The apparent rotational lag of the thick disc of ESO533-4 is compatible with that expected from the combinations of two effects: differential asymmetric drift and the projection effects arising from studying a disc a few degrees (2-3) away from edge-on. Thus, ESO533-4 contains little or no retrograde material. This is compatible with three formation scenarii: the secular heating of an initially thin disc, the formation of the thick disc at high redshift in a turbulent disc phase, and its creation in a major merger event. If happening in all galaxies, this last mechanism would cause retrograde thick discs in half of them. Retrograde discs have not been observed in the five massive disc galaxies (v_c>120km s^-1) for which thick disc kinematics are known. The populations of the thin and the thick discs are separated in the Age-log(Z/Z_Sun) plane. Thus, the thin and thick discs are made of two distinct stellar populations. Although the stellar population results are not conclusive due to the high dust extinction in ESO533-4, they do not favour a secular evolution origin for the thick disc. Hence, we suggest that the thick disc of ESO533-4 formed in a relatively short event (Abridged)., Comment: Accepted for publication in A&A

We use mid-infrared 3.6 and 4.5microns imaging of NGC 3906 from the Spitzer Survey of Stellar Structure in Galaxies (S4G) to understand the nature of an unusual offset between its stellar bar and the photometric center of an otherwise regular, circular outer stellar disk. We measure an offset of ~720 pc between the center of the stellar bar and photometric center of the stellar disk; the bar center coincides with the kinematic center of the disk determined from previous HI observations. Although the undisturbed shape of the disk suggests that NGC 3906 has not undergone a significant merger event in its recent history, the most plausible explanation for the observed offset is an interaction. Given the relatively isolated nature of NGC 3906 this interaction could be with dark matter sub structure in the galaxy's halo or from a recent interaction with a fast moving neighbor which remains to be identified. Simulations aimed at reproducing the observed offset between the stellar bar / kinematic center of the system and the photometric center of the disk are necessary to confirm this hypothesis and constrain the interaction history of the galaxy., Comment: 18 pages, 8 figures, ApJ in press

The Spitzer Survey of Stellar Structure in Galaxies (S4G) is a volume, magnitude, and size-limited survey of 2352 nearby galaxies with deep imaging at 3.6 and 4.5um. In this paper we describe our surface photometry pipeline and showcase the associated data products that we have released to the community. We also identify the physical mechanisms leading to different levels of central stellar mass concentration for galaxies with the same total stellar mass. Finally, we derive the local stellar mass-size relation at 3.6um for galaxies of different morphologies. Our radial profiles reach stellar mass surface densities below 1 Msun pc-2. Given the negligible impact of dust and the almost constant mass-to-light ratio at these wavelengths, these profiles constitute an accurate inventory of the radial distribution of stellar mass in nearby galaxies. From these profiles we have also derived global properties such as asymptotic magnitudes (and the corresponding stellar masses), isophotal sizes and shapes, and concentration indices. These and other data products from our various pipelines (science-ready mosaics, object masks, 2D image decompositions, and stellar mass maps), can be publicly accessed at IRSA (http://irsa.ipac.caltech.edu/data/SPITZER/S4G/)., Comment: Accepted for publication in ApJS

The aim of our study is to use dynamical simulations to explore the influence of two important dynamical bar parameters, bar strength and bar pattern speed, on the shape of the bar dust lanes. To quantify the shape of the dust lanes we have developed a new systematic method to measure the dust lane curvature. Previous numerical simulations have compared the curvature of bar dust lanes with the bar strength, predicting a relation between both parameters which has been supported by observational studies but with a large spread. We take into account the bar pattern speed to explore, simultaneously, the effect of both parameters on the dust lane shape. To that end, we separate our galactic bars in fast bars $\left(1 < \mathcal{R} < 1.4 \right)$ and slow bars $\left(\mathcal{R} > 1.4 \right)$, obtaining, as previous simulations, an inverse relation between the dust lane curvature and the bar strength for fast bars. For the first time, we extend the study to slow bars, finding a constant curvature as a function of the bar strength. As a result, we conclude that weak bars with straight dust lanes are candidates for slow bars. Finally, we have analysed a pilot sample of ten S$^4$G galaxies, obtaining dust lane curvatures lying within the range covered by the simulations., Comment: 7 pages, 3 figures, accepted for publication in MNRAS 2015 April 7

The Spitzer Survey of Stellar Structure in Galaxies (S$^4$G, Sheth et. al. 2010) is a deep 3.6 and 4.5 $\mu$m imaging survey of 2352 nearby ($< 40$ Mpc) galaxies. We describe the S$^4$G data analysis pipeline 4, which is dedicated to 2-dimensional structural surface brightness decompositions of 3.6 $\mu$m images, using GALFIT3.0 \citep{peng2010}. Besides automatic 1-component S\'ersic fits, and 2-component S\'ersic bulge + exponential disk fits, we present human supervised multi-component decompositions, which include, when judged appropriate, a central point source, bulge, disk, and bar components. Comparison of the fitted parameters indicates that multi-component models are needed to obtain reliable estimates for the bulge S\'ersic index and bulge-to-total light ratio ($B/T$), confirming earlier results \citep{laurikainen2007, gadotti2008, weinzirl2009}. In this first paper, we describe the preparations of input data done for decompositions, give examples of our decomposition strategy, and describe the data products released via IRSA and via our web page ({\bf \tt www.oulu.fi/astronomy/S4G\_PIPELINE4/MAIN}). These products include all the input data and decomposition files in electronic form, making it easy to extend the decompositions to suit specific science purposes. We also provide our IDL-based visualization tools (GALFIDL) developed for displaying/running GALFIT-decompositions, as well as our mask editing procedure (MASK\_EDIT) used in data preparation. In the second paper we will present a detailed analysis of the bulge, disk, and bar parameter derived from multi-component decompositions., Comment: Accepted to ApJS