Farnik Nikakhtar
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Farnik Nikakhtar has a total of 464 co-authors. Here are the top co-authors by number of shared publications:
- Ravi K. Sheth
University of Pennsylvania
11 shared publications - Nikhil Padmanabhan
Yale University
7 shared publications - Robyn E. Sanderson
University of Pennsylvania
5 shared publications - Bruno Lévy
Centre Inria de Saclay
4 shared publications Ravi K. Sheth
Center for Particle Cosmology University of Pennsylvania 19104 Philadelphia PA USA
4 shared publicationsRoya Mohayaee
Sorbonne Université CNRS Institut d'Astrophysique de Paris 98bis Bld Arago 75014 Paris France
4 shared publications
Latest Publications
- 2 October 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 114 (8) , 083502
Optimal transport reconstruction of biased tracers in primordial non-Gaussian fields
Optimal transport provides an efficient method to infer the displacement of objects by mapping their initial positions to their present-day locations over cosmic time; equivalently, it enables the reconstruction of initial positions from measurements taken at later times. The method has been shown to be accurate even if positions for only a biased subset of the particles are measured, provided that the initial displacement field was Gaussian. The method does not rely on the assumption of a Gaussian displacement field and thus may be extended to the reconstruction of non-Gaussian initial conditions. Here, we demonstrate how this is achieved for a class of ``local'' primordial non-Gaussian fields of current interest in cosmology. For these models, there is a distinctive signature in the large scale clustering of biased tracers that depends on the product of the primordial amplitude ${f}_{\mathrm{NL}}$ and the nature of the tracers ${b}_{\ensuremath{\phi}}$. Our method exploits the fact that this signature is not present in the full field; it is only present in biased fields. Therefore, the mass that is not in the biased subset, what we call the ``dust,'' also has a characteristic scale dependence, albeit of a different amplitude. We show that the quality of the optimal transport reconstruction improves as the model for this dust becomes more realistic. Importantly, the dust model requires only the combination ${f}_{\mathrm{NL}}{b}_{\ensuremath{\phi}}$, and not ${f}_{\mathrm{NL}}$ itself: this is the same quantity that the scale-dependent bias signal measures directly, so the reconstruction is conditioned on an observable rather than on prior knowledge of the primordial amplitude.- 2 October 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 114 (8) , 083502
Optimal transport reconstruction of biased tracers in primordial non-Gaussian fields
Optimal transport provides an efficient method to infer the displacement of objects by mapping their initial positions to their present-day locations over cosmic time; equivalently, it enables the reconstruction of initial positions from measurements taken at later times. The method has been shown to be accurate even if positions for only a biased subset of the particles are measured, provided that the initial displacement field was Gaussian. The method does not rely on the assumption of a Gaussian displacement field and thus may be extended to the reconstruction of non-Gaussian initial conditions. Here, we demonstrate how this is achieved for a class of ``local'' primordial non-Gaussian fields of current interest in cosmology. For these models, there is a distinctive signature in the large scale clustering of biased tracers that depends on the product of the primordial amplitude ${f}_{\mathrm{NL}}$ and the nature of the tracers ${b}_{\ensuremath{\phi}}$. Our method exploits the fact that this signature is not present in the full field; it is only present in biased fields. Therefore, the mass that is not in the biased subset, what we call the ``dust,'' also has a characteristic scale dependence, albeit of a different amplitude. We show that the quality of the optimal transport reconstruction improves as the model for this dust becomes more realistic. Importantly, the dust model requires only the combination ${f}_{\mathrm{NL}}{b}_{\ensuremath{\phi}}$, and not ${f}_{\mathrm{NL}}$ itself: this is the same quantity that the scale-dependent bias signal measures directly, so the reconstruction is conditioned on an observable rather than on prior knowledge of the primordial amplitude.- Open Access
- 1 October 2026
- Published by IOP Publishing in Journal of Cosmology and Astroparticle Physics
- Vol. 2026 (10)
The linear point standard ruler with DESI DR1 and DR2 data
The linear point s LP , a purely geometric feature in the monopole of the two-point correlation function, has been proposed as an alternative standard ruler that, unlike the peak in the correlation function, is more robust to late-time non-linear effects at the percent level. In light of improved simulations and high-quality data, we revisit this claimed robustness and test whether it persists at the sub-percent precision demanded by current surveys, using the linear point as an alternative to the template-based fitting approaches employed in standard BAO analyses. We present the first linear point measurements on galaxy samples from the first and second data releases (DR1 and DR2) of the Dark Energy Spectroscopic Instrument (DESI). We convert the linear point to a dimensionless parameter α iso,LP , defined as the ratio of the fiducial linear point to the observed value, analogous to the isotropic BAO scaling parameter α iso,BAO obtained from template-based BAO fits. Using the second generation of AbacusSummit mock catalogs for DESI DR1 ( Abacus - 2 ), we find that post-reconstruction linear point measurements are more precise than pre-reconstruction, with 15–60% smaller uncertainties. We compare the mean α iso,LP across the 25 Abacus - 2 mocks for each tracer to the mean α iso,BAO from the DESI template-based fits, and find a systematic offset which we attribute to the smearing of the correlation function, and thereby the linear point, in the non-linear regime. This shift is small relative to the method's original percent-level robustness, but becomes significant at the sub-percent level. We propose a sample-dependent correction to α iso,LP , as a function of the isotropic damping parameter Σ iso , that mitigates this offset; the resulting α iso,LP corrected values show excellent agreement with α iso,BAO in both mocks and data, particularly post-reconstruction. This agreement comes at the price of the model-independence that motivated the linear point in the first place, since both reconstruction and the proposed correction are cosmology-dependent. We conclude that, at present precision, the linear point is best regarded as a complementary geometric cross-check to BAO analyses rather than a strict alternative to template-based fitting approaches — one that grows more valuable as DESI hints at time-evolving dark energy, where independent geometric probes of the expansion history help test its robustness.
- 1 October 2026
- Published by IOP Publishing in Journal of Cosmology and Astroparticle Physics
- Vol. 2026 (10)
Open AccessThe linear point standard ruler with DESI DR1 and DR2 data
The linear point s LP , a purely geometric feature in the monopole of the two-point correlation function, has been proposed as an alternative standard ruler that, unlike the peak in the correlation function, is more robust to late-time non-linear effects at the percent level. In light of improved simulations and high-quality data, we revisit this claimed robustness and test whether it persists at the sub-percent precision demanded by current surveys, using the linear point as an alternative to the template-based fitting approaches employed in standard BAO analyses. We present the first linear point measurements on galaxy samples from the first and second data releases (DR1 and DR2) of the Dark Energy Spectroscopic Instrument (DESI). We convert the linear point to a dimensionless parameter α iso,LP , defined as the ratio of the fiducial linear point to the observed value, analogous to the isotropic BAO scaling parameter α iso,BAO obtained from template-based BAO fits. Using the second generation of AbacusSummit mock catalogs for DESI DR1 ( Abacus - 2 ), we find that post-reconstruction linear point measurements are more precise than pre-reconstruction, with 15–60% smaller uncertainties. We compare the mean α iso,LP across the 25 Abacus - 2 mocks for each tracer to the mean α iso,BAO from the DESI template-based fits, and find a systematic offset which we attribute to the smearing of the correlation function, and thereby the linear point, in the non-linear regime. This shift is small relative to the method's original percent-level robustness, but becomes significant at the sub-percent level. We propose a sample-dependent correction to α iso,LP , as a function of the isotropic damping parameter Σ iso , that mitigates this offset; the resulting α iso,LP corrected values show excellent agreement with α iso,BAO in both mocks and data, particularly post-reconstruction. This agreement comes at the price of the model-independence that motivated the linear point in the first place, since both reconstruction and the proposed correction are cosmology-dependent. We conclude that, at present precision, the linear point is best regarded as a complementary geometric cross-check to BAO analyses rather than a strict alternative to template-based fitting approaches — one that grows more valuable as DESI hints at time-evolving dark energy, where independent geometric probes of the expansion history help test its robustness.
- Open Access
- 5 February 2026
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 998 (1) , 130
Reducing the Dimensions of Active Galactic Nuclei Light-curve Manifolds
The active galactic nuclei (AGNs) glossary is vast and complex. Depending on the selection method, observing wavelength, and brightness, AGNs are assigned distinct labels, yet the relationship between different selection methods and the diversity of time-domain behavior within and across classes remains difficult to characterize in a unified framework. Changing-look AGNs (CLAGNs), which transition between classifications over time, further complicate this picture. In this work, we learn a data-driven, low-dimensional representation of multiwavelength photometric light curves of AGNs, in which the structure of the projected manifold correlates with AGN class and independent spectroscopic properties. Using the NASA Fornax Science Platform, we assemble light curves from Zwicky Transient Facility, Pan-STARRS, Gaia, and Wide-field Infrared Survey Explorer/NEOWISE for two samples: (1) a heterogeneous set of ∼2000 AGNs spanning z ≲ 1, including Sloan Digital Sky Survey quasars, variability-selected sources, and CLAGNs; and (2) a homogeneous sample of ∼65,000 narrow-line AGNs at z ≈ 0.1 with well-characterized optical emission-line measurements. Without using class labels during training, the learned manifolds organize variability-selected AGNs into coherent regions of the low-dimensional space, distinguish between turn-on and turn-off CLAGNs, and place tidal disruption events in distinct regions. Manifold coordinates correlate with key spectroscopic and host-galaxy properties—including stellar mass, [O III ] luminosity, and D n (4000)—demonstrating that heterogeneous multiband variability can be combined in a purely data-driven manner to recover correlations with independent physical diagnostics, without requiring explicit physical modeling. These results show that manifold learning offers a practical, assumption-light approach for integrating time-domain surveys and prioritizing spectroscopic follow-up.
- 5 February 2026
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 998 (1) , 130
Open AccessReducing the Dimensions of Active Galactic Nuclei Light-curve Manifolds
The active galactic nuclei (AGNs) glossary is vast and complex. Depending on the selection method, observing wavelength, and brightness, AGNs are assigned distinct labels, yet the relationship between different selection methods and the diversity of time-domain behavior within and across classes remains difficult to characterize in a unified framework. Changing-look AGNs (CLAGNs), which transition between classifications over time, further complicate this picture. In this work, we learn a data-driven, low-dimensional representation of multiwavelength photometric light curves of AGNs, in which the structure of the projected manifold correlates with AGN class and independent spectroscopic properties. Using the NASA Fornax Science Platform, we assemble light curves from Zwicky Transient Facility, Pan-STARRS, Gaia, and Wide-field Infrared Survey Explorer/NEOWISE for two samples: (1) a heterogeneous set of ∼2000 AGNs spanning z ≲ 1, including Sloan Digital Sky Survey quasars, variability-selected sources, and CLAGNs; and (2) a homogeneous sample of ∼65,000 narrow-line AGNs at z ≈ 0.1 with well-characterized optical emission-line measurements. Without using class labels during training, the learned manifolds organize variability-selected AGNs into coherent regions of the low-dimensional space, distinguish between turn-on and turn-off CLAGNs, and place tidal disruption events in distinct regions. Manifold coordinates correlate with key spectroscopic and host-galaxy properties—including stellar mass, [O III ] luminosity, and D n (4000)—demonstrating that heterogeneous multiband variability can be combined in a purely data-driven manner to recover correlations with independent physical diagnostics, without requiring explicit physical modeling. These results show that manifold learning offers a practical, assumption-light approach for integrating time-domain surveys and prioritizing spectroscopic follow-up.
- Open Access
The Linear Point Standard Ruler with DESI DR1 and DR2 Data
The linear point, a purely geometric feature in the monopole of the two-point correlation function, has been proposed as an alternative standard ruler. Compared to the peak in the correlation function, it is more robust to late-time nonlinear effects at the percent level. In light of improved simulations and high quality data, we revisit the robustness of the linear point and use it as an alternative to template-based fitting approaches typically used in BAO analyses. We present the linear point measurements on galaxy samples from the first and second data releases (DR1 and DR2) of the DESI survey. We convert the linear point into a dimensionless parameter $α_{iso,LP}$, defined as the ratio of the linear point in the fiducial cosmology and the observed value, analogous to the isotropic BAO scaling parameter $α_{iso}$ used in previous BAO measurements. Using the 2nd generation of AbacusSummit mock catalogs, we find that linear point measurements are more precise when calculated in the post-reconstruction regime with 15-60% smaller uncertainties than those pre-reconstruction. We find a systematic shift in the linear point measurements compared against the isotropic BAO measurements in mocks; we attribute this to the isotropic damping parameter responsible for smearing the linear point in the nonlinear regime. We propose a sample-dependent correction that mitigates the impact of late-time nonlinear effects. While this introduces a cosmology dependence in an otherwise model-independent measurement, this is necessary given the sub-percent precision dictated by current cosmological surveys. Comparing $α_{iso,LP}$ with isotropic BAO measurements made on the DESI DR1 and DR2 galaxy samples, we find excellent agreement after applying this correction, particularly post-reconstruction. We discuss future scope regarding cosmological inference with linear point measurements.Open AccessThe Linear Point Standard Ruler with DESI DR1 and DR2 Data
The linear point, a purely geometric feature in the monopole of the two-point correlation function, has been proposed as an alternative standard ruler. Compared to the peak in the correlation function, it is more robust to late-time nonlinear effects at the percent level. In light of improved simulations and high quality data, we revisit the robustness of the linear point and use it as an alternative to template-based fitting approaches typically used in BAO analyses. We present the linear point measurements on galaxy samples from the first and second data releases (DR1 and DR2) of the DESI survey. We convert the linear point into a dimensionless parameter $α_{iso,LP}$, defined as the ratio of the linear point in the fiducial cosmology and the observed value, analogous to the isotropic BAO scaling parameter $α_{iso}$ used in previous BAO measurements. Using the 2nd generation of AbacusSummit mock catalogs, we find that linear point measurements are more precise when calculated in the post-reconstruction regime with 15-60% smaller uncertainties than those pre-reconstruction. We find a systematic shift in the linear point measurements compared against the isotropic BAO measurements in mocks; we attribute this to the isotropic damping parameter responsible for smearing the linear point in the nonlinear regime. We propose a sample-dependent correction that mitigates the impact of late-time nonlinear effects. While this introduces a cosmology dependence in an otherwise model-independent measurement, this is necessary given the sub-percent precision dictated by current cosmological surveys. Comparing $α_{iso,LP}$ with isotropic BAO measurements made on the DESI DR1 and DR2 galaxy samples, we find excellent agreement after applying this correction, particularly post-reconstruction. We discuss future scope regarding cosmological inference with linear point measurements. - 12 November 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 110 (10) , 103515
Eigen-decomposition of covariance matrices: An application to the BAO linear point
The baryon acoustic oscillations (BAO) feature in the two-point correlation function (TPCF) of discrete tracers such as galaxies is an accurate standard ruler. The covariance matrix of the TPCF plays an important role in determining how the precision of this ruler depends on the number density and clustering strength of the tracers, as well as the survey volume. An eigen-decomposition of this matrix provides an objective way to separate the contributions of cosmic variance from those of shot-noise to the statistical uncertainties. For the signal-to-noise levels that are expected in ongoing and next-generation surveys, the cosmic variance eigen-modes dominate. These modes are smooth functions of scale, meaning that they are insensitive to the modest changes in binning that are allowed if one wishes to resolve the BAO feature in the TPCF; they provide a good description of the correlated residuals which result from fitting smooth functional forms to the measured TPCF; they motivate a simple but accurate approximation for the uncertainty on the linear point (LP) estimate of the BAO distance scale. This approximation allows one to quantify the precision of the BAO distance scale estimate without having to generate a large ensemble of mock catalogs and explains why: the uncertainty on the LP does not depend on the functional form fitted to the TPCF or the binning used; the LP is more constraining than the peak or dip scales in the TPCF; the evolved TPCF is less constraining than the initial one, so that reconstruction schemes can yield significant gains in precision.
- 12 November 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 110 (10) , 103515
Eigen-decomposition of covariance matrices: An application to the BAO linear point
The baryon acoustic oscillations (BAO) feature in the two-point correlation function (TPCF) of discrete tracers such as galaxies is an accurate standard ruler. The covariance matrix of the TPCF plays an important role in determining how the precision of this ruler depends on the number density and clustering strength of the tracers, as well as the survey volume. An eigen-decomposition of this matrix provides an objective way to separate the contributions of cosmic variance from those of shot-noise to the statistical uncertainties. For the signal-to-noise levels that are expected in ongoing and next-generation surveys, the cosmic variance eigen-modes dominate. These modes are smooth functions of scale, meaning that they are insensitive to the modest changes in binning that are allowed if one wishes to resolve the BAO feature in the TPCF; they provide a good description of the correlated residuals which result from fitting smooth functional forms to the measured TPCF; they motivate a simple but accurate approximation for the uncertainty on the linear point (LP) estimate of the BAO distance scale. This approximation allows one to quantify the precision of the BAO distance scale estimate without having to generate a large ensemble of mock catalogs and explains why: the uncertainty on the LP does not depend on the functional form fitted to the TPCF or the binning used; the LP is more constraining than the peak or dip scales in the TPCF; the evolved TPCF is less constraining than the initial one, so that reconstruction schemes can yield significant gains in precision.
- 12 November 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 110 (10)
Leveraging protohalos and scale-dependent bias to calibrate the BAO scale in real space
- 12 November 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 110 (10)
Leveraging protohalos and scale-dependent bias to calibrate the BAO scale in real space
- Open Access
- 23 October 2024
- Published by The Open Journal in Journal of Open Source Software
- Vol. 9 (102) , 6234
Generating synthetic star catalogs from simulated data for next-gen observatories with py-ananke
We find ourselves on the brink of an exciting era in observational astrophysics, driven by groundbreaking facilities like JWST, Euclid, Rubin, Roman, SKA, or ELT. Simultaneously, computational astrophysics has shown significant strides, yielding highly realistic galaxy formation simulations, thanks to both hardware and software enhancements. Bridging the gap between simulations and observations has become paramount for meaningful comparisons.
- 23 October 2024
- Published by The Open Journal in Journal of Open Source Software
- Vol. 9 (102) , 6234
Open AccessGenerating synthetic star catalogs from simulated data for next-gen observatories with py-ananke
We find ourselves on the brink of an exciting era in observational astrophysics, driven by groundbreaking facilities like JWST, Euclid, Rubin, Roman, SKA, or ELT. Simultaneously, computational astrophysics has shown significant strides, yielding highly realistic galaxy formation simulations, thanks to both hardware and software enhancements. Bridging the gap between simulations and observations has become paramount for meaningful comparisons.
- Open Access
Leveraging protohalos and scale-dependent bias to calibrate the BAO scale in real space
The location of the baryon acoustic oscillation (BAO) feature in the two-point correlation function (2PCF) of matter produces a standard ruler that is useful for the measurement of the expansion history of the Universe. Inspired by the possibility of reconstructing the positions of protohalos in the initial density field with a novel method rooted in optimal transport theory, we revisit the BAO signal in the protohalo correlation function. Our work examines the performance of a template 2PCF built on a tracer bias relation that includes scale dependence -- a term that can be motivated by peaks theory or a general bias expansion. Working in protohalos, halos, and the linear combination of the protohalo and matter fields that is motivated by the continuity equation, we demonstrate that this model accurately captures the shape of the BAO feature and improves the precision of the BAO scale measurement relative to a model that does not include scale-dependent bias by 47% in protohalos, 15% in halos, and 14% in the linear combination of the protohalo and matter fields. Allowing for scale dependence does not appear to introduce any shift in the BAO feature. The precision of the BAO distance scale estimate is highest with the linear combination of the protohalo and matter fields, which offers a factor of 3.5 improvement over Eulerian-space measurements and a factor of 4-8 improvement over the estimate made with protohalos alone.
Open AccessLeveraging protohalos and scale-dependent bias to calibrate the BAO scale in real space
The location of the baryon acoustic oscillation (BAO) feature in the two-point correlation function (2PCF) of matter produces a standard ruler that is useful for the measurement of the expansion history of the Universe. Inspired by the possibility of reconstructing the positions of protohalos in the initial density field with a novel method rooted in optimal transport theory, we revisit the BAO signal in the protohalo correlation function. Our work examines the performance of a template 2PCF built on a tracer bias relation that includes scale dependence -- a term that can be motivated by peaks theory or a general bias expansion. Working in protohalos, halos, and the linear combination of the protohalo and matter fields that is motivated by the continuity equation, we demonstrate that this model accurately captures the shape of the BAO feature and improves the precision of the BAO scale measurement relative to a model that does not include scale-dependent bias by 47% in protohalos, 15% in halos, and 14% in the linear combination of the protohalo and matter fields. Allowing for scale dependence does not appear to introduce any shift in the BAO feature. The precision of the BAO distance scale estimate is highest with the linear combination of the protohalo and matter fields, which offers a factor of 3.5 improvement over Eulerian-space measurements and a factor of 4-8 improvement over the estimate made with protohalos alone.
- 6 June 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 109 (12)
Displacement field analysis via optimal transport: Multitracer approach to cosmological reconstruction
- 6 June 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 109 (12)
Displacement field analysis via optimal transport: Multitracer approach to cosmological reconstruction
- Open Access
Displacement Field Analysis via Optimal Transport: Multi-Tracer Approach to Cosmological Reconstruction
We demonstrate the effectiveness of one of the many multi-tracer analyses enabled by Optimal Transport (OT) reconstruction. Leveraging a semi-discrete OT algorithm, we determine the displacements between initial and observed positions of biased tracers and the remaining matter field. With only redshift-space distorted final positions of biased tracers and a simple premise for the remaining mass distribution as input, OT solves the displacement field. This extracted field, assuming asymptotically uniform density and a gradient flow displacement, enables reconstruction of the initial overdensity fluctuation field. We show that the divergence of the OT displacement field is a good proxy of the linear density field, even though the method never assumes the linear theory growth. Additionally, this divergence field can be combined with the reconstructed protohalos to provide a higher signal-to-noise measurement of the BAO standard ruler than was possible with either measurement individually.
Open AccessDisplacement Field Analysis via Optimal Transport: Multi-Tracer Approach to Cosmological Reconstruction
We demonstrate the effectiveness of one of the many multi-tracer analyses enabled by Optimal Transport (OT) reconstruction. Leveraging a semi-discrete OT algorithm, we determine the displacements between initial and observed positions of biased tracers and the remaining matter field. With only redshift-space distorted final positions of biased tracers and a simple premise for the remaining mass distribution as input, OT solves the displacement field. This extracted field, assuming asymptotically uniform density and a gradient flow displacement, enables reconstruction of the initial overdensity fluctuation field. We show that the divergence of the OT displacement field is a good proxy of the linear density field, even though the method never assumes the linear theory growth. Additionally, this divergence field can be combined with the reconstructed protohalos to provide a higher signal-to-noise measurement of the BAO standard ruler than was possible with either measurement individually.