{'fraction_non_alphanumeric': 0.03242769500438212, 'fraction_numerical': 0.025381244522348818, 'mean_word_length': 4.762828282828282, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A Search for Technosignatures toward the Galactic Centre at 150 MHz

This paper is the fourth in a series of low-frequency searches for technosignatures. Using the Murchison Widefield Array over two nights we integrate 7 hours of data toward the Galactic Centre (centred on the position of Sagittarius A*) with a total field-of-view of 200 deg^2. We present a targeted search toward 144 exoplanetary systems, at our best yet angular resolution (75 arc seconds). This is the first technosignature search at a frequency of 155 MHz toward the Galactic Centre (our previous central frequencies have been lower). A blind search toward in excess of 3 million stars toward the Galactic Centre and Galactic bulge is also completed, placing an equivalent isotropic power limit<1.1X10^19W at the distance to the Galactic Centre. No plausible technosignatures are detected.


INTRODUCTION
The prevalence of life beyond Earth is a central and unanswered question within astrobiology. The search for extraterrestrial intelligence (SETI) seeks to answer this question via detection of "technosignatures", artificial signals that indicate the existence of technologicallycapable societies (see review by Tarter, 2001). On Earth, low-frequency radio signals, like those used by FM radio, are a ubiquitous choice for communications. Many astrophysical processes give rise to low-frequency radio emission, and as such numerous large and sensitive low-frequency radio telescopes have been built, including the current-generation Murchison Widefield Array (MWA, Tingay et al., 2013;Wayth et al., 2018), Long Wavelength Array (LWA, Ellingson et al., 2009), Low-Frequency Array (LOFAR, van Haarlem et al., 2013) and Giant Metrewave Radio Telescope (GMRT, Gupta et al., 2017). The existence of both powerful transmitters and sensitive receivers at low frequencies-both of which emerged early in the history of radio engineeringmotivates low-frequency technosignature searches by providing an example class of engineered signals to search for, and instruments with which to do so. This paper is the fourth in a series of papers detailing * astrochenoa@gmail.com SETI observations with the Murchison Widefield Array (MWA, Tingay et al., 2013;Wayth et al., 2018), the details of which are summarised in Table 1. The MWA offers two advantages over other SETI searches; its large field of view and the low frequency range. These searches of ∼400-600 square degrees, are some of the largest published surveys, although no candidate technosignature signals were detected above the detection limits. Both Garrett et al. (2017) and more recently Houston et al. (2021) have discussed the benefits of using aperture arrays like MWA for efficiently completing an all-sky SETI survey. Houston et al. (2021) outlines strategies of SETI searches from past, present, and future and suggests that if a receiver and transmitter are aligned in "...space, time and frequency, with adequate receive power, a detection can occur." They suggest that unless there is a compelling reason to only search stellar regions, wide-field searches of any signal of unknown origin are required. However, before we get to all-sky technosignature searches there are a number of computational challenges to overcome and these surveys have provided insight on how to accomplish this goal with an aperture array. While each of the MWA SETI publications follows a similar processing and search approach, our data analysis techniques have been gradually and significantly improved. The observations toward Orion represented an improvement in imaging techniques and source finding. In the observations toward Vela, the data were collected with an updated "Phase II" array, increasing the spatial resolution by more than a third (3 arc minutes down to 1 arc minute).
In addition to the large field-of-view offered by the MWA, our surveys also represent the first published low radio frequency searches (see Figure 1). Since we don't know what frequency another technologically advanced civilisation might broadcast or operate at, there is no reason to ignore available search space. There are additional motivations for low-frequency as well. Sullivan et al. (1978) suggested that the FM radio broadcasting stations of the world represents the greatest power per hertz in the radio band and this was further explored by Loeb & Zaldarriaga (2007). Overall, there is growing support for broadening the frequency range searched for technosignatures to lower frequencies. For example, Houston et al. (2021) find that, according to several detection optimisation metrics, SETI surveys should be undertaken down into the hundreds of megahertz frequency range.
In this survey we utilise the procedures developed in our search toward Vela to search 200 deg 2 toward the centre of our Galaxy but at a higher frequency of 155 MHz. This survey maintains the higher spatial resolution we obtained toward Vela, but encompasses the largest population of known exoplanets than across our previous surveys.

The Galactic Centre
The Galactic Centre (GC) is a prime SETI target as the line of sight toward the GC has the largest integrated count of Galactic stars for any direction. Here, we outline some arguments favouring and disfavouring the GC as a region where intelligent life may reside; however, we note that direct observational evidence remains the only method capable of proving the existence of life beyond Earth.
The high density of stars within the GC means that cataclysmic events such as stellar supernovae and magnetar flares are more likely to impact exoplanets within the GC, potentially destroying any life on their surface. Based on these factors, Lineweaver et al. (2004) identify a 'Galactic habitable zone' as an annular region between 7-9 kpc from the GC. Interactions due to close stellar flybys-more common within the dense GC-are also expected to damage planetary disks (Lineweaver et al., 2004;Jiménez-Torres et al., 2013); however, on longer timescales this could be advantageous to habitability. Despite these factors, modelling by Gowanlock et al. (2011) finds a majority of planets that may support complex life are found toward the inner Galaxy (less than 1 kpc from the Galactic Centre). Morrison & Gowanlock (2015) extend this model to include intelligent life, and also find higher probability within the inner Galaxy. Gajjar et al. (2021) update the model of Gowanlock et al. (2011) to include the galactic bulge within the radius R < 2.5 kpc, again finding the fraction of stars with a habitable planet is greatest in the inner regions of the Galaxy. Modelling by Cai et al. (2021), which includes factors such as abiogenesis (the idea that life arose from nonlife), evolutionary timescales, and self-annihilation also find a higher likelihood of intelligent life emerging in the Galactic inner disk (defined here as R ≤ 8 kpc), with peak likelihood at an annulus 4 kpc from the GC.
Separate to the propensity of life to emerge, the GC's high stellar density may be advantageous to the growth of advanced space-faring societies. Such societies are likely to be capable of producing technosignatures detectable across large distances. First discussed in Newman & Sagan (1981), diffusion of advanced societies across the Milky way is accelerated in areas of high stellar density. Di Stefano & Ray (2016) suggest that the close proximity of stars within dense globular clusters is favourable to the diffusion of space-faring societies. Modelling by Carroll-Nellenback et al. (2019) also suggests that high stellar density provides a front for settlements to rapidly expand from. Nevertheless, dust grains, gas and cosmic rays all pose hazards to interstellar travels, making the GC harder to traverse (Lacki, 2021). The GC can also be considered a natural cynosure, or 'Schelling point' within the Milky Way: an optimal location to place a transmitter to maximise the chances of its detection (Gajjar et al., 2021).
Despite arguments that motivate SETI searches of the GC, only a handful of observations have been published. In 1981, the Westerbork Synthesis Radio Telescope was used to search for narrow band pulsating beacons in the GC (Shostak & Tarter, 1985). Our previous observations (Paper 1, Tingay et al., 2016) searched toward 38 known exoplanets at a 10 kHz resolution in 2 hours of data. As part of the Breakthrough Listen search for intelligent life, the GC is being surveyed using the Robert C. Byrd Green Bank Telescope (GBT) and Parkes 64 m telescope to cover 0.7-93GHz (Worden et al., 2017;Gajjar et al., 2021). Our observations, which are at lower frequency, complement the Breakthrough Listen search. Both Parkes and Green Bank are large single-dish telescopes with wide-bandwidth receivers. Interferometers, such as the MWA, offer wider fields of view and much better spatial resolution, giving improved signal localisation and radio frequency interference (RFI) rejection. However, compared to single-dish telescopes, interferometers have higher ingest data rates and require more demanding signal processing systems; consequently, their instantaneous bandwidth and spectral resolution are generally more constrained than single-dish instruments. As such, single-dish and interferometric approaches to technosignatures are complementary.
As shown in Figure 2 (Section 3.1), the known exo-  planets in this survey are distributed along either side of the Galactic plane (plus and minus longitude centred along b=0). Microlensing is the most common method of detection and, as such, this population represents a population of planets significantly more distant than some of our previous surveys.

OBSERVATIONS
The Murchison Widefield Array (MWA; Tingay et al. 2013) is a low-frequency interferometer operating between 70 and 300 MHz at the Murchison Radioastronomy Observatory in Western Australia. In 2018, the telescope was upgraded to the "Phase II" array (Wayth et al., 2018), doubling the number of aperture array tiles from 128 to 256 and approximately doubling the maximum baseline from 3 km to 5.5 km, of which 128 tiles are correlated at one time. These observations use the extended "long base-line array" with baselines between ∼22 m and 5.5 km.
Observations of the Galactic Centre, centred on J2000 coordinates 17:45:40.04 −29:00:28, were taken during 5-7 September 2020 for a total of 10.5 hours and are summarised in Table 2. The observations were taken when the Galactic Centre was low on the horizon 1 , so only 7 hours of the most sensitive observations were processed. The data were calibrated from 2-minute observations of Hercules A each night and processed as described in Tremblay et al. (2020). To image these observations we used w s c l e a n version 2.9.2 with Briggs weighting of −1. Due to the large volume of these data, we did not process using other weightings. This weighting improves our point-source sensitivity to reduce the potential effects of beam dilution over previous publications. These settings produced a cube with a field-of-view (FOV) of 200 square degrees and a synthesised beam of 75×67 arc seconds. Based on the bandpass solutions, three of the 128 tiles were flagged for every observation.
The MWA uses a two-stage polyphase filter-bank which channelises the data into 24 × 1.28 MHz "coarse" frequency channels which are then further divided into 128 × 10 kHz "fine" frequency channels resulting in 3072 (10 kHz-wide) spectral channels. Each of the 85 5-minute observations were imaged independently using W S C L E A N (Offringa & Smirnov, 2017) and a Stokes I primary beam model was created using the simulation-based Full Embedded Element (FEE) model of Sokolowski et al. (2017) for each of 24 coarse channels. The creation of the model, taking into account the phase centre and other individual observational parameters 1 The MWA observations are obtained through competitive time allocation and when the observations happen are based on scientific priorities and work around maintenance schedules. Therefore preferential observational time when the Galactic Centre was closer to zenith was not available. in the metadata, was created using code by Morgan & Galvin (2021). The primary beam model was applied to the 100 fine frequency channels imaged 2 within each coarse channel to create primary beam corrected spectral channel images. The images were then stacked together into a three-dimensional data cube for each observation, which is then time-averaged together using inverse variance weighting. After correction, the flux density error was determined by comparison to the Molonglo Reference Catalogue (MRC; Large et al. 1981) which was scaled down to our frequency using a spectral index of −0.83, as described in detail in Hurley-Walker et al. (2017). The error across the field was 9±20 mJy.
Radio Frequency Interference (RFI) was flagged in two stages in the u, v, w (visibility) domain. First, each 5-minute observation was flagged using a statistical technique invoked by ao f l ag g e r (Offringa et al., 2015) and is a standard option for data downloaded through the All-Sky Virtual Observatory. Secondly, a baseline flagging algorithm which flagged amplitudes in u, v, w of baselines that were higher than three times the standard deviation 3 . This flagged <1-3% of the baselines. In previous work , strong FM band interference required a third-stage image-based filtering to flag channels where RFI was detected in a large fraction of pixels; this third-stage flagging was not required for these data. An example spectrum is shown in Figure 4.
The ionosphere can create an astrometric shift in the sources which is more pronounced at lower frequencies.
The estimated source position error is estimated in the final composite image by comparing to the MRC. After correction, the systematic spatial error in a fully integrated continuum image was +5±8 arcseconds in right ascension and +22±5 arcseconds in declination. The ratio of point source flux density and peak intensity is   Tables 3, 4, 5, and 6 shown in Galactic Coordinates. As shown there are more sources distributed along the negative latitude which is likely a selection effect from the exoplanet survey axis fields.
1.1±0.2. This suggests that the ionospheric distortions are corrected for in these observations.

Known Exoplanets
The three-dimensional data cubes were each searched in frequency for 10 kHz (frequency resolution of the MWA) signals of non-astrophysical origin toward the 144 known exoplanets in our field. This list of exoplanets from the The Exoplanet Encyclopaedia 4 , represents an accumulated list of exoplanets with a mass limit of 60 M J . In Tables 3, 4, and 5 we provide information on the planet designation, position, distance, planet mass, detection method, sensitivity limit as spectral variation in our data cube toward that position, and the equivalent isotropic power limit (EIRP). In Table 6 we provide the information of exoplanets in our sample which have no known distance, and therefore a limit on the EIRP could not be obtained.
Similar to our survey toward Vela (Tremblay & Tingay, 2020), we calculated the EIRP upper limit from the equation: where S rms is the RMS intensity value in Jy beam −1 along the spectral frequency axis toward the source position and R is the distance to the stellar system in pc. Similar to our previous work, we assume a signal width of 10 kHz to match the frequency resolution of the MWA, and discuss in Section 4.1 about potential for spectral broadening when higher resolution is assumed. In these upper limits we assume that the signal optimally fills the 10 kHz channel width. If an artificial signal was an unresolved signal with a width between 1-10 kHz but constrained within the channel, the estimated EIRP would be underestimated by 10kHz/∆νt, where ∆νt is the transmission bandwidth. We do not account for additional edge cases where a signal spans the channel gap between two fine frequency channels. We note that for the EIRP Limit toward the known sources we use a one sigma value for direct comparison with previous work and a 6 sigma limit for the blind search. As shown in Figure 5, the majority of the stars containing exoplanets are near or below solar mass. It is also shown that the sample of known exoplanets have Flagged channels and the edge channels which suffer from instrumental aliasing (see Section 2), are blanked out in the spectrum. The horizontal grey shaded region represents the ±1σ RMS value used in Table 3.
a mean distance of 4.8 kpc, which is much farther than our sample of stellar objects we studied toward Vela or Orion.

Blind Signal Search
Each of the 2400 continuum-subtracted fine-channel (10 kHz) images, each containing 2000×2000 pixels, are independently searched using the source-finding software A e g e a n (Hancock et al., 2018). This is done using the function "slice", to set which channel in the cube is searched, and setting a "seed clip" value of 5, in order to search the image for pixels with a peak intensity value greater than 5 σ (where σ is set from an input RMS image). A e g e a n works by fitting Gaussian's to the pixel data and applies a correction for the background 5 to calculate the flux density for potential sources. This source-finding threshold has the goal of detecting all signals >6 σ, which may not be pixel centred. Following this search, we found two signals over the >10 9 voxels searched. This is consistent with expected number of spurious signals within the data set. We therefore, do not detect any unknown emission sources.
In Tremblay & Tingay (2020) we used the Gaia catalogue to approximate the number of stars in the field toward Vela. However, due to the high dust extinction toward the Galactic Centre and Galactic bulge, optical surveys have difficulty detecting sources.      Southern Observatory VLT HAWK instrument, they determined the precise photometric properties of 3.3×10 6 stars with an assumed distance of 8178 ± 13 stat pc (Gravity Collaboration, 2019). The GNS infrared survey only covers 0.15% of our field and accounts for in excess of 3.3 million stars. So our blind search will likely cover billions of stars (and background galaxies). Using the distance of 8.1 kpc, we place a limit on the EIRP of putative transmitters within our observing band of <1.1×10 19 W.

DISCUSSION
The median distance for the 144 known exoplanet systems toward the Galactic Centre is 5585 pc with the shortest distance at 39 pc. This is compared to the 28 pc for the 6 exoplanets in Tremblay & Tingay (2020), 50 pc for the 22 exoplanets examined by Tingay et al. (2018) and ≈2000 pc for the 45 exoplanets examined by Tingay et al. (2016). The EIRP upper limits in this paper represent some of the highest limits in our surveys to date, but as this is the first published survey in this frequency range we find this an important starting point.
There are several metrics which allow for context of the transmitter values and the potential for detecting other societies. An extrapolation of terrestrial technology from the Kardashev (Kardashev, 1964) scale, where a Type I civilisation with technology close to Earth is predicted to be able to emit an non-isotropically radiated signal at ∼ 10 17 ,W, is one suggested method. We note that such a civilisation would only need to be able to build a big dish. Sullivan & Knowles (1985) suggest that with an Arecibo-like single dish we could detect passive radio signals (or radio leakage) out to 9 pc and with a large array of smaller dishes the distance could be more than 10 times that. Overall, they concluded the most delectable signal from another planet would be a powerful military satellite.
Using examples on Earth, especially in the frequency range of the MWA, we can look at the Air Force Space Surveillance System known as "Space Fence" which operated up until 2013. This system was a 1 MW continuous wave (0.1 Hz BW) system operating at 216 MHz illuminating a 120 deg × 1.5 arc min field at an EIRP of approximately 1.5×10 10 W (Sullivan & Knowles, 1985). If we compare this with our upper limits for the closest source, we are still three orders of magnitude from detecting a signal of this strength.
Similar to the metric used by Sheikh et al. (2020), we could measure our limits against the Arecibo Planetary Radar experiment (∼20×10 12 W transmitted through a 305m parabolic reflector; L A ; Siemion et al. 2013). Their largest value for the EIRP limit for stars in 7-143 pc is an upper limit in L A of 0.88. However, with our frequency resolution of 10 kHz and the distances of 39 pc or greater, the smallest L A upper limit is ∼13. This is significantly larger than our results from our other MWA surveys, which represented a population of closer stars.
Currently the MWA is undergoing an upgrade which will be pertinent to future SETI advances. Over the next 6 months, the MWA will be upgraded to increase frequency and time resolution, which will better match the kHz and sub-kHz searches completed by other facilities. An alternative approach to narrow-band SETI, first suggested by Cole & Ekers (1979), is to search for impulsive wide-band technosignatures; however comparatively few searches have been conducted. Gajjar et al. (2021) search for artificial transient signals (0.7-194 ms duration) from the GC, using non-physical dispersion as a discriminant. Our interferometric data does not have sufficient time resolution for such a search, and at low frequencies signals are severely broadened on sight-lines (as discussed in Section 4.1) with high electron density, such as toward the GC. Nevertheless, the upgraded in higher time resolution modes which could be employed to search for impulsive low-frequency technosignatures toward the GC or other target fields.

Spectral broadening due to intervening media
Electron density fluctuations in the interstellar medium (ISM), solar wind and interplanetary medium (IPM) cause spectral broadening by scattering signals as they are propagating through them. This is of particular concern for low-frequency observations of the GC, as the effect is stronger at low frequencies, and the large electron density fluctuations in the GC cause strong scattering. Although the MWA capabilities at the time of these observations have a relatively wide frequency resolution of 10 kHz, future upgrades planned for 2021 will allow for sub-kHz resolution. We therefore explore here the potential impact of future observations in this section.
The effect of the ISM on narrowband signals is detailed in Cordes & Lazio (1991) and Cordes et al. (1997). In the strong scattering regime a narrowband sinusoid with frequency ν GHz in GHz will be broadened by ∆ν broad = 0.14 Hz ν −6/5 GHz where ∆ν broad is the spectral broadening in Hz, V ⊥ is the transverse velocity of the source in km s −1 , and SM is the scattering measure along the line of sight. The SM is an integrated measure of electron density fluctuations C 2 ne along the line of sight to a distance L, (Rickett, 1990), defined as Estimates of SM toward the GC can be found using the NE2001 Galactic electron density model (Cordes & Lazio, 2002). To estimate the spectral broadening of narrowband signals at MWA frequencies, we use the P y G E D M code (Price et al., 2021) to generate SM from NE2001 as a function of distance toward the GC (Fig. 6, left). To avoid a discontinuity at the GC (l=0, b=0), we add a 0.5 deg offset in Galactic latitude when calculating SM. The highest radial velocities within the Galaxy occur toward the GC, reaching |V ⊥ | ≈280 km s −1 (Dame et al., 2001). Applying Eq. 2, we find spectral broadening of under 7.5 Hz at an observational frequency of 150 MHz: well below our channel resolution of 10 kHz. Following Siemion et al. (2013), spectral broadening due to the IPM can be estimated as ∆ν broad = 300 Hz ν where R is the radius of the Sun, and R is the solar impact distance. The R −9/5 dependence is detailed in Woo (2007), based on an empirical fit to spectral broadening measured from the Pioneer 10, 11 and Helios 1 space mission telemetry, and corroborated by radar reflections from Venus at superior conjunction (Harmon & Coles, 1983); this power law holds from ∼ 10R to ∼ 200R . Solar activity and coronal mass ejections can nonetheless cause significant variations. Fig. 6 (right) shows the estimated effect of the IPM from Eq. 4 for MWA frequencies. The effect of the IPM can be more significant than the ISM, particularly when the solar impact distance is small (i.e., the Sun passes close to the GC with respect to Earth).
Electron density fluctuations also impart scintillation on propagating signals, the characteristic timescale of which is given by If ∆t d is much shorter than the observation length, scintillation effects will be quenched; however, if ∆t d is much longer than the observation length, the signal may be effectively amplified or attenuated by scintillation during the observation, simplifying or confounding detection above our sensitivity threshold (Cordes & Lazio, 1991). Only the two nearest known exoplanets in the survey field, HD 164604 b and HD 165155 b, have ∆t d estimates significantly larger than our 7-hr observation length (using SM derived from NE2001). These sources could be reobserved at a later date, in case of scintillation-induced attenuation.
For our observations at 150 MHz, neither ISM nor IPM-induced broadening are significant effects. However, these effects should be considered in ∼Hz frequency resolution observations at low frequencies, in particular for proposed Moon-based SETI missions at frequencies below 1 MHz (e.g. Michaud et al., 2021) and future upgrades to the MWA. Signals from exoplanets that are occulted by their host star (with respect to Earth) can have large stellar impact distances and consequently large spectral broadening. The space weather around an exoplanet's host star will also affect broadening and scintillation; we note this as a potential field for future research.

CONCLUSIONS
We searched sightlines toward 144 known explanatory systems for artificial signals, finding no plausible technosignatures above an EIRP min from 10 13 W to 10 19 W (depending on distance to the source). We also conducted a blind survey across a 200 deg 2 field toward the GC, which covers billions of stars, including 3.3 million stars within 6000 pc 2 (0. (103-133 MHz) and higher-frequency (0.7-93 GHz) technosignature searches currently underway with the Parkes and Green Bank telescopes (Gajjar et al., 2021). Combined, these observations will allow statistical limits on the prevalence of putative transmitters at radio wavelengths in the GC. Future GC observations with the MWA and other instruments could allow full frequency coverage from 80 MHz to 93 GHz. Planned improvements to the MWA, due to be online in late 2021, will provide data with finer frequency and time resolution, which will improve the MWA's sensitivity to narrowband technosignatures and pulsating signals. Continual improvement of telescope capabilities, when combined with methodical observational campaigns, provides a means to explore the vast parameter space within which signs of technologically-capable life may be waiting to be found.

ACKNOWLEDGEMENTS
We thank A. Siemion and V. Gajjar for discussion on spectral broadening. We also thank the anonymous reviewers for their helpful comments and suggestions to improve the manuscript.

Facilities
This scientific work makes use of the Murchison Radioastronomy Observatory, operated by CSIRO. We acknowledge the Wajarri Yamatji people as the traditional owners of the Observatory site. Support for the operation of the MWA is provided by the Australian Government (NCRIS), under a contract to Curtin University administered by Astronomy Australia Limited.

Computer Services
We acknowledge the Pawsey Supercomputing Centre which is supported by the Western Australian and Australian Governments. Access to Pawsey Data Storage Services is governed by a Data Storage and Management Policy (DSMP). The All-Sky Virtual Observatory (ASVO) has received funding from the Australian Commonwealth Government through the National eResearch Collaboration Tools and Resources (NeCTAR) Project, the Australian National Data Service (ANDS), and the National Collaborative Research Infrastructure Strategy. This research has made use of NASA's Astrophysics Data System Bibliographic Services.

Software
The following software was used in the creation of the data cubes: