1. Introduction
A special place among the studied astronomical phenomena belongs to non-stationary ones. In a sense, all astronomy is the science of change (evolution) of the Universe with time, both as a whole and its parts of various scales, from meteors to clusters of galaxies. It is no coincidence that in recent years a new direction has been formed and is rapidly developing—astronomy in the time domain (“Time domain astronomy”). It combines methods, tools and ideas focused on the study of non-stationary phenomena in the Universe at different time and space scales. The website of the IAU working group (URL:
http://timedomainastronomy.net/resources.html#master (accessed on 1 December 2023), Master list of Time Domain Surveys section) provides information on 62 instruments that investigate (or have investigated) variable objects. As long as the list of telescopes is, so is the range of specific tasks for which these instruments are intended. This includes observations of meteors, comets and asteroids, artificial satellites of the Earth, space debris, searches for the effects of microlensing and transits of exoplanets, the study of variable stars, studies of optical afterglows of gamma-ray bursts and searches for optical flares synchronous with these bursts, searches for supernovae, and the study of brightness variations of galactic nuclei, search and characterization of optical counterparts of fast radio bursts and gravitational wave signals. The latter direction is all the more important and exceptionally interesting in connection with the discovery of the gravitational wave event GW170817 [
1], caused by the merger of two neutron stars and the associated short gamma-ray burst GRB 170817A [
2,
3,
4], the kilonova emission [
5,
6] and subsequently with a long multi-wavelength afterglow [
7,
8,
9,
10,
11,
12]. At the same time, there are considerations that precursors of short gamma-ray bursts can precede gravitational-wave pulses [
13], thereby marking the region of their localization, which gives hope for searching for electromagnetic radiation synchronously with the pulse itself. When searching for variability, as a rule, two modes of observation are implemented (we note a certain conventionality of such a division): survey and monitoring. In the first case, the radiation of relatively large areas of the sky is periodically recorded, while the ratio between the time of their exposure (determining the depth of the survey) and the size of the fraction of the celestial sphere observed during the night is determined by the nature of the required variability. Some of these projects and tools are SDSS [
14]; Sky Mapper [
15]; QUEST-La Silla survey [
16]; ASAS-SN [
17]; Catalina Surveys [
18]; ZTF [
19], LSST under construction [
20]. In the second observation mode, we are talking about the longest possible registration of the radiation of a certain object in anticipation of changes in its characteristics (intensity, spectrum, polarization, etc.). An outstanding example of this kind of TAOS II program is the search for trans-Neptunian planets by their occultations of 10,000 stars, whose long-term simultaneous monitoring of high temporal resolution (exposure—50 ms) in individual subapertures is carried out using three telescopes with a diameter of 1.3 m [
21]. Another classic example is the monitoring of UV Cet-type stars to study their flare activity using the largest telescopes. In particular, subsecond polarized spikes of synchrotron origin were first detected during a powerful UV Cet flare with the 6-meter telescope of SAO RAS [
22]. In essence, both modes consist of monitoring, spatiotemporal in the first case, and temporal in the second.
A special place among the non-stationary phenomena studied within the framework of Time Domain Astronomy is occupied by unexpectedly appearing (and also suddenly disappearing) objects, whose localization in space and/or time is not known in advance, and whose duration is rather short—the so-called transients (transient sources). They manifest themselves as stochastic fluctuations of electromagnetic radiation of different frequencies (from radio to gamma-rays), low (tens of MeV) and high (TeV-PeV) energy neutrino events, cosmic rays, and gravitational waves. It should be emphasized that despite the variety of forms of energy release in these events, the possibilities of understanding their nature and constructing their models are determined to a certain extent by the detection and study of their manifestations in the optical range—optical transients. These include non-stationary phenomena in the near-Earth space and the Earth’s (exoplanet’s?) atmosphere, non-moving phenomena such as auroras, transient luminous events (TLE)—elves, sprites, jets [
23], and moving phenomena—comets, asteroids, meteors, satellites, space debris. Finally, undoubtedly, such transients can be optical flashes (including periodic) of artificial origin—signals (transmissions) of extraterrestrial civilizations [
24,
25]. The characteristic duration of these phenomena ranges from minutes-hours to weeks-months (transits of exoplanets, outbursts of novae, supernovae, variable stars, microlensing effects) and from milliseconds to tens of seconds (gamma and radio bursts, gravitational waves, TLE, flyby meteors and satellites). In essence, we are talking about two different types of transients—long and short—the detection and study of which require the use of different methods and tools. At the same time, within the framework of the prevailing ideas and conceptual apparatus, events of the first type are often referred to as fast transients. For example, one of the most extensive modern programs “Wider Faster Deeper” [
26], which unites about 40 ground and space instruments and is dedicated to the study of all of the fastest transients, uses the follow-up mode in the optical range (repointing of the instruments to the area of the already detected fast transients). With the use of the 4-meter Blanco telescope and DECam, at an exposure of 20 s, information on the subsecond fast transients is completely lost, but their environment and slow transients are successfully studied—gamma-ray burst afterglows, localization regions of fast radio bursts and gravitational wave events, bursts of novae, supernovae, red dwarfs [
27].
There is a real contradiction between the need to use instruments with an extremely high temporal resolution in the search and study of fast optical transients and the real characteristics of telescopes, both existing and under construction.
Figure 1 and
Table A1 from
Appendix A demonstrate this perfectly. We present the characteristics of 30 survey instruments focused on the search and study of non-stationary objects of various types, whose fields of view exceed 4 square degrees.
The minimum exposure duration exceeds 10 s for the vast majority of these telescopes, which, given the possibility of summing up a sequence of individual frames, fully provides a depth sufficient for detecting and studying variability up to magnitudes 19–25. At the same time, these limits for systems with close fields of view may differ significantly due to the difference in the diameters of the mirrors. Note, that most of the telescopes in our list are multi-element systems. Indeed, to implement wide-angle spatial-temporal monitoring with a single aperture, a combination of initially mutually contradictory conditions is necessary—a sufficiently high detection limit (large diameter of the supply optics) and a wide field of view (short focus). This circumstance significantly limits the set of possible optical schemes of instruments that determine the optimum of such a combination. So, with an aperture diameter of D > 1–2 m, the size of the field of view is 2–3 degrees, and with D < 1 it can reach 5–10 degrees [
28,
29]. This implies the need to move to multi-telescope systems with fields of view of hundreds and thousands of square degrees, which also have a number of other advantages—relative cheapness, the ability to change the configuration, and (most importantly), since the dimensions of the detectors can be quite small, the use of high-temporal resolution instruments. [
29,
30]. However, as can be seen from
Figure 1 and
Table A1 of
Appendix A, only one multi-aperture system is currently being developed—Argus [
31], which has a time resolution of 0.05 s. It consists of 38 (as a starting point, with at least 900 planned in the future) 20 cm telescopes with a total field of view of 344 (8000 in the final version) sq. degrees. Unfortunately, Argus is not able to change the configuration of the field of view, orienting all telescopes to one region of the sky, which will not allow it to be used as a single telescope with an effective diameter of about 5 m. The remaining three instruments have the limiting dimensions of the fields of view (3–5 degrees) for a monolithic configuration in the Schmidt scheme and are equipped with sufficiently fast receivers. Nevertheless, it is clear that due to the small size of the fields of view when searching for fast transients, they can function only in the alert (follow-up) mode, like many large telescopes with standard fields of view of <1 degree (see, for example, [
32] and references therein). At the same time, this mode does not allow one to detect optical radiation synchronously with the fastest transients—gamma-ray bursts, fast radio bursts, and gravitational waves—for this, it is necessary to continuously monitor the celestial sphere using instruments that have the widest possible fields of view of hundreds and thousands of square degrees. In other words, it is necessary to independently detect and study fast optical transients and only afterward look for their connection with events in other spectral ranges and of a different nature, comparing the time of occurrence and the position of both. We emphasize that the understanding of their origin, and the choice of models that describe them, will largely be determined by the solution of this most complex problem of modern practical astrophysics. The areas of monitoring of gravitational waves, gamma and radio telescopes when searching for transients cover almost half of the area of the celestial sphere (the exception is the BAT detector of the Swift instrument with a field of view of about 5000 square degrees). This means that even with a relatively high accuracy (arc minutes–degrees) of determining the coordinates of an event for gamma detectors [
32], positioning optical telescopes according to an alert will take tens of seconds–minutes. At the same time, the duration of 90% of gamma-ray bursts lies in the 20 ms–100 s interval, and about 30% of the events last less than 2 s [
33]. The latter is the result of the merger of two neutron stars or a black hole and a neutron star and form a class of short gamma-ray bursts (SGRB), in contrast to long-duration GRBs (LGRB), caused by the collapse of some massive stars to form black holes [
34,
35]. Thus, in the alert mode, it is impossible to detect the optical companion of the short gamma-ray burst itself—gamma radiation is no longer detected by the time the telescope is pointed at the event localization region. When a long burst is detected, it is sometimes possible to begin optical observations of its location zone until the gamma-ray emission disappears. As a rule, in both cases, in the alert mode, it is possible to detect and study only the so-called afterglow—optical radiation resulting from the interaction of the plasma ejected during the generation of a gamma-ray burst with interstellar gas [
36].
It should be emphasized that so far, despite numerous works on observations and modeling of the phenomenon of gamma-ray bursts, no self-consistent theory is yet available. It is generally accepted that the study of joint optical radiation with gamma emission (prompt emission), the comparison of their characteristics and temporal structure can provide the key to solving this problem [
34,
37,
38,
39,
40]. At the same time, out of 900 optically identified gamma-ray bursts, only 21 have scant optical information (only 1–3 brightness measurements) simultaneous with gamma-ray emission, and only six events have truly informative optical light curves obtained during the same period, usually in its final phase [
38]. Finally, for the single GRB 080319B (Naked Eye Burst) event, we detected a joint optical burst with a temporal resolution of 0.13 s during the entire GRB by the wide-field monitoring independent of gamma-ray telescope data [
41,
42]. These unique observations with the wide-angle (600 sq. degrees) high-temporal resolution camera TORTORA [
43] made it possible to establish the essential features of the mechanisms of gamma-ray burst generation [
44]. The prototype of this camera is our similar instrument FAVOR (FAst Variability Optical Registrator), which was used for observations in 2003–2009 [
45,
46].
Further, in the course of developing the program for wide-angle high temporal resolution monitoring, SAO RAS (together with KFU and OOO Parallax) created a 9-channel system with a field of view of 900 square degrees and a temporal resolution of 0.1 s, Mini-MegaTORTORA (MMT) [
47]. It was used to detect optical flares synchronous with the bursts GRB 160625B [
48] (however, with a time resolution of 30 s) and GRB 210619B. In the latter case, radiation was recorded in 5 channels simultaneously with exposures of 1, 5, 10, and 30 s in white light (3 channels) and in filters B and V. The implementation of such a research mode is a fundamental feature of our approach, it allowed us to compare the structures of light curves in the gamma and optical ranges, build spectra in these ranges, and finally show that in this case, the burst emission is due to a backward shock wave propagating in a relativistic jet [
49].
Note, that the MMT system is also equipped with polaroids with different orientations for measuring the linear polarization of optical flares, which should provide unique information about the physical properties of the bursts, namely, the structure of the emitting regions, the characteristics of the magnetic field, and the details of the mechanisms for generating radiation of different energies [
34,
50,
51]. Space-telescope measurements of the polarization of hard burst radiation have low accuracy (30–50%) and temporal resolution [
52], which does not allow one to choose between models, providing information only about the most general features of the phenomenon [
37]. Because of this circumstance, polarization studies of optical flashes accompanying gamma-ray emission are of particular importance. Such observations are practically absent—a small number of works on this topic contain the results of measurements of the polarization of optical afterglows, at the level of 10–20%, but minutes–hours after the burst itself ([
53] and references therein). In [
54] an upper limit for linear polarization at a level of 12% has been obtained for the second episode of the GRB 140439A burst 2.5 min after the alert. Thus, the polarization of optical radiation simultaneous with the bursts, as well as its spectral characteristics, have yet to be investigated in high-resolution wide-angle surveys.
At present, one of the unsolved (and most important) problems of observational astrophysics is the detection of optical companions of fast radio bursts—the choice of adequate models of these phenomena from a large number of the ones considered requires the study of their possible optical counterparts, which, unlike gamma-ray bursts, are the subject of only theoretical constructions. Some of the several dozen models predict the generation of sufficiently bright optical flares [
55,
56,
57]. On the other hand, the frequency of occurrence of fast radio bursts throughout the sky can reach several thousand per day [
58], and, consequently, the optical flares accompanying them can be quite frequent. At the same time, as we know, unlike radio observations, optical surveys to search for these events with a high (1 ms) time resolution are not carried out at the moment. As a rule, studies are focused on searching for stationary optical burst hosts [
59], or searching for millisecond optical flashes in repetitive bursts. About 20 such objects, which erupted about 200 times, were found among the almost 800 bursts proper. Their coordinates were determined, which made it possible to observe the regions of localization of these sources using large telescopes (see, for example, [
60,
61]), in which optical flares were not detected—the same result was obtained by us with a 6-meter telescope [
62]. Gravitational wave events can also be accompanied by optical flares ranging from a few milliseconds to a few seconds, especially when they come from merging neutron stars, like the recently discovered event GW 170817 [
63,
64]. At the same time, since the accuracy of localization of gravitational waves is hundreds of degrees, such flares can only be detected by monitoring the celestial sphere independent of the data of gravitational detectors using instruments with comparable fields of view and millisecond time resolution. Moreover, in the constantly updated catalog of fast radio burst models [
65], probable deep connections between them, the gravitational waves and gamma-ray bursts are clearly visible, based on the universal cause of these phenomena—the interaction and evolution of relativistic objects. It is all the more important to look for signatures of similar phenomena in the optical range using wide-angle monitoring of high temporal resolution. Based on the above considerations, and on the experience of developing and using the FAVOR, TORTORA, Mini- MegaTORTORA tools, we proposed the concept of a multi-telescope complex with high temporal resolution for wide-angle monitoring of the celestial sphere—SAINT (Small Aperture Imaging Network Telescope) [
62,
66,
67,
68] using telescopes with a diameter of 30–50 cm, with the limiting magnitude at the level of single- component instruments (see
Figure 1 and
Table A1 in the
Appendix A) and capable of detecting and investigating non-stationary objects at times from milliseconds to months-years. This paper is devoted to the development of this project using modern devices, the choice of the design of the complex, its components, mode of operation, and the assessment of its parameters.
2. Small Aperture Imaging Network Telescope (Saint) Project
The performed analysis demonstrates the obvious need to create a multi-channel optical monitoring system that would be capable of conducting complex searches and studies of non-stationary objects and phenomena within the framework of universal instrumental and methodological approaches and would not have the disadvantages of existing systems.
The features of this kind of instrument should be as follows:
a large field of view of several hundred degrees, which requires a multi-channel design;
high temporal resolution in the range of tenths to hundredths of a second;
detection limit of 18–20 mag at time scales of 20–30 s to be at the level of the limits of existing survey telescopes (
Figure 1);
a combination of monitoring (wide-field) and follow-up (narrow-field) modes with their rapid change;
the possibility of obtaining maximum information (spectral, polarimetric, photometric) about non-stationary objects in follow-up mode;
processing of accumulated information in real time to detect, characterize and classify transient phenomena, and make decisions on the transition to the follow-up mode;
preservation of all raw and reduced data, as well as maintenance of databases obtained as a result of its a posteriori analysis;
complete robotization of the complex operation using information from external sources (meteo station, network, other instruments) and a system to control its condition.
As a practical implementation of these principles, we propose the project of SAINT (Small Aperture Imaging Network Telescope)—12-telescope complex with high temporal resolution, built mostly from off-the-shelf commercially available components that may be implemented in about
$2 M cost (plus additionally about
$500 K for labor costs). Its overall parameters are listed in
Table 1, schematic view of a single channel is shown in
Figure 2, and details on individual components and details of its operation modes are given below.
2.1. Telescope
The basic component of the complex is GENON Max—a catadioptric two-mirror Schenker–Terebizh telescope with five corrective lenses [
29] with a diameter of 30 cm (F/1.5) and a field of view of about 40 square degrees (
Table 2). In its Cassegrain focus a multimode photopolarimeter of a rather complex design is installed, capable of operating both in wide-field monitoring and in the study of a single object (
Section 2.4).
2.2. Mount
The telescopes are mounted on ASA DDM100 mounts based on direct drive technology (see
Table 3). This design has a number of important advantages listed below, which are especially significant for the proposed complex.
High acceleration and, as a result, high pointing speed up to 50/s, which will allow changing the field of view of the complex and the observation mode in several seconds.
Due to the absence of drive belts and gears, the rotation of the axes is very uniform and silent, which ensures the accuracy of the guidance and the absence of vibrations, which determine the high stability of the position of sources in the focal plane of telescopes.
Due to the absence of gearboxes, the absence of mechanical backlash and hysteresis is ensured. For the same reason, moving parts wear out little, which increases the period of maintaining the stability of the mount’s characteristics.
All movable structural elements are equipped with 28-bit encoders and feedback systems, which ensures the accuracy of rotation angles and their determination up to 0″.004.
Maintaining the stability of the functioning of all elements of the mount and its control is provided along with feedback systems by a high-speed controller online with a computer.
2.3. Detector
Every channel of SAINT is equipped with a photopolarimeter (
Section 2.4), which is the primary instrument in the SAINT complex. The detector used in it is Andor sCMOS Balor 17F-12 camera, whose characteristics are listed in
Table 4.
This detector currently has the best combination of a large size of 16.9 Mpix (70 mm diagonal), high temporal resolution when reading a full-frame (about 20 ms) and a fairly low readout noise (2.9 e
−), which is unattainable in CCD matrices for several data transmission channels [
69]. A sufficiently large pixel size (12 µm) allows the use of the maximum temporal resolution in observations even with relatively low image quality (>1
), nevertheless concentrated in one pixel, which keeps the readout noise to a minimum. At the same time, a high degree of signal response linearity (>99.7%) is combined with a large electron well depth (80,000) and a special way of amplifying and digitizing the signal, which ensures its continuous recording in the maximum dynamic range. These features make it possible to detect sources with an intensity from the read-out noise level to the saturation limit without distortion in a single camera frame, which is fundamental in wide-field observations. Moreover, reducing the read-out area (Region of Interest, ROI) (
Table 5) allows for improving the temporal resolution down to 2.5 or 0.6 ms which may be used in the regime of targeted observations of individual objects.
The microlens array on top of the sensor in the Balor camera allows lossless operation for the relative apertures of up to F/0.3 at cone angles up to 110, which makes it possible to avoid light losses even at the outer margins of the telescope field of view at aperture ratio of F/1.5. An important point is what the camera allows to timestamp and synchronize the acquired frames using the signal from the GPS receiver which is crucial for the combination of data from different channels pointed toward a single target. On the other hand, in contrast to traditional CCD sensors, the question of long-term spatio-temporal stability of individual pixels of CMOSes poses a potentially significant problem for the co-addition of images from such detectors, as it introduces additional per-pixel noise that cannot be mitigated by classical calibration methods (bias and dark frame subtraction, non-linearity correction and flat fielding). The limits imposed by these low-frequency effects on frame co-addition are still poorly studied and have to be thoroughly investigated prior to the decision on an exact strategy of frame co-addition in the operation of the complex.
2.4. Photopolarimeter and Different Modes of Operation
The primary instrument of every SAINT channel, and the only custom-built part of it, is a photopolarimeter which provides observations in both monitoring and research modes. It is attached to the Cassegain focus of the telescope and consists of two parts—an optical-mechanical unit and a camera unit, its general view with installation on the telescope is shown in
Figure 2.
The camera unit houses the detector, which is mounted on sliding rails that allow it to move along the optical axis, thereby focusing images on the sensor. The hoses of the water cooling system of the chamber, not shown in the figure, are inserted into the rear end of the casing of the unit.
The structure of the photopolarimeter is shown in the diagrams of
Figure 3, its internal view is in
Figure 4, its longitudinal view is in
Figure 5. Control electronics units, power supply, motors that ensure the movement of moving elements of the device, and connecting cables are not shown here.
The front part of the photopolarimeter, located in the gap between the back focus of the telescope and the entrance of the camera, contains the following components: a turret with filters, an empty window, and a narrow-field block containing mirrors 1 and 4, as well as a field-of-view diaphragm (the optical axis of the telescope passes through the centers of the turret holes); input (COMPUTAR M2514-MP2 with a focal length of 25 mm and aperture ratio of 1:1.4) and output (COMPUTAR V5014-MP with the focal length of 50 mm) lenses of the collimator; mirror 3.
The side part of the photopolarimeter is located outside the camera casing. It contains a turret with narrow-field filters and an empty window, as well as a turret coaxial with it with three polaroids of different orientations (with polarization plane rotated by 120 in respect to each other) and an empty window; mirror 2.
During observations in the wide-field (monitoring, or primary survey) mode, an empty turret window 1 or one of its filters is installed on the optical axis of all 12 telescopes (depending on the specific tasks of the survey, observations of different areas with different filters, etc. are possible). The cameras acquire a sequence of full-frame images (4128 × 4104 pixels, or 6.3× 6.3 each) with a rate of 54 fps, and the entire complex monitors the overall field of view at 470 sq. degrees.
In the narrow-field (research) mode, a block with mirrors 1, 4 and a field diaphragm measuring 11
× 11
(1.45 × 1.45 mm) of turret 1 (see
Figure 3) is installed on the optical axis, while mirror 1 focuses the redirected beam on the diaphragm placed at the focus of the input lens of the collimator. The parallel beam constructed by it passes through one of the
filters (or empty window) of turret 2 and one of the polaroids (empty window) of turret 3. After being reflected by mirrors 2 and 3, the parallel beam is focused by the output lens of the collimator and redirected by mirror 4 to the sensor of the camera. As a result, the light from the 11
× 11
sky region with a linear size of 3 × 3 mm (240 × 240 pixels) is imaged with a scale of 2.75
/pix, twice the original (5.5
/pix)— the scale is determined by the focal length of the collimator output lens.
This image can be recorded with a time resolution of 20 to 2.5 ms, and its central 6
× 6
part (128 × 128 pixels)—with a resolution of 0.6 ms (see
Table 5). The transition from monitoring to research mode takes a few seconds and is determined by both the rotation speed of the photopolarimeter turrets and the repointing speed of the mount. In particular, when an optical transient is detected during real-time monitoring, the following operations are performed:
its characteristics are determined (initial brightness, characteristic duration, structure of the light curve);
rough initial classification is performed (as one of noise, meteor, satellite, new or already known astrophysical object classes);
the follow-up mode for every channel is selected (photometry and/or polarimetry, filters to use, exposure time);
in parallel, all telescopes are repointed to the source, placing it in a small central field (with the exception of the instrument that detected the transient, it retains the original mode used for detection, in order to keep an uninterrupted sequence of data);
observations start.
Table 6 shows the estimated detection limits of the complex in different modes and with different time resolutions. Estimates were made for a single exposure in white light at a characteristic wavelength of 5500 Å, a sky background of 21.5 mag/sq.arcsec, an optical and atmospheric transmission of 0.25, the characteristics of the telescope and detector from
Table 2 and
Table 4, and assuming that the object flux is fully contained inside a single pixel (i.e., significantly undersampled PSF). It should be noted that the real limits when using filters and polaroids will be 0.5–1.0 magnitude worse.
2.5. Shelter and Infrastructure of the Complex
Twelve telescopes of the SAINT complex are installed in two cylindrical shelters housing six channels each (see
Figure 7). They consist of a reinforced concrete supporting slab mount and a steel truss frame sheathed with sandwich panels on which it is mounted. Its free part serves as a flyover to accommodate the cylindrical roof of the shelter in the working state of the complex, which, moving along rail guides, closes the plate with the telescopes in the non-working state. The frame is fixed on a reinforced concrete base covered with metal-plastic cladding, and installed on a concrete pad dug into the ground.
The internal volumes of the frame are the technical compartments of the SAINT complex. They house the following equipment and infrastructure components:
power supply input unit of the complex network, including power drives of all its mechanical elements, electronic components of the telescopes and instruments, computer cluster, and uninterruptible power supplies;
system for providing water cooling of detectors;
climate control system based on dehumidifiers COTES CR300.
Both parts of the shelter are equipped with video cameras for monitoring the internal and external space. The complex is also equipped with a meteorological system, a fish-eye all-sky camera and an IR cloud sensor, which allows to automatically monitor weather conditions and make a decision on the start and end of observations.
The movable roofs of the shelters are covered with a sun-protection coating, and flat screens with uniform illumination are placed on their inner surface to calibrate the sensitivity of the detectors.
3. Data Acquisition, Processing and Storage
During the primary monitoring mode of operation, the complex will acquire the data for every sky field (470 square degrees) for 20 min in a row, thus covering up to 12,700 square degrees in a typical 9-hour observational night. The data will be processed in real time, and a number of data products (transients on various time scales, co-added images with different temporal resolutions, etc.) will be derived and stored.
As both the computational hardware and data analysis algorithms evolve very quickly, it is not possible (and not necessary) to specify the IT infrastructure of the SAINT complex here in full detail as we did for the telescope itself. Therefore, below we will outline its generic structure, requirements and possible approaches for handling an enormous amount of data the complex will produce.
In general, the operation of the complex will be carried out using the software that is based on that we used in observations with the FAVOR, TORTORA and Mini-MegaTORTORA instruments (see Beskin et al. [
47] and references therein). On the shortest time scales, where timing constraints are the most demanding, transient detection will be performed in real time using fast image subtraction methods that are based on sufficient local stability of image PSF, as described e.g., in Karpov et al. [
70]. On longer time scales, where PSF stability cannot be ensured, the methods developed for image subtraction with PSF matching developed over the last 20 years will be used, such as Alard and Lupton [
71] family of methods, ZOGY [
72] algorithm, or SFFT [
73]. These methods generally allow parallelization using modern GPU hardware and thus may be expected to run fast enough to be used for analyzing the data in real time on minutes to hours time scale. An alternative to image subtraction methods may also be implemented using convolutional neural networks (CNN, see e.g., Sedaghat and Mahabal [
74] or Acero-Cuellar et al. [
75]). CNNs will also be employed in order to reliably separate the image subtraction artifacts, cosmic rays, as well as, e.g., the effects of stellar scintillations on short time scales from bona fide transients using the methods developed in e.g., Cabrera-Vives et al. [
76], Gieseke et al. [
77], Carrasco-Davis et al. [
78], Makhlouf et al. [
79].
Image co-addition will also allow simultaneously achieving the highest possible temporal resolution (on individual frames) and going for deeper objects (in running coadds of various lengths, corresponding to, e.g., every 100 or 1000 original images). Unfortunately, without a thorough laboratory study of the stability of the detector, it is impossible to define the limits for such co-addition, as at some point the effects of pixel parameter drift specific for CMOSes will be more significant than the further sensitivity gain. However, we may foresee that implementing some tailored observational strategy in the data sequence—e.g., dithering using some pre-defined pattern—will allow us to mitigate such effects to some degree.
As a minimal architecture for the data processing cluster for SAINT, we propose a hierarchical scheme with two identical computers per telescope channel—one for data acquisition from the CMOS and real-time image processing, second for the channel hardware and high-level controls, intermediate data storage and image processing—plus a central computer that handles the overall operation of the complex, including survey scheduling and reaction to the detected transients, communication with external networks, as well as various databases. Such architecture also allows for the use of all cluster machines for more computationally-intensive data processing tasks during the day time, when no observations are performed and no real-time data processing is necessary.
Let’s outline the basic requirements for data processing and storage infrastructure. They are defined by the extremely large amount of data being acquired from several large format detectors (16.9 million pixels each for the Balor camera) routinely operating at frame rates of up to 54 frames per second (or even up to 1684 fps in follow-up regime), thus amounting to 1.8 gigabytes of data per second for every channel, or up to 800 Terabytes per night for the whole complex. As a minimal realistic configuration using present-day commercially available hardware, we may use a camera connected to the PC through a four-channel CoaXPress interface, and a RAID-0 array consisting of five SATA-III hard drives 20 terabytes each. Due to parallel writing, a speed of up to 4 gigabytes per second may be achieved for such RAID configuration, which is sufficient for handling the real-time data flow. Longer-term data storage will require datacenter-class hardware which is also commercially available.
4. Expected Results
The proposed complex is able to operate in still poorly studied region of the parameter space shown in
Figure 1, and thus may provide a number of results related to previously unknown classes of rapid optical transients, as well as a vast amount of data on known ones. Below we will briefly outline some of the fields where such instruments may provide important developments.
Flaring stars are now being routinely discovered in modern time-domain sky surveys like Kepler [
80], TESS [
81,
82], EvryScope [
83], NGTS [
84], ZTF [
85] or Tomo-e Gozen [
86] what perform high cadence continuous observations of the fixed sky regions. However, most of these experiments do not provide temporal resolution better than half a minute, or in rare cases—better than a second, and are being performed in a single photometric filter that somehow reduces the scientific content of collected statistical results. Thus, using a specifically tailored monitoring mode with SAINT which is possible due to its multi-channel architecture—i.e., simultaneous multicolor observations—would allow us to augment the statistics of white-light flares with at least some estimations of their temperature. Moreover, employing the polarimetric mode for targeted (or follow-up, for the flares detected by real-time processing pipeline) observations would also allow placing statistically reliable upper limits on (or even obtain a detection of) the polarized components of stellar flares of various durations, including the shorter ones not typically resolved by existing surveys (see also
Figure 8 for details).
Aizawa et al. [
86] reported 22 flares detected in 40 h of 1 s cadence observations with Tomo-e Gozen camera (21 square degrees simultaneous field of view). While the depth of the SAINT survey is not as good, due to the superior field of view we may expect up to 12 flares per hour in a wide-field mode, thus obtaining a significantly larger amount of flashes, covered with better temporal resolution.
The observations of faint meteors with the Mini-MegaTORTORA system [
47] typically resulted in the detection of several hundreds of events per observational night, with maximum brightness down to 8–10 magnitudes, which is 2–4 magnitudes deeper than the limits of most modern meteor detection experiments (see, for example, Weryk et al. [
87], Jenniskens [
88]) and comparable to the capabilities of the Tomo-e Gozen camera [
89]. Such faint, and even fainter meteor events are caused by micrometeoroids—particles of interplanetary (interstellar?) dust with masses of 0.1
g–0.1 g [
90]—that, according to direct transatmospheric studies, make up the majority of the cosmic matter falling to the Earth, 100–200 tons/day [
91]. On the other hand, observations of meteors give an order of magnitude smaller volume [
92]—this discrepancy is most probably due to an underestimation of the role of ablation and fragmentation processes of micrometeoroids, as well as inaccuracy in determining their velocities [
92]. These problems are largely associated with the difficulty of interpreting radar observation data, which is the main means of studying such faint meteors with optical magnitudes of 10–14 [
93]. Thus, carrying out a long-term survey of a population of faint meteors, covering the entire range of velocities from 12.3 to 71.9 km/s, will allow us to study the features of their optical emission, taking into account possible fragmentation, and compare it with dynamic models of micrometeoroids. Another direction of meteor study might be the analysis of their mass distribution both in the showers and in the sporadic component, as well as massive colorimetry of events during the peaks of different meteor showers in order to compare their spectral properties.
UV Ceti type flare stars belong to the main sequence red dwarfs, making up 70–75% of the population of the Galaxy [
94]. The duration of their stay on the main sequence exceeds the age of the Universe, and during this time their main physical characteristics remain unchanged [
95]. These features determine great interest in red dwarfs as possible regions where life arises [
96]. The probabilities of the occurrence of Earth-like planets orbiting these stars in habitable zones with a characteristic size of about 0.2 AU [
97], estimated from various observations, are in the range of 30–50% [
98,
99]. This gives about 100 habitable Earth-like planets in the vicinity of 10 pc, and several billion in the entire Galaxy [
96]. The small masses and sizes, as well as the low luminosity of red dwarfs, as a result of which the habitable zones are located close to the star, determine the relatively high probability of detecting Earth-like planets during transits. Their depth for earth-like planets ranges from 0.002 to 0.0084 magnitudes with the radius of the host star from 0.6 to 0.1 solar, and the probability of a successful “edge-on” orientation of the orbit of the star-planet system for an observer is 0.5–1% [
100]. The photometric accuracy required for transit detection is close to the limit for ground-based instruments. Nevertheless, it is achieved in various programs, thanks to the optimal strategy of observations and selection of objects, the use of effective methods of data analysis and methods of modeling the processes of obtaining them. Thus, when studying several dozen red dwarfs of 12–15 magnitudes with systems of distributed telescopes with diameters of 0.2 to 2 m at exposures of 30–120 s, measurement accuracies in the range of 0.001–0.004 magnitudes were achieved, leading to the detection of effects with a depth of 0.005–0.007 magnitudes for planets of 1–3 Earth masses located in habitable zones, with a probability of 2.5–8% [
101,
102,
103]. These results were obtained from large data sets using BLS [
104,
105] and TLS [
106] folding algorithms in the space of periods and filling factors.
In this context, it seems that such studies using the SAINT complex should be very effective. Instead of simultaneously monitoring up to 10 red dwarfs (according to the number of instruments used) within the framework of the mentioned programs, SAINT will be able to observe more than a thousand stars. With a minimum brightness of objects with potentially habitable Earth-like planets at a level of 16 magnitude, they will be located at a distance of up to 100 pc, their number in the northern sky will be about 100,000 [
96], and transits of these planets may be detected at several hundred red dwarfs. Note, that when observing each star about 50 times over a visibility period of 100 days and a transit duty cycle of 0.001–0.005 [
100] the accuracy of intensity determination in phased light curves increases by tens of times relative to its initial estimates in individual exposures. Thus, as a result of monitoring the sky with the SAINT complex, it is possible to increase the sample of Earth-like planets in the habitable zones by 4–6 times relative to its current volume (See Habitable Exoplanet Catalogue at
http://phl.upr.edu/projects/habitable-exoplanets-catalog (accessed on 1 December 2023)).
Similar results within the framework of the strategy for detecting periodic signals in wide-angle surveys can be obtained when searching for transits of Earth-like (and not only) planets around white dwarfs. These objects may also host planetary systems and Earth-like planets located in habitable zones. Indeed, with a luminosity of
and a mass of 0.5–0.6 solar, their radii are close to those of the Earth, and the temperature lies in the range of 4000–9000 K, which leads to the existence of a habitable zone in the 0.005–0.02 AU range of distances from stars over several billion years [
97,
107,
108]. On the other hand, there is evidence of the existence of various proto-planetary and post-planetary structures around white dwarfs—disks, planetesimals, asteroids [
109,
110]. Finally, four giant planets were discovered in systems with white dwarfs [
111], one of them by the transit method [
112]. Thus, one can hope also for the presence of Earth-like (rocky) planets with a mass of 1–2 of Earth one in the habitable zones of white dwarfs, as established for 63 G-M dwarfs (according to Habitable Exoplanet Catalogue). The depth of their transits, repeated with periods in the interval of 4–30 h (corresponding to the size of the habitable zone), is 10–100%, the duration is 1–2 min, and the probability of the best “edge-on” orientation of the orbital plane relative to the observer is close to 1–2% for planet radii 1–2 Earth’s [
107,
108,
113]. In ground-based observations, even with small-diameter telescopes and standard photometric accuracy, it is possible to search for such effects. The primary difficulty in such studies is the needed relatively high temporal resolution at the level of several seconds. In particular, Faedi et al. [
113] searched for transits in long-term monitoring of 194 white dwarfs of 9–15 mag with a time resolution of 30 s (an 8-channel SuperWASP telescope was used—see
Figure 1 and
Table A1) and obtained only upper limits at 10% for the frequency of occurrence of planets and brown dwarfs. The authors come to the conclusion that it is necessary to increase the studied sample, and the duration of observations, and to improve the temporal resolution. Similar programs have been planned, and are being partially implemented using space-borne instruments [
114,
115], ground-based telescopes [
116,
117] and their combinations [
118].
In the northern sky, according to various estimates, there are 10–15 thousand white dwarfs brighter than 18 magnitude (see, for example, Kleinman et al. [
119] or Gentile Fusillo et al. [
120]), and for approximately 100–150 ones the transits with periods between 3 h and 50 days, with duty cycles of 0.03–0.0006 and depths 0.1–1, may be detected by SAINT. Note, however, that this estimate is solely for the objects where it is possible to detect planets if they are there—and the probability of it is unknown to us. On the other hand, in the absence of an effect, it will be possible to obtain an estimate of the upper limit for this probability.
Finally, there are several other (high risk) directions placed in the area of “Terra incognita” in
Figure 8 where SAINT observations may be potentially crucial:
detection of optical flashes accompanying gamma-ray bursts (about five events per year?),
registration of optical companions of fast radio bursts, in particular, the ones with repeating activity, where the same approach as for exoplanets may be used to improve the performance,
registration of optical companions of gravitational wave events or setting upper limits on their electromagnetic counterparts’ energy,
search for signals of terrestrial civilizations.