Exploring nine simultaneously occurring transients on April 12th 1950

Nine point sources appeared within half an hour on a region within \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼ 10 arcmin of a red-sensitive photographic plate taken in April 1950 as part of the historic Palomar Sky Survey. All nine sources are absent on both previous and later photographic images, and absent in modern surveys with CCD detectors which go several magnitudes deeper. We present deep CCD images with the 10.4-m Gran Telescopio Canarias, reaching brightness \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r \sim 26$$\end{document}r∼26 mag, that reveal possible optical counterparts, although these counterparts could equally well be just chance projections. The incidence of transients in the investigated photographic plate is far higher than expected from known detection rates of optical counterparts to e.g. flaring dwarf stars, Fast Radio Bursts, Gamma Ray Bursts or microlensing events. One possible explanation is that the plates have been subjected to an unknown type of contamination producing mainly point sources with of varying intensities along with some mechanism of concentration within a radius of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼ 10 arcmin on the plate. If contamination as an explanation can be fully excluded, another possibility is fast (t \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<0.5$$\end{document}<0.5 s) solar reflections from objects near geosynchronous orbits. An alternative route to confirm the latter scenario is by looking for images from the First Palomar Sky Survey where multiple transients follow a line.

www.nature.com/scientificreports/ the POSS-I E red plates. Most of them could only be found in one image. Since then, we have further examined the VASCO candidates and find that several images show at least two or three transients. In this Scientific Report, we report a special case where we find nine objects present at the POSS-I red plate from April 12th 1950 at the Palomar Observatory, which reaches a depth of r ∼ 20 mag. The plate has plate number XE 325 in the DSS Plate Finder. Whether the effect is a real astrophysical phenomenon or some type of instrumental artifact, the effect is crucial to understand the remaining ∼ 100 red transients, as some of them may be caused by similar phenomena. On one of the POSS red plates (XE 325), eight star-like objects and one with possible elongation, appear within a 10 * 10 arcmin 2 field-of-view centered at (RA,DEC) = (212.9291,+26.8311) (eq. J2000.0). We shall be referring to these objects as "nine simultaneous transients", despite that we are unsure on whether the object with slight elongation has the same origin as the others. The sources are not seen in the blue exposure 30 minutes earlier, nor in a subsequent red plate (XE 324) obtained six days later (see Fig. 1), nor in any subsequent archival images, including the POSS-II red imaging survey. Even in Pan-STARRS (limiting magnitude r ∼ 23.4) and SDSS (limiting magnitude r ∼ 23) nothing is seen at the positions of the transients. We investigate possible causes for the nine transients and present new, deep observations of the same field-of-view.

Method and preliminaries
In order to extract the coordinates, FWHMs and r magnitude of each object on the plate, we ran the photometric procedure described by Andruk & Villarroel et al. (to be submitted). This procedure builds on the methods of Andruk et al. (1995Andruk et al. ( , 2017Andruk et al. ( , 2019 [7][8][9] , and is adapted for POSS-I red emulsions specifically. We also used the procedure to check for other strong transients on the same plate (see Supplementary Information).
In order to gain a better understanding of what we see, we performed observations of the given field in the r and g filters with the GTC 10.4 meter telescope and the OSIRIS instrument on the night from the 25th to the 26th of May 2020. The astrometry was improved with Astrometry.net. Further, we used the Terapix swarp procedure to improve the astrometry, using SDSS images as a reference field for zero-point calculations in the astrometrical solutions. Exact details are in the Supp. Info, Sections A.1 and A.2. Out of the 9 sources marked in Fig. 1, only 6 fall into the GTC field of view -the two lower right ones fall outside the field of view, as well as the most upper left transient. We WCS-matched the POSS-I and astrometry-corrected GTC images.

Results
We manually identify the GTC counterparts to the transients using DS9. For each given transient, we look for counterparts residing within a circle with 5", 3" and 1.7" radius. We perform photometry for these counterparts (Supp. Info, Section A.4). Table 1 summarizes the POSS-I and GTC photometry for the six sources covered by the GTC/OSIRIS observations. The transients and their counterparts are presented in Figs. 2 and 3. Proper motion may make it difficult to be sure if one has assigned the right counterpart to the right object. For example, in Fig. 2 we see four transients, including the slightly elongated one (nr 5). All four have a counterpart within 5 arcseconds, but only three of these are easily distinguishable (the fourth being fainter). Only one of the counterparts are within a radius of 1.7 arcseconds. In the POSS-1 image we see a number of objects that cannot be subsequently found, marked with green circles. Purple circles are artifacts during the scanning process. About 9 objects are present in the POSS-I E image (left) from the 12th of April, but not in the POSS-2 image (right) from 1996. One slightly larger circle host two transients. In addition, the 9 objects are neither visible in the blue POSS-1 taken half an hour earlier, nor in a second POSS-1 red image taken six days later on April 18th (see additional images in Supp. Info. Fig. 2.) The 9 transients are not caused by a difference in depth or spectral sensitivity. The images are based on the DSS digitizations of the Palomar plates. www.nature.com/scientificreports/ We find the following: within a 5 arcseconds radius, 6 out of 6 transients have at least one counterpart. Within a radius of 3 arcseconds, 3 transients have one or more possible counterparts. Within a radius of 1.7 arcseconds (corresponding to the POSS-1 red images pixel size), 2 transients have counterparts. Given the density of objects within the given GTC field, we need to know whether the probability of finding this many counterparts is significant or not. We used a Sextractor procedure with a threshold of 3 σ (corresponding to the visual detectability of the counterparts), discounting objects right at the edge of the plate, to calculate the number of objects (2157) within the field of view of the red image (CCD1), that is 241 × 515 arcsec 2 in size. The average density is thus ρ = 0.0173 arcsec −2 . To calculate the probability of finding at least one counterpart within an area of A, Hence, the corrected probabilities to have at least one object within 1.7", 3.0", 5.0" are For calculating the probability of having at least m counterparts in n transients, with a probability P 1 , we applied the binomial distribu- 169 , which cannot reject the hypothesis of pure chance overlaps. It is not unthinkable, that eventual counterparts to the transients may have moved outside the 5 arcsecond circles during 70 years.
We then examined the colours of the possible counterparts, listed in Table 1. M dwarfs, which are generally very red, usually have significantly redder colours g − r > 1.5 . The colours of the possible counterparts are oppositely bluer with g − r < 1 . This excludes the M dwarf hypothesis for the simultaneous transients, and indicates that the counterparts are either star-forming galaxies or F, G or K main sequence stars.

Discussion
A number of instrumental issues should be commented upon. The glass cover used during the scanning process can produce many small, false stars. Adding the independently scanned SuperCosmos digitization, we can identify the artifacts in the image (see Supp. Info. A.7.2). Ghost images do not leave a PSF shape on the photographic plate. One may further speculate about the observers having exposed the same plate twice with a partially or completely open shutter. If so, one would expect a high density of extra stars all over the plate.
Non-astronomical effects apart from instrumental errors may arise from elementary particles. We can exclude muons and natural radioactivity, as these particles hit the photographic plates at a steady rate, which disagrees with a few plates having extra stars (had muons and natural radioactivity caused extra stars, the photographic plates would be full of them). We discuss cosmic rays and rare cosmic ray events in detail in Section A.7.5. (Supp. Info.). On the other hand, alpha, beta or gamma particles can certainly leave traces in a photographic plate 12-14 . An elementary particle entering the emulsion at some incoming angle will expose silver halide crystals within a cylinder that is slanted by that same angle. That slanted cylinder projects to the surface of the plate as an elongated oval shape. Such particles that hit the 103aE emulsion at 30 degrees can, in principle, result in an elongation that is so small (30.7 microns) compared to the FWHM of the seeing PSF that the image might appear nearly symmetric and "star-like" (see Supp Info, Section A.7.6).
For example, it has been shown 15 that the detonation of an atomic bomb in Alamogordo, New Mexico in July 1945 left radioactive contaminants due to the isotope Ce-141, seen as fogged spots, on photographic sensitive x-ray film produced by the Eastman Kodak Company. It was the same company that produced the 103aE and 103aO plates during POSS-I. The radioactive contaminants were transported to the Kodak factory by the wind. Following this incident Kodak developed containers that protected the photographic material from radioactive particles. As Mount Palomar is not too far away from New Mexico, it is not impossible that particles from atomic bomb tests contaminated some of the photographic plates, either while storing them or during observation.
A human source may be to blame. Small local events in an observatory have been known give rise to anomalous detections, e.g. the potassium flare stars 16 that were spectra contaminated by French matches, the discovery of unusual radio bursts, that were caused by microwave ovens 17 . In a recent search for laser signals from Proxima Table 1. Coordinates of simultaneous transients. We display the coordinates and magnitudes of the multiple transients that are covered by the GTC/OSIRIS images, including their counterparts. We see that 6 out of the originally 9 points are covered by the images. When two stars, like object 2 and 3, are too close to each other, the photometry code (that calculates Gaia r magnitudes) fails to detect one of the transients in the POSS-I image, which means we have no photometric information. These are marked with a star (*). Removing outliers with FWHMs larger than 10 pixel in the sample, the average seeing is about 4.8 + / − 0.8 pixel for all objects in the field. Finally, we show the g and r magnitudes (AB) for the closest counterpart.  www.nature.com/scientificreports/ plate. An observer suffering from hay fever during the spring time in California, might have produced aerosols that expose the photographic plates. But sprayed microdroplets of saliva, peaking around ∼ 100 micron in size, do not produce circular shapes on a flat surface, but elongated, wobbly and clustered driplets 19 , in sharp contrast with the point sources we see. Pure water driplets can be of sizes ranging from a few tens of microns to almost a millimeter in size, but also have a broader distribution of sizes than what is found in the POSS-1 plates. But maybe, another type of unanticipated contamination or degradation that coincidentally only creates objects with star-like brightness profiles is the guilty party. We investigate the other cases of multiple transients among the ∼ 100 VASCO transients and note that some very faint transients can occur in apparently contaminated regions, see Supplementary Information A.7.7. The glass cover during the plate scanning process is a possible source of contamination. For the more luminous transients among the simultaneous transients, the point sources are more difficult to reconcile with the hypothesis of contamination. The best way to exclude the possibility of contamination causing the simultaneous transients is by examining the original photographic plates with a microscope 11 . Unfortunately, we have no access to the original POSS-I plates. We next examine the possibility of naturally occurring transients. The total density of naturally occurring fast transients such as flaring stars, GRBs, or FRBs is estimated as about 1.63 events per square degree per day, up to a limiting magnitude of 23 or 24.7 mag in the g band 20 . The Large Synoptic Survey Telescope expects to find 100,000 short-lived transients per day, which corresponds to about 2.5 events per square degree per day 21 . Given that we have 9 transients within 30 minutes inside a 10 × 10 arcmin 2 box, no well-known natural transients can explain this phenomenon.
Another possibility are fragmenting asteroids. However, the field is far away from the ecliptic. Could the fragmentation of an asteroid sufficiently far from Earth result in a fragmented shower? Let us, for example, consider an asteroid at the main belt at ∼ 1.5 AU from the Earth and 2.5 AU from the Sun. The Earths motion ( ∼29.9 km/s) dominates over the motion of asteroids. During this time, the fragments must have moved at least 50 * 60 * 29.9 kilometers during the 50 minutes of exposure, resulting in 82 arcseconds long streaks in the image. At a distance of 30 AU (Kuiper belt), the displacement of fragments on the image is around 4 arcsec. At a distance on 0.1 AU from Earth and 1 AU from the Sun, the fragments would have moved roughly 20 arc minutes accross www.nature.com/scientificreports/ the field. Given the constant illumination of the Sun, long streaks are expected due to the orbital motion. Also, had a fragmentation occured, we' d see a parent body in the blue exposure taken 30 minutes earlier. We see no traces of a parent body in the previous blue exposure, taken half an hour earlier.
But, what if faint and small meteoride fragments burn up right upon atmospheric entry and shortly "glint" as they catch fire? A dust grain as small as 1 gram, can reach visible peak magnitudes around 0 -2 mag, as it burns up during the course of about 10 seconds. An object entering the Earth's atmosphere at 100 km can have a wide range of velocities. If the object hits the Earth directly "face-on", it can maximum reach a speed of 72 km/s. Or if the meteoride orbits the Earth, it could instead spiral down at speeds as low as 7.9 km/s. At 100 km above the surface of the Earth, such burning fragments will travel over a large fraction of the sky and extend its trajectory far beyond the observed 10 * 10 arcmin 2 , while reaching elongated shapes from the burning. Assuming a typical velocity of 12 km/s, an object travels 36,000 km during the exposure time, i.e. nearly the circumference of the Earth at ∼ 100 km. Such a burning process needs to last well below a fraction of a millisecond to be seen as a point source and otherwise shows up as a elongated streak accross the image. Given the extremely short time scale required to see a glint on the POSS-I images under these circumstances, the expected flux dilution is on the orders of ∼ 10 7 , which corresponds to a magnitude dilution of at least 19 magnitudes. Given that the POSS-I plates limiting magnitude is around ∼ 20 mag, such an event will not be observed. The time scale is also inconsistent with the observed, longer ablation times for meteorites. We therefore conclude that asteroids, meteorites and similar objects cannot explain the simultaneous transients.
Airplane strobes can flash fast, but they are also accompanied by continuous lights that leave streaks. Nine high-altitude fireflies will doubtfully all flash synchronized in time and motion. However, rotating geosynchronous satellites (or debris) may produce short sub-second glints with the observed amplitudes. In fact, most of the very short-lived transients are satellite glints 22 . One of the key signatures of small metallic objects with flat, reflective surfaces orbiting our Earth in geosynchronous orbits, is the presence of multiple, very fast and bright glints within the same field-of-view within a few minutes from each other 22 . Many objects only glint once or twice, other can leave several glints along a line 23 . The only problem with this scenario is that no satellites are known to have existed prior to the Soviet-made Sputnik in 1957, seven years after the appearance of the transients in the 1950 POSS-I image.
Finally, we know that the fact that we have 9 simultaneous transients is inconsistent with any well-known natural explanation. This points towards two possible, but not unique, explanations. The first is high-energy contamination of the photographic plates, most likely by radioactive particles. In an epoch of time when many atomic bomb tests were carried out in the United States and in the Soviet Union, this is not an unreasonable explanation, as the Palomar Observatory in California is not located far from the testing sites in Nevada. However, since no official atomic bomb tests were carried between 1949 and 1951 (which means during the exposure of the plate and during its storage time, since plates were normally developed the same night as the observation), we may also consider an alternative, astronomical explanation.
We first consider the possibility of transients from a high-redshift cluster of galaxies. Recently, the discovery of a possible, but extremely unlikely, GRB in a z ∼ 11 galaxy 24 with an apparent mag ∼ 26 (as faint as our possible counterparts), suggested a high rate of GRBs in the early universe 25 . A cluster of high redshift galaxies cannot explain the 9 simultaneous transients because they should continue to flash after the 1950s, which they do not. A follow-up critique showed that the z ∼ 11 transient is more likely explained by a Solar system satellite of artificial or natural origin 26 . If the transient point sources represent real phenomena in the sky, one must also consider how the simultaneous transients can be so well timed when extended over such a large region of the sky, at least covering 10 × 10 arcmin 2 . In order for the transients to have a some synchronization in turning on or off, the furthest transients must be physically separated by less than 30 light minutes, so that their maximum separation is d s < 5.4 * 10 11 m. The distance r to this phenomenon, can then be calculated as r < d s /ω , where ω is expressed in radians. For ω = 10 ′ , this corresponds to r = 1.9 × 10 14 meters or 0.02 light years which is inside our Solar System. Small, flat and highly reflective objects at, or near, geosynchronous orbits around the Earth could produce multiple, fast glints. For an object at a geosynchronous orbit moving with a typical speed of 15"/s, the produced glint must be faster than 0.5 s to look like a transient in POSS-I given the average FWHM ∼ 7", meaning the flux dilution of the glint is ∼ 9.4 magnitudes during the 50 minutes exposure time. Correcting the magnitudes in Table 1 for flux dilution, they reach peak magnitudes ∼ 9 − 11 mags, which is consistent of what is observed for many reflective pieces that are not larger than a few centimeters in size 22 . If these are truly glints, we estimate the rate of cases of multiple glints to be ∼ 0.07 h −1 sky −1 i.e. negligible in comparison to the current number of satellite glints an astronomer at a ground-based telescope near the equator sees today, ∼ 1800 hour −1 sky −1 (see Supp. Info., Sec. A.8.1.). Such glints are so rare that they are likely to go unnoticed. As a curiosity, we note that six of the nine transients (including the most central ones), can be explained by two straight lines each holding three transients together. If on top of this already extraordinary scenario, we see the GTC counterparts in an "off " state, some of the objects need to follow (artificial) trajectories mimicking stars, in order to compensate for the huge expected proper motion and thus moving at least tens of AU every year. In comparison to this, one may find the first proposed scenario with radioactive bomb particles more down-to-Earth. As this case of simultaneous transients is not unique, contamination cannot be excluded as an explanation.
We acknowledge the interesting surplus of star-like objects in the POSS-I catalogue. It is essential to figure out what is causing them, as several of the general VASCO transients belong to the same category. It is also a showcase of the types of anomalies that may be found in the VASCO citizen science project 27 . We believe the mystery of the simultaneous transients is a detective story worth of the attention of the astronomical community.