Galaxy NGC Research Articles

Uncovering the geometry of the hot X-ray corona in the Seyfert galaxy NGC 4151 with IXPE

ABSTRACT We present an X-ray spectropolarimetric analysis of the bright Seyfert galaxy NGC 4151. The source has been observed with the Imaging X-ray Polarimetry Explorer (IXPE) for 700 ks, complemented with simultaneous XMM–Newton (50 ks) and NuSTAR (100 ks) pointings. A polarization degree Π = 4.9 ± 1.1 per cent and angle Ψ = 86° ± 7° east of north (68 per cent confidence level) are measured in the 2–8 keV energy range. The spectropolarimetric analysis shows that the polarization could be entirely due to reflection. Given the low reflection flux in the IXPE band, this requires, however, a reflection with a very large (>38 per cent) polarization degree. Assuming more reasonable values, a polarization degree of the hot corona ranging from ∼4 to ∼8 per cent is found. The observed polarization degree excludes a ‘spherical’ lamppost geometry for the corona, suggesting instead a slab-like geometry, possibly a wedge, as determined via Monte Carlo simulations. This is further confirmed by the X-ray polarization angle, which coincides with the direction of the extended radio emission in this source, supposed to match the disc axis. NGC 4151 is the first active galactic nucleus with an X-ray polarization measure for the corona, illustrating the capabilities of X-ray polarimetry and IXPE in unveiling its geometry.

Lightning: An X-Ray to Submillimeter Galaxy SED-fitting Code with Physically Motivated Stellar, Dust, and AGN Models

We present an updated version of Lightning, a galaxy spectral energy distribution (SED) fitting code that can model X-ray to submillimeter observations. The models in Lightning include the options to contain contributions from stellar populations, dust attenuation and emission, and active galactic nuclei (AGNs). X-ray emission, when utilized, can be modeled as originating from stellar compact binary populations with the option to include emission from AGNs. We have also included a variety of algorithms to fit the models to observations and sample parameter posteriors; these include adaptive Markov Chain Monte Carlo (MCMC), affine-invariant MCMC, and Levenberg–Marquardt gradient-descent (MPFIT) algorithms. To demonstrate some of the capabilities of Lightning, we present several examples using a variety of observational data. These examples include (1) deriving the spatially resolved stellar properties of the nearby galaxy M81, (2) demonstrating how X-ray emission can provide constraints on the properties of the supermassive black hole of a distant AGN, (3) exploring how to rectify the attenuation effects of inclination on the derived the star formation rate of the edge-on galaxy NGC 4631, (4) comparing the performance of Lightning to similar Bayesian SED-fitting codes when deriving physical properties of the star-forming galaxy NGC 628, and (5) comparing the derived X-ray and UV-to-IR AGN properties from Lightning and CIGALE for a distant AGN. Lightning is an open-source application developed in IDL and is available at https://github.com/rafaeleufrasio/lightning.

NGC 1068 constraints on neutrino-dark matter scattering

The IceCube collaboration has observed the first steady-state point source of high-energy neutrinos, coming from the active galaxy NGC 1068. If neutrinos interacted strongly enough with dark matter, the emitted neutrinos would have been impeded by the dense spike of dark matter surrounding the supermassive black hole at the galactic center, which powers the emission. We derive a stringent upper limit on the scattering cross section between neutrinos and dark matter based on the observed events and theoretical models of the dark matter spike. The bound can be stronger than that obtained by the single IceCube neutrino event from the blazar TXS 0506+056 for some spike models.

Star Formation in the Dwarf Seyfert Galaxy NGC 4395: Evidence for Both AGN and SN Feedback?

We present a detailed multiwavelength study of star formation in the dwarf galaxy NGC 4395, which hosts an active galactic nucleus (AGN). From our observations with the Ultra-Violet Imaging Telescope, we have compiled a catalog of 284 star-forming (SF) regions, out of which we could detect 120 SF regions in Hα observations. Across the entire galaxy, we found the extinction-corrected star formation rate (SFR) in the far ultraviolet (UV) to range from 2.0 × 10−5 M ⊙ yr−1 to 1.5 × 10−2 M ⊙ yr−1 with a median of 3.0 × 10−4 M ⊙ yr−1, and the age to lie in the range of ∼1 to 98 Myr with a median of 14 Myr. In Hα we found the SFR to range from 7.2 × 10−6 M ⊙ yr−1 to 2.7 × 10−2 M ⊙ yr−1 with a median of 1.7 × 10−4 M ⊙ yr−1, and the age to lie between 3 to 6 Myr with a median of 5 Myr. The stellar ages derived from Hα show a gradual decline with galactocentric distance. We found three SF regions close to the center of NGC 4395 with high SFRs both from Hα and UV, which could be attributed to feedback effects from the AGN. We also found six other SF regions in one of the spiral arms having higher SFRs. These are very close to supernova remnants, which could have enhanced the SFR locally. We obtained a specific SFR (SFR per unit mass) for the whole galaxy of 4.64 × 10−10 yr−1.

The Isotropic Center of NGC 5419—A Core in Formation?

With its cored surface brightness profile, the elliptical galaxy NGC 5419 appears as a typical high-mass early-type galaxy (ETG). However, the galaxy hosts two distinct nuclei in its center. We use high-signal MUSE (Multi-unit Spectroscopic Explorer (Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO program 099.B-0193(A).)) spectral observations and novel triaxial dynamical orbit models to reveal a surprisingly isotropic central orbit distribution in NGC 5419. Recent collisionless simulations of merging massive ETGs suggest a two-phase core formation model, in which the low-density stellar core forms rapidly by supermassive black holes (SMBHs) sinking into the center due to dynamical friction. Only afterwards do the SMBHs form a hard binary, and the black hole scouring process slowly changes the central orbit distribution from isotropic to tangential. The observed cored density profile, the double nucleus, and the isotropic center of NGC 5419 together thus point to an intermediate evolutionary state where the first phase of core formation has taken place, yet the scouring process is only beginning. This implies that the double nucleus is an SMBH binary. Our triaxial dynamical models indicate a total mass of the two SMBHs in the center of NGC 5419 of M BH = (1.0 ± 0.08) × 1010 M ⊙. Moreover, we find that NGC 5419's complex kinematically distinct core can be explained by a coherent flip of the direction of orbital rotation of stars on tube orbits at ∼3 kpc distance from the galaxy center together with projection effects. This is also in agreement with merger simulations hosting SMBHs in the same mass regime.

The Effect of Superpositions on the Planetary Nebula Luminosity Function

Planetary nebula (PN) surveys in systems beyond ∼10 Mpc often find high-excitation, point-like sources with [O iii] λ5007 fluxes greater than the apparent bright-end cutoff of the planetary nebula luminosity function (PNLF). Here we identify PN superpositions as one likely cause for the phenomenon and describe the proper procedures for deriving PNLF distances when object blends are a possibility. We apply our technique to two objects: a model Virgo-distance elliptical galaxy observed through a narrowband interference filter, and the Fornax lenticular galaxy NGC 1380 surveyed with the MUSE integral-field unit spectrograph. Our analyses show that even when the most likely distance to a galaxy is unaffected by the possible presence of PN superpositions, the resultant value will still be biased toward too small a distance due to the asymmetrical nature of the error bars. We discuss the future of the PNLF in an era where current ground-based instrumentation can push the technique to distances beyond ∼35 Mpc.

The Lively Accretion Disk in NGC 2992. III. Tentative Evidence of Rapid Ultrafast Outflow Variability

We report on the 2019 XMM-Newton+NuSTAR monitoring campaign of the Seyfert galaxy NGC 2992, observed at one of its highest flux levels in the X-rays. The time-averaged spectra of the two XMM-Newton orbits show ultrafast outflows (UFOs) absorbing structures above 9 keV with >3σ significance. A detailed investigation of the temporal evolution on a ∼5 ks timescale reveals UFO absorption lines at a confidence level >95% (2σ) in 8 out of 50 XMM-Newton segments, estimated via Monte Carlo simulations. We observe a wind variability corresponding to a length scale of 5 Schwarzschild radii r S . Adopting the novel Wind in the Ionized Nuclear Environment model, we estimate the outflowing gas velocity (v = 0.21–0.45c), column density (N H = 4–8 × 1024 cm−2) and ionization state (), taking into account geometrical and special relativity corrections. These parameters lead to instantaneous mass outflow rates of M ⊙ yr−1, with associated outflow momentum rates L Bol/c and kinetic energy rates L Bol. We estimate a wind duty cycle of ≈12% and a total mechanical power of ≈2 times the active galactic nuclei (AGN) bolometric luminosity, suggesting that the wind may drive significant feedback effects between the AGN and the host galaxy. Notably, we also provide an estimate for the wind launching radius and density of ≈5r S , 1011 cm−3, respectively.

The fragility of thin discs in galaxies – II. Thin discs as tracers of the assembly history of galaxies

ABSTRACT Thin galactic discs and nuclear stellar discs (NSDs) are fragile structures that can be easily disturbed by merger events. By studying the age of the stellar populations in present-day discs, we can learn about the assembly history of galaxies and place constraints on their past merger events. Following on the steps of our initial work, we explore the fragility of such disc structures in intermediate-mass-ratio dry encounters using the previously constructed N-body model of the Fornax galaxy NGC 1381 (FCC 170), which hosts both a thin galactic disc and an NSD. We dismiss major and minor encounters, as the former were previously shown to easily destroy thin-disc structures, whereas the latter take several Hubble times to complete in the specific case of FCC 170. The kinematics and structure of the thin galactic disc are dramatically altered by the mergers, whereas the NSD shows a remarkable resilience, exhibiting only a smooth increase of its size when compared to the model evolved in isolation. Our results suggest that thin galactic discs are better tracers for intermediate-mass-ratio mergers, while NSDs may be more useful for major encounters. Based on our simulations and previous analysis of the stellar populations, we concluded that FCC 170 has not experienced any intermediate-mass-ratio dry encounters for at least ∼10 Gyr, as indicated by the age of its thin-disc stellar populations.

Very long baseline interferometry imaging of H2O maser emission in the nearby radio galaxy NGC 4261

Abstract We report dual-frequency very long baseline interferometry (VLBI) observations at 22 and 43 GHz toward the nucleus of a nearby radio galaxy NGC 4261. In particular, we present a VLBI image of the 22 GHz H2O maser line and its location in the circumnuclear region of NGC 4261. H2O maser emission is marginally detected above the three times the image rms level at channels within a velocity range of approximately 2250–2450 km s−1, slightly red-shifted with respect to the systemic velocity. H2O maser emission is located approximately 1 milliarcsecond east of the brightest continuum component at 22 GHz, where the continuum spectrum is optically thick, that is, at the free–free absorbed receding jet by ionized gas. A positional coincidence between H2O maser emission and an ionized gas disk implies that the H2O maser emission arises from the near side of the disk, amplifying continuum emission from the background receding jet. If the disk axis is oriented 64° relative to the line of sight, the H2O maser emission is expected to be at a mean radius of 0.3 pc in the disk. The broad line width of the H2O maser emission can be attributed to complex kinematics in the immediate vicinity of the supermassive black hole (SMBH), including ongoing gas infall on to the SMBH, turbulence, and outflow. This is analogous to the multi-phase circumnuclear torus model in the nearest radio-loud H2O megamaser source NGC 1052. An alternative explanation for H2O maser association is the shock region between the jet and ambient molecular clouds. However, this explanation fails to describe the explicit association of H2O maser emission only with the free–free absorbed receding jet.

The period–luminosity relation for Mira variables in the Milky Way using Gaia DR3: a further distance anchor for H0

ABSTRACT Gaia DR3 parallaxes are used to calibrate preliminary period–luminosity relations of O-rich Mira variables in the 2MASS J, H, and Ks bands using a probabilistic model accounting for variations in the parallax zero-point and underestimation of the parallax uncertainties. The derived relations are compared to those measured for the Large and Small Magellanic Clouds, the Sagittarius dwarf spheroidal galaxy, globular cluster members, and the subset of Milky Way Mira variables with VLBI parallaxes. The Milky Way linear JHKs relations are slightly steeper and thus fainter at short period than the corresponding LMC relations, suggesting population effects in the near-infrared are perhaps larger than previous observational works have claimed. Models of the Gaia astrometry for the Mira variables suggest that, despite the intrinsic photocentre wobble and use of mean photometry in the astrometric solution of the current data reduction, the recovered parallaxes should be on average unbiased but with underestimated uncertainties for the nearest stars. The recommended Gaia EDR3 parallax zero-point corrections evaluated at $\nu _\mathrm{eff}=1.25\, \mu \mathrm{m}^{-1}$ require minimal ($\lesssim 5\, \mu \mathrm{as}$) corrections for redder five-parameter sources, but overcorrect the parallaxes for redder six-parameter sources, and the parallax uncertainties are underestimated at most by a factor ∼1.6 at $G\approx 12.5\, \mathrm{mag}$. The derived period–luminosity relations are used as anchors for the Mira variables in the Type Ia host galaxy NGC 1559 to find $H_0=(73.7\pm 4.4)\, \mathrm{km\, s}^{-1}\, \mathrm{Mpc}^{-1}$.

Fuelling the nuclear ring of NGC 1097

ABSTRACT Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar dust lanes, which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes. Measuring the mass inflow rates is therefore important to understanding the mass/energy budget and evolution of galactic nuclei. In this work, we use CO datacubes from the PHANGS-ALMA survey and a simple geometrical method to measure the bar-driven mass inflow rate on to the nuclear ring of the barred galaxy NGC 1097. The method assumes that the gas velocity in the bar lanes is parallel to the lanes in the frame co-rotating with the bar, and allows one to derive the inflow rates from sufficiently sensitive and resolved position–position–velocity diagrams if the bar pattern speed and galaxy orientations are known. We find an inflow rate of $\dot{M}=(3.0 \pm 2.1)\, \rm M_\odot \, yr^{-1}$ averaged over a time span of 40 Myr, which varies by a factor of a few over time-scales of ∼10 Myr. Most of the inflow appears to be consumed by star formation in the ring, which is currently occurring at a star formation rate (SFR) of $\simeq\!1.8\!-\!2 \, \rm M_\odot \, yr^{-1}$, suggesting that the inflow is causally controlling the SFR in the ring as a function of time.

NGC 1436: the making of a lenticular galaxy in the Fornax Cluster

ABSTRACT We study the evolutionary path of the Fornax Cluster galaxy NGC 1436, which is known to be currently transitioning from a spiral into a lenticular morphology. This galaxy hosts an inner star-forming disc and an outer quiescent disc, and we analyse data from the MeerKAT Fornax Survey, ALMA, and the Fornax 3D survey to study the interstellar medium and the stellar populations of both disc components. Thanks to the combination of high resolution and sensitivity of the MeerKAT data, we find that the H i is entirely confined within the inner star-forming disc, and that its kinematics is coincident with that of the CO. The cold gas disc is now well settled, which suggests that the galaxy has not been affected by any environmental interactions in the last ∼1 Gyr. The star-formation history derived from the Fornax 3D data shows that both the inner and outer discs experienced a burst of star formation ∼5 Gyr ago, followed by rapid quenching in the outer disc and by slow quenching in the inner disc, which continues forming stars to this day. We claim that NGC 1436 has begun to effectively interact with the cluster environment 5 Gyr ago, when a combination of gravitational and hydrodynamical interactions caused the temporary enhancement of the star-formation rate. Furthermore, due to the weaker gravitational binding H i was stripped from the outer disc, causing its rapid quenching. At the same time, accretion of gas onto the inner disc stopped, causing slow quenching in this region.

Exploring the Magnetic Field Geometry in NGC 891 with SOFIA/HAWC+

Stratospheric Observatory For Infrared Astronomy/High-resolution Airborne Wideband Camera-plus 154 μm Far-Infrared polarimetry observations of the well-studied edge-on galaxy NGC 891 are analyzed and compared to simple disk models with ordered (planar) and turbulent magnetic fields. The overall low magnitude and the narrow dispersion of fractional polarization observed in the disk require significant turbulence and a large number of turbulent decorrelation cells along the line of sight through the plane. Higher surface brightness regions along the major axis to either side of the nucleus show a further reduction in polarization and are consistent with a view tangent to a spiral feature in our disk models. The nucleus also has a similar low polarization, and this is inconsistent with our model spiral galaxy where the ordered magnetic field component would be nearly perpendicular to the line of sight through the nucleus on an edge-on view. A model with a barred spiral morphology with a magnetic field geometry derived from radio synchrotron observations of face-on barred spirals fits the data much better. There is clear evidence for a vertical field extending into the halo from one location in the disk coincident with a polarization null point seen in near-infrared polarimetry, probably due to a blowout caused by star formation. Although our observations were capable of detecting a vertical magnetic field geometry elsewhere in the halo, no clear signature was found. A reduced polarization due to a mix of planar and vertical fields in the dusty regions of the halo best explains our observations, but unusually significant turbulence cannot be ruled out.

WISDOM project – XIV. SMBH mass in the early-type galaxies NGC 0612, NGC 1574, and NGC 4261 from CO dynamical modelling

ABSTRACT We present a CO dynamical estimate of the mass of the supermassive black hole (SMBH) in three nearby early-type galaxies: NGC 0612, NGC 1574, and NGC 4261. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3–6 observations of the 12CO(2–1) emission line with spatial resolutions of 14–58 pc (0.01″–0.26″). We detect disc-like CO distributions on scales from ≲ 200 pc (NGC 1574 and NGC 4261) to ≈10 kpc (NGC 0612). In NGC 0612 and NGC 1574 the bulk of the gas is regularly rotating. The data also provide evidence for the presence of a massive dark object at the centre of NGC 1574, allowing us to obtain the first measure of its mass, MBH = (1.0 ± 0.2) × 108 M⊙ (1σ uncertainty). In NGC 4261, the CO kinematics is clearly dominated by the SMBH gravitational influence, allowing us to determine an accurate black hole mass of (1.62 ± 0.04) × 109 M⊙ (1σ uncertainty). This is fully consistent with a previous CO dynamical estimate obtained using a different modelling technique. Signs of non-circular gas motions (likely outflow) are also identified in the inner regions of NGC 4261. In NGC 0612, we are only able to obtain a (conservative) upper limit of MBH ≲ 3.2 × 109 M⊙. This has likely to be ascribed to the presence of a central CO hole (with a radius much larger than that of the SMBH sphere of influence), combined with the inability of obtaining a robust prediction for the CO velocity curve. The three SMBH mass estimates are overall in agreement with predictions from the MBH − σ* relation.

WISDOM Project – XV. Giant molecular clouds in the central region of the barred spiral galaxy NGC 5806

ABSTRACT We present high spatial resolution (≈24 pc) Atacama Large Millimeter/sub-millimeter Array 12CO(2-1) observations of the central region of the nearby barred spiral galaxy NGC 5806. NGC 5806 has a highly structured molecular gas distribution with a clear nucleus, a nuclear ring, and offset dust lanes. We identify 170 spatially and spectrally resolved giant molecular clouds (GMCs). These clouds have comparable sizes (Rc) and larger gas masses, observed linewidths (σobs, los), and gas mass surface densities than those of clouds in the Milky Way disc. The size–linewidth relation of the clouds is one of the steepest reported so far ($\sigma _{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{1.20}$ ), the clouds are on average only marginally bound (with a mean virial parameter 〈αvir〉 ≈ 2), and high velocity dispersions are observed in the nuclear ring. These behaviours are likely due to bar-driven gas shocks and inflows along the offset dust lanes, and we infer an inflow velocity of ≈120 km s−1 and a total molecular gas mass inflow rate of ≈5 M⊙ yr−1 into the nuclear ring. The observed internal velocity gradients of the clouds are consistent with internal turbulence. The number of clouds in the nuclear ring decreases with azimuthal angle downstream from the dust lanes without clear variation of cloud properties. This is likely due to the estimated short lifetime of the clouds (≈6 Myr), which appears to be mainly regulated by cloud–cloud collision and/or shear processes. Overall, it thus seems that the presence of the large-scale bar and gas inflows to the centre of NGC 5806 affect cloud properties.

Regular Chains of Star Formation Regions in Spiral Arms and Rings of Disk Galaxies

The regularity in the distribution of young stellar groups along the spiral arms of galaxies, first discovered by Bruce and Debra Elmegreen in 1983, was considered a rather rare phenomenon. However, recent studies of the spatial regularities in the distribution of the young stellar populations along the arms of the spiral galaxies NGC 628, NGC 895, NGC 4321, NGC 5474, NGC 6946, as well as along the rings of the spiral galaxy NGC 6217 and the lenticular galaxy NGC 4324, have revealed that this spatial (quasi) regularity and/or the presence of regular chains of star-forming regions is a fairly common phenomenon. Across all galaxies, the characteristic regularity scale is 350–500 pc or a multiple thereof. It should be noted that theoretical models predict an instability scale of a stellar-gas disk on the order of a few kpc, which is several times larger than what has been observed. The paper is partly based on the report presented at the Modern Stellar Astronomy 2022 Conference held at the Caucasian Mountain Observatory of the Sternberg Astronomical Institute, Moscow State University, on November 8–10, 2022.

Counter-Rotating Gaseous Disk and Star Formation in the S0 Galaxy NGC 934

Counter-Rotating Gaseous Disk and Star Formation in the S0 Galaxy NGC 934

The Mass of the Black Hole in NGC 5273 from Stellar Dynamical Modeling

We present a new constraint on the mass of the black hole in the active S0 galaxy NGC 5273. Due to the proximity of the galaxy at 16.6 ± 2.1 Mpc, we were able to resolve and extract the bulk motions of stars near the central black hole using adaptive-optics-assisted observations with the Gemini Near-infrared Integral Field Spectrograph, as well as constrain the large-scale kinematics using archival Spectroscopic Areal Unit for Research and Optical Nebulae spectroscopy. High-resolution Hubble Space Telescope imaging allowed us to generate a surface-brightness decomposition, determine approximate mass-to-light ratios for the bulge and disk, and obtain an estimate for the disk inclination. We constructed an extensive library of dynamical models using the Schwarzschild orbit-superposition code FORSTAND, exploring a range of disk and bulge shapes, halo masses, etc. We determined a black hole mass of M • = [0.5–2] × 107 M ⊙, where the low side of the range is in agreement with the reverberation mapping measurement of M • = [4.7 ± 1.6] × 106 M ⊙. NGC 5273 is one of the few nearby galaxies that hosts a broad-lined active galactic nucleus, allowing a crucial comparison of black hole masses derived from independent mass-measurement techniques.

A Luminous Dust-obscured Tidal Disruption Event Candidate in a Star-forming Galaxy at 42 Mpc

While the vast majority of tidal disruption events (TDEs) have been identified by wide-field sky surveys in the optical and X-ray bands, recent studies indicate that a considerable fraction of TDEs may be dust obscured and thus preferentially detected in the infrared (IR) wave bands. In this Letter, we present the discovery of a luminous mid-IR nuclear flare (termed WTP14adbjsh), identified in a systematic transient search of archival images from the NEOWISE mid-IR survey. The source reached a peak luminosity of L ≃ 1043 erg s−1 at 4.6 μm in 2015 before fading in the IR with a TDE-like F ∝ t −5/3 decline, radiating a total of more than 3 × 1051 erg in the last 7 yr. The transient event took place in the nearby galaxy NGC 7392, at a distance of around 42 Mpc; yet, no optical or X-ray flare is detected. We interpret the transient as the nearest TDE candidate detected in the last decade, which was missed at other wavelengths due to dust obscuration, hinting at the existence of TDEs that have been historically overlooked. Unlike most previously detected TDEs, the transient was discovered in a star-forming galaxy, corroborating earlier suggestions that dust obscuration suppresses significantly the detection of TDEs in these environments. Our results demonstrate that the study of IR-detected TDEs is critical in order to obtain a complete understanding of the physics of TDEs and to conclude whether TDEs occur preferentially in a particular class of galaxies.

Integral Field Spectroscopy of the cometary starburst galaxy NGC 4861

ABSTRACT Using the PMAS Integral Field Unit on the Calar Alto 3.5 m telescope, we observed the southern component (Markarian 59) of the ‘cometary’ starburst galaxy NGC 4861. Mrk 59 is centred on a giant nebula and concentration of stars 1 kpc in diameter. Strong $\rm H\,\alpha$ emission points to a star-formation rate (SFR) at least 0.47 $\rm M_{\odot }\, yr^{-1}$. Mrk 59 has a very high [O iii] $\rm \lambda 5007/H\,\beta$ ratio, reaching 7.35 in the central nebula, with a second peak at a star-forming hotspot further north. Fast outflows are not detected but nebular motion and galaxy rotation produce relative velocities up to 40 km $\rm s^{-1}$. Spectral analysis of different regions with ‘Fitting Analysis using Differential evolution Optimization’ (FADO) finds that the stars in the central and ‘spur’ nebulae are very young, $\rm \le 125~Myr$ with a large $\rm \lt 10~Myr$ contribution. Older stars ($\rm \sim 1~Gyr$) make up the northern disk component, while the other regions show mixtures of 1 Gyr age with very young stars. This and the high specific SFR $\rm \sim 3.5~Gyr^{-1}$ imply a bimodal star formation history, with Mrk 59 formed in ongoing starbursts fuelled by a huge gas inflow, turning the galaxy into an asymmetric ‘green pea’ or blue compact dwarf. We map the He ii λ4686 emission, and identify a broad component from the central nebula, consistent with the emission of ∼300 Wolf–Rayet stars. About a third of the He ii λ4686 flux is a narrow line emitted from a more extended area covering the central and spur nebulae, and may have a different origin.