Next Article in Journal
New Insight into the FS CMa System MWC 645 from Near-Infrared and Optical Spectroscopy
Previous Article in Journal
Hawking Radiation and Lifetime of Primordial Black Holes in Braneworld
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating Gravitationally Lensed Quasars Observable by Nancy Grace Roman Space Telescope

by
Lindita Hamolli
1,*,†,
Mimoza Hafizi
1,†,
Francesco De Paolis
2,3,† and
Esmeralda Guliqani
4
1
Department of Physics, University of Tirana, Boulevard “Zogu I”, 1001 Tirana, Albania
2
Department of Mathematics and Physics “Ennio De Giorgi”, University of Salento, Via Arnesano, 73100 Lecce, Italy
3
INFN, Sezione di Lecce, Via Arnesano, 73100 Lecce, Italy
4
Department of Mathematics and Physics, University “Fan S. Noli”, 7001 Korca, Albania
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Galaxies 2023, 11(3), 71; https://doi.org/10.3390/galaxies11030071
Submission received: 22 April 2023 / Revised: 17 May 2023 / Accepted: 22 May 2023 / Published: 1 June 2023

Abstract

:
In this work, we investigate the possibility of observing quasars, particularly lensed quasars, by the Nancy Grace Roman Space Telescope (Roman). To this aim, based on the capabilities of the Roman Space Telescope and the results from the quasar luminosity function (QLF) in the infrared band of the Spitzer Space Telescope imaging survey, we calculated the number of quasars expected to be in its field of view. In order to estimate the number of lensed quasars, we develop a Monte Carlo simulation to estimate the probability that a quasar is lensed once or more times by foreground galaxies. Using the mass–luminosity distribution function of galaxies and the redshift distributions of galaxies and quasars, we find that 1 per 180 observed quasars will be lensed by foreground galaxies. Further on, adopting a singular isothermal sphere (SIS) model for lens galaxies, we calculate the time delay between lensed images for single and multiple lensing systems and present their distributions. We emphasize that detailed studies of these lensing systems will provide a powerful probe of the physical properties of quasars and may allow testing the mass distribution models of galaxies in addition to being extremely helpful for constraining the cosmological parameters.

1. Introduction

Gravitational lensing refers to the deflection of light rays from a background source due to the presence of a massive object close to the line of sight. In the case of strong lensing, a source (quasar) will appear multiply imaged, magnified, or deformed into an arc-like shape around a lens (foreground galaxy). Since gravitational lensing depends only on the distribution and geometry of foreground matter, it is an ideal tool for mapping the distribution of dark matter and testing the values of cosmological parameters, such as Hubble constant H 0 [1].
Today, there is a well-known discrepancy in the reported measurements of H 0 . From one hand, the Planck telescope finds H 0 = 67.4 ± 0.5 km s 1 Mpc 1 by observing the cosmic microwave background (CMB) [2], which is in good agreement with the measurements from galaxy clustering and weak lensing ( H 0 = 67.4 ± 1.2 km s 1 Mpc 1 [3]). On the other hand, the S H 0 E S team, using the distance ladder method with type I a supernovae and Cepheids, reported a value H 0 = 74.0 ± 1.4 km s 1 Mpc 1 [4], which substantially differs from the Planck ( Λ C D M ) result. Another possibility for determining H 0 is provided by gravitational lensing through the time-delay method. The H0LiCOW collaboration, measuring the time delay of six lensed quasars, found H 0 = 73 . 3 1.8 + 1.7 km s 1 Mpc 1 , with an uncertainty of only 2.4 % [5]. This method is completely independent of both the supernovae and C M B analyses. Different ground- and space-based surveys, such as the Cosmic Lens All Sky Survey (CLASS [6]), the Sloan Digital Sky Survey (SDSS [7]) and Hubble Space Telescope (HST [8]) provide huge amount of photometric and spectroscopic data, which greatly helped to find gravitationally lensed quasars. These systems can be double systems, such as the first gravitational lens system, Q0957+561, discovered accidentally by [9] in 1979, triple systems or even richer systems. To date, more than 200 lensed quasars are known1, and this number is expected to grow rapidly by forthcoming surveys, including the Vera Rubin Observatory’s Legacy Survey of Space and Time (LSST), Euclid [10] and Nancy Grace Roman Space Telescope [11].
The Roman Space Telescope, formerly known as the Wide Field InfraRed Survey (WFIRST [12]), is a NASA infrared space telescope in development. It is planned to be launched in May 2027, and will be positioned at the second Lagrange point ( L 2 ) of the Earth–Sun system. Its primary mission lifetime is set for 5 years, with the possibility of extending it for an additional 5 years. One of the primary mission objectives is to investigate the nature of dark energy with a variety of methods. This will be achieved through two instruments, the wide field instrument (WFI) and the coronagraph instrument. The WFI field of view is 0.281 square degrees and carries 8 science filters ( F 062 , F 087 , F 106 , F 129 , F 158 , F 184 , F 146 , and F 213 ) with overlapping band passes spanning 0.48–2.3 µm. The Roman Space Telescope survey is intended to cover 2000 square degrees over a five-year period with NIR imaging and spectroscopy. It has an angular resolution of about 0.11 arcsec.
The goal of the present work is to examine the capability of the Roman Space Telescope to observe quasars and, in particular, quasars that are lensed by foreground galaxies. Quasars are very powerful and bright distant active galactic nuclei (AGNs). Being very far and with extremely high intrinsic luminosity, they are considered to be important instruments for understanding the content and the history of the universe. Currently, the theoretical work is focused on systems lensed by one galaxy, that is, characterized by a single lens plane, but also the gravitational lensing by two foreground galaxies2 located at different redshifts and separation is studied by many authors (see [13,14,15] for details). In these papers, by modeling the lensing galaxies through the singular isothermal sphere ( S I S ) law, the number of images and their position, as well as the formation of the Einstein rings, were analyzed. Moreover, it was found that about 5 % of lensed quasars are caused by multiple systems.

2. Quasars Observed by Roman Space Telescope

In this section, we examine the possibility of the Roman Space Telescope to observe quasars. Roman WFI will observe in infrared wavelengths (0.48–2.3 µm), and its survey area is planned to be about 2000 deg 2 . We use the quasar luminosity function (QLF), which is the comoving number density of quasars as a function of their luminosity and redshift, to find the observable number of quasars. They are luminous in almost all accessible bands, but here, we consider the QFL in the IR band, derived by [16]. Based on the Spitzer Space Telescope imaging surveys, a double power-law function of the form
ϕ ( L , z ) = d ϕ d l o g L = ϕ * L L * γ 1 + L L * γ 2 ,
was extracted, where ϕ is the comoving space density of the AGN, ϕ * is the characteristic space density (both in units of comoving Mpc 3 ), L is the rest-frame luminosity at 5 µm, and L * is the break luminosity at 5 µm, both in units of erg s 1 Hz 1 . The evolution in L * is a cubic expression (see, e.g., [17]):
l o g 10 L * ( z ) = l o g 10 L 0 * + k 1 ϵ + k 2 ϵ 2 + k 3 ϵ 3 ,
where ϵ = l o g 10 ( ( z + 1 ) / ( z + z r e f ) ) , L 0 * is a free parameter in the fit, which corresponds to the break luminosity at redshift zero, γ 1 is the faint-end slope, γ 2 is the bright-end slope and z r e f is fixed at 2.5 . The coefficients k 1 , k 2 , k 3 , γ 1 , γ 2 , ϕ * , and L 0 * are given in Table 1 by Lacy et al. [16].
The number of quasars in the redshift range z 1 z 2 is obviously given by
N = z 1 z 2 L m i n L m a x ϕ ( L , z ) d L d V c
where d V c = D H ( 1 + z ) 2 D A 2 E ( z ) d Ω d z is the comoving volume element in the solid angle d Ω and redshift interval d z , D H = c H 0 is the Hubble distance, and D A is the angular diameter distance at redshift z (see [18] for details).
The Roman Space Telescope is planned to have eight filters, but in our calculation we consider the F213 filter (the band 1.95–2.30 µm with midpoint at 2.125 µm), which is closer to the band for which the QFL parameters are extracted. In Table 1, the AB magnitude limits achieved at 5 σ in 1 h or 55 s of integration, both for point sources and for compact galaxies with half-light radius of 0.3″, are given3. Using the Vega - AB Magnitude Conversion tool, we convert the AB magnitude (e.g., m A B = 26.2 ) to the apparent magnitude ( m = 24.35 ) and then find L m i n = L 10 20.42 l o g 23.04 × 10 24 D L 2 , where D L is the luminosity distance of the source. Assuming L m a x = 10 47 erg s 1 in Equation (3), we find the number of quasars expected to be observed by the Roman Space Telescope (survey area, 2000 deg 2 ). Our results for different time intervals (1 h or 55 s) are shown in Table 1. Moreover, we give the limit of the redshift z m a x for each case. As one can see, for a simple exposure, the number of quasars is about 20 times smaller than for deep surveys.

3. Review of Strong Lensing

In strong gravitational lensing, the distances to the lens ( D L ) and source ( D s ) are very large. By considering the “thin” lens model, where the mass of the lens is represented as a two-dimensional mass sheet lying perpendicular to the line of sight (known as the lens plane), the lens equation is defined as
β = θ α ( θ ) ,
where β is the angular position of source, α ( θ ) = D L S D S α ^ ( θ ) is the scaled deflection angle and α ^ is the deflection angle4 (see [19] for details). The lens equation is, in general, a non-linear equation and can be characterized by multiple solutions of θ for a given source position β . When β = 0 , the lens and the source are perfectly aligned, and the image is given by the Einstein ring radius θ E , which results to be
θ E = 4 G M ( θ ) c 2 D L S D S D L ,
where D S , D L , and D L S are the observer–source, observer–lens and lens–source angular diameter distances, respectively, which can be calculated by the expressions [18],
D L = c H 0 ( 1 + z L ) 0 z L d z E ( z ) ; D S = c H 0 ( 1 + z S ) 0 z S d z E ( z ) ;
D L S = D S 1 + z L 1 + z S D L ; w i t h E ( z ) = Ω m ( 1 + z ) 3 + Ω k ( 1 + z ) 2 + Ω Λ .
Here, H 0 is the Hubble constant and Ω m , Ω k , and Ω Λ are the dimensionless density parameters, i.e., the sum of the cold dark matter and baryonic matter, the space curvature, and the dark energy, respectively.
The scaled deflection angle is related to the scaled lens potential by α ( θ ) = Ψ ( θ ) , and the lens Equation (4) takes the form
( θ β ) θ Ψ = 0 .
Another important quantity in gravitational lensing is the scaled arrival time defined as [19]
τ ( θ ) = 1 2 ( θ β ) 2 Ψ ( θ ) .
Light rays emitted by a distant source will take some time to reach the observer. This amount of time depends on the path followed by the light rays and will be different if the source is gravitationally lensed or not. Light rays are delayed because of two effects. First, a light ray that is bent is subject to a longer path and therefore needs more time to propagate. This gives rise to a geometrical time delay Δ t g e o m , which is given by
Δ t g e o m = 1 + z L c D L D S D L S [ 1 2 ( θ β ) 2 ] .
Second, when light rays propagate through a gravitational potential, this results, in general relativity, in an induced time-dilation effect known as the Shapiro effect [20]. The gravitational time delay caused by the lens with redshift z L is given by
Δ t g r a v = 1 + z L c Ψ ^ ( ξ ) + c o n s t .
Here, ξ is the position where the light rays cross the lens plane, which corresponds to the two-dimensional angular position θ = ξ / D L . The total time delay, Δ t , is the sum of both of them i.e.,
Δ t ( θ ) = 1 + z L c D L D S D L S 1 2 ( θ β ) 2 Ψ ( θ ) + c o n s t .
In fact, the unlensed source is not observable, so the time delay Δ t cannot be measured. However, if the source is time variable, an observable quantity is the time delay between images, which, for two images, labeled as ( 1 ) and ( 2 ) , can be written as
Δ t 1 , 2 = 1 + z L c D L D S D L S 1 2 ( θ 1 β ) 2 1 2 ( θ 2 β ) 2 Ψ ( θ 1 ) + Ψ ( θ 2 ) .
The factor outside the square bracket, noted as 1 + z L c D L D S D L S = D Δ t c (where D Δ t is the so called time-delay distance) is inversely proportional to H 0 (see Equations (6) and (7)), while the term in the square bracket depends on the lens configuration and its mass distribution. Of course, observations may allow defining only the positions of image angles ( θ 1 , θ 2 ), whereas the lens potential Ψ ( θ ) and the source position β must be inferred via lens modeling.
In the subsequent Section 3.1 and Section 3.2, we summarize the gravitational lensing principles in both single and multiple lensing systems, considering a singular isothermal sphere (SIS) model for lensing galaxies. This model corresponds to a distribution of self-gravitating particles, with a Maxwellian velocity distribution at all radii and one-dimensional velocity dispersion σ S I S . The density distribution is described by [19]
ρ ( r ) = σ S I S 2 2 π G r 2 ,
and the mass of the lensing galaxy inside a radius θ E from the galactic center in the galactic plane is given by
M ( θ E ) = π G σ S I S 2 θ E D L .
Replacing the value above in the lens Equation (5), one can obtain
θ E = 4 π σ S I S 2 c 2 D L S D S .
There are some conclusions about the relation between σ S I S and the stellar velocity distribution, σ . We remark that Ofek et al. [21] predicted σ S I S that is not exactly equal to the observed stellar velocity dispersion ( σ ) , which measures the stellar kinematics related to the gravitational potential of the system, including both the stellar and dark matter components. They proposed that σ S I S = f e σ , with ( 0.8 ) 1 / 2 < f e < ( 1.2 ) 1 / 2 . Additionally, Bolton et al. [22] conducted a study involving 53 massive early-type strong gravitational lens galaxies and found that the ratio between the stellar velocity dispersion and the velocity dispersion obtained from lensing, assuming an isothermal halo model, is approximately one. So, our assumption to calculate the Einstein radius from Equation (16), by substituting the value of σ for σ S I S , appears acceptable.

3.1. Single Lensing Systems

Adopting the SIS model to describe the lens galaxy in the single lensing systems, one can find that α ( θ ) is independent of θ and is given by the Einstein radius θ E . So, in terms of the Einstein angle θ E , one can write α ( θ ) = θ E | | θ | | θ . Then, the lens Equation (4) takes the form
θ β = θ E | | θ | | θ
Multiple images are obtained only if the source lies inside the Einstein ring, i.e., if β < θ E . When this condition is satisfied, the lens equation has two solutions:
θ 1 = β + θ E , θ 2 = β θ E .
The images at θ 1 and θ 2 , are on opposite sides of the lens galaxy and together with the source lie on a straight line. The separation between the formed images always remains equal to twice the Einstein ring’s radius, Δ θ = 2 θ E . Since α ( θ ) = Ψ ( θ ) , the scaled lens potential can be written as
Ψ ( θ ) = D L S D S 4 π σ 2 c 2 | | θ | | = θ E | | θ | | ,
and total time delay can be given by
Δ t ( θ ) = 1 + z L c D L D S D L S 1 2 ( θ β ) 2 θ E | | θ | | + c o n s t .
For a time-variable source, the time delay between two lensed images is given by
Δ t 1 , 2 = 1 + z L 2 c D L D S D L S ( θ 1 2 θ 2 2 ) .

3.2. Multiple Lensing Systems

In multiple lensing systems, there are two foreground galaxies along the line of sight towards the source. Their redshifts positioned at z 1 < z 2 are smaller than the source redshift, z s . Defining as the “optical axis” the line that passes through the mass center of the foreground galaxy, the lens equation can be written as
x S = x 1 α 1 ( x 1 ) β α 2 ( x 2 ) ,
where x S is the unlensed angular position of the source. α 1 ( x 1 ) and α 2 ( x 2 ) are the physical deflections of a ray when it passes through x 1 and x 2 , respectively. The angular position on plane of the second galaxy is found by x 2 = x 1 β α 1 ( x 1 ) , and β is defined as β D 12 D 3 D 2 D 2 S 5. Here, D 1 , D 2 , and D S are the angular diameter distances between the observer and first lens, second lens and source, respectively, while D 12 and D 2 S are the angular diameter distances between both galaxies and between the second galaxy and source (see [23] for detail). Considering the SIS model for lenses, Equation (22) indicates that up to six images of the background source can potentially be produced. In multiple lensing systems, the time delay for the light ray following a deflected light path relative to the undeflected path is given by the sum of delays during the two planes,
Δ t = 1 + z 1 c D 1 D 2 D 12 1 2 ( x 1 x 2 ) 2 β Ψ ( x 1 ) + 1 + z 2 c D 2 D S D 2 S 1 2 ( x 2 x S ) 2 Ψ ( x 2 )
where z 1 and z 2 are, respectively, the redshifts of the foreground and background lenses [24].

4. Quasar Lensing with Roman Space Telescope

In this section, we investigate lensed quasars by foreground galaxies. These events are very important for investigations in cosmology. Firstly, we describe an algorithm to estimate the number of lensed quasars that would be observed by the Roman Space Telescope and then give our results.

4.1. Simulations

In order to calculate the probability to detect the lensed quasars, we use the Monte Carlo method (see [25] for a more detailed description of our algorithm). For each event, we generate the redshift of the quasar [26] and the redshift of the galaxy [27], provided that it is smaller than the quasar redshift, extract the galaxy mass using the stellar mass function [28], and find its velocity dispersion from the relation between the galaxy’s stellar mass and stellar velocity distribution [29]. Using the SIS model for galaxies, we define the Einstein angle θ E for each quasar/galaxy pair. For a single galaxy, the probability to reside inside the Einstein angle about the observer–quasar direction would scale as θ E 2 / 4 . Since [27] give the redshift distribution of 7000 galaxies, we normalize the probability considering the whole number of galaxies, 200 billions (0.25 galaxies per arcsec 2 , [30]), and find to be 10 8 θ E 2 / 14 . We compare this result with a number n, uniformly distributed in the interval ( 0 , 1 ) , which is extracted by the Monte Carlo code. We keep this pair when its probability is smaller; otherwise, we reject it. The procedure is repeated for all expected quasars to be observed by the Roman Space Telescope. For each aligned system, we calculate the positions of the images and then the time delay between them. In order to estimate the probability that a quasar is lensed by one or more foreground galaxies, we compare the Einstein angle of the system with the angular accuracy of the considered telescope.
Let us give some more details about the distributions mentioned above:
– The quasar redshift distribution is taken by Ref. [26], where the authors determined the redshift distribution of 46,420 quasars in the SDSS Data Catalog (Third Data Release), covering an area of 4188 deg 2 . Most quasars have z < 2 , with a median value of 1.47 and a mode of 1.85 . There are 520 quasars at z > 4 , including 17 at z > 5 , and 69 with z < 0.15 .
– The galaxy redshift distribution was found by [27]. The authors analyzed a sample of 7000 distant galaxies observed near the South Galactic Pole by the FORS Deep Field (ESO VLT), covering a sky area of about 7 × 7 arcmin 2 . This distribution was also confirmed by the analysis of Davison et al. [28].
– The galaxy mass distribution obtained by Davidzon et al. [28] is used to extract the mass of the lensing galaxies. They utilized the COSMOS2015 catalog to provide a comprehensive view of the galaxy stellar mass assembly in the redshift range between z = 0.1 and 6. Measurements were fitted with a double Schechter function up to z = 3 and a single Schechter function beyond that, given by
Φ ( M ) d M = Φ 1 * M M * α 1 + Φ 2 * M M * α 2 exp M M * d M M * .
Schechter parameters ( M * , α 1 , Φ 1 * , α 2 , Φ 2 * ) of the COSMOS2015 galaxy stellar mass function (SMF) are found in Table 1 of Ref. [28].
– The stellar velocity dispersion  σ is found by the relation proposed by Zahid et al. [29]. Analyzing data from the SDSS and SHELS, they established the relation between the central stellar velocity dispersion σ and the stellar mass of galaxies, given by
σ ( M ) = σ b M M b α 1 f o r M M b σ ( M ) = σ b M M b α 2 f o r M > M b .
The fit parameters are l o g ( M b / M ) = 10.26 , l o g ( σ b ) = 2.073 , α 1 = 0.403 and α 2 = 0.293 . The two indices α 1 and α 2 define the power law below and above the break point, respectively.
In the next section, we present our results considering observations with the Roman Space Telescope. In our calculations, we use the following cosmological parameters: Ω m = 0.30 , Ω k = 0 , Ω Λ = 0.70 and H 0 = 70 km s 1 Mpc 1 .

4.2. Results

We generate a sample of 101,0006 events adopting the redshift distributions of both quasars and galaxies, the galactic mass distribution and stellar velocity dispersion presented in Section 4.1. Assuming the SIS model for the galaxies, we find the probability that a quasar is lensed by one or more foreground galaxies. Considering the Roman Space Telescope, with an angular resolution of 0.11 a r c s e c , we find that about one quasar in 180 in the field of view of Roman is expected to be lensed by foreground galaxies. In the fifth column of Table 1, we give the number of lensed quasars expected to be observed by the Roman Space Telescope for each imaging sensitivity.
Furthermore, our analysis reveals that approximately 10 % of configurations can be attributed to the presence of two or more foreground galaxies, and that the probability of a quasar being lensed by more than two galaxies is very low as anticipated.
In Figure 1, the redshift distribution of galaxies and quasars in the case of the single lensing systems is plotted. In these events, two images are formed at the two sides of the quasar. As one can note, the redshift of the galaxy is up to two, while the redshift of the quasars goes up to five. In Figure 2, we present the redshift distributions in the case of a quasar lensed by two galaxies, that is, in the multiple lensing systems that are expected to be observed by the Roman Space Telescope. As can be seen, the redshift of the first galaxy is up to 1.5, and the redshift of the second galaxy is up to 4, while the quasar redshift goes up to 5.
The multiple lensing systems can exhibit up to six distinct images. Our findings regarding the distribution of the number of images are visually depicted in Figure 3.
The formation of doubly lensed quasars can occur in lensing systems comprising one or two lenses. Figure 4 displays the distribution of time delays between two images as a function of the Einstein radius. It is evident from Equations (13) and (19) that these two quantities are strongly correlated. The Einstein radii in our study range from approximately 0.11 arcsec (which is the angular accuracy of the Roman Space Telescope) to over 5 arcsec, while the time delays span from around 0.2 h to 2900 days, with a median value of approximately 10 days.

5. Conclusions

In this paper, we estimate the number of quasars that would be observed by the Roman Space Telescope. Our calculation is based on the quasar luminosity function (QLF) obtained from the Spitzer Space Telescope imaging survey [16]. We remark that we used the same QLF for both instruments, even if there is some discrepancy in their band observation. Roman F213 filter will observe (1.95–2.30 µm) (with midpoint in 2.125 µm), whereas the QLF is established for the waveband of 5 µm. We emphasize that the predicted quasar number depends on the sensitivity of the filter. For the deep survey, the number is increased about 20 times (see Table 1).
Moreover, we investigated the case of quasars lensed by one or two foreground galaxies in the case of the Roman Space Telescope. Using the redshift distributions of both quasars and galaxies, together with the galactic mass distribution and velocity dispersion of galaxies, we find that one quasar among 180 observed by Roman Space Telescope would be lensed by foreground galaxies.
Multiple lensing events are a rather rare phenomenon in astronomy, but it is very likely that the first such events will be discovered soon by the next generation of space telescopes. We find that about 10 % of the lensed quasars observable by the Roman Space Telescope are expected to be due to multiple lensing systems. Lensing by more than one galaxy along the line of sight can lead to interesting image configurations. Such systems are expected to be very important, especially in the future, both for constraining lens models of individual systems and for statistical analysis in gravitational lensing.
As discussed in the previous section, if the quasar is time variable, the time delay between the lensed images can be measured. The theoretical expression of the time delay between two lensed images is given by Equation (13). As one can see, it is the product of the time-delay distance D Δ t c , which depends on cosmology (is inversely proportional to H 0 ), and a term that depends on the lens configuration and its mass distribution. Among these quantities, only the image positions can be defined directly by observations, whereas the lens potential Ψ ( θ ) and the source position β must be inferred via lens modeling. Although the time delay can provide valuable constraints on the cosmological parameters, there are several limitations and challenges that one would need to take into account: (1) modeling uncertainties (small errors in the modeling of lensing system can lead to significant biases in the inferred cosmological parameters); (2) substructure and environment effects (the presence of substructure within the lensing object and the presence of other nearby galaxies can affect the lensing potential introducing additional time delays or perturbations in the images); (3) systematic errors (sources of systematic errors in strong gravitational lensing include telescope calibration, analysis pipeline, parametric model choice, and assumptions about the lensing object and cosmological model); (4) small sample size (few observed strong gravitational lensing systems with measured time delays limit the statistical power and robustness of the inferred cosmological parameters); (5) degeneracies with other parameters (the Hubble constant measurement from time delay cosmography is affected by degeneracies with respect to other cosmological parameters, which can lead to a broader range of allowed values for H 0 ) [31].
The Roman Space Telescope observations may allow addressing these limitations with its large field of view and high sensitivity, enabling the detection of more strong lensing events and increasing the sample size for improved statistical power. It has superior angular resolution and image quality, allowing for more accurate measurements of the lensing geometry and better identification of substructures. Additionally, operating at near-infrared wavelengths reduces uncertainties caused by dust and other obscuring materials.
We also mention that in the present work, we focused only on the SIS (singular isothermal sphere) model for lensing galaxies. However, we are currently exploring other models and how they impact the estimated time delay between lensed images. The investigation of these additional models and their influence will be addressed in a separate, more comprehensive paper.

Author Contributions

Conceptualization, L.H; M.H. and F.D.P.; methodology, L.H; M.H. and F.D.P.; software, L.H. and E.G; validation, L.H., M.H. and F.D.P.; writing—original draft preparation, L.H.; writing—review and editing, F.D.P., M.H. and E.G.; supervision, M.H. and F.D.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The APC was funded by the Faculty of Natural Science, University of Tirana, Albania.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the results presented in this paper are derived from rigorous numerical calculations based on the equations provided in the paper. No new data was created.

Acknowledgments

This article is based upon work from COST Action CA21136—Addressing observational tensions in cosmology with systematics and fundamental physics (CosmoVerse), supported by COST (European Cooperation in Science and Technology). FDP acknowledge the TAsP and Euclid projects of INFN.

Conflicts of Interest

The authors declare no conflict of interest.

Notes

1
See to https://research.ast.cam.ac.uk/lensedquasars/ (accessed on 21 May 2023).
2
In this work, the lens systems consisting of two lenses are called the multiple lensing systems.
3
4
In the case of a point-mass lens, by solving the Einstein field equations, the deflection angle is α ^ = 4 G M c 2 r m i n , where r m i n is the closest approach distance
5
This definition is only in multiple systems, whereas for its definition in single systems see Equation (4).
6
It is the sum of many smaller event samples generated separately due to the rather long time necessary for the Monte Carlo simulation to run.

References

  1. Schechter, P. The Hubble constant from gravitational lens time delays. Proc. Int. Astron. Union 2004, 2004, 281–296. [Google Scholar] [CrossRef]
  2. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.; Barreiro, R.; Bartolo, N.; Basak, S.; et al. Planck 2018 results-VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6. [Google Scholar]
  3. Abbott, T.; Abdalla, F.; Annis, J.; Bechtol, K.; Blazek, J.; Benson, B.; Bernstein, R.; Bernstein, G.; Bertin, E.; Brooks, D.; et al. Dark energy survey year 1 results: A precise H0 estimate from DES Y1, BAO, and D/H data. Mon. Not. R. Astron. Soc. 2018, 480, 3879–3888. [Google Scholar] [CrossRef]
  4. Riess, A.; Casertano, S.; Yuan, W.; Macri, L.; Scolnic, D. Large Magellanic Cloud Cepheid standards provide a 1% foundation for the determination of the Hubble constant and stronger evidence for physics beyond ΛCDM. Astrophys. J. 2019, 876, 85. [Google Scholar] [CrossRef]
  5. Wong, K.; Suyu, S.; Chen, G.; Rusu, C.; Millon, M.; Sluse, D.; Bonvin, V.; Fassnacht, C.; Taubenberger, S.; Auger, M.; et al. H0LiCOW–XIII. A 2.4 per cent measurement of H 0 from lensed quasars: 5.3 σ tension between early-and late-Universe probes. Mon. Not. R. Astron. Soc. 2020, 498, 1420–1439. [Google Scholar] [CrossRef]
  6. Myers, S.; Jackson, N.; Browne, I.; De Bruyn, A.; Pearson, T.; Readhead, A.; Wilkinson, P.; Biggs, A.; Blandford, R.; Fassnacht, C.; et al. the cosmic lens all-sky survey-I. source selection and observations. Mon. Not. R. Astron. Soc. 2003, 341, 1–12. [Google Scholar] [CrossRef]
  7. York, D.; Adelman, J.; Anderson, J., Jr.; Anderson, S.; Annis, J.; Bahcall, N.; Bakken, J.; Barkhouser, R.; Bastian, S.; Berman, E.; et al. The sloan digital sky survey: Technical summary. Astron. J. 2000, 120, 1579. [Google Scholar] [CrossRef]
  8. Keeton, C.; Falco, E.; Impey, C.; Kochanek, C.; Lehár, J.; McLeod, B.; Rix, H.; Munoz, J.; Peng, C. The host galaxy of the lensed quasar Q0957+561. Astrophys. J. 2000, 542, 74. [Google Scholar] [CrossRef]
  9. Walsh, D.; Carswell, R.; Weymann, R. 0957+ 561 A, B: Twin quasistellar objects or gravitational lens? Nature 1979, 279, 381–384. [Google Scholar] [CrossRef]
  10. Laureijs, R.; Amiaux, J.; Arduini, S.; Augueres, J.; Brinchmann, J.; Cole, R.; Cropper, M.; Dabin, C.; Duvet, L.; Ealet, A.; et al. Euclid definition study report. arXiv 2011, arXiv:1110.3193. [Google Scholar]
  11. Beaton, R.; Al-Kowsi, H.; Bellini, A.; Casertano, S.; Christian, C.; Cosentino, R.; De Rosa, G.; Desjardins, T.; Girard, J.; Gomez, S.; et al. The Roman Space Telescope Science Operations Center: Overview and Progress. Am. Astron. Soc. Meet. Abstr. 2023, 55, 207.01. [Google Scholar]
  12. Spergel, D.; Gehrels, N.; Baltay, C.; Bennett, D.; Breckinridge, J.; Donahue, M.; Dressler, A.; Gaudi, B.; Greene, T.; Guyon, O.; et al. Wide-field infrarred survey telescope-astrophysics focused telescope assets WFIRST-AFTA 2015 report. arXiv 2015, arXiv:1503.03757. [Google Scholar]
  13. Kochanek, C.; Apostolakis, J. The two-screen gravitational lens. Mon. Not. R. Astron. Soc. 1988, 235, 1073–1109. [Google Scholar] [CrossRef]
  14. Möller, O.; Blain, A. Strong gravitational lensing by multiple galaxies. Mon. Not. R. Astron. Soc. 2001, 327, 339–349. [Google Scholar] [CrossRef]
  15. Werner, M.; An, J.; Evans, N. On multiple Einstein rings. Mon. Not. R. Astron. Soc. 2008, 391, 668–674. [Google Scholar] [CrossRef]
  16. Lacy, M.; Ridgway, S.; Sajina, A.; Petric, A.; Gates, E.; Urrutia, T.; Storrie-Lombardi, L. The Spitzer mid-infrared AGN Survey. II. The demographics and cosmic evolution of the AGN population. Astrophys. J. 2015, 802, 102. [Google Scholar] [CrossRef]
  17. Hopkins, P.; Richards, G.; Hernquist, L. An observational determination of the bolometric quasar luminosity function. Astrophys. J. 2007, 654, 731. [Google Scholar] [CrossRef]
  18. Liao, K. Measuring the distances to quasars at high redshifts with strong lensing. Astrophys. J. 2019, 883, 3. [Google Scholar] [CrossRef]
  19. Schneider, P.; Ehlers, J.; Falco, E. Gravitational Lenses; Astronomy and Astrophysics Library; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1999; ISSN 0941-7834. [Google Scholar]
  20. Shapiro, I. Fourth test of general relativity. Phys. Rev. Lett. 1964, 13, 789. [Google Scholar] [CrossRef]
  21. Ofek, E.O.; Rix, H.W.; Maoz, D. The redshift distribution of gravitational lenses revisited: Constraints on galaxy mass evolution. Mon. Not. R. Astron. Soc. 2003, 343, 639–652. [Google Scholar] [CrossRef]
  22. Bolton, A.S.; Treu, T.; Koopmans, L.V.E.; Gavazzi, R.; Moustakas, L.A.; Burles, S.; Schlegel, D.J.; Wayth, R. The Sloan Lens ACS Survey VII. Elliptical galaxy scaling laws from direct observational mass measurements. Astrophys. J. 2008, 684, 248. [Google Scholar]
  23. Collett, T.; Bacon, D. Compound lensing: Einstein zig-zags and high-multiplicity lensed images. Mon. Not. R. Astron. Soc. 2016, 456, 2210–2220. [Google Scholar] [CrossRef]
  24. Chae, K.; Mao, S.; Augusto, P. Modelling the first probable two-plane lens system B2114+022: Reproducing two compact radio cores A and D. Mon. Not. R. Astron. Soc. 2001, 326, 1015–1026. [Google Scholar] [CrossRef]
  25. Hamolli, L.; Hafizi, M.; De Paolis, F.; Nucita, A. Exploiting the IRT-THESEUS capability to observe lensed Quasars. Galaxies 2021, 9, 35. [Google Scholar] [CrossRef]
  26. Schneider, D.; Hall, P.; Richards, G.; Berk, D.; Anderson, S.; Fan, X.; Jester, S.; Stoughton, C.; Strauss, M.; SubbaRao, M.; et al. The sloan digital sky survey quasar catalog. III. Third data release. Astron. J. 2005, 130, 367. [Google Scholar] [CrossRef]
  27. Appenzeller, I.; Bender, R.; Böhm, A.; Frank, S.; Fricke, K.; Gabasch, A.; Heidt, J.; Hopp, U.; Jäger, K.; Mehlert, D.; et al. Exploring Cosmic Evolution with the FORS Deep Field. Messenger 2004, 116, 18. [Google Scholar]
  28. Davidzon, I.; Ilbert, O.; Laigle, C.; Coupon, J.; McCracken, H.; Delvecchio, I.; Masters, D.; Capak, P.; Hsieh, B.; Le Fèvre, O.; et al. The COSMOS2015 galaxy stellar mass function-Thirteen billion years of stellar mass assembly in ten snapshots. Astron. Astrophys. 2017, 605, A70. [Google Scholar] [CrossRef]
  29. Zahid, H.; Geller, M.; Fabricant, D.; Hwang, H. The scaling of stellar mass and central stellar velocity dispersion for quiescent galaxies at z < 0.7. Astrophys. J. 2016, 832, 203. [Google Scholar]
  30. Beckwith, S.; Stiavelli, M.; Koekemoer, A.; Caldwell, J.; Ferguson, H.; Hook, R.; Lucas, R.; Bergeron, L.; Corbin, M.; Jogee, S.; et al. The Hubble ultra deep field. Astron. J. 2006, 132, 1729. [Google Scholar] [CrossRef]
  31. Suyu, S.H.; Treu, T.; Hilbert, S.; Sonnenfeld, A.; Auger, M.W.; Blandford, R.D.; Collett, T.; Courbin, F.; Fassnacht, C.D.; Koopmans, L.V.E.; et al. Cosmology from gravitational lens time delays and Planck data. Astrophys. J. Lett. 2014, 788, L35. [Google Scholar] [CrossRef]
Figure 1. Redshift distributions in single lensing events expected to be observed by Roman Space Telescope. The red line shows the redshift distribution of the galaxy and the gray shadow shows the redshift distribution of quasars.
Figure 1. Redshift distributions in single lensing events expected to be observed by Roman Space Telescope. The red line shows the redshift distribution of the galaxy and the gray shadow shows the redshift distribution of quasars.
Galaxies 11 00071 g001
Figure 2. Redshift distributions in multiple lensing events expected to be observed by the Roman Space Telescope. The red line shows the redshift distribution of the first galaxy, the dashed black line shows the redshift distribution of the second galaxy and gray shadow shows the redshift distribution of quasars.
Figure 2. Redshift distributions in multiple lensing events expected to be observed by the Roman Space Telescope. The red line shows the redshift distribution of the first galaxy, the dashed black line shows the redshift distribution of the second galaxy and gray shadow shows the redshift distribution of quasars.
Galaxies 11 00071 g002
Figure 3. Distribution of number of images for strong lensing events caused by multiple lenses. As one can see, up to six images can form (see text for details).
Figure 3. Distribution of number of images for strong lensing events caused by multiple lenses. As one can see, up to six images can form (see text for details).
Galaxies 11 00071 g003
Figure 4. Time delay distribution between two lensed images by the Einstein radius.
Figure 4. Time delay distribution between two lensed images by the Einstein radius.
Galaxies 11 00071 g004
Table 1. Expected number of quasars (column 3) observed by Roman Space Telescope in filter F213 for simple expose and deep survey, both for point-like sources and compact galaxies with half-light radius of 0.3″. In the fourth column is the redshift limit for each imaging sensitivity, while in the fifth column, the number of the lensed quasars are given.
Table 1. Expected number of quasars (column 3) observed by Roman Space Telescope in filter F213 for simple expose and deep survey, both for point-like sources and compact galaxies with half-light radius of 0.3″. In the fourth column is the redshift limit for each imaging sensitivity, while in the fifth column, the number of the lensed quasars are given.
F213 m AB Number of Quasars ( × 10 6 ) z max Number of Lensed Quasars
1 h, Point26.218.433.6102,000
1 h, r 50 = 0.3 25.26.322.135,000
55 s, Point23.71.0211.95660
55 s, r 50 = 0.3 22.70.2457.91360
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hamolli, L.; Hafizi, M.; De Paolis, F.; Guliqani, E. Investigating Gravitationally Lensed Quasars Observable by Nancy Grace Roman Space Telescope. Galaxies 2023, 11, 71. https://doi.org/10.3390/galaxies11030071

AMA Style

Hamolli L, Hafizi M, De Paolis F, Guliqani E. Investigating Gravitationally Lensed Quasars Observable by Nancy Grace Roman Space Telescope. Galaxies. 2023; 11(3):71. https://doi.org/10.3390/galaxies11030071

Chicago/Turabian Style

Hamolli, Lindita, Mimoza Hafizi, Francesco De Paolis, and Esmeralda Guliqani. 2023. "Investigating Gravitationally Lensed Quasars Observable by Nancy Grace Roman Space Telescope" Galaxies 11, no. 3: 71. https://doi.org/10.3390/galaxies11030071

APA Style

Hamolli, L., Hafizi, M., De Paolis, F., & Guliqani, E. (2023). Investigating Gravitationally Lensed Quasars Observable by Nancy Grace Roman Space Telescope. Galaxies, 11(3), 71. https://doi.org/10.3390/galaxies11030071

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop