Source author record

A. Mangalam

A. Mangalam 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

20works
6topics
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

20 published item(s)

preprint2020arXiv

A relativistic orbit model for temporal properties of AGN

We present a unified model for X-ray quasi-periodic oscillations (QPOs) seen in Narrow-line Seyfert 1 (NLSy1) galaxies, $γ$-ray and optical band QPOs that are seen in Blazars. The origin of these QPOs is attributed to the plasma motion in corona or jets of these AGN. In the case of X-ray QPOs, we applied the general relativistic precession model for the two simultaneous QPOs seen in NLSy1 1H 0707-945 and deduce orbital parameters, such the radius of the emission region, and spin parameter $a$ for a circular orbit; we obtained the Carter's constant $Q$, $a$, and the radius in the case of a spherical orbit solution. In other cases where only one X-ray QPO is seen, we localized the orbital parameters for NLSy1 galaxies REJ 1034+396, 2XMM J123103.2+110648, MS 2254.9-3712, Mrk 766, and MCG-06-30-15. By applying the lighthouse model, we found that a kinematic origin of the jet based $γ$-ray and optical QPOs, in a relativistic MHD framework, is possible. Based on the inbuilt Hamiltonian formulation with a power-law distribution in the orbital energy of the plasma consisting of only circular or spherical trajectories, we show that the resulting Fourier power spectral density (PSD) has a break corresponding to the energy at ISCO. Further, we derive connection formulae between the slopes in the PSD and that of the energy distribution. Overall, given the preliminary but promising results of these relativistic orbit models to match the observed QPO frequencies and PSD at diverse scales in the inner corona and the jet, it motivates us to build detailed models, including a transfer function for the energy spectrum in the corona and relativistic MHD jet models for plasma flow and its polarization properties.

preprint2020arXiv

Cosmic spin and mass evolution of black holes and its impact

We build an evolution model of the central black hole that depends on the processes of gas accretion, the capture of stars, mergers as well as electromagnetic torque. In case of gas accretion in the presence of cooling sources, the flow is momentum-driven, after which the black hole reaches a saturated mass; subsequently, it grows only by stellar capture and mergers. We model the evolution of the mass and spin with the initial seed mass and spin in $Λ$CDM cosmology. For stellar capture, we have assumed a power-law density profile for the stellar cusp in a framework of relativistic loss cone theory that include the effects of black hole spin, Carter's constant, loss cone angular momentum, and capture radius. Based on this, the predicted capture rates of $10^{-5}$--$10^{-6}$ yr$^{-1}$ are closer to the observed range. We have considered the merger activity to be effective for $z \lesssim 4$, and we self-consistently include the Blandford-Znajek torque. We calculate these effects on the black hole growth individually and in combination, for deriving the evolution. Before saturation, accretion dominates the black hole growth ($\sim 95\%$ of the final mass), and subsequently, stellar capture and mergers take over with roughly equal contribution. The simulations of the evolution of the $M_{\bullet} - σ$ relation using these effects are consistent with available observations. We run our model backward in time and retrodict the parameters at formation. Our model will provide useful inputs for building demographics of the black holes and in formation scenarios involving stellar capture.

preprint2020arXiv

Dynamics of accretion and winds in tidal disruption events

We have constructed self-similar models of a time-dependent accretion disk in both sub and super-Eddington phases with wind outflows for tidal disruption events (TDEs). The physical input parameters are the black hole (BH) mass $M_{\bullet}$, specific orbital energy $E$ and angular momentum $J$, star mass $M_{\star}$ and radius $R_{\star}$. We consider the sub-Eddington phase to be total pressure (model A1) and gas pressure (model A2) dominated. In contrast, the super-Eddington phase is dominated by radiation pressure (model B) with Thomson opacity. We derive the viscosity prescribed by the stress tensor, $Π_{rϕ}\propto Σ_d^b r^d$ where $Σ_d$ is the surface density of the disk, $r$ is the radius and $b$ and $d$ are constants. The specific choice of radiative or $α$ viscosity is motivated, and its parameters are decided by the expected disk luminosity and evolution time scale being in the observed range. The disk evolves due to mass loss by accretion onto the black hole and outflowing wind, and mass gain by fallback of the debris; this results in an increasing outer radius. We have simulated the luminosity profile for both sub and super-Eddington disks. As an illustrative example, we fit our models to the observations in X-ray, UV, and Optical of four TDE events and deduce the physical parameters above.

preprint2018arXiv

Black hole demographics from TDE modeling

The occurrence rate of tidal disruption events (TDEs) by survey missions depend on the black hole mass function of the galaxies, properties of the stellar cusp and mass of the central black hole. Using a power law density profile with Kroupa mass function, we solve the steady state Fokker-Planck to calculate the theoretical capture rate of stars by the black hole. Using a steady accretion model, the Schechter black hole mass function (BHMF) and the cosmological parameters, we calculate the detection rate of TDEs for various surveys which is then fit with the observed TDE rates to extract the Schechter parameters. The rate tension between observation ($\sim 10^{-5}~{\rm yr^{-1}}$) and theory ($\sim 10^{-4}~{\rm yr^{-1}}$) for individual galaxies is explained by the statistical average over the BHMF.

preprint2018arXiv

Evolution of the $M_{\bullet} - σ$ relation

Black holes at the centers of the galaxies grow mainly by the processes of accretion, mergers, and consumption of stars. In the case of gas accretion with cooling sources, the flow is momentum driven, after which the black hole reaches a saturated mass, and subsequently, it grows only by consumption of stars. In addition, we include the effect of mergers on the growth of black hole spin and mass and study its evolution as a function of redshift in a $Λ$CDM cosmology using an initial seed mass and spin distribution functions that we have derived. For the stellar ingestion, we have assumed a power-law density profile for the galaxy in our framework of a new relativistic loss cone theory that includes the effect of the black hole spin. We predict the impact of the evolution on the $M_{\bullet} - σ$ relation and compare it with available observations.

preprint2017arXiv

The $M_{\bullet} - σ$ relation in spherical systems

To investigate the $M_\bullet -σ$ relation, we consider realistic elliptical galaxy profiles that are taken to follow a single power law density profile given by $ρ(r) = ρ_{0}(r/ r_{0})^{-γ}$ or the Nuker intensity profile. We calculate the density using Abel's formula in the latter case by employing the derived stellar potential in both cases, we derive the distribution function $f(E)$ of the stars in presence of the supermassive black hole (SMBH) at the center and hence compute the line of sight (LOS) velocity dispersion as a function of radius. For the typical range of values for masses of SMBH, we obtain $M_{\bullet} \propto σ^{p}$ for different profiles. An analytical relation $p = (2γ+ 6)/(2 + γ)$ is found which is in reasonable agreement with observations (for $γ$ = 0.75 - 1.4, $p$ = 3.6 - 5.3). Assuming that a proportionality relation holds between the black hole mass and bulge mass, $ M_{\bullet} =f M_b$, and applying this to several galaxies we find the individual best fit values of $p$ as a function of $f$; also by minimizing $χ^{2}$, we find the best fit global $p$ and $f$. For Nuker profiles we find that $p$ = $3.81 \pm 0.004$ and $f$ = $(1.23 \pm 0.09)\times 10^{-3}$ which are consistent with the observed ranges.

preprint2016arXiv

Parsec-scale jet properties of the quasar PG 1302$-$102

The quasar PG 1302$-$102 is believed to harbour a supermassive binary black hole (SMBBH) system. Using the available 15 GHz and $2-8$ GHz, multi-epoch Very Long Baseline Array data, we constrain the pc-scale jet properties based on the inferred mean proper motion, including a bulk Lorentz factor $\geq 5.1 \pm 0.8$, jet inclination angle $\leq (11.4 \pm 1.7)$ degrees, projected position angle $= 31.8$ degrees, intrinsic half opening angle $\leq (0.9 \pm 0.1)$ degrees and a mean $2-8$ GHz spectral index of 0.31. A general relativistic helical jet model is presented and applied to predict quasi-periodic oscillations of $\sim$ 10 days, power law power spectrum shape and a contribution of up to $\sim$ 53 percent to the observed variable core flux density. The model is used to make a case for high resolution, moderately sampled, long duration radio interferometric observations to reveal signatures due to helical knots and distinguish them from those due to SMBBH orbital activity including a phase difference $\sim π$ and an amplitude ratio (helical light curve amplitude/SMBBH light curve amplitude) of $0.2-3.3$. The prescription can be used to identify helical kinematic signatures from quasars, providing possible candidates for further studies with polarization measurements. It can also be used to infer promising SMBBH candidates for the study of gravitational waves if there are systematic deviations from helical signatures.

preprint2015arXiv

A global galactic dynamo with a corona constrained by relative helicity

We present a model for a global axisymmetric turbulent dynamo operating in a galaxy with a corona which treats the supernovae (SNe) and magneto-rotational instability (MRI) driven turbulence parameters under a common formalism. The nonlinear quenching of the dynamo is alleviated by inclusion of small-scale advective and diffusive magnetic helicity fluxes, which allow the gauge invariant magnetic helicity to be transferred outside the disk and consequently build up a corona during the course of dynamo action. The time-dependent dynamo equations are expressed in a separable form and solved through an eigenvector expansion constructed using the steady-state solutions of the dynamo equation. The parametric evolution of the dynamo solution allows us to estimate the final structure of the global magnetic field and the saturated value of the turbulence parameter $α_m$, even before solving the dynamical equations for evolution of magnetic fields in the disk and the corona, along with $α$-quenching. We then solve these equations simultaneously to study the saturation of large-scale magnetic field, its dependence on the small-scale magnetic helicity fluxes and corresponding evolution of the force-free field in the corona. The quadrupolar large-scale magnetic field in the disk is found to reach equipartition strength within a timescale of 1 Gyr. The large-scale magnetic field in the corona obtained is much weaker in strength compared to the field inside the disk and has only a weak impact on the dynamo operation.

preprint2015arXiv

Frequency dependent core shifts and parameter estimation for the blazar 3C 454.3

We study the core shift effect in the parsec scale jet of the blazar 3C 454.3 using the 4.8 GHz - 36.8 GHz radio light curves obtained from three decades of continuous monitoring. From a piecewise Gaussian fit to each flare, time lags $Δt$ between the observation frequencies $ν$ and spectral indices $α$ based on peak amplitudes $A$ are determined. From the fit $Δt \propto ν^{1/k_r}$, $k_r = 1.10 \pm 0.18$ indicating equipartition between the magnetic field energy density and the particle energy density. From the fit $A \propto ν^α$, $α$ is in the range $-0.24$ to $1.52$. A mean magnetic field strength at 1 pc, $B_1 = 0.5 \pm 0.2$ G, and at the core, $B_{\rm core} = 46 \pm 16$ mG, are inferred, consistent with previous estimates. The measure of core position offset is $Ω_{rν} = 6.4 \pm 2.8$ pc GHz$^{1/k_r}$ when averaged over all frequency pairs. Based on the statistical trend shown by the measured core radius $r_{\rm core}$ as a function of $ν$, we infer that the synchrotron opacity model may not be valid for all cases. A Fourier periodogram analysis yields power law slopes in the range $-1.6$ to $-3.5$ describing the power spectral density shape and gives bend timescales in the range $0.52 - 0.66~$yr. This result, and both positive and negative $α$, indicate that the flares originate from multiple shocks in a small region. Important objectives met in our study include: the demonstration of the computational efficiency and statistical basis of the piecewise Gaussian fit; consistency with previously reported results; evidence for the core shift dependence on observation frequency and its utility in jet diagnostics in the region close to the resolving limit of very long baseline interferometry observations.

preprint2015arXiv

Kinematics of and emission from helically orbiting blobs in a relativistic magnetized jet

We present a general relativistic (GR) model of jet variability in active galactic nuclei due to orbiting blobs in helical motion along a funnel or cone shaped magnetic surface anchored to the accretion disk near the black hole. Considering a radiation pressure driven flow in the inner region, we find that it stabilizes the flow, yielding Lorentz factors ranging between 1.1 and 7 at small radii for reasonable initial conditions. Assuming these as inputs, simulated light curves (LCs) for the funnel model include Doppler and gravitational shifts, aberration, light bending, and time delay. These LCs are studied for quasi-periodic oscillations (QPOs) and the power spectral density (PSD) shape and yield an increased amplitude ($\sim$ 12 %); a beamed portion and a systematic phase shift with respect to that from a previous special relativistic model. The results strongly justify implementing a realistic magnetic surface geometry in Schwarzschild geometry to describe effects on emission from orbital features in the jet close to the horizon radius. A power law shaped PSD with a typical slope of $-2$ and QPOs with timescales in the range of $(1.37 - 130.7)$ days consistent with optical variability in Blazars, emerges from the simulations for black hole masses $M_{\bullet} = (0.5 - 5) \times 10^8 M_{\odot}$ and initial Lorentz factors $γ_{jet,i} = 2 - 10$. The models presented here can be applied to explain radio, optical, and X-ray variability from a range of jetted sources including active galactic nuclei, X-ray binaries and neutron stars.

preprint2015arXiv

Stellar and gas dynamical model for tidal disruption events in a quiescent galaxy

A detailed model of the tidal disruption events (TDEs) has been constructed using stellar dynamical and gas dynamical inputs that include black hole (BH) mass $M_{\bullet}$, specific orbital energy $E$ and angular momentum $J$, star mass $M_{\star}$ and radius $R_{\star}$, and the pericenter of the star orbit $r_{p}(E,\hspace{1mm}J,\hspace{1mm}M_{\bullet})$. We solved the steady state Fokker--Planck equation using the standard loss cone theory for the galactic density profile $ρ(r) \propto r^{-γ}$ and stellar mass function $ξ(m) $ where $m=M_{\star}/M_{\odot}$ and obtained the feeding rate of stars to the BH integrated over the phase space as $\dot{N}_{t} \propto M_{\bullet}^β$, where $β= -0.3\pm 0.01$ for $M_{\bullet}>10^7 M_{\odot}$ and $\sim 6.8 \hspace{1mm} \times 10^{-5}$ Yr$^{-1}$ for $γ=0.7$. We use this to model the in-fall rate of the disrupted debris, $\dot{M}(E,\hspace{1mm}J,\hspace{1mm}m,\hspace{1mm}t)$, and discuss the conditions for the disk formation, finding that the accretion disk is almost always formed for the fiduciary range of the physical parameters. We also find the conditions under which the disk formed from the tidal debris of a given star with a super Eddington accretion phase. We have simulated the light curve profiles in the relevant optical g band and soft X-rays for both super and sub-Eddington accretion disks as a function of $\dot{M}(E,\hspace{1mm}J,\hspace{1mm}t)$. Using this, standard cosmological parameters, and mission instrument details, we predict the detectable TDE rates for various forthcoming surveys finally as a function of $γ$.

preprint2014arXiv

Models of force-free spheres and applications to solar active regions

Here we present a systematic study of force-free field equation for simple axisymmetric configurations in spherical geometry. The condition of separability of solutions in radial and angular variables leads to two classes of solutions: linear and non-linear force-free fields. We have studied these linear solutions Chandrasekhar (1956) and extended the non-linear solutions given in Low \& Lou (1990) to the irreducible rational form $n= p/q$, which is allowed for all cases of odd $p$ and to cases of $q>p$ for even $p$. We have further calculated their energies and relative helicities for magnetic field configurations in finite and infinite shell geometries. We demonstrate here a method here to be used to fit observed magnetograms as well as to provide good exact input fields for testing other numerical codes used in reconstruction on the non-linear force-free fields.

preprint2014arXiv

Orbital signatures from observed light curves of blazars

Variability in active galactic nuclei is observed in ultraviolet to X-ray emission based light curves. This could be attributed to orbital signatures of the plasma that constitutes the accretion flow on the putative disk or in the developing jet close to the inner region of the central black hole. We discuss some theoretical models which build on this view. These models include general relativistic effects such as light bending, aberration effects, gravitational and Doppler redshifts. The novel aspects relate to the treatment of helical flow in cylindrical and conical geometries in the vicinity of a Schwarzschild black hole that leads to amplitude and frequency modulations of simulated light curves as well as the inclusion of beaming effects in these idealized geometries. We then present a suite of time series analysis techniques applicable to data with varied properties which can extract detailed information from them for their use in theoretical models.

preprint2014arXiv

Parametric models of the periodogram

The maximum likelihood estimator is used to determine fit parameters for various parametric models of the Fourier periodogram followed by the selection of the best fit model amongst competing models using the Akaike information criteria. This analysis, when applied to light curves of active galactic nuclei can be used to infer the presence of quasi-periodicity and break or knee frequencies. The extracted information can be used to place constraints on the mass, spin and other properties of the putative central black hole and the region surrounding it through theoretical models involving disk and jet physics.

preprint2014arXiv

Separable solutions of force-free spheres and applications to solar active regions

In this paper, we present a systematic study of the force-free field equation for simple axisymmetric configurations in spherical geometry and apply it to the solar active regions. The condition of separability of solutions in the radial and angular variables leads to two classes of solutions: linear and nonlinear force-free fields. We have studied these linear solutions and extended the nonlinear solutions for the radial power law index to the irreducible rational form $n= p/q$, which is allowed for all cases of odd $p$ and cases of $q>p$ for even $p$, where the poloidal flux $ψ\propto1/r^n$ and field $\mathbf{B}\propto 1/r^{n+2}$). We apply these solutions to simulate photospheric vector magnetograms obtained using the spectropolarimeter on board Hinode. The effectiveness of our search strategy is first demonstrated on test inputs of dipolar, axisymmetric,and non axisymmetric linear force-free fields. Using the best-fit to these magnetograms, we build three-dimensional axisymmetric field configurations and calculate the energy and relative helicity with two independent methods, which are in agreement. We have analyzed five magnetograms for AR 10930 spanning a period of three days during which two X-class flares occurred which allowed us to find the free energy and relative helicity of the active region before and after the flare; our analysis indicates a peak in these quantities before the flare events which is consistent with the results mentioned in literature. We also analyzed single-polarity regions AR 10923 and 10933, which showed very good fits with potential fields. This method can provide useful reconstruction of the nonlinear force-free (NLFF) fields as well as reasonably good input fields for other numerical techniques.

preprint2014arXiv

X-Ray Variability and the Inner Region in Active Galactic Nuclei

We present theoretical models of X-ray variability attributable to orbital signatures from an accretion disk including emission region size, quasi-periodic oscillations (QPOs) and its quality factor $Q$, and the emergence of a break frequency in the power spectral density shape. We find a fractional variability amplitude of $F_{var}\propto M^{-0.4}_{\bullet}$. We conduct a time series analysis on X-ray light curves ($0.3-10$ keV) of a sample of active galactic nuclei (AGNs). A statistically significant bend frequency is inferred in 9 of 58 light curves (16 %) from 3 AGNs for which the break timescale is consistent with the reported BH spin but not with the reported BH mass. Upper limits of $2.85 \times 10^7 M_\odot$ in NGC 4051, $8.02 \times 10^7 M_\odot$ in MRK 766 and $4.68 \times 10^7 M_\odot$ in MCG-6-30-15 are inferred for maximally spinning BHs. For REJ 1034+396, where a QPO at 3733 s was reported, we obtain an emission region size of $(6 - 6.5) M$ and a BH spin $a\lesssim$ 0.08. The relativistic inner region of a thin disk, dominated by radiation pressure and electron scattering is likely to host the orbital features as the simulated $Q$ ranges from $6.3 \times 10^{-2}$ to $4.25 \times 10^6$, containing the observed $Q$. The derived value of $Q \sim$ 32 for REJ 1034+396 therefore suggests that the AGN hosts a thin disk.

preprint2012arXiv

Multiwavelength Intraday Variability of the BL Lac S5 0716+714

We report results from a 1 week multi-wavelength campaign to monitor the BL Lac object S5 0716+714 (on December 9-16, 2009). In the radio bands the source shows rapid (~ (0.5-1.5) day) intra-day variability with peak amplitudes of up to ~ 10 %. The variability at 2.8 cm leads by about 1 day the variability at 6 cm and 11 cm. This time lag and more rapid variations suggests an intrinsic contribution to the source's intraday variability at 2.8 cm, while at 6 cm and 11 cm interstellar scintillation (ISS) seems to predominate. Large and quasi-sinusoidal variations of ~ 0.8 mag were detected in the V, R and I-bands. The X-ray data (0.2-10 keV) do not reveal significant variability on a 4 day time scale, favoring reprocessed inverse-Compton over synchrotron radiation in this band. The characteristic variability time scales in radio and optical bands are similar. A quasi-periodic variation (QPO) of 0.9 - 1.1 days in the optical data may be present, but if so it is marginal and limited to 2.2 cycles. Cross-correlations between radio and optical are discussed. The lack of a strong radio-optical correlation indicates different physical causes of variability (ISS at long radio wavelengths, source intrinsic origin in the optical), and is consistent with a high jet opacity and a compact synchrotron component peaking at ~= 100 GHz in an ongoing very prominent flux density outburst. For the campaign period, we construct a quasi-simultaneous spectral energy distribution (SED), including gamma-ray data from the FERMI satellite. We obtain lower limits for the relativistic Doppler-boosting of delta >= 12-26, which for a BL\,Lac type object, is remarkably high.

preprint2011arXiv

On the injection spectrum of relativistic electrons in high-redshift radio galaxies

We point out that the remarkable linearity of the ultra-steep radio spectra of high redshift radio galaxies reflects a previously reported general trend for powerful radio galaxies, according to which the spectral curvature is lesser for sources having steeper spectra (measured near rest-frame 1 GHz). We argue based on existing theoretical and observational evidence that it is premature to conclude that the particle acceleration mechanism in sources having straight, ultra-steep radio spectra gives rise to an ultra-steep injection spectrum of the radiating electrons. In empirical support to this we show that the estimated injection spectral indices, available for a representative sample of 35 compact steep spectrum (CSS) radio sources are not correlated with their rest-frame (intrinsic) rotation measures, which are known to be typically large, indicating a dense environment, as is also the case for high-$z$ radio galaxies.

preprint2009arXiv

The changing ISM of massive elliptical galaxies and cosmic evolution of radio galaxies and quasars

The recently discovered apparent dramatic expansion in the effective radii of massive elliptical galaxies from $z \simeq 2$ to $z \simeq 0.1$ has been interpreted in terms of either galaxy mergers or the rapid loss of cold gas due to AGN feedback. In examining the latter case we have quantified the extent of the expansion, which is uncertain observationally, in terms of the star formation parameters and time of the expulsion of the cold gas. In either case, the large global decrease in stellar density should translate into a major drop in the ISM density and pressure, and a much steeper radial decline in those quantities with cosmic epoch. These cosmological changes are expected to have a major influence on the gas accretion mode, which will shift from `cold' thin disk accretion at high redshifts toward `hot' Bondi fed ADAF accretion at low redshifts. The decline of angular momentum inflow would then lead to a spin down of the black hole, for which we have calculated more precise time scales; a value of about 0.2 Gyr is typical for a $10^9 M_\odot$ central black hole. These results have implications for the different cosmological evolutionary patterns found for the luminosity functions of powerful and weak radio galaxies.

preprint2008arXiv

Stellar Disruption and the Quasar Radio Dichotomy

The origin of the dichotomy of radio loudness among quasars can be explained using recent findings that the mass of the central supermassive black hole (SMBH) in extended radio-loud quasars is systematically a few times that of their counterparts in radio-quiet quasars. This sensitive dependence of radio jet ejection upon SMBH mass probably arises from the blockage of jets by the presence of substantial quantities of gas tidally stripped from stars by the central BH. This disruptive gas, however, will only be available around BHs with masses less than $M_c ~\gtrsim ~10^8\Msun$, for which the tidal disruption radius lies outside the SMBH's event horizon. Consequently, we find that AGN with $M_{BH} > M_c$ can successfully launch jets with a wide range of powers, thus producing radio-loud quasars. The great majority of jets launched by less massive BHs, however, will be truncated in the vicinity of the SMBH due to mass loading from this stellar debris. This scenario also can naturally explain the remarkable dearth of extended radio structures in quasars showing broad absorption line spectra.