Source author record

Mohammadtaher Safarzadeh

Mohammadtaher Safarzadeh appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

24 published item(s)

preprint2022arXiv

A Statistical Detection of Wide Binary Systems in the Ultra-Faint Dwarf Galaxy Reticulum II

Binary stars can inflate the observed velocity dispersion of stars in dark matter dominated systems such as ultra-faint dwarf galaxies (UFDs). However, the population of binaries in UFDs is poorly constrained by observations, with preferred binary fractions for individual galaxies ranging from a few percent to nearly unity. Searching for wide binaries through nearest neighbor (NN) statistics (or the two-point correlation function) has been suggested in the literature, and we apply this method for the first time to detect wide binaries in a UFD. By analyzing the positions of stars in Reticulum II (Ret II) from Hubble Space Telescope images, we search for angularly resolved wide binaries in Ret II. We find that the distribution of their NN distances shows an enhancement at projected separations of $<1$ arcsec relative to a model containing no binaries. We show that such an enhancement can be explained by a wide~binary fraction of $f_b\approx0.007^{+0.008}_{-0.003}$ at separations of more than 3000 AU. Under the assumption that the binary separation distribution is similar to that in the Milky Way, the total binary fraction in Ret II may be on the order of 50%. We also use the observed magnitude distribution of stars in Ret II to constrain the initial mass function over the mass range $0.34-0.78~M_{\odot}$, finding that a shallow power-law slope of $1.01 \le α\le 1.15$ matches the data.

preprint2022arXiv

Interpreting a Machine Learning Model for Detecting Gravitational Waves

We describe a case study of translational research, applying interpretability techniques developed for computer vision to machine learning models used to search for and find gravitational waves. The models we study are trained to detect black hole merger events in non-Gaussian and non-stationary advanced Laser Interferometer Gravitational-wave Observatory (LIGO) data. We produced visualizations of the response of machine learning models when they process advanced LIGO data that contains real gravitational wave signals, noise anomalies, and pure advanced LIGO noise. Our findings shed light on the responses of individual neurons in these machine learning models. Further analysis suggests that different parts of the network appear to specialize in local versus global features, and that this difference appears to be rooted in the branched architecture of the network as well as noise characteristics of the LIGO detectors. We believe efforts to whiten these "black box" models can suggest future avenues for research and help inform the design of interpretable machine learning models for gravitational wave astrophysics.

preprint2021arXiv

How can LISA probe a population of GW190425-like binary neutron stars in the Milky Way?

The nature of GW190425, a presumed binary neutron star (BNS) merger detected by the LIGO/Virgo Scientific Collaboration (LVC) with a total mass of $3.4^{+0.3}_{-0.1}$ M$_{\odot}$, remains a mystery. With such a large total mass, GW190425 stands at five standard deviations away from the total mass distribution of Galactic BNSs of $2.66\pm 0.12$ M$_{\odot}$. LVC suggested that this system could be a BNS formed from a fast-merging channel rendering its non-detection at radio wavelengths due to selection effects. BNSs with orbital periods less than a few hours - progenitors of LIGO/Virgo mergers - are prime target candidates for the future Laser Interferometer Space Antenna (LISA). If GW190425-like binaries exist in the Milky Way, LISA will detect them within the volume of our Galaxy and will measure their chirp masses to better than 10% for those binaries with gravitational wave frequencies larger than 2 mHz. This work explores how we can probe a population of Galactic GW190425-like BNSs with LISA and investigate their origin. We assume that the Milky Way's BNS population consists of two distinct sub-populations: a fraction $w_1$ that follows the observed Galactic BNS chirp mass distribution and $w_2$ that resembles chirp mass of GW190425. We show that LISA's accuracy on recovering the fraction of GW190425-like binaries depends on the BNS merger rate. For the merger rates reported in the literature, $21 - 212\,$Myr$^{-1}$, the error on the recovered fractions varies between $\sim 30 - 5$%.

preprint2020arXiv

A trend in the effective spin distribution of LIGO binary black holes with mass

Binary black holes (BBHs) detected by gravitational wave (GW) observations could be broadly divided into two formation channels: those formed through field binary evolution and those assembled dynamically in dense stellar systems. Each of these formation channels, and their sub-channels, populate a distinct region in the effective spin-mass ($χ_{\rm eff}-M$) plane. Depending on the branching ratio of different channels, an ensemble of BBHs could show a trend in this plane. Here we fit a mass-dependent distribution for $χ_{\rm eff}$ to the GWTC-1 BBHs from the first and second observing runs of Advanced LIGO and Advanced Virgo. We find a negative correlation between mass and the mean effective spin ($\barχ_{\mathrm{eff}}$), and positive correlation with its dispersion ($σ_{χ_\mathrm{eff}}$) at 75\% and 80\% confidence. This trend is robust against the choice of mass variable, but most pronounced when the mass variable is taken to be the chirp mass of the binary. The result is consistent with significant contributions from both dynamically assembled and field binaries in the GWTC-1 catalog. The upcoming LIGO O3a data release will critically test this interpretation.

preprint2020arXiv

Being Careful with the Field Formation Interpretation of GW190412

The LIGO/Virgo Scientific Collaboration recently announced the detection of a compact object binary merger, GW190412, as the first asymmetric binary black hole (BBH) merger with mass ratio $q\approx0.25$. Other than the mass ratio, this BBH has shown to have a positive effective spin of around $χ_{\rm eff}\approx0.28$. Assuming a field formation channel, associating this effective spin to either the primary or the secondary BH each has its implications: If the spin of the BBH comes form the primary BH, it has consequences for the efficiency of angular momentum transport in the formation of the BH. If, on the other hand, the spin is due to the secondary BH through tidal spin-up processes, one has to note that (i) such processes have very short delay-times, and (ii) subsequently, their local merger rate is determined by local star formation rate at assumed formation metallicity of the BBH. We show that the predicted merger rate density from this channel is $\lesssim 0.3~\rm Gpc^{-3} yr^{-1}$ and in tension with the rather high local merger rate of such systems which we estimate from this single event to be $\sim 1.7^{+2.5}_{-1.4}~\rm Gpc^{-3} yr^{-1}$ (90\% confidence interval, and assuming 50 days of observing time). Large natal kicks ($v\gtrsim 500\,{\rm km/s}$) would be required to get such BBHs with an in-plane spin component to account for the marginal detection of precession in GW190412. However, this would only exacerbate the tension as the estimated local merger rate would be further decreased. Similarly, the formation of such systems through the dynamical assembly is exceedingly rare, leaving this system a dilemma hard to account for with the currently accepted paradigms of BBH formation.

preprint2020arXiv

Confronting the Magnetar Interpretation of Fast Radio Bursts Through Their Host Galaxy Demographics

We explore the prompt magnetar progenitor scenario in the context of fast radio burst (FRB) host galaxies demographics and offset distributions. Magnetars are neutron stars with strong magnetic fields on the order of $10^{15}$ G with a short decay lifetime of less than $10^4$ years. Due to their extremely short lifetimes, magnetars should follow the demographics of galaxies according to their current star-formation rate (SFR). Moreover, we hypothesize that magnetars should follow the SFR profile within galaxies, which we assume to follow an exponential profile. We construct a simple model for the host galaxies of magnetars assuming these events track SFR in all galaxies and compare it to observed properties from a sample of \nsecure\ secure FRB hosts. We find the distribution of observed SFRs is inconsistent with the model at $>95\%$ c.l. The offset distribution is consistent with this scenario; however, this could be due to the limited sample size and the seeing limited estimates for the effective radii of the FRB host galaxies. Despite the recent association of an FRB with a magnetar in the Milky Way, magnetars may not be the only source of FRBs in the universe, yet any other successful model must account for the demographics of the FRB host in SFR and their observed galactocentric offsets.

preprint2020arXiv

Does GW190425 require an alternative formation pathway than a fast-merging channel?

The LIGO/Virgo Scientific Collaboration (LVC) recently announced the detection of a compact object binary merger, GW190425, with a total mass of $3.4^{+0.3}_{-0.1}$ M$_{\odot}$, and individual component masses in the range of about 1.1 to 2.5 $M_{\odot}$. If the constituent compact objects are neutron stars, then the total mass is five standard deviations higher than the mean of $2.66\pm 0.12$ M$_{\odot}$ for Galactic binary neutron stars (BNSs). LVC suggests such massive BNS systems are born from a fast-merging channel, and therefore, their non-detection in the Galaxy to be due to a selection effect. However, we are unable to reconcile the inferred formation efficiency from the reported merger rate, $\mathcal{R}_{\rm GW190425}=460^{+1050}_{-390}$ yr$^{-1}$ Gpc$^{-3}$, with predictions from our own study for fast-merging BNS systems. Moreover, the comparable merger rates of GW190425, and GW170817 are possibly in tension with our results for two reasons: (i) more massive systems are expected to have a lower formation rate, and (ii) fast merging channels should constitute $\lesssim 10\%$ of the total BNS systems if case BB unstable mass transfer is permitted to take place as a formation pathway. We argue that to account for the high merger rate of GW190425 as a BNS system requires: (i) revisiting our understanding of NS formation in supernova explosions, or (ii): that more massive NSs need to be preferentially born with either very weak or very high magnetic fields so that they would be undetectable in the radio surveys. Perhaps massive NSs detected in NS-white dwarf binaries are our clues to the formation path of GW190425 systems.

preprint2020arXiv

Formation of mass gap objects in highly asymmetric mergers

The LIGO/Virgo Collaboration (LVC) recently reported the detection of GW190814, a merger of a $23^{+1.0}_{-0.9}~M_{\odot}$ primary black hole (BH), and a $2.6^{+0.08}_{-0.08}~M_{\odot}$ secondary. The secondary's mass falls into the mass-gap regime, which refers to the scarcity of compact objects in the mass range of 2-5 $M_{\odot}$. The first clue to the formation of the GW190814 lies in the fact that the primary is a very massive BH. We suggest that the secondary was born as a neutron star (NS) where a significant amount of the supernova ejecta mass from its formation remained bound to the binary due to the presence of the massive BH companion. The bound mass forms a circumbinary accretion disk, and its accretion onto the NS created a mass-gap object. In this scenario, LIGO/Virgo will only detect mass-gap objects in binary mergers with an extreme mass ratio. We also predict a correlation between the mass of the secondary and the mass of the primary in such asymmetric mergers. Our model can be tested with future data from the LVC's third-observing run.

preprint2020arXiv

Modeling Photoionized Turbulent Material in the Circumgalactic Medium II: Effect of Turbulence within a Stratified Medium

The circumgalactic medium (CGM) of nearby star-forming galaxies shows clear indications of OVI absorption accompanied by little to no detectable NV absorption. This unusual spectral signature, accompanied by highly non-uniform absorption from lower ionization state species, indicates that the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. Motivated by previous isotropic turbulent simulations, we carry out chemodynamical simulations of stratified media in a Navarro-Frenk-White (NFW) gravitational potential with a total mass of $10^{12}$ solar masses and turbulence that decreases radially. The simulations assume a metallicity of 0.3 solar, a redshift zero metagalatic UV background, and they track ionizations, recombinations, and species-by-species radiative cooling using the MAIHEM package. We compare a suite of ionic column densities with the COS-Halos sample of low-redshift star-forming galaxies. Turbulence with an average one-dimensional velocity dispersion approximately 40 km/s, corresponding to an energy injection rate of approximately $10^{49}$ erg/yr, produces a CGM that matches many of the observed ionic column densities and ratios. In this simulation, the NVI to OVI ratio is suppressed from its equilibrium value due to a combination of radiative cooling and cooling from turbulent mixing. This level of turbulence is consistent with expectations from observations of better constrained, higher-mass systems, and could be sustained by energy input from supernovae, gas inflows, and dynamical friction from dark matter subhalos. We also conduct a higher resolution run which yields smaller-scale structures, but remains in agreement with observations.

preprint2020arXiv

Numerical Relativity Simulations of the Neutron Star Merger GW170817: Long-Term Remnant Evolutions, Winds, Remnant Disks, and Nucleosynthesis

We present a systematic numerical-relativity study of the dynamical ejecta, winds and nucleosynthesis in neutron star merger remnants. Binaries with the chirp mass compatible with GW170817, different mass ratios, and five microphysical equations of state (EOS) are simulated with an approximate neutrino transport and a subgrid model for magnetohydrodynamics turbulence up to 100 milliseconds postmerger. Spiral density waves propagating from the neutron star remnant to the disk trigger a wind with mass flux ${\sim}0.1{-}0.5\,{\rm M_\odot/s}$ persisting for the entire simulation as long as the remnant does not collapse to black hole. This wind has average electron fraction $\gtrsim 0.3$ and average velocity ${\sim}0.1-0.17\,$c and thus is a site for the production of weak $r$-process elements (mass number $A<195$). Disks around long-lived remnants have masses ${\sim}0.1{-}0.2\,{\rm M_\odot}$, temperatures peaking at $\lesssim10\,$MeV near the inner edge, and a characteristic double-peak distribution in entropy resulting from shocks propagating through the disk. The dynamical and spiral-wave ejecta computed in our targeted simulations are not compatible with those inferred from AT2017gfo using two-components kilonova models. Rather, they indicate that multi-component kilonova models including disk winds are necessary to interpret AT2017gfo. The nucleosynthesis in the combined dynamical ejecta and spiral-wave wind in the comparable-mass long-lived mergers robustly accounts for all the $r$-process peaks, from mass number ${\sim}75$ to actinides in terms of solar abundances. Total abundandes are weakly dependent on the EOS, while the mass ratio affect the production of first peak elements.

preprint2020arXiv

The branching ratio of LIGO binary black holes

Formation of binary black holes (BBHs) detected by gravitational-wave (GW) observations could be broadly divided into two categories: those formed through field binary evolution and those assembled dynamically in dense stellar systems. The branching ratio of the BBHs refers to the contribution of each channel. The dynamical assembly channel would predict a symmetric distribution in the effective spins of the BBHs while field formation predicts BBHs to have positive effective spins. By modeling these two populations based on their effective spin distribution we show that in the 10 BBHs detected by LIGO/Virgo the contribution of the dynamically assembled BBHs to be more than about 50\% with 90\% confidence. This result is based on the assumption that the field binaries are born with positive effective spins not restricted to have small values.

preprint2020arXiv

The nearest discovered black hole is likely not in a triple configuration

HR 6819 was recently claimed to be a hierarchical triple system of a Be star in a wide orbit around an inner binary system of a black hole (BH) and a B III type star. We argue that this system is unlikely to be a hierarchical triple due to three reasons: (i) Given that this system is discovered in a magnitude limited Bright Star Catalog, the expected number of such systems in the Milky Way amounts to about $10^4$ while the estimate for the MW budget for such systems is between $10^2-10^3$ systems under generous assumptions. Such a large gap cannot be reconciled as it would otherwise likely overflow the MW budget for BHs; (ii) The dynamical stability of this system sets lower bounds on the orbital separation of the outer Be star, while it not being resolved by Gaia places an upper limit on its projected sky separation. We show that these two constraints would imply a narrow range for the outer orbit without resorting to geometrical fine-tuning; (iii) The triple system should have survived the stellar evolution prior to the formation of the BH in the inner binary. We perform numerical simulations starting with conservative initial conditions of this system and show that a small parameter space for BH progenitor star's mass loss, BH natal kicks, and initial orbital separation can reproduce HR 6819. Therefore, we propose this system is a chance superposition of a Be star with a binary.

preprint2020arXiv

Ultra-light Dark Matter is Incompatible with the Milky Way's Dwarf Satellites

The density profiles of dwarf galaxies are a highly varied set. If the dark matter is an ultra-light particle such as axions, then simulations predict a distinctive and unique profile. If the axion mass is large enough to fit the ultra-faint dwarf (UFD) satellites($m\gtrapprox 10^{-21}$ eV), then the models do not fit the density profile of Fornax and Sculptor and are ruled out by more than $3-σ$ confidence. If the axion mass is in the mass range that can fit mass profiles of Fornax and Sculptor dwarf spheroidals, then its extended profile implies enormous masses ($\approx10^{11}-10^{12}M_{\odot}$) for the UFDs. These large masses for the UFDS are ruled out by more than $3-σ$ confidence by dynamical friction arguments. The tension would increase further considering star formation histories and stellar masses of the UFDs. Unless future ultra-light dark matter (ULDM) simulations with baryonic feedback show a significant change in the density structure of the halos, the current data is incompatible with the ULDM scenario. Relaxing the slope constraint from classical dwarf galaxies would lead to excluding ULDM with mass less than $6\times10^{-22}$ eV.

preprint2020arXiv

Was GW190412 born from a hierarchical 3+1 quadruple configuration?

The gravitational wave source GW190412 is a binary black hole (BBH) merger with three unique properties: i) its mass ratio is about 0.28, the lowest found so far, ii) it has a relatively high positive effective spin parameter χ_eff=0.25, and iii) it is observed to be precessing due to in-plane projected spin of the binary with an in-plane precession parameter χ_p=0.3. The two main formation channels of BBH formation fail to account for GW190412: field formation scenarios cannot explain the observed precession unless by invoking large natal kicks, and dynamical assembly in dense stellar systems is inefficient in producing such low mass-ratio BBH mergers. Here, we investigate whether "double mergers" in wide hierarchical quadruple systems in the "3+1" configuration could explain the unique properties of GW190412. In this scenario, a compact object quadruple system experiences two mergers: first, two compact objects in the innermost orbit merge due to secular chaotic evolution. At a later time, the merged compact object coalesces with another compact object due to secular Lidov-Kozai oscillations. We find that our scenario is consistent with GW190412. In particular, we find a preferential projected spin around χ_p=0.2. However, the likelihood of a double merger is small and the formation efficiency of these systems is uncertain. If GW190412 originated from a double merger in a 3+1 quadruple, we find a strong constraint that the first merger likely occurred between roughly equal-mass BHs in the innermost orbit, since the recoil velocity from unequal-mass BHs would otherwise have disrupted the system.

preprint2019arXiv

Formation and merging of Mass Gap Black Holes in Gravitational Wave Merger Events from Wide Hierarchical Quadruple Systems

We investigate secular evolution in hierarchical quadruple systems as a formation channel of mass-gap black holes (with masses of about $3-5$ $M_{\odot}$) in systems that will eventually lead to binary black hole mergers detectable by ground-based gravitational wave detectors (LIGO/Virgo). We show that in a 3+1 hierarchical system, two episodes of induced mergers would first cause two neutron stars to merge and form a mass-gap black hole, which will subsequently merge with another (more massive) black hole through a second induced merger. We demonstrate that such systems are stable to flybys, and their formation would predict a high mass ratio and eccentric merger of a mass-gap black hole with a more massive black hole companion. Such a formation channel may explain observed gravitational wave events such as the recently-discovered LIGO/Virgo events S190814bv and S190924h.

preprint2019arXiv

LISA and the Existence of a Fast-Merging Double Neutron Star Formation Channel

Using a Milky Way double neutron star (DNS) merger rate of 210 Myr$^{-1}$, as derived by the Laser Interferometer Gravitational-Wave Observatory (LIGO), we demonstrate that the Laser Interferometer Space Antenna (LISA) will detect on average 240 (330) DNSs within the Milky Way for a 4-year (8-year) mission with a signal-to-noise ratio greater than 7. Even adopting a more pessimistic rate of 42 Myr$^{-1}$, as derived by the population of Galactic DNSs, we find a significant detection of 46 (65) Milky Way DNSs. These DNSs can be leveraged to constrain formation scenarios. In particular, traditional NS-discovery methods using radio telescopes are unable to detect DNSs with $P_{\rm orb}$ $\lesssim$1 hour (merger times $\lesssim$10 Myr). If a fast-merging channel exists that forms DNSs at these short orbital periods, LISA affords, perhaps, the only opportunity to observationally characterize these systems; we show that toy models for possible formation scenarios leave unique imprints on DNS orbital eccentricities, which may be measured by LISA for values as small as $\sim$10$^{-2}$.

preprint2016arXiv

The Impact of Unresolved Turbulence on the Escape Fraction of Lyman Continuum Photons

We investigate the relation between the turbulent Mach number (\mach) and the escape fraction of Lyman continuum photons ($f_{\rm esc}$) in high-redshift galaxies. Approximating the turbulence as isothermal and isotropic, we show that the increase in the variance in column densities from $\mathcal{M}=1$ to $\mathcal{M}=10$ causes $f_{\rm esc}$ to increase by $\approx 25$\%, and the increase from $\mathcal{M}=1$ to $\mathcal{M}=20$ causes $f_{\rm esc}$ to increases by $\approx 50$\% for a medium with opacity $τ\approx1$. At a fixed Mach number, the correction factor for escape fraction relative to a constant column density case scales exponentially with the opacity in the cell, which has a large impact for simulated star forming regions. Furthermore, in simulations of isotropic turbulence with full atomic/ionic cooling and chemistry, the fraction of HI drops by a factor of $\approx 2.5$ at $\mathcal{M}\approx10$ even when the mean temperature is $\approx5\times10^3 K$. If turbulence is unresolved, these effects together enhance $f_{\rm esc}$ by a factor $>3$ at Mach numbers above 10. Such Mach numbers are common at high-redshifts where vigorous turbulence is driven by supernovae, gravitational instabilities, and merger activity, as shown both by numerical simulations and observations. These results, if implemented in the current hydrodynamical cosmological simulations to account for unresolved turbulence, can boost the theoretical predictions of the Lyman Continuum photon escape fraction and further constrain the sources of reionization.

preprint2016arXiv

UV Absorption Line Ratios in Circumgalactic Medium at Low Redshift in Realistic Cosmological Hydrodynamic Simulations

Utilizing high-resolution cosmological hydrodynamic simulations we investigate various ultra-violet absorption lines in the circumgalactic medium of star forming galaxies at low redshift, in hopes of checking and alleviating the claimed observational conundrum of the ratio of NV to OVI absorbers, among others. We find a satisfactory agreement between simulations and extant observational data with respect to the ratios of the following four line pairs examined, NV/OVI, SiIV/OVI, NIII/OVI and NII/OVI. For the pairs involving nitrogen lines, we examine two cases of nitrogen abundance, one with constant N/O ratio and the other with varying N/O ratio, with the latter motivated by theoretical considerations of two different synthetic sources of nitrogen that is empirically verified independently. Along a separate vector, for all line pairs, we examine two cases of radiation field, one with the Haardt-Madau background radiation field and the other with an additional local radiation field sourced by hot gas in the host galaxy. In all cases, two-sample Kolmogorov-Smirnov tests indicate excellent agreements. We find that the apparent agreements between simulations and observations will be strongly tested, if the bulk of current upper limits of various line ratios are turned into actual detections. We show that an increase in observational sensitivity by 0.2 dex will already start to significantly constrain the models.

preprint2015arXiv

Testing Dark Matter Halo Models of Quasars With Thermal Sunyaev-Zeldovich Effect

A statistical analysis of stacked Compton$-y$ maps of quasar hosts with a median redshift of $1.5$ using Millennium Simulation is performed to address two issues, one on the feedback energy from quasars and the other on testing dark matter halo models for quasar hosts. On the first, we find that, at the resolution of FWHM=$10$ arcmin obtained by Planck data, the observed thermal Sunyaev-Zeldovich (tSZ) effect can be entirely accounted for and explained by the thermal energy of halos sourced by gravitational collapse of halos, without a need to invoke additional, large energy sources, such as quasar or stellar feedback. Allowing for uncertainties of dust temperature in the calibration of observed Comton$-y$ maps, the maximum additional feedback energy is $\sim 25\%$ of that previously suggested. Second, we show that, with FWHM=$1$ arcmin beam, tSZ measurements will provide a potentially powerful test of quasar-hosting dark matter halo models, limited only by possible observational systematic uncertainties, not by statistical ones, even in the presence of possible quasar feedback.

preprint2015arXiv

What shapes the far-infrared spectral energy distributions of galaxies?

To explore the connection between the global physical properties of galaxies and their far-infrared (FIR) spectral energy distributions (SEDs), we study the variation in the FIR SEDs of a set of hydrodynamically simulated galaxies that are generated by performing dust radiative transfer in post-processing. Our sample includes both isolated and merging systems at various stages of the merging process and covers infrared (IR) luminosities and dust masses that are representative of both low- and high-redshift galaxies. We study the FIR SEDs using principle component analysis (PCA) and find that 97\% of the variance in the sample can be explained by two principle components (PCs). The first PC characterizes the wavelength of the peak of the FIR SED, and the second encodes the breadth of the SED. We find that the coefficients of both PCs can be predicted well using a double power law in terms of the IR luminosity and dust mass, which suggests that these two physical properties are the primary determinants of galaxies' FIR SED shapes. Incorporating galaxy sizes does not significantly improve our ability to predict the FIR SEDs. Our results suggest that the observed redshift evolution in the effective dust temperature at fixed IR luminosity is not driven by geometry: the SEDs of $z \sim 2-3$ ultraluminous IR galaxies (ULIRGs) are cooler than those of local ULIRGs not because the high-redshift galaxies are more extended but rather because they have higher dust masses at fixed IR luminosity. Finally, based on our simulations, we introduce a two-parameter set of SED templates that depend on both IR luminosity and dust mass.

preprint2014arXiv

A New Model for Dark Matter Halos Hosting Quasars

A new model for quasar-hosting dark matter halos, meeting two physical conditions, is put forth. First, significant interactions are taken into consideration to trigger quasar activities. Second, satellites in very massive halos at low redshift are removed from consideration, due to their deficiency of cold gas. We analyze the {\em Millennium Simulation} to find halos that meet these two conditions and simultaneously match two-point auto-correlation functions of quasars and cross-correlation functions between quasars and galaxies at $z=0.5-3.2$. %The found halos have some distinct properties worth noting. The masses of found quasar hosts decrease with decreasing redshift, with the mass thresholds being $[(2-5)\times 10^{12}, (2-5)\times 10^{11}, (1-3)\times 10^{11}]\msun$ for median luminosities of $\sim[10^{46}, 10^{46}, 10^{45}]$erg/s at $z=(3.2, 1.4, 0.53)$, respectively, an order of magnitude lower than those inferred based on halo occupation distribution modeling. In this model quasar hosts are primarily massive central halos at $z\ge 2-3$ but increasingly dominated by lower mass satellite halos experiencing major interactions towards lower redshift. But below $z=1$ satellite halos in groups more massive than $\sim 2\times 10^{13}\msun$ do not host quasars. Whether for central or satellite halos, imposing the condition of significant interactions substantially boosts the clustering strength compared to the total population with the same mass cut. The inferred lifetimes of quasars at $z=0.5-3.2$ of $3-30$Myr are in agreement with observations. Quasars at $z\sim 2$ would be hosted by halos of mass $\sim 5\times 10^{11}\msun$ in this model, compared to $\sim 3\times 10^{12}\msun$ previously thought, which would help reconcile with the observed, otherwise puzzling high covering fractions for Lyman limit systems around quasars.

preprint2014arXiv

De-Confusing blended field images using graphs and bayesian priors

We present a new technique for overcoming confusion noise in deep far-infrared \Herschel space telescope images making use of prior information from shorter $λ<2$\micron wavelengths. For the deepest images obtained by \Herschels, the flux limit due to source confusion is about a factor of three brighter than the flux limit due to instrumental noise and (smooth) sky background. We have investigated the possibility of de-confusing simulated \Herschel PACS-160\micron images by using strong Bayesian priors on the positions and weak priors on the flux of sources. We find the blended sources and group them together and simultaneously fit their fluxes. We derive the posterior probability distribution function of fluxes subject to these priors through Monte Carlo Markov Chain (MCMC) sampling by fitting the image. Assuming we can predict FIR flux of sources based on ultraviolet-optical part of their SEDs to within an order of magnitude, the simulations show that we can obtain reliable fluxes and uncertainties at least a factor of three fainter than the confusion noise limit of $3σ_{c} $=2.7 mJy in our simulated PACS-160 image. This technique could in principle be used to mitigate the effects of source confusion in any situation where one has prior information of positions and plausible fluxes of blended sources. For \Herschel, application of this technique will improve our ability to constrain the dust content in normal galaxies at high redshift.

preprint2013arXiv

A Comparison between Semi-Analytic Model Predictions for the CANDELS Survey

We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED

preprint2011arXiv

Spinning up black holes with super-critical accretion flows

We study the process of spinning up black holes by accretion from slim disks in a wide range of accretion rates. We show that for super-Eddington accretion rates and low values of the viscosity parameter alpha (\lesssim 0.01) the limiting value of the dimensionless spin parameter a_* can reach values higher than a_*=0.9978 inferred by Thorne (1974) in his seminal study. For \mdot=10\mdot_{Edd} and alpha=0.01 spin equilibrium is reached at a_*=0.9994. We show that the equilibrium spin value depends strongly on the assumed value of alpha. We also prove that for high accretion rates the impact of captured radiation on spin evolution is negligible.