C Catalán-Torrecilla

Á Castillo-Morales

A Gil De Paz

J Gallego

E Carrasco

J Iglesias-Páramo

R Cedazo

M Chamorro-Cazorla

S Pascual

M L García-Vargas

Á Castillo

N Cardiel

P Gómez-Alvarez

A Pérez-Calpena

I Martínez-Delgado

B T Dullo

P Coelho

G Bruzual

S Charlot

SPATIALLY-RESOLVED ANALYSIS OF NEUTRAL WINDS, STARS AND IONIZED GAS KINEMATICS WITH MEGARA/GTC: NEW INSIGHTS ON THE NEARBY GALAXY UGC 10205

Keywords: galaxies: evolution, galaxies: spiral, galaxies: star formation, galactic winds, techniques: spectroscopic

whether they are published or not. The documents may come

INTRODUCTION

Being a fundamental physical process for the galaxy evolution models, a vast observational evidence have probed the ubiquity of galaxy-scale outflows across all cosmic epochs (see reviews by Veilleux et al. 2005;Rupke 2018). Stellar winds, supernovae or active galactic nuclei (AGNs) have been postulated as the potential mechanisms to drive these phenomena to which we will broadly refer as Galactic Winds (GWs).

GWs are candidates to contribute to the metal enrichment of the intergalactic medium (IGM) (Nath & Trentham 1997;Dalcanton 2007). The gas entrained in the wind may escape the halo potential (Heckman et al. 2000) or it could fall back down on to the galaxy in a process referred to as galactic fountain (Shapiro & Field 1976). GWs are essential for understanding the feedback mechanisms and their impact on the host galaxy. While they might produce a positive effect in which the compression of the insterstellar and circumgalactic medium can potentially boost the star formation (Nayakshin & Zubovas 2012;Zubovas et al. 2013), they can also be considered as a plausible mechanism to suppress star forma-tion within galaxies (Granato et al. 2004;King & Pounds 2015).

Given their intrinsic multiphase nature, different tracers (X-rays, optical nebular emission lines, ISM absorption lines, CO transitions) ought to be used to characterize their different phases (hot, warm, cold, molecular). While most of the studies have traditionally focused on the study of GWs on extreme objects such as luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs), mergers and strong starburst (Martin 2005;Rupke et al. 2005;Rupke & Veilleux 2005;Arribas et al. 2014;Cazzoli et al. 2014Cazzoli et al. , 2016)), the presence of relatively strong winds in less active galaxies have remain mostly unexplored (Castillo-Morales et al. 2007;Jiménez-Vicente et al. 2007;Chen et al. 2010;Roberts-Borsani & Saintonge 2019). Within this scenario, the use of new instrumental capabilities to explore the existence, frequency and strength of winds in relatively quiet galaxies seems the necessary step to follow.

The advent of the Multi-Espectrógrafo en GTC de Alta Resolución para Astronomía (MEGARA, Gil de Paz et al. in prep.) will improve the detectability of the GWs helping to separate them from their hosts with exquisite detail. These high-quality Integral Field Spectroscopy (IFS) data are fundamental to kinematically discriminate both the ionized medium and the cold interstellar gas phase of winds and their geometry. That is why, the characterization of GWs in the central regions of nearby galaxies is one of the main scientific objectives of the MEGARA instrument [START_REF] Carrasco | Society of Photo-Optical Instrumentation Engineers[END_REF].

In this context, UGC 10205 was selected as part of the night-time commissioning observations due to the presence of residual Na I D interestelar absorption previously identified in the CALIFA survey (Sánchez et al. 2012). In the optical range, neutral gas motions are commonly studied using this interstellar medium (ISM) Na I D absorption line doublet at 5889.95 Å and 5895.92 Å. The lower value of the ionization potential for the Na I D (5.14 eV) in comparison with the one for the hydrogen, makes it an ideal candidate to trace the cool (T < 10 4 K) neutral gas.

Besides, by examining the multiple components of the ISM in galaxies such as the case of UGC 10205, where mergers might have played an important role (Reshetnikov & Evstigneeva 1999), one could provide important clues on the impact of these events on their evolution. Past long-slit spectroscopy studies limited to a few directions only (see Vega Beltran et al. 1997) have revealed the complexity of this object. The use of high-resolution IFS data is essential to provide a comprehensive view on the dynamics of this galaxy. In this work, we perform a reconstruction of the stellar kinematics and gas morphology in the central regions of this nearby galaxy. Thus, UGC 10205 is intended to provide a exploratory study to define the final selection for the optimal exploitation of the MEGARA capabilities both in terms of the stellar and gas kinematics.

The structure of this paper is as follows: Section 2 describes the general properties of the UGC 10205 galaxy, Section 3 presents MEGARA observations and data reduction, Section 4 outlines the various steps in our analysis, Section 5 describes the main results of this study including details on the kinematics of the ionized gas and the stellar component and the properties of the cold gas. Finally, we discuss the proposed scenario for this object in Section 6 and we summarize our conclusions in Section 7. Unless otherwise stated, throughout this paper we adopt a a flat ΛCDM cosmology with H 0 = 70 km s -1 Mpc -1 , Ω m = 0.3 and Ω Λ = 0.7.

GENERAL PROPERTIES OF UGC 10205

UGC 10205 is an edge-on disk galaxy classified as an Sa/S0 galaxy by [START_REF] Nilson | Uppsala general catalogue of galaxies Noordermeer, E[END_REF]; Tarenghi et al. (1994) and Nair & Abraham (2010), respectively, with an inclination of i = 84 • [START_REF] Rubin | [END_REF]. A SDSS composite image for this galaxy appears in Figure 1 showing the presence of dust lanes roughly parallel to the major axis of the galaxy. A shell structure in the outer regions of this galaxy was also discovered by the early studies of [START_REF] Rubin | IAU Symposium[END_REF].

The main properties of UGC 10205 are summarized in Table 1 and briefly described below. A systemic heliocentric radial velocity of 6556 km -1 was calculated in Theureau et al. (1998) using measurements of the 21 cm neutral hydrogen line carried out with the meridian transit Nançay radiotelescope. Walcher et al. (2014) used Bruzual & Charlot (2003) stellar population models with a Chabrier (2003) stellar IMF to construct UV to NIR SED to estimate a total stellar mass of log(M /M ) = 10.997 ± 0.111. A set of panchromatic luminosity measurements can be also found in the literature for this galaxy. In particular, a far-infrared luminosity of log(L F IR /L ) = 10.10 was reported by Vlahakis et al. (2005) using the Submillimetre Common-User Bolometer Array (SCUBA) data and a log(L T IR /L ) = 10.35 by Willmer et al. (2009) combining Spitzer Space Telescope observations with SCUBA 850 µm data. UGC 10205 was observed as part of the CALIFA survey (Sánchez et al. 2012). In this context, Catalán-Torrecilla et al. (2015) derived a global L Hα,obs = (10.35 ± 0.57) × 10 40 erg s -1 and a L F U V,obs = (2.71 ± 0.67) × 10 42 erg s -1 from the Galaxy Evolution Explorer (GALEX, Martin 2005).

Questions arise about the complexity of this object as has been already manifested by several authors. The numerical simulations run by Reshetnikov & Evstigneeva (1999) show the capture and tidal disruption of a companion (a small type E/S0). Other authors have considered this galaxy as a possible polar ring candidate (Whitmore et al. 1990;[START_REF] Van Driel | [END_REF]. A profound study done by Vega Beltran et al. (1997) using long-slit spectroscopy data revealed the presence of up to three kinematically distinct gaseous components. The use of IFS data will help to discriminate the spatial distribution of these components and to overcome limitations associated with the use of one-dimensional velocity curves. At a distance of 150 Mpc (H 0 = 75.2 ± 3.3 km s -1 Mpc -1 , Sorce et al. 2014), the central 5.7 × 5.1 kpc 2 of UGC 10205 can be mapped with MEGARA. These observations have the potential to shed some light on the multi-phase components of the ISM to understand the current state and future evolution of this galaxy. The orientation is indicated in the bottom-right corner. Bottom panel: Distribution of the continuum emission of UGC 10205 derived from the MEGARA IFS datacube in the spectral range 5600-5700 Å. The flux units are erg s -1 cm -2 pixel -2 Å-1 . Pixel size is 0.4" in this case. The area of these synthetic spaxels was chosen as the best compromise between the seeing of the observing nights and the routine used to create the datacubes (see more details in Section 4.2).

gral Field Unit (IFU) and a Multi-Object Spectroscopy (MOS) mode. For the purpose of this work, we used only IFU data taken in 2017 June 30th and July 1st as part of the science verification of its commissioning phase. The IFU unit provides a field-of-view (FoV) of 12.5" x 11.3" using 567 hexagonal spaxels of 0.62" in size which makes it ideal for studying the central regions of nearby galaxies. Due to the reasons explained later in Section 4.2, the analysis of the data will be performed using squared spaxels. A recent example showing the potential of MEGARA data can be found in Dullo et al. (2019).

The values of the seeing during the observations ranged between approximately 0.93" and 1.11" as recorded by the differential image motion monitor (DIMM) close to the Telescopio Nazionale Galileo (TNG). Four 1200 s science exposures were combined for the LR-V data set while in the HR-R set-up a total integration time of 2200s were made in two exposures. These data provide a unique view into the cold and warm phases of the interstellar medium. In particular, the LR-V grating is optimal to analyze the ISM Na I λλ 5890,5896 doublet (Na I D) as it yields a wavelength coverage from 5143 Å to 6164 Å with a spectral resolution of FWHM = 0.937 Å at central wavelength λ c = 5695 Å. The HR-R grating provides a wavelength coverage from 6445 Å to 6837 Å (FWHM = 0.355 Å at λ c = 6646 Å) that allow us to characterize the ionized gas by means of the Hα emission line. The reciprocal dispersion is 0.27 Å/pix and 0.097 Å/pix respectively for LR-V and HR-R.

Additional images were acquired to help in the data reduction process. In particular, the spectrophotometric standard star BD+33 2642 was observed and used to create a response function to calibrate the absolute flux scale. Arc calibration lamp frames and twilight sky flat-field/lamp flat-field were obtained at the beginning and at the end of each night, respectively. Finally, the reduction of the raw data comprises the standard procedures of bias subtraction, cosmic-ray removal, flat fielding, tracing and extraction of the spectra, arc calibration solutions, sky-subtraction and flux calibration. It was performed applying the standard MEGARA pipeline (Pascual et al. 2018;Cardiel & Pascual 2018).

ANALYSIS

Galaxy continuum subtraction and interstellar

Na I D detection In this section, we describe the main characteristics of the process followed to model the galaxy continuum for the case of the LR-V grating. As commented before, the LR-V covers the range from 5143 Å to 6164 Å with a spectral resolution of FWHM = 0.937 Å at λ c = 5695 Å.

Firstly, a Voronoi binning algorithm (Cappellari & Copin 2003) is applied to increase the signal-to-noise ratio (S/N) of the spectra in the individual spaxels, reaching a S/N ∼ 20 over the continuum. This process allows a more homogeneous distribution of the S/N within the datacube. Then, we model the galaxy continuum in the spatially-binned MEGARA datacube using a Penalized Pixel-Fitting method (pPXF, Cappellari & Emsellem 2004; Cappellari 2017). The aim of pPXF is to build the stellar template that will best reproduce the observed spectrum in each Voronoi cell. For that purpose, we used the latest theoretical stellar library by Coelho (2014) and P. Coelho, G. Bruzual & S. Charlot in preparation, which represents an expansion of the previous stellar population modelling in Coelho et al. (2005Coelho et al. ( , 2007) ) and Walcher et al. (2009). The resolution of these models, R = 20000, is optimal for all the MEGARA setups and especially for the highest resolution VPHs (R ≈ 20000). We selected a total of 660 models with different ages ranging from 30 Myr to 14 Gyr and 3 different metallicities ([Fe/H] = -1.0,0.0,0.2) at both scaled-solar and α-enhanced mixtures. The stellar continuum was fitted over the spectral range (5150 -6160) Å which contains multiple stellar absorption lines. We carefully masked the ISM features as they are the result of both the stellar and the ISM contribution. Also, the He I emission line at 5875.67 Å was masked due to its proximity to the Na I D line. Finally, the best stellar model was sub-Fig. 2.-Spectrum for the LR-V MEGARA datacube is shown in black while its best fit from pPXF appears in red. The green spectrum at the bottom represents the fit residuals (note the different vertical scale). The blue shaded areas correspond to the masked regions (ISM features and sky line residuals). The inset shows the position of the Voronoi cell corresponding to this spectrum on the same continuum image of the galaxy as in Figure 1.

tracted from the observed spectra to detect the presence of Na I D interstellar absorption. To assess the robustness of the pPXF fits within the region measured and its impact on the results, we have applied a standard bootstrap analysis (see Wu 1986) using 100 simulations to each of the Voronoi spectra. First, a redistribution of the best-fit model residuals is done. Then, we fitted the resampled spectrum using pPXF. For the two first LOSVD moments (V and σ) we obtained dispersions of up to 9 km s -1 and 10 km s -1 after running the 100 repetitions, respectively. We also found mean differences of only 7 km s -1 (5 km s -1 ) between the original velocities (stellar velocity dispersions) values and the ones obtained after applying the bootstrapping technique.

Figure 2 shows a representative fit for one of the Voronoi cells. The position of this particular Voronoi cell within the galaxy can be seen in the inset of the figure. The black line is the observed spectrum while the red line is the fit for the stellar component. The vertical blue shaded regions show the ISM features and sky-lines that have been masked in the fitting process. The residuals of the fitting appear in green color at the bottom. There is a clear detection of ISM Na I D absorption in this spectrum.

The analysis done along this Section allows us (i) to derive fundamental parameters of the galaxy kinematics such as the first four order moments of the line of sight velocity distribution (V, σ, h 3 and h 4 ) for each Voronoi cell of the MEGARA datacube (Section 5.2) and (ii) to explore the presence of interstellar Na I D absorption in UGC 10205 and its connection with neutral gas outflows (Section 5.5).

Recovering the Hα emission

In the following paragraphs we describe the analysis performed to obtain a reliable measure of the Hα emission flux for the case of the HR-R grating. This set-up is characterized by a wavelength coverage from 6445 Å to 6837 Å and a FWHM = 0.355 Å at λ c = 6646 Å. The HR-R observations were taken on different nights, June 30th and July 1st (2017), and the pointing was slightly offset between them (we remind the reader that the data were taken as part of the Commissioning Time). Thus, the first step was to combine both observations to create a new datacube taking into account the relative offsets between them. The second step was to remove the stellar continuum of the galaxy. However, due to the lack of absorption features in the wavelength range covered by the HR-R setup, the use of codes such as pPXF is not recommended. For that reason, the method applied for the HR-R case slightly differs from the one used in the LR-V grating (Section 4.1). The approach followed here to estimate the Hα absorption-corrected emission flux is to derive (i) the equivalent width of the stellar Hα absorption (EW Hα,abs ) and (ii) the stellar velocity corresponding to each spaxel in the datacube by anchoring the stellar populations models to the ones obtained for the LR-V grating. The latest calculation requires to reference the HR-R observations to the LR-V ones. Due to pointing effects, the spatial coverage of the galaxy is slightly different in both set-ups so we have created datacubes with a pixel scale of 0.4 arcsec in each of the spatial axes. The area of these new synthetic spaxels, on which we will perform the subsequent analysis, will therefore be 0.4×0.4 arcsec 2 and will be simply referred as spaxels. The value of the spaxel size was chosen as the best compromise between the seeing of the observing nights and the routine used to create the datacubes. To estimate the EW Hα,abs we used the CALIFA spectrum (3700-7000 Å) integrated in the same area as the one in our MEGARA FoV datacube. To recover the best model of the stellar continuum, we run the pPXF code. Finally, we have obtained the gaussian profile corresponding to the stellar absorption that is added to the Hα emission flux in the HR-R datacube. The main assumption in this approach is that the model that best reproduce the stellar contin-uum is the same for all the spaxels in the MEGARA datacube. As the mean value of the recovered stellar flux at the Hα wavelength is around 9 % of the total Hα flux, we applied this correction assuming that it represents a small second order correction for the estimation of the total Hα flux. As done in the case of the LR-V, we also increase the S/N ratio making a spatial binning of the datacube. In this case, we have defined the S/N as the ratio between the Hα emission line and the nearest continuum and imposed a S/N > 15.

The high spectral resolution of MEGARA allow us to distinguish double and triple-peaked Hα emission line profiles in this galaxy. The finding of up to three kinematically distinct gaseous components is in accordance with the results from Vega Beltran et al. (1997). Some examples showing the complexity of the Hα emission line profiles will be shown and examined in the following sections.

RESULTS

Kinematics of the Emitting Gas as revealed by the

Hα emission line To investigate the gas kinematics, we make use of the Hα emission line. A previous long-slit spectroscopic study of the kinematics of the gas in UGC 10205 was performed by Vega Beltran et al. (1997) using the Intermediate Dispersion Spectrograph (IDS) at the Isaac Newton Telescope (INT) obtaining a mean value of FWHM = 0.86 Å (σ ∼ 17 km s -1 ). These authors found line splitting for three kinematically distinct gaseous components in the position-velocity diagram along the major axis of the galaxy in the inner ± 13". A double-peaked Hα emission appeared in the region -13"≤ r < 3" and a third Hα component is found at distances in the range 3" ≤ r ≤ 13". Later, García-Lorenzo et al. (2015) presented a comprehensive analysis of the gas velocity fields for 177 galaxies among which UGC 10205 is present. The galaxies were observed as part of the CALIFA survey using the Potsdam Multi-Aperture Spectrophotometer (PMAS, Roth et al. 2005) in the PPak mode (Kelz et al. 2006) which provides a nominal resolution (R ∼ 850 at λ ∼ 5000 Å) that is not sufficient to probe the presence of multiple components in this object.

Here, the high-resolution MEGARA IFS data improves the ability to map the kinematics of the line-emitting gas in a spatially-resolved manner, combining the information on the morphology and the kinematics simultaneously. In the next section, we attempt an interpretation for the kinematic properties of the ionized-gas phase by fitting a rotating exponential thin disk model to the observations.

Thin disk model

Both García-Lorenzo et al. (2015) and Vega Beltran et al. (1997) derived monotonically rising rotation velocity curves for the inner regions of UGC 10205. Following these results, we adopt an arctangent parameterization with a smooth turn-over to recover the shape of the rotation curve, V r = (2/π) v max arctan(r/r t ) (Courteau 1997), where v max is the asymptotic maximum rotation speed and r t is the transition radius between the rising and flat part of the rotation curve.

A suite of 2D kinematic models is built to recover the 2D ionized velocity field of the purely rotating galactic disk. The models contain six free parameters: v max , r t , position angle (PA), inclination (i) and the position of the dynamical center (x0,y0). To define the parameters range we use reported measurements and allow variations consistent with our data. The corresponding ranges are: v max = (250 -350) km s -1 , r t = (0.05 -0.5) kpc, i = (20 -89) • , PA = (40 -170) • . For the dynamical center, we use the photometric center as our reference value. The photometric center is defined as the peak of the continuum emission in the region (6540 -6620) Å and we allow a variation of 2 arcsecs around it. A systemic velocity of 6556 km s -1 is taken from NED.

We utilize a minimization based on a mean-squared error to find the best fit between the model and the data. As the Hα channels show a somehow clumpy distribution (see Figure 3), we apply the minimization on the voronois that are more relevant to the shape of the rotation curve. The model parameters that best describe the data are v max = 292 km s -1 , r t = 0.12 kpc, i = 75 • and PA = 130 • . In the case of the inclination, a wealth of values can be found in the literature. García-Lorenzo et al. (2015) derived a photometric inclination of 54 • using the whole galaxy while [START_REF] Rubin | [END_REF] reported i = 84 • . Due to the recent interactions suffered by the galaxy (as suggested in Reshetnikov & Evstigneeva 1999), the gas disk of the outer parts might be mostly affected by them. Thus, the inclination for the innermost parts might be different. We have obtained i = 75 • for these regions. The orientation of the line of nodes is in reasonably good agreement with the value inferred by Barrera-Ballesteros et al. (2015). They derived a morphological PA = (132.4 ± 1.6) • by fitting an ellipse to an isophote at radius 15 arcsec in the r-band SDSS image. They also measured kinematic PA for the ionized component and derived PA = (120.3 ± 2.7) • and PA = (129.4 ± 1.5) • for the approaching and receding sides, respectively.

Figure 3 shows the Hα channel maps of the HR-R datacube. Each map represents the line emission at a particular wavelength channel (from 6705.0 Å to 6709.5 Å) with a channel step of ∆λ = 0.5 Å. From the analysis of this figure, we infer the presence of a highly-inclined gasrich disk located in the central regions of this galaxy. The isovelocity contours are well reproduced by our model for the central parts (see the white solid lines in the channels maps from 6707.0 Å to 6708.5 Å of Figure 3). However, some of the ionized gas do not follow the general rotation pattern as indicated by the isovelocity curves (see green and orange crosses in the channel maps from λ = 6704.5 Å to λ = 6707.5 Å in Figure 3). This component will be described later.

To recover the velocity maps of the ionized gas components, we need to disentangle the different emission lines in the Hα profile. We fit multiple Gaussian components and a polynomial continuum to separate up to three kinematically distinct gaseous components. The corresponding ionized gas kinematic maps are shown in the top panel of Figure 4.

The first component (component 1) is the one associated with the thin-disk model. We trace this component based on the velocity associated to the centroid of each Voronoi cell in the disk model, i. e., we defined the Gaus- There is also a third component (component 3) with a velocity close to the systemic, ranging from -60 km s -1 to 85 km s -1 relative to the systemic. This last component is spatially coincident with the regions previously characterized as deviations from the general rotation pattern (green and orange crosses in Figure 3).

Finally, we select a few characteristic Voronoi cells to show the complexity of the emission line profiles together with the Gaussian fits of the three kinematically distinct components. The position of these Voronoi cells is shown with blue and orange colours in the middle panel of the top row in Figure 5. Blue, red and orange Gaussian fits in the spectra of this figure correspond to components 1, 2 and 3 in Figure 4, respectively. The blue dashed line is the velocity of the thin-disk model while the orange dashed line represents the systemic velocity. Also, a significant rotational broadening for the Hα profiles in the central regions can be seen in this figure (see panel corresponding to voronoi = 3). As we will describe later in Section 5.4, this effect is due to the steep rise of the rotational velocity curve in the innermost regions of this galaxy. We would expect to better constraint the internal dynamics of these regions in nearly face-on or less inclined systems where the σ values are expected to be less broaden by the effect of rotation. In spite of the complex structure that follows the ionized gas distribution, our thin disk model is able to reproduce the kinematic pattern found in this object.

Deviations from the pure thin disk model

As commented before, the distribution of the ionized gas in the low velocity channels (from 6704.5 Å to 6707.5 Å) reveals the presence of a kinematically distinct component that do not follow the rotation pattern of the gas-rich disk. To explain the gas velocity field in those particular regions, two possible interpretations arise.

On one hand, this dynamical component seems to be spatially coincident with the position of the dust lanes that might be driving the motions of the gas (see positions marked as green and orange crosses in Figure 3). Dust lanes are in most cases morphologically associated with a gas dense component even in early-type galaxies (Finkelman et al. 2010(Finkelman et al. , 2012)). Judging from the largescale images available for UGC 10205 (see upper panel in Figure 1 as an example), the well-defined and prominent dust lanes of this object are highly inclined. Due to the reduced size of the MEGARA FoV, we are sampling a relatively small angle on the sky in comparison with the large-scale size of the object (∼ 1 arcmin). Assuming a close to edge-on orientation for the dust lanes, we expect to recover angles (φ) around 90 • and 270 • within this structure. Should that be the case, the difference between the radial velocity component of the gas in each Voronoi cell and the systemic velocity of the thin-disk model (i = 75 • ) should be nearly zero. Possible values for the inclination of the dust lane were investigated to check whether or not the velocities in these regions might be coherent with the structure of a disk model with higher inclination than the main thin disk model (i = 75 • ). An optimal value of i = 83 • was found.

Thus, the bottom panel in Figure 6 shows the distribution of the distances and angles probed by the Voronoi cells in the plane of the dust lanes assuming i = 83 • for this plane. The SE (approaching) side of the dust lane is represented by the orange points while the SW-to-NW (receding) side is marked by the green points. In particular, the points represent the centroid of each Voronoi cell associated with the regions where the line-of-sight passes through the plane of the dust lanes (positions also marked as green and orange crosses in Figure 3).

Additionally, the difference between the velocity of this gas component and the systemic velocity of the thin-disk model (i = 75 • ) computed for each Voronoi cell is shown in the upper panel of Figure 6 as green/orange points. As expected, this difference in velocity (∆v) is nearly zero for φ ∼ 90 • and φ ∼ 270 • . As an example of a particular Voronoi cell, we include the spectrum corresponding to voronoi = 7 (Figure 5). In this spectrum, two different features can be clearly seen: (i) a component almost at the systemic velocity of the model and (ii) the component of the thin-disk model. The former one would be the gas component associated with the dust-lane structure.

At the galactocentric distances corresponding to the positions of the receding side of the dust lane (∼ 30 arcsec), we expect to reach the maximum velocity for the thin disk model (v max = 292 km s -1 , green solid line). i = 83 • is the inclination of the dust-lane that minimizes the distance between the model and ∆v. For the case of 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 3 3 3 9 9 9 9 9 9 2 2 2 2 2 2 2 2 2 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 5 5 5 5 8 8 8 8 8 8 4 4 4 the SE side of the dust lane (orange points), the galactocentric distances are ∼ 10 arcsec so the velocity curve has still not reach its maximum velocity. A velocity of 130 km s -1 (orange dashed line) is the best value consistent with the inclination derived previously. As a conclusion, having a slightly more inclined dust-lane (i = 83 • ) with the same galaxy's potential than the thin disk model (i = 75 • ) seems to be a reasonable good explanation for these regions.

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 3 3 3
On the other hand, the origin of this component could be associated with the merging scenario that creates a localized enhancement of the star formation with a distinct kinematic pattern and the gas might be tracing non-circular motions in those regions. The numerical model presented by Reshetnikov & Evstigneeva (1999) relied on the tidal disruption of a small E/SO object by the massive early-type galaxy. Non-circular gas motions in the nuclear region of UGC 10205 were suggested by these authors. However, the nice match between both the ∆v values and its azimuthal variation with respect to the one of highly-inclined dust lanes seems to favor the former scenario.

Resolved stellar kinematics

The MEGARA absorption-line kinematics for the central regions of the galaxy UGC 10205 are presented in Figure 7. This figure shows both the stellar velocity field and the stellar velocity dispersion maps derived from our pPXF analysis. pPXF yields the values of the stellar mean velocity and stellar velocity dispersion for each Voronoi cell applying a S/N rejection criterion of 20. The reader is referred to Section 4.1 where the main procedure followed to create these maps is explained with detail.

The stellar velocity map (upper panel in Figure 7)

shows a clear stellar rotation pattern with a symmetric velocity field respect to the galaxy center and with a maximum radial velocity of ∼ 6628 km s -1 . Typical errors in our stellar velocity measurements range from 4 to 23 km s -1 with mean errors of 10 km s -1 . To check the accuracy of our kinematic analysis a visual inspection is performed using the results presented in Falcón-Barroso et al. ( 2017). These authors found a similar stellar rotation pattern using the mid-resolution CALIFA data (V1200 with R ∼ 1650) in the wavelength range 3850 -4600 Å.

The stellar velocity dispersion map (bottom panel in Figure 7) displays values ranging from 124 to 230 km s -1 , with values in the nuclear region ∼ 190 km s -1 . These values are already corrected for the instrumental resolution, which in our case is σ ∼ 21 km s -1 . Although no signs of kinematically decoupled components appeared to be present, we have performed the analysis of the Gauss-Hermite moment h3 (skewness of the velocity profile) to examine this effect. The h3 map (not shown) indicates that we can safely discard the presence of an anticorrelation between these values and the stellar velocities ones.

Estimating the asymmetric drift correction

By applying an asymmetric drift correction (ADC) it is possible to account for the effect of the stellar random motions within galaxies. Basically, the mean tangential velocity of a stellar population lags the actual circular velocity as defined by the gravitational potential of the galaxy. ADC accounts for this difference. The effect of the ADC phenomenon will be more significant in those cases where the velocity dispersion of different stellar populations is higher. Although this effect was first noted in the solar neighborhood (Strömberg 1924(Strömberg , 1925)), asymmetric drift corrections have been extensively used in the literature to reconcile the gas and stellar rotation curves observed in external galaxies.

The Stellar Velocity Ellipsoid (SVE) is described in cylindrical coordinates by its radial, tangential and vertical components: σ R , σ φ and σ Z . The shape of the SVE could be parameterized by the axial ratios σ Z /σ R , σ φ /σ R and σ Z /σ φ . For the purpose of this work, we focus our attention on the α = σ Z /σ R and β = σ φ /σ R parameters. In particular, Mogotsi & Romeo (2018) found a global vertical-to-radial velocity dispersion ratio of α = 0.97 ± 0.07 for UGC 10205. Similar values of the radial and vertical dispersions point out to an isotropic heating.

Here, we assume that the stellar azimuthal anisotropy (β φ ) can be approximated using the epicycle anisotropy, i. e., the stellar orbits are considered nearly circular (Gerssen et al. 1997(Gerssen et al. , 2000;;[START_REF] Sersic | Atlas de Galaxias Australes[END_REF]Ciardullo et al. 2004;Westfall et al. 2011;Gerssen & Shapiro Griffin 2012;Blanc et al. 2013;Gentile et al. 2015). Thus, we measure β EA (hereafter β for simplicity) using the following expression:

β = 1 2 ∂ln(v φ ) ∂ln(R) + 1 (1)
where v φ is the tangential speed of the stars. We found a mean value of β = 0.86 ± 0.12. Due to the small dispersion found and in order to simplify the calculations, we assume the previous value for all the Voronoi cells.

Once α and β are known, we infer σ R in each Voronoi cell using the values of σ LOS previously obtained (bottom panel of Figure 7):

σ 2 LOS = σ 2 R sin φ 2 sin i 2 + σ 2 φ cos φ 2 sin i 2 + σ 2 Z cos i 2 (2)
where φ is the angle of the central position of each Voronoi cell in the plane of the galaxy measured from the line of nodes and i refers to the inclination. We assume the same inclination for each Voronoi cell, i. e., in this case we adopt the same inclination for the most inner regions of the galaxy. Due to the agreement between the ionized gas rotation axis and the stellar rotation axis with a difference of (6 ± 3) • , we use an inclination of 75 • for the ADC.

To estimate the ADC in each Voronoi cell, we apply the following expression as reproduced from equation (4-33) of [START_REF] Binney | Galactic dynamics Blanc[END_REF]:

V 2 C = v 2 φ -σ 2 R ∂ln(νR) ∂ln(R) + ∂ln(σ 2 R ) ∂ln(R) + 1 - σ 2 φ σ 2 R + R σ 2 R ∂(vRvZ ) ∂Z (3)
This expression shows the difference between the true velocity (i.e., V C is the circular velocity of a test particle in the potential) and the observed circular velocity. If the SVE is aligned with the cylindrical coordinate system (R,φ,z), the last term in equation 3 could be neglected.

The term ν R in equation 3 is proportional to the stellar mass surface density (Σ) under the assumption that Σ is well traced by the surface brightness profile of the galaxy. In this case, to trace the radial dependence of ν we apply a one dimensional (1D) photometric decomposition over the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS, Chambers et al. 2016) y-band light profile. Among all the Pan-STARRS bands, the near-infrared one is considered the best tracer of the stellar mass as the light is mainly dominated by old MS and RGB stars. Also, problems associated with the interstellar dust (as shown by the presence of prominent dust lanes in UGC 10205) are partly avoided. The fit was performed following the procedure explained in Dullo et al. (2019). We discard non-symmetric features due to bars or spiral arms in the fitting process. The bulge component is parameterized using a Sérsic profile [START_REF] Sersic | Atlas de Galaxias Australes[END_REF]) with a best-fitting Sérsic index n = 1.48 ± 0.22 and a half-light radius r e = 4.10 ± 0.81 arcsec. The galaxy disk component is described with an exponential profile with scale length h = 10.8 ± 1.3 arcsec.

Asymmetric drift correction: comparison with

previous data in the literature Figure 8 shows the results obtained after applying the ADC to the stellar component. For comparison, the predictions for the rotation curves of Sérsic bulges obtained by Noordermeer (2008) are also shown. These curves have been scaled up to take into account the mass and the effective radius of the bulge component. The effective radius, obtained directly from our 1D decomposition, is r e = 4.10 ± 0.81 arcsec (r e ∼ 1.97 kpc). For the stellar mass derivation, we apply the bulge-to-total flux ratio (B/T = 0.19 ± 0.3) obtained as a result of the 1D decomposition and the total stellar mass of the galaxy derived in Walcher et al. (2014) (log(M ) = 10.997 ± 0.111 M ). As the 1D decomposition yields n = 1.48 ± 0.22, we represent the rotation curve for a fiducial bulge with a Sérsic concentration parameter n = 2 as the best approximation in our case. For the same n = 2, blue dashed line curves correspond to variations of the axis ratio q = 0.2 (top) and q = 1 (bottom) in Figure 8. It is expected that the rotation curves from Noordermeer (2008) are a lower limit as they only consider the bulge mass component but our V C measurements also trace the mass in the inner disk of the galaxy.

We have also included in Figure 8 some observational data found in the literature for this galaxy. In particular, the stellar rotation curve derived in Kalinova et al. (2017) is shown in orange. In the new classification proposed by these authors, the rotation curve of UGC 10205 is classified as round-peaked meaning that the circular velocity rises steeply and has a round peak, gradually changing -Stellar (blue circles) rotation curve for UGC 10205. The ionized gas rotation curves (red squares) after applying the ADC under the assumption of an inclination i = 75 • . The light-blue dashed lines are the predictions for the rotation curves of a Sérsic bulge with n = 2 and different values of the axis ratio q = 0.2 (top) and q = 1 (bottom) derived by Noordermeer (2008) and scaled up to match the mass and the effective radius of our bulge component. The best rotation thin disk model is represented by the red solid line. The observed velocity curve values from Vega Beltran et al. (1997) and García-Lorenzo et al. (2015) are shown in dotted gray lines and dotted-dashed green lines, respectively. The velocity curve profile measured by Kalinova et al. (2017) is also shown (orange line).

to the flat part of the circular velocity curve (CVC) at radius ∼ 0.5 R e , where R e is the value they derived for the whole galaxy. They also found that this class is the most common one among early-type galaxies (E4-S0a) and spiral galaxies classified as Sa-Sbc. The typical value expected for the average amplitude of the CVC is around 300 km s -1 , similar to the v max value derived in Section 5.1 for the thin-disk model.

On the other hand, measurements of the ionized gas velocity rotation curve are also included. As the thindisk model derived in Section 5.1 was recovered using mainly the inner regions of the galaxy, we include here other literature data that cover the whole galaxy. The velocity curve based on the Hα + [NII] emission-line data obtained by García-Lorenzo et al. (2015) within CAL-IFA appears as the green dot-dashed lines in Figure 8. The UGC 10205 major-axis kinematic curves derived by Vega Beltran et al. (1997) for the multiple Hα components found on this object using long-slit spectroscopy data are shown by the gray dotted lines. Both García-Lorenzo et al. (2015) and Vega Beltran et al. (1997) derived monotonically rising rotation velocity curves for the inner regions of UGC 10205. These data have been homogenized to i = 75 • . Our stellar data points once corrected by ADC effects are in good agreement with the values obtained by the previous authors up to a galactocentric distance of R ∼ 17 arcsec. For R > 17 arcsec, our measurements match the ones by Vega Beltran et al. (1997) along the NW side of the major axis of the galaxy. These results suggest that assuming an inclination of 75 • for the stellar component in the central part of this galaxy was a reasonable approach. Also, the adoption of a tangential parameterization with a smooth turn-over to recover the shape of the thin disk model in Section 5.1 is justified by the shape of the previous rotation curves.

Average values of the ADC in local disk galaxies are usually 10 % -20 % in late-type systems (see Ciardullo et al. 2004;Merrett et al. 2006;Noordermeer et al. 2008;Herrmann & Ciardullo 2009;Westfall et al. 2011;Martinsson et al. 2013). Higher values for the ADC are expected for E/S0 and early-type spirals. Here, we are analyzing the central part of a S0/Sa galaxy. We found a mean value of V rot (stars)/V rot (gas) = 0.75 ± 0.08. This result is consistent with the one by Cortese et al. (2014). These authors explored only the central parts of a sample mainly composed by early-type spirals finding an average ratio of V rot (stars)/V rot (gas) = 0.75. These results point out the necessity of having high values of the ADC in order to reconcile the discrepancy between the ionized and the stellar velocity rotation curves for these morphological types.

5.5. Nature of the absorbing gas through ISM Na I D measurements Galactic winds are an essential element that determines the evolution of galaxies through cosmic time. In particular, as explained in the introduction, galactic outflows are thought to play a key role in self-regulating star formation in galaxies as well as they constitute necessary mechanisms to explain the chemical evolution and metal enrichment of the circumgalactic medium (Oppenheimer & Davé 2006[START_REF] Nilson | Uppsala general catalogue of galaxies Noordermeer, E[END_REF].

At optical wavelengths, the ISM Na I D absorption doublet is commonly used to trace the properties of neutral atomic gas inflows and outflows in the interstellar medium (Heckman et al. 2000;Schwartz & Martin 2004;Veilleux et al. 2005;Rupke et al. 2005;Martin 2005Martin , 2006;;Chen et al. 2010;Sarzi et al. 2016;Roberts-Borsani & Saintonge 2019). Due to the high inclination of this object the Na I D absorption line profiles are confused with multiple components along the line of sight. Searching for signatures of inflowing/outflowing material in highly inclined systems require having high spatial and spectralresolution observations to break up the absorption profile into different components. Here, the observations provided by MEGARA show the variety of the profiles and the distinct kinematic components for the Na I D excess. Some examples of the ISM absorption Na I D lines that will be explained later (Sections 5.5.1 and 5.5.2) are provided in Figure 9.

For the discussion about inflows/outflows that appear in the following sections, we assume that foreground gas is the dominant source of the absorptions found in the Na I D profiles. In this scenario, blue-shifted absorption is a signature of gas moving towards the observer (outflow) while red-shifted absorption is indicative of gas moving towards the galaxy (inflow). We do not detect reemission of absorbed photons from the background gas. For more details about this topic, the reader is referred to the detailed discussion found in Roberts-Borsani & Saintonge (2019). The analysis of the previous profiles allows us to capture the presence of a red-shifted high-velocity component (see left panel in Figure 9) that might be associated with material moving towards the galaxy center. In particular, the ∆V distribution of the red-shifted narrow Na I D continuum absorption features shows a positive mean value of ∆V = (108.35 ± 18.05) km s -1 relative to the systemic velocity of the galaxy. This red-shifted component is located all over the MEGARA FoV. Also, the small variation of the ∆V distribution among all the Voronoi cells is an indication that the infalling neutral gas extends beyond the spatial coverage of the FoV.

Infalling neutral gas

The previous value suggests a scenario for the ISM Na I D absorption lines in which the red-shifted Na I D interstellar gas is moving inward toward the central regions of UGC 10205. Thus, the absorbing neutral material would be located in a shell-like spatial distribution region moving towards the galaxy as an inflowing gas. We consider this structure to be a shell-like type rather than a stream based mainly in the following two arguments. From a morphological point, there is a complete 2D coverage of the FoV. In the case of the shell a 2D coverage is expected while the stream structure would be seen as a strip of stars or gas. From a kinematical point of view, in the current scenario where the bulk velocity of the inflowing gas presents almost a constant value (approximately 108 km s -1 relative to the systemic velocity of the galaxy) instead of showing a velocity gradient, we consider that the most simple explanation would be to have roughly spherically symmetric gas inflows. The current stage of this galaxy should be highly affected by the previous minor merger suffered in the past as proposed in Reshetnikov & Evstigneeva (1999). The presence of these gas clouds that are inflowing toward the center of the galaxy might be a consequence of such scenario.

Presence of outflow signatures: mass and outflow rate of the wind

Major galaxy mergers and ULIRGS in the local universe exhibit outflows with relatively large velocities and a cone shape (Rupke et al. 2002;Rupke & Veilleux 2005;Westmoquette et al. 2012;Soto et al. 2012;Cazzoli et al. 2014;Perna et al. 2019). Here, we present results from a less active star-forming galaxy. The right panel in Figure 9 shows that clear blue-shifted Na I D absorption signatures are present in this galaxy. The spectra clearly show the presence of double-peaked components in the Na I D absorption line profiles. High spectral resolution is needed to distinct these components. The analysis of the spatially-resolved optical spectroscopy of the blue-shifted NaI D suggests the presence of a wideangle conical outflow emerging from the central regions of UGC 10205 and perpendicular to the galactic plane (Figure 10). The outflow shows a maximum projected velocity of up to -87 km s -1 , blue-shifted with respect to the systemic.

To explore to what extent this neutral gas outflow might remove relatively large amounts of gas from the galaxy, we estimate both the total mass in the wind and the outflow mass rate applying the thin-shell free wind model (equations 13 and 14 in Rupke et al. 2005). Since we have information on the spatial distribution of the wind thanks to the IFS data, we follow the approach of Shih & Rupke (2010); Rupke & Veilleux (2013); Cazzoli et al. (2014Cazzoli et al. ( , 2016) ) using a Voronoi-by-Voronoi basis and assuming a series of thin shells to make it a valid approach for this case. Following the previous prescription, the outflow mass rate is proportional to the deprojected velocity, the HI column density and the deprojected solid angle measured in each Voronoi cell. The prescriptions applied to convert Na column densities into HI column densities are largely affected by the assumptions made on the ionization fraction, the depletion onto dust grains and the Na abundance. To avoid the previous uncertainties in the estimation of the column density of HI atoms, we adopt here the same approach as in Cazzoli et al. (2014Cazzoli et al. ( , 2016)). These authors make use of both the relation between the HI and the color excess proposed by Bohlin et al. (1978) (N HI = 4.8 × E (B-V ) × 10 21 cm -2 ) and the relation between the color excess and the equivalent width as in [START_REF] Turatto | From Twilight to Highlight: The Physics of Supernovae[END_REF] (E (B-V ) = -0.04 + 0.51 × EW N aD ). We calculate a mean value of 1.65 ± 0.07 × 10 21 cm -2 for the HI column density which is in the range of the expected values for low-z galaxies as reported in Heckman et al. (2015) and Martín-Fernández et al. (2016). The cold gas mass outflow rate is 0.78 ± 0.03 M yr -1 and the total mass in the wind is equal to 4.55 ± 0.06 × 10 7 M . These values are consistent with the ones derived using the formalism of the thin-shell model in galaxies with lower star formation rates as in Roberts-Borsani & Saintonge (2019) where the authors used galaxies with SFR in the range 1.45 < SFR [M yr -1 ] < 16.98 obtaining outflow mass rates ranging from 0.15 < Ṁ out [M yr -1 ] < 1.74. A significant mass loss has been postulated to regulate the star formation processes within galaxies. The mass loading factor, expressed as the ratio η = Ṁ out / SFR, is used to describe the strength of the winds. This factor has important consequences for the predictions in the theoretical models with some findings suggesting that η increases with redshift (Barai et al. 2015;Muratov et al. 2015;Sugahara et al. 2017). To calculate the SFR, we have derived dust-corrected Hα growth curves using CALIFA data (following the methodology explained in detail in Catalán-Torrecilla et al. 2017). In the region covered by the MEGARA FoV and applying the conversion SFR(Hα) = 5.5 × 10 -42 L(Hα corr ) (Kennicutt et al. 2009), we obtained SFR(Hα corr ) = 0.49 M yr -1 and η = 1.59. Due to the large amount of dust present in this galaxy, the attenuation correction based on Balmer Decrement measurements might be underestimating the SFR and affecting the real value of η. To mitigate this problem, we make use of hybrid calibrations combin-ing luminosities measured directly with that of the light emitted by dust after being heated by young massive stars (Gordon et al. 2000;Inoue et al. 2001;Iglesias-Páramo et al. 2006;Calzetti et al. 2007;Hao et al. 2011;Kennicutt & Evans 2012;Catalán-Torrecilla et al. 2015). We adopt the Hα obs , FUV obs , 22µm and TIR luminosities calculated for this object in Catalán-Torrecilla et al. (2015) as our best estimations. Using the hybrid SFR prescriptions for the different combinations of Hα + a×IR and FUV + b×IR, we derived a mean SFR of SFR = 1.32 M yr -1 and a loading factor of η = 0.59. As these SFR measurements are for the whole galaxy while the effect of the outflow has been computed at a restricted distance limited to our current FoV, the estimated value of η serves here as a lower limit. Mass loading factors below unity are typically found on galaxies similar to UGC 10205 (see Roberts-Borsani & Saintonge 2019). Large-scale observations that trace the wind extension outside the disk are needed to explore the effect in the neighboring ISM and circumgalactic medium. UGC 10205 was proposed as a possible candidate for being a polar-ring galaxy (Whitmore et al. 1990;[START_REF] Van Driel | [END_REF]. Further observations made with CALIFA revealed that there is not a clear misalignment between the velocity maps of the gas (García-Lorenzo et al. 2015) and the stellar components (Falcón-Barroso et al. 2017). As seen in the 2D stellar kinematics maps of Figure 7, the approaching radial velocities are in the SE side of the galaxy which means the rotation of the disk is counterclockwise.

Dust lanes appeared roughly parallel to the major axis of the galaxy. In particular, the one that is crossing the South region is more clearly visible in front of the bulge (i. e., on the near side of the galaxy) while the one in the North seems to be hidden behind the bulge. This location within the galaxy is consistent with the orientation of the reported outflow, which is seen over the near side of the galaxy. As the cold neutral gas traced by Na I D needs to be observed against the stellar background, the outflow is detected in the approaching cone while the receding part of the cone is not observed. The dust lane that is crossing the galaxy in the South position would prevent us for seeing the outflow component arising in this region which suggest that the dust lane is located on the near side of the galaxy from our position. On the contrary, the dust lane in the North would be located on the far side of the galaxy.

To elucidate the possible mechanism that originated the outflow, it is important to discriminate between different ionization sources in the nuclear region. In Catalán-Torrecilla et al. (2015) we applied the standard Baldwin, Phillips & Terlevich (BPT) diagram (Baldwin et al. 1981) with the demarcation lines of Kauffmann et al. (2003) and Kewley et al. (2001). In particular, we estimated the [OIII]/Hβ and the [NII]/Hα line flux ratio for the central 3" spectra using the CALIFA IFS data. The results pointed out that the center of this galaxy could be classified as a mix SF/AGN region, i. e., the points lied on the composite area of the diagram. This analysis discards the presence of a strong AGN contribution. This is consistent with the studies by Sarzi et al. (2016) where the authors examined a sample composed mostly of massive early-type galaxies suggesting that the cold-gas outflows typically occur in objects dominated either by central starbursts or composite AGN/star formation activity. The spatial distribution of the Hα emission shows that there is a strong resolved Hα contribution in the nuclear region (see the channel maps presented in Figure 3). This might favor a scenario where star formation activity could be the potential mechanism triggering the galactic wind. The presence of the outflow is compatible with being created at the very center of the galaxy where most of the Hα emission is found. Also, AGNdriven outflows are expected to display higher velocities than the ones observed here (Rupke & Veilleux 2013;Rupke et al. 2017). Although further investigations need to be done, it seems plausible to discard the AGN-driven as the potential mechanism that generate the outflow in this case.

Finally, the analysis of the ionized and the cold phase of the ISM shows that they have different locations within the galaxy. The Hα emission is confined to a disk structure supported by rotation while the Na I D is associated with (i) a shell-like structure moving infalling towards the galaxy center and (ii) a blue-shifted Na I D absorption in the near side of the galaxy that reflects the presence of a galactic outflow along the NE semi-minor axis. The ionized gas in the central region shares the kinematic pattern of the thin-disk model proposed to explain the motion of the gas. Thus, the analysis of the ionized Hα emission line shows the presence of three different dynamical components, with one of them associated with a highly-inclined (i = 75 • ) rotating disk. Also, a second dynamical component associated with the position of the dust lanes is present. The gas is rotating in this dusty structure surrounding the central regions with a higher inclination, i = 83 • .

SUMMARY AND CONCLUSIONS

In this paper, we present high resolution optical IFS observations of the nearby galaxy UGC 10205, made with MEGARA on the GTC 10.4m telescope at La Palma Observatory. The combination of the LR-V and HR-R set-ups provides a wealth of information on the state of the multi-phase interstellar medium and the stellar kinematics in UGC 10205 and allows us to study the kinematics of the cool and the ionized gas in detail. MEGARA IFS observations have shown an intriguing and complex structure in the central regions of UGC 10205. The main results obtained in this study can be summarized as the following:

• We present fundamental parameters of the galaxy kinematics by means of creating 2D stellar kinematics maps for the first two order moments of the line of sight velocity distribution (V and σ).

• In an attempt to reconciliate the ionized and the stellar velocity rotation curves, ADC corrections are applied. Having high quality IFS data is fundamental to properly estimate the ADC in each Voronoi cell and to recover the shape of the SVE. A mean value of V rot (stars)/V rot (gas) = 0.75 ± 0.08 is found. Applying this correction, there is a significant agreement between the derived stellar velocity rotation curve and the gaseous ionized Hα one previously derived in the literature by both Vega Beltran et al. (1997) using long-slit spectroscopy data and García-Lorenzo et al. (2015) with the CALIFA data. These results suggest that the innermost regions of early-type spirals needed large values of the asymmetric drift correction as proposed by Cortese et al. (2014).

• We also explore the presence of interstellar Na I λλ 5890,5896 doublet absorption that is commonly used to trace the properties of neutral gas inflows/outflows. The complexity of the neutral gas kinematics is revealed by the presence of multiple narrow components. An overall shift of ∆V = (114 ± 5) km s -1 for the foreground redshifted component of the neutral gas is detected through the Na I D absorption. This suggests the presence of mild inflows towards the inner regions. The inflowing gas detection might be triggered by a past minor merger suffered by the galaxy (Reshetnikov & Evstigneeva 1999).

• The characterization of outflows is essential to place constraints on the models of galaxy formation and evolution. The majority of works in the literature has focused on understanding the complexity of outflows in extreme objects such as ULIRGS. Here, we point out the necessity to extend these studies including less active galaxies. To elucidate the role that this outflow might play in removing gas from UGC 10205, we calculate physical properties of the cold gas entrained in the wind applying the formalism of the thinshell model. We obtained a total mass in the wind M out = 4.55 ± 0.06 × 10 7 M and a cold gas mass outflow rate Ṁ out = 0.78 ± 0.03 M yr -1 . The calculation of the mass loading factor yield values of 0.59 < η < 1.59, although more accuracy SFR measurements are needed to constraint this value as well as large-scale observations to explore the effect of the wind in the neighboring ISM and circumgalactic medium.

• Finally, assuming that the gas describes circular orbits in the galaxy's plane, we model the Hα gas distribution in the inner regions using a thin-disk model with (i = 75 • ). The primary goal of our modeling is to characterize the kinematics and the ionized gas geometry. Despite the complexity of the Hα emission line profile where up to three kinematically distinct gaseous components appeared, our simple thin-disk model still can reproduce most of the kinematical features seen in the central regions of UGC 10205. Also, a distinct kinematic component of the gas with a higher inclination (i = 83 • ) pattern is detected in a position coincident with the dust lane structures.

Finally, we emphasize the potential of MEGARA to make meaningful comparisons between the neutral and the ionized gas phase components in the central regions of nearby galaxies.

UGC 10205 is a remarkable example of an object that has revealed an active current phase, reinforcing the im-portance of spatially-resolved high-resolution analyses in less active galaxies. A systematic exploration of galactic winds in normal star-forming objects and in its multiple phases is required to provide a complete picture on the role that these phenomena might play to shape galaxies and its influence on the surrounding environments.

Authors acknowledge financial support from the Spanish Ministry of Economy and Competitiveness (MINECO) under grant number AYA2016-75808-R which is partly funded by the European Regional Development Fund (ERDF), AYA2017-90589-REDT, RTI2018-096188-B-I00, S2018/NMT-4291 and ESP2017-83197P projects. C.C.-T. gratefully acknowledges the support of the European Youth Employment Initiative (YEI) by means of the Postdoctoral Fellowship Program. GB acknowledges financial support from DGAPA-UNAM through PAPIIT project IG100319.

  [START_REF] Nilson | Uppsala general catalogue of galaxies Noordermeer, E[END_REF], (2): Tarenghi et al. (1994), (3): Nair & Abraham (2010) , (4): NASA/IPAC Extra-galactic Database , (5): Rubin et al. (1985), (6): Walcher et al. (2014), (7): Willmer et al. (2009), (8): Catalán-Torrecilla et al. (2015),

  Fig. 1.-Top panel: SDSS composite image of UGC 10205.Overlaid on the image is the footprint of the MEGARA IFU field-of-view (red rectangle) covering a region of 5.1 × 5.7 kpc 2 . The orientation is indicated in the bottom-right corner. Bottom panel: Distribution of the continuum emission of UGC 10205 derived from the MEGARA IFS datacube in the spectral range 5600-5700 Å. The flux units are erg s -1 cm -2 pixel -2 Å-1 . Pixel size is 0.4" in this case. The area of these synthetic spaxels was chosen as the best compromise between the seeing of the observing nights and the routine used to create the datacubes (see more details in Section 4.2).

  Fig. 3.-Hα velocity channel maps for the central part of the galaxy UGC 10205 in the velocity range 6363.14 -6797.40 km s -1 . A velocity bin of 22.86 km s -1 corresponding to a ∆λ = 0.5 Å is used to slice the datacube. The wavelength of each channel is shown in the lower right corner. White lines represent the isovelocity curves for a thin-disk model with inclination of 75 • . Green and orange crosses represent the positions of the Voronoi cells that are best reproduced by a model with inclination of 83 • as discussed in Section 5.1.2. Continuum flux contours from the LR-V datacube are shown as black solid lines.

  Fig. 4.-Top panel: velocity maps of the three kinematically distinct gaseous components: (a) component defined by the thin-disk model (arctangent function), (b) second component with a similar rotation pattern than the previous one, (c) component with a velocity close to the systemic value. These components are defined by the blue, red and orange Gaussian fits in the spectra of Figure 5, respectively. Continuum flux contours from the LR-V datacube are shown as gray solid lines. Bottom panel: (d) velocity map for the ionized gas component defined by the thin-disk model, (e) the best thin-disk model for the previous component, (f) the residuals after subtracting the best thin-disk model.

Fig. 5 .

 5 Fig. 5.-Left top panel: Hα flux contours together with the rotating thin disk model derived for the Hα velocity field are shown as blue, grey and red colours for ranges [6703.5 -6704.5] Å, [6707.0 -6708.0] Å and [6710.5 -6711.5] Å, respectively. Middle top panel: Isovelocity contours for the thin disk model (i = 75 • )showing the location of some specific Voronoi cells within the datacube. The color coding assigned to the Voronoi cells is chosen for convenience, allowing a better visual inspection of their spatial extent within the map. Rest panels: Examples showing the presence of kinematically distinct gaseous Hα components for the Voronoi cells indicated in the middle top panel. The zero velocity is set at the systemic redshift of the galaxy (orange dashed line) while the velocity of the thin-disk model for each particular Voronoi cell is shown by the blue dashed line. Gaussian fits represent the component associated with the thin-disk model (blue), the component with a similar rotation pattern to the previous one but restricted to a smaller region (red) and the component with a velocity close to the systemic (orange), respectively. The units of the horizontal axis are in km s -1 relative to the systemic velocity of the galaxy.

  Fig. 6.-Top panel: Velocity distribution in the plane of the dust lanes (i = 83 • ). Orange (approaching side) and green (receding side) points represent the difference between the radial velocity component of the gas and the systemic velocity of the thin-disk model (i = 75 • ) computed for each Voronoi cell (∆v) within the dust lane. Green solid line and orange dashed line are the velocity models with vmax = 292 km s -1 and vmax = 130 km s -1 , respectively. Bottom panel: Distribution of the galactocentric distances and angles for the different Voronoi cells. Orange and green points represent the values of the ones located in the dust lane.

  Fig. 7.-Top panel: Stellar velocity field for UGC 10205. Note that the galaxy systemic velocity is reported to be 6556 km s -1 . Bottom panel: Stellar velocity dispersion. Contours of the continuum emission (black) are overlaid on both maps. Color scales and units are indicated on the right part of each map.

  Fig.8.-Stellar (blue circles) rotation curve for UGC 10205. The ionized gas rotation curves (red squares) after applying the ADC under the assumption of an inclination i = 75 • . The light-blue dashed lines are the predictions for the rotation curves of a Sérsic bulge with n = 2 and different values of the axis ratio q = 0.2 (top) and q = 1 (bottom) derived byNoordermeer (2008) and scaled up to match the mass and the effective radius of our bulge component. The best rotation thin disk model is represented by the red solid line. The observed velocity curve values from VegaBeltran et al. (1997) andGarcía-Lorenzo et al. (2015) are shown in dotted gray lines and dotted-dashed green lines, respectively. The velocity curve profile measured byKalinova et al. (2017) is also shown (orange line).

Fig. 9 .

 9 Fig. 9.-Examples of the ISM absorption Na I D line profiles for UGC 10205. The insets show the position of the Voronoi cells used to obtain each spectrum. The vertical dashed lines represent the wavelengths corresponding to the Na I λλ 5890, 5896 at the systemic velocity of the galaxy. In the right panel, the thick blue line shows the best fit to the individual narrow Gaussian components for the blue-shifted component of the Na I D doublet that traces the outflow. The red line is the fit to the red-shifted Na I D absorption signatures that appears as an indicator of gas moving towards the galaxy as an infalling gas. In the left panel only the presence of the red-shifted Na I D absorption could be seen.

  Fig. 10.-Neutral gas kinematics derived for the blue-shifted absorption component of the Na I D doublet with respect to the systemic velocity of the galaxy as inferred from the MEGARA data. Contours indicate the distribution of the continuum flux from the LR-V datacube.

  6. DISCUSSION: PROPOSED SCENARIO TO UNDERSTAND THE COMPLEX KINEMATICS IN UGC 10205