# The evolution of superluminous supernova LSQ14mo and its interacting host galaxy system

T.-W. Chen, M. Nicholl, S.J. Smartt, P.A. Mazzali, R.M. Yates, T.J. Moriya, C. Inserra, N. Langer, T. Krühler, Y.-C. Pan, et al.

► To cite this version:

T.-W. Chen, M. Nicholl, S.J. Smartt, P.A. Mazzali, R.M. Yates, et al. The evolution of superluminous supernova LSQ14mo and its interacting host galaxy system. Astronomy & Astrophysics - A&A, 2017, 602, pp.A9. 10.1051/0004-6361/201630163. hal-01554685

HAL Id: hal-01554685

https://hal.science/hal-01554685v1

Submitted on 14 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The evolution of superluminous supernova LSQ14mo and its interacting host galaxy system*

T.-W. Chen<sup>1,2</sup>, M. Nicholl<sup>3</sup>, S. J. Smartt<sup>4</sup>, P. A. Mazzali<sup>5,6</sup>, R. M. Yates<sup>1</sup>, T. J. Moriya<sup>7</sup>, C. Inserra<sup>4</sup>, N. Langer<sup>2</sup>, T. Krühler<sup>1</sup>, Y.-C. Pan<sup>8</sup>, R. Kotak<sup>4</sup>, L. Galbany<sup>9,10</sup>, P. Schady<sup>1</sup>, P. Wiseman<sup>1</sup>, J. Greiner<sup>1</sup>, S. Schulze<sup>11</sup>, A. W. S. Man<sup>12</sup>, A. Jerkstrand<sup>6</sup>, K. W. Smith<sup>4</sup>, M. Dennefeld<sup>13</sup>, C. Baltay<sup>14</sup>, J. Bolmer<sup>1,15</sup>, E. Kankare<sup>4</sup>, F. Knust<sup>1</sup>, K. Maguire<sup>4</sup>, D. Rabinowitz<sup>14</sup>, S. Rostami<sup>14</sup>, M. Sullivan<sup>16</sup>, and D. R. Young<sup>4</sup>

<sup>1</sup> Max-Planck-Institut für Extraterrestrische Physik, Giessenbachstraße 1, 85748 Garching, Germany
e-mail: jchen@mpe.mpg.de

<sup>2</sup> Argelander Institute for Astronomy, University of Bonn, Auf dem Hügel 71, 53121 Bonn, Germany

<sup>3</sup> Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA

<sup>4</sup> Astrophysics Research Centre, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, UK

<sup>5</sup> Astrophysics Research Institute, Liverpool John Moores University, IC2, Liverpool Science Park, 146 Brownlow Hill, Liverpool L3 5RF, UK

<sup>6</sup> Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, 85748 Garching-bei-München, Germany

<sup>7</sup> Division of Theoretical Astronomy, National Astronomical Observatory of Japan, National Institutes of Natural Sciences, 2-21-1 Osawa, Mitaka, 181-8588 Tokyo, Japan

<sup>8</sup> Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064, USA

<sup>9</sup> Pittsburgh Particle Physics, Astrophysics, and Cosmology Center (PITT PACC), 3941 Pittsburgh, USA

<sup>10</sup> Physics and Astronomy Department, University of Pittsburgh, Pittsburgh, PA 15260, USA

<sup>11</sup> Department of Particle Physics and Astrophysics, Weizmann Institute of Science, 76100 Rehovot, Israel

<sup>12</sup> ESO, Karl-Schwarzschild-Strasse 2, 85748 Garching-bei-München, Germany

<sup>13</sup> Institut d'Astrophysique de Paris, CNRS, and Université Pierre et Marie Curie, 98bis boulevard Arago, 75014 Paris, France

<sup>14</sup> Physics Department, Yale University, 217 Prospect Street, New Haven, CT 06511-8499, USA

<sup>15</sup> European Southern Observatory, Alonso de Córdova 3107, Vitacura, Casilla 19001, Santiago 19, Chile

<sup>16</sup> Department of Physics and Astronomy, University of Southampton, Southampton, SO17 1BJ, UK

Received 29 November 2016 / Accepted 25 January 2017

## ABSTRACT

We present and analyse an extensive dataset of the superluminous supernova (SLSN) LSQ14mo ($z = 0.256$), consisting of a multi-colour light curve from $-30\,d$ to $+70\,d$ in the rest-frame (relative to maximum light) and a series of six spectra from PESSTO covering $-7\,d$ to $+50\,d$. This is among the densest spectroscopic coverage, and best-constrained rising light curve, for a fast-declining hydrogen-poor SLSN. The bolometric light curve can be reproduced with a millisecond magnetar model with $\sim 4\,M_{\odot}$ ejecta mass, and the temperature and velocity evolution is also suggestive of a magnetar as the power source. Spectral modelling indicates that the SN ejected $\sim 6\,M_{\odot}$ of CO-rich material with a kinetic energy of $\sim 7 \times 10^{51}$ erg, and suggests a partially thermalised additional source of luminosity between $-2\,d$ and $+22\,d$. This may be due to interaction with a shell of material originating from pre-explosion mass loss. We further present a detailed analysis of the host galaxy system of LSQ14mo. PESSTO and GROND imaging show three spatially resolved bright regions, and we used the VLT and FORS2 to obtain a deep (five-hour exposure) spectra of the SN position and the three star-forming regions, which are at a similar redshift. The FORS2 spectrum at $+300$ days shows no trace of SN emission lines and we place limits on the strength of [OI] from comparisons with other Ic supernovae. The deep spectra provides a unique chance to investigate spatial variations in the host star-formation activity and metallicity. The specific star-formation rate is similar in all three components, as is the presence of a young stellar population. However, the position of LSQ14mo exhibits a lower metallicity, with $12 + \log(O/H) = 8.2$ in both the $R_{23}$ and N2 scales (corresponding to $\sim 0.3\,Z_{\odot}$). We propose that the three bright regions in the host system are interacting, which could induce gas flows triggering star formation in low-metallicity regions.

**Key words.** supernovae: general – supernovae: individual: LSQ14mo – galaxies: abundances – galaxies: dwarf – galaxies: interactions

## 1. Introduction

Superluminous supernovae (SLSNe) are 10–100 times brighter than normal core-collapse SNe (CCSNe) and can reach absolute magnitudes of $\sim -21$ (see Gal-Yam 2012 for a review). They are divided into two distinct groups from the optical

spectral features around the peak brightness. The first is that of SLSNe, which do not generally exhibit hydrogen or helium lines in spectra and show no spectral signatures of interaction between fast moving ejecta and circumstellar shells, and thus are classified as SLSNe type I or type Ic (e.g. Pastorello et al. 2010; Quimby et al. 2011; Chomiuk et al. 2011; Inserra et al. 2013; Howell et al. 2013; Nicholl et al. 2013). The light curves of SLSNe type I span a wide range of rise ($\sim 15–40\,d$) and

<footnote>* Based on observations at ESO, Program IDs: 191.D-0935, 094.D-0645, 096.A-9099.</footnote>
decline timescales (~30–100 d). These may form two separate subclasses of slowly- and fast-declining objects (with a paucity of events at the midpoints of these ranges), but, with the current small sample size, they are also consistent with a continuous distribution (see Nicholl et al. 2015a). The second group is that of hydrogen-rich SLSNe (SLSNe II). A small number show broad Hα features during the photospheric phase and do not exhibit any clear sign of interaction in their early spectra (see Inserra et al. 2017). These intrinsically differ from the strongly interacting SNe such as SN 2006gy (e.g. Smith et al. 2007; Ofek et al. 2007), which are also referred to as SLSNe IIn. Hereafter in this paper, we only address the SLSN I types.

Three main competing theoretical models have been proposed to explain the extreme luminosity of SLSNe of type I or Ic. These are (i) a central engine scenario, such as millisecond magnetar spin-down (e.g. Woosley 2010; Kasen & Bildsten 2010; Dessart et al. 2012) or black hole accretion (e.g. Dexter & Kasen 2013); (ii) pair-instability SNe (PISNe; e.g. Heger & Woosley 2002) for the broad light curve SLSNe; and (iii) the interaction of SN ejecta with dense and massive circumstellar medium (CSM) shells (e.g. Chevalier & Irwin 2011; Chatzopoulos et al. 2012; Sorokina et al. 2016).

Light curve modelling is an important diagnostic for current SLSNe studies. Firstly, the millisecond magnetar model can reproduce a wide range of light curves, fitting both fast and slow decliners (e.g. Inserra et al. 2013; Nicholl et al. 2013); Chen et al. (2015) found that the slowly-fading late-time light curve also fits a magnetar model well if the escape of high-energy gamma rays is taken into account (for time-varying leakage see Wang et al. 2015). In contrast, the PISN model can only explain slowly-fading light curve events, such as SN 2007bi, which was initially suggested to be a PISN (Gal-Yam et al. 2009) based on a $^{56}$Co decay-like tail. However, the relatively rapid rise time of similar SLSNe PTF12dam and PS1-11ap (Nicholl et al. 2013; McCrum et al. 2014) argues against this interpretation. Recently, Kozyreva et al. (2017) demonstrated that radiative-transfer simulations have essential scatter depending on the input ingredients, a relatively fast rising time was found for the helium PISN model (Kasen et al. 2011). Finally, the CSM model is very flexible due to having many free parameters (e.g. Nicholl et al. 2014). However, the mechanism that could generate the dense and massive CSM (>few $M_\odot$) inferred from light curve fitting is still unclear. The hydrogen detected in the late-time spectra of iPTF13ehe was interpreted as CSM interaction features (Yan et al. 2015) but Moriya et al. (2015) showed that interaction with a binary component was equally plausible.

From the spectral point of view, the ionic line identifications and investigation of the excitation processes have been recently provided by spectral modelling of Mazzali et al. (2016) of early, photospheric phase spectra. The late-time nebular spectra allows us to investigate the composition of the SN ejecta. Dessart et al. (2013), Jerkstrand et al. (2016, 2017) simulated the nebular spectra expected in PISN models, which show little emission in the blue region (<6000 Å) of the optical spectrum because this region is fully blocked by the optically thick ejecta. In the red, strong [Ca II] λ7300 and [Fe I] lines are expected along with neutral silicon and sulphur in the near infra-red. These model nebular spectra look markedly different from that of several slow-fading SLSNe of Jerkstrand et al. (2017) providing no support for a PISN interpretation. Instead, the similarity of the late-time nebular spectra of a number of energetic SNe Ic (e.g. SN 1998bw) and those of SLSNe has been identified by Jerkstrand et al. (2017), Nicholl et al. (2016a). Jerkstrand et al. (2017) find that at least 10 $M_\odot$ of oxygen is required in the ejecta

to produce the strong [O I] lines and that the material must be highly clumped. Moreover, Inserra et al. (2016) presented the first spectropolarimetric observations of a SLSN, showing that SN 2015bn has an axisymmetric geometry, which is similar to those SNe that are connected with long-duration gamma-ray bursts (LGRBs; e.g. Galama et al. 1998).

The study of the host galaxies of SLSNe provide a strong constraint for understanding the stellar progenitors of SLSNe given that the distances of SLSNe ($0.1 < z < 4$; Cooke et al. 2012) are too far away to detect their progenitors directly. The host galaxies of SLSNe are generally faint dwarf galaxies (Neill et al. 2011) which tend to have low metallicity and low mass (Stoll et al. 2011; Chen et al. 2013; Lunnan et al. 2013). They share properties with LGRB host galaxies (Lunnan et al. 2014; Japelj et al. 2016), although SLSN host galaxies appear more extreme (Vreeswijk et al. 2014; Chen et al. 2015; Leloudas et al. 2015b; Angus et al. 2016). Leloudas et al. (2015b) found that the equivalent widths of [O III] are much higher in SLSN hosts than in LGRB hosts, which they argue may imply that the stellar population is, on average, younger in SLSN than in LGRB host galaxies, and by implication that the progenitors of SLSNe are more massive than GRBs. However, this is not corroborated by the *Hubble* Space Telescope study of hosts by Lunnan et al. (2015). There is no significant difference in host environments of fast-declining and slowly-fading SLSNe (Chen et al. 2015; Leloudas et al. 2015b; Perley et al. 2016). Recently, Perley et al. (2016) and Chen et al. (2016) both suggested that a sub-solar metallicity is required for SLSN progenitors.

The SLSN LSQ14mo was discovered by the La Silla-QUEST (LSQ; Baltay et al. 2013) on 2014 January 30, located at RA = 10:22:41.53, Dec = −16:55:14.4 (J2000.0). The Public ESO Spectroscopic Survey of Transient Objects (PESSTO; Smartt et al. 2015) took a spectrum on 2014 January 31 and identified it as a SLSN Ic (Leloudas et al. 2014). The spectrum showed a blue continuum with strong O II absorption features around 4200–4600 Å, mimicking PTF09cnd at a phase of ~1 week before maximum light. Narrow Mg II λλ 2796, 2803 ISM absorption lines indicated a redshift of $z = 0.253$. However, these lines were not resolved in the low-resolution spectrum and thus it only provided a rough redshift measurement (or the ISM is slightly in the foreground with respect to the SN location). In this paper, we re-identify the redshift of LSQ14mo to be $z = 0.2563$ using the narrow host emission lines (more details see Sect. 2.3). Observations of LSQ14mo were first presented by Nicholl et al. (2015a), who included the *r*-band and pseudo bolo-metric light curve in their statistical sample of SLSNe. Furthermore, in Leloudas et al. (2015a), as part of the first polarimetric study of a SLSN, they found there is no evidence for significant deviations from spherical symmetry of LSQ14mo.

This paper is organised as follows: in Sect. 2, we give details on the photometric follow-up and spectroscopic observations of LSQ14mo and its host galaxy. We construct the bolometric light curve of LSQ14mo and fit models in Sect. 3, and here we also compare it to the light curves and spectral evolution of the other SLSNe. In Sect. 4, we apply the new spectral synthesis model from Mazzali et al. (2016) for LSQ14mo. In Sect. 5, we present the various host galaxy properties. In the discussion in Sects. 6 and 7, we search for other SLSNe that have spectral features indicative of a thin shell interaction; we argue that the entire host system is an interacting dwarf galaxy. Finally, we conclude in Sect. 8. The appendix contains log and magnitude tables and provides more details on data reduction. In this paper, we use a cosmology of $H_0 = 70$, $\Omega_{\Lambda} = 0.73$, $\Omega_M = 0.27$. All magnitudes
<table><caption>**Table 1.** Log of spectroscopic observations of LSQ14mo and its host galaxy.</caption><thead><tr><th>Date</th><th>MJD</th><th>Phase (day)</th><th>Telescope</th><th>Instrument</th><th>Grism</th><th>Exp. time (s)</th><th>Slit (′′)</th><th>Resolution (Å)</th><th>Range (Å)</th></tr></thead><tbody><tr><td rowspan="2">2014 Jan. 31</td><td rowspan="2">56 688.16</td><td rowspan="2">−7.0</td><td rowspan="2">NTT</td><td rowspan="2">EFOSC2</td><td>Gr#13</td><td>1800 × 1</td><td>1.0</td><td>18.2</td><td>3668–9269</td></tr><tr><td>Gr#11</td><td>1800 × 2</td><td>1.0</td><td>13.8</td><td>3357–7486</td></tr><tr><td rowspan="2">2014 Feb. 6</td><td rowspan="2">56 694.19</td><td rowspan="2">−2.2</td><td rowspan="2">NTT</td><td rowspan="2">EFOSC2</td><td>Gr#11</td><td>1800 × 3</td><td>1.0</td><td>13.8</td><td>3357–7486</td></tr><tr><td>Gr#16</td><td>1800 × 3</td><td>1.0</td><td>13.4</td><td>6008–10 008</td></tr><tr><td>2014 Feb. 20</td><td>56 708.06</td><td>8.8</td><td>NTT</td><td>EFOSC2</td><td>Gr#13</td><td>1800 × 3</td><td>1.0</td><td>18.2</td><td>3668–9269</td></tr><tr><td>2014 Feb. 28</td><td>56 716.13</td><td>15.2</td><td>NTT</td><td>EFOSC2</td><td>Gr#13</td><td>1800 × 3</td><td>1.0</td><td>18.2</td><td>3668–9269</td></tr><tr><td>2014 Mar. 8</td><td>56 724.26</td><td>21.7</td><td>NTT</td><td>EFOSC2</td><td>Gr#13</td><td>2400 × 3</td><td>1.0</td><td>18.2</td><td>3668–9269</td></tr><tr><td>2014 Apr. 21</td><td>56 768.04</td><td>56.5</td><td>NTT</td><td>EFOSC2</td><td>Gr#13</td><td>2400 × 4</td><td>1.0</td><td>18.2</td><td>3668–9269</td></tr><tr><td>2015 Jan. 23</td><td>57 045.18</td><td>277.2</td><td>VLT</td><td>FORS2</td><td>GRIS_300V</td><td>1230 × 10</td><td>1.0</td><td>10.2</td><td>4300–9600</td></tr><tr><td>2015 Feb. 21</td><td>57 074.28</td><td>300.4</td><td>VLT</td><td>FORS2</td><td>GRIS_300V</td><td>1230 × 3</td><td>1.0</td><td>10.2</td><td>4300–9600</td></tr><tr><td>2015 Feb. 24</td><td>57 077.27</td><td>302.8</td><td>VLT</td><td>FORS2</td><td>GRIS_300V</td><td>1230 × 2</td><td>1.0</td><td>10.2</td><td>4300–9600</td></tr></tbody></table>

**Notes.** The phase (day) has been corrected for time dilation ($z = 0.256$) and relative to the SN r-band maximum on MJD 56 697. Resolution of NTT + EFOSC2 is adopted from Smartt et al. (2015), as measured from the sky lines. For VLT + FORS2 with a filter GG435+81 and a 1′′0 slit, the resolving power is 440 at the central wavelength 5900 Å, and the dispersion is 3.27 Å/pixel.

reported are in the AB system. We assumed a Chabrier (2003) initial mass function (IMF) of the host galaxy.

## 2. Observational data

### 2.1. Supernova data

The multi-band optical light curves of LSQ14mo were obtained using several telescope and instrumental configurations listed in Tables C.1. These are the QUEST instrument on the ESO 1.0-m Schmidt telescope dedicated to the LSQ survey; the 2.0-m Liverpool Telescope (LT; Steele et al. 2004) using the IO:O imager; the 6.5-m Magellan telescope and the Inamori-Magellan Areal Camera & Spectrograph (IMACS); the 3.58-m New Technology Telescope (NTT) using the ESO Faint Object Spectrograph and Camera (EFOSC2) in the framework of the PESSTO program (Smartt et al. 2015). We discuss the reduction and photometric calibration of these images in detail in the Appendix.

Additionally, we have ultraviolet data from the Ultraviolet and Optical Telescope (UVOT) on the Swift satellite, with (uvw2, uvm2, uvw1, and u) observations; see Poole et al. (2008) covering the spectral range ~1600–3800 Å. Those observations were acquired only within a few days of maximum light. The data were reduced using the HEASARC<fnref n="1" /> software. All measured magnitudes and detection limits are listed in Tables C.1 and C.2.

We obtained a series of spectra of LSQ14mo from −7.0d to +56.5d using NTT+EFOSC2 within the PESSTO program. The observational log is reported in Table 1. The NTT + EFOSC2 data were reduced in a standard fashion using the PESSTO pipeline (Smartt et al. 2015). This applies bias-subtraction and flat-fielding with halogen lamp frames. The spectra were extracted using the pipeline and wavelength-calibrated by identifying lines of He and Ar lamps, and flux-calibrated using sensitivity curves obtained using spectroscopic standard star observations on the same nights. The final data products can be found in the ESO Science Archive Facility as part of PESSTO SSDR2², and all spectra will also be available through WISEREP<fnref n="3" /> (Yaron & Gal-Yam 2012), and along with photometry through the Open SN Catalog<fnref n="4" /> (Guillochon et al. 2017).

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Left: VLT + FORS2 slit angle overlaid on NTT + EFOSC2 (PESSTO) i-band image of the host galaxy system taken on 2015 February 12 (+293.1 d). We have labelled the LSQ14mo position (C), bright H II region (B) and PSO J155.6730-16.9216 (A). North is up and east is left. Right: portion of the 2D spectrum from VLT + FORS2 showing the emission lines Hβ and [O III] λλ4959,5007. Three spatially resolved components, A, B, and C, also offset in velocity, are clearly seen.</figcaption></figure>

### 2.2. Host galaxy photometry

We gathered deep, high-spatial-resolution images of LSQ14mo and its host galaxy system after +250 d with PESSTO (see Sect. A). We chose the best-seeing (0.8′′) i-band image taken at +293.1 d to determine that there are three components in the host galaxy system PSO J155.6730-16.9216. We have named these positions A, B, and C from south to north, respectively (see Fig. 1) and use these names throughout this paper. The SN is located coincident with the flux at position C. We selected appropriate aperture sizes for positions A, B, and C to limit the contamination from each corresponding component (2.5′′, 0.75′′, and 1′′ for positions A, B, and C, respectively). Aperture photometry was then carried out within IRAF/DAOPHOT, and we used the same aperture size to measure the flux of local secondary standards in the field (see Sect. A) to set the zero-point. We found $i = 20.98 \pm 0.04$ mag for component A and $23.75 \pm 0.10$ mag and $23.44 \pm 0.10$ mag for components B and C, respectively. We obtained pure host images on 2015 December 17 (+539.2 d) using the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND; Greiner et al. 2008). We did not detect a 4000 Å break for positions B and C, expected to have $g - r > 0.3$ mag from the synthetic photometry of the spectra. All photometry measurements are listed in Table 2.

The brightest component, position A, was detected in Pan-STARRS1 $3\pi$ images as PSO J155.6730-16.9216. Those stack

<footnote n="1">¹ NASA's High Energy Astrophysics Science Archive Research Center.</footnote>

<footnote n="2">² Data access described on www.pessto.org</footnote>

<footnote n="3">³ http://wiserep.weizmann.ac.il</footnote>

<footnote n="4">⁴ https://sne.space</footnote>
<table><caption>**Table 2.** Photometry of each component of the host galaxy system of LSQ14mo using GROND.</caption><thead><tr><th>Position (aperture)</th><th>g</th><th>r</th><th>i</th><th>z</th><th>J</th><th>H</th><th>K</th></tr></thead><tbody><tr><td>A (1.9'')</td><td>21.85 (0.15)</td><td>21.27 (0.02)</td><td>21.11 (0.05)</td><td>21.07 (0.04)</td><td>21.18 (0.09)</td><td>21.31 (0.17)</td><td>20.65 (0.14)</td></tr><tr><td>B (0.6'')</td><td>23.67 (0.15)</td><td>23.59 (0.06)</td><td>23.59 (0.12)</td><td>23.67 (0.18)</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C (0.8'')</td><td>23.82 (0.15)</td><td>23.61 (0.06)</td><td>23.50 (0.12)</td><td>23.56 (0.18)</td><td>23.61 (0.62)</td><td>—</td><td>—</td></tr><tr><td>A + B (3.2'')</td><td>21.49 (0.16)</td><td>21.00 (0.03)</td><td>20.97 (0.05)</td><td>20.92 (0.06)</td><td>20.79 (0.09)</td><td>—</td><td>—</td></tr></tbody></table>

**Notes.** Data were taken on 2015 December 17 (MJD = 57374.28), at +539.2 d. All magnitudes are in the AB system and the errors include statistical and systematic errors. The aperture reported was measured from the *r'*-band image.

images (from February 2010 to 2013 May) were taken before the SN explosion. The position A was also marginally detected in Galaxy Evolution Explorer (GALEX) images and designated as GALEX J102241.6–165517, with a magnitude of NUV = 22.87 ± 0.39 (~2271 Å) in the AB system (GR6 catalogue<fnref n="5" />).

## 2.3. Host galaxy spectroscopy

The deep host spectra were taken with the 8.2 m Very Large Telescope (VLT) + FOcal Reducer and low dispersion Spectrograph (FORS2; Appenzeller et al. 1998) during 2015 January and February (approximately +294 d, listed in Table 1). The slit angle is 13.066 degrees rotated from the north (see Fig. 1 for the slit position) to cover the C (SN location) and A positions. We shifted the slit position by ~40 pixels for each individual frame, allowing us to remove the sky emission lines by two-dimensional (2D) image subtraction.

The 2D images were reduced using standard techniques within IRAF for bias subtraction and then removed cosmic rays using L.A.COSMIC (van Dokkum 2001), an algorithm for robust cosmic ray identification and rejection. We finally stacked those 15 frames using an inverse-variance weighted average. The stacked image (Fig. 1) clearly shows emission lines at three resolved spatial positions. The separations between those emission lines are consistent with the brightest regions at positions A, B, and C in the deep NTT + EFOSC2 *i*-band image.

We extracted 1D spectra from the stacked science frame using the IRAF/LONGSLIT package. We took a spectrophotometric standard LTT 2415 for the flux calibration and used daytime He+HgCd+Ar arcs for the wavelength calibration. More details are given in the Appendix. Subsequently, we corrected observed spectra for a Milky Way extinction of $A_V$ = 0.20 mag, which, for the assumed $R_V$ = 3.1, corresponds to $E(B-V)$ = 0.06 mag (Schlaflay & Finkbeiner 2011). We then shifted the spectra to the rest-frame: $z$ = 0.2556 for spectrum A, $z$ = 0.2553 for spectrum B and $z$ = 0.2563 for spectrum C. The redshift was estimated based on the observed wavelength of $H\beta$ and [O III] $\lambda\lambda$4959, 5007 lines. The slightly different redshifts of three components are clearly indicated from the emission line positions, see Fig. 1. Figure 2 shows the final spectra corresponding to positions A, B, and C.

We estimated radial velocity offsets from their relative redshifts: ~−90 km s⁻¹ from position A to B, ~300 km s⁻¹ from B to C, and ~210 km s⁻¹ from A to C. For comparison, the rotational velocity of the Milky Way is approximately 220 km s⁻¹ at the position of the Sun. The radial velocities are approximately 270 km s⁻¹ for the Large Magellanic Cloud (LMC) and 150 km s⁻¹ for the Small Magellanic Cloud (SMC; see Fig. 3 of Brüns et al. 2005). Hence, these three components may be kinematically associated, although it is not possible to

distinguish whether the three components are part of a single galaxy with B and C being tidal tails of A, or if C is a satellite galaxy of A. Therefore, we refer to them as a “host galaxy system” in this paper. For further discussion of the nature of this host system, we refer to Sect. 7, where we propose that the component B is a new-born H II region and the component C may be a satellite or merging with the component A.

## 3. Supernova results

### 3.1. Light curve

Figure 3 shows the observed light curves of LSQ14mo. We have ground-based photometry of LSQ14mo in *gri* bands covering the period −20 d to +70 d, and late-time (>250 d) detection limits. We also have a few epochs around maximum light from the *Swift* UV observations. The comparison of our data with the early-phase photometry of Leloudas et al. (2015a) shows excellent agreement. The two-day cadence survey strategy employed by LSQ provides a good constraint on the explosion epoch. The earliest detections of LSQ14mo appear to show an initial decline, which could be consistent with the double-peaked structure exhibited by LSQ14bdq (Nicholl et al. 2015b), SN 2006oz (Leloudas et al. 2012), DES14X3taz (Smith et al. 2016) and PTF12dam (Vreeswijk et al. 2017), and has been proposed to be common in SLSNe Ic (Nicholl & Smartt 2016). However, the photometric errors on these points are relatively large, and an early bump cannot be definitively identified in the case of LSQ14mo. The inset panel of Fig. 3 shows the absolute *g*-band pre-maximum light curves of LSQ14mo and LSQ14bdq with a stretch in time divided by 2.5.

Our best-sampled photometry is in the *r* band, which, at redshift $z$ = 0.256, approximately corresponds to the rest-frame *g* band. We converted to absolute magnitudes using a distance modulus from the cosmology calculator of Wright (2006), corrected for a Milky Way extinction ($A_V$ = 0.201), and applied *K*-corrections. We employed the PYSYNPHOT module within PYTHON to calculate synthetic photometric quantities for spectra at the given *gri*-band pass, then we shifted the spectra to the rest-frame ($z$ = 0.256) and calculated the synthetic photometry again. The differences between these were adopted as *K*-correction values, which were double-checked using the SNAKE code (Inserra et al. 2017). These values are given in Table C.4. We interpolated the corrections linearly to all epochs with imaging. Finally, we considered that the internal dust extinction from the host galaxy is negligible given by the ratio of $H\alpha$ and $H\beta$ lines.

In Fig. 4, we compare the rest-frame *g*-band light curve to a sample of other SLSNe. The peak luminosity, $M_g$ ≈ −21, is typical of these fast-declining events and fainter than those slowly-fading objects (for a larger sample comparison, we refer to Nicholl et al. 2015a). The rise time and overall

<footnote n="5">⁵ http://galex.stsci.edu/GR6/</footnote>
<figure><img src="image_3.png" /><figcaption>(a) Position A.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>(b) Position B.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>(c) Position C.</figcaption></figure>

**Fig. 2.** VLT + FORS2 spectra of the host galaxy system extracted at the three spatially resolved positions (A, B, and C; shown in red) and fit with the stellar population models from STARLIGHT (black lines overlapping observed spectra; see Sect. B). The blue spectra are subtracted residuals for the models. For position A, the *upper inset panels* show strong stellar absorption features of Balmer lines, while the *bottom inset panels* show detections of a weak auroral λ4363 line and clear detections of [N II] lines around Hα. For position C, the *inset panels* zoom in to WR feature regions, and no broad feature was detected. The atmospheric telluric absorption features are marked with the Earth symbol.

<figure><img src="image_6.png" /><figcaption>Rest-frame phase since r-band maximum (day)</figcaption></figure>

**Fig. 3.** Photometry of LSQ14mo from the UV to the optical. The early- and late-time $3\sigma$ detection limits are shown in empty symbols. The phase (day) has been corrected for time dilation ($z = 0.256$) and relative to the SN r-band maximum on MJD 56697. The period from 100 to 200 days is skipped due to no observations. The *inset panel* shows the absolute g-band pre-maximum light curves of LSQ14mo and LSQ14bdq with a stretch in time divided by 2.5, an early bump of LSQ14mo is plausibly detected. For comparison, the optical photometry given by Leloudas et al. (2015a), which is plotted as dashed lines, is in a good agreement with our measurements.

<figure><img src="image_7.png" /></figure>

**Fig. 4.** Absolute rest-frame g-band light curve of LSQ14mo compared with other type I SLSNe. The peak brightness of LSQ14mo of -21 mag is typical for fast-declining events and fainter than slowly-fading SLSNe. We note that after +50 d, LSQ14mo and SN 2011ke show a shallower decline rate, which is similar to a $^{56}$Co decay tail. Data for comparison are taken from Pastorello et al. (2010), Quimby et al. (2011), Leloudas et al. (2012), Inserra et al. (2013), Nicholl et al. (2013, 2015b, 2016b), Vreeswijk et al. (2014).

light curve width are very similar to SN 2011ke (Inserra et al. 2013), SN 2006oz (Leloudas et al. 2012), and SN 2010gx (Pastorello et al. 2010) and much narrower than other events, such as PTF12dam (Nicholl et al. 2013; Chen et al. 2015) and LSQ14bdq (Nicholl et al. 2015b). LSQ14mo appears to transition to a shallower decline rate after ≈30 d from peak. Similar transitions are apparent in other SLSNe, notably SN 2011ke, but the tail of LSQ14mo is brighter than SN 2011ke, and by +70 d it is much brighter than SN 2010gx, highlighting the rather diverse late-time behaviour of SLSN light curves.
## 3.2. Spectroscopic evolution

The NTT+EFOSC2 (of the PESSTO program) spectra of LSQ14mo are shown in Fig. 5. They have been corrected for Galactic extinction and shifted to the rest-frame at z = 0.256. The pre-maximum spectra are dominated by a blue continuum with a blackbody colour temperature of ≈13 000–15 000 K, and a distinctive W-shaped absorption feature between 4100 and 4500 Å. These lines have been attributed to O II by many authors since they were first identified by Quimby et al. (2011), and are generally seen in SLSNe Ic at this phase. The presence of ionised oxygen is consistent with the high temperatures around the peak. Recently, Mazzali et al. (2016) suggested that these lines must be excited non-thermally. One scenario that could explain this is X-ray injection from a magnetar.

As LSQ14mo ages, its spectrum resembles various other well-observed SLSNe, such as PTF09cnd (Quimby et al. 2011), SN 2013dg (Nicholl et al. 2014), and SN 2010gx (Pastorello et al. 2010). The continuum becomes substantially cooler between one and two weeks after maximum and the black-body fitting temperature drops from ~14 000 K at peak to ~7400 K at +15.2 d. At this phase, it is dominated by Fe II and intermediate-mass elements (bottom panel of Fig. 5). In Sect. 4, we examine the evolution of the lines and continuum in more detail using spectral synthesis models (Mazzali et al. 2016).

## 3.3. Bolometric light curve and model fits

To build the bolometric light curve of LSQ14mo, we used our photometry and spectroscopy to construct the spectral energy distribution (SED) as a function of time, and built a bolometric light curve with a code SUPERBOL<fnref n="6" />. As we lack coverage in the NIR and have Swift UV photometry only at maximum light, we estimated the UV and NIR contribution at each epoch with spectroscopy via blackbody fits. Since blackbody fits do not take into account metal absorption lines in the UV (e.g. Chomiuk et al. 2011), we used our Swift measurements to calculate the degree of absorption. Around -5d, we found that the flux in the UVW2 (filter central wavelength at 1928 Å) is 48% of the blackbody prediction, and 56% in the UVW1 (filter central wavelength at 2600 Å). We assumed that this suppression ratio is constant at all times.

We then linearly interpolated the K-corrected g- and i- band – and predicted UV and NIR – fluxes to all epochs with r-band photometry, corrected for the distance to the SN, and numerically integrated the SED assuming that the flux goes to zero outside of the Swift uvw2 and K bands. We also added additional 4% flux from the UV contribution for the last two epochs, because, as discussed by Chen et al. (2015), UV radiation is not zero in the late time. The fraction of energy emitted in the optical, UV, and NIR for LSQ14mo is very similar to other fast-declining SLSNe (Inserra et al. 2013). At -20 d, the fraction of UV is 41%, optical is 53%, and NIR is 6%; at +20 d, those fractions evolved to 12% UV, 69% optical, and 19% NIR; and at +70 d, these fractions are 4% UV, 60% optical, and 36% NIR. This gives us confidence in the method used to derive the fluxes from UVW2 through to K band, but does not account for flux beyond these limits. For comparison, we also present the optical pseudo-bolometric light curve as well, based on real ugri photometry.

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 5.** Top: spectral evolution of LSQ14mo taken by NTT+EFOSC2 of the PESSTO program. Blackbody curve fits (grey) show that the temperature before maximum light is ≈13 000–15 000 K and cooling rapidly to ~7000 K by 20 days after peak. The *inset* shows the Mg II absorption features from the host ISM taken at -2.2d. Bottom: comparison with other SLSN spectra. The prominent lines and their profiles are virtually identical to other SLSNe.</figcaption></figure>

The pseudo bolometric light curves and full bolometric light curves of LSQ14mo are plotted in Fig. 6 and the data points are listed for reference in Table C.5. Uncertainties are accounted for in the error bars in Table C.5. These include the error in photometric measurements and an uncertainty arising from the extrapolation methods for the missing filter data, using constant colour evolution or linear interpolation. The luminosity limits derived from our late non-detections are unfortunately not deep enough to provide constraints to the model, hence we do not plot these on the figures.

As shown by Fig. 6, the magnetar spin-down models (Inserra et al. 2013)<fnref n="7" /> provide an excellent fit to the data. We assumed that the magnetar spin-down energy is fully trapped in the ejecta and obtained a best fit (reduced $\chi^2$ of 1.40) with a magnetic field strength of $5.1 \times 10^{14}$ G and initial spin period of 3.9 ms with $M_{ej} = 3.9 M_{\odot}$ and $E_k = 3.14 \times 10^{51}$ erg. Figure 7 shows the evolution of temperature and velocity of LSQ14mo.

<footnote n="6">⁶ https://github.com/mnicholl/astro-scripts/blob/master/superbol.py</footnote>

<footnote n="7">⁷ Magnetar light curve fitting code is available at A. Jerkstrand's webpage: https://star.pst.qub.ac.uk/webdav/public/ajerkstrand/Codes/Genericarnett</footnote>
<figure><img src="image_10.png" /><figcaption>**Fig. 6.** Bolometric light curves of LSQ14mo and alternative model fits. Black empty circles show the pseudo bolometric light curve from *ugri* photometry and blue squares show the final bolometric light curve integrated from *UVW2* to *K*. The fully trapped magnetar model shown with a red curve indicates a magnetar with a magnetic field of $5.1 \times 10^{14}$ G and an initial spin period of 3.9 ms and an SN ejecta mass of $M_{ej} = 3.9 M_{\odot}$. The Ni-powered model requires $7 M_{\odot}$ of $^{56}$Ni, which is ~90% of the SN ejecta mass. The interacting model gives a SN ejecta of $M_{ej} = 3.3 M_{\odot}$ and a dense ($\rho_{CSM} = 7 \times 10^{-13} \text{ g cm}^{-3}$) CSM of $M_{CSM} = 2.0 M_{\odot}$.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 7.** Temperature and velocity evolution of LSQ14mo. Red and green lines are given from the magnetar modelling results. For observational data, the velocity was measured from the Fe II λ5169 line, and temperature was measured by fitting a black-body curve to spectra. The consistency of those two independent parameters is suggestive of a magnetar as the power source.</figcaption></figure>

The results from the magnetar model fit are compared with observational data: the velocity was measured from the minimum of Fe II λ5169 absorption and the temperature was measured by fitting a black-body curve to spectra. Those independent parameters from the model and from the observation got a good agreement, which supports a magnetar as the underlying power source of LSQ14mo.

Alternatively, we can fit the data with an interaction-powered model (Chatzopoulos et al. 2012). We assume a uniform CSM shell with an inner radius of $10^{15}$ cm (for details see Nicholl et al. 2014). Our best-fit model is shown in Fig. 6, and reproduces the data very well. The parameters for this fit are: $M_{ej} = 3.3 M_{\odot}$,

$M_{CSM} = 2.0 M_{\odot}$, $E_k = 6.6 \times 10^{50}$ erg, and $\rho_{CSM} = 7 \times 10^{-13} \text{ g cm}^{-3}$. Although this gives a satisfactory fit to the light curve, our spectroscopic modelling (Sect. 4) suggests that most of the luminosity is emitted by rapidly expanding material, therefore we disfavour this model relative to the magnetar scenario. However, as we discuss below, there may still be a (subdominant) contribution to the light curve from interaction.

Finally, we show a bolometric light curve model powered by the nuclear decay of $^{56}$Ni in Fig. 6. The $^{56}$Ni-powered LC model is derived in the same way as in Inserra et al. (2013). The model has a $^{56}$Ni mass of 7 $M_{\odot}$ with $M_{ej} = 8 M_{\odot}$ and $E_k = 5 \times 10^{52}$ erg. The unrealistically high fraction of $^{56}$Ni to ejecta disfavours the $^{56}$Ni-powered model. Figure 6 shows the magnetar model can mimic the $^{56}$Ni-powered light curve very well. As discussed in Moriya et al. (2017), they found that most $^{56}$Ni-powered light curves can be reproduced by magnetars that require the magnetar spin-down to be by almost pure dipole radiation with the breaking index close to 3.

### 3.4. Late-time spectrum of LSQ14mo compared to SN 1998bw

The late-time spectrum of LSQ14mo taken with VLT around +300d is dominated by the host galaxy. We estimated the detection limit of [O I] 6300 Å, the strongest feature of type Ic SNe in nebular phase, by comparing the spectrum of LSQ14mo to the spectrum of the bright broad-line type Ic SN 1998bw at +337d (Patat et al. 2001). We did this in two ways.

Firstly, we scaled the spectrum of SN 1998bw to match the photometry and placed it at the same distance as LSQ14mo. We compensated for the difference in phase according to the $^{56}$Co decay rate and γ-ray trapping of SN 1998bw. We then shifted the SN 1998bw spectrum to match the red part of the continuum in LSQ14mo, and estimated the detection limit of [O I] 6300 Å in the spectrum of LSQ14mo at +300d to be $1.4 \times 10^{-18} \text{ erg s}^{-1} \text{ cm}^{-2} \text{ Å}^{-1}$. This is similar to the comparison routines employed in Inserra et al. (2017).

Alternatively, we took the LSQ14mo FORS2 spectrum, subtracted the stellar population model at position C (see Fig. 2), and corrected the spectrum to rest-frame luminosity ($L_{\text{rest}}(\lambda) = 4\pi D_1^2(1+z)f_{\text{obs}}(\lambda)$). Then, we applied the same to the SN 1998bw spectrum, scaled it by five, and added it to the LSQ14mo spectrum. This illustrates that we would have detected the [O I] 6300 Å doublet if it had been five times more luminous than in SN 1998bw, or if the luminosity was $L_{6300} \ge 2.3 \times 10^{39} \text{ erg s}^{-1}$ (see Fig. 8).

The peak bolometric luminosity of LSQ14mo is $10^{44.04} \text{ erg s}^{-1}$, which is ten times brighter than that of SN 1998bw (~$10^{43.05} \text{ erg s}^{-1}$). Jerkstrand et al. (2017), Nicholl et al. (2016a) pointed out that the late-time spectrum of the SLSN SN 2015bn is very similar to SN 1998bw. Since our FORS2 spectrum shows that LSQ14mo is less than five times more luminous than SN 1998bw at late phases, this suggests an upper limit for the nickel mass of five times the nickel mass in SN 1998bw. If we take 0.4 $M_{\odot}$ as the representative $^{56}$Ni mass for SN 1998bw (Maeda et al. 2006), then the flux limit of LSQ14mo suggests an upper limit of approximately 2 $M_{\odot}$ of $^{56}$Ni ejected by LSQ14mo. This is much lower than the 7 $M_{\odot}$ $^{56}$Ni mass needed to power the luminosity of LSQ14mo at maximum light (see Sect. 3.3). Therefore, we rule out nickel-power as the underlying energy source in LSQ14mo. This spectroscopic late-time detection limit result supports the finding from Chen et al. (2013), who
<figure><img src="image_12.png" /></figure>

**Fig. 8.** VLT+FORS2 spectrum of LSQ14mo (position C) in black, with the stellar population model subtracted and corrected to restframe luminosity. The restframe spectrum of the broad-line type Ic SN 1998bw at +337d (Patat et al. 2001), scaled to the distance of LSQ14mo, is shown in red, indicating the strength of the [O I] λλ6300, 6364 doublet. The green spectrum is SN 1998bw multiplied by a factor 5. The blue spectrum shows LSQ14mo with the SN 1998bw ×5 spectrum added, illustrating an upper limit to the detection of the [O I] λλ6300, 6364 doublet. A broad feature of this strength would have been visible in the LSQ14mo spectrum. The atmospheric telluric absorption features are marked with the Earth symbol.

used a similar method. They found the upper limit of $ ^{56}\text{Ni} $ production for the fast-declining SLSN SN 2010gx to be similar to or below that of SN 1998bw (0.4 $M_{\odot}$). This small amount of $ ^{56}\text{Ni} $ cannot reproduce such bright peak luminosity. However, the late-time (+253 d) photometric detection limit of LSQ14mo ($9.31 \times 10^{37}$ erg s$^{-1}$) is not deep enough to compare with that of SN 1998bw ($1.35 \times 10^{37}$ erg s$^{-1}$).

# 4. Synthetic spectrum modelling of LSQ14mo

## 4.1. Theoretical modelling

The spectra of LSQ14mo were modelled following the approach described in Mazzali et al. (2016). These models use a Monte Carlo supernova radiative transfer code (Mazzali & Lucy 1993; Mazzali et al. 2000; Lucy 2000), and have been used for SNe Ia (e.g. Mazzali et al. 2014) as well as SNe Ib/c (Mazzali et al. 2013), and have recently been applied to SLSNe (Mazzali et al. 2016)$^8$.

The code assumes that the SN luminosity is emitted at a sharp inner boundary and follows the propagation of energy packets through the expanding SN envelope, which has a characteristic SN density distribution. In the case of SLSNe, Mazzali et al. (2016) find that a steep density gradient yields good fits to the absorption-line spectra of different SLSNe. They find that thermal conditions in SLSNe explain some of the lines, but O II lines near maximum and He I lines after maximum require non-thermal excitation/ionisation, which they claim may be the result of X-rays leaking from a magnetar-interaction.

## 4.2. Modelling of LSQ14mo

We adopt a model similar to those used in Mazzali et al. (2016). We prefer a modest mass for the ejecta, ~6 $M_{\odot}$. The kinetic energy is proportionally high, ~7 × 10$^{51}$ ergs, which gives a E/M (kinetic energy in units of 10$^{51}$ ergs and ejecta mass in log $M_{\odot}$) ratio of ~1. This mass is slightly higher than, but similar to, that of our magnetar light curve fit of 4 $M_{\odot}$ (Sect. 3.3), and the E/M ratio is also similar. Overall, the scenario we consider here is not inconsistent with the magnetar model inferred by the light curve fitting. The spectral synthesis models for the LSQ14mo spectral time series are presented in Fig. 9 and discussed individually below.

The first spectrum at −7.0 d has formally an epoch of 18 days after explosion (fits obtained from the model). This is well within the reconstruction of the detection limit and light curve of LSQ14mo. The −7.0 d spectrum is characterised by a steep blue continuum and absorption lines of mostly O II. This spectrum has a luminosity log L = 44.10 (erg s$^{-1}$) and a photospheric velocity of 17 000 km s$^{-1}$. The composition is (relatively) dominated by O (0.51) and C (0.31). We assume He (0.10) and lower abundances of Mg (0.01), Si (0.04), Ne (0.02) and other intermediate-mass elements, including S, Ca, Ti and Cr. The Fe abundance in the layers above the photosphere is approximately 1/3 solar, consistent with the host metallicity measurement. The strongest features reproduced are the O II near 3200, 3600, and 3900 Å and the distinctive 4200 Å (blended with Si III) and 4400 Å blends of O II transitions. The spectrum does appear to turn over at 3200, which is probably the Si III line near 3000 Å. Other very strong lines are predicted in the UV (further blue than our rest-frame spectra) including Mg II near 2600 Å. The C II line at 6300 Å is predicted to be shallow and weak, and is not visible at the signal-to-noise of our spectra in this region.

The second spectrum at −2.2 d has log L = 44.04 (erg s$^{-1}$) and a photospheric velocity of 15 000 km s$^{-1}$. This is a small change in L but a 2000 km s$^{-1}$ change in velocity over five days. The observed spectrum is very similar to that of −7.0 d. We used a very similar composition as for −7.0 d, with slightly less Mg. The abundances of several elements would be best set if the flux further in the UV was known, which is unfortunately not the case. Our model reproduces the overall spectral shape and the quality and strength of most lines. There may be reprocessed flux from the UV via metal lines in the observed spectrum, but it is unlikely that line fluorescence processes lead to such a smooth enhancement of the flux in such a limited spectral region. We may instead be witnessing the beginning of a short-lived additional energy source, as we discuss in the following section.

The third spectrum, taken at +8.8 d, is less blue and the O II lines have disappeared. This intermediate spectrum has a luminosity log L = 43.83 (erg s$^{-1}$) and a photospheric velocity of 13 750 km s$^{-1}$. Other lines in the model are similar to the later epoch, but Fe II and Mg II are weaker. Fe III lines are still present at 4700–5000 Å. The line at 4300 Å could be matched with He I. Near 3700 Å, Ca II is also very weak, and the line may be fit with He I. However, the fit to this spectrum is relatively poor, and some of the strongest predicted helium lines are not observed (e.g. 5876 and 6678 Å). Therefore, the assumption of 10% He may be too high, suggesting a highly-stripped core. It appears that while the luminosity is still high and the continuum blue ($T_{\text{eff}}$ ≈ 10 000 K), most lines are compatible with a cool spectrum such as that observed at +21.7 d. The spectrum of LSQ14mo changes quite dramatically between the previous epoch (−2.2 d) and this one at +8.8 d.

<footnote>$^8$ The luminosity of PTF09cnd is one order of magnitude smaller than what indicated in their Table 3. This error only affects the table, not the actual calculation.</footnote>
<figure><img src="image_13.png" /><figcaption>**Fig. 9.** Spectral synthesis model (in black) from Mazzali et al. (2016) applied for LSQ14mo spectra. Each spectrum is labelled with its rest-frame time from the *r*-band light curve maximum. The −7.0 d spectrum is dominated by O and C, other intermediate-mass elements also shown and labelled. The −2.2 d spectrum is similar to −7.0 d in terms of element identification, temperature, luminosity and velocity. The +8.8 d spectrum, with a poor model fitting, shows a hot continuum but rather cool lines. This implies an additional luminosity may come from the SN ejecta interacting with a H- and He-poor shell. The +15.2 d spectrum is similar to +21.7 d and the strongest features are intermediate elements such as Si II, Fe II, Mg II and Ca II. All identified lines are labelled.</figcaption></figure>

The fourth spectrum at +15.2 d is very similar to that at +21.7 d, but it also has similar lines as that at +8.8 d. Our model has log L = 43.68 (erg s⁻¹) and a photospheric velocity of 13 250 km s⁻¹. This means that a significant drop in luminosity is not accompanied by a drop in velocity. The model spectrum is redder than that at +8.8 d, but it still does not reproduce the observed lines very well because the temperature is too high. The composition is the same as at +21.7 d. Excess flux seems to be present from the NUV all the way to ~6000 Å.

The fifth spectrum at +21.7 d is relatively red, and it shows similar lines to the ones in the +8.8 d spectrum. It has a luminosity log L = 43.53 (erg s⁻¹) and a photospheric velocity of 12 500 km s⁻¹. The composition is similar to the earlier epoch. The strongest features in the model are Si II near 6100 Å (a noisy region in the observed spectrum, but not clearly observed), Fe II multiplet 48 lines near 4700–4800 Å, Fe II and Mg II near 4300 Å, Ca II blended with Si II near 3700 Å, and Ti II and Cr II

near 3100 Å. Heavy line blocking suppresses the NUV flux. Compared to iPTF13ajg (Mazzali et al. 2016, Fig. 4), the luminosity is a lot lower, singly ionised species are now present, and many of the features that were attributed to He I for iPTF13ajg spectra are fit by Fe II or Si II in LSQ14mo.

It should be noted that LSQ14mo is significantly less luminous than iPTF13ajg, and that it evolves much more rapidly than both iPTF13ajg and PTF09cnd (see Fig. 4), and therefore, direct comparisons must be treated with caution. Near maximum, the O II-rich spectrum of LSQ14mo may be explained by non-thermal excitation. The later absence of He lines may be interpreted either as an indication that the He abundance is small, or that non-thermal processes are no longer at play, which is less likely. The modelling of iPTF13ajg (Mazzali et al. 2016) indicated the presence of He I lines, but these were not required in PTF09cnd. Overall, we do not see an obvious requirement for helium in LSQ14mo.
<figure><img src="image_14.png" /><figcaption>**Fig. 10.** Subtracted continuum (assumed black-body fitting) spectra of LSQ14mo at explosion epochs +8.8 d, +15.2 d, and +21.7 d. The identification of spectral element features are the same at those three epochs. This supports the idea that the blue continuum observed in the +8.8 d spectrum could be due to an additional energy source that likely comes from the interaction between the SN ejecta and a thin CO shell.</figcaption></figure>

### 4.3. Interaction contribution seen in spectra

The spectrum after maximum (+8.8 d) shows a hot continuum, but rather cool lines. This is difficult to reproduce if the SN ejecta are in radiative equilibrium, and is actually the opposite situation with respect to the epochs when O II lines are present (in that case, the excitation temperature of the lines is above the continuum temperature). One option is that the SN absorption-line spectrum is indeed a cool one, and that additional luminosity is provided in a form that does not affect the lines. In Fig. 10, we present the spectra of LSQ14mo at +8.8 d, +15.2 d, and +21.7 d with the continuum subtracted (assuming a black-body). The subtracted spectra have a similar profile; in fact, the line identifications are the same at those three epochs. This could arise if the ejecta were interacting with some porous or clumpy, H- and He-poor shell, which might add a hot black-body component to the flux; for example, SN 2009dc had its luminosity augmented by interaction of the ejecta with a H-/He-poor circumstellar medium (Hachinger et al. 2012), as opposed to a H-rich one (e.g. Deng et al. 2004).

We suggest that we may be witnessing the contribution of a partially thermalised, additional source of luminosity, starting after −2.2 d and ending before +21.7 d. The spectrum at +8.8 d is the most heavily affected. Because the nature of the lines changes between −2.2 d and +8.8 d (from O II to Fe II and Mg II), while the continuum remains blue, it appears that this additional component may not come from the same physical location where the lines are formed.

## 5. Host galaxy properties

### 5.1. Galaxy size, extinction corrections, and luminosity

The profile of the host galaxy is noticeably broader than the stellar PSF, hence we can determine the physical diameter of the extended source, assuming the relation (galaxy observed $FWHM)^2 = (\text{PSF FWHM})^2 + (\text{intrinsic galaxy FWHM})^2$. We measured the FWHM at both the host positions A (1.74'') and C (1.03'') and the average FWHM (0.80'') of seven reference stars within a 2' radius around the host on the NTT i-band image

taken on 2015 February 12, which had the best seeing conditions (0.8''). We were unable to fit the FWHM of a Gaussian profile at position B because of its irregular profile. This provides a physical diameter of 6.17 kpc for position A at the angular size distance of 824.4 Mpc, and a physical diameter of 2.63 kpc for position C at 826.0 Mpc. For comparison, the size of the LMC is ~4.3 kpc in diameter, and 2.1 kpc for the SMC.

The angular separation between positions A and C is ~4'', measured from the brightest region of each galaxy component, which corresponds to ~15.2 kpc at the distance of the LSQ14mo host system. The separation between positions A and B is ~10.5 kpc, and between positions B and C it is ~4.7 kpc. Those separations were consistent with that measured from the emission line positions in three spectral traces.

We used the Balmer decrements to estimate the internal dust extinction of the host galaxy system. The intrinsic line ratio is expected to be $H\alpha/H\beta = 2.86$, assuming case B recombination for $T_e = 10\,000$ K and $n_e = 100\,\text{cm}^{-3}$ (Osterbrock 1989). The measured fluxes (Table C.3) of $H\alpha$ and $H\beta$ give a ratio of 3.10 at position A, which indicates an internal dust extinction of approximately $A_V = 0.26$ in the rest-frame, assuming $R_V = 3.1$. For positions B and C, the $H\alpha/H\beta$ value is slightly lower than 2.86, corresponding to 2.57 and 2.74, respectively. Therefore, we assumed that the internal dust reddening is negligible in positions B and C, that is, $A_V = 0$.

Finally, we obtained the absolute *g*-band magnitudes ($M_g$) from the observed *r*-band magnitude ($m_r$) after applying *K*-correction (Table C.4), foreground and internal dust extincions, using the formula: $M_g = m_r - A_r$, MilkyWay + $K_{r\to g} - A_{g, \text{host}}$. This gives us $M_g = -19.04 \pm 0.15$ mag, $-16.94 \pm 0.15$ mag, and $-16.85 \pm 0.15$ mag at the three components A, B, and C, respectively.

### 5.2. Stellar population ages in the host galaxy system

There are multiple stellar populations inside the galaxy. We used the EW of $H\alpha$ line as a tracer for ongoing star formation and as an age indicator to probe a young stellar population of H II region. Such methods have also been employed recently to analyse the spectra of the environments of local SNe (e.g. Leloudas et al. 2011; Kuncarayakti et al. 2013). We found the rest-frame EWs of $H\alpha$ to be 54 Å, 175 Å, and 85 Å at positions A, B, and C, respectively. To be more consistent with our metallicity measurement of the host system of LSQ14mo, we chose the new simple stellar population (SSP) model with low-metallicity (0.2 $Z_\odot$) from Kuncarayakti et al. (2016), which gave slightly older ages than that given by solar metallicity. We further considered the SSP with a Kroupa IMF (Levesque & Leitherer 2013), which is similar to our assumption of a Chabrier IMF. Our EW range of the $H\alpha$ line is much smaller than the continuous star-formation models (>500 Å) and therefore we chose the SSP model of an instantaneous burst of star formation. This resulted in a stellar population age between 6.3 and 9 Myr for position A (as its metallicity of 0.4-0.8 $Z_\odot$ is between the SSP models of 0.2-1 $Z_\odot$), 6.3 Myr at position B and 7.7 Myr at position C. The youngest stellar population is located at the B position, while older, relatively similar stellar populations are present at positions A and C (SN location). The largest EW of $H\alpha$ line at position B supports our interpretation that the position B is a new-born bright H II region from the interacting activities of galaxies A and C.

Galaxies certainly contain multiple stellar populations rather than a single one. We also reported the *r*-band light-weighted
<table><caption>**Table 3.** Temperature, photospheric velocity and luminosity evolution of LSQ14mo from the spectral modelling results.</caption><thead><tr><th>Date</th><th>Phase (day)</th><th>Exp. epoch (day)</th><th>Temp<sub>R</sub> (K)</th><th>Temp<sub>*</sub> (K)</th><th>log L (erg s<sup>-1</sup>)</th><th>Velocity (km s<sup>-1</sup>)</th><th>Temp<sub>BB</sub> (K)</th><th>v<sub>5169</sub> (km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>2014 Jan. 31</td><td>-7.0</td><td>18</td><td>14 885</td><td>13 272</td><td>44.10</td><td>17 000</td><td>14 641</td><td>-</td></tr><tr><td>2014 Feb. 6</td><td>-2.2</td><td>23</td><td>13 971</td><td>12 351</td><td>44.04</td><td>15 000</td><td>13 636</td><td>-</td></tr><tr><td>2014 Feb. 20</td><td>+8.8</td><td>34</td><td>9674</td><td>8961</td><td>43.83</td><td>13 750</td><td>10 083</td><td>10 500 ± 1500</td></tr><tr><td>2014 Feb. 28</td><td>+15.2</td><td>40</td><td>7698</td><td>7465</td><td>43.68</td><td>13 250</td><td>7421</td><td>10 400 ± 500</td></tr><tr><td>2014 Mar. 8</td><td>+21.7</td><td>47</td><td>6965</td><td>6700</td><td>43.53</td><td>12 500</td><td>6488</td><td>10 200 ± 1500</td></tr></tbody></table>

**Notes.** The phase (day) has been corrected for time dilation ($z = 0.256$) and relative to the SN r-band maximum on MJD 56 697. The explosion epoch is given from the spectral model. Temp<sub>R</sub> refers to radiation temperature and Temp<sub>*</sub> is effective temperature. Additionally, we list a black body temperature (Temp<sub>BB</sub>) estimated from fitting the spectrum assuming a black body; this value is similar with the radiation temperature. $v_{5169}$ measured from minimum of Fe II λ5169 absorption is also shown.

stellar population age given by MAGPHYS. The mean age of position A is slightly older at 212 Myr, than the 94 Myr and 118 Myr at positions B and C, respectively. Alternatively, the fitting from spectral continuum using STARLIGHT also implies both young (3-7 Myr) and old (500 Myr) populations at positions A and C. However, those values are model dependent and there are not only those populations inside the galaxy. We can only argue that a young stellar population is present as Hα emission line is detected. For comparison, to date, Thöne et al. (2015) reported the largest EW of Hα line (>800 Å) among the SLSN host galaxies that has been seen in the host of PTF12dam, and thus suggested a very young stellar population at the SN site of ~3 Myr.

## 5.3. Stellar masses and star-formation rates in the host galaxy system

To obtain stellar masses of the LSQ14mo host galaxy system, we fit the available photometry with the stellar population models in MAGPHYS (da Cunha et al. 2008) using Bruzual & Charlot (2003) templates. The best fit is shown in Fig. 11, and returns masses (in unit of log $M/M_⊙$) of $8.8_{-0.2}^{+0.1}$ for position A, $7.6_{-0.2}^{+0.3}$ for position B, and $7.7_{-0.2}^{+0.2}$ for position C. Here, the stellar mass is the median of the probability density function (PDF) over a range of models and the errors correspond to the $1\sigma$ credible intervals of the PDF. For comparison, the stellar mass of the SN position (C) given by Schulze et al. (2016) of $7.89_{-0.19}^{+0.15}$ log $M/M_⊙$ from the SED fitting is shown, which is consistent with our estimate.

After correcting the Hα line luminosity for Milky Way and internal host dust reddening, we applied the conversion of SFR = $L(Hα) \times 7.9 \times 10^{-42}$ to calculate the star-formation rate (SFR; Kennicutt 1998), and then divided it by 1.6 assuming a Chabrier IMF. The SFR is 0.52 $M_⊙$ yr$^{-1}$ for component A, 0.08 $M_⊙$ yr$^{-1}$ for component B and 0.06 $M_⊙$ yr$^{-1}$ for component C. Together with their stellar masses, we determined a specific SFR (sSFR) of 0.82 Gyr$^{-1}$, 1.89 Gyr$^{-1}$, and 1.16 Gyr$^{-1}$ for components A, B, and C, respectively.

## 5.4. Gas-phase metallicity of the host galaxy

We carry out an abundance analysis based on both the strong line methods and the “direct” measurement, derived from the weak [O III] λ4363 auroral line detected at position A. The emission line fluxes (Table C.3) have been corrected for Milky Way and internal dust extinction (Sect. 5.1).

First of all, we checked the locations of the host system of LSQ14mo on the BPT diagram (Baldwin et al. 1981). The $log([O III]/Hβ) = 0.39, 0.58, and 0.52$ and

<figure><img src="image_15.png" /><figcaption>**Fig. 11.** Best-fit SED models of the host system of LSQ14mo from MAGPHYS, for three components, A, B, and C, are shown in blue, green, and red lines. The colour circles show the different bands’ photometry (shifted to the rest-frame) for the three components of host system where it was possible to have a clear detection. For instance, for component A, the colour circles from left to right are from NUV, griz to JHK-bands photometry. The median and a $1\sigma$ range of stellar masses of components A, B, and C are: $8.8_{-0.2}^{+0.1}$, $7.6_{-0.2}^{+0.3}$ and $7.7_{-0.2}^{+0.2}$ log $M/M_⊙$, respectively.</figcaption></figure>

$log([N II]/Hα) = -0.91, -1.47, and -1.24$ for positions A, B, and C, separately. They are all located at the star-forming galaxy region, while components B and C are near the high ionised boundary with $log([O III]/Hβ) > 0.5$. This characteristic of SLSN host galaxies has been noticed in Leloudas et al. (2015b), and the host galaxy of LSQ14mo (position C) also follows this trend.

We used the open-source PYTHON code PYMCZ (Bianco et al. 2016) to calculate the oxygen abundance in different diagnostics. We adopted the Kobulnicky & Kewley (2004, hereafter KK04) calibration of the $R_{23}$ method, which uses the $([O III]λ4959, 5007 + [O II]λ3727)/Hβ$ ratio. We found an oxygen abundance of $12 + \log(O/H) = 8.62 \pm 0.03$, $8.28 \pm 0.03$, and $8.18 \pm 0.02$ for positions A, B, and C, respectively (where the Hα and [N II] λ6584 ratio was used to break the $R_{23}$ degeneracy). We also used the N2 method of Pettini & Pagel (2004, hereafter PP04), which uses the $log([N II]λ6583/Hα)$ ratio, giving $12 + \log(O/H) = 8.31 \pm 0.01$, $8.09 \pm 0.05$, and $8.18 \pm 0.05$ for positions A, B, and C, respectively. In addition, we considered a new metallicity diagnostic proposed by
<table><caption>**Table 4.** Main properties of the host galaxy system of LSQ14mo.</caption><thead><tr><th></th><th>Position A</th><th>Position B</th><th>Position C</th></tr></thead><tbody><tr><td>RA (J2000)</td><td>10:22:41.55</td><td>10:22:41.52</td><td>10:22:41.48</td></tr><tr><td>Dec (J2000)</td><td>-16:55:17.9</td><td>-16:55:15.2</td><td>-16:55:14.2</td></tr><tr><td>Redshift</td><td>0.2556</td><td>0.2553</td><td>0.2563</td></tr><tr><td>Projected separation (kpc)</td><td>–</td><td>10.5 (A to B)</td><td>15.2 (A to C)</td></tr><tr><td>Physical diameter (kpc)</td><td>6.17</td><td>–</td><td>2.63</td></tr><tr><td>Radial velocity offset (km s<sup>-1</sup>)</td><td>–</td><td>-90 (A to B)</td><td>210 (A to C)</td></tr><tr><td>Apparent r (mag)</td><td>21.26 ± 0.02</td><td>23.58 ± 0.06</td><td>23.60 ± 0.06</td></tr><tr><td>Galactic extinction A<sub>V</sub> (mag)</td><td>0.20</td><td>0.20</td><td>0.20</td></tr><tr><td>Internal extinction A<sub>V</sub> (mag)</td><td>~0.26</td><td>~0.0</td><td>~0.0</td></tr><tr><td>Absolute g (mag)</td><td>-19.04 ± 0.15</td><td>-16.94 ± 0.15</td><td>-16.85 ± 0.15</td></tr><tr><td>GALEX NUV (mag)</td><td>22.87 ± 0.39</td><td>–</td><td>–</td></tr><tr><td>Hα luminosity (erg s<sup>-1</sup>)</td><td>1.04 × 10<sup>41</sup></td><td>1.62 × 10<sup>40</sup></td><td>1.30 × 10<sup>40</sup></td></tr><tr><td>Hα EW<sub>rest</sub> (Å)</td><td>54.10</td><td>174.67</td><td>84.94</td></tr><tr><td>SFR from Hα (M<sub>☉</sub> yr<sup>-1</sup>)</td><td>0.52</td><td>0.08</td><td>0.06</td></tr><tr><td>log stellar mass (M<sub>☉</sub>)</td><td>8.8<sup>+0.1</sup><sub>-0.2</sub></td><td>7.6<sup>+0.3</sup><sub>-0.2</sub></td><td>7.7<sup>+0.2</sup><sub>-0.2</sub></td></tr><tr><td>sSFR (Gyr<sup>-1</sup>)</td><td>0.82</td><td>1.89</td><td>1.16</td></tr><tr><td>12 + log (O/H) (T<sub>e</sub>)</td><td>8.00 ± 0.31</td><td>–</td><td>–</td></tr><tr><td>12 + log (O/H) (KK04 R<sub>23</sub>)</td><td>8.62 ± 0.03</td><td>8.28 ± 0.03</td><td>8.18 ± 0.02</td></tr><tr><td>12 + log (O/H) (PP04 N2)</td><td>8.31 ± 0.01</td><td>8.09 ± 0.05</td><td>8.18 ± 0.05</td></tr><tr><td>12 + log (O/H) (D16)</td><td>8.02 ± 0.03</td><td>–</td><td>7.92 ± 0.16</td></tr><tr><td>Youngest stellar population from Hα (Myr)</td><td>6.3–9</td><td>6.3</td><td>7.7</td></tr></tbody></table>

**Notes.** The SN site is located at the position C.

Dopita et al. (2016, hereafter D16), which uses [N II] λ6583, Hα, and [S II] λ66717, 6731 lines. We found 12 + log(O/H) = 8.02 ± 0.03 for position A and 7.92 ± 0.16 for position C. We note that, for position B, the line ratio combination required for the D16 diagnostic returns a value ($y = -1.12$) at the edge of where their calibration is applicable. Therefore, we did not apply this metallicity diagnostic at position B.

We also attempted an estimate of the metallicity at position A from direct measurements of the [O III] electron temperature. The auroral [O III] λ4363 line is intrinsically weak and difficult to detect since all SLSNe have been found at distances greater than $z = 0.1$ (see Nicholl et al. 2015a). The marginal detection ($S/N = 3.6$) of the auroral line seen in the spectrum of position A gives an electron temperature of $T_e = 13\,000$ K assuming an electron density of $100\,\text{cm}^{-3}$. This indicates a metallicity of $12+\log(\text{O/H}) = 8.00 \pm 0.31$.

The metallicities obtained from each of these diagnostics are reported in Table 4. To conclude, the SN location (position C) exhibits the lowest metallicity when using the KK04 and D16 diagnostics, even when accounting for the ~0.1 dex systematic uncertainty in these diagnostics. This is perhaps unsurprising, given the know correlation between $M_*$ and metallicity (e.g. Tremonti et al. 2004), as position C also has the lowest stellar mass. It is known that there is a systematic offset of 0.2–0.4 dex (and sometimes up to 0.6 dex) between the KK04 method and both the PP04 and direct $T_e$ method (Liang et al. 2007; Bresolin 2011), with the former giving systematically higher values. As discussed in Kewley & Ellison (2008), for example, position A with a stellar mass of $8.80 \log M/M_☉$, a typical offset between KK04 R$_{23}$ and PP04 N2 is 0.35 dex, which is consistent with the values we measured.

We further obtained a nitrogen-to-oxygen ratio of $\log(\text{N/O}) = -1.35 \pm 0.24$ at position A from the R$_{23}$ calibration method, assuming that $\text{N/O} = \text{N}^+/\text{O}^+$. We measured $\log(\text{N/O}) = -1.58 \pm 0.86$ at position B, and $-1.36 \pm 1.60$ at position C. This places those three components close

to the “plateau” seen at $\log(\text{N/O}) \approx -1.5$ on the N/O versus O/H diagram of galaxies (Pilyugin et al. 2003). For comparison, the mean values of oxygen and nitrogen abundances in 21 H II regions of the Small Magellanic Cloud are $12 + \log(\text{O/H}) = 8.07 \pm 0.07$ ($T_e$ method) and $\log(\text{N/O}) = -1.55 \pm 0.08$ (from literature values summarised in Pilyugin et al. 2003). Hence, the metallicity environments in the host system of LSQ14mo are similar to (or perhaps somewhat more extreme than) that of the SMC.

## 6. Discussion I: superluminous supernovae with interacting features

As the hot continuum has been seen in our $+8.8\,\text{d}$ spectrum, we propose that this additional energy is from the interaction between the SN ejecta and a thin shell or some clumps inside the CSM. This material likely originates from a mass loss before the SN explosion. Here we roughly calculated an upper limit on the energy input from such as interaction, and the place of this interacting shell.

Firstly, we assumed $+21.7\,\text{d}$ spectrum as a baseline, and that any additional luminosities at $+8.8\,\text{d}$ and $+15.2\,\text{d}$ spectra came from the interaction. This is of course an upper limit on the interaction energy, since the SN ejecta also cools down from $+8.8\,\text{d}$ to $+21.7\,\text{d}$. The upper limit on the luminosity from the interaction is up to 60% at $+8.8\,\text{d}$ and 36% at $+15.2\,\text{d}$. From our synthetic spectrum modelling results, the two early-phase spectra ($-7.0\,\text{d}$ and $-2.2\,\text{d}$) were well fitted by the model (see Fig. 9). Hence we propose that the interaction only happened after $-2.2\,\text{d}$ and that its energy contribution peaked at $+8.8\,\text{d}$. We assumed that the interaction plays a role of 60% at $+8.8\,\text{d}$ and 36% at $+15.2\,\text{d}$, and keep 30% after $+15.2\,\text{d}$ while 0% around and before the SN peak. This is clearly a coarse assumption, but we use this simple approach to illustrate the essential points. Figure 12 demonstrates the assumed bolometric light curve after subtracting the interaction contribution. We assumed the remaining luminosity
<figure><img src="image_16.png" /></figure>

**Fig. 12.** Combination of magnetar-powered and interaction contribution of bolometric luminosity of LSQ14mo. The symbol “S” marks the time when spectra were taken, and the green “S” highlights the epoch with a hot continuum seen in the spectrum. Blue squares show the final bolometric light curve integrated from UVW2 to K and fitted with a fully trapped magnetar model in red; black empty circles show the lower limit of assumed magnetar-powered light curve and fitted in blue. The luminosity difference between red and blue magnetar models is an additional energy source in green, which is speculated to come from the interaction between the SN ejecta and a thin shell or some clumps inside the CSM likely originating from a massive star mass loss before the SN explosion.

is powered purely by a magnetar. Fitting this light curve (reduced $\chi^2 = 9.53$), we obtain a magnetic field strength of $5.5 \times 10^{14}$ G and an initial spin period of 4.8 ms with $M_{ej} = 2.7 M_{\odot}$. We subtract this new magnetar model fit from the model for the full bolometric light curve (Sect. 3.3) and plot the difference in Fig. 12 to give an estimate of the light curve component from interaction.

Secondly, we estimated the location of this interacting material. We determined an initial radius ($R_0$) of the SN photosphere at $-7.0$ d from the equation of $L = 4\pi R_0^2 \sigma T^4$, where $\sigma$ is the Stefan-Boltzmann constant, and we used the effective temperature and luminosity listed in Table 3. Thus we got the $R_0 = 2.39 \times 10^{15}$ cm. We approximated the radius of the next epoch to $R = R_0 + v \times \delta_t$, where $v$ is the photospheric velocity, which we adopted to be the average between two epochs from $R_0$ to $R$. Then we moved on to the next epoch of +8.8 d, and obtained the radius of $4.51 \times 10^{15}$ cm of this interacting thin shell.

LSQ14mo is the first event for which interacting features in fast-evolving SLSN spectra have been reported, by seeing a hot-ter continuum and cooler lines. We therefore searched for other SLSNe with similar behaviour and found several fast-declining SLSNe at a similar phase to LSQ14mo, with interaction around 10 days after the maximum light. They are: SN 2013dg at +9 d, SN 2011ke at +9 d, and LSQ12dlf at +21 d, shown in Fig. 13. A further larger sample study will be interesting to investigate.

Overall, LSQ14mo is one of the fastest evolving SLSNe of type I which has a well constrained rise and decline (see Fig. 4). This fact leads to a relatively low ejecta mass of approximately $4 M_{\odot}$ from the magnetar powering fits. This mass is on the lower end of the distribution of low redshift SLSNe Ic presented by Nicholl et al. (2015a). For such low masses, magnetar powering inherently implies formation of a neutron star remnant and hence

<figure><img src="image_17.png" /></figure>

**Fig. 13.** Spectral comparison between LSQ14mo and other fast-declining SLSNe at a similar observational phase. The interacting feature of a blue continuum but rather cool lines is also shown in SN 2013dg at +9d, SN 2011ke at +9d, and LSQ12dlf at +21d. Data for comparison are taken from Inserra et al. (2013), Nicholl et al. (2014).

a progenitor mass of $\sim 6 M_{\odot}$. Such a carbon-oxygen star would be relatively low in mass for the final mass of a Wolf-Rayet star (Crowther 2007).

# 7. Discussion II: interacting host galaxy

## 7.1. Star-formation activity

In Fig. 14, we plot the relation between $M_*$ and specific SFR ($sSFR = SFR/M_*$) for systems A, B, and C, along with a selection of local galaxies. SFRs for the Lee et al. (2004), Kennicutt et al. (2008), Izotov et al. (2011), and Thuan et al. (2016) samples are measured via the dust-corrected $H\alpha$ flux or narrow-band luminosity, using the $L(H\alpha)$-SFR relation from Kennicutt (1998) and the Cardelli et al. (1989) extinction law. Total SFRs for the Guseva et al. (2009), Yates et al. (2012), and Berg et al. (2012) samples were obtained from the SDSS-DR7 spectroscopic catalogue (Brinchmann et al. 2004). All SFRs and stellar masses have been corrected to a Chabrier IMF.

We note that the Kennicutt et al. (2008) and Berg et al. (2012) data (orange and green points) contain very local galaxies ($z < 0.002$), the Izotov et al. (2011) data (dark green points) comprise “green pea” galaxies with particularly high specific SFRs, and the other data represent a statistically significant sample of $H\alpha$-bright dwarf galaxies out to a redshift of $z = 0.3$. Therefore, these samples combined constitute a fair representation of the range of sSFR found in the local dwarf galaxy population.

We can see from Fig. 14 that the three LSQ14mo systems each have typical sSFR compared to the general star-forming dwarf-galaxy population out to $z = 0.3$. This is consistent with the idea that SLSN hosts, such as system C, need to have a high sSFR in an absolute sense (i.e. $\log(sSFR/yr^{-1}) \ge -10.3$) (Leloudas et al. 2015b; Chen et al. 2016), but not necessarily higher than other galaxies of the same stellar mass.

## 7.2. Interacting galaxy system

All three systems, A, B, and C, are in close proximity to each other. Their projected separations are $r_{proj,AB} \sim 10.5$ kpc,
<figure><img src="image_18.png" /><figcaption>**Fig. 14.** Relation between stellar mass ($M_*$) and specific SFR (sSFR) for positions A, B, and C (red points) from LSQ14mo system compared to a selection of local dwarf galaxies ($z < 0.3$). The contours represent the star-forming sample of Yates et al. (2012), showing the 68th, 95th, and 99.5th percentiles of the distribution. All three LSQ14mo systems have SFRs (and sSFRs) slightly higher than most galaxies of their mass.</figcaption></figure>

$r_{\text{proj,AC}} \sim 15.2$ kpc, and $r_{\text{proj,BC}} \sim 4.7$ kpc, and their line-of-sight velocity ($v_{\text{los}} \equiv c z$) differences are $\Delta v_{\text{los,AB}} = 89.94$ km s$^{-1}$, $\Delta v_{\text{los,AC}} = -209.86$ km s$^{-1}$, and $\Delta v_{\text{los,BC}} = -299.80$ km s$^{-1}$, respectively. This means that all three systems are easily close enough to be considered interacting under typical selection criteria (e.g. Kewley et al. 2006; Ellison et al. 2008; Michel-Dansac et al. 2008).

Interacting systems are known to exhibit enhanced SFRs and lowered gas-phase metallicities compared to isolated galaxies of the same mass (e.g. Lee et al. 2004; Peeples et al. 2009). This is also seen for other SLSN type I host galaxies with disturbed morphologies (Lunnan et al. 2015; Perley et al. 2016). However, none of the three systems in LSQ14mo have particularly high sSFR or low $Z_g$ for their stellar mass. Integral field unit (IFU) observations are therefore required to investigate this host system more carefully, in order to better determine its interacting state.

The maximum escape velocity of system C from A (i.e. assuming that their true separation is simply their projected separation, $r_{\text{proj,AC}}$) can also be estimated. Firstly, assuming system A is a central galaxy, we can use the multi-epoch abundance matching technique of Moster et al. (2013) to obtain an estimate of its dark matter (DM) halo mass. Abundance matching statistically determines the relationship between stellar mass and DM halo mass by fitting the observed galaxy stellar mass function to simulations of dark matter halo growth (e.g. Moster et al. 2010). For the stellar mass and redshift of system A, this technique estimates a DM halo mass of $\log(M_{\text{DM}}/M_{\odot}) \sim 11.0$. In order to estimate the amount of this DM that is within $r_{\text{proj,AC}}$, we further assume the standard NFW (Navarro et al. 1997) DM density profile and a DM concentration parameter of 16, both of which are suitable for DM haloes of this size (Sawala et al. 2015; Bullock et al. 2001). This implies that ~29 per cent of the DM is within the 15.2 kpc projected radius of system C, and that the maximum escape velocity of system C from A is $v_{\text{esc}} = (2GM_{\text{tot}}/r_{\text{proj,AC}})^{1/2} \sim 130$ km s$^{-1}$. This is only ~62 per cent of the (orthogonal) relative line-of-sight velocity of $\Delta v_{\text{los,AC}} = -209.86$ km s$^{-1}$, suggesting that the system may not be gravitationally bound and that system C could be in the process of a fly-by past system A.

System B has the largest $H\alpha$ line EW, indicating that it contains the youngest stellar population. In addition, its weak stellar continuum suggests that system B is not a galaxy but a new-born bright H II region, possibly formed due to the interaction between galaxies A and C. This is similar to the Antennae galaxies (NGC 4038/NGC 4039) where many H II regions with high star-formation are distributed around the off-nuclear regions (Wang et al. 2004). This kind of starbursting clump is also seen in other galaxy merger systems. For example, in tidal features and very young stellar populations (e.g. Hancock et al. 2009; Levesque & Leitherer 2013).

### 7.3. The stellar population of complex C

We found a young stellar population of ~7 Myr at the SN explosion site (position C). Thöne et al. (2015) found a very young stellar population of ~3 Myr at the SN site of PTF12dam, and furthermore linked this to a quantitative estimate of 60 $M_{\odot}$ for the SN progenitor. However, these data and analysis do not directly prove that the progenitor was of this age and mass. As discussed in Chen et al. (2015), the bulk of the stellar population in the host of PTF12dam spans ages between 10 and 40 Myr. There is certainly evidence of a younger population of massive stars (a few Myr old) existing from the emission lines, but this does not indicate that the progenitor was one of them. The same argument should apply to the host system of LSQ14mo. The spatial resolution of the spectra and images are seeing dominated and $1''$ corresponds to 4 kpc at $z = 0.256$. Therefore a direct age estimate is not reliable from ground-based spectra, and a quantitative mass is not supported as a unique turn off mass to apply for the progenitor mass of the SLSN.

The very deep spectrum of the host system of LSQ14mo allows us to search for Wolf-Rayet (WR) star features, which are a good tracer of recent massive star formation (Kunth & Sargent 1981). There are two main WR features (see e.g. Brinchmann et al. 2008) seen in optical spectra: the blue (4600–4680 Å) and red (5650–5800 Å) bumps. However, we do not find these features in our spectra for components A, B, or C. We followed the method used in Thöne et al. (2015) to determine the upper limit of WR star populations. We measured
the $3\sigma$ detection limit of He\,{\sc ii} λ4686 and C IV λ5808 lines assuming a full width at half maximum (FWHM) between 1000 and 3000 km s$^{-1}$ (Brinchmann et al. 2008). The luminosity of He\,{\sc ii} line is <9.29 × 10$^{38}$–2.80 × 10$^{39}$ erg s$^{-1}$, and that of C IV line is <1.29 × 10$^{39}$–3.86 × 10$^{39}$ erg s$^{-1}$. We also considered the metallicity-luminosity relation for WR stars using Eqs. (7) and (8) in López-Sánchez & Esteban (2010) and the host metallicity of 12 + log(O/H) ~ 8.2 (i.e. 0.3 solar). We obtained the upper limit of <760–2280 of late-type WN stars, and <670–2000 of early-type WC stars.

The lack of WR star features may indicate that either there is a very young stellar population environment, which is younger than 2 Myr before WR stars evolved, or it may imply an older environment, such that WR stars have already died within 5 Myr after a star-forming episode. The Hα EW indication of a 7 Myr old stellar population at LSQ14mo host galaxy appears consistent with the latter explanation. In addition, there is no narrow galaxy He\,{\sc ii} λ4686 line detected in the entire host system of LSQ14mo. To date, only the host galaxy of PTF12dam has been reported to have this line detection. However, unambiguous WR features are not been detected (Chen et al. 2015; Thöne et al. 2015).

## 8. Conclusions

LSQ14mo is a spectroscopically normal type Ic SLSN at $z = 0.2563$, with a fast-declining light curve. We have presented its photometric and spectroscopic evolution from the early rising phase through to approximately +70 d after maximum light. Combining our multicolour data and K-correction from spectra, we were able to build the bolometric light curve, which we fitted with magnetar and CSM interaction-powered models. As with other SLSNe, both scenarios can satisfactorily reproduce the observations. Using SN 1998bw as a template, we find that our upper limit of $2 M_{\odot}^{56}$Ni mass derived for LSQ14mo is not sufficient to power the peak luminosity of LSQ14mo, and thus we ruled out the purely $^{56}$Ni-powered mechanism. We fit our spectra of LSQ14mo with a synthetic spectral model (Mazzali et al. 2016). The early spectra can be reproduced assuming a magnetar-like energy source, but residuals from the later fits suggest that there may be a contribution from a partially thermalised, additional source of luminosity, starting after day −7 d and ending before +22 d. The spectrum at +9 d is the most heavily affected, which shows a hot continuum but rather cool lines. We propose that this additional component could be produced by the interaction between the SN ejecta and a thin shell of material that was likely released during a pre-explosion mass-loss event of the massive progenitor star.

Late-time imaging of the SN revealed a host galaxy system with three resolved components. This was confirmed by deep VLT+FORS2 spectroscopy, which showed emission lines at three positions that are all spatially distinct and at slightly different velocities. LSQ14mo is coincident with the position we designated “C”, which is the furthest from the brightest component in the system (designated “A”), at a velocity offset of ≈200 km s$^{-1}$. This is comparable to the relative velocities of the Magellanic Clouds relative to the Milky Way. The physical diameter of component C (A) was found to be 2.6 kpc (6.2 kpc), and stellar mass from the SED fitting was found to be $10^{7.7} M_{\odot} (10^{8.8} M_{\odot})$, which is also comparable to the Magellanic Clouds. By resolving the SN location in spectroscopy and imaging, we were able to carry out a detailed analysis comparing the SN site to the global properties of the host system. We found that the SN location had the lowest oxygen abundance of

the three resolved components, with $12+\log(O/H) = 7.92–8.18$, depending on the metallicity diagnostics used. This is ~0.4 dex lower than the metallicity at position A. On the other hand, we measured a fairly constant specific SFR across the host system, and the specific SFR is typical compared to the general star-forming dwarf galaxies at similar redshift. Finally, we derive a young stellar population (6–9 Myr) at all three positions as traced by the Hα line width. We propose that components C and A are an interacting galaxy system that has triggered star-formation, and that component B is a new-born H II region, which thus shows the strongest EW of Hα line and the weakest stellar continuum.

The low metallicity constraint is favoured by the magnetar model as the power source for LSQ14mo. The spin-down period (3.9 ms) and the host galaxy metallicity ($12+\log(O/H) = 8.18 \pm 0.05$, N2 scale) of LSQ14mo is in agreement with the magnetar spin-period and host metallicity relation proposed by Chen et al. (2016).

**Acknowledgements.** T.-W. Chen thanks the organisers and participants of MIAPP workshop “Supernovae: The Outliers” for stimulating discussions: You-Hua Chu for interacting galaxy, Maryam Modjaz for metallicity and Ken'ichi Nomoto for CSM model. T.-W. Chen thanks Stefan Taubenberger for bolometric light curve comparison and participants in the Garching SN meeting for useful advise. T.W.C, T.K., R.Y. and P.W. acknowledge the support through the Sofia Kovalevskaja Award to P. Schady from the Alexander von Humboldt Foundation of Germany. S.J.S. acknowledges ERC grant 291222 and STFC grants ST/IO01123/1 and ST/L000709/1. L.G. was supported in part by the US National Science Foundation under Grant AST-1311862. A.J. acknowledges funding by the European Union's Framework Programme for Research and Innovation Horizon 2020 under Marie Skłodowska-Curie grant agreement No. 702538. M.S. acknowledges support from EU/FP7-ERC grant 615929. K.M. acknowledges support from the STFC through an Ernest Rutherford Fellowship. The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias with financial support from the UK Science and Technology Facilities Council. Part of the funding for GROND (both hardware as well as personnel) was generously granted from the Leibniz-Prize to Prof. G. Hasinger (DFG grant HA 1850/28-1). We acknowledge use of the Pan-STARRS1 photometric reference catalogue: the Pan-STARRS1 Surveys (PS1) have been made possible through contributions of the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, Queen's University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation under Grant No. AST-1238877, the University of Maryland, and Eotvos Lorand University (ELTE) and the Los Alamos National Laboratory.

## References

* Angus, C. R., Levan, A. J., Perley, D. A., et al. 2016, MNRAS, 458, 84

* Appenzeller, I., Fricke, K., Fürtig, W., et al. 1998, The Messenger, 94, 1

* Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5

* Baltay, C., Rabinowitz, D., Hadjiyska, E., et al. 2013, PASP, 125, 683

* Berg, D. A., Skillman, E. D., Marble, A. R., et al. 2012, ApJ, 754, 98

* Bianco, F. B., Modjaz, M., Oh, S. M., et al. 2016, Astron. Comput., 16, 54

* Bresolin, F. 2011, ApJ, 729, 56

* Brinchmann, J., Charlot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151

* Brinchmann, J., Kunth, D., & Durret, F. 2008, A&A, 485, 657

* Brüns, C., Kerp, J., Staveley-Smith, L., et al. 2005, A&A, 432, 45

* Bruzual, G. 2007, in From Stars to Galaxies: Building the Pieces to Build Up the Universe, eds. A. Vallenari, R. Tantalo, L. Portinari, & A. Moretti, ASP Conf. Ser., 374, 303

* Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000

* Bullock, J. S., Kolatt, T. S., Sigad, Y., et al. 2001, MNRAS, 321, 559

* Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245

* Chabrier, G. 2003, PASP, 115, 763
Chatzopoulos, E., Wheeler, J. C., & Vinko, J. 2012, ApJ, 746, 121

Chen, T.-W., Smartt, S. J., Bresolin, F., et al. 2013, ApJ, 763, L28

Chen, T.-W., Smartt, S. J., Jerkstrand, A., et al. 2015, MNRAS, 452, 1567

Chen, T.-W., Smartt, S. J., Yates, R. M., et al. 2016, J. Comput. Phys., 327, 1

Chevalier, R. A., & Irwin, C. M. 2011, ApJ, 729, L6

Chomiuk, L., Chornock, R., Soderberg, A. M., et al. 2011, ApJ, 743, 114

Cid Fernandes, R., Mateus, A., Sodré, L., Stasińska, G., & Gomes, J. M. 2005, MNRAS, 358, 363

Cid Fernandes, R., Schoenell, W., Gomes, J. M., et al. 2009, Rev. Mex. Astron. Astrofis., 35, 127

Cooke, J., Sullivan, M., Gal-Yam, A., et al. 2012, Nature, 491, 228

Crowther, P. A. 2007, ARA&A, 45, 177

da Cunha, E., Charlot, S., & Elbaz, D. 2008, MNRAS, 388, 1595

Deng, J., Kawabata, K. S., Ohyama, Y., et al. 2004, ApJ, 605, L37

Dessart, L., Hillier, D. J., Waldman, R., Livne, E., & Blondin, S. 2012, MNRAS, 426, L76

Dessart, L., Waldman, R., Livne, E., Hillier, D. J., & Blondin, S. 2013, MNRAS, 428, 3227

Dexter, J., & Kasen, D. 2013, ApJ, 772, 30

Dopita, M. A., Kewley, L. J., Sutherland, R. S., & Nicholls, D. C. 2016, Ap&SS, 361, 61

Ellison, S. L., Patton, D. R., Simard, L., & McConnachie, A. W. 2008, AJ, 135, 1877

Gal-Yam, A. 2012, Science, 337, 927

Gal-Yam, A., Mazzali, P., Ofek, E. O., et al. 2009, Nature, 462, 624

Galama, T. J., Vreeswijk, P. M., van Paradijs, J., et al. 1998, Nature, 395, 670

Galbany, L., Stanishev, V., Mourão, A. M., et al. 2014, A&A, 572, A38

Galbany, L., Anderson, J. P., Rosales-Ortega, F. F., et al. 2016, MNRAS, 455, 4087

Gonzalez-Delgado, R. M., Perez, E., Tenorio-Tagle, G., et al. 1994, ApJ, 437, 239

Greiner, J., Bornemann, W., Clemens, C., et al. 2008, PASP, 120, 405

Guillochon, J., Parrent, J., Kelley, L. Z., & Margutti, R. 2017, ApJ, 835, 64

Guseva, N. G., Bulladeros, P., Meyer, H. T., Izotov, Y. I., & Fricke, K. J. 2009, A&A, 505, 63

Hachinger, S., Mazzali, P. A., Taubenberger, S., et al. 2012, MNRAS, 427, 2057

Hancock, M., Smith, B. J., Struck, C., Giroux, M. L., & Hurlock, S. 2009, AJ, 137, 4643

Heger, A., & Woosley, S. E. 2002, ApJ, 567, 532

Howell, D. A., Kasen, D., Lidman, C., et al. 2013, ApJ, 779, 98

Inserra, C., & Smartt, S. J. 2014, ApJ, 796, 87

Inserra, C., Smartt, S. J., Jerkstrand, A., et al. 2013, ApJ, 770, 128

Inserra, C., Bulla, M., Sim, S. A., & Smartt, S. J. 2016, ApJ, 831, 79

Inserra, C., Smartt, S. J., Gall, E. E. E., et al. 2017, ApJ, submitted [arXiv:1604.01226]

Izotov, Y. I., Guseva, N. G., & Thuan, T. X. 2011, ApJ, 728, 161

Japelj, J., Vergani, S. D., Salvaterra, R., Hunt, L. K., & Mannucci, F. 2016, A&A, 593, A115

Jerkstrand, A., Smartt, S. J., & Heger, A. 2016, MNRAS, 455, 3207

Jerkstrand, A., Smartt, S. J., Inserra, C., et al. 2017, ApJ, 835, 13

Kasen, D., & Bildsten, L. 2010, ApJ, 717, 245

Kasen, D., Woosley, S. E., & Heger, A. 2011, ApJ, 734, 102

Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189

Kennicutt, Jr., R. C., Lee, J. C., Funes, J. G., et al. 2008, ApJS, 178, 247

Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183

Kewley, L. J., Geller, M. J., & Barton, E. J. 2006, AJ, 131, 204

Kobulnicky, H. A., & Kewley, L. J. 2004, ApJ, 617, 240

Kozyreva, A., Gilmer, M., Hirschi, R., et al. 2017, MNRAS, 464, 2854

Kuncarayakti, H., Doi, M., Aldering, G., et al. 2013, AJ, 146, 30

Kuncarayakti, H., Galbany, L., Anderson, J. P., Krühler, T., & Hamuy, M. 2016, A&A, 593, A78

Kunth, D., & Sargent, W. L. W. 1981, A&A, 101, L5

Lee, J. C., Salzer, J. J., & Melbourne, J. 2004, ApJ, 616, 752

Leloudas, G., Gallazzi, A., Sollerman, J., et al. 2011, A&A, 530, A95

Leloudas, G., Chatzopoulos, E., Dilday, B., et al. 2012, A&A, 541, A129

Leloudas, G., Ergon, M., Taddia, F., et al. 2014, ATel, 5839

Leloudas, G., Patat, F., Maund, J. R., et al. 2015a, ApJ, 815, L10

Leloudas, G., Schulze, S., Krühler, T., et al. 2015b, MNRAS, 449, 917

Levesque, E. M., & Leitherer, C. 2013, ApJ, 779, 170

Liang, Y. C., Hammer, F., Yin, S. Y., et al. 2007, A&A, 473, 411

López-Sánchez, Á. R., & Esteban, C. 2010, A&A, 516, A104

Lucy, L. B. 2000, MNRAS, 318, 92

Lunnan, R., Chornock, R., Berger, E., et al. 2013, ApJ, 771, 97

Lunnan, R., Chornock, R., Berger, E., et al. 2014, ApJ, 787, 138

Lunnan, R., Chornock, R., Berger, E., et al. 2015, ApJ, 804, 90

Maeda, K., Nomoto, K., Mazzali, P. A., & Deng, J. 2006, ApJ, 640, 854

Magnier, E. A., Schlafly, E., Finkbeiner, D., et al. 2013, ApJS, 205, 20

Mazzali, P. A., & Lucy, L. B. 1993, A&A, 279, 447

Mazzali, P. A., Iwamoto, K., & Nomoto, K. 2000, ApJ, 545, 407

Mazzali, P. A., Walker, E. S., Pian, E., et al. 2013, MNRAS, 432, 2463

Mazzali, P. A., Sullivan, M., Hachinger, S., et al. 2014, MNRAS, 439, 1959

Mazzali, P. A., Sullivan, M., Pian, E., Greiner, J., & Kann, D. A. 2016, MNRAS, 458, 3455

McCrum, M., Smartt, S. J., Kotak, R., et al. 2014, MNRAS, 437, 656

Michel-Dansac, L., Lambas, D. G., Alonso, M. S., & Tissera, P. 2008, MNRAS, 386, L82

Moriya, T. J., Liu, Z.-W., Mackey, J., Chen, T.-W., & Langer, N. 2015, A&A, 584, L5

Moriya, T. J., Chen, T.-W., & Langer, N. 2017, ApJ, 835, 177

Moster, B. P., Somerville, R. S., Maulbetsch, C., et al. 2010, ApJ, 710, 903

Moster, B. P., Naab, T., & White, S. D. M. 2013, MNRAS, 428, 3121

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493

Neill, J. D., Sullivan, M., Gal-Yam, A., et al. 2011, ApJ, 727, 15

Nicholl, M., & Smartt, S. J. 2016, MNRAS, 457, L79

Nicholl, M., Smartt, S. J., Jerkstrand, A., et al. 2013, Nature, 502, 346

Nicholl, M., Smartt, S. J., Jerkstrand, A., et al. 2014, MNRAS, 444, 2096

Nicholl, M., Smartt, S. J., Jerkstrand, A., et al. 2015a, MNRAS, 452, 3869

Nicholl, M., Smartt, S. J., Jerkstrand, A., et al. 2015b, ApJ, 807, L18

Nicholl, M., Berger, E., Margutti, R., et al. 2016a, ApJ, 828, L18

Nicholl, M., Berger, E., Smartt, S. J., et al. 2016b, ApJ, 826, 39

Ofek, E. O., Cameron, P. B., Kasliwal, M. M., et al. 2007, ApJ, 659, L13

Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and active galactic nuclei (Mill Valley, CA: University Science Books)

Pastorello, A., Smartt, S. J., Botticella, M. T., et al. 2010, ApJ, 724, L16

Patat, F., Cappellaro, E., Danziger, J., et al. 2001, ApJ, 555, 900

Peeples, M. S., Pogge, R. W., & Stanek, K. Z. 2009, ApJ, 695, 259

Perley, D. A., Quimby, R. M., Yan, L., et al. 2016, ApJ, 830, 13

Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59

Pilyugin, L. S., Thuan, T. X., & Vílchez, J. M. 2003, A&A, 397, 487

Poole, T. S., Breeveld, A. A., Page, M. J., et al. 2008, MNRAS, 383, 627

Quimby, R. M., Kulkarni, S. R., Kasliwal, M. M., et al. 2011, Nature, 474, 487

Sawala, T., Frenk, C. S., Fattahi, A., et al. 2015, MNRAS, 448, 2941

Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103

Schlafly, E. F., Finkbeiner, D. P., Jurić, M., et al. 2012, ApJ, 756, 158

Schulze, S., Krühler, T., Leloudas, G., et al. 2016, ArXiv e-prints [arXiv:1612.05978]

Smartt, S. J., Valenti, S., Fraser, M., et al. 2015, A&A, 579, A40

Smith, M., Sullivan, M., D'Andrea, C. B., et al. 2016, ApJ, 818, L8

Smith, N., Li, W., Foley, R. J., et al. 2007, ApJ, 666, 1116

Sorokina, E., Blinnikov, S., Nomoto, K., Quimby, R., & Tolstov, A. 2016, ApJ, 829, 17

Stanishev, V., Rodrigues, M., Mourão, A., & Flores, H. 2012, A&A, 545, A58

Steele, I. A., Smith, R. J., Rees, P. C., et al. 2004, in Ground-based Telescopes, ed. J. M. Oschmann, Jr., SPIE Conf. Ser., 5489, 679

Stoll, R., Prieto, J. L., Stanek, K. Z., et al. 2011, ApJ, 730, 34

Thöne, C. C., de Ugarte Postigo, A., García-Benito, R., et al. 2015, MNRAS, 451, L65

Thuan, T. X., Goehring, K. M., Hibbard, J. E., Izotov, Y. I., & Hunt, L. K. 2016, MNRAS, 463, 4268

Tonry, J. L., Stubbs, C. W., Lykke, K. R., et al. 2012, ApJ, 750, 99

Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898

van Dokkum, P. G. 2001, PASP, 113, 1420

Vreeswijk, P. M., Savaglio, S., Gal-Yam, A., et al. 2014, ApJ, 797, 24

Vreeswijk, P. M., Leloudas, G., Gal-Yam, A., et al. 2017, ApJ, 835, 58

Wang, S. Q., Wang, L. J., Dai, Z. G., & Wu, X. F. 2015, ApJ, 799, 107

Wang, Z., Fazio, G. G., Ashby, M. L. N., et al. 2004, ApJS, 154, 193

Woosley, S. E. 2010, ApJ, 719, L204

Wright, E. L. 2006, PASP, 118, 1711

Yan, L., Quimby, R., Ofek, E., et al. 2015, ApJ, 814, 108

Yaron, O., & Gal-Yam, A. 2012, PASP, 124, 668

Yates, R. M., Kauffmann, G., & Guo, Q. 2012, MNRAS, 422, 215
# Appendix A: Additional details of data reduction

LSQ detected the SN during the early rising phase (−7.8 rest-frame days relative to the *r*-band maximum light on MJD = 56697). We investigated pre-discovery LSQ images and found even earlier detections of this object at −22.2 d. Although the LSQ pipeline reports magnitudes of the SN at each epoch at which it is detected, we opted to carry out our own point-spread function (PSF)-fitting photometry on every image in order to maintain consistency between epochs at which the pipeline did and did not detect SN flux automatically. The LSQ survey observations are carried out in a wide optical filter, effectively *g* + *r*. We calibrated our PSF photometry to *r* band using a sequence of local field stars. We calculated the LSQ−*r* colour on each night with a pipeline detection, and checked this against the expected LSQ−*r* colour derived from synthetic photometry of the spectra of LSQ14mo, finding the two to be consistent. The calibrated photometry in the LSQ system was found to be very close to *r* band, with LSQ−*r* < 0.02 mag.

The LT images were reduced automatically by the LT pipeline, which applies de-bias and flat-field corrections. We extracted the SN flux using PSF-fitting, and calibrated the magnitudes using local sequence stars, which were themselves calibrated against standard fields observed with IO:O on photometric nights. We also obtained one epoch from Magellan + IMACS, reduced with bias subtraction and flat-fielding. The NTT + EFOSC2 data were reduced using the PESSTO pipeline (Smartt et al. 2015) with bias subtraction and flat-fielding using twilight sky frames. Zeropoints were determined with nine Pan-STARRS1 catalogue sequence stars in the field (Schlaflay et al. 2012; Magnier et al. 2013) and converted to SDSS photometric system (Tonry et al. 2012), and PSF photometry was carried out using the DAOPHOT task within IRAF. For the last two epochs of *g* band from the NTT + EFOSC2, we have applied a band-pass correction of +0.31 mag from the *Gunn* *g* to the SDSS *g* (Inserra & Smartt 2014).

After +250 d, PESSTO took deep images to search for any late SN signal (e.g. as in Inserra et al. 2013). We used *i*-band data taken at +339.3 d and *r*-band data taken at +321.7 d as templates for image subtraction from the late epochs. The *R*-filter was misused for two epochs, and thus we used *R*-band data taken at +339.3 d as a template to subtract the *R*-band image taken at +252.6 d. We aligned the images using the GEOMAP and GEO-TRAN tasks within IRAF and employed the High Order Transform of PSF And Template Subtraction (HOTPANTS)<fnref n="9" /> package to carry out the subtraction. We did not detect a point source at the SN position from +252.6 d to +321.7 d, which indicates that either there is no SN signal or the SN does not change its brightness during these epochs. We assumed that the templates do not contain SN signal, as the measured host magnitude remains stable over these late epochs and no SN features are detected in deep (~5 h of VLT exposure) VLT spectra obtained during the same phase (see Sect. 2.3). We therefore measured 3σ limiting magnitudes on the late-time images, by adding ten artificial stars at ~0.5' around the SN position, and re-measuring their magnitudes until we obtained a standard deviation of 0.3 mag.

We obtained pure host images on 2015 December 17 (+539.2 d) using the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND; Greiner et al. 2008), a 7-channel imager that collects multi-band photometry simultaneously with *g'*r'*i'*z'JHK bands, mounted at the 2.2 m MPG telescope at ESO La Silla Observatory in Chile. We used different aperture sizes

to measure each component in each band. The *griz* photometry was calibrated against Pan-STARRS1 sequence stars and converted to the SDSS photometric system; while *JHK* was calibrated against the 2MASS catalogue and converted to the AB system. We noticed that the *g'*-band image has a higher light dispersion and the FWHM of stars are larger than other filters. This affected our flux measurement of the host components since they are not resolved clearly. We compared GROND photometry to the host magnitude measurement from the NTT + EFOSC *i*-band image, because different seeing conditions cause differences in the PSF, and thus different aperture sizes were needed. The two datasets are in agreement with a systematic uncertainty of approximately 0.1 mag.

For our host spectra, FORS2, with the standard resolution collimator, provides a pixel scale of 0.25"/pixel. We used an aperture size of 1.25" for trace A and a smaller aperture of 0.5" for the spatially-close traces B and C to avoid including fluxes from one to the other. Traces A and C are clearly identified from the emission lines and the stellar continuum. No continuum trace is visible at position B, but the nebular emission lines are clear; we therefore used the continuum trace of C as a reference to define the trace for position B. We employed a module of PYSYN-PHOT from PYTHON to calculate synthetic photometric quantities for the extracted, flux-calibrated spectra A, B, and C in the NTT + EFOSC2 *i*-band bandpass, and then scaled the fluxes of spectra to match their *i*-band photometry (Sect. 2.2), by a scaling factor of 2.72, 1.46, and 1.86 for spectra A, B and C, respectively.

# Appendix B: Emission line flux measurement

The spectrum A has strong stellar absorptions superimposed on the Balmer lines. Spectra B and C suffer from the same effect, but to a lesser extent (see Fig. 2). Hence, we employed the spectral synthesis code STARLIGHT (Cid Fernandes et al. 2005, 2009) to fit our observed spectra. We used the SSP model base of Bruzual (2007) and selected 66 models from it. This selection base has been applied in Stanishev et al. (2012), Galbany et al. (2014, 2016), which consists of 66 components with 17 different ages (from 1 Myr to 18 Gyr) and four metallicities (0.2, 0.4, 1.0 and 2.5 $Z_\odot$, where $Z_\odot = 0.02$) coming from a slightly modified version of the models of Bruzual & Charlot (2003) and uses a Chabrier (2003) IMF. In general, the model fittings are good for three components with a $\chi^2$ of 2.22, 0.03, and 0.08 for positions A, B, and C, respectively. The modelled spectra A, B, and C were subtracted from our observations and thus recover the true Balmer emission line fluxes (Fig. 2).

The spectrum A has the highest signal-to-noise ratio (S/N), and the main optical emission lines were all clearly detected; especially an intrinsic weak auroral [O III] $\lambda$4363 line is marginally detected at 3.6σ. The spectrum B has the weakest continuum but the largest emission line equivalent widths. The spectrum C, where the SN is located, shows no sign of a broad [O I] $\lambda$6300 line, which would generally be the strongest type Ic SN feature in the nebular phase. Therefore, there is no evidence for supernova flux from LSQ14mo at this epoch. No He II $\lambda$4686 line is detected at any position in the host system.

Line flux measurements were made after subtraction of the STARLIGHT synthesis models. We fitted Gaussian line profiles using the QUB custom built PROCSPEC environment within IDL. The FWHM for each line in a given spectrum was fixed to that of the [O III] $\lambda$4959 line. The line width is fairly closely matched with the instrumental resolution (~9 Å in the observed-frame) measured from the skylines. We further normalised the spectrum and determined the equivalent widths (EWs) of all identified

<footnote n="9">⁹ http://www.astro.washington.edu/users/becker/v2.0/hotpants.html</footnote>
lines. Uncertainties were estimated from the EW and the root mean square (rms) of the continuum according to the equation from Gonzalez-Delgado et al. (1994):

$$ \sigma_i = \sigma_c \sqrt{N + \frac{W_i}{\Delta}}, \qquad (B.1) $$

where $\sigma_i$ is the error on the flux of the emission line, $\sigma_c$ is the rms measured from the local continuum, $N$ is the measured line profile width in pixels, $W_i$ is the absolute value of the line EW in Å, and $\Delta$ is the spectral dispersion in Å/pixel. The flux measurements and related parameters are listed in Table C.3.
## Appendix C: Observations log

<table><caption>Table C.1. Log of optical imaging of LSQ14mo.</caption><thead><tr><th>Date</th><th>MJD</th><th>Phase</th><th>g (error)</th><th>r (error)</th><th>i (error)</th><th>Telescope</th><th>Instrument</th></tr></thead><tbody><tr><td>2014 Jan. 02</td><td>56 659.23</td><td>-30.1</td><td>-</td><td>&gt;20.94</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 04</td><td>56 661.22</td><td>-28.5</td><td>-</td><td>&gt;21.43</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 06</td><td>56 663.16</td><td>-26.9</td><td>-</td><td>&gt;20.28</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 08</td><td>56 665.16</td><td>-25.4</td><td>-</td><td>&gt;20.36</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 10</td><td>56 667.16</td><td>-23.8</td><td>-</td><td>&gt;20.79</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 12</td><td>56 669.15</td><td>-22.2</td><td>-</td><td>20.82 (0.22)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 14</td><td>56 671.14</td><td>-20.5</td><td>-</td><td>21.04 (0.29)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 16</td><td>56 673.14</td><td>-19.0</td><td>-</td><td>20.64 (0.22)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 18</td><td>56 675.13</td><td>-17.4</td><td>-</td><td>20.41 (0.26)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 20</td><td>56 677.12</td><td>-15.8</td><td>-</td><td>19.95 (0.20)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 22</td><td>56 679.12</td><td>-14.2</td><td>-</td><td>19.99 (0.12)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 24</td><td>56 681.14</td><td>-12.6</td><td>-</td><td>19.77 (0.21)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 28</td><td>56 685.12</td><td>-9.4</td><td>-</td><td>19.92 (0.16)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Jan. 30</td><td>56 687.17</td><td>-7.8</td><td>-</td><td>19.49 (0.14)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Feb. 1</td><td>56 689.09</td><td>-6.3</td><td>-</td><td>19.40 (0.24)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Feb. 3</td><td>56 691.09</td><td>-4.7</td><td>-</td><td>19.49 (0.11)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Feb. 9</td><td>56 697.01</td><td>0.0</td><td>19.48 (0.04)</td><td>19.42 (0.05)</td><td>19.57 (0.06)</td><td>LT</td><td>IO</td></tr><tr><td>2014 Feb. 9</td><td>56 697.07</td><td>0.1</td><td>-</td><td>19.42 (0.19)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Feb. 10</td><td>56 698.03</td><td>0.8</td><td>19.52 (0.03)</td><td>19.50 (0.07)</td><td>19.53 (0.07)</td><td>LT</td><td>IO</td></tr><tr><td>2014 Feb. 10</td><td>56 698.98</td><td>1.6</td><td>19.54 (0.04)</td><td>19.42 (0.05)</td><td>19.49 (0.06)</td><td>LT</td><td>IO</td></tr><tr><td>2014 Feb. 12</td><td>56 700.01</td><td>2.4</td><td>19.52 (0.06)</td><td>19.50 (0.10)</td><td>19.50 (0.12)</td><td>LT</td><td>IO</td></tr><tr><td>2014 Feb. 23</td><td>56 711.07</td><td>11.2</td><td>20.13 (0.07)</td><td>19.74 (0.08)</td><td>19.74 (0.09)</td><td>LT</td><td>IO</td></tr><tr><td>2014 Feb. 26</td><td>56 714.07</td><td>13.6</td><td>-</td><td>19.78 (0.16)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Feb. 28</td><td>56 716.08</td><td>15.2</td><td>-</td><td>19.86 (0.15)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Mar. 02</td><td>56 718.13</td><td>16.8</td><td>-</td><td>19.85 (0.21)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Mar. 04</td><td>56 720.03</td><td>18.3</td><td>-</td><td>20.03 (0.19)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Mar. 06</td><td>56 722.06</td><td>20.0</td><td>-</td><td>20.32 (0.15)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Mar. 07</td><td>56 723.91</td><td>21.4</td><td>21.27 (0.20)</td><td>20.44 (0.07)</td><td>20.11 (0.11)</td><td>LT</td><td>IO</td></tr><tr><td>2014 Mar. 08</td><td>56 724.06</td><td>21.7</td><td>-</td><td>20.46 (0.16)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Mar. 10</td><td>56 726.14</td><td>23.2</td><td>-</td><td>20.65 (0.20)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Mar. 12</td><td>56 728.05</td><td>24.7</td><td>-</td><td>20.57 (0.19)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Mar. 12</td><td>56 728.91</td><td>25.4</td><td>21.48 (0.27)</td><td>20.80 (0.17)</td><td>20.41 (0.13)</td><td>LT</td><td>IO</td></tr><tr><td>2014 Mar. 18</td><td>56 734.13</td><td>29.4</td><td>-</td><td>21.21 (0.35)</td><td>-</td><td>ESO 1.0-m Schmidt</td><td>QUEST</td></tr><tr><td>2014 Mar. 20</td><td>56 736.88</td><td>31.8</td><td>21.78 (0.20)</td><td>20.97 (0.24)</td><td>20.54 (0.16)</td><td>LT</td><td>IO</td></tr><tr><td>2014 Mar. 24</td><td>56 740.23</td><td>34.4</td><td>-</td><td>21.11 (0.10)</td><td>-</td><td>Magellan</td><td>IMACS</td></tr><tr><td>2014 Apr. 21</td><td>56 768.02</td><td>56.5</td><td>23.77 (0.10)</td><td>22.42 (0.13)</td><td>22.07 (0.16)</td><td>NTT</td><td>EFOSC2</td></tr><tr><td>2014 May 6</td><td>56 783.01</td><td>68.5</td><td>24.14 (0.14)</td><td>22.85 (0.15)</td><td>22.32 (0.13)</td><td>NTT</td><td>EFOSC2</td></tr><tr><td>2014 Dec. 23</td><td>57 014.31</td><td>252.6</td><td>-</td><td>&gt;23.92</td><td>&gt;24.46</td><td>NTT</td><td>EFOSC2</td></tr><tr><td>2014 Dec. 23</td><td>57 014.31</td><td>252.6</td><td>-</td><td>&gt;23.79*</td><td>-</td><td>NTT</td><td>EFOSC2</td></tr><tr><td>2015 Jan. 21</td><td>57 043.31</td><td>275.7</td><td>-</td><td>&gt;23.96</td><td>&gt;24.04</td><td>NTT</td><td>EFOSC2</td></tr><tr><td>2015 Feb. 12</td><td>57 065.18</td><td>293.1</td><td>-</td><td>&gt;24.03</td><td>&gt;24.03</td><td>NTT</td><td>EFOSC2</td></tr><tr><td>2015 Mar. 20</td><td>57 101.07</td><td>321.7</td><td>-</td><td>&gt;23.63</td><td>&gt;23.34</td><td>NTT</td><td>EFOSC2</td></tr><tr><td>2015 Apr. 11</td><td>57 123.18</td><td>339.3</td><td>-</td><td>&gt;23.22*</td><td>&gt;23.24</td><td>NTT</td><td>EFOSC2</td></tr></tbody></table>

Notes. The phase (day) has been corrected for time dilation ($z = 0.256$) and is relative to the SN r-band maximum on MJD 56 697. The ">" denotes the $3\sigma$ detection limit. Optical gri magnitudes are in the AB system and parenthesized with statistical errors. Note: R-band magnitudes are marked with a * symbol.
<table><caption>**Table C.2.** Ultraviolet photometry of LSQ14mo in the *Swift* UVOT bands.</caption><thead><tr><th>Date</th><th>MJD</th><th>Phase</th><th>uvw2</th><th>uvm2</th><th>uvw1</th><th>u</th></tr></thead><tbody><tr><td>2014 Jan. 31</td><td>56 688.76</td><td>-6.6</td><td>-</td><td>21.49 (0.17)</td><td>20.80 (0.16)</td><td>-</td></tr><tr><td>2014 Feb. 02</td><td>56 690.03</td><td>-5.5</td><td>21.81 (0.26)</td><td>-</td><td>-</td><td>19.62 (0.17)</td></tr><tr><td>2014 Feb. 04</td><td>56 692.03</td><td>-4.0</td><td>-</td><td>-</td><td>-</td><td>19.41 (0.15)</td></tr><tr><td>2014 Feb. 05</td><td>56 693.05</td><td>-3.1</td><td>-</td><td>-</td><td>-</td><td>19.92 (0.21)</td></tr><tr><td>2014 Feb. 09</td><td>56 697.97</td><td>0.8</td><td>-</td><td>-</td><td>-</td><td>19.96 (0.20)</td></tr><tr><td>2014 Feb. 11</td><td>56 699.72</td><td>2.2</td><td>-</td><td>-</td><td>-</td><td>19.94 (0.20)</td></tr></tbody></table>

**Notes.** The phase (day) has been corrected for time dilation ($z = 0.256$) and is relative to the SN r-band maximum on MJD 56 697. *Swift* magnitudes are in the AB system, with the uncertainty including statistical and systematic errors.
<table><caption>**Table C.3.** Emission line measurements of the host galaxy system of LSQ14mo in the observed frame.</caption><thead><tr><th colspan="7">Position A</th></tr><tr><th>Line</th><th>λ (Å)</th><th>Flux</th><th>Error</th><th>rms</th><th>EW (Å)</th><th>Luminosity</th></tr></thead><tbody><tr><td>[O II]</td><td>3727</td><td>498</td><td>3</td><td>0.7</td><td>58</td><td>101</td></tr><tr><td>[Ne III]</td><td>3868</td><td>39.2</td><td>3.5</td><td>1.2</td><td>4.6</td><td>7.9</td></tr><tr><td>Hδ</td><td>4102</td><td>23.1</td><td>2.3</td><td>0.8</td><td>2.7</td><td>4.7</td></tr><tr><td>Hγ</td><td>4340</td><td>53</td><td>2</td><td>0.6</td><td>3</td><td>11</td></tr><tr><td>[O III]</td><td>4363</td><td>4.63</td><td>1.28</td><td>0.47</td><td>1.1</td><td>0.94</td></tr><tr><td>Hβ</td><td>4861</td><td>137</td><td>1</td><td>0.4</td><td>11</td><td>28</td></tr><tr><td>[O III]</td><td>4959</td><td>121</td><td>1</td><td>0.3</td><td>12</td><td>25</td></tr><tr><td>[O III]</td><td>5007</td><td>332</td><td>8</td><td>2.1</td><td>32</td><td>67</td></tr><tr><td>[O I]</td><td>6300</td><td>35</td><td>2</td><td>0.6</td><td>5</td><td>7</td></tr><tr><td>[N II]</td><td>6548</td><td>11.6</td><td>1.6</td><td>0.6</td><td>1.7</td><td>2.4</td></tr><tr><td>Hα</td><td>6563</td><td>425</td><td>3</td><td>0.6</td><td>68</td><td>86</td></tr><tr><td>[N II]</td><td>6583</td><td>48</td><td>2</td><td>0.6</td><td>8</td><td>10</td></tr><tr><td>[S II]</td><td>6717</td><td>93</td><td>3</td><td>0.9</td><td>15</td><td>19</td></tr><tr><td>[S II]</td><td>6731</td><td>60</td><td>3</td><td>0.9</td><td>10</td><td>12</td></tr></tbody></table>

<table><thead><tr><th colspan="7">Position B</th></tr><tr><th>Line</th><th>λ (Å)</th><th>Flux</th><th>Error</th><th>rms</th><th>EW (Å)</th><th>Luminosity</th></tr></thead><tbody><tr><td>[O II]</td><td>3727</td><td>74</td><td>0.8</td><td>0.13</td><td>90</td><td>15</td></tr><tr><td>[Ne III]</td><td>3868</td><td>12</td><td>0.4</td><td>0.13</td><td>9</td><td>2.4</td></tr><tr><td>Hδ</td><td>4102</td><td>7.4</td><td>0.33</td><td>0.11</td><td>4</td><td>1.5</td></tr><tr><td>Hγ</td><td>4340</td><td>13</td><td>0.2</td><td>0.06</td><td>14</td><td>2.7</td></tr><tr><td>[O III]</td><td>4363</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Hβ</td><td>4861</td><td>31</td><td>0.4</td><td>0.08</td><td>40</td><td>6.3</td></tr><tr><td>[O III]</td><td>4959</td><td>39</td><td>0.3</td><td>0.08</td><td>45</td><td>8.0</td></tr><tr><td>[O III]</td><td>5007</td><td>119</td><td>3.5</td><td>0.49</td><td>144</td><td>24</td></tr><tr><td>[O I]</td><td>6300</td><td>6.1</td><td>0.63</td><td>0.20</td><td>11</td><td>1.2</td></tr><tr><td>[N II]</td><td>6548</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Hα</td><td>6563</td><td>80</td><td>1.1</td><td>0.13</td><td>219</td><td>16</td></tr><tr><td>[N II]</td><td>6583</td><td>2.7</td><td>0.41</td><td>0.13</td><td>9</td><td>0.6</td></tr><tr><td>[S II]</td><td>6717</td><td>7.8</td><td>0.79</td><td>0.21</td><td>23</td><td>1.6</td></tr><tr><td>[S II]</td><td>6731</td><td>6.9</td><td>0.81</td><td>0.21</td><td>25</td><td>1.4</td></tr></tbody></table>

<table><thead><tr><th colspan="7">Position C</th></tr><tr><th>Line</th><th>λ (Å)</th><th>Flux</th><th>Error</th><th>rms</th><th>EW (Å)</th><th>Luminosity</th></tr></thead><tbody><tr><td>[O II]</td><td>3727</td><td>53</td><td>1.0</td><td>0.19</td><td>60</td><td>10.8</td></tr><tr><td>[Ne III]</td><td>3868</td><td>8.8</td><td>0.54</td><td>0.15</td><td>9</td><td>1.8</td></tr><tr><td>Hδ</td><td>4102</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Hγ</td><td>4340</td><td>9.3</td><td>0.43</td><td>0.13</td><td>4</td><td>1.9</td></tr><tr><td>[O III]</td><td>4363</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Hβ</td><td>4861</td><td>23</td><td>0.5</td><td>0.12</td><td>21</td><td>4.8</td></tr><tr><td>[O III]</td><td>4959</td><td>26</td><td>0.7</td><td>0.18</td><td>23</td><td>5.2</td></tr><tr><td>[O III]</td><td>5007</td><td>77</td><td>6.2</td><td>1.11</td><td>70</td><td>15.7</td></tr><tr><td>[O I]</td><td>6300</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>[N II]</td><td>6548</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Hα</td><td>6563</td><td>64</td><td>1.5</td><td>0.24</td><td>107</td><td>13.0</td></tr><tr><td>[N II]</td><td>6583</td><td>3.7</td><td>0.83</td><td>0.24</td><td>8</td><td>0.7</td></tr><tr><td>[S II]</td><td>6717</td><td>6.8</td><td>0.81</td><td>0.22</td><td>13</td><td>1.4</td></tr><tr><td>[S II]</td><td>6731</td><td>5.0</td><td>0.80</td><td>0.22</td><td>11</td><td>1.0</td></tr></tbody></table>

**Notes.** The spectra have been corrected for Milky Way reddening corresponding to a dust extinction of $A_V = 0.20$. We measured the emission line fluxes after subtracting the stellar continuum from the STARLIGHT spectral synthesis models. The error is derived by the Eq. (B.1). We fixed the FWHM widths at the values derived from the [O III] λ4959 line, which correspond to the FORS2 resolution (9.05, 9.06, and 9.21 Å for positions A, B, and C, respectively). The relativity large rms and error of [O III] λ5007 is due to a strong sky line nearby. For luminosity, we have not applied host galaxy dust reddening corrections. The internal host extinction at position A is $A_V = 0.26$ derived by the intrinsic ratio of $H\alpha$ and $H\beta$, and negligible at positions B and C. The unit of Flux, Error, rms of $\times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$, Luminosity of $\times 10^{39}$ erg s$^{-1}$.

<table><caption>**Table C.4.** K-corrections for LSQ14mo and its host galaxy system derived from the spectroscopic observations using SNAKE (Inserra et al. 2017).</caption><thead><tr><th>Date</th><th>MJD</th><th>Phase</th><th>K<sub>g→u</sub></th><th>K<sub>r→g</sub></th><th>K<sub>i→r</sub></th></tr></thead><tbody><tr><td colspan="6"><strong>LSQ14mo</strong></td></tr><tr><td>2014 Jan. 31</td><td>56 688.16</td><td>-7.0</td><td>0.20</td><td>0.26</td><td>0.29</td></tr><tr><td>2014 Feb. 6</td><td>56 694.19</td><td>-2.2</td><td>0.22</td><td>0.24</td><td>0.28</td></tr><tr><td>2014 Feb. 20</td><td>56 708.06</td><td>8.8</td><td>0.28</td><td>0.27</td><td>0.26</td></tr><tr><td>2014 Feb. 28</td><td>56 716.13</td><td>15.2</td><td>0.37</td><td>0.34</td><td>0.26</td></tr><tr><td>2014 Mar. 8</td><td>56 724.26</td><td>21.7</td><td>0.41</td><td>0.38</td><td>0.25</td></tr><tr><td>2014 Apr. 21</td><td>56 768.04</td><td>56.5</td><td>0.38</td><td>0.32</td><td>0.31</td></tr><tr><td colspan="6"><strong>Host galaxy</strong></td></tr><tr><td>2015 Jan.-Feb.</td><td>57 065.58</td><td>293.5</td><td></td><td></td><td></td></tr><tr><td>A</td><td></td><td></td><td>0.57</td><td>0.33</td><td>0.21</td></tr><tr><td>B</td><td></td><td></td><td>0.44</td><td>0.35</td><td>0.15</td></tr><tr><td>C</td><td></td><td></td><td>0.58</td><td>0.37</td><td>0.19</td></tr></tbody></table>

**Notes.** The phase (day) has been corrected for time dilation ($z = 0.256$) and is relative to the SN r-band maximum on MJD 56 697. The MJD and phase of the host galaxy was adopted from the middle time of three VLT observations.
<table><caption>**Table C.5.** Final and optical pseudo bolometric luminosities of LSQ14mo with their uncertainties.</caption><thead><tr><th>Phase<br>(day)</th><th>log L<sub>bol</sub><br>(erg/s)</th><th>log L<sub>pseudo</sub><br>(erg/s)</th></tr></thead><tbody><tr><td>-22.2</td><td>43.47 ± 0.05</td><td>43.16 ± 0.04</td></tr><tr><td>-20.5</td><td>43.38 ± 0.05</td><td>43.08 ± 0.04</td></tr><tr><td>-19.0</td><td>43.55 ± 0.04</td><td>43.24 ± 0.04</td></tr><tr><td>-17.4</td><td>43.64 ± 0.04</td><td>43.33 ± 0.04</td></tr><tr><td>-15.8</td><td>43.83 ± 0.04</td><td>43.52 ± 0.03</td></tr><tr><td>-14.2</td><td>43.81 ± 0.04</td><td>43.50 ± 0.03</td></tr><tr><td>-12.6</td><td>43.90 ± 0.04</td><td>43.59 ± 0.03</td></tr><tr><td>-9.4</td><td>43.84 ± 0.04</td><td>43.54 ± 0.03</td></tr><tr><td>-7.8</td><td>44.01 ± 0.04</td><td>43.71 ± 0.03</td></tr><tr><td>-7.0</td><td>44.03 ± 0.03</td><td>43.73 ± 0.03</td></tr><tr><td>-6.3</td><td>44.04 ± 0.03</td><td>43.75 ± 0.03</td></tr><tr><td>-4.7</td><td>44.03 ± 0.03</td><td>43.74 ± 0.03</td></tr><tr><td>-2.2</td><td>44.02 ± 0.03</td><td>43.74 ± 0.03</td></tr><tr><td>0.0</td><td>43.99 ± 0.03</td><td>43.72 ± 0.02</td></tr><tr><td>0.8</td><td>43.99 ± 0.03</td><td>43.72 ± 0.02</td></tr><tr><td>1.6</td><td>43.98 ± 0.03</td><td>43.73 ± 0.02</td></tr><tr><td>2.4</td><td>43.97 ± 0.02</td><td>43.72 ± 0.02</td></tr><tr><td>8.8</td><td>43.86 ± 0.03</td><td>43.64 ± 0.03</td></tr><tr><td>11.2</td><td>43.79 ± 0.02</td><td>43.61 ± 0.02</td></tr><tr><td>13.6</td><td>43.74 ± 0.02</td><td>43.58 ± 0.03</td></tr><tr><td>15.2</td><td>43.69 ± 0.02</td><td>43.56 ± 0.04</td></tr><tr><td>16.8</td><td>43.67 ± 0.02</td><td>43.55 ± 0.03</td></tr><tr><td>18.3</td><td>43.62 ± 0.02</td><td>43.50 ± 0.04</td></tr><tr><td>20.0</td><td>43.57 ± 0.02</td><td>43.44 ± 0.04</td></tr><tr><td>21.4</td><td>43.52 ± 0.02</td><td>43.39 ± 0.04</td></tr><tr><td>21.7</td><td>43.52 ± 0.02</td><td>43.38 ± 0.04</td></tr><tr><td>23.2</td><td>43.46 ± 0.02</td><td>43.34 ± 0.04</td></tr><tr><td>25.4</td><td>43.39 ± 0.03</td><td>43.27 ± 0.05</td></tr><tr><td>31.8</td><td>43.28 ± 0.02</td><td>43.17 ± 0.05</td></tr><tr><td>34.4</td><td>43.19 ± 0.02</td><td>43.07 ± 0.05</td></tr><tr><td>56.5</td><td>42.83 ± 0.03</td><td>42.60 ± 0.05</td></tr><tr><td>68.5</td><td>42.67 ± 0.04</td><td>42.41 ± 0.04</td></tr><tr><td>252.6</td><td>&lt;41.86</td><td>&lt;41.44</td></tr><tr><td>275.7</td><td>&lt;41.88</td><td>&lt;41.47</td></tr><tr><td>293.1</td><td>&lt;41.85</td><td>&lt;41.45</td></tr><tr><td>321.7</td><td>&lt;42.04</td><td>&lt;41.65</td></tr></tbody></table>

**Notes.** The SN phase is reported with respect to the time of observed maximum r-band light in the rest-frame. Note: The epoch of bolometric light curve corresponds to photometric epoch, except those epochs -7.0, -2.2 and +8.8, which were derived from the spectroscopy, and epochs 0.1 and 29.4 were dropped due to large photometry errors.