Radio frequency emissions from dark-matter-candidate magnetized quark nuggets interacting with matter

Quark nuggets are theoretical objects composed of approximately equal numbers of up, down, and strange quarks. They are also called strangelets, nuclearites, AQNs, slets, Macros, and MQNs. Quark nuggets are a candidate for dark matter, which has been a mystery for decades despite constituting ~ 85% of the universe’s mass. Most previous models of quark nuggets have assumed no intrinsic magnetic field; however, Tatsumi found that quark nuggets may exist in magnetars as a ferromagnetic liquid with a magnetic field BS = 1012±1 T. We apply that result to quark nuggets, a dark-matter candidate consistent with the Standard Model, and report results of analytic calculations and simulations that show they spin up and emit electromagnetic radiation at ~ 104 to ~ 109 Hz after passage through planetary environments. The results depend strongly on the value of Bo, which is a parameter to guide and interpret observations. A proposed sensor system with three satellites at 51,000 km altitude illustrates the feasibility of using radio-frequency emissions to detect 0.003 to 1,600 MQNs, depending on Bo, during a 5 year mission.

1. self-magnetic field aggregates MQNs with baryon number A = 1 into MQNs with a broad mass distribution 4 that is characterized by the value of B o and typically has A between ~ 10 3 and 10 37 , 2. aggregation dominates decay by weak interaction so massive MQNs can form and remain magnetically stabilized in the early universe even though they have not been observed in particle accelerators 4 , 3. the self-magnetic field forms a magnetopause that strongly enhances the interaction cross section of a MQN with a surrounding plasma 16 that is primarily sustained by radiation and electron impact ionization in the high temperature plasma formed by normal matter stagnating against the magnetopause, 4. B o > 3 × 10 12 T is excluded 4 by the lack of observed deeply penetrating impacts that deposit > Megaton-TNT equivalent energy per km, so Tatsumi's range of B S is reduced to 1 × 10 11 T ≤ B o ≤ 3 × 10 12 T, 5. MQNs satisfy criteria for dark matter even with baryon's interacting with MQNs through the self-magnetic field 4 , and 6. the unexcluded range of B o , the low (~ 7 × 10 -22 kg m −3 ) density of local dark matter, net incident velocity of ~ 250 km/s, and the high average mass of MQNs constrain the flux of MQNs of all masses to between 10 -7 and 4 × 10 -15 m −2 y −1 sr −1 and constrain the flux for MQN masses > 1 kg to between 3 × 10 -14 and 2 × 10 -17 m −2 y −1 sr −1 , so very large area detectors or very long observation times are required 4 to detect MQNs.
In this paper, we investigate the possibility of detecting MQNs interacting with the largest accessible target: Earth, with its magnetosphere to 10 Earth radii. We show MQNs necessarily experience a net torque while passing through matter, can spin up to kHz to GHz frequencies, and emit sufficient narrow-band electromagnetic radiation that could be detected in additional tests of the MQN hypothesis for dark matter.
Detection is also complicated by the magnetopause plasma (hot ionized gas) shielding some radio-frequency (RF) emissions. Spacecraft re-entering the atmosphere cannot communicate with ground stations because the surrounding plasma shields the emissions until the spacecraft slows down. The same blackout effect prevents RF emissions from MQNs from being detected during the high-velocity interaction of MQMs with matter in the troposphere or ionosphere.
However, MQNs that exit the atmosphere can be detected by their narrow-band, time-varying RF radiation, with frequency equal to the rotation frequency. In addition, MQN detection rate should be strongly correlated with the direction of dark-matter flux into the detector aperture. That preferred direction is determined by Earth's velocity about the galactic center and, consequently, through the dark-matter halo. These characteristics should permit their detection and differentiation from background RF.

Results
In this section, we show (1) when MQNs interact with plasma, their rotational velocity increases as their translational velocity decreases, (2) MQNs passing through Earth's atmosphere on a fly-by trajectory produce sufficient RF emissions to be detected by satellites but not ground-based sensors, and (3) the frequency, RF power, and event rate depend strongly on the value of the B o parameter in Eq. (1).
Rotational spin-up. Plasma is an ionized state of matter with sufficient electron number density n e at temperature T e for the electrons and ions to behave as a quasi-neutral fluid. Quantitatively, that means there is at least one electron-ion pair in a spherical volume of radius equal to a Debye shielding length λ D , in which the permittivity of free space ε o = 8.854 × 10 -12 F/m, the Boltzmann constant k B = 1.38 × 10 -23 J/K, and electron charge e = 1.6 × 10 -19 C. The Debye length is also the distance over which the thermal energy of the electrons permits charge separation to occur and defines electrical quasi-neutrality of a plasma 36 versus an ensemble of electrons and ions which behave as particles. The distinction is important for MQNs. A plasma impacting the magnetic field of a MQN acts as a conducting fluid and compresses the magnetic field to form a magnetopause 16 . A co-moving electron-ion pair impacting the magnetic field will separate because the forces from their opposite charges deflect them in opposite directions and will create a local electric field between them. The combined force from their electric E and magnetic B fields causes them to move in the same E × B direction at the same www.nature.com/scientificreports/ velocity 36 and through the magnetic field, instead of strongly compressing it. In contrast, conducting plasma shorts out the electric field and prevents particle penetration. A MQN moving through a plasma experiences a greatly-enhanced slowing down force 16 through its magnetopause, which is the magnetic structure formed by particle pressure from a plasma stream balancing magnetic field pressure around a magnetic dipole. For example, the solar wind forms a magnetopause with Earth's magnetic field. Since the particles' mean free path for collisions is much larger than the Larmor radius in the magnetic field, the physics of Earth's magnetopause is collisionless and applicable to the very small-scale lengths of a quark-nugget's magnetopause.
As derived in Ref. 16 , the cross section σ m for momentum transfer by the magnetopause effect is for magnetopause radius r m , MQN radius r QN , MQN speed v, and mass density of surrounding matter ρ x . The total force F o exerted by the plasma on the quark nugget is approximately.
Equations (3) and (4) let us calculate the decelerating force on the MQN and its translational velocity during passage through matter 16 .
In this paper, we extend the dynamics to calculate the torque and rotational velocity of the MQN passing through matter. Papagiannis 37 showed that the solar wind, which has mass density ρ x ≈ 10 -20 kg/m 3 and velocity v ≈ 3.5 × 10 5 m/s, exerts a torque T (N⋅m) on Earth (radius r o = 6.37 × 10 6 m and magnetic field B o = 3. × 10 -5 T) as a function of the angle χ between the magnetic axis and the normal to both the magnetic axis and the direction of the solar wind. His semi-empirical result is expressed in MKS units as in which C 2 = 1,400 with units of Ns kg −0.5 m −1.5 T −1 . Papagiannis validated the expression for the angles χ within 0.61 radians of 0 and within 0.61 radians of π (i.e. -0.61 ≤ χ ≤ + 0.61, -3.14 ≤ χ ≤ -2.53, and + 2.53 ≤ χ ≤ + 3.14, as illustrated in Fig. 1). By symmetry of the magnetic field, the torque is 0 at χ = 0 and χ = ± π/2. Since Papagiannis was only considering Earth, he limited his calculations to within 0.61 radians of the normal to the plasma velocity. Rigorously reproducing and extending his computational results to the larger angles required for MQN rotation is beyond the scope of this paper. Therefore, we extend his results by observing the amplitude of the torque is symmetric about χ = -3π/4, -π/4, π/4 and 3π/4, as shown in Fig. 1 by solid blue lines, and approximate the rest of the torque function by extrapolation of the tan(χ) function in Eq. (5), as shown with dotted blue lines in Fig. 1. The resulting functions F χ and T, shown in Eq. (6), replace Eq. (5). Equations (2)-(7) were solved for the angular velocity versus time. The interaction produces a velocitydependent and angle-dependent torque that causes MQNs to oscillate initially about an equilibrium. Since the quark nugget slows down as it passes through ionized matter, the decreasing forward velocity decreases the torque with time, so the time-averaged torque in one half-cycle is greater than the opposing time-averaged torque in the next half-cycle. The amplitude of the oscillation necessarily grows, as shown in Fig. 2. Once the angular momentum is sufficient to give continuous rotation, the net torque continually accelerates the angular motion to produce a rapidly-rotating quark nugget. As shown in Fig. 2, MHz frequencies are quickly achieved even with a 0.1 kg quark nugget moving through 1 kg/m 3 density air at 250 km/s. For smaller or larger masses, the resulting angular acceleration and velocity are respectively larger or smaller.
Several other approximations for the torque in the intervals shown with dotted lines in Fig. 1 gave the same frequency within 5%, which is within the uncertainty of the magnetic field parameter B o . equilibrium frequency and radiated power. Rotating magnetic dipoles emit electromagnetic radiation in the far field with power per steradian 38  The spin-up process strongly depends on the details of the surrounding material mass and MQN velocity along the path of the MQN, the MQN mass, and surface magnetic field B o . In spite of these complexities, we find that the spin-up time is very much less than the MQN transit time through the region of highest torque, and the lower limit of final rotation frequency can be adequately estimated by assuming the energy gained per cycle equals the energy radiated per cycle: in which the torque T is given by Eq. (6) and the radiated power P is given by Eq. (9). Combining Eqs. (6)-(10) gives 2π ω 0 Tdt = 2πP ω 2 , www.nature.com/scientificreports/ In the simulations discussed below, we solve Eqs. (3) and (4) for position and velocity v(t) calculated along a trajectory through the atmosphere to the position of maximum density ρ x , where we solve Eq. (11) for the lower-limit to the maximum frequency ω max . Attenuation of Rf power by magnetopause plasma. The surrounding magnetopause plasma has a characteristic plasma frequency in which n e = the local electron number density, e = the electron charge, and m e = the electron mass. The characteristic e-fold length for attenuating the radiated power for frequency ω < ω pe is k −1 atten : In practice, the scale length 0.5c/ω pe is approximately 0.5% of the ion Larmor radius, which is approximately the minimum thickness of the magnetopause boundary. Therefore, magnetopause plasma strongly absorbs RF energy for frequencies less than the plasma frequency.
Since plasmas strongly absorb electromagnetic radiation with frequency less than the plasma frequency, which varies with solar activity, practical detection of MQNs by their RF emissions is limited to ≥ 0.03 MHz in the solar-wind plasma near Earth orbit, to ≥ 0.4 MHz in the magnetosphere, and to ≥ 40 MHz in the ionosphere. The equilibrium frequency of the MQN spin-up process and the high density of the troposphere means, in practice, all RF emissions in the troposphere are strongly shielded. Therefore, we will focus on detecting MQNs in the magnetosphere after they have transited Earth's atmosphere, as illustrated in Fig. 3.
The middle MQN trajectory in Fig. 3 represents a direct impact 16 . The bottom trajectory represents a MQN that is gravitationally captured and does not exit into the magnetopause. The top trajectory represents a MQN that spins up during transit and is detected by satellite S1. Satellite S2 would not detect these MQNs since they have not passed through sufficient matter to spin up. Therefore, appropriate differences in event rates as a function of satellite position would support detection of dark matter.
Trajectories of quark nuggets are from a preferred direction in Fig. 3 because the velocity of the solar system about the galactic center and through the halo of dark matter nearly dominates the random velocity of quark nuggets, as discussed in a subsequent section. Figure 3. Near Earth environment with MQN trajectories (red lines). Three concentric circles represent MQN interaction volume: highly ionized and low density magnetosphere and ionosphere (purple); weakly ionized or neutral, low-density troposphere (blue), and neutral, high-density planet (gray). Satellites S1, S2 and S3 are shown in orbits that let them monitor narrow-band RF emissions above background as described below. The vertical ring (light blue) is a detection-area element for simulations of MQNs interacting with magnetosphere, ionosphere and troposphere, as discussed below. www.nature.com/scientificreports/ Radiofrequency background from MQns distributed throughout the galaxy. Dark-matter is distributed throughout the galaxy 2,3 . Therefore, MQN dark matter in interstellar space might interact with the local plasma density and emit RF radiation that fills the universe over billions of years. The resulting RF might be detectable as a galactic background near Earth. Olbers' Paradox on why the night sky is dark 39 addresses the same phenomenon for photons from stars. Evaluating the potential for detecting and interpreting this radiation is complex. A detailed analysis is beyond the scope of this paper, which focuses on detection of near-Earth MQNs. However, the following preliminary analysis shows that the MQN hypothesis cannot be tested by observations of galactic RF background near Earth. The plasma density and temperature in interstellar space vary greatly: 1. a minimum of ~ 10 -4 to ~ 10 -2 particles/cm 3 in Hot Ionized Media (HIM) at T e ~ 10 6 to 10 7 K composing 20% to 70% of interstellar space, 2. to ~ 0.2 to 0.5 particles/cm 3 in Warm Ionized Media (WIM) at T e ~ 8,000 K composing 20-50% of space, 3. to higher mass densities represented by ~ 10 2 to 10 6 particles/cm 3 in molecular clouds at T e ~ 10 to 20 K composing < 1% of space 40 .
As discussed above, matter has to act like a fluid plasma (instead of an ensemble of isolated charged particles) to form a magnetopause. Quantitatively, the scale length λ D for charge separation in Eq. (2) has to be less than the magnetopause radius r m defined in Eq. (3), so For example, using the mid-range value of B o = 2.0 × 10 12 T and ρ QN = 1 × 10 18 kg/m 3 , the MQN mass has to be greater than 10 7 kg for a magnetopause to form in HIM. The equilibrium rotation frequency for such massive MQNs is much less than 30 MHz and the corresponding RF is absorbed by the solar-wind plasma near Earth.
The corresponding threshold MQN mass for forming a magnetopause in WIM is 10 kg. However, the higher density plasma in the WIM slows those MQNs, so their RF is still less than 30 kHz, so it is also absorbed in the solar-wind plasma near Earth.
In general, the larger mass density of the rest of interstellar media produce even lower frequency RF, which is absorbed even further from Earth. We analyzed these effects for the full range of interstellar plasma conditions 40 and MQN mass distributions 4 as a function of the B o parameter. We found wherever MQN dark matter can produce RF emissions, the emissions are either absorbed in the surrounding plasma or in the solar-wind plasma near Earth. The RF is not detectable near Earth and cannot be used to test the MQN hypothesis. Therefore, we focus on detecting MQNs passing very near Earth and radiating well above the 30-kHz cutoff frequency of the solar-wind plasma.
Discriminating MQn events from background. Theoretical profiles of dark matter halos are guided by astrophysical observations, which are consistent with a mass density of dark matter near Earth of about 7 × 10 -22 kg/m 3 ± 70% 3,14 . This extremely low mass density and the very broad mass distributions 4 mean that the flux of MQNs is very low. Detecting them over the full range of not-excluded values of B o requires interaction volumes of at least planetary size and requires a means of reliably subtracting background. The orientation of the sensed MQN flux with respect to the direction of inflowing dark matter provides one such opportunity.
Models differ in their degree of self-interaction. Computer simulations 41 covering a substantial portion of these models indicate dark matter occupies a halo within and around the galaxy's ordinary matter and has a Maxwellian-like, isotropic velocity distribution. Although simulations predict the velocity distribution of dark matter varies somewhat with the self-interaction model, the most probable, isotropic speed is ~ 220 km/s with a full-width-at-half-maximum of ~ 275 km/s.
The solar system moves through this high-speed dark-matter halo in its ~ 250 km/s motion about the galactic center. The direction of this motion is towards the star Vega, which has celestial coordinates: right ascension 18 h 36 m 56.33635 s, declination + 38° 47′ 01.2802″. In addition, Earth moves around the Sun at ~ 30 km/s. The vector sum of these two velocities gives the net velocity of Earth through dark matter, and the negative of this vector sum is the velocity of dark matter relative to Earth, shown in Fig. 4 for the position of Earth on the first day of each month.
A sensor on Earth should detect the most events per hour when it is sensitive to the flux of dark matter from the direction of Vega and much less when Earth shields the detector from the flux. A satellite in orbit about Earth would encounter a higher flux of MQNs when it is not shielded by Earth and when it is within range of MQNs that have transited through enough matter to spin up, as illustrated by S1 in Fig. 3.
The dark-matter velocity distribution has a streaming component U s relative to the sensor and an isotropic component U iso relative to the galactic center. The isotropic component smooths the transition between the directly exposed and Earth-shielded conditions. U iso can be adequately approximated for our purposes as the highest probability thermal speed U iso = 2kT M of a 3D Maxwellian velocity distribution of dark matter with mass M and with temperature T, where k is the Boltzman constant. Simulations 41 indicate that U s ≈ U iso . Therefore, we calculated the detection rate as a function of the sensor's orientation on Earth with respect to Vega and S = U s /U iso using the method developed by Cai et al. 42 . The results are shown in Fig. 5.   Fig. 3 is appropriately and systematically less than the event rate for S1, the result would be convincing evidence of having detected MQN dark matter.
including MQns from all directions. The flux in Fig. 5 is normalized to F θ=0 , the total number of events m −2 y −1 sr −1 for a detection surface facing directly into the streaming velocity U s , and was approximated from the mass distributions in Ref. 4 , assuming the mean incoming velocity U s = 2.5 × 10 5 m/s and assuming the effect of U iso ≠ 0 on the flux is negligible to first order for θ = 0. Assuming cylindrical symmetry in azimuthal direction φ and integrating the curve in   to a tangent to Earth, as illustrated by the three trajectories in Fig. 3. Satellites sense MQNs after they have transited the highest density matter along their trajectories and experienced the corresponding torque, as described in Eq. (6), to produce the maximum frequency and radiated power. After they pass into the magnetosphere, they can be detected since their emissions above ~ 0.4 MHz are no longer shielded by the higher-density plasma of the ionosphere. In our simulations, Earth's atmosphere is divided into increments ∆h of altitude h. MQN trajectories are characterized by their minimum altitude for 0 < h ≤ 9 r e , for r e = 6.378 × 10 6 m, which is one Earth radius. For each increment ∆h, test MQNs with masses consistent with the mass distributions of Ref. 4 are injected from the right in Fig. 3 with initial velocity in the x direction v x = -U s = -2.5 × 10 5 m/s. Their positions and velocities are calculated under the combined effects of gravity and magnetopause interaction. For each test particle, the maximum torque encountered in its trajectory, i.e. the torque in Eq. (6) at the maximum value of the product of total velocity and the square root of the local mass density, is used to calculate the equilibrium frequency from Eq. (11) and the corresponding RF power from Eq. (9).
Characteristic times τ up = ω max I mom /T max for spin up and τ down = 0.5 I mom ω max 2 /P max for spin down by radiation loss are calculated. Since we find τ up is much less than transit time through the atmosphere, ω max is the lower limit to the maximum frequency, from Eq. (11). The value of τ down helps determine the detectability of each representative MQN. Characteristic e-fold times for spin up and spin down are included in Supplementary Results: Representative Data Tables for Sensor Design. Values of τ down vary from a minimum of 2.4 × 10 3 s to a maximum of 1.9 × 10 6 s and provide adequate time for detection. The cylindrically symmetric cross sectional area A h associated with the altitude increment ∆h and altitude h above Earth radius r e is illustrated in Fig. 3

and is given by
MQNs with final velocity exceeding low-Earth orbital velocity ≥ 7,400 m/s escape the RF-absorbing ionosphere and will be recorded by a satellite-based sensor.
Atmospheric density as a function of altitude h = r -r e for MQNs at radius r and Earth radius r e were derived from the literature and fit with the following equations: For radius r below the magnetosphere 43 , i.e. 0 > r -r e > 2.873 × 10 5 m: and for radius r in the magnetosphere 44 , i.e. 2.873 × 10 5 m > r -r e > 10 r e : Mass density in Earth's magnetosphere depends strongly on solar activity and varies greatly. The data in Ref. 44 were averaged to produce Eq. (18), which should be adequate to estimate the annual event rate.
Our simulations show MQNs radiating between 10 -28 W and 10 +13 W and at frequencies between 0.35 MHz and ~ 2 GHz. Very high frequencies are associated with negligible RF power, and very high powers are associated with frequencies that are shielded by the magnetopause plasma. Results for RF power as a function of frequency are shown in Fig. 6 for events radiating at more than 1 μW and at frequencies more than 100 kHz. Results are shown for four representative and non-excluded values of B o .
As shown in Fig. 6, highest power emissions occur at the lowest frequencies and highest values of B o . These originate from the most massive MQNs penetrating the troposphere, but there are very few of them. The map associated with B o = 1.5 × 10 12 T is common to the maps of all B o values. The differences represent aggregation run-away as discussed in Ref. 4 . Figure 6 represents events by frequency and RF power but does not indicate the expected number of events per year. For each test MQN, the effective target area, given by Eq. (16), was multiplied by the corresponding (16) A h = 2π(r e + h)(�h). Detailed results for MQNs with RF power greater than 1 nW and with sufficient flux to be in the most probable 80% of events are provided in Supplementary Results: Representative Data Tables for Sensor Design. The information should be useful for designing sensors for detecting MQNs.
Baseline sensor system. A realistic sensor system is, of course, essential for testing the MQN dark-matter hypothesis. A convenient coincidence of frequency range and emerging technology enable a practical sensor. We outline one such baseline system and compute the corresponding event rate in this section.
As shown in Fig. 3, the system has three satellites equally spaced in a circular orbit with inclination 38.783° (to match MQN flux) at 51,000 km altitude where (1) the background plasma density is sufficiently low to permit good RF propagation in the intended detection band of 10 5 to 10 6 Hz, (2) radiation damage from electrons in the outer Van Allen belt is minimized, and (3) coverage by three satellites is acceptable. The Interplanetary Monitoring Platform IMP-6 (Explorer 43) spacecraft 45,46 is the reference architecture for the sensor system's spacecraft. IMP-6 was a 16-sided drum, 1.8-m long by 1.35-m diameter, having four 46-m long monopole antennas operating in pairs as 91 m long dipoles, and spinning at 5.4 revolutions per minute to scan space. Our baseline design is the same architecture but with 350 m long dipoles.
The sensor is based on the standard radar equation, formulated as transmit-receive equation: where P r is the received power, f is the frequency in Hz, c is the speed of light, G r is the antenna gain, P s is the source power, R is the distance, and G s is the source gain, which will be taken to be unity to represent the time averaged value. Lossless dipoles up to a half-wavelength long have gains G r < 1.6 and directivities D r from 1.5 to 1.6. The weak signals require amplification in the receiver. The noise power P N and the signal to noise ratio S/N for an amplified receiver are www.nature.com/scientificreports/ where Boltzmann's constant k B = 1.38 × 10 -23 J/K, T noise is the absolute radiation noise temperature, T a is the preamplifier noise temperature, and △ f is the receiver resolution bandwidth. Solving Eq. (21) for the range R gives Since many measurements will be made on each MQN as it transits through detection range, we can operate at signal/noise ratio S/N = 1 and use signal averaging and pattern recognition to reliably detect the signal.
The background noise temperature T noise is a major factor in determining sensor performance. IMP-6 measured the galactic background noise in the magnetosphere between 354 km and 206,000 km altitude and between 130 and 2,600 kHz. To determine noise temperatures from the noise power measurements of Brown 45 , we equate his values for spectral brightness B to the brightness in the Rayleigh-Jeans law for blackbody radiation and solve for the noise temperature T noise The results are shown in Table 1. The frequency range of interest encompasses a temperature maximum of 24 million degrees (MK) near 0.8 MHz. Ground based measurements at the poles during solar minima by Cane 47 are also shown in Table 1 for comparison and are consistent with our analysis of Brown's 45 data.
Researchers accustomed to working at much higher frequencies, e.g. in satellite communications or radio astronomy, may find these values of T noise to be unreasonably high. However, T noise is decreasing from 0.8 MHz to 3.0 MHz in Table 1; it continues to decrease with increasing frequency to agree with T noise observed in those disciplines.
Fitting the data with the minimum of two polynomials reproduces the table for 0.2 MHz ≤ f ≤ 3.0 MHz within ± 1% and provides a useful function for estimating sensor performance: For these low frequencies f (MHz), the noise temperature is millions of degrees Kelvin. From the denominator under the square-root term in Eq. (22), we see that variation of range with antenna gain is small if Since noise temperature T a can be less than 100 K (which does not require a cryogenic amplifier) for the frequencies of interest, 25 MK > T noise > 3 MK in Table 1, and D r ~ 1.5, Eq. (23) gives 6 × 10 -6 < G r < 5 × 10 -5 before the range becomes very sensitive to antenna gain. For source powers of 1·w, antenna gains of 0.001 to 1, radiation noise temperatures of 10 million K for the present case, amplifier noise temperatures of 100 K, receiver bandwidths of 1 Hz, and signal to noise ratios of unity, one gets ranges of 14,954 km and 15,740 km at 0.5 MHz. Therefore, range is very insensitive to antenna gain, and antenna gains that are well below unity are acceptable when the radiation temperature is very high.
The antenna gain G r depends on the dimensions and the termination impedance of the dipole antenna. Tang, Tieng, and Guna's theory 48 was used to design the antenna subsystem. The antenna gain G r as a function of frequency is shown in Fig. 8 for four terminating impedances.
For a wide range of terminating impedances, G r easily satisfies Eq. (25) and the simple dipole is suitable for the baseline system. The optimal dipole design will likely incorporate resistive and/or possibly reactive loading at discrete positions along its length to achieve the best bandwidth and gain response if the additional complexity does not preclude deployment in space.
The other variables in Eq. (22) are source gain G s and signal-to-noise ratio (S/N). Both are set to 1.0. For multiple measurements like the satellites will be recording, detections with S/N = 1 are routine, and detections have been demonstrated for S/N as low as 0.1.
We found two options that apparently meet key requirements for the receiver: Δf = 1 Hz and a million simultaneously recorded frequencies. Assuming that such compact and low-power, ground-based electronic hardware can be made into space-qualified hardware, they provide proof-of-principle of receiver feasibility. The Tektronix RSA 306B Real Time Spectrum Analyzer system with Tektronix proprietary software on a separate computer  www.nature.com/scientificreports/ monitors 10 6 frequencies every 2 s with the required bandwidth Δf = 1 Hz and outputs the results to a laptop computer. In addition, the N210 USRP unit from National Instruments, which we have been using for autonomous data collection with our custom software for 254 frequencies every 25 μs, can meet the requirements by replacing the Field Programmable Gate Arrays (FPGAs) currently used for processing fast Fourier transforms with General Purpose Computing on Graphics Processing Units (GPGPUs). In both cases, processing the data on 10 6 simultaneous frequencies to find the unique signature of an MQN and extract the desired information will be a significant challenge. A combination of on-board processing with GPGPUs to select data for download and ground-based post-processing should be sufficient but has yet to be demonstrated. estimated MQn detection rates. Although the non-excluded range for B o at the time of this publication is 1 × 10 11 T ≤ B o ≤ 3 × 10 12 T, recent results being prepared for publication reduce the non-excluded range to 1.5 × 10 12 T ≤ B o ≤ 3 × 10 12 T. If B o < 1.5 × 10 12 T is included, the predicted MQN detection rate will be even larger and will shift to higher frequencies.
The number of events expected with the three- Since quark nuggets are also baryons, their magnetic field may be similar to that of protons, which have the equivalent magnetic field B o between 0.9 × 10 12 T and 2.2 × 10 12 T.
The case for B o = 1.5 × 10 12 T is quite different because the corresponding MQN mass is much smaller, so the flux is much larger and the frequencies are much higher. Approximately 7% of events are between 0.1 and 1.0 MHz, 26% are between 1 and 3 MHz, 54% are between 3 and 10 MHz, and 13% are from 10 to 200 MHz. The system should be optimized for 0.1 to 1.0 MHz. If the RF background permits some additional measurements between 1 and 200 MHz would explore the special case of B o ~ × 10 12 T, unless additional observations in the near future exclude that case. obtainable measurements and information. MQNs transiting through the magnetosphere, as shown in Fig. 3, provides a frequency source f 0 (t) moving at constant non-relativistic velocity v in a straight line. A stationary observer at distance r(t) records the observed frequency versus time. The closest approach occurs at time t 0 and at distance r 0 . The observed frequency f (t) will be Total Doppler shifts are expected to be about that 1.5% of f 0 . The source frequency is expected to decrease with time. During the brief periods of detectability, the source frequency should change linearly with time by amounts comparable to or less than the Doppler shifts. www.nature.com/scientificreports/ Figure 9 shows what can be extracted from measurements of frequency versus time from the moving source. The velocity v 0 is determined from the asymptotic limit lines. The slope of these lines determines the rate of change of frequency of the source. The location of the inflection point determines f 0 and t 0 , and the slope of the central asymptote determines (v 0 /r 0 )f 0 -df 0 /dt from which r 0 is determined with the other measured variables.
Because signal strengths are not likely to be measurable in the far asymptotic regions, it is important that good data be observed in the regions of maximum d 2 f/dt 2 , which occur at distances less than 2r 0 . Optimal fits of data to Eq. (26) will provide estimates of the variables of interest.
Thus, from the measurement it is possible to estimate t 0 , f 0 , df 0 /dt, r 0 , v 0 , and dv 0 /dt and analyze the data to get two measurements of MQN mass, one from f 0 and df 0 /dt and one from v 0 and dv 0 /dt. With enough data the MQN mass distribution can be constructed and compared with the predictions in Ref. 4 and be used to refine Tatsumi's value of the surface magnetic field in magnetar cores.

Discussion
We have shown that MQNs transiting through ionized matter experience a torque on their magnetopause and spin up to high frequencies, ranging from kHz to GHz, depending on MQN mass and velocity, mass density of the surrounding matter, and the B o parameter. Rotating MQNs radiate. If the radiation is above the plasma frequency of the surrounding matter, the RF radiation will propagate and can, in principle, be used to detect MQN dark matter at a substantial distance.
We did examine the possibility of using ground-based sensors to detect MQNs transiting the magnetosphere and radiating at > 40 MHz, the cut-off frequency of the ionosphere. The event rate is at most 0.3, 0.02, 0.0004, and 0.00007 per year even if 2π steradians solid angle can be observed. A space-based system is more promising.
Our results have identified requirements for MQN detection systems and a baseline system of three satellites that should test the MQN hypothesis for dark matter by detecting between 1,600 and 0.003 MQNs for B o between 1.5 × 10 12 T and 3.0 × 10 12 T, respectively. Recording their Doppler shifted frequencies during their passage by a sensor provides two different measurements of the MQN mass, so a mass distribution can be obtained in time.
The pattern of trajectories with respect to the direction of Vega, frequency, and rate of change of frequency would provide strong evidence of MQNs and eventually characterize the actual mass distribution of MQNs to compare with the predictions of Ref. 4 . Very low mass MQNs will have frequencies of up to 2 GHz, but their negligible RF power precludes their being detected. Very high mass MQNs will emit RF power of up to 10 13 W, but their negligible flux makes their detection unlikely in a year of observation. As shown in Supplementary Results: Representative Data Tables for Sensor Design, frequencies between 100 and 800 kHz are the most likely frequencies to be observed from MQNs.
Discerning quark-nugget signals from all human-caused and naturally occurring, non-MQN RF is facilitated by their characteristics: (1) narrow-band RF, unlike RF from lightning and other discharges, (2) continuous emissions, unlike pulsing radars, (3) relatively unmodulated in frequency and amplitude, unlike communication RF, (4) moving at ~ 200 km/s, unlike all human sources, and (5) initially increasing in frequency to a maximum and then slowly decreasing, unlike magnetosonic waves.
Finally, we note that solar and planetary atmospheres can be even larger-area, but less accessible, targets for RF-emitting MQNs.

Data availability
All final analyzed data generated during this study are included in this published article with its Supplements.
Received: 10 May 2020; Accepted: 3 August 2020 Figure 9. Observed Doppler shifted frequency versus time for source radiating at f 0 = unshifted frequency, moving with v 0 = relative velocity, passing at r 0 = distance of closest approach at time t 0 . Dotted lines show asymptotic rate of change of frequency df 0 /dt well before and after closest approach.