We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 seconds on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications.

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.

Interferometers, such as the Event Horizon Telescope (EHT), do not directly observe the images of sources but rather measure their Fourier components at discrete spatial frequencies up to a maximum value set by the longest baseline in the array. Construction of images from the Fourier components or analysis of them with high-resolution models requires careful treatment of fine source structure nominally beyond the array resolution. The primary EHT targets, Sgr A* and M87, are expected to have black-hole shadows with sharp edges and strongly filamentary emission from the surrounding plasma on scales much smaller than those probed by the currently largest baselines. We show that for aliasing not to affect images reconstructed with regularized maximum likelihood methods and model images that are directly compared to the data, the sampling of these images (i.e., their pixel spacing) needs to be significantly finer than the scale probed by the largest baseline in the array. Using GRMHD simulations of black-hole images, we estimate the maximum allowable pixel spacing to be approximately equal to (1/8)GM/c^2; for both of the primary EHT targets, this corresponds to an angular pixel size of <0.5 microarcseconds. With aliasing under control, we then advocate use of the second-order Butterworth filter with a cut-off scale equal to the maximum array baseline as optimal for visualizing the reconstructed images. In contrast to the traditional Gaussian filters, this Butterworth filter retains most of the power at the scales probed by the array while suppressing the fine image details for which no data exist.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The Reduced Relativistic Gas (RRG) is a simplified version of the ideal relativistic gas, which assumes that all particles have the same momentum magnitude. Although this is a very idealized situation, the resulting model preserves the phenomenology of Maxwell-Boltzmann distribution and, in some situations, can be described as a perfect fluid, without introducing large errors in both cosmological background and first-order perturbations. The perfect fluid description of RRG model was already used to study the warmness of dark matter, massive neutrinos and interaction of baryons and photons before recombination, showing very good agreement with previous works based on the full Einstein-Boltzmann system of equations. In order to understand these results and construct a more general and formal framework for RRG, we develop a theoretical description of first-order cosmological perturbations of RRG, based on a distribution function which encodes the simplifying assumption that all particles have the same momentum magnitude. The full set of Einstein-Boltzmann equations for RRG distribution are derived and quantities beyond the perfect fluid approximation are studied. Using RRG to describe warm dark matter, we show that, for particles with $m \sim \text{keV}$, the perfect fluid approximation is valid on scales $k < 10\, \text{h}/\text{Mpc}$, for most of the universe evolution. We also determine initial conditions for RRG in the early universe and study the evolution of potential in a toy model of universe composed only by RRG.

An extension of the Starobinsky model is proposed. Besides the usual Starobinsky Lagrangian, a term proportional to the derivative of the scalar curvature, $\nabla_{\mu}R\nabla^{\mu}R$, is considered. The analyzis is done in the Einstein frame with the introduction of a scalar field and a vector field. We show that inflation is attainable in our model, allowing for a graceful exit. We also build the cosmological perturbations and obtain the leading-order curvature power spectrum, scalar and tensor tilts and tensor-to-scalar ratio. The tensor and curvature power spectrums are compared to the most recent observations from BICEP2/Keck collaboration. We verify that the scalar-to-tensor rate $r$ can be expected to be up to three times the values predicted by Starobinsky model.

We study cosmological perturbations in a model of unified dark matter and dark energy with a sharp transition in the late-time universe. The dark sector is described by a dark fluid which evolves from an early stage at redshifts $z > z_C$ when it behaves as cold dark matter (CDM) to a late time dark energy (DE) phase ($z < z_C$) when the equation of state parameter is $w = -1 + \epsilon$, with a constant $\epsilon$ which must be in the range $0 < \epsilon < 2/3$. We show that fluctuations in the dark energy phase suffer from an exponential instability, the mode functions growing both as a function of comoving momentum $k$ and of conformal time $\eta$. In order that this exponential instability does not lead to distortions of the energy density power spectrum on scales for which we have good observational results, the redshift $z_C$ of transition between the two phases is constrained to be so close to zero that the model is unable to explain the supernova data.

Synthesis imaging of the black hole in the center of the Milky Way, Sgr A*, with the Event Horizon Telescope (EHT) rests on the assumption of a stationary image. We explore the limitations of this assumption using high-cadence GRMHD simulations of Sgr A*. We employ analytic models that capture the basic characteristics of the images to understand the origin of the variability in the simulated visibility amplitudes. We find that, in all simulations, the visibility amplitudes for baselines oriented perpendicular to the spin axis of the black hole typically decrease smoothly over baseline lengths that are comparable to those of the EHT. On the other hand, the visibility amplitudes for baselines oriented parallel to the spin axis show significant structure with one or more minima. This suggests that fitting EHT observations with geometric models will lead to reasonably accurate determination of the orientation of the black-hole on the plane of the sky. However, in the disk-dominated models, the locations and depths of the minima in the visibility amplitudes depend primarily on the width and asymmetry of the crescent-like images and are highly variable. In the jet-dominated models, the locations of the minima are determined by the separation of the two image components but their depths depend primarily on the relative brightness of the two components and are also variable. This suggests that using time-independent models to infer additional black-hole parameters, such as the shadow size or the spin magnitude, will be severely affected by the variability of the accretion flow.

The equivalence between theories depending on the derivatives of $R$, i.e. $f\left( R,\nabla R,...,\nabla^{n}R\right) $, and scalar-multi-tensorial theories is verified. The analysis is done in both metric and Palatini formalisms. It is shown that $f\left( R,\nabla R,...,\nabla^{n}R\right) $ theories are equivalent to scalar-multi-tensorial ones resembling Brans-Dicke theories with kinetic terms $\omega_{0}=0$ and $\omega_{0}= - \frac{3}{2}$ for metric and Palatini formalisms respectively. This result is analogous to what happens for $f(R)$ theories. It is worthy emphasizing that the scalar-multi-tensorial theories obtained here differ from Brans-Dicke ones due to the presence of multiple tensorial fields absent in the last. Furthermore, sufficient conditions are established for $f\left( R,\nabla R,...,\nabla^{n}R\right) $ theories to be written as scalar-multi-tensorial theories. Finally, some examples are studied and the comparison of $f\left( R,\nabla R,...,\nabla^{n}R\right) $ theories to $f\left( R,\Box R,...\Box^{n}R\right) $ theories is performed.

Scattering in the ionized interstellar medium is commonly observed to be anisotropic, with theories of magnetohydrodynamic (MHD) turbulence explaining the anisotropy through a preferred magnetic field direction throughout the scattering regions. In particular, the line of sight to the Galactic Center supermassive black hole, Sgr A*, exhibits strong and anisotropic scattering, which dominates its observed size at wavelengths of a few millimeters and longer. Therefore, inferences of the intrinsic structure of \sgra\ at these wavelengths are sensitive to the assumed scattering model. In addition, extrapolations of the scattering model from long wavelengths, at which its parameters are usually estimated, to 1.3 mm, where the Event Horizon Telescope (EHT) seeks to image Sgr A* on Schwarzschild-radius scales, are also sensitive to the assumed scattering model. Past studies of Sgr A* have relied on simple Gaussian models for the scattering kernel that effectively presume an inner scale of turbulence far greater than the diffractive scale; this assumption is likely violated for Sgr A* at 1.3 mm. We develop a physically motivated model for anisotropic scattering, using a simplified model for MHD turbulence with a finite inner scale and a wandering transverse magnetic field direction. We explore several explicit analytic models for this wandering and derive the expected observational properties --- scatter broadening and refractive scintillation --- for each. For expected values of the inner scale, the scattering kernel for all models is markedly non-Gaussian at 1.3 mm but is straightforward to calculate and depends only weakly on the assumed model for the wandering of the magnetic field direction. On the other hand, in all models, the refractive substructure depends strongly on the wandering model and may be an important consideration in imaging Sgr A* with the EHT.

We explore the use of principal component analysis (PCA) to characterize high-fidelity simulations and interferometric observations of the millimeter emission that originates near the horizons of accreting black holes. We show mathematically that the Fourier transforms of eigenimages derived from PCA applied to an ensemble of images in the spatial-domain are identical to the eigenvectors of PCA applied to the ensemble of the Fourier transforms of the images, which suggests that this approach may be applied to modeling the sparse interferometric Fourier-visibilities produced by an array such as the Event Horizon Telescope (EHT). We also show that the simulations in the spatial domain themselves can be compactly represented with a PCA-derived basis of eigenimages allowing for detailed comparisons between variable observations and time-dependent models, as well as for detection of outliers or rare events within a time series of images. Furthermore, we demonstrate that the spectrum of PCA eigenvalues is a diagnostic of the power spectrum of the structure and, hence, of the underlying physical processes in the simulated and observed images.

Fast and accurate integration of geodesics in Kerr spacetimes is an important tool in modeling the orbits of stars and the transport of radiation in the vicinities of black holes. Most existing integration algorithms employ Boyer-Lindquist coordinates, which have coordinate singularities at the event horizon and along the poles. Handling the singularities requires special numerical treatment in these regions, often slows down the calculations, and may lead to inaccurate geodesics. We present here a new general-purpose geodesic integrator, GRay2, that overcomes these issues by employing the Cartesian form of Kerr-Schild coordinates. By performing particular mathematical manipulations of the geodesic equations and several optimizations, we develop an implementation of the Cartesian Kerr-Schild coordinates that outperforms calculations that use the seemingly simpler equations in Boyer-Lindquist coordinates. We also employ the OpenCL framework, which allows GRay2 to run on multi-core CPUs as well as on a wide range of GPU hardware accelerators, making the algorithm more versatile. We report numerous convergence tests and benchmark results for GRay2 for both time-like (particle) and null (photon) geodesics.

We investigate the influence of massive photons on the evolution of the expanding universe. Two particular models for generalized electrodynamics are considered, namely de Broglie-Proca and Bopp-Podolsky electrodynamics. We obtain the equation of state (EOS) $P=P(\varepsilon)$ for each case using dispersion relations derived from both theories. The EOS are inputted into the Friedmann equations of a homogeneous and isotropic space-time to determine the cosmic scale factor $a(t)$. It is shown that the photon non-null mass does not significantly alter the result $a\propto t^{1/2}$ valid for a massless photon gas; this is true either in de Broglie-Proca's case (where the photon mass $m$ is extremely small) or in Bopp-Podolsky theory (for which $m$ is extremely large).

The propagation of non-linear electromagnetic waves is carefully analyzed on a curved spacetime created by static spherically symmetric mass and charge distribution. We compute how non-linear electrodynamics affects the geodesic deviation and the redshift of photons propagating near this massive charged object. In the first order approximation, the effects of electromagnetic self-interaction can be distinguished from the usual Reissner-Nordstr\"om terms. In the particular case of Euler-Heisenberg effective Lagrangian, we find that these self-interaction effects might be important near extremal compact charged objects.

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*.

The observable quantities that carry the most information regarding the structures of the images of black holes in the interferometric observations with the Event Horizon Telescope are the closure phases along different baseline triangles. We use long time span, high cadence, GRMHD+radiative transfer models of Sgr A$^*$ to investigate the expected variability of closure phases in such observations. We find that, in general, closure phases along small baseline triangles show little variability, except in the cases when one of the triangle vertices crosses one of a small regions of low visibility amplitude. The closure phase variability increases with the size of the baseline triangle, as larger baselines probe the small-scale structures of the images, which are highly variable. On average, the jet-dominated MAD models show less closure phase variability than the disk-dominated SANE models, even in the large baseline triangles, because the images from the latter are more sensitive to the turbulence in the accretion flow. Our results suggest that image reconstruction techniques need to explicitly take into account the closure phase variability, especially if the quality and quantity of data allow for a detailed characterization of the nature of variability. This also implies that, if image reconstruction techniques that rely on the assumption of a static image are utilized, regions of the $u-v$ space that show a high level of variability will need to be identified and excised.

We propose a phenomenological unified model for dark matter and dark energy based on an equation of state parameter $w$ that scales with the $\arctan$ of the redshift. The free parameters of the model are three constants: $\Omega_{b0}$, $\alpha$ and $\beta$. Parameter $\alpha$ dictates the transition rate between the matter dominated era and the accelerated expansion period. The ratio $\beta / \alpha$ gives the redshift of the equivalence between both regimes. Cosmological parameters are fixed by observational data from Primordial Nucleosynthesis (PN), Supernovae of the type Ia (SNIa), Gamma-Ray Bursts (GRB) and Baryon Acoustic Oscillations (BAO). The calibration of the 138 GRBs events is performed using the 580 SNIa of the Union2.1 data set and a new set of 79 high-redshift GRBs is obtained. The various sets of data are used in different combinations to constraint the parameters through statistical analysis. The unified model is compared to the $\Lambda$CDM model and their differences are emphasized.

We explore the variability properties of long, high cadence GRMHD simulations across the electromagnetic spectrum using an efficient, GPU-based radiative transfer algorithm. We focus on both disk- and jet-dominated simulations with parameters that successfully reproduce the time-averaged spectral properties of Sgr A* and the size of its image at 1.3mm. We find that the disk-dominated models produce short timescale variability with amplitudes and power spectra that closely resemble those inferred observationally. In contrast, jet-dominated models generate only slow variability, at lower flux levels. Neither set of models show any X-ray flares, which most likely indicate that additional physics, such as particle acceleration mechanisms, need to be incorporated into the GRMHD simulations to account for them. The disk-dominated models show strong, short-lived mm/IR flares, with short (<~ 1hr) time lags between the mm and IR wavelengths, that arise from strong-field gravitational lensing of magnetic flux tubes near the horizon. Such events provide a natural explanation for the observed IR flares with no X-ray counterparts.

This paper analyses the cosmological consequences of a modified theory of gravity whose action integral is built from a linear combination of the Ricci scalar $R$ and a quadratic term in the covariant derivative of $R$. The resulting Friedmann equations are of the fifth-order in the Hubble function. These equations are solved numerically for a flat space section geometry and pressureless matter. The cosmological parameters of the higher-order model are fit using SN Ia data and X-ray gas mass fraction in galaxy clusters. The best-fit present-day $t_{0}$ values for the deceleration parameter, jerk and snap are given. The coupling constant $\beta$\ of the model is not univocally determined by the data fit, but partially constrained by it. Density parameter $\Omega_{m0}$ is also determined and shows weak correlation with the other parameters. The model allows for two possible future scenarios: there may be either a premature Big Rip or a Rebouncing event depending on the set of values in the space of parameters. The analysis towards the past performed with the best-fit parameters shows that the model is not able to accommodate a matter-dominated stage required to the formation of structure.

We have detected bright HC7N J = 21-20 emission toward multiple locations in the Serpens South cluster-forming region using the K-Band Focal Plane Array at the Robert C. Byrd Green Bank Telescope. HC7N is seen primarily toward cold filamentary structures that have yet to form stars, largely avoiding the dense gas associated with small protostellar groups and the main central cluster of Serpens South. Where detected, the HC7N abundances are similar to those found in other nearby star forming regions. Toward some HC7N `clumps', we find consistent variations in the line centroids relative to NH3 (1,1) emission, as well as systematic increases in the HC7N non-thermal line widths, which we argue reveal infall motions onto dense filaments within Serpens South with minimum mass accretion rates of M ~ 2-5 M_sun Myr^-1. The relative abundance of NH3 to HC7N suggests that the HC7N is tracing gas that has been at densities n ~ 10^4 cm^-3, for timescales t < 1-2 x 10^5 yr. Since HC7N emission peaks are rarely co-located with those of either NH3 or continuum, it is likely that Serpens South is not particularly remarkable in its abundance of HC7N, but instead the serendipitous mapping of HC7N simultaneously with NH3 has allowed us to detect HC7N at low abundances in regions where it otherwise may not have been looked for. This result extends the known star-forming regions containing significant HC7N emission from typically quiescent regions, like the Taurus molecular cloud, to more complex, active environments.

This work presents a complete cyclic cosmological scenario based on nonlinear magnetic field. It is constructed a model composed by five fluids namely baryonic matter, dark matter, radiation, neutrinos and a cosmological magnetic field. The first four fluids are treated in the standard way and the fifth fluid, the magnetic field, is described by a nonlinear electrodynamics. The free parameters are fitted by observational data (SNIa, CMB, extragalactic magnetic fields, etc) and by simple theoretical considerations. As result arises a cyclic cosmological model which preserves the main successes of standard big bang model and solve some other problems like the initial singularity, the present acceleration and the Big Rip.

We construct a phenomenological theory of gravitation based on a second order gauge formulation for the Lorentz group. The model presents a long-range modification for the gravitational field leading to a cosmological model provided with an accelerated expansion at recent times. We estimate the model parameters using observational data and verify that our estimative for the age of the Universe is of the same magnitude than the one predicted by the standard model. The transition from the decelerated expansion regime to the accelerated one occurs recently (at $\sim9.3\;Gyr$).

Early universe equations of state including realistic interactions between constituents are built up. Under certain reasonable assumptions, these equations are able to generate an inflationary regime prior to the nucleosynthesis period. The resulting accelerated expansion is intense enough to solve the flatness and horizon problems. In the cases of curvature parameter \kappa equal to 0 or +1, the model is able to avoid the initial singularity and offers a natural explanation for why the universe is in expansion.

This work is an application of the second order gauge theory for the Lorentz group, where a description of the gravitational interaction is obtained which includes derivatives of the curvature. We analyze the form of the second field strenght, $G=\partial F +fAF$, in terms of geometrical variables. All possible independent Lagrangians constructed with quadratic contractions of $F$ and quadratic contractions of $G$ are analyzed. The equations of motion for a particular Lagrangian, which is analogous to Podolsky's term of his Generalized Electrodynamics, are calculated. The static isotropic solution in the linear approximation was found, exhibiting the regular Newtonian behaviour at short distances as well as a meso-large distance modification.

The high precision attained by cosmological data in the last few years has increased the interest in exact solutions. Analytic expressions for solutions in the Standard Model are presented here for all combinations of $\Lambda = 0$, $\Lambda \ne 0$, $\kappa = 0$ and $\kappa \ne 0$, in the presence and absence of radiation and nonrelativistic matter. The most complete case (here called the $\Lambda \gamma CDM$ Model) has $\Lambda \ne 0, \kappa \ne 0$, and supposes the presence of radiation and dust. It exhibits clearly the recent onset of acceleration. The treatment includes particular models of interest such as the $\Lambda$CDM Model (which includes the cosmological constant plus cold dark matter as source constituents).