\section{Introduction}
Mrk\thinspace1239\ is a Narrow Line Seyfert 1  galaxy (NLS1) located at a redshift of $z=0.01993$ \citep{Beers+1995}. The optical spectra of Mrk\thinspace1239\ was first reported by \cite{Rafanelli+1984} and later used by \cite{Osterbrock+1985} to define the subclass of galaxies known as NLS1s. They suggest a classification criterion for NLS1 to have  \hbeta\ line widths less than 2000\,km~s$^{-1}$. This is now commonly adopted by the community \cite[see][for a complete review]{Komossa+2008}. These active galaxies are thought to be low mass counterparts of broad line Seyfert 1s  (BLS1). They show intense X-ray variability on both long-term and rapid time scales (see \citealt{Gallo+2018,Wilkins+2017, Wilkins+2015B}). These objects are seen through low column densities and often exhibit steep powerlaw continuum \citep{Waddell+2020} and strong soft excesses \citep{Waddell+2019}.

Mrk\thinspace1239\ is a radio quite galaxy that shows evidence of non-stellar radio emission \citep{Doi+2015}. This indicates the possible presence of a radio jet originating from the AGN in this galaxy. Mrk\thinspace1239\ is also a highly polarized in the broad \halpha\ and \hbeta\ emission lines, along with the continuum \citep{Pan+2019, Pan+2021}. The high levels of polarization in Mrk\thinspace1239\ have previously been observed by \cite{Goodrich+1989} where Mrk\thinspace1239\ was the highest polarized source in the sample. This suggests some amount of scattered emission in this source. Mrk\thinspace1239\ is classified as a polar-scattered Seyfert 1 galaxy, meaning the line of sight is through the upper layers of the tours \citep{Smith+2004,Jiang+2019}.
As well \cite{Goodrich+1989} found that the Balmer lines had a higher polarization than the forbidden \oiii\ lines. This is strong evidence that the sources of these emissions lines is physically distinct.

Mrk\thinspace1239\ was first observed in  X-ray using {\it ROSAT}\ where a steep photon index was found  \citep[$\Gamma \simeq 3$,][]{Rush+1996,Jiang+2021}.  Mrk\thinspace1239\ was a serendipitous source in a 2001 observation by {\it XMM-Newton}\ and analysed by \cite{Grupe+2004}. They found the continuum was best fit using a powerlaw with two intrinsic absorbers: the first a direct, highly absorbed path and the second a less absorbed scattered path. They concluded that this was consistent with the distinct emitting/scattering regions found in \cite{Goodrich+1989}.  Finally, \cite{Grupe+2004} found an excess of emission at 0.91 \kev, which was attributed to a blend of the \neix\ triplet.  However, the \ovii\ triplet was not observed in the spectrum suggesting a super-solar ratio of \NeOx\ in Mrk\thinspace1239.

Mrk\thinspace1239\ was revisited by \cite{Buhariwalla+2020} with archival X-ray data from {\it Suzaku}\ and \textit{NuSTAR}, obtained in the subsequent years since \cite{Grupe+2004}. With this broadband multi-epoch analysis  some truly unique behaviour  was revealed. First and foremost, across the 18 years of data analysed the shape and flux of the spectra below 3\,\kev\ remain remarkably similar. Meanwhile, the spectrum above 3\kev\ retained the characteristic NLS1 variability. This immediately indicated the presence of a strong absorbing feature that blotted out most of the emission from the central engine below 3\kev. The 0.91 \kev\ feature was present in the 2007 {\it Suzaku}\ observation of Mrk\thinspace1239.

Continuum modelling of this source between 0.3 and 30 \kev\ saw a ionized partial covering model compared to a blurred reflection model. The data were insufficient to significantly distinguish between the two models but their similarities proved to be more interesting. Both models required high levels of intrinsic absorption and both models required the presence of a collisionally ionized component to fit the spectrum below 3\kev. \cite{Buhariwalla+2020} attributed this collisionally ionized emission (CIE) to the presence of star forming activity in this object. They proposed that the 0.91 \kev\ feature found originally by \cite{Grupe+2004} was not an overabundance of Ne, but a blend of Fe-L lines from the regions of star formation. Using a relationship between the $L_{2-10\kev}$ of the CIE and star formation rate (SFR) \citep{Franceschini+2003}, the SFR of Mrk\thinspace1239\ was found to be $\simeq4-6\Msun \, {\rm yr^{-1}}$. This is comparable with the SFR of $<7.5\Msun \, {\rm yr^{-1}}$measured using PAH signatures \citep{RuschelDutra+2017},  $\sim3.5\Msun \, {\rm yr^{-1}}$ found using SED fitting \citep{Gruppioni+2016}, and $\sim2.1\Msun \, {\rm yr^{-1}}$ found using IR measurements \citep{Smirnova+2022}.  To fully probe the soft band of this galaxy and determine the origin of the 0.91 \kev\ feature conclusively deeper observations of Mrk\thinspace1239\ were needed.

This paper presents the first high-resolution X-ray spectrum of Mrk\thinspace1239\ in the $7-31$\,\AA\ band with the {\it XMM-Newton}\  RGS spectra. Limited MOS/PN spectra of a deep observation of Mrk\thinspace1239\ and a {\it Chandra}\ image is also presented. In Section \ref{sec:data}, the observations and data reduction techniques are summarized. Section \ref{sec:results} examines the {\it Chandra}\ image,  emission line identification,  possible outflow velocities, plasma diagnostics and the possibility of a \NeOx\  and in Section \ref{sec:spectral} the spectra are analyzed. A discussion of the results is given in Section \ref{sec:disc}, and conclusions are drawn in Section \ref{sec:conclusion}.
\section{Observations and Data Reduction}
\label{sec:data}
A snapshot observation of Mrk\thinspace1239\ was taken on 2021-04-11 with {\it Chandra}\ as part of a sample study exploring duel AGN \citep{Foord+2020}. Several months later a deep observation of Mrk\thinspace1239\ was taken on 2021-11-04 with {\it XMM-Newton}\ and \textit{NuSTAR}, with each obtaining $\approx100$ \ks\ of on source time. During the second half of the \textit{NuSTAR}\ observation Mrk\thinspace1239\ entered a flaring state. To fully explore all facets of this object analysis of the \textit{NuSTAR}\ and broad band {\it XMM-Newton}\ data is relegated to a second paper.  The data for analysis are listed in  Table \ref{tab:obs}. This section describes the observations and data reduction.
\begin{table*}
	\begin{tabular}{c c c c c c c c}
		\hline
		(1) 				  & (2) 			 			& (3) 		 						  & (4) 		  				& (5) 	 			& (6) 				& (7) 	 		& (8)				\\
		Observatory & Observation ID   &  Instrument Name 	 & Start date  			& Duration 	& Exposure 		& Counts & Energy range		\\
								 & 				   				  & 			 							& (yyyy-mm-dd)	& [s]      			& [s]      				&    	 		& 			 		\\
		\hline
		{\it XMM-Newton}			& 0891070101	  & PN 									 &2021-11-04 		 & 105000 	& 79175				& 1512	 &	$0.2-2.0$\,\kev	\\
		{\it XMM-Newton} 	 		& 				    			  & RGS 		 					   &2021-11-04  		& 105000 & $193963^{\star}$	&$5660$	 &  $7-31$\,\AA 	\\
		\hline
		{\it Chandra}     & 23693		            & ACIS-S 							& 2021-04-11 		& 2100 		& 2021					& 50 	 & $0.2-2$ \kev  	\\
		\hline
	\end{tabular}
	\caption{Observations log for Mrk 1239. The observations and  instruments used for analysis are listed in column (1). The observation ID and labels used in this work are given in columns (2) and (3), respectively. The start date of each observation is given in column (4). The duration of each observation, total exposure time and total counts for each observation are given in columns (5), (6), and  (7), respectively. The energy/wavelength each observation were fit over is given in column (8). $^{\star}$For RGS the combined exposure and counts for RGS1 and RGS2 are reported.}
	\label{tab:obs}
\end{table*}
\subsection{{\it XMM-Newton}}
Mrk\thinspace1239\ was observed with  {\it XMM-Newton}\ \citep{Jansen+2001} for 105 \ks\ starting on November $4^{th}$ 2021. The {\it XMM-Newton}\ Observation Data Files (ODF) were processed to produce calibrated event lists using the {\it XMM-Newton}\ Science Analysis System, {\sc sas}   V17.0.0.
The first order RGS data was extracted from the ODF using {\sc rgsproc}. Upon analysis of the light curve a background flaring event was seen within the first 15 \ks. {\sc tabgtigen} was used to filter all instances where the \cps\ was greater than or equal to 0.1. The spectra was then re-extracted using {\sc rgsproc}. RGS1 and RGS2 were then combined using  {\sc rgscombine}, which produced a combined source spectra, background spectra, and response matrix. The source and background spectra were then optimally binned using {\sc ftgrouppha} \citep{optbin}, the spectra were loaded into {\sc xspec} where the background was modelled. Fits were evaluated  using Cash-statistics \citep[C-stat,][]{Cash+1979}.

For the PN spectra an event list was created from the ODF using {\sc epproc}. The lightcurve was generated using {\sc evselect} and used to examine for background flaring. The same flare seen in the RGS data was present in the PN data, thus {\sc tabgtigen} was used to filter all times when the \cps\ were greater than or equal to 10. The spectra was also examined for pileup and was found to be satisfactory. The source spectra was extracted from a circular region 35"  centred on Mrk\thinspace1239. A background region was extracted from a circular region with a radius of 50" on the same chip. {\sc rmfgen} and {\sc arfgen} were used to generate response files for the observation.

\subsection{{\it Chandra}}
Mrk\thinspace1239\ was observed with {\it Chandra}\ in April 2021 for 2.1 \ks, and the data entered the archive in April 2022. The image was captured using  {\it Chandra}\ ACIS-S instrument and the data was processed using {\sc ciao 4.14}\footnote{\url{https://cxc.cfa.harvard.edu/ciao/}}. Initially a spectrum was extracted using {\sc dmextract}. The spectrum was then loaded into {\sc Sherpa} where a rough fit was obtained and saved over the energy region of interest, 0.2-2.0 \kev. The spectral shape obtained by this rough fit was imputed in {\sc ChaRT }\footnote{\url{https://cxc.cfa.harvard.edu/ciao/PSFs/chart2/runchart.html}} along with the nominal RA and DEC of the observation to generate a simulated PSF. The PSF was then projected in the detector plane using {\sc MARX} via {\sc simulate\_psf}  with a bin size of 1 pixel and a blur factor of 0.25. Next the image was exposure corrected using {\sc fluximage}, where an exposure corrected map centred around the brightest pixel in the image with a 600 pixel box, with bin size of 1 pixel. The image was extracted with a minimum energy of 0.2 \kev, a maximum energy of 2 \kev, and an effective energy of 1 \kev. Finally the image was deconvolved with the PSF generated  previously using {\sc arestore} with  the Lucy-Richardson deconvolution algorithm and 100 iterations. The finial image  presented is smoothed using a Gaussian interpolation method.
\section{Analysis}
\label{sec:results}
\subsection{Extended emission in 0.2-2 keV}
The 2.1\ks\ image of Mrk\thinspace1239\ taken with {\it Chandra}\ ACIS-S  is presented in Figure \ref{fig:june7_softChandra}. The radio contours are taken from \cite{Jarvela+2022} and are based on JVLA observations centred at 5.2 GHz with a 2 GHz bandwidth.  They  give a  sense of size and shape of the extended radio emission region.

Most of the soft X-ray emission originates from the very central region of the galaxy, suggesting at least part of it is originating from the AGN. Previous studies have indicated that there is star formation present within the central 400 pc of nucleus of Mrk\thinspace1239\ \citep{RuschelDutra+2017}. This appears comparable to the scale of the extended soft X-ray emission in the {\it Chandra}\ image  thus some of the soft emission seen in the central region may be the result of this star formation.
There appears to be asymmetry extending to the Southeast quadrant of Mrk\thinspace1239.  Detailed analysis is limited by the short exposure time, and deeper observations would be needed to determine the exact location of the soft emission.
\begin{figure}
	\centering
	\includegraphics[width=\linewidth]{aug10ChandraV9-105x185.png}
	\caption{Snap-shot image of Mrk\thinspace1239\ taken with {\it Chandra}-ACIS-S instrument.The JVLA 5.2 GHz radio contours are provided by  \protect\cite{Jarvela+2022}.   }
	\label{fig:june7_softChandra}
\end{figure}
\subsection{Emission Line Search}
\begin{figure*}
	\centering
	\includegraphics[width=\linewidth]{oct19_dataV1_mnras}
	\caption{\textit{Top panel:} the RGS spectrum is shown in red, while the background model is shown in black. \textit{Middle panel:} the ratio of the background model to the background data. \textit{Bottom panel:} The ratio of the source data to the background model is shown in red. In blue is the maximum fit improvement (Max $\dec$) found by the line search (right axis). The black line is at a ratio of 1  on the left axis and at a maximum fit improvement of 0 C-stat on the right axis. The blue line details where the source is above the background.}
	\label{fig:aug11_dataV1_mnras}
\end{figure*}
Mrk\thinspace1239\ is not a particularly bright source, as such we cannot confidently identify the continuum level. Instead,  we observe  emission lines poking above the background (Figure \ref{fig:aug11_dataV1_mnras} top panel). Due to the high background level it was imperative to find all regions where the data were significant over the background. This was achieved by modelling the background (see Figure \ref{fig:aug11_dataV1_mnras} middle panel) and applying the model to the source spectrum.

A Gaussian profile with a width of $\sigma=0.01$\,\AA\ was then applied beginning at 7\,\AA, The normalization was stepped through twenty values space evenly in log space between $5 \times 10^{-7}$ and $5 \times 10^{-5}$ \anorm\ using the {\sc xspec} command {\sc steppar}. For each normalization step the best fit improvement was recorded. When complete the Gaussian profile was moved 0.005\,\AA\ (an order of magnitude finer than  the instrument resolution of RGS ) and the process was repeated, covering all wavelengths between 7 and 35\,\AA.  This allowed for the spectra to be oversampled and provide a clear picture of where the source was detected above the background.

The maximum fit improvement was taken at each wavelength regardless of the normalisation required.  The maximum fit improvement (Max $\dec$) verse wavelength was plotted as a visualization tool to more clearly identify the locations of possible lines. This max fit improvement line can be seen in the bottom panel of Figure \ref{fig:aug11_dataV1_mnras}.  Forty-five line candidates were identified, all having a fit improvement of at least $\dec=5$.

Next a Gaussian profile  with a width of $0.01$\,\AA\ was applied for each of the 45 line candidates, the normalizations were fit and subsequently frozen. The error on the line energies was then calculated using the {\sc xspec} error command. The uncertainty in the line energy, plus an outflow correction factor (see Sec. \ref{sec:out}), was used as the search range up until a maximum search range of 0.1\,\AA. With the line wavelength and search range  imputed in AtomDB, the most probable line candidate was identified.

Forty line candidates were recovered, they are detailed in Table \ref{tab:BLS}, while the data and the search line are shown in Figure \ref{fig:april18figure1}. No radiative recombination continuum (RRC) features were identified through this process as a continuum was not confidently identified. The search was preformed over the wavelength range $7-35$\,\AA, however lines were only significantly detected between $9-26$\,\AA.
\renewcommand{\arraystretch}{1.25}
\begin{figure*}
	\centering
	\includegraphics[width=\linewidth]{aug11_lnsrchV9_mnras.png}
	\caption{Results of the blind line search with candidate lines labelled. Red denotes the ratio of the RGS data to the background model. The blue line is the maximum fit improvement (Max $\dec$) at any wavelength. The shaded green and gray regions denote the search radius used in AtomDB and the label given in each shaded region is the most probable line ID that produced that emission feature. The different colour bands are simply used to distinguish nearby lines. As a reminder, anywhere the RGS data drops below 1 is a region where the background was above the source.}
	\label{fig:april18figure1}
\end{figure*}
\begin{table}
	\centering
	\begin{tabular}{r c c c c c }
\hline
\multicolumn{1}{c}{(1)} 			& (2) 						& (3) 	& (4)       &(5)       	\\
\multicolumn{1}{c}{ }				& Observed					&Search &  Rest     &	       	\\
\multicolumn{1}{c}{Ion ID}	  		& wavelength	  			& range	&wavelength &$\Delta$C 	\\
									&[\AA]		  				&[\AA]  &[\AA]		& 			\\
\hline
		Fe\,{\sc xx} 	 			& $9.05 \pm 0.04$			& 0.059	&	9.066	&	17	  	\\
		Mg\,{\sc xi} 	 			& $9.2\pm 0.2 $  	 		& 0.100	&	9.169	&	10		\\
		Fe\,{\sc xxi} 	  			& $9.60\pm 0.08$	      	& 0.100	&	9.504	&	7		\\
		Fe\,{\sc xxi} 				& $9.8\pm 0.2 $		  		& 0.100	&	9.819	&	11		\\
		Fe\,{\sc xix} 	 			& $9.92^{+0.01}_{-0.02}$ 	& 0.038	&	9.945	&	30		\\
		Ne\,{\sc x} 	 			& $10.22 \pm 0.04$			& 0.057	&	10.238	&	16		\\
		Fe\,{\sc xviii}$^{\star}$	& $10.32^{+0.06}_{-0.05}$ 	& 0.072	&	10.360	&	15		\\
		Fe\,{\sc xvii} 	  			& $10.47^{+0.03}_{-0.04}$ 	& 0.053	&	10.504	&	20		\\
		Fe\,{\sc xxiv} 	  			& $10.7^{+0.1}_{-0.2}$		& 0.100	&	10.663	&	8		\\
		Fe\,{\sc xxiii}	  			& $10.96^{+0.06}_{-0.04}$ 	& 0.069	&	10.981	&	12		\\
		Fe\,{\sc xvii} 	 			& $11.26^{+0.04}_{-0.13}$ 	& 0.100	&	11.254	&	12		\\
		Fe\,{\sc xxiv} 	 			& $11.49\pm 0.05$			& 0.071	&	11.432	&	18		\\
		Fe\,{\sc xxii} 	 			& $12.0\pm 0.2$		  		& 0.100	&	11.977	&	9		\\
		Ne\,{\sc x} ($\alpha$)  	& $12.14^{+0.03}_{-0.04}$ 	& 0.059	&	12.132	&	26		\\
		Fe\,{\sc xxi} 				& $12.34^{+0.06}_{-0.04}$ 	& 0.072	&	12.284	&	22		\\
		Fe\,{\sc xx} 	 			& $12.531^{+0.1}_{-0.2}$ 	& 0.100	&	12.576	&	14		\\
		Fe\,{\sc xx} 				& $13.1^{+0.6}_{-0.3}$		& 0.100	&	12.864	&	9		\\
		Ne\,{\sc ix} ($r$)  		& $13.47^{0.05}_{-0.03}$ 	& 0.067	&	13.447	&	25		\\
		Ne\,{\sc ix} ($f$)  		& $13.68^{+0.04}_{-0.03}$ 	& 0.061	&	13.699	&	37		\\
		Fe\,{\sc xxi} 	  			& $13.98^{+0.03}_{-0.02}$ 	& 0.053	&	14.008	&	34		\\
		Fe\,{\sc xviii} 			& $14.17^{+0.03}_{-0.02}$ 	& 0.055	&	14.208	&	45		\\
		Fe\,{\sc xviii}${\dag}$  	& $14.31^{+0.03}_{-0.12}$ 	& 0.100	&	14.373	&	35		\\
		O\,{\sc viii} 				& $14.82^{+0.05}_{-0.03}$ 	& 0.070	&	14.821	&	17		\\
		Fe\,{\sc xvii} 				& $15.01^{+0.02}_{-0.04}$ 	& 0.060	&	15.014	&	37		\\
		Fe\,{\sc xix} 			  	& $15.13^{+0.06}_{-0.04}$ 	& 0.081	&	15.079	&	33		\\
		Fe\,{\sc xvii}				& $15.44^{+0.04}_{-0.08}$ 	& 0.092	&	15.453	&	18		\\
		Fe\,{\sc xviii}	 			& $15.93^{+0.02}_{-0.03}$ 	& 0.057	&	15.931	&	25		\\
		Fe\,{\sc xviii}	 			& $16.05^{+0.04}_{-0.2}$ 	& 0.100	&	16.071	&	15		\\
		Fe\,{\sc xvi}				& $16.60^{+0.05}_{-0.03}$ 	& 0.072	&	16.567	&	29		\\
		Fe\,{\sc xvii}			 	& $16.70^{+0.04}_{-0.09}$ 	& 0.099	&	16.780	&	32		\\
		Fe\,{\sc xvii}				& $17.03\pm 0.02$			& 0.052	&	17.051	&	64		\\
		O\,{\sc vii} 				& $18.65^{+0.2}_{-0.03}$ 	& 0.100	&	18.627	&	5		\\
		O\,{\sc viii} ($\alpha$) 	& $18.93\pm 0.02$			& 0.054	&	18.967	&	73		\\
		N\,{\sc vii} 				& $19.88^{+0.05}_{-0.06}$ 	& 0.094	&	19.826	&	6		\\
	  	O\,{\sc vii} ($r$)			& $21.55^{+0.04}_{-0.03}$ 	& 0.079	&	21.602	&	6		\\
	  	O\,{\sc vii} ($i$)			& $21.82^{+0.06}_{-0.04}$ 	& 0.096	&	21.804	&	6		\\
	  	O\,{\sc vii} ($f$)			& $22.06^{+0.01}_{-0.02}$ 	& 0.059	&	22.098	&	56		\\
		Fe\,{\sc xxiv} 				& $22.3^{+0.1}_{-0.2}$		& 0.100	&	22.249	&	6		\\
		N\,{\sc vii} 				& $24.74^{+0.02}_{-0.01}$ 	& 0.067	&	24.779	&	97		\\
		N\,{\sc vii} 				& $24.94^{+0.02}_{-0.04}$ 	& 0.078	&	24.898	&	34		\\
\hline
	\end{tabular}
	\label{tab:BLS}
	\caption{ (1) the emission Line ID as determined by the highest probability transition within the search range given by AtomDB. (2) The wavelength detected by the Gaussian line profile, given in Angstroms. (3) The search range used to identify the lines, it was determined by the size of the uncertainty in wavelength placement plus a z =0.002 inflow/outflow allowance up to a maximum of  0.1\,\AA. (4) The rest wavelength as reported from AtomDB. (5) The fit improvement of including this line. ($\star$) This line was fit independently of all others due to its close proximity of the Ne\,{\sc x} line at 10.238\,\AA. ($\dagger$) This line was initially frozen and then fit will all other lines free}
\end{table}
\subsection{\ovii\ Outflow velocity}
\label{sec:out}
\begin{figure}
	\centering
	\includegraphics[width=1\linewidth]{jan5_fig3.png}
	\caption{The velocity profile for the forbidden line of \ovii\ with a Gaussian line overlaid. The forbidden \ovii\ line is the strongest single feature in the RGS spectra, and thus is the best line candidate to explore the  outflow velocity. Based on this line the outflow velocity is  $-500^{+100}_{-300}$ km~s$^{-1}$.  }
	\label{fig:july28oviiprofilr}
\end{figure}
While examining possible line candidates it was found that several emission features exhibited slight blue shifts in their line energy. Figure \ref{fig:july28oviiprofilr} shows the velocity profile of  the forbidden \ovii\ line. This was the strongest isolated line. The rest frame wavelength of this line is $22.06^{+0.01}_{-0.02}$\,\AA\ while the laboratory measurement is  22.098\,\AA. The resultant outflow velocity is $-500^{+100}_{-300}$ km~s$^{-1}$. The measured velocity  is comparable to other outflows in Seyfert galaxies (see NGC\;5548 \citealt{Mao+2018}; NGC\;4151 \citealt{Armentrout+2007}; NGC\;1068 \citealt{Grafton+2021}). Outflow velocities from other emission lines were not as well constrained.  As well,  Mrk\thinspace1239\ does show evidence for an outflow in the \oiii\ line of $\sim2000$ km~s$^{-1}$, located between 10-100 pc of the SMBH \citep{Pan+2021}.
\subsection{Plasma Diagnostics}
He-like ions produce very specific emission features dependant on the density, temperature and processes by which the atoms were ionized. A common method for determining the density and temperature of the emitting plasma is by using the $G$ and $R$ ratios defined in Equations \ref{eq:R} and \ref{eq:G}. For these diagnostics $z$ is the  forbidden line ($f$), $x+y$ are the inter-combination lines ($i$), and $w$ is the resonance line ($r$) \citep{Porquet+2000}.
\begin{equation}
\centering
R = \frac{z}{x+y} = \frac{f}{i}
\label{eq:R}
\end{equation}
\begin{equation}
G  = \frac{x+y+z}{w} = \frac{f+i}{r}
\label{eq:G}
\end{equation}
For Mrk\thinspace1239\ only  the He-like triplet of  \ovii\ is strongly detected. There is some evidence of Mg\,{\sc xi} and \neix, but the \ovii\  triplet is clearly seen. This may be due to the blend of lines around Mg and Ne originating from other species and the relative lack of other lines around the \ovii\  triplet. For this work both the \ovii\ and \neix\ He-like triplets are explored.

Due to the weak intercombination lines in this spectra it was difficult to ascertain the uncertainties of the $G$ and $R$ ratios. To overcome this we implemented a model in {\sc xspec} so that the values and their uncertainties might be directly outputted. This was achieved by re-writing the $G$ and $R$ ratio such that $f$ and $r$ could be expressed in terms of $G$, $R$ and $i$. To do this the first step is a simple rearrangement of the $R$ ratio:
\begin{equation}
f = R\times i \;.
\label{Eq:f}
\end{equation}
Next the $G$ ratio:
\begin{equation}
r  = \frac{f+i}{G} \;.
\end{equation}
Substituting in what we know about $f$ from Equation \ref{Eq:f};
\begin{equation}
r  = \frac{R\times i+i}{G}
\end{equation}
The subtlety of this method is revealed when we consider how to implement it directly in {\sc xspec}. The expression would be: $(constant_G\times zagauss_r + zagauss_i + constant_R \times zagauss_f)$, where the subscripts denote the parameter the component represents. The correct normalization for $zagauss_r$ is $\frac{R\times i+i}{G}$ but the component is multiplied by $G$ therefore the normalization must be $\frac{R\times i+i}{G^2}$. Similarly the correct normalization on f is  $R\times i$, because the component is already multiplied by $R$ we simply link the two components. This leaves $G$, $R$ and $i$ free to fit the triplet.  The added complication here is the outflow we see on \ovii($f$), thus we have an addition component, {\sc vmshift}, to account for the outflow. The outflow velocities are linked for the $r$, $i$ and  $f$ lines. The $G$ and $R$ ratio for \ovii\ and \neix\ as well as their outflow velocities are displayed in Table \ref{tab:GR}.
\begin{table*}
	\begin{tabular}{c c c c c}
		\hline
		(1)								  & (2)									 &     (3)		&   (4)				&	(5)					\\
		Model Component   & Model Parameters	& Short name	& \ovii\ 			& \neix\			\\
		\hline
		Resonance Line  	 & Wavelength 				   & $\lambda_r$	& 21.602\,\AA		& 13.447\,\AA	\\
											 & Normalization	 	& $r$			& $=(i+i*R)/G^2$	& $=(i+i*R)/G^2$ \\
		\hline
		Intercombination Line	& Wavelength 		& $\lambda_i$	& 21.804\,\AA 		& 13.553\,\AA	 \\
											& Normalization	 	& $i$			& $10^{-5}$			& $3\times10^{-6}$\\
		\hline
		Forbidden Line		 & Wavelength  		& $\lambda_f$	& 22.098\,\AA		& 13.699\,\AA 	 \\
											& Normalization	 	& $f$			& $=i$				&$=i$					\\
		\hline
		G 								  &			 		& $G$			&$5^{+15}_{-3}$		&$1.5^{+3.5}_{-0.9}$			\\
		{\sc constant} 		  &					&				&					&						\\
		\hline
		R 								  & 					& 	$R$			&	$>2.5$			&	$>1$				\\
		{\sc constant} 																							\\
		\hline
		Outflow Velocity    & 					& 				&	$-500\pm200$	km~s$^{-1}$	&	$-300^{+800}_{-600}$ km~s$^{-1}$				\\
		{\sc vmshift} 																							\\
		\hline
	\end{tabular}
\label{tab:GR}
\caption{G and R diagnostic ratios for \ovii\ and \neix. The model components are listed in column (1), the model parameters are listed in column (2) and their short form names used in the text are listed in column (3). The values for each parameter for the \ovii\ and \neix\ triplets are given in columns (4) and (5), respectively.  }
\end{table*}

Comparing the strength of the resonance line to the strength of the Ly${\alpha}$ line reveals that the ratio of H-like/He-like atoms for both the \ovii\ and \neix\ atoms is of order unity. This allows us to use the G ratio to find a temperature of the order $10^6$ K \citep{Porquet+2000}. This estimation is not with out its caveats, as it assumes that the regions where the H-like ions are formed is the same as where the He-like ions are formed which is not necessarily the case according to \cite{Porquet+2000}.

Due to the weak intercombination line in both triplets the R ratio was not well constrained and we can only provide lower limits on it and the density estimations. Comparing the calculated values to the $R(n_e)$ curves produced by \cite{Porquet+2000} we gain limited information. For \ovii, which has $R>2.5$, this tells us that the plasma that produced these emission lines has a density less than $10^{11}$\cm$^{-3}$. The plasma that produced the \neix\ emission lines has a density less than  $10^{12}$\cm$^{-3}$.  Moreover, the outflow velocity of the \ovii\ triplet is consistent with the velocities reported in Section \ref{sec:out}. The outflow velocity of the \neix\ triplet  agrees with that of the  \ovii\ triplet, and it also agrees with being at rest with the galaxy as a whole.

The finial piece of information we can deduce from the He-like ions is from the relative strength of the resonance line.  In purely photoionized plasmas the resonance line is weaker compared to the intercombination and forbidden lines. Whereas, in mixed plasmas, containing both photoionized and collisionally ionized emission, the resonance line is relatively strong. If we look at the \ovii\ triplet we see an incredibly strong forbidden line and weak intercombination and resonance lines. Assuming we are in the region where $R$ is very sensitive to density (which was allowed by the calculated $R$ ratio) than this is indicative to a photoionized plasma.

\subsection{Neon Abundance}
\label{sec:noNeon}
A neon overabundance was proposed in \cite{Grupe+2004} to explain a feature they observed at 0.91 $\kev$ in the 5\ks\ PN data. The feature was attributed to the \neix\ triplet. They did not detect the \ovii\ triplet located at $\sim$ 0.56 \kev\ leading them to conclude that there was an overabundance of Ne in Mrk\thinspace1239. Subsequent works have used this overabundance while fitting the X-ray spectra (see \citealt{Jiang+2021}).

The analysis of \cite{Buhariwalla+2020} suggested the excess emission at 0.91 \kev\ was in fact due to a blend of emission lines including the \neix\ triplet and Fe-L emission lines. In Figure \ref{fig:Neon}  a comparison between the RGS spectra and the EPIC PN spectra is made. Here we see both the RGS (in red) and PN spectra (in black) compared from $11-16$\,\AA. Both spectra have been normalized by their mean in this region to ease the comparison.

The red RGS data points show a plethora of emission features while the black PN data show an excess emission between 0.85 and 1 \kev, similar to what was seen in the original 2001 observation by \cite{Grupe+2004}. Even with the deep 2021 PN observation, the \neix\ emission lines remain unresolved. The RGS spectrum is required to resolve the \neix\ triplet and the surrounding features to  state that the excess emission seen in this object at 0.91 \kev\ is not due solely to the \neix\ triplet but instead a it is a blend of many emission features in this region.

Furthermore the overabundance described by \cite{Grupe+2004} was reliant on the non-detection of the \ovii\ triplet in addition to the enhancement of the \neix\ line by Fe-L emission. In the bottom panel of Figure \ref{fig:april18figure1} we can clearly see the \ovii\ present in our RGS spectra. The detection of the \ovii\ triplet and the confirmation that the \neix\ triplet was being enhanced by surrounding emission lines the shows a Ne overabundance is not needed in this source.
\begin{figure}
	\centering
	\includegraphics[width=\linewidth]{oct19_ne_rgs_vs_pn}
	\caption{Comparison of the spectral resolution of the EPIC PN instrument to the RGS instrument. Red shows the RGS data, blue is the maximum fit improvement line (same as Figure \ref{fig:aug11_dataV1_mnras}) black denotes the PN data and green shows the location of the \neix\ triplet. The PN and RGS spectra have been normalized so that they may be presented on the same scale.}
	\label{fig:Neon}
\end{figure}
\section{Spectral Modelling}
\label{sec:spectral}
The RGS data were background modelled in {\sc xspec}, and C-statistics were used \citep{Cash+1979} to evaluate the fit quality throughout. Errors were calculated at the 90 \% confidence level using the {\sc xspec} error command. Based on the results of the plasma diagnostics we proceeded with the assumption that there were at least two ionized regions present in the spectra.

We tested collisionally ionized plasma with  {\sc apec}, version 3.0.9, \citep{Smith+2001} and a photoionized plasma using  {\sc xstar} grids, version 2.31 \citep{Kallman+2001}. The model  {\sc mekal}\ \citep{Liedahl+1995} was tested in this analysis as it was what was originally used to model the data in \cite{Buhariwalla+2020}. Ultimately {\sc apec}\ was used as it contains more up to date atomic data, and produced better fit statistics. All fits were absorbed  with a Galactic  column density of {$N_{\rm H}$}$=4.43\times10^{20}\; \text{atoms/cm}^2$ \citep{Willingale+2013}.

{\sc xstar} grids were calculated with a covering fraction of 1, luminosity of $10^{42}$\ergps\ and the turbulent velocity was fixed at 300 km~s$^{-1}$. The ionizing continuum was described by a powerlaw with $\Gamma=2.2$.  The column density was logarithmically sampled at five points between $10^{20}$ and $10^{24}$ cm$^{-2}$. The density of the emitting material was logarithmically sampled at six points between $10^9$ and $10^{14}$ cm$^{-3}$  The ionization parameter is defined as $\xi = L/nr^2$ \citep{Tarter+1969, Kallman+2001} where $L$ is the ionizing luminosity, $n$ is the density of the material and $r$ is the distance between the ionizing source and the emitting material. The log ionization was sampled linearly at 6 points between $\text{log} \xi$ 1 and 4 erg cm s$^{-1}$. The column density, ionization and density were all interpolated for a total of a 300 step grid.

Seven subsequent grids were created with density at a fixed value ranging from $10^4$ to $10^{10}$ cm$^{-3}$, \, log$\xi$ was linearly samples at 10 points between 0 and 5, the column density was logarithmically sampled at ten points between $10^{20}$ and $10^{24}$ cm$^{-2}$. All other parameters were kept the same. These were 100 step grids. The $10^9$ and $10^{10}$ cm$^{-3}$ grids were tested against the 300 step grid for consistency, they produced the same fit statistic and ionization parameters. All abundances were fixed at solar.

For  {\sc apec}\  the plasma temperature ($kT$) and normalization remained free, while the abundances remained fixed at solar. The powerlaw, when included, had free photon index and normalization.

Initially, emission from only one plasma component ({\sc apec}\ or  {\sc xstar})  was fitted to the RGS spectrum between $7-31$\,\AA. With each individual plasma, the addition of a powerlaw component was considered to account for the underlying continuum.
Fitting each of these components individually was able to account for select emission lines, but not all the emission features. The fits were $\cdof=1319/814$ for {\sc apec}\ and $\cdof=1360/813$ for {\sc xstar}. The addition of a powerlaw ({\sc po}) component had fit improvements of   $\dec=22$ for  {\sc apec}+{\sc po} for two additional free parameters. The addition of the powerlaw to {\sc xstar}\ had a drastic fit improvement of $\dec=83$ for two additional free parameters, but this was due to the poor initial fit of the {\sc xstar}\ model {\sc xstar}+{\sc po} results in a fit of $\cdof=1277/811$.

Next two component plasmas were tested where both components were of the same type, this was done to explore the possibility of a single plasma with a gradient of temperature and/or ionization values. This produced slightly better fits than the addition of the powerlaw, with the  fit improvement for {\sc apec}+{\sc apec}\  was $\dec=27$ the single component counterpart, for two additional free parameters. {\sc xstar}+{\sc xstar}\ resulted in a fit improvement of $\dec=69$, producing a worse fit than a single component+powerlaw model. \cite{Jiang+2021} fit the the soft band of Mrk\thinspace1239\ with two collisionally ionized components. However their CCD spectra lacked the data quality to  distinguish between a collisionally ionized and photoionized component.

\subsection{Testing for Two Distinct Plasmas }
\label{subsec:two_plasma}
A mixed plasma was tested such that it contained a collisionally ionized ({\sc apec}) and a photoionized  ({\sc xstar}) component.  This produced a fit of  \cdof=1201/811. The {\sc apec}\ component produces a $\dec \sim15$ better than the {\sc mekal}\ component for the same degrees of freedom. The \ovii$(f)$ line, which was outflowing with a velocity of $\sim-500$ km~s$^{-1}$, was produced by the photoionized plasma. Thus we test the {\sc xstar}+{\sc apec}\ models with the photoionized component allowed to outflow.  The {\sc vmshift}\,$\times$\,{\sc xstar}\,+\,{\sc apec}\ model produces a superior fit of \cdof=1156/810, $\dec=45$ for one additional free parameter. Similar fit improvements were seen in the {\sc xstar}+{\sc xstar}\ and {\sc apec}+{\sc apec}\ fits when one component is allowed to outflow. No significant fit improvement were seen in these fits when both components are allowed to outflow. The best fit model can be found in Figure \ref{fig:aug8fig3},  the best fit parameters can be found in Table \ref{tab:xstar}.
\begin{figure}
	\centering
	\includegraphics[width=1\linewidth]{Aug8_fig3.png}
	\caption{\textit{Top panel:} the ratio of the RGS data and background model is shown in red, the fit here is \cdof= 1716/816. All data above ratio$=1$  is above the background level. The best fit two component plasma  is shown in green ({\sc apec}) and blue ({\sc xstar}). The individual models are plotted to show which features are fit by each. \textit{Bottom panel:} the red data are the same as above. The black data are the ratio of red data with the best fit model folded through the response (\cdof=1158/810).  }
	\label{fig:aug8fig3}
\end{figure}

\begin{table}
	\centering
	\label{tab:xstar}
	\begin{tabular}{c c c }
		\hline
		(1)		&(2)								&(3)					\\
Model Component	& Model Parameter					& Value 				\\
		\hline
		{\sc apec} 	& kT [keV]										& $0.74_{-0.07}^{+0.06}$\\
		 		& Abundance 										& $1^f$					\\
		 		& Normalization  [$\times10^{-5}$]  & 		$3.8\pm 0.6$			\\
		{\sc xstar} 	& Column density [cm$^{-2}$]& $<\num{2.5e22}$	  	\\
		 		& Log Ionization [erg cm s$^{-1}$]	& $1.7_{-0.2}^{+0.3}$ 	\\
		 		& Density 	[cm$^{-3}$]						& $>\num{3.2e10}$     	\\
		 		& Redshift												& $0.01993^f$		  	\\
				& Normalization	[$\times10^{-6}$]	 & 		$2.6_{-0.4}^{+0.5}$ \\
  {\sc vmshift} & Velocity	[km/s]					& 	$-660\pm150$	  	\\
		\hline
	\end{tabular}
	\caption{Best fit model parameters for the two component plasma model. Column (1) indicates the model component,  Column (2) gives the Model Parameter and  Column (3) gives the value of each parameter. All parameters with the superscript `$f$' are kept fixed at quoted values. }
\end{table}
We can see that the photoionized material is outflowing with a velocity of $-660\pm150$ km~s$^{-1}$\ which is consistent with the outflow velocity that we saw with \ovii($f$) emission line. The density of the emitting material is $>\num{3.2e10}$ cm$^{-3}$, which is consistent with the density estimation found using the \ovii\ $R$ ratio. The temperature of the emitting plasma is consistent with the temperature found in \cite{Buhariwalla+2020}, despite the difference in model component ({\sc apec}\ vs {\sc mekal}) and the difference in instrument (RGS vs MOS). There appears to be an excess of emission below 12\,\AA\ and the $\sim25$\,\AA\ feature attributed to N{\sc vii} appears under fit.

The best fit two component plasma gave an upper limit on density of less than $\num{3.2e10}$ cm$^{-3}$, the {\sc xstar}\ grid that gave this result was limited to densities above $10^9$ cm$^{-3}$. Thus the lower limit of density was estimated using previously made grids with the same parameters.   Figure \ref{fig:sept29_Figure_2} (left axis) shows the C-stat as a function of density, the minimum C-stat is 1156, given by the best fit model. 1, 2, and 3 $\sigma$ lines are drawn to show the allowed density of the photoionized plasma.  The allowed density is then between $10^5 - 10^{10}$ cm$^{-3}$. Ionization as a function of density (Figure \ref{fig:sept29_Figure_2} right axis) will be discussed in section \ref{subsec:PIP}.
\begin{figure}
	\centering
	\includegraphics[width=1\linewidth]{sept29_Figure_2.png}
	\caption{ Density as a function of C-Stat (left axis) for the best fit two component plasma at different densities, and measured ionization at different density values (right axis).}
	\label{fig:sept29_Figure_2}
\end{figure}

{\sc xtar}  grids with ionizing continuum described by a powerlaw with  $\Gamma=1.7$ and $\Gamma=2.5$ were generated to test the effects of  the continuum shape on the spectral appearance. They produced similar fit statistics as the best fit model, with $\Gamma=2.2$. The measured properties of the {\sc xtar}  and {\sc apec} components  were consistent, within error, to those quoted in Table \ref{tab:xstar}.  This led us to conclude that the shape of the continuum  has a minimal effect on the properties of the photoionized emitter. This does not seem unreasonable in Mrk\thinspace1239\ due to the modest data quality and small number of emission lines fit by the {\sc xstar } component.

{\sc xstar}\ grids with variable abundances were generated to explore the possibility of an overabundance in Mrk\thinspace1239. To ease computational time the density was fixed at $10^6$ cm$^{-3}$ and column density was fixed at $10^{23}$ cm$^{-2}$. All other parameters remained the same. Grids with variable Ne, O, Fe, and N were generated. Variable Ne and Fe  had no affect on the fit of the RGS spectrum.  Variable O and N produced best fit models with super solar ratios of \NOx. The best fit grid produced a relative abundance of $\NOx =4\pm2$ $A_{\odot}$, with a $\dec=16$. After investigation it was found that the driving mechanism for this super solar abundance is the 25\,\AA\ N{\sc vii} feature that was under fit in Figure \ref{fig:aug8fig3}. When the variable N grids are fit to the data below 25\,\AA\ no overabundance or fit improvement is found.   However  when the data above 25\,\AA\ was included the overabundance is present. \cite{Mao+2019} showed that the relative \NOx\ abundance could be super-solar in the ICM due to enrichment from low and intermediate mass stars. Their measured \NOx\ agree with the value found here. However it is not immediately clear if this enrichment effect  could also be present in the photoionized material of  Mrk\thinspace1239, or if some other processes is driving the excess emission at 25\,\AA.

The model under fits the data below 12\,\AA,  and no combination of collisionally or photoionized plasma could sufficiently fit the data in this region.  This excess emission is also seen in the PN spectrum from this observation (Buhariwalla {et al.\ }\ in prep.) indicating that the excess emission is probably not from statistical fluctuations. Instead it may be that the AGN continuum in  Mrk\thinspace1239\ is becoming visible at wavelengths less then 12\,\AA\ (energies greater than 1\;\kev). The RGS data are insufficient to make any inferences about the continuum. Thus this feature will be followed up in the future work on Mrk\thinspace1239\ and its broadband spectra.
\section{Discussion}
\label{sec:disc}
\subsection{The Origin of the Photoionized Plasma}
\label{subsec:PIP}
\begin{figure}
	\centering
	\includegraphics[width=1\linewidth]{jan6_fig1.png}
	\caption{Density as a function of distance outlying the density-distance combinations that would produce the correct ionization for a given luminosity. The boxes outline the density-distance combinations for known components in  AGN,  they include the narrow line region (NLR), the broad line region (BLR) and the torus. See Section \ref{subsec:PIP} for further details.  }
	\label{fig:sept29_Figure_1}
\end{figure}
To explore the origin of the photoionized emitter we examined properties that could result in the measured ionization parameter. Recall the  definition of the ionization parameter $\xi = L/nr^2$ (see Sec. \ref{sec:spectral}). This  can be rearranged as $r(n) = \sqrt{L/n \xi}$, allowing us to plot distance as a function of density. The ionization parameter ($\xi$) is degenerate with the density of the emitting material. This degeneracy is mitigated with the use of multiple {\sc xstar}\ grids generated at different plasma densities. The right panel (right axis) of Figure \ref{fig:sept29_Figure_1} shows $\text{log} \xi$ as a function of density.

The ionizing luminosity ($L$) was estimated using the models in \cite{Buhariwalla+2020}, however due to the uncertainty in the location of the emission material it is undetermined if the photoionized emitter sees the absorbed or unabsorbed luminosity, thus both luminosities are tested. Although the shape of the ionizing continuum would change between the absorbed and unabsorbed scenarios we showed in Section \ref{sec:spectral} where dramatic changes in the continuum have a minimal effect on the properties of the photoionized emitter. Thus for this exercise both luminosities can be used to estimate the properties of the photoionized emitter. We caution that the solution for the unabsorbed luminosity is to be taken as  approximate, as we did not re-run the models using the absorbed continuum to illuminate the emitting gas at all densities. A single {\sc xstar}\ grid with density $10^{10}$ cm$^{-3}$ was constructed using an absorbed ionizing continuum based on the absorbed continuum from \cite{Buhariwalla+2020}. It produced an ionization parameter greater than the unabsorbed scenario, however the difference was small on the distance-density plot. It resulted in a slight shift downward of the allowed regions for the absorbed luminosity.

The blue band in Figure \ref{fig:sept29_Figure_1} (left panel) shows where allowed combinations of distance and density lie for an unabsorbed source. The  red band shows the same for an absorbed source. The variable $\xi$ has some influence on the $r(n)$ plot, causing a small wave-like effect on the bands in Figure \ref{fig:sept29_Figure_1}.  To illustrate what kind of material would be able to produce these emission lines density and distance combinations of known AGN components are drawn in this plane.

The outer radius of the NLR was estimated using   a $\textrm{R}_{\textrm{NLR}} - \textrm{L}_{\textrm{\oiii}}$ relationship derived for a sample of quasars and Seyfert galaxies \citep{Bennert+2002}.  The value of $\sim3$ kpc is consistent with the radius of the extended NLR in Mrk\thinspace1239\ found by \cite{Husemann+2022}. We take the inner radius of the NLR to be 30 pc, obtained from multi-density photoionization modelling of the NLR of Seyfert galaxies that reproduces a large range of NLR emission line intensities \citep{Komossa+1997}. This is consistent with a sample of Seyfert 2s and intermediate Seyfert ($\sim$\ Sy1.5) galaxies \citep{Vaona+2012}, and estimates using a second $\textrm{R}_{\textrm{NLR}} - \textrm{L}_{\textrm{\oiii}}$ relationship for the effective radius of the NLR in Seyfert galaxies given by \cite{Schmitt+2003}. For Mrk\thinspace1239, comparable values of L$_{\textrm{\oiii}}$ are  given by \cite{Rafanelli+1984} and \cite{Malkan+2017}.  The density is estimated measurements of intermediate Seyfert, and  Seyfert 2 galaxies \citep{Bennert+2006, Vaona+2012}.  We note that the density and inner radius of the NLR may extend to lower values, but this provides a good order-of-magnitude estimate for our work.

The BLR size was estimated using R$_{\textrm{BLR}} - \lambda \textrm{F}_{\lambda}(5100 \, \mbox{\normalfont\AA})$ \citep{Kaspi+2005}, with $\lambda \textrm{F}_{\lambda}(5100\, \mbox{\normalfont\AA})$ for Mrk\thinspace1239\ taken from \cite{Grupe+2004} and \cite{Pan+2021}. The density of the BLR  is given by \cite{Netzer+2013, Arav+1998}.

The outer radius of the torus in Mrk\thinspace1239\ is measured to be 3 \pc\ based on 12 $\mu m$ observations \citep{Tristram+2011}. The inner torus radius  (0.1-0.01 pc) is estimated  by the dust sublimation radius with the luminosity estimated using models from \cite{Buhariwalla+2020}. The density of the torus is taken from \cite{Netzer+2013}.

With all these components placed in the density - distance plane we begin to see what kind on material is producing the photoionized emission in Mrk\thinspace1239. A tours like structure produces the correct density - distance combination in both the absorbed and unabsorbed scenario. As discussed in section \ref{subsec:two_plasma}  the allowed density is  between $10^5 - 10^{10}$ cm$^{-3}$, which agrees with a tours like structure again. The farthest allowed emitting region would be dependant on the lowest allowed density ($10^5 \; \textrm{cm}^{-3}$) and it would be at most a few pc, around the measured outer edge of the torus in Mrk\thinspace1239\ \citep{Tristram+2011}.

If we analyses the full width half max (FWHM) of the \ovii\ ($f$) line we fine that FWHM$_{\textrm{\ovii}\;(f)}$ <520 km~s$^{-1}$. The Keplerian velocity of an object orbiting around the SMBH at the inner radius of the torus is $\sim$ 500 km~s$^{-1}$. We can then conclude that the photoionized emitting material cannot be located any closer than the inner radius of the torus.

\begin{figure*}
	\centering
	\includegraphics[width=1\linewidth]{jan5_fig1.png}
	\caption{ \textit{Top panel:} 100\;ks \textit{Athena-XIFU} simulation with the best fit {\sc xstar}+{\sc apec}\ model. The {\sc xstar}\ grid was illuminated by a powerlaw continuum with $\Gamma=2.2$.   The approximate  \textit{XIFU} background is shown in red. \textit{Bottom panel:} the Athena simulation fit with {\sc apec}\ +{\sc xstar}, for four different {\sc xstar}\ grids. The four grids have the same input parameter (covering fraction, luminosity, turbulent velocity, density, column density, ionization and abundances) but are illuminated by different ionizing continuum. The ionizing continuum are a mathematical powerlaw  with $\Gamma=2.2$, $\Gamma=2.5$ and $\Gamma=1.5$, as well as an empirical  continuum model based on the absorbed continuum from the broad band spectra of Mrk\thinspace1239\ \protect\citep{Buhariwalla+2020}. The ratios have been binned up for visual clarity.  }
	\label{fig:dec19_fig1}
\end{figure*}
In summary,  the RGS data have greatly increased our knowledge of the photoionized plasma in Mrk\thinspace1239, allowing us to determine its ionization, density and location. We now know that the material exists  within several parsec of the central engine. Previous to these results it was not abundantly clear that there were emission line features in Mrk\thinspace1239, let alone to the extent detected here. However, several questions about the physical descriptions of the  few parsec surrounding the central engine of Mrk\thinspace1239\ remain. Namely, is this photoionized material related to the absorbing material we know is present in this source? Modelling the broadband X-ray continuum of Mrk\thinspace1239\ requires heavy absorption from neutral and warm absorbers \citep{Buhariwalla+2020}. It is easy to image that a large cloud could be illuminated by the central engine obscuring our line of sight to the central engine, while other regions of the cloud are  photoionized, producing the photoionized spectra we see here.

Broadband analysis will   shed some light onto the physical description of this region. For future X-ray telescope missions this exercise may be trivial. The top panel of Figure \ref{fig:dec19_fig1} shows a simulated \textit{Athena-XIFU}\footnote{\textit{XIFU} response files used were XIFU\_CC\_BASELINECONF\_2018\_10\_10, retrieved from \url{http://x-ifu-resources.irap.omp.eu/PUBLIC/RESPONSES/CC_CONFIGURATION/}} 100\,\ks\ observation \citep{Barret+2022}, using the best fit {\sc xstar}+{\sc apec}\ model. Note that unlike with RGS, the \textit{XIFU} spectrum is several orders of magnitude   above the expected background. Most notably, \textit{XIFU} data will  distinguish between different ionizing continuum for the photoionized {\sc xstar}\ grids. The data were simulated with an ionizing continuum of $\Gamma=2.2$, equivalent to the continuum used for the majority of this work.  The four grids fit to the \textit{XIFU} data were ionized with: $\Gamma=2.2$; $\Gamma=2.5$; $\Gamma=1.5$; and an empirical continuum based on the best fit absorbed continuum  given by  \cite{Buhariwalla+2020}.  The $\dec$ value between the best fit $\Gamma=2.2$ grid and the $\Gamma=2.5$, $\Gamma=1.5$, and absorbed continuum were  $\dec=391$, $\dec=688$ and  $\dec=2476$, respectively.  Thus with \textit{XIFU} we  easily  place restrictions on the ionizing continuum for the photoionized material. This is not even taking into account that \textit{XIFU} will have high resolution spectroscopy up to $\sim$10 \kev\  allowing for direct testing of the link between the absorber and the photoionized emitter.

\subsection{The Collisional Plasma in Mrk 1239}
When testing spectral models no combination of photoionized emitters could fit the $13-17$\,\AA\ region well. Only the {\sc apec}\ and {\sc mekal}\ models could adequately fit the data in this region. Figure \ref{fig:aug8fig3} shows that this is the region of the RGS spectra that is dominated by the {\sc apec}\ component. The luminosity of each component was calculated  using the XSPEC\ command  {\sc clumin}. The {\sc apec}\ component had a luminosity of  $L_{0.3-2.5 \kev}=41.02^{+0.06}_{-0.07}$ \ergps, the {\sc xstar}\ had $L_{0.3-2.5 \kev}=41.05^{+0.06}_{-0.08}$ \ergps. The brightness of each feature was approximately equal.

Emission from  \fexvii\ at  16.777\,\AA\  (3F), 17.050\,\AA\ (3G), and 17.097\,\AA\ (M2)  are commonly seen in collisionally ionized plasma. The 3G and M2 lines are  blended  due to instrument resolution. They are often labelled $I$(17.076\,\AA) as they appear as one emission line feature at that wavelength. The \ovii\ RRC is found at 16.769\,\AA\ which blends with the 3F line much the same way 3G and M2 blend ($I$(16.777\,\AA)). \cite{Phillips+1997} measured the ratios of 3G/3F and M2/3F lines in the spectrum of the solar corona. Using their data the  average value of ratio (3F)/(3G + M2) was calculated to be $0.48 \pm 0.2$. This suggests that if the  ratio of $I$(16.777\,\AA)/$I$(17.076\,\AA) is measured to be above 1 then than 50\% of the flux in the $I$(16.777\,\AA) feature is contributed by the \ovii\ RRC component.

\cite{Bianchi+2010} found that the {\sc apec} model produced a ratio of $\approx 0.44$ at varying plasma temperatures. They then measured the $I$(16.777\,\AA)/$I$(17.076\,\AA) ratio in Mrk 573 to be $1.2_{-0.5}^{+0.7}$, concluding that the \ovii\ RRC feature was contaminated by the 3F emission line.
Similarly \cite{Marinucci+2011} measured the ratio in  NGC 424 to be $2.6_{-1.9}^{+1.8}$, and stated that the RRC component dominated over the \fexvii\ lines. In Mrk\thinspace1239, the ratio is measured at $0.9_{-0.4}^{+0.7}$ indicating that both the \fexvii\ emission and \ovii\ RRC are contributing to $I$(16.777\,\AA).

\cite{Guainazzi+2007} studied obscured AGN (Sy 1.5 or greater) and found  16 of the 69 objects they investigated  showed evidence of the \ovii\ RRC feature. They took a ratio of 0.6 to indicate the RRC feature was dominate of the 3F emission line. Of those 16 galaxies, 15 had a ratio >0.6, and 13 were >1.2.  Mrk\thinspace1239\ ($0.9_{-0.4}^{+0.7}$) fits in the lower end of obscured AGN distribution and would be classified as having a dominate \ovii\ RRC feature over the 3F line by \cite{Guainazzi+2007}.

When the best-fit {\sc xstar} model was applied to the data without the {\sc apec} component and the (3F)/(3G +M2) ratio  was measure  it  resulted in a value <0.5. This indicated that the photoionized emission was accounting for the \ovii\ RRC feature and the remaining flux in this region was due to \fexvii\ emission from collisionally ionized material.

Collisionally ionized emitters have been used to model the X-ray spectra of star forming regions in galaxies.  \cite{Franceschini+2003}  approximated the SFR from the  $L_{2-10 \kev}$ in ULIRGs without dominant AGN components. It is given by:
\begin{equation}
{\rm SFR^{ULIRG}_{X-ray}} \approx \frac{L_{2-10 \kev}}{10^{39} {\rm erg\; s^{-1}}} \Msun {\rm yr^{-1}}\;.
\end{equation}
Using just the {\sc apec}\ component allows us to estimate the non-AGN contribution of the  2-10\kev\ luminosity. We measure the {\sc apec}\ contribution to be $L_{2-10 \kev} \approx \num{3e39}$ \ergps, resulting in a SFR of 3 $\Msun \, {\rm yr^{-1}}$.  This is consistent with  previous measurements SFR in Mrk\thinspace1239\ such as X-ray measurements \citep[$ 3.7-5.8\Msun \, {\rm yr^{-1}}$;][]{Buhariwalla+2020},  PAH measurements \citep[$\lesssim 7.5\Msun \, {\rm yr^{-1}}$;][]{RuschelDutra+2017},    SED fitting \citep[$3.47\pm 0.26\Msun \, {\rm yr^{-1}}$;][]{Gruppioni+2016}, and IR measurements \citep[$2.1^{+0.5}_{-0.4}\Msun \, {\rm yr^{-1}}$;][]{Smirnova+2022} . In the sample by \cite{RuschelDutra+2017}, the luminosity from the collisionally ionized component matched those galaxies whose X-ray spectra  were dominated by starburst component rather than  those dominated by AGN.

\subsection{Mrk\,1239 Compared to Seyfert Galaxies}
In a sample of 21 bright {\it ROSAT}\ selected AGN  viewed with {\it XMM-Newton}-\textit{EPIC}, only 2 of the galaxies studied preferred a soft excess with a line emitting region  over a smooth blackbody  \citep{Gallo+2006C}. The sample included 19 Seyfert galaxies (from Seyfert 1s to Seyfert 1.9s) and 2  LINERs and it was the two LINERs that required the CI component ({\sc mekal}).

The RGS spectrum of Seyfert 1 galaxies is typically dominated by features from  photoionized plasma. This primarily includes emission from He-like ion triplets (e.g. \ovii, \neix, N\,{\sc vi})  and RRC features (\oviii, \neix, C\,{\sc vi}).   NGC 5548 is among the first Seyfert galaxies studied \citep{Seyfert+1943}, it possesses  narrow and broad emission features in both the optical and X-ray bands \citep{Cappi+2016,Whewell+2015}. The RGS spectra of this object shows RRC features from \neix, \ovii, \oviii, C\,{\sc v}, and C\,{\sc vi}. In addition the spectra show a strong \ovii$(f)$ line compared to the $i$ and $r$ lines. This galaxy shows strong evidence of several photoionized emitters and no evidence of any collisionally ionized material \citep{Mao+2018}.

NGC 4151 is another well known Seyfert 1.5 galaxy. It shows many of the same features that NGC 5548 does\ including RRC and He-like triplets. The $R$ and $G$ ratios straddle the boarder between photoionized and collisional ionized plasma. However, the strength and shape of the RRC features in this galaxy point to a purely photoionized plasma \citep{Armentrout+2007}.

The NLS1, Mrk\,335 has been extensively studied with {\it XMM-Newton}. While early RGS spectra with low signal-to-noise supported a CI plasma \citep{Grupe+2008}, theses were not confirmed in later works (e.g. \citealt{Longinotti+2013,Longinotti+2019}). \cite{Parker+2019} examined the low continuum flux state when the emission lines were strongest and attributed them to photoionized emission.

Mrk\thinspace1239\ shows strong evidence of collisionally ionized plasma, putting it at odds with the many of Seyfert 1 observations that favour photoionized emission.  In Mrk\thinspace1239\ the measured  ratio of the forbidden \ovii\ line to the \oviii\ Ly$_{\alpha}$ emission line is $\approx 1$ as is  expected for an obscured AGN (Sy 1.5 or greater) compared to the lower ratio for starburst galaxy \citep{Guainazzi+2007}.

NGC 1068 is a prototypical Seyfert 2 galaxy, and it shows a myriad of emission features stemming from multiple photoionized components. Based on visual inspection this galaxy has a spectrum that looks similar to Mrk\thinspace1239. The He-like triplets in NGC\,1068 have similar appearances, with the  \neix\ triplet showing equally strong $r$ and $f$ lines and no $i$ line matching Mrk\thinspace1239\ completely \citep{Grafton+2021}. However the \neix\ triplet is not as blended with other lines in NGC 1068 as in Mrk\thinspace1239.

Original analysis of NGC\,1068 showed that the maximum strength of a collisionally ionized plasma would be an order of magnitude less than the photoionized plasma \citep{Kinkhabwala+2002}. However, the photoionization models used to describe the RGS spectrum failed to properly account for \fexvii\ at 15\,\AA\ and 17\,\AA, thus \cite{Grafton+2021} used a collisionally ionized component to model these lines. They stressed that this was most likely not caused by starburst activity. Instead it may be due to incomplete photoexcitation information in {\sc spex},  causing {\sc pion}\ to under predict emission in these lines.  As for Mrk\thinspace1239\ this does not seem to be the case as no matter the number of photoionized components applied (in addition to powerlaw components) the spectra is still under fit between $13-17$\,\AA\ without the addition of {\sc apec}.

 Mrk\thinspace1239\ is a polar scattered Seyfert 1 galaxy, meaning our line of sight (LOS) passes through the upper layers of the torus polarizing the optical emission    \citep{Smith+2004,Jiang+2021}. ESO 323-G77, another polar scattered Seyfert1,   presents an X-ray spectra similar  to Mrk\thinspace1239. This object shows high levels of variable absorption, that leave the soft band remarkably consistent below 1.5 \kev\ (Mrk\thinspace1239\ is consistent below 3 \kev) and the harder bands show  typical Seyfert variability (the same in Mrk\thinspace1239). In the soft band of ESO 323-G77 excess emission is seen at 0.9 \kev, and can be modelled using a collisionally ionized plasma. This plasma has been interpreted as a region of star formation activity, and a SFR can be extracted in much the same way as is done here in Mrk\thinspace1239\   \citep{Miniutti+2014}.

ESO 323-G77 has been interpreted to be obscured by a clumpy torus, BLR clouds and a warm outflowing medium  \citep{Miniutti+2014}.
The soft X-ray  band of this object has been modelled using  two {\sc apec}\ components, one with similar temperature to Mrk\thinspace1239\ ($kT=0.74$\kev) , the other much cooler, $kT=0.09$\kev. This is a difference in the modelling of these  two objects, as the  cool {\sc apec}\ component in ESO 323-G77 takes the place of a photoionized emitter. This cool {\sc apec}\ component fits  a slight excess emission around $\sim0.5$ \kev, however it could be argued that the this excess emission could be originating from photoionized plasma  (\ovii\ triplet/N\,{\sc vii}  emission lines), but the CCD spectra lacks the resolution to distinguish between the components. Mrk\thinspace1239\ could have a geometry similar to this object  where  the LOS is through the torus, creating a polarized BLR and absorbing the X-ray continuum, allowing for distant star forming regions to be visible in the  soft X-ray band. ESO 323-G77  may exhibit a similar RGS spectrum as we have seen here in Mrk\thinspace1239.

If the soft X-ray band is Mrk\thinspace1239\ is obscured in a similar manor as it is in ESO 323-G77, then an interesting avenue for analysis would be to investigate the connections with the photoionized emitter producing the \ovii\ triplet and the obscuring torus. Both sources would benefit from much deeper observations to explore this connection.

\section{Conclusions}
\label{sec:conclusion}
In this work we present the first deep RGS spectrum of Mrk\thinspace1239, it  shows a myriad of ionized emission lines originating from a blend of photoionized  and collisionally ionized plasma. The collisionally ionized plasma dominates the RGS spectra below 17\,\AA, while the photoionized plasma dominates above 17\,\AA. The \ovii\ triplet is detected for the very first time in this source. Our main conclusions are as follows:
\begin{itemize}
	\item[i.] Based on a strong 0.91 \kev\ feature interpreted as \neix\ triplet and the lack of \ovii\ triplet,  \cite{Grupe+2004} concluded that there was a super solar abundance of  Ne/O. We show that the \neix\ triplet is enhanced in a CCD spectra  by Fe-L transitions, as proposed by \cite{Buhariwalla+2020}, and that the \ovii$(f)$ triplet is present. All fits done in this work use a solar abundance of Ne/O, allowing us to conclude that Mrk\thinspace1239\ does not have an over abundance of Ne.
	\item[ii.]  Based on the line ratios of the \ovii\ triplet and the variable density {\sc xstar}\ grids, we can infer that the density of the photoionized material is between $10^5-10^{10} \cm^{-3}$. Using the definition of ionization and luminosities  estimated from \cite{Buhariwalla+2020} we can estimate that the photoionized emitter is located within a few \pc\ of the SMBH. Based in on the FWHM of the \ovii$(f)$ emission line, the emitting material cannot be closer than the inner radius of the torus.
	\item[iii.] The collisionally ionized material is interpenetrated as a region of star forming activity. Based on relations from \cite{Franceschini+2003} we estimate the SFR to be $\sim3 \Msun \, {\rm yr^{-1}}$. This is consistent with the the SFR found by \cite{Buhariwalla+2020} using the same method, and using PAH measurement of Mrk\thinspace1239\ \citep{RuschelDutra+2017}.
	\item[iv.] Seyfert 1 galaxies typically show evidence of photoionized emission, it is uncommon for them to show evidence of collisionally ionized plasma in their X-ray spectra, placing Mrk\thinspace1239\ in the minority along with ESO 323-G77.
\end{itemize}
The broadband spectrum Mrk\thinspace1239\ will be presented in a future work. We will investigate any connection between the photoionized emitter and the highly absorbed continuum.
\section*{Acknowledgements}
This work was based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. This research has also made use of data obtained from the Chandra Data Archive and the Chandra Source Catalog, and software provided by the Chandra X-ray Center (CXC) in the application packages CIAO and Sherpa.
We thank the referee for comments that helped clarify the work. LCG acknowledges financial support from the Natural Sciences and Engineering Research Council of Canada (NSERC) and from the Canadian Space Agency (CSA). JJ acknowledges support from the Leverhulme Trust, Isaac Newton Trust, and St. Edmund's College University of Cambridge.

\section*{Data Availability}
All data will be available through the HEASARC  archive after the proprietary period has ended.

\bibliographystyle{mnras}


\section{Introduction}
One of the main challenges in theoretical physics is to find a theory that consistently incorporates quantum and gravitational physics. The attempts to construct such a theory can be grouped into two categories: first, there are theories proposed from first principles that incorporate some quantum aspects into gravitation, or vice versa. This strategy includes theories like string theory, loop quantum gravity, and other remarkable ideas (for a review of most quantum gravity approaches see Ref.~\cite{Oriti}). In contrast, another strategy is to consider general relativity as an effective low-energy description of the fundamental theory and to look for clues on the theory of quantum gravity by looking for small deviations from general relativity. To search for these deviations, it is often useful to have a parametrisation that, in turn, can be generated by assuming violations of some fundamental principles of the current paradigm.

One of the most important principles of current physics is Lorentz invariance. This principle states that all inertial observers are equivalent. Lorentz invariance is an action symmetry and is closely related to the proposal that all fields in the action must be dynamical~\cite{Cristobal1, Cristobal2}, thereby avoiding  preferred spacetime directions. Relevantly, there are works within several promising quantum gravity candidates that argue that Lorentz invariance is broken at energies close to the Planck scale \cite{KosteleckySamuel, GambiniPullin}.

Moreover, there exists a model-independent framework to study (explicit) violations of Lorentz invariance, from hereon referred to as Lorentz violation. The idea, pioneered by Kostelecky and collaborators, is to create a parametrisation in the effective field theory framework~\cite{Rob}; the result is known as the Standard Model Extension (SME)~\cite{SME1, SME2, Kostelecky2004}. As such, the SME action contains the action terms of conventional physics, that is, the Standard Model and general relativity, plus all corrections that can be written in terms of the fields of this conventional part, but that includes terms that are not Lorentz invariant. The Lorentz-violating terms are weighted by nondynamical tensor fields, dubbed the SME coefficients, and are expected to be ranked in terms of their mass dimensions.

There is vast literature exploring a variety of Lorentz-violating effects (a list of the empirical bounds on the SME coefficients by many different experiments can be found in Ref.~\cite{DataTables}). There are also interesting theoretical considerations~\cite{PhysRevD.63.065008}, particularly when considering dynamical spacetimes or spontaneous Lorentz violation \cite{SMEtheo1,SMEtheo2,SMEtheo3,SMEtheo4,SMEtheo5,SMEtheo6,SMEtheo7,SMEtheo8}. In this letter, we focus on a Lorentz-violating quantised field. In these regards, there is a paper where radiative corrections and the running of the SME coefficients are studied~\cite{MauroRalfRob}. However, the symmetries of the vacuum, which are the main topic of this paper, are not analysed in Ref.~\cite{MauroRalfRob}. We also refer the reader to Ref.~\cite{Husain16}, where the authors study how a detector may respond to a quantum field with modified dispersion relations, where the modifications are suppressed by a large mass. Here, we quantise a free Lorentz-violating scalar coupled to a dimensionless SME-like coefficient, which is the simplest Lorentz-violating theory, and examine its vacuum states to verify whether they are unitarily equivalent. We find that the vacuum state in a particular Lorentz frame is not a vacuum in a frame connected by a Lorentz transformation. This raises serious questions on the phenomenological viability of Lorentz-violating theories, even if they are merely seen as low-energy effective theories.

To gain a deeper insight, we study whether different inertial observers can feel such effects. We equip our inertial observers with Unruh-De Witt detectors to probe the particle content of an inertial vacuum in different inertial frames. We find that, despite the initial concerns, inertial detectors remain silent upon the change of frames. Following Unruh's motto that `a particle is what a particle detector detects,' in this simple Lorentz-violating model, inertial observers experience the vacuum associated with one fixed inertial frame as devoid of particles in all inertial frames.

Given that it may be interesting to extend the above results to perturbative interacting theory, we complement the aforementioned results with a discussion on the Hadamard property in this Lorentz-violating theory. This condition, as it turns out, does not hold strictly and needs to be adapted. Nevertheless, we see no obstruction in using this `modified' Hadamard property to construct interacting observables (including time-ordered products of non-linear observables) and explore these questions in full perturbative interacting theory.

Throughout the paper, Latin indices from the beginning of the alphabet represent abstract indices, namely, they stand for open slots in the associated tensors and are thus not associated with any coordinates. In contrast, Greek indexes represent the components of the corresponding tensors in a given coordinate basis. For simplicity, the spacetime under consideration is the flat $4$-dimensional spacetime (without torsion or nonmetricity), whose metric (inverse metric) is denoted by $\eta_{ab}$ ($\eta^{ab}$), and it has signature $(-+++)$. The covariant derivative is denoted by $\nabla_a$.

\section{The model}

Although our model is not part of the SME, which deals with fields present in the Standard Model, it is instructive to study the simplest possible setting: that of a quantum scalar field $\phi$. This example already contains valuable lessons for other types of matter fields. In analogy with SME spinors, one can propose the Lorentz-violating Lagrangian
\begin{equation}
\mathcal L=-\frac{1}{2}\left[\left(\eta^{ab}+c^{ab}\right)\nabla_a\nabla_b\phi^2+m^2\phi^2\right],
    \label{scalarLagrangian}
\end{equation}
where $c^{ab}$ is the SME coefficient, which generates an effective metric for $\phi$. A particular $c^{ab}$, known in the SME literature as an isotropic coefficient, happens when there is a global set of Minkowskian coordinates $(t,x,y,z)$ for which the components of $c^{ab}$ take the simple form
\begin{equation}
c^{\mu\nu}=\mathrm{diag}(C,0,0,0),
\label{particular}
\end{equation}
 where $C$ is a constant. In these `privileged' coordinates, the Lagragian~(\ref{scalarLagrangian}) is
\begin{equation}
\mathcal L=-\frac{1}{2}\left[|\nabla\phi|^2-(1-C)(\partial_t\phi)^2+m^2\phi^2\right].
    \label{privLag}
\end{equation}

In the case of the scalar field above, one writes the Klein-Gordon product between two solutions $\phi$ and $\psi$ in the privileged coordinates as
\begin{equation}
    (\phi,\psi)=-\mathrm{i}(1-C)\int(\psi\partial_t\bar\phi-\bar\phi\partial_t\psi)\ \mathrm{d}^3x,
    \label{kg}
\end{equation}
where we integrate over a constant $t$ surface and the bar stands for complex conjugation. In fact, for technical convenience, we quantise the fields in a toroidal box of length $L$, so that integration is taken over the cube $[-L/2,L/2]^3$. Yet, observable quantities must be independent of $L$.

We can readily check that the modes
\begin{equation}
    \phi_{\vec k}=[L^32\omega(1-C)]^{-1/2}\exp{\mathrm{i}(\vec k\cdot \vec x-\omega t)},
    \label{privmodes}
\end{equation}
with
\begin{equation}
\omega:=\sqrt{\frac{m^2+\vec k^2}{1-C}}
    \label{privdisp}
\end{equation}
for $k_i=2\pi j_i/L$, $j_i\in\mathbb Z, i=1,2,3$. In addition, the dot in Eq.~\eqref{privmodes} stands for the standard $3$-dimensional scalar product. Note that these solutions, together with their complex conjugates, comprise a complete set of solutions to the field equations derived from the Lagrangian~(\ref{privLag}). The modes are normalised so that $\left(\phi_{\vec k_1},\phi_{\vec k_2}\right)=\delta_{\vec k_1\vec k_2}$. Canonical quantisation can be carried out in the standard way, as we shall see in the next section.

\section{Inequivalence of inertial vacua}
\label{Sec:Vacua}

In standard quantum field theory, the existence of a Lorentz-invariant vacuum is of fundamental importance, so much so that the existence of such a unique Poincar\'e-invariant state was adopted as an \emph{axiom} in the axiomatic formulation of quantum field theory in Minkowski spacetime~\cite{Streater64}. This property is lost in Lorentz-violating theories.

The modes~(\ref{privmodes}) are adapted to the privileged coordinates in the sense that the Lie derivative of $\phi_{\vec k}$ along $\partial_t$ vanishes. One can likewise adopt a set of solutions that are Lie-dragged according to $\partial_{t^\prime}$, which is obtained from $\partial_t$ by a (global) Lorentz boost. These are
\begin{equation}
\varphi_{\vec k^\prime}=\mathcal N_{\vec k^\prime}\exp{\mathrm{i}(\vec k^\prime\cdot\vec x^\prime-\omega^\prime t^\prime),}
    \label{unpmodes}
\end{equation}
where $\mathcal N_{\vec k^\prime}$ are normalisation constants. Crucially, reflecting the nonequivalent equation of motion, the unprivileged dispersion relation is altered, it solves
\begin{equation}
    (\omega^\prime)^2-(\vec k^\prime)^2-m^2+C[-\gamma^2(\omega^\prime)^2-2v\gamma k_z^\prime\omega^\prime-\gamma^2v^2(k_z^\prime)^2]=0,
    \label{undisp}
\end{equation}
where $\gamma=1/\sqrt{1-v^2}$ and the boost direction was chosen to be, without loss of generality, to be $+z$.

To find $\mathcal N_{\vec k}$, first, we note that, from the Lagrangian written in terms of the primed coordinates, the Klein-Gordon product reads
\begin{multline}
    \left(\phi,\psi\right)^\prime=-\mathrm{i}\int\mathrm{d}^3\vec x^\prime\left\{\psi[(1-C-v^2)\gamma^2\partial_{t^\prime}\bar\phi+\gamma^2vC\partial_{z^\prime}\bar\phi]\right.\\
    -\left.\bar\phi[(1-C-v^2)\gamma^2\partial_{t^\prime}\psi+\gamma^2vC\partial_{z^\prime}\psi]\right\},
    \label{unkg}
\end{multline}
where integration now is on a $t^\prime$ surface. With this product we get
\begin{equation}
\left(\varphi_{\vec k_1^\prime},\varphi_{\vec k_2^\prime}\right)^\prime=2|\mathcal N_{\vec k1^\prime}|^2\frac{L^3}{\gamma}
\left[(1-C-v^2)\gamma^2\omega_1^\prime-C\gamma^2vk_{z1}^\prime\right]\delta_{\vec k_1^\prime \vec k_2^\prime}.
    \label{unpnorm}
\end{equation}

Since both sets of solutions $\{\phi_{\vec k},\bar\phi_{\vec k}\}$ and $\{\varphi_{\vec k^\prime},\bar\varphi_{\vec k^\prime}\}$ are complete, we may write
\begin{equation}
\phi_{\vec k}=\sum_{\vec k^\prime}\left(\alpha_{\vec k\vec k^\prime}\varphi_{\vec k^\prime}+\beta_{\vec k\vec k^\prime}\bar\varphi_{\vec k^\prime}\right).
    \label{bogo}
\end{equation}
Let $\boldsymbol\Phi$ denote the field operator. We write
\begin{equation}
    \boldsymbol\Phi=\sum_{\vec k}\left(\phi_{\vec k}a_{\vec k}+ \mathrm{h.c.}\right)=\sum_{\vec k^\prime}\left(\varphi_{\vec k^\prime}a^\prime_{\vec k^\prime}+\mathrm{h.c.}\right),
    \label{modeexp}
\end{equation}
from which we find the corresponding Bogoliubov coefficients (see e.g. Ref. \cite{Bire82})
\begin{equation}
\alpha_{\vec k\vec k^\prime}=\left(\varphi_{\vec k^\prime},\phi_{\vec k}\right)^\prime,\qquad\beta_{\vec k\vec k^\prime}=-\left(\bar\varphi_{\vec k^\prime},\phi_{\vec k}\right)^\prime.
    \label{alphabeta}
\end{equation}
Moreover, from Eqs.~(\ref{bogo}) and (\ref{alphabeta}), we can
\begin{equation}
a_{\vec k}^\prime=\sum_{\vec k^\prime}\ \left(\alpha_{\vec k\vec k^\prime}a_{\vec k^\prime}+\overline{\beta_{\vec k\vec k^\prime}}a^\dagger_{\vec k^\prime}\right).
    \label{annihilation}
\end{equation}

Clearly, the vacuum annihilated by $a_{\vec k}^\prime$ does not coincide with the one annihilated by $a_{\vec k}$ when at least one $\beta_{\vec k\vec k^\prime}\neq0$. Furthermore, denoting by $|0\rangle$ the state for which $a_{\vec k}|0\rangle=0\ \forall\vec k$, we see that $a^\prime_{\vec k}|0\rangle\neq0$, proving the vacuum state is not Lorentz-invariant, unlike in conventional theories. We can go further and notice that the number operator $N_{\vec k}^\prime=(a^\prime_{\vec k})^\dagger a^\prime_{\vec k}$ satisfies
\begin{equation}
\langle 0|N_{\vec k}^\prime|0\rangle=\sum_{\vec k^\prime}\left|\beta_{\vec k\vec k^\prime}\right|^2 \left(\phi_{\vec k^\prime},\phi_{\vec k^\prime}\right),
    \label{number}
\end{equation}
where, in the second step, we have applied Eq.~(\ref{annihilation}) and the commutation relation
\begin{equation}
    \left[a_{\vec k^\prime_1}, a^\dagger_{\vec{k_2}^\prime}\right]=\left(\phi_{\vec k_1^\prime},\phi_{\vec k_2^\prime}\right)\mathbf 1.
\label{commutation}
\end{equation}

Again, whenever the coefficients $\beta_{\vec k\vec{k^\prime}}$ are nonzero, the preferred vacuum $|0\rangle$ has a nonzero number of particles in the modes adapted to the unprivileged frame. However, since the quantisations in the two frames are independent, one must evaluate how such particles interact with other matter fields. Indeed, as we shall see in Sec.~\ref{Sec:Detector}, these particles cannot excite a detector under reasonable hypotheses.

From Eqs.~(\ref{unpmodes}), (\ref{undisp}), (\ref{alphabeta}), (\ref{number}), and choosing the normalisation constant in Eq.~(\ref{unpnorm}) so that the resulting solution basis is orthonormal, it can be shown that
\begin{multline}
\beta_{\vec k\vec k^\prime}=\frac{8\mathrm{e}^{-\mathrm{i}(\omega^\prime-\gamma(\omega+v k_z))t^\prime}}{2L^3\sqrt{\omega(1-C)/\gamma\  [(1-C-v^2)\gamma^2\omega^\prime-C\gamma^2v k_z^\prime]}}\times\\
\frac{\sin(k_x-k_x^\prime)L}{k_x-k_x^\prime}\frac{\sin(k_y-k_y^\prime)L}{k_y-k_y^\prime}\frac{\sin(k_z^\prime-\gamma(k_z+v\omega))L/\gamma}{k_z^\prime-\gamma(k_z+v\omega)}\\
\times\left[(1-C-v^2)\gamma^2(\omega^\prime+\gamma(k_zv-\omega))+C\gamma^2v(\gamma(k_z-\omega v)-k_z^\prime)\right].
    \label{beta}
\end{multline}
After imposing the discrete values for $\vec k$ and $\vec k^\prime$, the last three factors on the first line of Eq.~(\ref{beta}) become
$L^3\delta_{k_xk^\prime_x}\delta_{k_yk^\prime_y}\delta_{k_z,\gamma(k_z^\prime+v\omega^\prime)}.$ Thus,
\begin{equation}
\langle 0|N^\prime_{\vec k}|0\rangle=\frac{16[(1-C-v^2)\gamma^2(\omega^\prime-\gamma(\omega-k_z v))-C\gamma^2v(k_z^\prime-\gamma(k_z-\omega v))]^2}{\omega(1-C)[(1-C-v^2)\gamma\omega^\prime-C\gamma vk_z^\prime]},
    \label{number2}
\end{equation}
with a small abuse of notation, where $k_z^\prime$ is not free, but rather obtained from $k_z$ from the corresponding Lorentz transformation, and, as before, $\omega$ and $\omega^\prime$ are calculated from $\vec k$ from their corresponding dispersion relations.

Importantly, from Eq.~(\ref{undisp}), we can readily verify that, when $C=0$, then $\omega^\prime=\gamma(\omega-vk_z)$ and $k_z^\prime=\gamma(k_z-v\omega)$, which are the conventional Lorentz-invariant results. Moreover, since, in this case, the factor in the square brackets in Eq.~(\ref{beta}) vanishes for the diagonal elements, all the $\beta_{\vec k\vec k^\prime}$ vanish identically, reflecting the equivalence of all inertial vacua for theories that preserve the Lorentz symmetry.

Notice that expression~(\ref{number2}) is ill-behaved for $C>1$ when the `effective' metric $\tilde\eta_{ab} = \eta_{ab} + c_{ab}$ ceases to be Lorentzian, and the field equations to be hyperbolic. It also diverges when $C<0$ and the phase velocity in the `primed frame', $\omega^\prime/k_z^\prime$ reaches $Cv/(1-C-v^2)$. These effects are closely related to vacuum \v{C}erenkov radiation, which is well understood in the context of the SME \cite{VacCerenkov1, VacCerenkov2}.

Another interesting feature of Eq.~(\ref{number2}) is that, for large values of $k_z$, it approaches a constant. Figure~\ref{fig:my_label} shows a plot of Eq.~(\ref{number2}) for the particular case $k_x=k_y=0$. Thus, one can see that the number operator contains excitations of arbitrarily large momentum. We turn now to study the interaction of the Lorentz-violating scalar field with inertial Unruh-De Witt detectors.

\begin{figure}
    \centering
    \includegraphics[width=.65\textwidth]{Nkk4.png}
    \caption{$\langle0|N^\prime_{\vec k}|0\rangle$ as a function of $k_z$ for $m=1$ (blue curve) and $m=0$ (red curve), and with $k_x=k_y=0$. We used the parameters $C=0.1$ and $v=0.8$. Note that $\langle0|N^\prime_{\vec k}|0\rangle$ approaches nonzero constants as $k_z\to\pm\infty$.}
    \label{fig:my_label}
\end{figure}

\section{Inertial Unruh-De Witt detectors}
\label{Sec:Detector}

If one wishes to extract observable consequences of a theory, one must allow for interactions. The simplest model that captures the physics is an Unruh-De Witt detector \cite{Bire82}: a pointlike two-level quantum system that follows a spacetime worldline as it interacts with the quantum field. The energy eigenvectors of the detector form a basis for its two-level Hilbert space. We denote them by $|E_0\rangle $ and $|E_1 \rangle$ with energies $E_0$ and $E_1>E_0$, respectively.

If we assume that the detector and the field are constantly coupled for an arbitrarily large proper time of the detector, we can write the interaction Hamiltonian as
\begin{equation}
H_I=\lambda\mu(\tau)\otimes\phi(x(\tau)),
    \label{UDWscalar}
\end{equation}
where $\mu$ denotes an internal degree of freedom of the detector, usually known as its monopole moment, $\lambda$ is a coupling constant, and $\phi(x(\tau))$ denotes pullback of the filed to the worldline of the detector parametrised by its proper time $\tau$.

In first-order perturbation theory, the transition amplitude from the detector's ground state $|E_0\rangle$ to the excited state $|E_1\rangle$ is
\[ A^{(1)}:= -\mathrm{i}\lambda\langle E_1|\langle\boldsymbol\Psi|\int_{-\infty}^\infty\mathrm{d}\tau\ \mu(\tau)\phi(x(\tau))\ |0\rangle|E_0\rangle\]
for each field state $|\boldsymbol\Psi\rangle$. As before, $|0\rangle$ denotes the `privileged' vacuum. This last equation factorises to
\begin{equation}
  A^{(1)}=  -\mathrm{i}\lambda\langle E_1|\mu(0)|E_0\rangle\int\mathrm{d}\tau\ \mathrm{e}^{\mathrm{i}(E_1-E_0)\tau}\langle\boldsymbol\Psi|\phi(x(\tau))|0\rangle
   \label{Tamplitude}
\end{equation}
so the only nonzero matrix elements are
\[M_k:=\langle 0|a_k\phi(x(\tau))|0\rangle=[L^32\omega(1-C)]^{-1/2}\mathrm{e}^{\mathrm{i}(\omega t-\vec k\cdot\vec x)},\]
where, in the second step, we used the expansion in modes given in Eq.~(\ref{privmodes}). Here, $t$ and $\vec x$ are coordinates along the detector's worldline, $x^\mu(\tau)=t\delta^\mu_0+vt\delta^\mu_z$, with $t=\gamma\tau$.

It can be straightforwardly verified that the other terms in the Dyson series~\cite{Reed75} share the form
\begin{multline}
 A^{(n)}:=   \frac{(-\mathrm{i}\lambda)^n}{n!}\langle E_1|\langle\boldsymbol\Psi|\int\mathrm{d}\tau_1\int\mathrm{d}\tau_2\ldots\int\mathrm{d}\tau_n\ \theta(\tau_{n-1}-\tau_n)\ldots\theta(\tau_2-\tau_1)\\
 H_I(x(\tau_1)\cdot\ldots\cdot H_I(x(\tau_n))|0\rangle|E_0\rangle.
    \label{odern}
\end{multline}

To order $n$, all the matrix elements for which $|\boldsymbol{\Psi}\rangle$ is not the product of exactly $n$ operators $a^\dagger_k$ acting on $|0\rangle$ are trivially zero. The other matrix elements are of the form
\begin{multline}
  A^{(n)}=      \frac{(-\mathrm{i}\lambda)^n}{n!}\langle E_1|\mu(0)^n|E_0\rangle
        \int_{-\infty}^\infty\mathrm{d}\xi\ \mathrm{e}^{\mathrm{i}[(E_1-E_0)+\frac{\gamma}{n}\sum_{j=1}^n(\omega_j-\vec k_j\cdot\vec v)]\xi}\times\\
        \int_0^\infty\ldots\int_0^\infty\prod_{i=1}^{n-1}\mathrm{d}\eta_i\ \frac{(-1)^{n+1}}{n}\ \sum_{\sigma,i,l}(-1)^{\sigma}\sigma\left(\frac{\exp\left\{-\mathrm{i}\gamma[\sum_{r=1}^{n-1}k_lx(P^{-1}_{i,r+1}\eta_r)]\right\}}{\sqrt{L^32\omega_i(1-C)}}\right),
        \label{interm}
\end{multline}
on which we introduce the variables $\xi=\sum_{i=1}^n\tau_i$ and $\eta_i=\tau_i-\tau_{i+1}$ for $i=1,2,\ldots,n-1$. Also, $\sigma$ denotes a permutation of $i$ and $l$, and $P$ denotes the matrix implementing the transformation form $(\xi,\eta_i)$ to $\tau_i$. We also use the equalities $\det P=(-1)^{n+1}n$ and $P_{i1}^{-1}=1/n\ \forall i$.

Equation~(\ref{interm}) assumes $n\geq2$. For $n=1$, only the integral over $\xi$ remains. In any case, it yields a factor \[\delta\left(E_1-E_0+\frac{\gamma}{n}\sum_{j=1}^n(\omega_j-\vec k_j\cdot\vec v)\right).\]
Given that $E_1>E_0$, if the detector moves more slowly than the phase velocity of each field mode (which coincides with the group velocity for linear dispersion relations),
\begin{equation}
v<\omega_j/k_j \quad\forall j,
\label{bound}
\end{equation}
it follows that the argument of the delta function is strictly positive. Consequently, the inertial detector cannot detect particles in any order in the Dyson series.

This result can be generalised to any field satisfying the homogeneous Dirac or Klein-Gordon equation and couplings other than Eq.~(\ref{UDWscalar}). A particular realisation of condition~(\ref{bound}) is when the phase velocities are bounded by the speed of light, in which case Eq.~(\ref{bound}) holds for all $v<1$. In our example, the necessary and sufficient condition for this realisation is $0\leq C<1$\footnote{Indeed, let $\alpha=1/\sqrt{1-C}$. The condition we wish to satisfy is $\alpha\sqrt{\vec k^2+m^2}\geq|\vec k|\ \forall |\vec k|\geq 0$. Since the left-hand side is a convex function which is larger than $|\vec k|$ for $|\vec k|=0$, we need the equation $\alpha\sqrt{\vec k^2+m^2}=|\vec k|$ to have no real solutions in this range. Writing the discriminant of the square of the equation, we see that this happens for $\alpha\geq 1$.}.

There are, of course, limitations in using Unruh-De Witt detectors. Nevertheless, they function as a technically simpler tool for identifying cases of interest. For example, if the Lorentz-violating sector responds to an effective metric $\tilde\eta_{ab}$ that is not flat, the transition amplitudes~(\ref{Tamplitude}) could be nonzero. Likewise, the question of how noninertial detectors respond in Lorentz-violating situations could also have phenomenological implications in the future, once we have the experimental capacity of measuring the Unruh and related effects of noninertial motion~\cite{Good:2020hav} (see, e.g., Ref.~\cite{Gooding:2020scc} for a proposal in this direction).

An alternative procedure to look for observable consequences of Lorentz violation consists in violating the bound~(\ref{bound}). A direct calculation of the Wightman function written in terms of the privileged coordinates in the massless case reveals a divergence when $v=1/\sqrt{1-C}$. In the detector model, such divergence can influence the rate of spontaneous de-excitation~\cite{Benito17}, but this is a regime where other effects like vacuum \v{C}erenkov radiation~\cite{VacCerenkov1, VacCerenkov2} are also expected to occur.

This section shows that, despite an `unprivileged' inertial detector describing the vacuum $|0\rangle$ as a sea of particles following the distribution~(\ref{number2}), it does not interact with them, except for detectors at least as fast as the propagation of the field modes. We turn now to analyze the singular structure of the two-point function in the theory at hand.

\section{The Hadamard property and Lorentz-violating theories}

A property of central importance for free and perturbatively-interacting quantum field theories is the Hadamard condition, which characterises the distributional singularities of the two-point function for physically reasonable states \cite{Fewster:2013lqa}. In this section, we examine how the Hadamard condition becomes altered for the Lorentz-violating theory under consideration. The Hadamard property is typically formulated in terms of the Wightman function, which, for a state vector $\Psi$ in a concrete representation, is $\omega_2(x,x') = \langle \Psi | \bf \Phi(x) \bf \Phi(x') | \Psi \rangle$.

For the theory defined by the Lagrangian~(\ref{scalarLagrangian}), the Wightman function is a bi-solution of the field equation
\begin{equation}
P \otimes 1 \omega_2 = 1 \otimes P \omega_2 = 0,
\label{P-op-bi}
\end{equation}
where the differential operator, $P$, is defined as
\begin{equation}
P = \left(\eta^{ab}+c^{ab}\right)\nabla_a\nabla_b-m^2 = \tilde \eta_{ab} \nabla_a \nabla_b - m^2.
\label{P-op}
\end{equation}
Direct inspection of the characteristic set of the differential operator $P$ indicates that the standard Hadamard property does not hold in Lorentz-violating theories. Indeed,
the principal symbol of $P$ is given by $p(x,k) = \tilde \eta^{ab}(x) k_a k_b$ and hence the characteristic set,
\begin{equation}
\mathcal{N} = {\rm Char} P = \{ (x,k) \in \mathbb{R}^4 \times \mathbb{R}^{4 *}\setminus{0}: \eta^{ab}(x) k_a k_b = 0  \},
\end{equation}
lies along nonvanishing null covectors with respect to the effective metric $\tilde \eta_{ab}$, and not along the standard (nonvanishing) null vectors of Minkowski spacetime.

The above observations suggest that it is convenient to study the Hadamard property from the microlocal viewpoint introduced by Radzikowski \cite{Radzikowski}, which might be adapted more straightforwardly to the Lorentz-violating case. Also, it appears to us that this has the advantage that the result can be directly generalised to a large class of Lorentz-violating theories. By a standard result in microlocal analysis \cite{Hormander}, any distribution defined on smooth functions of compact support, $\varphi: C^\infty_0 (\mathbb{R}^4) \to \mathbb{R}$, that satisfies $P \varphi = 0$, must have a wavefront set\footnote{The wavefront set of a distribution is defined as follows: If $\varphi$ is a distribution as above, we call $(x, k)$ with $k \neq 0$ a \emph{regular direction} of $\varphi$ if one can find $u \in C_0^\infty(\mathbb{R}^4)$ with support including $x$, such that the Fourier transform of the pointwise product $\varphi u$ decays rapidly in a conic neighbourhood of $k$. The wavefront set of $\varphi$, denoted by ${\rm WF}(\varphi)$, is then defined as the set of pairs $(x,k)$ with $k \neq 0$, such that $(x,k)$ are not regular directions of $\varphi$.} satisfying
\begin{align}
    {\rm WF} (\varphi) \subset \mathcal{N}.
    \label{WFset1}
\end{align}
In fact, this analysis is valid for any second-order hyperbolic differential operator in Lorentz-violating theories whose principal symbol is given in terms of an effective Lorentz-violating metric, analogous to $\tilde \eta$ in Eq.~(\ref{P-op}). Thus, we can continue our discussion in full generality.

From Eq.~(\ref{P-op}), it can be verified that any distributional singularities of a distributional solution of the Lorentz-violating Klein-Gordon equation propagate along the `effective' null cones of $\tilde \eta_{ab}$, where we work under the reasonable assumption that $c_{ab}$ does not alter the Lorentzian character of $\tilde \eta_{ab}$. This result can be extended straightforwardly to distributional bi-solutions, such as $\omega_2$. If $\varpi$ is a distributional bi-solution, one is now led to a wavefront set condition satisfying that
\begin{align}
{\rm WF}(\varpi) \subset {\cal N}_0 \times {\cal N}_0,
\end{align}
where ${\cal N}_0$ is equal to ${\cal N}$, but including now in addition zero co-vectors, i.e., of the form $(x,0)$.

With this in mind, a direct generalisation of the Lorentz-invariant case suggests that the correct microlocal spectrum condition for the (positive frequency) Wightman function, $\omega_2$, should further satisfy
\begin{equation}
{\rm WF}(\omega_2) \subset \cal N^+ \times \cal N^-,
\end{equation}
where $\cal N^+$ is the set of covectors in $\cal N$ with future-pointing momenta and $\cal N^-$ is the set of covectors in $\cal N$ with past-pointing momenta \cite{Radzikowski}.

As it can be seen, so far, all the results for characterising the wavefront set of $\omega_2$ are as in the Lorentz invariant case, but replacing the Minkowski metric by the effective metric $\tilde \eta$. We can therefore directly write down the microlocal spectrum condition \cite{Radzikowski} for Lorentz-violating theories,
\begin{multline}
{\rm WF}(\omega_2)  = \{ ((x_1, k_1),(x_2 k_2)) \in (\mathbb{R}^4 \times \mathbb{R}^{4 *}\setminus{0})
 \times (\mathbb{R}^4 \times \mathbb{R}^{4 *}\setminus{0}):\\
 (x_1,k_1) \sim (x_2, -k_2),
 \text{ and } k_1 \text{ is future pointing}\},
\label{MLSC}
\end{multline}
where the notation here is that $(x,k) \sim (y,p)$ if there is a curve $\gamma$ whose image contains the points $x$ and $y$, such that the tangent covector field to $\gamma$ satisfies that $\tilde \eta^{ab} v_a v_b $ (note that here we use the effective metric and not the Minkowski metric) at every point along the curve, and that $k$ and $p$ are parallel transports along $\gamma$.

In short, this section shows explicitly that the singular structure of the two-point function is that of the Lorentz-invariant case, but associated with the causal structure of the effective metric $\tilde \eta^{ab}$.

\section{Discussion}

We analysed the quantisation of a very simple Lorentz-violating theory. We found that the vacuum, as perceived by an inertial observer, is populated from the viewpoint of a different inertial frame, refuting the idea that the vacuum is always Lorentz invariant. This result could have had profound empirical consequences, had it had physical effects on observers moving in different frames. To this end, we examined the response of an observer equipped with an Unruh-DeWitt detector and showed however that it is insensitive---to any order in perturbation theory---to the excitations of this scalar field. Following the maxim that `a particle is what a particle detector detects' the vacuum appears to be empty in all Lorentz frames, after all. So, the answers to the questions in the title are negative, at least for the model we explicitly solved.

In addition to all the former, we discussed how the Hadamard property needs to be adapted in any theory whose effective light cone does not coincide with that of the Minkowski metric, finding {\it a priori} no obstructions to treating Lorentz-violating fields in interacting quantum theory or for the construction of renormalised non-linear observables. This is an interesting open question that can be examined in the future as a test of the internal mathematical consistency of Lorentz-violating theories.

It should be straightforward to generalise this study to other maximally symmetric spacetimes, where there are preferred vacua. More interestingly, we believe that the presented results can be generalised to more realistic theories without finding major conceptual or technical difficulties. In this sense, we think that nonequivalent inertial vacua are present in the actual SME. Notice that more general SME coefficients (e.g., spacetime-dependent coefficients) could have detectable signals of the form considered here, which could lead to promising experimental tests. In this sense, we urge the proponents of particular Lorentz-violating models to consider the corresponding interacting quantum field theory and to verify their mathematical and empirical adequacy in light of the lessons of this work.

\section*{Acknowledgments}
The authors wish to thank Lukas Nellen for his helpful ideas. This research was partially supported by UNAM-DGAPA-PAPIIT Grant No. IG100120, CONACyT FORDECYT-PRONACES Grant No. 140630. BAJ-A is funded by a CONACYT Postdoctoral Fellowship.

\bibliographystyle{unsrt}