A fraction of Galactic stars have compact companions that could be white dwarfs WDs neutron stars NSs or stellarmass black holes SBHs. In a detached and edgeon binary system including a mainsequence star and a compact object denoted by WD mainsequence WDMS NS mainsequence NSMS and SBH mainsequence BHMS systems the stellar brightness can change periodically due to selflensing or eclipsing features. The shape of a selflensing signal is a degenerate function of stellar radius and the compact objects mass because the selflensing peak strongly depends on the projected source radius normalized to the Einstein radius. Increasing the inclination angle i changes the selflensing shape from a strict tophat model to one with slowly increasing edges. We simulate stellar light curves due to these binary systems that are observed by NASAs Transiting Exoplanet Survey Satellite TESS telescope and evaluate the efficiencies to detect their periodic signatures using two sets of criteria i signaltonoise ratio SNR gt3 and N SUBtranSUB gt 1 low confidence LC and ii SNR gt5 and N SUBtranSUB gt 2 high confidence HC. The HC efficiencies for detecting WDMS NSMS and BHMS systems with the inclination angle i lt 20 during different time spans are 57 4.56 and 45 respectively. Detecting lensinginduced features is possible in only 3 and 33 of detectable WDMS and NSMS events. The detection efficiencies for closer source stars with higher priorities are higher and drop to zero for b  R SUBSUB where inlineformula inlineformula is the impact parameter a is the semimajor axis. We predict the numbers of WDs NSs and SBHs that are discovered from the TESS Candidate Target List stars are 1518 67 and lt1.
['Compact objects', 'Gravitational lensing', 'Compact binary stars', 'Stellar remnants', 'Transient detection', 'Eclipsing binary stars', 'Astronomical simulations', 'Space telescopes', '288', '670', '283', '1627', '1957', '444', '1857', '1547', 'Astrophysics - Solar and Stellar Astrophysics', 'Astrophysics - Earth and Planetary Astrophysics', 'Astrophysics - Astrophysics of Galaxies', 'Astrophysics - Instrumentation and Methods for Astrophysics']
Simulating Selflensing and Eclipsing Signals due to Detached Compact Objects in the TESS Light Curves
{'Simulating Self-Lensing and Eclipsing Signals due to Detached Compact Objects in the TESS Light Curves': 'Sedighe Sajadian 1 and Niayesh Afshordi 2, 3, 4  \n1 Department of Physics, Isfahan University of Technology, Isfahan 84156-83111, Iran, s.sajadian@iut.ac.ir  \n2 Waterloo Centre for Astrophysics, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada  \n- 3 Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada  \n4 Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada', 'ABSTRACT': "A fraction of Galactic stars have compact companions which could be white dwarfs (WDs), neutron stars (NSs) or stellar-mass black holes (SBHs). In a detached and edge-on binary system including a main-sequence star and a compact object (denoted by WDMS, NSMS, and BHMS systems), the stellar brightness can change periodically due to self-lensing or eclipsing features. The shape of a self-lensing signals is a degenerate function of stellar radius and compact object's mass because the self-lensing peak strongly depends on the projected source radius normalized to Einstein radius. Increasing the inclination angle i changes the self-lensing shape from a strict top-hat model to one with slow-increasing edges. We simulate stellar light curves due to these binary systems which are observed by NASA's Transiting Exoplanet Survey Satellite (TESS) telescope and evaluate the efficiencies to detect their periodic signatures using two sets of criteria (i)SNR > 3 and N tran > 1 (Low-Confidence, LC), and (ii) SNR > 5 and N tran > 2 (High-Confidence, HC). The HC efficiencies for detecting WDMS, NSMS, and BHMS systems with the inclination angle i < 20 · during different time spans are 5-7%, 4 . 5-6%, and 4-5%, respectively. Detecting lensing-induced features is possible in only ≲ 3% and ≲ 33% of detectable WDMS and NSMS events. The detection efficiencies for closer source stars with higher priorities are higher and drop to zero for b ≳ R ⋆ , where b ≃ tan( i ) a is the impact parameter( a is the semi-major axis). We predict the numbers of WDs, NSs, and SBHs that are discovered from the TESS Candidate Target List stars are 15-18, 6-7, and < 1.  \nKeywords: Space telescope- Astronomical simulations - Eclipsing binary stars - Transient detection Stellar remnants - Compact binary stars - Gravitational lensing - Compact objects", '1. INTRODUCTION': "The NASA's Transiting Exoplanet Survey Satellite (TESS 1 , Ricker et al. (2014)) telescope observed (and observes) stars in the Candidate Target List (CTL, Stassun et al. (2018, 2019)) with a 2-min cadence and an accuracy better than 60 ppm on hourly timescales. Although its main goal from these observations is to discover Earth-size planets transiting bright stars in the solar neighborhood, its observing strategy is uniquely matched to capture other types of periodic and weak variations in stellar light curves.  \nFor instance, stellar light curves from edge-on and detached binary systems including main-sequence stars and compact objects, i.e., white dwarfs (WDs), neutron  \nstars (NSs), and stellar-mass black holes (SBHs), have potentially self-lensing, eclipsing (blocking the compact objects' brightness by their stellar companions), and occultation (blocking the light of stellar images by their compact companions) signals which all of them are periodic. Self-lensing refers to the lensing effect on the flux of the luminous object by its compact companion in edge-on systems (see, e.g., Gould 1995). In this regard, simulating and numerically studying selflensing/eclipsing/occultation features that are realizable in the TESS data have several advantages as (i) revealing the importance of searching these signals in the TESS data, (ii) learning the machines to extract the real events from a huge ensemble of the TESS data archive which is growing up with time, and (iii) evaluating the known models describing binary systems, NSs,  \nand SBHs through comparing their results with ones from the real observations.  \nIn this work, we therefore simulate these periodic features and study their properties and detectability in the TESS data. However, Masuda & Hotokezaka (2019) roughly evaluated the number of detectable BHs in the mass range [1 , 100] M ⊙ in tight and detached orbits around stars with the orbital period [0 . 3 , 30]days through either self-lensing signals or phase-curve variations of stellar light curves. They injected these photometric signals in spotted stellar light curves detected by the Kepler telescope and concluded that 10-100 BHs can be detected by searching 10 5 -10 6 stellar light curves. Also, Wiktorowicz et al. (2021) studied detection of BHs/NSs through their self-lensing signals in a synthetic ensemble of binary systems and reported a higher number for detectable BHs. Here, we extend their works by (i) considering binary systems including different types of compact objects (i.e., WDs, NSs, and BHs) separately, (ii) modeling eclipsing and occultation in addition to self-lensing signals, (iii) generating synthetic light curves and data points based on the real observations by the TESS telescope, and (iv) taking their companion stars from the TESS CTL targets.  \nHere, we first review the known properties of binary systems including main-sequence stars and compact objects in three following paragraphs.  \nWD main-sequence (WDMS) binaries: It is predicted that the number of WDMS binaries in our galaxy is 10 7 -10 8 which depends on their mass-ratio distribution. The number of detached WDMS binaries is more than that of interacting ones by more than one order of magnitude (e.g., see Willems & Kolb 2004), nevertheless detecting interacting WDMS binaries is easier and plausible via either spectroscopic observations or variability surveys. For instance, the Sloan Digital Sky Survey (SDSS, York et al. (2000)) telescope discovered more than 3200 WDMS binaries (see, e.g., Heller et al. 2009; Rebassa-Mansergas et al. 2016). Additionally, ∼ 320 , 000 (either candidate or confirmed) WDs in binary systems have been reported in the Galaxy Evolution Explorer (GALEX, Martin et al. (2005)) and the Gaia Data Release 3 (GDR3, Gaia Collaboration et al. (2023)) (Bianchi et al. 2011; Gentile Fusillo et al. 2021). For detecting detached WDMS binaries, a common method is spectroscopic observations because they can make a combined spectrum due to one star and one WD. If these systems are edge-on as seen by the observer, additionally eclipsing, self-lensing signals or variation in stellar radial velocities are realizable through precise photometric and spectroscopic observations (Korol et al. 2017).  \nNS main-sequence (NSMS) binaries: Another type of compact object is an NS which could be a companion for a main-sequence star in a binary system. A NS is formed after the collapse of a massive star with an initial mass higher than 8 M ⊙ , whereas more massive stars with initial masses higher than ∼ 20 M ⊙ are converted to SBHs after their gravitational collapse. Isolated NSs can be discovered through radio emissions, the so-called pulsars (e.g., Zhang et al. 2022). NSs and SBHs in interacting binary systems with mainsequence stars, the so-called X -ray binaries, are discernible through X -ray emissions owing to mass transferring from their stellar counterparts. Depending on the dominant mechanism for mass transferring toward the compact objects, these binaries are divided into two subclasses: (i) low-mass X -ray binaries (with RocheLobe overflow), and (ii) high-mass X -ray binaries (with wind-accreting) (see, e.g., Ogelman & Swank 1974; Pfahl et al. 2002; Casares et al. 2017). Around 4% of all discovered NSs are in binary systems (Tauris & van den Heuvel 2006), and our galaxy hosts up to one billion NSs, whereas only ∼ 4 , 000 NSs have been discovered up to now.  \nSBH main-sequence (BHMS) binaries: Black holes with masses ≲ 100 M ⊙ are the so-called stellarmass black holes. The lightest SBH was discovered up to now has the mass ≃ 3-3 . 3 M ⊙ which was located in the mass-gap between NSs and SBHs (Thompson et al. 2019; Jayasinghe et al. 2021; Ozel et al. 2010; Farr et al. 2011). Although the number of SBHs in our galaxy is predicted to be more than 10 million, up to now only 72 SBHs were confirmed mostly through X -ray transients from their accretion disks in interacting BHMS binaries 2 (Corral-Santana et al. 2016). 20 of these discovered SBHs were confirmed dynamically by discerning periodic variations in the radial velocities of their luminous companions.  \nIn addition to spectroscopic and variability surveys, and X -ray observations to capture compact objects in binary systems with main-sequence stars, there are other channels for discovering these objects, e.g., (i) based on their gravitational impacts which is the so-called selflensing effect, (ii) eclipsing signals, (iii) ellipsoidal variation of stellar brightness, (iv) Doppler boosting, etc (e.g., see, Masuda & Hotokezaka 2019; Sorabella et al. 2022). Two last effects happen for massive compact companions and small orbital radii. Unlike other types of lensing, including microlensing, strong lensing and weak lensing which all are not repeatable, a self-lensing  \nsignal is periodic and its period is exactly equal to the orbital period of the binary system.  \nUp to now, five self-lensing/eclipsing binary systems containing white dwarfs and main-sequence stars were discovered through photometric observations by the Kepler telescope (Kruse & Agol 2014; Kawahara et al. 2018; Masuda et al. 2019). All of these systems are wide binaries with orbital periods from 88 to 728 days. The next generation of the Kepler telescope is the TESS telescope which was launched on 18 April 2018. During its two-year primary mission, it covered 85% of the sky by dividing it into 26 sectors. The area of each sector is 24 × 90 square degrees, and sectors have some overlapped parts over the ecliptic poles. Each sector is observed during two 13 . 7-day observing periods with a one-day gap in the middle. The TESS pixel scale is 21 arc-second which leads to a high photometric accuracy better than 60 ppm (parts per million) for brightest stars on hourly time scales. This telescope could also re-cover self-lensing signals in a binary system originally discovered by the Kepler telescope, i.e., KIC 12254688 (Sorabella et al. 2024).  \nIn this work, we simulate possible self-lensing and eclipsing signals due to compact objects in stellar fluxes from detached binary systems as seen by the TESS telescope, to (i) study the characteristics of self-lensing signals and (ii) estimate the TESS efficiency for detecting them. We also investigate how this detection efficiency depends on the physical parameters of compact objects, source stars, and binary orbits. We finally estimate the numbers of WDs, NSs, and SBHs that are realizable through their self-lensing/eclipsing signals in the photometric data of the TESS CTL targets.  \nThe outline of this paper is as follows. In Section 2, we explain the details of simulating self-lensing, occultation, and eclipsing signals and then discuss on their characteristics as a function of models' parameters. In Section 3, we explain the details of Monte Carlo simulations from self-lensing, occultation, and eclipsing light curves due to detached and edge-on binary systems while we assume that these events are observed by the TESS telescope. We extract the detectable events based on two sets of criteria. In Section 4, we explain the results and conclusions.", '2. BINARY SYSTEMS CONTAIN COMPACT OBJECTS AND MAIN-SEQUENCE STARS': 'To simulate a stellar light curve from a binary system including a compact object and a main-sequence star, in the first step, we simulate a binary orbit as explained in Subsection 2.1. By assuming that its orbital plane is edge-on as seen by the observer, in the next step we cal-  \nculate its self-lensing, occultation, and eclipsing signals. All details are explained in Subsection 2.2. We then discuss the characteristics of self-lensing signals as functions of orbital properties and parameters of compact objects in Subsection 2.3.', '2.1. Simulating a Stellar Orbit Around a Compact Object': "Let's consider a detached binary system containing two companions: a star and a compact object with masses M ⋆ , and M c ( c can be either WD, or NS, or SBH). We assume this system is isolated and there is no external force. Hence, their center of mass (CM) moves with a constant velocity and both components rotate over elliptical orbits so that these orbits have a common barycenter on their CM's location. In the CM coordinate system, their CM is fixed and the star with respect to the compact object moves over an elliptical orbit. We assume that the semi-major axis of this elliptical orbit is a . The period of this elliptical orbit is given by the third Kepler's law as follows:  \nT = 2 π √ G ( M ⋆ + M c ) a 3 / 2 . (1)  \nThis elliptical orbit is characterized as (See, e.g., Dominik 1998):  \nx = a ( cos ξ -ϵ ) , y = a sin ξ √ 1 -ϵ 2 , (2)  \nwhere, ϵ is the orbital eccentricity. ξ , the so-called eccentric anomaly, is given by Kepler's equation, i.e., ϕ = ξ -ϵ sin ξ , where ϕ = ω ( t -t p ) is called the mean anomaly, and ω = 2 π / T is the angular velocity. Also, t p is a characteristic time which indicates the time of crossing the orbital periapsis point. Hence, ( x, y ) defines the orbital plane of the star around the compact object. In the simulation, we solve Kepler's equation numerically using:  \nξ = ϕ + ∞ ∑ i =1 2 π sin( iϕ ) B i ( iϵ ) , (3)  \nwhere, B i is the known Bessel function of i th order.  \nWe define the observer's coordinate system, ( x o , y o , z o ), where x o is toward the observer, and ( y o , z o ) defines the sky plane. We then convert the CM coordinate system to the observer coordinate system using two projection angles, i.e., θ and i . θ is the angle between the minor axis of this elliptical orbit and the sky plane, and i , the so-called inclination angle, is the angle between the line of sight toward the observer and  \nthe orbital plane. Hence, the orbital components in the observer coordinate system are given by:  \nx o =cos i ( -y sin θ + x cos θ ) , y o = y cos θ + x sin θ, z o = -sin i ( -y sin θ + x cos θ ) . (4)  \nFor stellar orbits with small inclination angles (the socalled edge-on ones) the stellar light is magnified by the compact object once an orbital period and whenever the star is passing behind it. We explain the details of calculating self-lensing, occultation and eclipsing signals in the following subsection.", '2.2. Simulating Self-Lensing, Occultation, and Eclipsing Signals': "In the stellar orbit around the compact object (the lens object), there is a phase angle Φ which is the angle between two lines of sight from the source star: one toward the compact object and another one toward the observer. It is calculated by cos Φ = -x o / D , where D = √ x 2 o + y 2 o + z 2 o is the distance between two components of the binary system versus time. When the source star is passing behind the compact component, this phase angle changes from 0 to 90 · , and when the source star is passing in front of the compact component this phase angle alters from 90 to 180 · . Hence, in the simulation when Φ < 90 · (or x o < 0), we calculate the lensing effect.  \nThe lensing magnification factor depends on the source radius projected on the lens plane and normalized to the Einstein radius ρ ⋆ , and the lens-source distance projected on the lens plane and normalized to the Einstein radius u . In our formalism, these two parameters are given by  \nρ ⋆ = R ⋆ D l R E ( D l -x o ) , R E = √ 4 GM c c 2 D l | x o | D l -x o , u = d p R E , d p = √ y 2 o + z 2 o , (5)  \nwhere, D l is the distance of the compact object from the observer, R ⋆ is the radius of the luminous star, G is the gravitational constant, c is the light speed, and d p is the projected distance between the compact object and the source star on the sky plane.  \nAnother factor which impacts the lensing magnification is the limb-darkening effect for stellar surface brightness. In the simulation, we consider a linear limb-darkening profile for the source star, as I = I 0 [ 1 -Γ(1 -µ ) ] , where I 0 is the stellar brightness at the center of the source disk, Γ is the so-called limb-darkening coefficient, and  \nµ = √ 1 -R 2 / R 2 ⋆ , where R is the radial distance over the source disk. We determine the magnification factor, A ( u, ρ ⋆ , Γ ) , using the public RT-model (Bozza 2010; Bozza et al. 2018). We ignore the General Relativistic (GR) effects on the magnification factor because this effect is significant only for a very small part of the source disk which is exactly behind and collinear with the compact object as seen by the observer.  \nIn self-lensing events, the Einstein radius (Eq. 5) is estimated using R E ≃ √ 4 GM c a/c , which is considerably smaller than the Einstein radius for common microlensing events toward the Galactic bulge. Hence, in self-lensing events, ρ ⋆ ≳ 1 and their magnification factors are estimated by the ratio of the area of the images' ring generated at the time of the complete alignment to the source area (e.g., Han 2016). The inner and outer radii of the images' ring are (respectively):  \nR in = 1 2 [√ R 2 ⋆, p +4 R 2 E -R ⋆, p ] , R out = 1 2 [√ R 2 ⋆, p +4 R 2 E + R ⋆, p ] , (6)  \nwhere, R ⋆, p = R ⋆ D l / ( D l -x o ) is the projected source radius on the lens plane. Accordingly, the magnification factor for a uniform source star during a self-lensing signal is given by (Maeder 1973; Gould & Gaucherel 1996; Agol 2003):  \nA ≃ R 2 out -R 2 in R 2 ⋆, p ≃ 1 + 2 ρ 2 ⋆ . (7)  \nWe note that when x o < 0 and during self-lensing signals, the compact object (the gravitational lens) can block some part of images' area, which is the socalled finite-lens effect or occultation of images' light by the compact object (see, e.g., Marsh 2001; Han 2016). Finite-lens effect is considerable in WDMS binaries. This effect decreases the magnification factor by O = R 2 c / R 2 ⋆, p when R in ≤ R c ≤ R out , and does not change the magnification factor ( O = 0) if R c < R in . Here, R c is the radius of the compact object. If the radius of the compact object is larger than the outer radius of the images' ring, the light of images is completely blocked by the compact object and O = A .  \nWhenever x o > 0 (or Φ ∈ [90 , 180 · ]) the source star is passing in front of its compact companion and it can block its luminosity, the so-called eclipsing effect (when the compact objects are either WDs or NSs with the masses M c ≲ 2 . 9 M ⊙ ). According to our formalism, eclipsing features happen when x o > 0 and d p ≤ ( R c + R ⋆, p ) . We calculate the fraction of the compact object's disk eclipsed by the stellar compan-  \nFigure 1. Maps of log 10 [∆ F L ] in detached and edge-on binary systems over the 2D space of log 10 [ M c ( M ⊙ )] -log 10 [ a ( R ⊙ )]. We consider four different source stars in binary systems include: (i) a red dwarf with M ⋆ = 0 . 3 M ⊙ and R ⋆ = 0 . 35 R ⊙ , (ii) a Sun-like star, (iii) a B-type star with M ⋆ = 2 M ⊙ , and R ⋆ = 1 . 8 R ⊙ , and (iv) a sub-giant star with M ⋆ = 2 M ⊙ , and R ⋆ = 3 . 6 R ⊙ . The contours display log 10 [ T (days)]. The interacting binaries with D RL < R ⋆ (see, Eq. 10) are excluded and covered with white color on maps.  \n<!-- image -->  \n, y o , z o ), numerically by:  \nE ( x o , y o , z o ) = 1 πR 2 c ∫ R c -R c dy ' ∫ √ R 2 c -y ' 2 -√ R 2 c -y ' 2 dz ' Θ [ R ⋆, p d ' p ] , (8)  \nwhere d ' p = √ ( y o -y ' ) 2 +( z o -z ' ) 2 is the distance of each element over the compact object's disk ( y ' , z ' ) from the source center projected on the sky plane. Θ is a step function which is one if its argument is larger than one and it is zero when the argument is less than one. We note that the whole disk of the compact object is eclipsed (i.e., E = 1) if d p ≤ ( R ⋆, p -R c ), and there is no eclipsing ( E = 0) when either x o < 0, or x o > 0 and d p ≥ ( R c + R ⋆, p ).  \nThe apparent magnitude of the source star by considering self-lensing, occultation (or finite-lens effect), and eclipsing signals in edge-on binary systems as measured by the observer is given by:  \nm o = m ⋆ -2 . 5 log 10 [ f b A ( u, ρ ⋆ , Γ) -O + FE 1 + F +1 -f b ] (9)  \nwhere, F is the ratio of the compact object's flux to the stellar flux, and f b is the blending factor which is the ratio of the source flux to the total flux received from the source star PSF (Point Spread Function), m ⋆ is the apparent magnitude of source star when it is isolated and without any companion.", '2.3. Characteristics of Self-Lensing Signals': "Here, we evaluate the peak amounts in self-lensing signals, i.e., ∆ F L ≃ 2 ρ -2 ⋆ , and consider four types of source stars: (i) a red dwarf with M ⋆ = 0 . 3 M ⊙ , and R ⋆ = 0 . 35 R ⊙ , (ii) a Sun-like star, (iii) a B-star with M ⋆ = 2 M ⊙ , and R ⋆ = 1 . 8 R ⊙ , and (iv) a sub-giant star with M ⋆ = 2 M ⊙ , and R ⋆ = 3 . 6 R ⊙ , as well as two wide ranges for the compact object's mass and the orbital semi-major axis as M c ( M ⊙ ) ∈ [0 . 1 , 50], and a ( R ⊙ ) ∈ [5 , 5000], respectively. For all of these binary systems, we calculate ∆ F L values and show their maps in Figure 1.  \n<!-- image -->  \n<!-- image -->  \n0, O(deg) = 0,  \n<!-- image -->  \nMc  \n<!-- image -->  \n<!-- image -->  \nFigure 2. Simulated self-lensing signals versus time (normalized to the orbital period) by considering different values for M c ( M ⊙ ), ϵ , i , R ⋆ ( R ⊙ ), θ , and Γ. For each panel, several parameters are fixed and reported at the top of that panel. The vertical axis is the flux of the source star normalized to its baseline value.  \n<!-- image -->  \nOver these maps, the contours of log 10 [ T (days)] are shown with black solid and dashed lines (corresponding to positive and negative values, respectively). In these plots, binary systems with source radii larger than the Roche-Lobe distance are excluded and covered with white colour on maps. For these systems, the mass transfers from the bright star to the compact object and there are some other sources for variability (e.g., ellipsoidal effect) in stellar fluxes which makes discerning self-lensing signals hard and impossible. Here, we determine the Roche-Lobe distance from the source center as (Paczy'nski 1971; Eggleton 1983):  \nD RL = a -a 0 . 49 q 2 / 3 0 . 6 q 2 / 3 +log 10 [ 1 + q 1 / 3 ] , (10)  \nwhere, q = M c /M ⋆ . According to these plots, generally more massive compact objects in wider orbits (with larger semi-major axes) have larger lensing-induced signals, so that increasing (a) the mass of compact objects from 0 . 1 M ⊙ to 50 M ⊙ , and (b) the semi-major axis by three orders of magnitude (from 5 R ⊙ to 5000 R ⊙ ) enlarge ∆ F L by ∼ 5 orders of magnitude. However, more massive compact objects in wider orbits have longer orbital periods. For that reason, less massive source stars (e.g., red dwarfs in comparison to B-stars) are more suitable for detecting their self-lensing signals generated by their compact companions because (i) they generate shorter orbital periods, and (ii) they have smaller radii, smaller ρ ⋆ values, and as a result higher ∆ F L values. We probe this point again in the next section and through Monte Carlo simulations.  \nGenerally, peaks, shapes and durations of self-lensing signals determine their detectability. They depend on several parameters including M c , ϵ , i , R ⋆ , θ , and Γ. We therefore simulate self-lensing signals by considering different values of these parameters as shown in different panels of Figure 2. In each panel, one parameter changes and other parameters are fixed and mentioned at the top of that panel. Accordingly, we summarize some points in follows.  \n- · The maximum enhancement in self-lensing signals increases with the lens mass as ∆ F L ∝ M c . The selflensing signals due to completely edge-on orbits are flattened (top-hat models).\n- · There is a degeneracy between M c and R ⋆ , so that small stellar radii make similar self-lensing signals to those due to more massive compact objects. Both of these parameters change the peaks of self-lensing from top-hat ones (but they do not change the self-lensing edges).  \n- · If the stellar orbits are eccentric, the resulting selflensing signals are asymmetric, unless the source star is passing from either periapsis or apoapsis point of its orbit while lensing (which are rare). This point can be found from the fifth panel of Figure 2.\n- · The inclination angle of stellar orbit with respect to the line of sight is the only factor which causes that self-lensing signals at the edges are not broken (not a strict top-hat model). By increasing the inclination angle self-lensing signals at edges are rather slowenhancing.\n- · In more eccentric stellar orbits, the resulting selflensing signals could be wider depending on the sourcelens distance while lensing (in our formalism this distance is determined by θ ).\n- · The limb-darkening effect has a very small impact on the width of self-lensing signals, whereas it changes the peak and shape of self-lensing signals.  \nIn the next step, we perform Monte Carlo simulations from all possible self-lensing signals due to different binary systems (including main-sequence stars and compact objects) and study their detectability in the TESS data.", '3. MONTE CARLO SIMULATIONS': "To simulate potential self-lensing, occultation, and eclipsing signals that can be detected in the TESS observations, we first take an ensemble of the Candidate Target List (CTL, Stassun et al. (2018, 2019)) from the Mikulski Archive for Space Telescopes (MAST) catalog (STScl 2022). The TESS CTL targets are relatively close and bright stars with pre-measured physical parameters which were (and are) observed by the TESS telescope with a 2-min cadence, based on their priority. In this ensemble, for CTL targets their mass M ⋆ , priority, blending factor f b , distance D l , source radius R ⋆ , effective surface temperature T ⋆ , and stellar apparent magnitude in the TESS filter m ⋆ are reported which all are used for simulating source stars in binary systems.  \nWe assume the TESS CTL targets live in detached binaries with compact companions. We take 1772 known WDs which were extracted from the SDSS data with the reported mass M WD , distance, radius R WD , and apparent magnitudes in the Gaia passbands, i.e., m G , m G BP , and m G RP at the distances closer than 100 parsec from Kilic et al. (2020). Using the same method which was explained in Sajadian et al. (2024), we convert their apparent magnitudes in the Gaia filters to their absolute magnitudes in the TESS passband.  \nFigure 3. Six examples of simulated stellar light curves due to detached and edge-on binary WDMS, NSMS, and BHMS systems. The synthetic data points are hypothetically taken by the TESS telescope. We fix the observational time to 27 . 4 days with a one-day gap in middle. The parameters used to make light curves are mentioned at the tops of panels. By applying HC detectability criteria, for three of these light curves the impacts of compact objects are realizable, as mentioned in their legends.  \n<!-- image -->  \nThe numbers of discovered NSs and SBHs up to now are low. Additionally, all necessary parameters for discovered ones were not measured. Therefore, in Monte Carlo simulations we generate their populations synthetically according to the known distribution functions for their physical parameters.  \nThe mass distribution function of NSs is a bimodal distribution with two peaks at 1 . 37 M ⊙ , and 1 . 73 M ⊙ whereas the second one is wider and more flattened (Valentim et al. 2011). So, we determine the mass of NSs using N ( 1 . 37 , 0 . 15 ) +0 . 5 N ( 1 . 73 , 0 . 3 ) from the range M NS ∈ [0 . 5 , 2 . 9] M ⊙ , where N ( µ, σ ) is a normal function with the mean value µ , and the width σ . The radius of NSs is a function of their mass and is in the range R NS ∈ [10 , 15] km (Lattimer 2019). The luminosity of NSs is a function of their ages and masses, because these objects cool down over time through thermal radiation. We determine the NSs' luminosity based on Figure (2) of Potekhin et al. (2020). We choose the NSs' ages ( A ) in the logarithmic scale uniformly from the range log 10 [ A (year)] ∈ [2 . 5 , 8]. For SBHs we choose their masses uniformly from the range M BH ∈ [3 . 3 , 50] M ⊙ (Sicilia et al. 2022; Sajadian & Sahu 2023).  \nIn each step of Monte Carlo simulations, we take one CTL target and one compact object as two components of a detached binary system. We take the semi-major axis of their orbit a from a log-uniform distribution and in the range [3 R ⋆ , 10 6 R ⋆ ] (see, e.g., Abt 1983), and put aside the interacting systems with the source radii larger than the Roche-Lobe distance ( D RL which is given by Eq. 10).  \nThe known period-eccentricity correlation indicates an upper limit on the orbital eccentricity ( ϵ max ) for a given orbital period T (Mazeh 2008). By considering this correlation, we uniformly choose the orbital eccentricity from the range ϵ ∈ [0 , ϵ max ]. We choose the inclination angle uniformly from the range i ∈ [0 , 20 · ], because for larger inclination angles neither self-lensing/occultation nor eclipsing signals happen. We choose the projection angle θ uniformly from the range [0 , 360 · ].  \nIn the next step, we generate synthetic data points taken by the TESS telescope. We fix the observing cadence to 2 minutes, and consider a one-day gap after each 13 . 7-day observation. The observing time span is calculated by T obs = No s × 27 . 4 days, where No s is the number of overlapping sectors and can be an integer number from one to thirteen. Therefore, the longest observing time is ∼ 360 days for the ecliptic poles during one year.  \nTo examine their detectability we calculate the signalto-noise ratio (SNR), which is defined for planetary transit events, as given by (see, e.g., Fatheddin & Sajadian  \n2024):  \nSNR = √ N tran ∆F × 10 6 CDPP , (11)  \nwhere, CDPP is the Combined Differential Photometric Precision (CDPP) metric which is a function of the stellar apparent magnitude in the TESS passband, N tran = T obs / T is the number of orbital period during an observing time span. ∆ F is the maximum variation in the stellar flux which is due to either selflensing/occultation or eclipsing signals (the maximum of ∆ F L , ∆ F O , and ∆ F E ). In the simulation, to extract the detectable events we consider two sets of criteria which are (i) SNR > 5, and N tran > 2, and (ii) SNR > 3, and N tran > 1. The first set extracts detectable events with a high confidence (HC), and the second one takes detectable events with a low confidence(LC).  \nIn Figure 3, we show six examples of simulated stellar light curves due to edge-on WDMS, NSMS, BHMS systems (with i ≤ 20 · ). Some of useful parameters to make light curves and their SNR values are reported at tops of plots. The magenta synthetic data points are taken by the TESS telescope with a 2-min cadence. Accordingly, three events are detectable (with a high confidence) in the TESS data and three others are not detectable.", '4. RESULTS: STATISTICS AND PROPERTIES': "We assume the maximum observational time span for a part of the sky during the TESS mission is 360 days. Although, the southern (northern) ecliptic hemisphere was re-observed in the third (forth and fifth) year(s) of the TESS mission again, discerning periodic variations with T (day) ≥ 360 (longer than the maximum value for the continuous TESS observing time span from a part of the sky) in stellar light curves is barely possible. Because, (i) there is a 1- up to 2-year gap in the middle, and (ii) 74% of stars are observed during only 27 . 4 days of a year (they are inside one sector) and the probability of occurring either self-lensing/occultation or eclipsing signals (when the orbital period is long) exactly during that 27 . 4-day observing time is low. We therefore simulate synthetic data points for the events with T < 360 days, which means we assume all events with T ≥ 360 days are not detectable in the TESS observations.  \nWe perform three Monte Carlo simulations from WDMS, NSMS, and BHMS binary systems by considering different observing time spans, which can be from  \nTable 1. The average parameters of the simulated WDMS binaries which have detectable WD-induced signals with a high confidence in their stellar light curves by considering different observing time spans.  \nNote -No s is the number of overlapping sectors. ε HC , and ε LC are the efficiencies for detecting WD-induced signals with high and low confidences, respectively.  \nTable 2. Same as Table 1, but for simulated NSMS binaries.  \n27 . 4 days to 356 days due to different numbers of overlapping sectors. In Table 1, the results from the Monte Carlo simulation of WDMS binary systems are reported. This table includes the average values of some orbital and WDs parameters (including M WD ( M ⊙ ), T (days), log 10 [ a/R ⋆ ], ϵ , i (deg), log 10 [ F ], D l (kpc), SNR, log 10 [ ρ ⋆ ], and b ( R ⋆ )) for HC detectable events during different observing time spans (given in the second column).  \nHere, b is the impact parameter in the lensing formalism, which is the minimum value of the projected distance d p .  \nIf this parameter is less than the source radius, eclipsing/occultation signals certainly happen. Also, this parameter determines the maximum magnification factor. The thirteenth column of this table determines (i) the fraction of simulated events which are detectable because of their self-lensing signals, f L , (ii) the fraction of simulated events which are detectable owing to their eclipsing signals f E , and (iii) the fraction of ones which are detectable due to their occultation (finite-lens effect) signals f O . Two last columns report the detection  \nTable 3. Same as Tables 1 and 2, but for simulated BHMS binaries.  \nefficiencies ε HC [%] and ε LC [%] which are the ratio of detectable events (with high and low confidences, respectively) to total simulated events.  \nThe results from Monte Carlo simulations of detached binary systems including main-sequence stars and either NSs or SBHs are reported in Tables 2, and 3, respectively. We note that for BHMS binary systems f E = 0, f O = 0, and F = 0.  \nGenerally, longer observing time windows have two positive effects on the detectability of compact objects' signatures. For longer observing time windows, the number of transits N tran is higher which (i) increases SNR values (see Eq. 11), and (ii) enhances N tran (the second detectability criterion). Enhancing N tran (and accordingly SNR) for longer observing times is beneficial for detecting fainter compact objects (WDs and NSs which could be even farther), in wider and more eccentric orbits. For that reason in these three tables for longer observing times detectable compact objects are on average fainter (with less log 10 [ F ]) and farther. We note that due to applying the eccentricity-orbital period correlation in Monte Carlo simulations binary systems with longer orbital periods are on average more eccentric. By increasing the observing time from 27 . 4 days to 356 . 2 days, the detection efficiency improves by ∼ 2.  \nFor ≲ 3% and ≲ 33% of detectable WDMS and NSMS binary systems, self-lensing signals are the mostdominant ones in stellar light curves. Indeed, for these detached WDMS and NSMS binary systems on average ρ ⋆ ∼ 145 , 100 which make very flattened self-lensing signals with the depths ∆ F L ∼ 2 ρ -2 ⋆ ∼ 10 -4 , 2 × 10 -4 , respectively, whereas their eclipsing signals, as given by ∆ F E ∼ F / (1 + F ) ∼ 1-2 × 10 -3 , 10 -5 -5 × 10 -4 , respec-  \nor detectable WDMS and NSMS binary systems we have ∆ F L ≲ ∆ F E . The occultation signal in WDMS binary systems is on average ∆ F O ∼ R 2 c /R 2 ⋆ ∼ 10 -4 which is in the same order of magnitude with the self-lensing ones. Here, we consider a common WD with the radius R WD ∼ 0 . 01 R ⊙ .  \nTo study what kinds of binary systems and compact objects are more detectable in the TESS observations, in Figure 4 we show the efficiency curves for detecting the impacts of compact objects with a high confidence, ε HC , in simulated stellar light curves versus nine parameters, which are D l (kpc), log 10 [Periority], log 10 [ a/R ⋆ ], M ⋆ ( M ⊙ ), log 10 [ ρ ⋆ ], i (deg), log 10 [ b/R ⋆ ], M NS ( M ⊙ ), and M BH ( M ⊙ ). We consider three amounts for T obs as mentioned in the first panel.  \nAccording to these plots, we conclude that closer source stars which have more priorities and higher detection efficiencies. Indeed, WDs and NSs in nearby binary systems are brighter with on average higher F values. Higher F values make deeper eclipsing signals.  \nClose binary systems with smaller semi-major axes have on average shorter orbital periods, which are more suitable to be detected. Because short-period binary systems have higher N tran and higher SNR values. We note that ρ ⋆ ∝ a -1 / 2 . Therefore, although the closer binary systems have shorter orbital periods (and higher SNR values), they have larger ρ ⋆ s, and as a result more flattened self-lensing signals. Hence, during longer observing times wider binary systems can be detected rather via self-lensing ( f L increases with T obs in Tables 1, and 2).  \nThe impact of stellar masses on the self-lensing signals has been shown in Figure 1. Accordingly, less massive  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 4. The HC efficiencies for detecting signals due to compact objects in stellar light curves (in the logarithmic scale log 10 [ ε HC (%)]) versus nine parameters which are D l (kpc), log 10 [Periority], log 10 [ a/R ⋆ ], M ⋆ ( M ⊙ ), log 10 [ ρ ⋆ ], i · , log 10 [ b/R ⋆ ], M NS ( M ⊙ ), and M BH ( M ⊙ ) as shown in different panels. Three efficiency curves with different line-styles are due to three values for T obs as mentioned in the first panel.  \n<!-- image -->  \nstars (with smaller radii) have higher self-lensing signals. Also, lower mass stars are on average fainter which results in higher F and deeper eclipsing signals.  \nStellar orbits with fewer inclination angles (more edgeon ones) have higher detection efficiencies because by decreasing the inclination angle the impact parameter reduces as well. The eclipsing/occultation signals can be detected only in binary systems with impact parameters less than source radii. The detection efficiency drops from ∼ 4 -5% to ≲ 1% when the inclination angle increases from 0 to 15 degrees.  \nWe also plot the detection efficiencies as a function of the mass of NSs and SBHs in two last panels of Figure 4. We note that the mass range for WDs is small, and the detection efficiency does not highly change with the mass of WDs. For NSs, we determine their luminosity according to their mass and age (based on Fig. (2) of Potekhin et al. 2020). Accordingly, by increas-  \ning the mass of NSs by ∼ 1 M ⊙ the luminosity of NSs decreases up to three orders of magnitude (when they are younger than 3 million years). Hence, on average less massive NSs are brighter with higher F values, and deeper eclipsing signals. For detached SBHs, the only method to detect them is self-lensing. In self-lensing formalism, we have ∆ F L ∼ 2 ρ -2 ⋆ ∝ M BH . Therefore, more massive SBHs make higher self-lensing signals with higher SNR values.  \nTop panel of Figure 5 shows the scatter plot of simulated WDMS binaries in 2D space log 10 [ F ] -log 10 [ b/R ⋆ ] with black circles, with two marginal and normalized 1D distributions. The binary systems with detectable WDinduced impacts due to eclipsing, self-lensing, and occultation effects are specified with red, blue, and green circles, respectively. Accordingly, most of WDMS binaries with detectable WDs impacts have b ≲ R ⋆ and F ≳ 10 -4 . The blue points (with detectable self-lensing  \nFigure 5. Top panel: The scatter plot of all simulated WDMS binaries (black points) and the ones with detectable WD-induced impacts (colored points) in 2D space log 10 [ F ] -log 10 [ b/R ⋆ ]. Their 1D and normalized distributions are shown at two sides of the plot. The red, blue, and green circles represent detectable events owing to their eclipsing, self-lensing, and occultation signals, respectively. Middle panel: Same as the previous one but resulted from the Monte Carlo simulation of NSMS binaries. Last panel: Same as two previous ones but resulted from simulating BHMS binaries and in 2D space log 10 [ T (days)] -log 10 [b / R ⋆ ].  \n<!-- image -->  \nsignals) have lower F up to 10 -5 . We therefore expect in the TESS data eclipsing-induced footprints due to WDs to be more realizable than their self-lensing/occultation impacts.  \nThe next panel of Figure 5 shows a scatter plot the same as the previous one but for NSMS binary systems. Considering this point that we chose the ages of NSs uniformly from the range log 10 [ A (year)] ∈ [2 . 5 , 8], several NSs in our simulation are cool with very low F values (ones older than ∼ 1 million years). For these systems with F ≲ 10 -4 , only lensing-induced impacts are detectable (and when b ≲ R ⋆ ). For brighter (younger and hotter) NSs, eclipsing signals are detectable in the events with b ≲ R ⋆ . We note that due to this considerable number of dim NSs in our simulation, only f L ∼ 15-33% of detectable events have dominant self-lensing signals.  \nThe last panel of Figure 5 shows the same plot as ones displayed in previous panels but for BHMS binary systems in 2D space log 10 [ T (days)] -log 10 [b / R ⋆ ]. The most important factor for the detectability of BHinduced lensing signals is the impact parameter b . Most of the simulated events with b ≲ R ⋆ have detectable self-lensing signals. We determine the impact parameter numerically from the simulation (the minimum value of d p ), nevertheless it can be estimated as b ≃ tan( i ) a . Therefore, BHMS binary systems with smaller semimajor axes have smaller b and shorter orbital periods which both impacts are beneficial for detecting SBHs through self-lensing signals.  \nHere, we estimate the number of WDs that the TESS telescope can detect through precise photometric observations from the CTL targets with a 2-min cadence during its mission. The TESS telescope is planned to detect 1 , 390 , 486 CTL targets during its mission 3 (STScl 2022). We calculated the fractions of these CTL targets which are observed during different T obs due to different No s values, i.e., N ⋆ ( No s ), as reported in the seventh column of Table (2) of Sajadian et al. (2024) and extracted from Fig (2) of Barclay et al. (2018). These numbers are N ⋆ ( No s ) = (1031 , 258, 38, 9, 5, 5, 3, 2, 1, 1, 1, 16, 13) × 1000.  \nThe number of WDMS binaries in our galaxy is predicted to be ∼ 100 million. Considering the total number of stars in our galaxy (which is 100 billion), one thousandth ( f 1 ) of the TESS CTL targets should be in binary systems with WDs. We note that in the simulation we limit the inclination angle to i ∈ [0 , 20 · ], whereas f 2 = 22% of all binary systems should have i ≤ 20 · . We assume that the inclination angle of the  \norbital plane is uniformly in the range [0 , 90 · ]. By assuming that the TESS CTL targets make a common sample of stars, the number of WDs detectable in the TESS observations from the CTL targets can be estimated by:  \nN WD ≃ f 1 × f 2 × 13 ∑ No s =1 N ⋆ ( No s ) × ε i ( No s ) , (12)  \nwhere, values of ε i ( No s ) (here, i = HC, and LC) are given in two last columns of Table 1. The estimated numbers of detectable WDs for different observing time windows and two levels of confidence are reported in two first rows of Table 4.  \nOverly, we expect ∼ 15 , 18 WDs in detached WDMS systems can be realized with high and low confidences through the TESS photometric data from the CTL targets. We note that some of these stars were (and will be) observed more than one year (two or three times) during the TESS mission. For instance, the TESS telescope observed each ecliptic hemisphere two times up to now. But, in the Monte Carlo simulations, we consider the TESS data for each CTL target taken for up to one year. Hence, these numbers are underestimations and by considering all data some more WDs in detached orbits around main-sequence stars with days will be discovered.  \nIn the same way, we estimate the numbers of NSs and SBHs expected to be discovered through detecting their self-lensing and eclipsing signals in the TESS CTL's light curves which are reported in Table 4. The fractions of the TESS CTL targets that are in binary systems with NSs and SBHs are f 1 ≃ 4 × 10 -4 , 4 × 10 -6 , respectively. Here, we assume the binarity fractions of NSs, and SBHs (with main-sequence stars) are 4% (Tauris & van den Heuvel 2006). The total numbers of NSs and SBHs will be discovered with two levels of confidence through photometric observations of the TESS CTL targets are ∼ 6 , 7, and less than one, respectively.  \nTherefore, changing detectability criteria from a high confidence to a low one does not alter the numbers of detectable compact objects significantly. Because, according to Tables 1, 2, and 3, HC and LC efficiencies for detecting impacts of compact objects in WDMS, NSMS, and BHMS systems are close to each other, and ε HC = 4 . 6 -7 . 0% , 4 . 5 -6 . 4% , 4 . 2 -5 . 3%, and ε LC = 5 . 6 -7 . 8% , 5 . 6 -7 . 0% , 4 . 3 -5 . 8%. In fact, the fractions of simulated WDMS, NSMS, and BHMS binary systems with T < 360 days are 35 . 3%, 37 . 3%, and 42 . 8%, respectively, and the corresponding fractions for systems with T < 180 days are 32 . 1%, 33 . 2%, and  \n39 . 2%. Therefore, more than ∼ 60% of simulated events have T > 360 days, and the numbers of binary systems with orbital periods ∈ [180 , 360] days are not significant.  \nAlthough the number of SBHs that can be detected through self-lensing is not promising, we expect ∼ 1518 WDs to be detected through their signatures in light curves of the TESS CTL targets. This number is three times larger than the number of WDs discovered from the Kepler data.", '5. CONCLUSIONS': "The TESS telescope observed (and observes) the CTL targets with a 2-min cadence and a great accuracy. Although its main goal from these observations is discovering Earth-size planets transiting bright stars in the solar neighborhood, its observing strategy is uniquely matched to capture any other types of periodic and weak variations in stellar light curves. Stellar light curves from edge-on and detached binary systems including main-sequence stars and compact objects have potentially self-lensing/occultation/eclipsing signals which all are periodic. Considering three types of compact objects (i.e., WDs, NSs, and SBHs), these binary systems usually are denoted by WDMS, NSMS, and BHMS, respectively. In this work, we studied statistics and properties of WDMS, NSMS, and BHMS binary systems with detectable signatures due to compact objects in the TESS observations of the CTL targets.  \nWe first modeled self-lensing signals due to detached and edge-on binary systems including main-sequence stars and compact objects. The self-lensing peak, which can be estimated by ∆ F L ∼ 2 ρ -2 ⋆ , is higher for low-mass and small stars. A self-lensing signal for a red-dwarf star with R ⋆ = 0 . 35 R ⊙ is higher (two orders of magnitude) than that for a sub-giant star with R ⋆ = 3 . 6 R ⊙ . Peaks of self-lensing signals are degenerate functions of two parameters (i) the mass of the compact object and (ii) the source radius. Increasing the mass of compact objects and decreasing the source radius have the same effects on peaks of self-lensing signals. The inclination angle of stellar orbits around the compact objects is the only parameter which changes the shape of self-lensing signals from strict top-hat models to ones with slow-increasing edges. Self-lensing signals from eccentric stellar orbits are asymmetric concerning their peaks unless the source star is passing from either apoapsis or periapsis point while lensing.  \nWe performed Monte-Carlo simulations from all possible stellar light curves due to detached and edge-on binary WDMS, NSMS, and BHMS systems, and assumed they are observed by the TESS telescope with a 2-min cadence. We chose source stars in these simu-  \n.  \nTable 4. Estimated numbers of WDs, NSs, and SBHs with detectable impacts in the TESS CTL's light curves, by considering different observing time spans, and two levels of confidence (HC, LC).  \nlations from the TESS CTL, WDs from an ensemble of 1772 discovered nearby ones through the SDSS observations, NSs, and SBHs from their known distributions. We made synthetic data points according to the TESS observing strategy for the CTL targets and extracted the ones with the detectable signatures due to compact objects based on two sets of criteria (i) SNR > 5 and N tran > 2, i.e., detecting with a high confidence(HC), and (ii) SNR > 3 and N tran > 1, i.e., detecting with a low confidence (LC).  \nThere are two issues for detecting (at least one) selflensing or eclipsing signal with the period longer than 360 days from the TESS data which are: (i) the longest (continuous) observing time span for a part of the sky is 360 days (happens for the ecliptic poles), and there is a gap (from one up to two years depending on the locations) between one-year missions, and (ii) around 74% of the CTL targets are observed during only 27 . 4 days of a one-year mission. We therefore assumed that the TESS detection efficiency for T ≥ 360 days is zero, although it is actually close to zero.  \nWe found that the probability that simulated binary WDMB, NSMS, and BHMS systems have T < 180(360) days were 32 . 1(35 . 3)%, 33 . 2(37 . 3)%, and 39 . 2(42 . 8)%, respectively. The HC (and LC) efficiencies for detecting periodic signatures from WDMS, NSMS, and BHMS binaries during different observing time spans of a one-year mission of TESS, (i.e., T obs = 27 . 4 × No s , where No s = 1 , 2 , 3 , ..., 13) are 4 . 6 -7 . 0(5 . 6 -7 . 8)%, 4 . 5 -6 . 4(5 . 6 -7 . 0)%, and 4 . 2 -5 . 3(4 . 3 -5 . 8)%, respectively. Increasing observing time spans improves the detection efficiency by ∼ 2% and fainter compact objects in wider and more eccentric orbits can be detected during longer observing time spans (see Tables 1, 2, and 3). The fractions of detectable WDMS and NSMS bi-  \ny events which their self-lensing signals are the mostdominant ones are ≲ 3%, and ≲ 33%, respectively.  \nWe found that the detection efficiency is higher for closer CTL targets with higher priorities and smaller radii. Most of binary systems with detectable periodic signals have the impact parameters b ≲ R ⋆ . The detection efficiencies for detecting more massive NSs, and SBHs is higher.  \nWe estimated the total number of WDs, NSs, and SBHs that can be discovered from the TESS CTL observations which are 15 -18, 6 -7, and less than one, respectively. The number of detectable WDs in the TESS data is three times higher than the number of WDs discovered in the Kepler data.  \nAll simulations that have been done for this paper are available at: https://github.com/SSajadian54/ SelfLensing Eclipsing simulator. The codes, and several examples of generated light curves can be found in the Zenodo repository(sajadian 2024).  \nIn Monte-Carlo simulations, we use the TESS CTL (with DOI number: doi:10.17909/fwdt-2x66) that are publicly available from the MAST catalog. Funding for the TESS mission is provided by NASA's Science Mission directorate. We acknowledge the use of TESS Alert data, which is currently in a beta test phase, from pipelines at the TESS Science Office and at the TESS Science Processing Operations Center. NA is supported by the University of Waterloo, Natural Sciences and Engineering Research Council of Canada (NSERC) and the Perimeter Institute for Theoretical Physics. Research at Perimeter Institute is supported in part by the Government of Canada through the Department of Innovation, Science and Economic Development Canada and by the Province of Ontario through the Ministry of Colleges and Universities.", 'REFERENCES': 'STScl. 2022, TESS Input Catalog and Candidate Target List, Version: 8.2, MAST, doi: 10.17909/fwdt-2x66  \n- Tauris, T. M., & van den Heuvel, E. P. J. 2006, in Compact stellar X-ray sources, Vol. 39, 623-665, doi: 10.48550/arXiv.astro-ph/0303456\n- Thompson, T. A., Kochanek, C. S., Stanek, K. Z., et al. 2019, Science, 366, 637, doi: 10.1126/science.aau4005  \nValentim, R., Rangel, E., & Horvath, J. E. 2011, MNRAS, 414, 1427, doi: 10.1111/j.1365-2966.2011.18477.x  \nWiktorowicz, G., Middleton, M., Khan, N., et al. 2021, MNRAS, 507, 374, doi: 10.1093/mnras/stab2135 Willems, B., & Kolb, U. 2004, A&A, 419, 1057, doi: 10.1051/0004-6361:20040085 York, D. G., Adelman, J., Anderson, John E., J., et al. 2000, AJ, 120, 1579, doi: 10.1086/301513  \n- Zhang, L., Ridolfi, A., Blumer, H., et al. 2022, ApJL, 934, L21, doi: 10.3847/2041-8213/ac81c3'}