I'm Boryana Hadzhiyska

a PhD Candidate


at the Center for Astrophysics | Harvard & Smithsonian

Contact Me

About Me

About Me

I'm Boryana Hadzhiyska, a PhD Candidate at the Center for Astrophysics | Harvard & Smithsonian

I originally come from Sofia, Bulgaria, and am currently a graduate student at the Center for Astrophysics (CfA), working with Prof. Daniel Eisenstein and Prof. Lars Hernquist on various concepts in (mostly) theoretical cosmology.

The research questions that keep me up at night can be divided into the following broad categories: how much can we learn about the galaxy-halo connection from state-of-the-art cosmological hydrodynamical simulations, how can we apply this knowledge in the most optimal and accurate way to N-body simulations, what innovative methods can we use to improve our inference on cosmological parameters from galaxy survey data, and how can we combine galaxy clustering and lensing data with early Universe probes such as the Cosmic Microwave Background (CMB) to place better constraints on cosmological parameters.

In my many hours of free time, I like biking to quiet places in nature, singing at the graduate student choir (and outside it, too), reading classic novels, learning new languages, supporting social justice issues as best as I can, and eating lots of desserts! I am also a member of the LGBTQ+ community and am open to talk to anyone about it!

  • Name:Boryana Hadzhiyska
  • Email:boryana.hadzhiyska@
  • Research Interests:Theoretical Cosmology, Large-scale Structure, N-Body Simulations, Galaxy-halo Connection
  • Advisors:Daniel Eisenstein, Lars Hernquist
  • Education:
    2022 - PhD in Astrophysics, Harvard University
    2018 - MASt, Cambridge University
    2017 - BA, Princeton University
Download CV





Apr 2020

Limitations to the `basic' HOD model and beyond

BH, Sownak Bose, Daniel Eisenstein, Lars Hernquist, David N. Spergel

We make use of the IllustrisTNG hydrodynamical simulations to test fundamental assumptions of the Halo Occupation Distribution (HOD) approach to modelling the galaxy-halo connection. By comparing the clustering of galaxies measured in TNG300 with that predicted by the standard HOD model, we find that, on average, the "basic" HOD model underpredicts the real-space correlation function by ∼15% on large scales, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. In our exploration of "secondary" parameters with which to augment the "basic" HOD, we find that the local environment of the halo and the velocity dispersion anisotropy help the most with reconciling the HOD clustering with that of TNG300.

Oct 2021

CompaSO: A new halo finder for competitive assignment to spherical overdensities

BH, Daniel Eisenstein, Sownak Bose, Lehman Garrison, Nina Maksimova

We describe a new method called CompaSO for identifying groups of particles in N-body simulations. This halo-finding algorithm is used in the AbacusSummit suite. CompaSO builds upon existing spherical overdensity algorithms by taking into consideration the tidal radius around a smaller halo before competitively assigning halo membership to the particles.

Feb 2021

The galaxy-halo connection of emission-line galaxies in IllustrisTNG

BH, Sandro Tacchella, Sownak Bose, Daniel Eisenstein

We employ the hydrodynamical simulation IllustrisTNG to explore the halo occupation distribution (HOD) and environmental dependence of luminous star-forming emission-line galaxies (ELGs) at z∼1. We demonstrate that the ELG populations are twice more likely to reside in lower-density regions (sheets) compared with the mass-selected populations and twice less likely to occupy the densest regions of the cosmic web (knots). We estimate the assembly bias signal at ∼4% for the DESI sample at z=0.8, which is less compared with that of LRG-like samples, but nevertheless, requiring attention.

Oct 2020

Analytic marginalization of N(z) uncertainties in tomographic galaxy surveys

BH, David Alonso, Andrina Nicola, Anže Slosar

We present a new method to marginalize over uncertainties in redshift distributions, N(z), within tomographic cosmological analyses applicable to current and upcoming photometric galaxy surveys. We allow for arbitrary deviations from the best-guess N(z) governed by a general covariance matrix describing the uncertainty in our knowledge of redshift distributions.

Feb 2021

Extensions to models of the galaxy-halo connection

BH, Sownak Bose, Daniel Eisenstein, Lars Hernquist

We explore two widely used empirical models for the galaxy-halo connection, subhalo abundance matching (SHAM) and the halo occupation distribution (HOD) and compare their predictions with the hydrodynamical simulation IllustrisTNG (TNG) for a range of statistics that quantify the galaxy distribution. We observe that in their most straightforward implementations, both models fail to reproduce the two-point clustering measured in TNG. We also consider a two-dimensional implementation of the HOD model augmented with environment and velocity anisotropy and find both samples to be consistent with full-physics across all statistical probes examined.

Oct 2021

The halo light cone catalogs of the AbacusSummit suite of simulations

BH, Daniel Eisenstein, Sownak Bose, Lehman Garrison, Nina Maksimova

We describe a method for generating halo light cones using the AbacusSummit suite. The halo light cone catalogs for base simulation (mass resolution of 2 10^9 Msun/h and box size of 2 Gpc/h) cover an octant of the sky to z = 0.8 and about 1800 deg^2, at z = 2.5. We perform various tests and to validate its accuracy and compare the clustering properties of light cone mock catalogs with that of snapshot catalogs.

July 2021

Galaxy assembly bias: a comparative study between IllustrisTNG and a semi-analytical model

BH, Sonya Liu, Rachel Somerville, Sownak Bose, Daniel Eisenstein, Lars Hernquist

We compare the hydrodynamical simulation IllustrisTNG with a semi-analytical model calibrated to match various galaxy properties. We study LRG- and ELG-like samples at redshifts z = 0 and z = 0.8, comparing them in terms of their large-scale structure properties such as halo occupation distribution, two-point clustering of the galaxies, void-galaxy correlation function, cumulative galaxy density, and galaxy-galaxy lensing. We further study the assembly bias signal in both TNG and the SAM and trace its origins iin both in an effort to explain the observed discrepancies.

March 2021

Hefty enhancement of cosmological constraints from the DES Y1 data using a Hybrid Effective Field Theory approach to galaxy bias

BH, Carlos Garcia-Garcia, David Alonso, Andrina Nicola, Anže Slosar

We present a re-analysis of the shear and galaxy clustering data from first-year Dark Energy Survey data (DES Y1), making use of a Hybrid Effective Field Theory (HEFT) approach to model the galaxy-matter relation on weakly non-linear scales, initially proposed in Modi et al. (2020). This allows us to explore the enhancement in cosmological constraints enabled by extending the galaxy clustering scale range typically used in projected large-scale structure analysis.








A-204, MS 10,
60 Garden Street,
Cambridge, MA 02138

(617) 384 9665


Follow Me

Get in Touch