This site uses cookies. By continuing to use this site you agree to our use of cookies. To find out more, see our Privacy and Cookies policy.

The following article is Open access

A Search for Radio Technosignatures at the Solar Gravitational Lens Targeting Alpha Centauri

, , , , , , , , , , , , , , , , , , , , and

Published 2022 August 31 © 2022. The Author(s). Published by the American Astronomical Society.
, , Citation Nick Tusay et al 2022 AJ 164 116 DOI 10.3847/1538-3881/ac8358

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

1538-3881/164/3/116

Abstract

Stars provide an enormous gain for interstellar communications at their gravitational focus, perhaps as part of an interstellar network. If the Sun is part of such a network, there should be probes at the gravitational foci of nearby stars. If there are probes within the solar system connected to such a network, we might detect them by intercepting transmissions from relays at these foci. Here, we demonstrate a search across a wide bandwidth for interstellar communication relays beyond the Sun's innermost gravitational focus at 550 au using the Green Bank Telescope (GBT) and Breakthrough Listen (BL) backend. As a first target, we searched for a relay at the focus of the Alpha Centauri AB system while correcting for the parallax due to Earth's orbit around the Sun. We searched for radio signals directed at the inner solar system from such a source in the L and S bands. Our analysis, utilizing the turboSETI software developed by BL, did not detect any signal indicative of a non-human-made artificial origin. Further analysis excluded false negatives and signals from the nearby target HD 13908. Assuming a conservative gain of 103 in the L band and roughly 4 times that in the S band, a ∼1 m directed transmitter would be detectable by our search above 7 W at 550 au or 23 W at 1000 au in the L band, and above 2 W at 550 au or 7 W at 1000 au in the S band. Finally, we discuss the application of this method to other frequencies and targets.

Export citation and abstract BibTeX RIS

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

1. Introduction

The Search for Extraterrestrial Intelligence (SETI) often focuses on looking for signs of technology around other stars. This, however, is just one way to approach the search for ETI. Bracewell (1960) posited that "superior communities" would send material probes to neighboring star systems for the purpose of communication. If this is the case, we should also be searching our own solar system for evidence of ETI. This idea is the basis of artifact SETI, or solar system SETI (also called the Search for Extraterrestrial Artifacts (SETA); Freitas 1983), the search for any material artifacts that may be sent by ETI during a long-term galactic exploration effort. Freitas (1983) proposed that a search for such objects in our solar system could be carried out using telescopes, radar, infrared radiation, direct probes, and other methods.

1.1. The Interstellar Communication Network Hypothesis

Communications between probes and their home systems across interstellar distances would require extraordinarily high directional gains or transmitting powers to ensure that a distant receiving dish can reliably reconstruct information packets. Eshleman (1979) proposed that the Sun's gravitation could be used as a lens to magnify the radiation to or from a distant source. A spacecraft beyond 550 au along the focal line could use this effect to observe and communicate at interstellar distances using instruments comparable to those that we currently use for interplanetary communication. One proposed application of this mechanism is to utilize the Solar Gravitational Lens (SGL) to provide huge magnification to a telescope along the Sun's focal line at a distance between 550 and 1000 au (e.g., the FOCAL mission; Maccone 1994; Maccone & Piantá 1997; Maccone 2010).

Gillon (2014) proposed a new SETI approach to monitor the focal regions of the SGL corresponding to nearby stars to search for communicative technology. Their hypothesis is based on two assertions about galactic exploration. First, "Von Neumann" (self-replicating) probes (Tipler 1980) exploring the Galaxy would require some communication among probes and back to the original system. Second, in order to communicate across these interstellar distances, the explorers would create a network of probes leveraging stars as gravitational lenses. These assumptions solve the problem of long-distance communication across the galaxy by eliminating long-range communication in favor of a locally multinodal network connected to their nearest neighbors (Gertz 2017, 2020); this strategy is commonly used for communication networks around the Earth, such as the Internet. Gillon (2014) acknowledged that such devices would not likely be detectable through imaging or stellar occultation, but that multiwavelength monitoring may detect leaked signals.

Hippke (2020a, 2020b, 2021) explored technical constraints involved in the configuration of an interstellar communication network using stellar gravitational lenses, in order to motivate searches at the proper locations and sizes for nodes and probes. Their work argued that optimal probe sizes would be on the order of 1 m, the optimal communication wavelength range would be from 100 μm to 1 nm, and that a distance of 1000 au from the Sun is preferred. Hippke (2020b) also notes that, in this proposed configuration, separate exploration probes would be required in order to study the inner solar system. Kerby & Wright (2021) explore the engineering requirements and sustainability for an SGL relay to remain in proper position, finding that such technology is feasible and noting that another observable property of such probes may be the by-products of station-keeping propulsion. They also argued that single stars (i.e., those with minimal gravitational perturbations from companions) were the best hosts for such probes because their station-keeping costs were smallest.

Gillon et al. (2021) searched for a probe on the Solar focal line in communication with the Wolf 359 system, the third nearest star to our Sun and one from which Earth's transit across the Sun would be observable, at a time when Earth could have intercepted the transmissions. Their search for optical signals from the solar system toward the star, as well as for an object with the position and motion hypothesized within the extent of Uranus' orbit (20 au), did not reliably identify any such probe.

1.2. Solar Gravitational Lensing for Interstellar Communication

Massive objects in the universe, such as black holes and stars, bend the trajectories of nearby photons. This process can create a lensing effect similar to a focusing element of a telescope (Einstein 1936). Gravitational lensing warps a source object's apparent shape into two distorted images. The distortion effect centers on a ring around the lens object, called the Einstein ring, which has an angular radius of:

Equation (1)

where G is the gravitational constant, M is the mass of the lens object, c is the speed of light, DL is the observer–lens distance, DS is the observer–source distance, and DLS is the lens–source distance. In a solar gravitational lens system, the Sun is the lens, the distant star is the source, and the probe is our origin/observer. The physical Einstein ring radius (RE = DL θE ) of the solar lens should be at least the radius of the Sun. This works out to a minimum DL of roughly 550 au, which we adopt as our minimum possible probe distance. Figure 1 shows the layout for a probe utilizing the SGL to send or receive messages to/from α Cen. We note that Gillon et al. (2021) argue that the use of out-of-focus nodes at lower separations may also be worthwhile.

Figure 1.

Figure 1. Visualization of an SGL probe opposite α Cen. The Sun is at the origin in our coordinate system, and Earth's orbit lies on the xy plane. The probe is indicated in green and its outgoing beams for various dish sizes in gray (see Section 1.2).

Standard image High-resolution image

It may be possible to detect signals from a relay spacecraft at a distance d ≳ 550 au from the Sun if the Earth is contained within the opening angle of the probe's outgoing transmission beam at any point in its orbit. For target stars that lie along the ecliptic plane, the Earth will always pass through the beam when it transits and eclipses the Sun. However, from the view of a relay off of the ecliptic plane, Earth will always have some nonzero separation from the Sun, which its beam may not encapsulate. From the point of view of a probe at a distance of 550 au, the maximum angular separation between the Earth and the Sun is 6farcm3, regardless of its ecliptic latitude. A probe with an ecliptic latitude of b sees a minimum angular separation

Equation (2)

The Earth's entire orbit would be in a beam for 10 cm wavelength signals for any transmitter with diameter ≲34 m. However, Earth's point of closest approach to the SGL line may fit into smaller beams. This is a unique feature of performing this search search at radio wavelengths. In optical wavelengths, the focus of prior searches, the beam sizes for dishes on the order of 10 m would be much less than 6farcm3. Thus, eavesdropping on optical SGL signals would only be feasible for target stars in the ecliptic plane. Following Equation (2), for a probe opposite α Cen (b ∼ 44°), the minimum visible beam size at the point of closest approach would be 4farcm4. This corresponds to transmitter diameter ≲50 m.

The directional gain for a transmitting relay is

Equation (3)

where rt is the radius of the transmitter, kt is its efficiency (which we assume to be of order unity), and λ is the emitting wavelength (Maccone 2011). Hippke (2020b) suggests ∼1 m transmitter size for submillimeter signals, but a larger transmitter would be required to produce similar gains at radio wavelengths. A 50 m transmitter could achieve directional gains on the order of 65 dB, and even a 20 m transmitter could achieve 55 dB.

From Equations (8) to (9) of Maccone (2011), we find that the gain resulting from a solar lens is:

Equation (4)

For λ ∼ 10 cm, the solar lens gain is roughly 60 dB. Combining the directional gain of the relay probe and the lensing of the Sun, a transmitting spacecraft could achieve a total gain of Gtot = G × GSC of over 120 dB, overcoming the difficulties of transmission across interstellar distances by focusing transmissions into tight parallel beams.

1.3. Searching for Technology at the SGL

1.3.1. What to Look for

There have been many suggestions in the literature of purposes for probes in the solar system. The SGL foci of nearby stars are special locations for probes because they allow access to a "Galactic Internet," but the purpose and functions of such probes is otherwise unconstrained. In contrast to classical radio SETI searches for intentional messages sent toward Earth from interstellar distances, this method primarily attempts to eavesdrop on communications between two technological structures. We offer here three suggestions from the literature on such functions and the resulting implications for our search, but emphasize that our search is not dependent on and does not assume the probes we seek serve such functions.

Direct interception of outgoing communication through the SGL is the obvious signal to search for under the interstellar communication network hypothesis, but it suffers from some technical limitations. For interstellar transmissions from a relay spacecraft to be detectable from the Earth, the stellar relay must be actively transmitting when we search, such searches can only be performed when the Earth happens to be in the beam, and/or the beam could be tightly focused on a portion of the Einstein ring the Earth does not transit. A constantly transmitting stellar relay is more likely to be detected from Earth than a relay that only intermittently sends signals to the target star or local probes.

When searching for a relay, one major obstacle is determining which wavelengths would be most appropriate for SGL lensing. Radio wavelengths are an obvious choice for interstellar communication, due to their low energy and extinction. However, the solar corona can interfere with long wavelength radio signals. Turyshev & Toth (2019) found that photons with λ ≥ 157 cm (ν ≤ 0.2 GHz) are completely blocked. At shorter wavelengths, plasma refraction steers some photons away from the focal line, decreasing signal strength (Hippke 2020b). Moving from 3 cm to 200 μm (10 GHz to 1.5 THz), the effect gradually decreases. At wavelengths below ∼100 μm (ν ≳ 3 THz), the decrease in solar lens gain is negligible. For wavelengths nonnegligibly but not fully blocked, the SGL could still provide significant gains, so these are not completely eliminated from consideration, but do make searches for intercepted transmissions at radio or microwave frequencies less well motivated.

On the other hand, shorter wavelengths produce smaller beam widths, which would make eavesdropping on an SGL signal impossible for reasonable transmitter sizes and targets not near the ecliptic plane.

In addition to sending messages via the SGL, a relay spacecraft might send signals to other probes in the inner solar system. Indeed, such a scheme is recommended by Maccone & Antonietti (2022) as a way to retrieve information from interstellar probes, and is reminiscent of relay schemes used by humans' interplanetary probes, for instance for communication with Mars landers. There is no reason such secondary communications might not happen at different wavelengths than the interstellar communications through the SGL, as coronal interference is not a problem for them. It is thus possible they would occur at radio frequencies and with beam sizes and directions that would allow them to be intercepted at Earth year-round. Figure 2 shows the view of the solar system from a probe at the Sun's focal point for communication with α Cen.

Figure 2.

Figure 2. View of the solar system from a probe opposite α Cen at 550 au on the night of our observations. The Sun lies at the origin, and the inner solar system planets' orbits are shown. Gray filled circles indicate beam widths for various transmitter diameters at 10 cm wavelengths. For a probe located at a distance d > 550 au, these orbital axes would scale down by a factor of d/(550 au).

Standard image High-resolution image

Bracewell (1960) suggested that the purpose of a solar system probe would be to serve as a beacon, and Freitas (1983) and others considered places where such beacons would reside. Success in intentional communication SETI requires a determination of how to find another group of beings without being able to communicate beforehand. This requires the determination of optimal places, times, frequencies, and other aspects of possible communication which may be chosen by the other group. In game theory, these optimal choices are referred to as "focal points" (Schelling 1960), but we use the term "Schelling points" in order to avoid confusion with the optical definition of focal points (Wright 2020). The SGL focal line is a recognizable location for both humans and the residents of nearby stars, and it is a relatively static location in the sky. Thus, it may be an optimal place to place a beacon to intentionally attempt to make contact with technological life in the Solar system.

Beacons sent to other star systems may be a preferred method of interstellar first contact, as it does not immediately reveal the location of the sender, reducing risk of aggressive retaliation (Gertz 2018). An interstellar communication relay only reveals the position of the next node in the web. By their nature, beacons are intended to be found, so a beacon at the SGL might transmit at any frequency its constructors thought potential recipients would guess, further justifying our choice of frequencies near the "water hole" (Oliver 1979) for this study.

For this work, our first observations of the SGL for nearby stars as a proof-of-concept, we have chosen to observe in L and S bands. We chose these bands for three reasons: (1) the size of beams at these frequencies allowed us to check them all in a single 1 hr session; (2) these are in or near the "water hole" (to search for beacons); and (3) these are or are near frequencies humans use for interplanetary communication (to search for communications with probes in the solar system).

1.3.2. Where to Look

To first order, we expect this type of probe to lie at the antipode of some nearby star, around which another interstellar communication network probe could reside. Initially, we neglect light travel time considerations and place the probe at the exact antipode of the target star as we currently see it in the sky. This position on the sky, viewed from Earth, will vary slightly based on the probe's assumed distance from the Sun because, being a small fraction of a parsec distant, the probes suffer significant parallax.

For observation time t, we consider the positions of the Sun and Earth in the barycentric frame, S(t) and E(t), respectively. We also define unit vector x(t) as the direction of the target star from the Sun. We consider Sun-probe distances, z, from 550 au to infinity. We then calculate the barycentric frame position of the probe as:

Equation (5)

Finally, on-sky coordinates of the probe are calculated from the Earth-probe vector PE = PE(t).

Finite light travel time requires some corrections to this expected position (Seto & Kashiyama 2020). We must consider the positions of our celestial objects at other points in time. This calculation can be performed for two probe types: receivers and transmitters.

For a probe transmitting signals through the SGL to a star, the signal must travel to the position the star will be one Sun-star light travel time in the future. As our naive calculation uses the apparent position of the star at our observation time, which is where the star was one travel time ago, we must advance the star's position by two light travel times.

In principle, we must also consider that the probe must aim to where the Sun will be when the signal arrives, and advance its position by one Sun-probe light travel time. We neglect this effect because the maximum change it has on the position of the probe is of order 0farcs1.

We thus calculate the transmitting probe position as:

Equation (6)

where c is the speed of light.

For a probe near the Sun receiving signals from a star, our naive calculation of the star's position is appropriate, because the signal will arrive from the apparent (retarded) position of the star, just as other light from the star does. But for a distant probe, we need to consider that we see the probe where it was one Earth-probe light travel time ago, when it was sitting in a position corresponding to the star's position one probe-Sun travel time before that. Formally, we should then retard the star's apparent position by a time (z + ∣PE∣)/c; however for computational simplicity we approximate ∣PE∣ ≈ z because the difference in times is small compared to the timescale over which the star's apparent position changes significantly.

Our receiving probe position is then:

Equation (7)

This light travel time complication means that the probe position does not converge at infinite distance. Here, we assume that such a probe would not be placed further than a tenth of the Sun-star distance.

For both types of probe we have a locus across the sky of possible locations for a probe at varying distance z from the Sun. In general, the full extent of this line across the sky should be observed in order to detect or rule out active communication from a probe at the target star's SGL location.

1.4. The SGL for α Cen

In this work, we demonstrate the SGL relay SETI search method described above with a search for a hypothetical probe in contact with α Centauri. Figure 3 shows the positions of transmitting and receiving probes for α Cen A (the two components of α Cen are separated by less than 10'', which is so much smaller than our beam that our search encompasses both components).

The observations are described in Section 2. We present the analysis of our data in Section 3. In Section 4, we discuss the results of our analysis and additional artifact SETI searches of this type. We conclude in Section 5.

2.  α Centauri Observations with GBT

We selected a set of positions corresponding to the possible locations of a relay probe in communication with α Centauri, considering both the receiver position and the transmitter position. We took observations with the Green Bank Telescope (GBT) in the L and S bands using the Breakthrough Listen (BL) backend (MacMahon et al. 2018). As shown in Figure 3, we observed the region of the sky directly opposite α Cen placing observations along a line to account for the parallax of a probe at finite distances from the Earth and Sun. A probe at infinite distance from the Sun with no accounting for light travel time will be positioned exactly opposite α Cen on the sky and have no parallax, but a probe at 550 au will be offset from the exact antipode by up to 6 arcmin, depending on the position of Earth in its orbit. We trace two of these lines, as receiver and transmitter probes require different considerations of the light travel time between the Sun and α Cen, as described in Section 1.3.2. The GBT pointings were selected to cover regions along the lines of possible positions on the sky, leveraging the different beamwidth of the various bands in each case.

Figure 3.

Figure 3. Observations of the possible probe positions for communication with α Cen on UT 2021 November 6. The blue and orange lines mark the positions for receiving and transmitting probes, respectively, with light travel time taken into account. Tics along the probe position lines show where a probe at certain distances from the Sun would lie. The FWHM beam positions of each observation are represented by black circles, labeled by filter.

Standard image High-resolution image

We observed our target area on UT 2021 November 6. The Earth is at its closest position to the vector connecting the Sun to the target point in early November. In addition to minimizing the beam size required for the hypothesized signal to be detectable from Earth, this also minimizes the number of pointings required to cover the entire line of possible probe locations.

Data were taken at 300 s per scan using an ABABAB pattern. The ABABAB sequence consists of three "nods" between our "ON" source target and an "OFF" target. The OFF target is used to identify and correct for radio frequency interference (RFI). We chose HD 13908 as our OFF target because it has the smallest on-sky separation from the ON target of any known planet-hosting star. The separation between the ON and OFF targets is 317', which is equivalent to 37 beam widths at the L band and 58 at the S band. Additionally, a scan of a pulsar was completed at the beginning of each night of observations to confirm that the telescope and instrument were performing as expected.

3. Analysis

3.1. TurboSETI

We used turboSETI (Enriquez & Price 2019) to analyze our observations by searching for signals with specific drift rates. In the rest frame of the Sun, the proposed communication relay is effectively stationary, so a monochromatic transmission viewed from Earth should only show Doppler shifts due to the orbital motion and the rotation of the Earth. To obtain a drift rate for a signal from the relay, we used barycorrpy (Kanodia & Wright 2018) to calculate the radial acceleration between GBT and the probe. Given this radial acceleration dvr /dt, we calculate a maximum drift rate ${\dot{f}}_{\max }$ in Hz s−1:

Equation (8)

where ${f}_{\max }$ is the maximum rest-frame frequency observed.

The magnitude of the drift rate depends on the frequency and relative acceleration at the time of observation. The largest contribution to the relative acceleration term is due to the rotation of the Earth. A signal from a probe at rest relative to the Sun should exhibit a negative drift rate to an observer on Earth (Sheikh et al. 2019).

We calculate drift rates for the hypothetical probe relative to the barycentric frame at the time of observation based on the pointing of each beam. These drift rates are around −0.124 to −0.127 Hz s−1 in the L band, and −0.18 Hz s−1 in the S band. To account for uncertainty in the probe's position within a generous margin, we multiply the calculated drift rate for each pointing by a small numerical factor (×1.3) and use these modified drift rates as the maximum drift rates in our search through the data. Without specifying a minimum, we employ turboSETI to search for hits up to this maximum drift rate.

As the maximum drift rate of the hypothetical probe is quite small, it is computationally feasible to conduct a high-sensitivity search. By default, turboSETI searches for signals above a signal-to-noise (S/N) threshold of 25. For our analysis, we lower this threshold to 10.

Following an initial search for signals, turboSETI performs an additional check on the preliminary detections using three different "filter thresholds" to remove RFI from the pool of candidates based on the ON/OFF observing cadence. The first filter threshold simply looks for signals above the specified S/N. The second threshold looks for signals above the S/N that appear in at least one on-pointing, but not in any off-pointings. The third looks for signals that appear consistently in all ON pointings and no OFF pointings.

It has been noted by the BL team (Enriquez et al. 2017) and others (Margot et al. 2021) that turboSETI has two potential weaknesses. First, frequency binning causes a significant drop in sensitivity for drift rates above 0.16 Hz s−1. Fortunately, for these observations the expected drift rates are of the same order as the drift rate where this drop in sensitivity occurs. Thus, we do not expect this issue to cause significant sensitivity loss in this case. Second, it has been suggested that the filtering algorithm within turboSETI may be prone to missing signals that overlap with easily identifiable RFI.

However, the narrow range of possible drift rates predicted by our theory allows for a much more sensitive search (Wright 2020). We were able to inspect all of the signals within this range by eye (3864 hits in the L band and 448 hits in the S band) to ensure that no signals of interest were discarded along with RFI during filtering.

3.2.  L Band

The bandpass filter used on GBT in the L band has sensitivity over the range of about 1.07 to 1.87 GHz, with a notch filter from 1.25 to 1.35 GHz. The RFI environment for the L band is notoriously "noisy," and the notch filter eliminates some of the powerful RFI that frequently dominates the band. 12 The "hits" plotted in Figure 4 show the crowded RFI environment. Hits outside regions of detector sensitivity, including the notch filter, have been included in grayed out regions to provide a more complete picture of the data set and RFI environment.

Figure 4.

Figure 4. All of the signals above an S/N of 10 detected at the first filter threshold (any signal above the S/N detected in any ON source pointing) of turboSETI in the L band. Each signal is sized by the log of the S/N. The drift rates of each signal are determined by the turboSETI algorithm. Hits passing only the first filter threshold (any signal above the S/N detected in any ON source pointing) are marked as blue dots. The purple triangles indicate hits passing the second filter threshold (any signal above the S/N in at least one ON and no OFF source windows). The golden stars indicate the hits that passed the third filter threshold (signal detection in all ONs and no OFFs). The orange dashed lines correspond to the barycentric drift rates, of −0.124, −0.125, and −0.127 Hz s−1, for each of the three ON source pointings. The faint blue vertical stripes correspond to signals within 1.87 MHz (0.1% of the maximum frequency in this band) of each other, likely indicating they are part of an RFI comb. Data outside regions of detector sensitivity as well as those falling within the notch filter (1.25–1.35 GHz) have been grayed out.

Standard image High-resolution image

The abundance of signals detected with multiple drift rates within narrow ranges of frequencies is characteristic of RFI. We refer to signals with the same morphology spanning multiple frequencies as combs. Although turboSETI is designed to not re-trigger on the same signal, it often triggers at multiple drift rates on the same comb when that comb spans a large enough frequency range and its morphology is sufficiently complex. The addition of the light blue vertical lines in Figure 4 illustrates the presence of these combs by highlighting hits closely grouped in frequency, separated by less than 0.1% of the maximum frequency in the given filter. The orange dashed lines show the expected drift rates for the ON source pointings as discussed in Section 3.1.

From the entire L-band data set, turboSETI identified 27,995 hits passing the first filter threshold, 23,306 of which were within the frequency ranges the bandpass filter is sensitive to, as described above and shown in the white regions of Figures 4 and 5. 9929 hits passed the second filter threshold, 9585 of which were within the frequency ranges the bandpass filter is sensitive to, and only three signals made it through the third filter threshold, all within the sensitivity range.

Figure 5 isolates all the hits within the barycentric drift rate window shown in Figure 4. Only the drift rates for the ON source targets are relevant, so we considered only those drift rates calculated by turboSETI near the expected values. The marker types indicate the highest filter threshold that each hit passed. Out of the three hits passing the third filter threshold, only one was within the drift rate range we were searching for.

Figure 5.

Figure 5. All of the signals above an S/N of 10 detected using the turboSETI algorithm around the barycentric drift rates, between −0.124 and −0.127 Hz s−1 in the ON target pointings, for these observations in the L band. This shows the S/N for all the hits around the dashed orange lines shown in Figure 4 within the detector sensitivity range. Hits passing only the first filter threshold (any signal above the S/N detected in any ON source pointing) are marked as blue dots. The purple triangles indicate hits passing the second filter threshold (any signal above the S/N in at least one ON and no OFF source windows). The lone golden star indicates the one hit within this drift rate window that passed the third filter threshold (signal detection in all ONs and no OFFs). Figure 6 shows the dynamic spectrum of this signal.

Standard image High-resolution image

Figure 6 shows the output plot of the dynamical spectrum for this best candidate and more generally illustrates the ON/OFF observing cadence. Figure 7 zooms out to show this signal over a wider frequency range, which clearly identifies both what turboSETI likely erroneously triggered on, and that the signal is part of a ∼80 kHz wide sweeping RFI signal. While we cannot be certain about the signal's origin, its frequency, brightness, and drift structure are consistent with what would be expected from an IRIDIUM satellite downlinking in its L-band allocation (Maine et al. 1995).

Figure 6.

Figure 6. The single L-band hit that passed turboSETI's strictest filtering threshold, intended to select for narrowband signals that appear in every ON source pointing and none of the OFF source pointings, in the target drift rate range near −0.124 Hz s−1. Each rectangular block labeled "ON_L" and "OFF_HD13908" stacked on top of each other shows the ON/OFF observation cadence. This is plotted over a narrow frequency range of 602 Hz centered at the detected signal frequency of 1626.076447 MHz. The red dashed line indicates turboSETI's best-fit drift rate for this signal. The color mapping indicates power, normalized to the maximum power in each panel. It is difficult to determine what turboSETI picked up from this plot alone, but Figure 7 shows that it is RFI with a bandwidth of ∼80 kHz. Many hits in turboSETI show similar features, where a zoomed-out look clarifies that the signal is drifting RFI.

Standard image High-resolution image
Figure 7.

Figure 7. Zoomed-out plot of Figure 6, covering more than 80 kHz in the L band (as compared to 602 Hz in the previous plot). The red box indicates the plot window in Figure 6. From this plot it is easier to see what the algorithm likely picked up as a potential signal. This shows our most promising result is simply an RFI signal.

Standard image High-resolution image

To be thorough, we concatenated all of the hits in the data from the ON sources and manually selected them by drift rate, over the range of −0.110 to −0.127 Hz s−1. We plotted all 3864 of these signals and examined them by eye to confirm that turboSETI only filtered out RFI, finding nothing of interest.

3.3.  S Band

The bandpass filter used on GBT in the S band is sensitive over the range of 1.80–2.80 GHz, with a notch filter from 2.30 to 2.36 GHz (Lebofsky et al. 2019). Figure 8 provides a look at the RFI environment in the S-band observations with all of the hits passing the first filter threshold of turboSETI according to their drift rate and frequency. Similar to Figure 4, regions where the detector has little to no sensitivity are grayed out. The RFI environment in the S band is not as crowded as in the L band, but turboSETI still detects multiple frequency combs, again highlighted in vertical bands.

Figure 8.

Figure 8. All of the signals above an S/N of 10 detected at the first filter threshold (any signal above the S/N detected in any ON source pointing) of turboSETI in the S band. Each signal is sized by the log of the S/N. The drift rates of each plotted point, corresponding to a signal above the S/N, are determined by the turboSETI algorithm. Hits passing only the first filter threshold (any signal above the S/N detected in any ON source pointing) are marked as blue dots. The purple triangles indicate hits passing the second filter threshold (any signal above the S/N in at least one ON and no OFF source windows). The golden stars indicate the hits that passed the third filter threshold (signal detection in all ONs and no OFFs). The orange dashed line corresponds to the barycentric drift rate, of −0.18 Hz s−1, for each of the three ON source pointings. The light blue vertical strips correspond to signals within 2.8 MHz (0.1% of the maximum frequency in this band) of each other, likely indicating they are part of an RFI comb. Data outside regions of detector sensitivity as well as those falling within the notch filter (2.30–2.36 GHz) have been grayed out.

Standard image High-resolution image

Figure 9 isolates all the hits within the barycentric drift rate window shown in Figure 8. Only the drift rates for the ON source targets are relevant, so we considered only those drift rates calculated by turboSETI near the expected values. The marker types indicate the highest filter threshold that each hit passed. Out of the six hits passing the third filter threshold, none were within the drift rate range expected for an SGL probe.

Figure 9.

Figure 9. All of the signals above an S/N of 10 detected using the turboSETI algorithm around the barycentric drift rates, of −0.18 Hz s−1, for these observations in the S band. This shows the S/N for all the hits around the dashed orange lines shown in Figure 8 within the detector sensitivity range. Hits passing only the first filter threshold (any signal above the S/N detected in any ON source pointing) are marked as blue dots. The purple triangles indicate hits passing the second filter threshold (any signal above the S/N in at least one ON and no OFF source windows). No hits within this drift rate window passed the third filter threshold (signal detection in all ONs and no OFFs).

Standard image High-resolution image

From the entire S-band data set, turboSETI identified 24,298 hits passing the first filter threshold, 23,563 of which were within the frequency ranges the bandpass filter is sensitive to, as described above and shown in the white regions of Figures 8 and 9. 2226 hits passed the second filter threshold, 2169 of which were within the frequency ranges the bandpass filter is sensitive to, and only 6 signals made it through the third filter threshold, all within the sensitivity range.

To be thorough, we concatenated all of the hits in the data from the ON sources and manually selected them by drift rate, over the range of −0.16 to −0.18 Hz s−1. We plotted all 448 of these signals and examined them by eye to be sure no signals were being filtered out with the RFI by turboSETI. We found no signals of interest within the frequency ranges of this bandpass.

4. Results and Discussion

Our analysis reveals no signals near the expected drift rate for a stellar relay opposite α Cen in L or S bands. Although turboSETI output three signals in the L band and six signals in the S band that passed all of its filters, only one of these was found around the expected drift rate, and we determined that all candidate signals are RFI. This was confirmed by visually inspecting the waterfall plots produced by turboSETI, paying particular attention to signals around the expected drift rate for each pointing and examining the RFI environment for each band. Figure 6 shows the only one of these candidates identified with a drift rate around what would be expected. Figure 7 shows this to be RFI.

To address concerns that a positive signal could be filtered out with the noise, we examined all of the signals detected around the barycentric drift rate for each band without applying turboSETI filter thresholds. Although not practical for many radio SETI searches, our narrow range of drift rates provided the opportunity to scrutinize the data more closely. The result of this search provided additional confidence that there were no true positives filtered out by the turboSETI algorithm.

4.1. HD 13908

In addition to our SGL search, we used the OFF source pointings to observe a nearby star, HD 13908. We ran another analysis up to the default drift rate of 10 Hz s−1 on these observations in the L band and S band at an S/N of 10. As this was not our primary target of interest, we did not address the limitations of turboSETI for this additional target. We report for completeness that no convincing ETI signals were found in these bands using turboSETI during the HD 13908 observations.

In the L band, turboSETI detected 4932 hits passing the first filter threshold, 1085 passing the second filter threshold, and 7 hits passing the third filter threshold. None of the hits that passed the third filter threshold were of interest. In the S band, turboSETI detected 1733 hits passing the first filter threshold, 392 passing the second filter threshold, and no hits passing the third filter threshold.

4.2. Detection Sensitivity

Sensitivity for radio SETI observations is often measured with the equivalent isotropic radiated power (EIRP; Tarter2001), which is a function of the observing instrument, the S/N threshold, the observing time, and the distance to the transmitter:

Equation (9)

where d is the distance to the transmitter, S/N is the signal-to-noise ratio, Δν is the frequency resolution of the detector (where here we assume the transmitter bandwidth is narrower than or equal to the resolution of the detector), and t is the on-source integration time. The system equivalent flux density (SEFD) is specific to the observing instrument in a given band. The SEFD is reported as 10 Jy for the L-band receiver and 12 Jy for the S band at GBT. 13 In both bands, the BL backend achieves a frequency resolution Δ ν of 2.79 Hz. The integration time for each of our pointings was 300 s, and we analyzed the data with a minimum S/N of 10.

At the shortest theoretical distance of 550 au, the minimum detectable powers of an isotropic transmitter in the L band and S band are EIRPL = 6.6 kW and EIRPS = 8.1 kW respectively. At the 1000 au suggested by Hippke 2020b, EIRPL = 22 kW and EIRPS = 27 kW. However, we anticipate a probe utilizing the SGL to be a directed transmitter. Assuming a conservative gain of 103 in the L band and roughly 4 times that in the S band (see Section 1.2), a ∼1 m transmitter would be detectable by our search above 7 W at 550 au or 23 W at 1000 au in the L band, and above 2 W at 550 au or 7 W at 1000 au in the S band. These transmitter powers are comparable to cell phones, which emit at 0.6 or 3 W, as well as CB radios at 4 W. A microbroadcasting FM transmitter puts out 7 W, and Voyager I 23 W.

Using GBT to make observations affords studies of this kind some of the most sensitive radio reception in the world, allowing us to search for relatively low-power transmitters. Conducting searches for more powerful transmitters may feasibly be carried out by less sensitive facilities down to the scale of amateur hobbyists.

4.3. Future Applications

Although α Cen belongs to our nearest neighboring star system, it may not be the best star to study in the SGL network context. Kerby & Wright (2021) found that close binary star systems are unsuitable for stellar gravitational lens relays due to the high delta-v cost required to maintain alignment. Kerby & Wright 2021 instead suggest that the best candidate stars for hosting relay systems are without close companions, without large gas giant planets, relatively low in mass, and not spinning rapidly enough to deviate from sphericity. However, a significant number of stars do not fit this criteria, and the difficulties may often be overcome with some extra delta-v or settling for lesser but still significant stellar lens gains. Hippke (2021) created a first attempt at a prioritized list of targets for these searches, based on factors including the lack of massive companions and the presence of confirmed planets.

Another consideration in selecting stellar relay probe search targets is whether our observations would be made from within the probe's beam, though we have shown that for relatively small spacecraft in radio wavelengths this is not of much concern. Focusing on stars in Earth's transit zone (Heller & Pudritz 2016; Sheikh et al. 2020; Kaltenegger & Faherty 2021) would allow for observations guaranteed to be in such a beam when Earth is eclipsed by the Sun from the star's viewpoint, even for large optical laser transmitters.

We observed the SGL position opposite α Cen over a single night. Such a limited search may miss extant probes if they are not constantly transmitting. A thorough search for stellar relay probes in the Solar System would involve the monitoring of the antipodes of many stars over longer periods of time. Some progress on this may, in fact, be possible with serendipitous archival observations that happen to include the antipodes of nearby stars (see M. L. Palumbo et al. 2022, in preparation).

Our search was limited to the region of sky that would be occupied by probes at or past the Sun's focal point of 550 au. However, Gillon et al. (2021) propose that searching for closer probes may still be worthwhile. Gillon et al. propose that an array of 1 m lasers may be placed much nearer the Sun, because even at 10 au, the loss in gain is only a factor of 33. Such transmitters would be subject to stronger perturbations due to planets, but they would also receive much more solar irradiation for power.

As discussed in Section 1.3.1, wavelengths longer than 100 μm are subject to signal strength loss due to the solar corona. Hippke (2020b) suggests that the optimal wavelength range for SGL signals would be between 100 μm and 1 nm. We observed at relatively long radio wavelengths, 1–4 GHz or 30–7.5 cm. Here, losses due to the solar corona are not prohibitive but are nevertheless significant.

The difficulty in using turboSETI to analyze observations at shorter wavelengths is, as mentioned in Section 3.1, that the 1-to-1 binning ratio is 0.16 Hz s−1. The drift rate of a signal is a function of the frequency of that signal, so higher frequencies (shorter wavelengths) increase the drift rate beyond the 0.16 Hz s−1 limit, which causes a significant reduction in sensitivity (Margot et al. 2021). Sensitivity loss is a trade-off that may be worthwhile for the opportunity to cover more parameter space, especially when specific power levels are not well motivated. It is also possible to recover some of the lost sensitivity for higher drift rates through "frequency scrunching," wherein frequency bins are effectively summed together. The next generation dedoppler and hit search algorithm from the BL team is in production with the intention to improve search sensitivity.

5. Summary and Conclusion

In this work, we have described the history and physics behind the proposed use of stars as gravitational lenses to create an interstellar communication network, which motivates an artifact SETI search for such probes in our own solar system. We presented our method for searching for radio transmission from a probe in the solar system in communication with nearby stars and demonstrated this method on α Cen. Our analysis found no signals near the expected drift rate during our observations. Therefore, if such a probe exists, it is not transmitting constantly in the 1–3 GHz range. We also present statistics regarding the RFI environment in this region of sky in the L and S bands. In the future, both archival data and new observations can use similar methods to further search for relays that make use of the Sun's lensing ability for interstellar communication.

N.T. and M.J.H. are co-first authors of this manuscript. This paper is a result of the class project for the 2020 graduate course in SETI at Penn State. We thank the members of the 2021 undergraduate Penn State SETI course for participating in these observations, and the support of the Green Bank Observatory staff for supporting these classes' visit. J.T.W. acknowledges useful conversations with Michael Hippke on the physics of the SGL.

N.T.'s contribution to this material is based upon work supported by the National Science Foundation Graduate Research Fellowship under grant No. DGE1255832. M.L.P. acknowledges the support of the Penn State Academic Computing Fellowship. J.L. acknowledges the support of the NASA Astrobiology NfoLD grant #80NSSC18K1140 and the NASA Pennsylvania Space Grant Consortium (#80NSSC20M0097 and #NNX15AK06H). Disclaimer: the findings and conclusions do not necessarily reflect the views of NASA. S.Z.S. acknowledges that this material is based upon work supported by the National Science Foundation MPS-Ascend Postdoctoral Research Fellowship under grant No. 2138147.

The Penn State Extraterrestrial Intelligence Center and the Penn State Center for Exoplanets and Habitable Worlds are supported by the Pennsylvania State University and the Eberly College of Science. Computations for this research were performed on Penn State's Institute for Computational and Data Sciences' Roar supercomputer.

Breakthrough Listen is managed by the Breakthrough Initiatives, sponsored by the Breakthrough Prize Foundation. This material is based upon work supported by the Green Bank Observatory which is a major facility funded by the National Science Foundation operated by Associated Universities, Inc.

Software: barycorrpy (Kanodia & Wright 2018), turboSETI (Enriquez & Price 2019), Skyfield (Rhodes 2019), Astropy (Astropy Collaboration et al. 2013, 2018), Matplotlib (Hunter 2007), NumPy (Harris et al. 2020).

Footnotes

Please wait… references are loading.
10.3847/1538-3881/ac8358