Source author record

Feryal Özel

Feryal Özel 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

13works
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

13 published item(s)

preprint2022arXiv

MeqSilhouette v2: Spectrally-resolved polarimetric synthetic data generation for the Event Horizon Telescope

We present MeqSilhouette v2.0 (MeqSv2), a fully polarimetric, time-and frequency-resolved synthetic data generation software for simulating millimetre (mm) wavelength very long baseline interferometry (VLBI) observations with heterogeneous arrays. Synthetic data are a critical component in understanding real observations, testing calibration and imaging algorithms, and predicting performance metrics of existing or proposed sites. MeqSv2 applies physics-based instrumental and atmospheric signal corruptions constrained by empirically-derived site and station parameters to the data. The new version is capable of applying instrumental polarization effects and various other spectrally-resolved effects using the Radio Interferometry Measurement Equation (RIME) formalism and produces synthetic data compatible with calibration pipelines designed to process real data. We demonstrate the various corruption capabilities of MeqSv2 using different arrays, with a focus on the effect of complex bandpass gains on closure quantities for the EHT at 230 GHz. We validate the frequency-dependent polarization leakage implementation by performing polarization self-calibration of synthetic EHT data using PolSolve. We also note the potential applications for cm-wavelength VLBI array analysis and design and future directions.

preprint2020arXiv

A Parametric model for the shapes of black-hole shadows in non-Kerr spacetimes

The Event Horizon Telescope (EHT) is taking the first images of black holes resolved at horizon scales to measure their shadows and probe accretion physics. A promising avenue for testing the hypothesis that astrophysical black holes are described by the Kerr solution to Einstein's equations is to compare the size and shape of the shadow a black hole casts on the surrounding emission to the predictions of the Kerr metric. We develop here an efficient parametric framework to perform this test. We carry out ray-tracing simulations for several parametrized non-Kerr metrics to create a large data set of non-Kerr shadows that probe the allowed parameter space for the free parameters of each metric. We then perform principal components analysis (PCA) on this set of shadows and show that only a small number of components are needed to accurately reconstruct all shadows within the set. We further show that the amplitude of the PCA components are smoothly related to the free parameters in the metrics, and therefore, that these PCA components can be fit to EHT observations in order to place constraints on the free parameters of these metrics that will help quantify any potential deviations from the Kerr solution.

preprint2020arXiv

Spectral analysis of the quiescent low-mass X-ray binary in the globular cluster M30

We present a recent Chandra observation of the quiescent low-mass X-ray binary containing a neutron star, located in the globular cluster M30. We fit the thermal emission from the neutron star to extract its mass and radius. We find no evidence of flux variability between the two observations taken in 2001 and 2017, nor between individual 2017 observations, so we analyse them together to increase the signal to noise. We perform simultaneous spectral fits using standard light-element composition atmosphere models (hydrogen or helium), including absorption by the interstellar medium, correction for pile-up of X-ray photons on the detector, and a power-law for count excesses at high photon energy. Using a Markov-chain Monte Carlo approach, we extract mass and radius credible intervals for both chemical compositions of the atmosphere: $R_{\textrm{NS}}=7.94^{+0.76}_{-1.21}$ km and $M_{\textrm{NS}}<1.19$ M$_{\odot}$ assuming pure hydrogen, and $R_{\textrm{NS}}=10.50^{+2.88}_{-2.03}$ km and $M_{\textrm{NS}}<1.78$ M$_{\odot}$ for helium, where the uncertainties represent the 90% credible regions. For H, the small radius is difficult to reconcile with most current nuclear physics models (especially for nucleonic equations of state) and with other measurements of neutron star radii, with recent preferred values generally in the 11-14 km range. Whereas for He, the measured radius is consistent with this range. We discuss possible sources of systematic uncertainty that may result in an underestimation of the radius, identifying the presence of surface temperature inhomogeneities as the most relevant bias. According to this, we conclude that either the atmosphere is composed of He, or it is a H atmosphere with a significant contribution of hot spots to the observed radiation.

preprint2016arXiv

Bayesian techniques for comparing time-dependent GRMHD simulations to variable Event Horizon Telescope observations

The Event Horizon Telescope (EHT) is a millimeter-wavelength, very-long-baseline interferometry (VLBI) experiment that is capable of observing black holes with horizon-scale resolution. Early observations have revealed variable horizon-scale emission in the Galactic Center black hole, Sagittarius A* (Sgr A*). Comparing such observations to time-dependent general relativistic magnetohydrodynamic (GRMHD) simulations requires statistical tools that explicitly consider the variability in both the data and the models. We develop here a Bayesian method to compare time-resolved simulation images to variable VLBI data, in order to infer model parameters and perform model comparisons. We use mock EHT data based on GRMHD simulations to explore the robustness of this Bayesian method and contrast it to approaches that do not consider the effects of variability. We find that time-independent models lead to offset values of the inferred parameters with artificially reduced uncertainties. Moreover, neglecting the variability in the data and the models often leads to erroneous model selections. We finally apply our method to the early EHT data on Sgr A*.

preprint2016arXiv

From Neutron Star Observables to the Equation of State. I. An Optimal Parametrization

The increasing number and precision of measurements of neutron star masses, radii, and, in the near future, moments of inertia offer the possibility of precisely determining the neutron star equation of state. One way to facilitate the mapping of observables to the equation of state is through a parametrization of the latter. We present here a generic method for optimizing the parametrization of any physically allowed EoS. We use mock equations of state that incorporate physically diverse and extreme behavior to test how well our parametrization reproduces the global properties of the stars, by minimizing the errors in the observables mass, radius, and the moment of inertia. We find that using piecewise polytropes and sampling the EoS with five fiducial densities between ~1-8 times the nuclear saturation density results in optimal errors for the smallest number of parameters. Specifically, it recreates the radii of the assumed EoS to within less than 0.5 km for the extreme mock equations of state and to within less than 0.12 km for 95% of a sample of 42 proposed, physically-motivated equations of state. Such a parametrization is also able to reproduce the maximum mass to within 0.04 M_sun and the moment of inertia of a 1.338 M_sun neutron star to within less than 10% for 95% of the proposed sample of equations of state.

preprint2016arXiv

Gravitational torque-driven black hole growth and feedback in cosmological simulations

We investigate black hole-host galaxy scaling relations in cosmological simulations with a self-consistent black hole growth and feedback model. The sub-grid accretion model captures the key scalings governing angular momentum transport from galactic scales down to parsec scales, while our kinetic feedback implementation enables the injection of outflows with properties chosen to match observed nuclear outflows. We show that "quasar mode" feedback can have a large impact on the thermal properties of the intergalactic medium and the growth of galaxies and massive black holes for kinetic feedback efficiencies as low as 0.1% relative to the bolometric luminosity. Nonetheless, our simulations suggest that the black hole-host scaling relations are only weakly dependent on the effects of black hole feedback on galactic scales, owing to feedback suppressing the growth of galaxies and massive black holes by a similar amount. In contrast, the rate at which gravitational torques feed the central black hole relative to the host galaxy star formation rate governs the slope and normalization of the black hole-host correlations. Our results suggest that a common gas supply regulated by gravitational torques is the primary driver of the observed co-evolution of black holes and galaxies.

preprint2016arXiv

Measuring the neutron star equation of state using X-ray timing

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

preprint2016arXiv

Neutron Star Mass-Radius Constraints of the Quiescent Low-mass X-ray Binaries X7 and X5 in the Globular Cluster 47 Tuc

We present Chandra ACIS-S sub-array observations of the quiescent neutron star low-mass X-ray binaries X7 and X5 in the globular cluster 47 Tuc. The large reduction in photon pile-up compared to previous deep exposures enables a substantial improvement in the spectroscopic determination of the neutron star radius and mass of these neutron stars. Modeling the thermal emission from the neutron star surface with a non-magnetized hydrogen atmosphere and accounting for numerous sources of uncertainties, we obtain for the neutron star in X7 a radius of $R=11.1^{+0.8}_{-0.7}$ km for an assumed stellar mass of $M=1.4$ M$_{\odot}$ (68% C.L.). We argue, based on astrophysical grounds, that the presence of a He atmosphere is unlikely for this source. Due to eclipses and variable absorption, the quiescent low-mass X-ray binary X5 provides less stringent constraints, leading to a radius of $R=9.6^{+0.9}_{-1.1}$ km, assuming a hydrogen atmosphere and a mass of $M=1.4$ M$_{\odot}$. When combined with all other existing spectroscopic radius measurements, these measurements strongly favor radii in the 9.9-11.2 km range for a ~1.5 M$_{\odot}$ neutron star and point to a dense matter equation of state that is somewhat softer than the nucleonic ones that are consistent with laboratory experiments at low densities.

preprint2016arXiv

Persistent Asymmetric Structure of Sagittarius A* on Event Horizon Scales

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

preprint2015arXiv

Torque-Limited Growth of Massive Black Holes in Galaxies Across Cosmic Time

We combine cosmological hydrodynamic simulations with analytic models to evaluate the role of galaxy-scale gravitational torques on the evolution of massive black holes at the centers of star-forming galaxies. We confirm and extend our earlier results to show that torque-limited growth yields black holes and host galaxies evolving on average along the Mbh-Mbulge relation from early times down to z = 0 and that convergence onto the scaling relation occurs independent of the initial conditions and with no need for mass averaging through mergers or additional self-regulation processes. Smooth accretion dominates the long-term evolution, with black hole mergers with mass ratios >1:5 representing typically a small fraction of the total growth. Winds from the accretion disk are required to eject significant mass to suppress black hole growth, but there is no need for coupling this wind to galactic-scale gas to regulate black holes in a non-linear feedback loop. Torque-limited growth yields a close-to-linear relation for the star formation rate and the black hole accretion rate averaged over galaxy evolution time scales. However, the SFR-AGN connection has significant scatter owing to strong variability of black hole accretion at all resolved time scales. Eddington ratios can be described by a broad lognormal distribution with median value evolving roughly as (1 + z)^1.9, suggesting a main sequence for black hole growth similar to the cosmic evolution of specific SFRs. Our results offer an attractive scenario consistent with available observations in which cosmological gas infall and transport of angular momentum in the galaxy by gravitational instabilities regulate the long-term co-evolution of black holes and star-forming galaxies.

preprint2014arXiv

Black Hole-Galaxy Correlations without Self-Regulation

Recent models of black hole growth in a cosmological context have forwarded a paradigm in which the growth is self-regulated by feedback from the black hole itself. Here we use cosmological zoom simulations of galaxy formation down to z = 2 to show that such strong self-regulation is required in the popular spherical Bondi accretion model, but that a plausible alternative model in which black hole growth is limited by galaxy-scale torques does not require self-regulation. Instead, this torque-limited accretion model yields black holes and galaxies evolving on average along the observed scaling relations by relying only on a fixed, 5% mass retention rate onto the black hole from the radius at which the accretion flow is fed. Feedback from the black hole may (and likely does) occur, but does not need to couple to galaxy-scale gas in order to regulate black hole growth. We show that this result is insensitive to variations in the initial black hole mass, stellar feedback, or other implementation details. The torque-limited model allows for high accretion rates at very early epochs (unlike the Bondi case), which if viable can help explain the rapid early growth of black holes, while by z = 2 it yields Eddington factors of 1%-10%. This model also yields a less direct correspondence between major merger events and rapid phases of black hole growth. Instead, growth is more closely tied to cosmological disk feeding, which may help explain observational studies showing that, at least at z > 1, active galaxies do not preferentially show merger signatures.

preprint2014arXiv

Cosmological Zoom Simulations of z = 2 Galaxies: The Impact of Galactic Outflows

We use high-resolution cosmological zoom simulations with ~200 pc resolution at z = 2 and various prescriptions for galactic outflows in order to explore the impact of winds on the morphological, dynamical, and structural properties of eight individual galaxies with halo masses ~ 10^11--2x10^12 Msun at z = 2. We present a detailed comparison to spatially and spectrally resolved Hα and other observations of z ~ 2 galaxies. We find that simulations without winds produce massive, compact galaxies with low gas fractions, super-solar metallicities, high bulge fractions, and much of the star formation concentrated within the inner kpc. Strong winds are required to maintain high gas fractions, redistribute star-forming gas over larger scales, and increase the velocity dispersion of simulated galaxies, more in agreement with the large, extended, turbulent disks typical of high-redshift star-forming galaxies. Winds also suppress early star formation to produce high-redshift cosmic star formation efficiencies in better agreement with observations. Sizes, rotation velocities, and velocity dispersions all scale with stellar mass in accord with observations. Our simulations produce a diversity of morphological characteristics - among our three most massive galaxies, we find a quiescent grand-design spiral, a very compact star-forming galaxy, and a clumpy disk undergoing a minor merger; the clumps are evident in Hα but not in the stars. Rotation curves are generally slowly rising, particularly when calculated using azimuthal velocities rather than enclosed mass. Our results are broadly resolution-converged. These results show that cosmological simulations including outflows can produce disk galaxies similar to those observed during the peak epoch of cosmic galaxy growth.

preprint2012arXiv

Modeling Collisional Cascades In Debris Disks: Steep Dust-Size Distributions

We explore the evolution of the mass distribution of dust in collision-dominated debris disks, using the collisional code introduced in our previous paper. We analyze the equilibrium distribution and its dependence on model parameters by evolving over 100 models to 10 Gyr. With our numerical models, we confirm that systems reach collisional equilibrium with a mass distribution that is steeper than the traditional solution by Dohnanyi (1969). Our model yields a quasi steady-state slope of n(m) ~ m^{-1.88} [n(a) ~ a^{-3.65}] as a robust solution for a wide range of possible model parameters. We also show that a simple power-law function can be an appropriate approximation for the mass distribution of particles in certain regimes. The steeper solution has observable effects in the submillimeter and millimeter wavelength regimes of the electromagnetic spectrum. We assemble data for nine debris disks that have been observed at these wavelengths and, using a simplified absorption efficiency model, show that the predicted slope of the particle mass distribution generates SEDs that are in agreement with the observed ones.