BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

S. H. Suyu

S. H. Suyu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COastro-ph.HEastro-ph.GAastro-ph.IMApplicationsastro-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2026arXiv

Rates of Strongly Lensed Tidal Disruption Events

In the coming years, surveys such as the Rubin Observatory's Legacy Survey of Space and Time (LSST) are expected to increase the number of observed Tidal Disruption Events (TDEs) substantially. We employ Monte Carlo integration to calculate the unlensed and lensed TDE rate as a function of limiting magnitude in $u$, $g$, $r$, and $i$-bands. We investigate the impact of multiple luminosity models, black hole mass functions (BHMFs), and flare temperatures on the TDE rate. Notably, this includes a semi-analytical model, which enables the determination of the TDE temperature in terms of black hole (BH) mass. We predict the highest unlensed TDE rate to be in $g$-band. It ranges from $16$ to $5,440\;\mathrm{yr}^{-1}\;(20,000\;\mathrm{deg}^2)^{-1}$ for the Zwicky Transient Facility, being more consistent with the observed rate at the low end. For LSST, we expect a rate in $g$-band between $3,580$ and $82,060\;\mathrm{yr}^{-1}\;(20,000\;\mathrm{deg}^2)^{-1}$. A higher theoretical prediction is understandable, as we do not consider observational effects such as completeness. The unlensed and lensed TDE rates are insensitive to the redshift evolution of the BHMF, even for LSST limiting magnitudes. The best band for detecting lensed TDEs is also $g$-band. Its predicted rates range from $0.43$ to $15\;\mathrm{yr}^{-1}\;(20,000\;\mathrm{deg}^2)^{-1}$ for LSST. The scatter of predicted rates reduces when we consider the fraction of lensed TDEs; that is, a few in ten thousand TDEs will be lensed. Despite the large scatter in the rates of lensed TDEs, our comprehensive considerations of multiple models suggest that lensed TDEs will occur in the $10$-year LSST lifetime, providing an exciting prospect for detecting such events. We expect the median redshift of a lensed TDE to be between $1.5$ and $2$. In this paper, we additionally report on lensed TDE properties, such as the BH mass and time delays.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

GIGA-Lens: Fast Bayesian Inference for Strong Gravitational Lens Modeling

We present GIGA-Lens: a gradient-informed, GPU-accelerated Bayesian framework for modeling strong gravitational lensing systems, implemented in TensorFlow and JAX. The three components, optimization using multi-start gradient descent, posterior covariance estimation with variational inference, and sampling via Hamiltonian Monte Carlo, all take advantage of gradient information through automatic differentiation and massive parallelization on graphics processing units (GPUs). We test our pipeline on a large set of simulated systems and demonstrate in detail its high level of performance. The average time to model a single system on four Nvidia A100 GPUs is 105 seconds. The robustness, speed, and scalability offered by this framework make it possible to model the large number of strong lenses found in current surveys and present a very promising prospect for the modeling of $\mathcal{O}(10^5)$ lensing systems expected to be discovered in the era of the Vera C. Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Galaxy cluster strong lensing cosmography: cosmological constraints from a sample of regular galaxy clusters

Cluster strong lensing cosmography is a promising probe of the background geometry of the Universe and several studies have emerged, thanks to the increased quality of observations using space and ground-based telescopes. For the first time, we use a sample of five cluster strong lenses to measure the values of cosmological parameters and combine them with those from classical probes. In order to assess the degeneracies and the effectiveness of strong-lensing cosmography in constraining the background geometry of the Universe, we adopt four cosmological scenarios. We find good constraining power on the total matter density of the Universe ($Ω_{\rm m}$) and the equation of state of the dark energy parameter $w$. For a flat $w$CDM cosmology, we find $Ω_{\rm m} = 0.30_{-0.11}^{+0.09}$ and $w=-1.12_{-0.32}^{+0.17}$ from strong lensing only. Interestingly, we show that the constraints from the Cosmic Microwave Background (CMB) are improved by factors of 2.5 and 4.0 on $Ω_{\rm m}$ and $w$, respectively, when combined with our posterior distributions in this cosmological model. In a scenario where the equation of state of dark energy evolves with redshift, the strong lensing constraints are compatible with a cosmological constant (i.e. $w=-1$). In a curved cosmology, our strong lensing analyses can accommodate a large range of values for the curvature of the Universe of $Ω_{\rm k}=0.28_{-0.21}^{+0.16}$. In all cosmological scenarios, we show that our strong lensing constraints are complementary and in good agreement with measurements from the CMB, baryon acoustic oscillations and Type Ia supernovae. Our results show that cluster strong lensing cosmography is a potentially powerful probe to be included in the cosmological analyses of future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

HOLISMOKES -- III. Achromatic Phase of Strongly Lensed Type Ia Supernovae

To use strongly lensed Type Ia supernovae (LSNe Ia) for cosmology, a time-delay measurement between the multiple supernova (SN) images is necessary. The sharp rise and decline of SN Ia light curves make them promising for measuring time delays, but microlensing can distort these light curves and therefore add large uncertainties to the measurements. An alternative approach is to use color curves where uncertainties due to microlensing are significantly reduced for a certain period of time known as the achromatic phase. In this work, we investigate in detail the achromatic phase, testing four different SN Ia models with various microlensing configurations. We find on average an achromatic phase of around three rest-frame weeks or longer for most color curves but the spread in the duration of the achromatic phase (due to different microlensing maps and filter combinations) is quite large and an achromatic phase of just a few days is also possible. Furthermore, the achromatic phase is longer for smoother microlensing maps, lower macro-magnifications and larger mean Einstein radii of microlenses. From our investigations, we do not find a strong dependency on the model or on asymmetries in the SN ejecta. Further, we find that three independent LSST color curves exhibit features such as extreme points or turning points within the achromatic phase, which make them promising for time-delay measurements. These curves contain combinations of rest-frame bands $u$, $g$, $r$, and $i$ and to observe them for typical LSN Ia redshifts, it would be ideal to cover (observer-frame) filters $r$, $i$, $z$, $y$, $J$, and $H$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

HOLISMOKES -- VII. Time-delay measurement of strongly lensed Type Ia supernovae using machine learning

The Hubble constant ($H_0$) is one of the fundamental parameters in cosmology, but there is a heated debate around the $>$4$σ$ tension between the local Cepheid distance ladder and the early Universe measurements. Strongly lensed Type Ia supernovae (LSNe Ia) are an independent and direct way to measure $H_0$, where a time-delay measurement between the multiple supernova (SN) images is required. In this work, we present two machine learning approaches for measuring time delays in LSNe Ia, namely, a fully connected neural network (FCNN) and a random forest (RF). For the training of the FCNN and the RF, we simulate mock LSNe Ia from theoretical SN Ia models that include observational noise and microlensing. We test the generalizability of the machine learning models by using a final test set based on empirical LSN Ia light curves not used in the training process, and we find that only the RF provides a low enough bias to achieve precision cosmology; as such, RF is therefore preferred over our FCNN approach for applications to real systems. For the RF with single-band photometry in the $i$ band, we obtain an accuracy better than 1\% in all investigated cases for time delays longer than 15 days, assuming follow-up observations with a 5$σ$ point-source depth of 24.7, a two day cadence with a few random gaps, and a detection of the LSNe Ia 8 to 10 days before peak in the observer frame. In terms of precision, we can achieve an approximately 1.5-day uncertainty for a typical source redshift of $\sim$0.8 on the $i$ band under the same assumptions. To improve the measurement, we find that using three bands, where we train a RF for each band separately and combine them afterward, helps to reduce the uncertainty to $\sim$1.0 day. We have publicly released the microlensed spectra and light curves used in this work.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

H0LiCOW XI. A weak lensing measurement of the external convergence in the field of the lensed quasar B1608+656 using HST and Subaru deep imaging

We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $κ_{\rm ext}$ at the position of B1608+656 is $κ_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

HOLISMOKES -- I. Highly Optimised Lensing Investigations of Supernovae, Microlensing Objects, and Kinematics of Ellipticals and Spirals

We present the HOLISMOKES programme on strong gravitational lensing of supernovae as a probe of supernova (SN) physics and cosmology. We investigate the effects of microlensing on early-phase SN Ia spectra using four different SN explosion models, and find that within 10 rest-frame days after SN explosion, distortions of SN Ia spectra due to microlensing are typically negligible ($<$1% distortion within the 1$σ$ spread, and $\lesssim$10% distortion within the 2$σ$ spread). This shows great prospects of using lensed SNe Ia to obtain intrinsic early-phase SN spectra for deciphering SN Ia progenitors. As a demonstration of the usefulness of lensed SNe Ia for cosmology, we simulate a sample of mock lensed SN Ia systems that are expected to have accurate and precise time-delay measurements in the era of the Rubin Observatory Legacy Survey of Space and Time (LSST). Adopting realistic yet conservative uncertainties on their time-delay distances and lens angular diameter distances (of 6.6% and 5%, respectively), we find that a sample of 20 lensed SNe Ia would allow a constraint on the Hubble constant ($H_0$) with 1.3% uncertainty in the flat $Λ$CDM cosmology. We find a similar constraint on $H_0$ in an open $Λ$CDM cosmology, while the constraint degrades to $3\%$ in a flat $w$CDM cosmology. We anticipate lensed SNe to be an independent and powerful probe of SN physics and cosmology in the upcoming LSST era.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

HOLISMOKES -- IV. Efficient mass modeling of strong lenses through deep learning

Modelling the mass distributions of strong gravitational lenses is often necessary to use them as astrophysical and cosmological probes. With the high number of lens systems ($>10^5$) expected from upcoming surveys, it is timely to explore efficient modeling approaches beyond traditional MCMC techniques that are time consuming. We train a CNN on images of galaxy-scale lenses to predict the parameters of the SIE mass model ($x,y,e_x,e_y$, and $θ_E$). To train the network, we simulate images based on real observations from the HSC Survey for the lens galaxies and from the HUDF as lensed galaxies. We tested different network architectures, the effect of different data sets, and using different input distributions of $θ_E$. We find that the CNN performs well and obtain with the network trained with a uniform distribution of $θ_E$ $>0.5&#34;$ the following median values with $1σ$ scatter: $Δx=(0.00^{+0.30}_{-0.30})&#34;$, $Δy=(0.00^{+0.30}_{-0.29})&#34; $, $Δθ_E=(0.07^{+0.29}_{-0.12})&#34;$, $Δe_x = -0.01^{+0.08}_{-0.09}$ and $Δe_y = 0.00^{+0.08}_{-0.09}$. The bias in $θ_E$ is driven by systems with small $θ_E$. Therefore, when we further predict the multiple lensed image positions and time delays based on the network output, we apply the network to the sample limited to $θ_E>0.8&#34;$. In this case, the offset between the predicted and input lensed image positions is $(0.00_{-0.29}^{+0.29})&#34;$ and $(0.00_{-0.31}^{+0.32})&#34;$ for $x$ and $y$, respectively. For the fractional difference between the predicted and true time delay, we obtain $0.04_{-0.05}^{+0.27}$. Our CNN is able to predict the SIE parameters in fractions of a second on a single CPU and with the output we can predict the image positions and time delays in an automated way, such that we are able to process efficiently the huge amount of expected lens detections in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Is every strong lens model unhappy in its own way? Uniform modelling of a sample of 13 quadruply+ imaged quasars

Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the accuracy of time-delay cosmography in the Frontier Fields Cluster MACS J1149.5+2223 with supernova Refsdal

We study possible systematic effects on the values of the cosmological parameters measured through strong lensing analyses of the HFF galaxy cluster MACS J1149.5+2223. We use the observed positions of a large set of spectroscopically selected multiple images, including those of supernova &#34;Refsdal&#34; with their estimated time delays. Starting from our reference model in a flat $Λ$CDM cosmology, published in Grillo et al. (2018), we confirm the relevance of the longest measurable time delay, between SX and S1, and an approximately linear relation between its value and that of $H_{0}$. We perform true blind tests by considering a range of time delays around its original estimate of $345 \pm 10$ days, as an accurate measurement of this time delay was not known at the time of analysis and writing. We investigate separately the impact of a constant sheet of mass at the cluster redshift, of a power-law profile for the mass density of the cluster main halo and of some scatter in the cluster member scaling relations. Remarkably, we find that these systematic effects do not introduce a significant bias on the inferred values of $H_{0}$ and $Ω_{\rm m}$ and that the statistical uncertainties dominate the total error budget: a 3% uncertainty on the time delay of image SX translates into approximately 6% and 40% (including both statistical and systematic $1σ$) uncertainties for $H_{0}$ and $Ω_{\rm m}$, respectively. Furthermore, our model accurately reproduces the extended surface brightness distribution of the supernova host, covering more than $3 \times 10^{4}$ $HST$ pixels. We also present the interesting possibility of measuring the value of the equation-of-state parameter $w$ of the dark energy density, currently with a 30% uncertainty. We conclude that time-delay cluster lenses have the potential to become soon an alternative and competitive cosmological probe.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

TDCOSMO. I. An exploration of systematic uncertainties in the inference of $H_0$ from time-delay cosmography

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the &#34;bulge-halo&#34; conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dissecting the Gravitational Lens B1608+656. I. Lens Potential Reconstruction

Strong gravitational lensing is a powerful technique for probing galaxy mass distributions and for measuring cosmological parameters. We present a pixelated approach to modeling simultaneously the lens potential and source intensity of strong gravitational lens systems with extended source-intensity distributions. For systems with sources of sufficient extent such that the separate lensed images are connected by intensity measurements, the accuracy in the reconstructed potential is solely limited by the quality of the data. We apply this potential reconstruction technique to deep HST observations of B1608+656, a four-image gravitational lens system formed by a pair of interacting lens galaxies. We present a comprehensive Bayesian analysis of the system that takes into account the extended source-intensity distribution, dust extinction, and the interacting lens galaxies. Our approach allows us to compare various models of the components of the lens system, which include the point-spread function (PSF), dust, lens galaxy light, source-intensity distribution, and lens potential. Using optimal combinations of the PSF, dust, and lens galaxy light models, we successfully reconstruct both the lens potential and the extended source-intensity distribution of B1608+656. The resulting reconstruction can be used as the basis of a measurement of the Hubble constant. We use our reconstruction of the gravitational potential to study the relative distribution of mass and light in the lensing galaxies. We find that the mass-to-light ratio for the primary lens galaxy is (2.0+/-0.2)h M_{\sun} L_{B,\sun}^{-1} within the Einstein radius 3.9 h^{-1} kpc, in agreement with what is found for noninteracting lens galaxies at the same scales. (Abridged)

0 citations0 reviewsNo code linkedOpen access