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 $\Lambda$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} - \sigma$ 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.

We construct two classes of magnetohydrostatic (MHS) equilibria for an axisymmetric vertical flux tube spanning from the photosphere to the lower part of the transition region within a realistic stratified solar atmosphere subject to solar gravity. We assume a general quadratic expression of the magnetic flux function for the gas pressure and poloidal current and solve the Grad-Shafranov equation analytically. The solution is a combination of a homogeneous and a particular part where the former is separable by a Coulomb function in $r$ and exponential in $z$, while the particular part is an open configuration that has no $z$ dependence. We also present another open field solution by using a self-similar formulation with two different profile functions and incorporating stratified solar gravity to maintain the magnetohydrostatic equilibria, which is a modification of earlier self-similar models with a twist. We study the admitted parameter space that is consistent with the conditions in the solar atmosphere and derive magnetic and the thermodynamic structures inside the flux tube that are reasonably consistent with the photospheric magnetic bright points (MBPs) for both open and closed field Coulomb function and self-similar models as estimated from observations and simulations. The obtained open and closed field flux tube solutions can be used as the background conditions for the numerical simulations for the study of the wave propagation through the flux tubes. The solutions can also be used to construct realistic magnetic canopies.

We derive alternate and new closed-form analytic solutions for the non-equatorial eccentric bound trajectories, $\{ \phi \left( r, \theta \right)$, $\ t \left( r, \theta \right),\ r \left( \theta \right) \}$, around a Kerr black hole by using the transformation $1/r=\mu \left(1+ e \cos \chi \right)$. The application of the solutions is straightforward and numerically fast. We obtain and implement translation relations between energy and angular momentum of the particle, ($E$, $L$), and eccentricity and inverse-latus rectum, ($e$, $\mu$), for a given spin, $a$, and Carter's constant, $Q$, to write the trajectory completely in the ($e$, $\mu$, $a$, $Q$) parameter space. The bound orbit conditions are obtained and implemented to select the allowed combination of parameters ($e$, $\mu$, $a$, $Q$). We also derive specialized formulae for spherical and separatrix orbits. A study of the non-equatorial analog of the previously studied equatorial separatrix orbits is carried out where a homoclinic orbit asymptotes to an energetically bound spherical orbit. Such orbits simultaneously represent an eccentric orbit and an unstable spherical orbit, both of which share the same $E$ and $L$ values. We present exact expressions for $e$ and $\mu$ as functions of the radius of the corresponding unstable spherical orbit, $r_s$, $a$, and $Q$, and their trajectories, for ($Q\neq0$) separatrix orbits; they are shown to reduce to the equatorial case. These formulae have applications to study the gravitational waveforms from EMRIs besides relativistic precession and phase space explorations. We obtain closed-form expressions of the fundamental frequencies of non-equatorial eccentric trajectories that are equivalent to the previously obtained quadrature forms and also numerically match with the equivalent formulae previously derived. We sketch several orbits and discuss their astrophysical applications.

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 $\Lambda$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} - \sigma$ relation and compare it with available observations.

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.

In this paper, we explore the possibility of accreting primordial black holes as the source of heating for the collapsing gas in the context of the direct collapse black hole scenario for the formation of super-massive black holes (SMBHs) at high redshifts, $z\sim 6-7$. One of the essential requirements for the direct collapse model to work is to maintain the temperature of the in falling gas at $\approx 10^4$ K. We show that even under the existing abundance limits, the primordial black holes of masses $\gtrsim 10^{-2} \ {\rm M}_\odot$, can heat the collapsing gas to an extent that the ${\rm H}_2$ formation is inhibited. The collapsing gas can maintain its temperature at $10^4$ K till the gas reaches a critical density $n_{crit} \approx 10^3$ cm$^{-3}$, at which the roto-vibrational states of ${\rm H}_2$ approaches local thermodynamic equilibrium and ${\rm H}_2$ cooling becomes inefficient. In the absence of ${\rm H}_2$ cooling the temperature of the collapsing gas stays at $\approx 10^4$ K even as it collapses further. We discuss scenarios of subsequent angular momentum removal and the route to find collapse through either a supermassive star or a supermassive disk.

Observational evidence of the braiding of magnetic field lines has been reported. The magnetic reconnection within the loop (nanoflares) and with other loops (microflares) disentangle the field. The coronal field then reorganizes itself to attain a force-free field configuration. We have evaluated the power law index of the energy distribution $f(E)=f_0 E^{-\alpha}$ by using a model of relaxation incorporating different profile functions of winding number distribution $f(w)$ based on braided topologies. We study the radio signatures that occur in the solar corona using the radio data obtained from the Gauribidanur Radio Observatory (IIA) and extract the power law index by using the Statistic-sensitive nonlinear iterative peak clipping (SNIP) algorithm. We see that the power law index obtained from the model is in good agreement with the calculated value from the radio data observation.

We construct two classes of the magnetohydrostatic equilibria of the axisymmetric flux tubes with twisted magnetic fields in the stratified solar atmosphere that span from the photosphere to the transition region. We built the models by incorporating specific forms of the gas pressure and poloidal current in the Grad-Shafranov equation. This model gives both closed and open field structure of the flux tube. The other open field model we construct is based on the self-similar formulation, where we have incorporated specific forms of the gas pressure, poloidal current and two different shape functions. We study the homology of the parameter space that is consistent with the solar atmosphere and find that the estimation of the magnetic structure inside the flux tubes is consistent with the observation and simulation results of the magnetic bright points.

We use our semi-analytic solution of the nonlinear force-free field equation to construct three-dimensional magnetic fields that are applicable to the solar corona and study their statistical properties for estimating the degree of braiding exhibited by these fields. We present a new formula for calculating the winding number and compare it with the formula for the crossing number. The comparison is shown for a toy model of two helices and for realistic cases of nonlinear force-free fields; conceptually the formulae are nearly the same but the resulting distributions calculated for a given topology can be different. We also calculate linkages, which are useful topological quantities that are independent measures of the contribution of magnetic braiding to the total free energy and relative helicity of the field. Finally, we derive new analytical bounds for the free energy and relative helicity for the field configurations in terms of the linking number. These bounds will be of utility in estimating the braided energy available for nano-flares or for eruptions.

We describe a model of tidal disruption events (TDEs) with input physical parameters that include the black hole (BH) mass $M_{\bullet}$, the specific orbital energy $E$, the angular momentum $J$, the star mass $M_{\star}$ and radius $R_{\star}$. We calculate the rise time of the TDEs, the peak bolometric luminosity in terms of these physical parameters and a typical light curve of TDEs for various All Sky Survey (ASS) and Deep Sky Survey (DSS) missions. We then derive the expected detection rates and discuss the follow up of TDEs through observations in various spectral bands from X-rays to radio wavelengths.

To investigate the $M_\bullet -\sigma$ relation, we consider realistic elliptical galaxy profiles that are taken to follow a single power law density profile given by $\rho(r) = \rho_{0}(r/ r_{0})^{-\gamma}$ 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 \sigma^{p}$ for different profiles. An analytical relation $p = (2\gamma + 6)/(2 + \gamma)$ is found which is in reasonable agreement with observations (for $\gamma$ = 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 $\chi^{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.

The blazar 3C454.3 exhibited a strong flare seen in gamma-rays, X-rays, and optical/NIR bands during 3--12 December 2009. Emission in the V and J bands rose more gradually than did the gamma-rays and soft X-rays, though all peaked at nearly the same time. Optical polarization measurements showed dramatic changes during the flare, with a strong anti-correlation between optical flux and degree of polarization (which rose from ~ 3% to ~ 20%) during the declining phase of the flare. The flare was accompanied by large rapid swings in polarization angle of ~ 170 degree. This combination of behaviors appear to be unique. We have cm-band radio data during the same period but they show no correlation with variations at higher frequencies. Such peculiar behavior may be explained using jet models incorporating fully relativistic effects with a dominant source region moving along a helical path or by a shock-in-jet model incorporating three-dimensional radiation transfer if there is a dominant helical magnetic field. We find that spectral energy distributions at different times during the flare can be fit using modified one-zone models where only the magnetic field strength and particle break frequencies and normalizations need change. An optical spectrum taken at nearly the same time provides an estimate for the central black hole mass of ~ 2.3 * 10^9 M_sun. We also consider two weaker flares seen during the $\sim 200$ d span over which multi-band data are available. In one of them, the V and J bands appear to lead the $\gamma$-ray and X-ray bands by a few days; in the other, all variations are simultaneous.

We build a single vertical straight magnetic fluxtube spanning the solar photosphere and the transition region which does not expand with height. We assume that the fluxtube containing twisted magnetic fields is in magnetohydrostatic equilibrium within a realistic stratified atmosphere subject to solar gravity. Incorporating specific forms of current density and gas pressure in the Grad--Shafranov equation, we solve the magnetic flux function, and find it to be separable with a Coulomb wave function in radial direction while the vertical part of the solution decreases exponentially. We employ improved fluxtube boundary conditions and take a realistic ambient external pressure for the photosphere to transition region, to derive a family of solutions for reasonable values of the fluxtube radius and magnetic field strength at the base of the axis that are the free parameters in our model. We find that our model estimates are consistent with the magnetic field strength and the radii of Magnetic bright points (MBPs) as estimated from observations. We also derive thermodynamic quantities inside the fluxtube.

We studied the pc-scale core shift effect using radio light curves for three blazars, S5 0716+714, 3C 279 and BL Lacertae, which were monitored at five frequencies ($\nu$) between 4.8 GHz and 36.8 GHz using the University of Michigan Radio Astronomical Observatory (UMRAO), the Crimean Astrophysical Observatory (CrAO), and Metsahovi Radio Observatory for over 40 years. Flares were Gaussian fitted to derive time delays between observed frequencies for each flare ($\Delta t$), peak amplitude ($A$), and their half width. Using $A \propto \nu^{\alpha}$ we infer $\alpha$ in the range $-$16.67 to 2.41 and using $\Delta t \propto \nu^{1/k_r}$, we infer $k_r \sim 1$, employed in the context of equipartition between magnetic and kinetic energy density for parameter estimation. From the estimated core position offset ($\Omega_{r \nu}$) and the core radius ($r_{\rm core}$), we infer that opacity model may not be valid in all cases. The mean magnetic field strength at 1 pc ($B_1$) and at the core ($B_{\rm core}$), are in agreement with previous estimates. We apply the magnetically arrested disk model to estimate black hole spins in the range $0.15-0.9$ for these blazars, indicating that the model is consistent with expected accretion mode in such sources. The power law shaped power spectral density has slopes $-$1.3 to $-$2.3 and is interpreted in terms of multiple shocks or magnetic instabilities.

Origin of magnetic fields, its structure and effects on dynamical processes in stars to galaxies are not well understood. Lack of a direct probe has hampered its study. The first phase of Square Kilometre Array (SKA-I), will have more than an order of magnitude higher sensitivity than existing radio telescopes. In this contribution, we discuss specific science cases that are of interest to the Indian community concerned with astrophysical turbulence and magnetic fields. The SKA-I will allow observations of a large number of background sources with detectable polarisation and measure their Faraday depths (FDs) through the Milky Way, other galaxies and their circum-galactic medium. This will probe line-of-sight magnetic fields in these objects well and provide field configurations. Detailed comparison of observational data with models which consider various processes giving rise to field amplification and maintenance will then be possible. Such observations will also provide the coherence scale of the fields and measure its random component. Measuring the random component is important to characterise turbulence in the medium. Observations of FDs with redshift will provide important information on magnetic field evolution as a function of redshift. The background sources could also be used to probe magnetic fields and its coherent scale in galaxy clusters and in bridges formed between interacting galaxies. Other than FDs, sensitive observations of synchrotron emission from galaxies will provide complimentary information on their magnetic field strengths in the sky plane. The core shift measurements of AGNs can provide more precise measurements of magnetic field very close (<pc) to the black hole and its evolution. The low band of SKA-I will also be useful to study circularly polarized emission from Sun and comparing various models of field configurations with observations.

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 \pi$ 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.

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 $\alpha_m$, even before solving the dynamical equations for evolution of magnetic fields in the disk and the corona, along with $\alpha$-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.

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 $\rho (r) \propto r^{-\gamma}$ and stellar mass function $\xi(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}^\beta$, where $\beta= -0.3\pm 0.01$ for $M_{\bullet}>10^7 M_{\odot}$ and $\sim 6.8 \hspace{1mm} \times 10^{-5}$ Yr$^{-1}$ for $\gamma=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 $\gamma$.

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 $\Delta t$ between the observation frequencies $\nu$ and spectral indices $\alpha$ based on peak amplitudes $A$ are determined. From the fit $\Delta t \propto \nu^{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 \nu^\alpha$, $\alpha$ 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 $\Omega_{r\nu} = 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 $\nu$, 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 $\alpha$, 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.

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 $\gamma_{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.

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.

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 $\psi\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.

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.

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.

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.

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.

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.

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.

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.

A general type of mathematical argument is described, which applies to all the cases in which dynamo maintenance of a steady magnetic field by motion in a uniform density is known to be impossible. Previous work has demonstrated that magnetic field decay is unavoidable under conditions of axisymmetry and in spherical or planar incompressible flows. These known results are encompassed by a calculation for flows described in terms of a generalized poloidal-toroidal representation of the magnetic field with respect to an arbitrary two dimensional surface. We show that when the velocity field is two dimensional, the dynamo growth, if any, that results, is linear in one of the projections of the field while the other projections remain constant. We also obtain criteria for the existence of and classification into two and three dimensional velocity results which are satisfied by a restricted set of geometries. In addition, we discuss the forms of spatial variation of the density and the resistivity that are allowed so that field decay still occurs for this set of geometries.

Supermassive disks are thought to be precursors of supermassive black holes that are believed to power quasars and exist at centers of galaxies. Formation scenarios of such disks are reviewed and it is argued that gas dynamical schemes are favourable compared to stellar dynamical schemes which could however be important feeding mechanisms for the growth of the black hole. A new self-similar model of a collapse of a self-gravitating disk due to radiation induced stresses applicable to two different situations of radiative viscosity and Compton drag is presented. The collapse timescale purely due to radiative viscosity is found to be a fraction of Hubble time, $\tau_\gamma \sim{\sigma_T c/ (m_p G)} ({L_{edd}/ L}) \simeq 6 \times 10^9$yrs is slow and probably magnetic fields play an important role before general relativistic effects take over. A model of self-gravitating disk collapsing due to Compton drag by the Cosmic Microwave Background is also presented which is found to be effective at redshifts $1400>z\gsim 300$. It is proposed that the small $\lsim 10^5 \Msun$ objects that form by this mechanism by $z \sim 20$ can merge and coalesce by dynamical friction to form the high redshift quasars seen. Supermassive stars which are systems (and could be end products of a supermassive disk phase) en route to the final collapse are also briefly reviewed.

We present a detailed model for the formation of massive objects at the centers of galaxies. The effects of supernovae heating and the conditions of gas loss are revisited. The escape time of the gas is compared with the cooling time, which provides an additional condition not previously considered. Its consequences for the allowed mass range of the halo is calculated. It is shown that sufficient gas is retained to form massive dark objects and quasars even for moderately massive halos but a decline is expected at low redshifts. Subsequently, a gaseous disk forms with a radial extent of a ~kpc, spun up by tidal torques and magnetized by supernovae fields with fields strengths of $10-100 \mu G$. In a model of a self-similar accretion flow in an initially dominant halo, it is shown that for typical halo parameters, about $10^8 M_\odot$ accretes via small magnetic stresses (or alternatively by self-gravity induced instability or by alpha viscosity) in $10^8$ years into a compact region. A model of a self-gravitating evolution of a compact magnetized disk, which is relevant when a significant fraction of the disk mass falls in, is presented, and it has a rapid collapse time scale of a million years. The two disk solutions, one for accretion in an imposed halo potential and the other for self-gravitating disk, obtained here, have general utility and can be adapted to other contexts like protostellar disks as well.

Violent relaxation during the collapse of a galaxy halo is known to be incomplete in realistic cases such as cosmological infall or mergers. We adopt a physical picture of strong but short lived interactions between potential fluctuations and particle orbits, using the broad framework outlined by Tremaine (1987) for incorporating incompleteness of the relaxation. We are guided by results from plasma physics, viz. the quasilinear theory of Landau damping, but allow for significant differences in our case. Crucially, wave particle scattering does not drive the system to an equilibrium distribution function of the exponential type, even in regions of phase space allowed by the constraints. The physical process is mixing without friction in ``action'' space, for which the simplest possible model is a constant phase space density modulated by the constraints. Our distribution function does not use the exponential functions of the energy prevalent in other work, which we regard as inappropriate to collisionless systems. The halo of the self-consistent, parameter-free solutions show an r^(-4) behavior in density at large r, an r^(1/4) surface brightness profile in the region 0.1-8 r_e, and a radially anisotropic velocity dispersion profile outside an isotropic core. The energy distribution seen in simulations N(E) singles out the pericenter cutoff model as the most realistic among the variants we have explored.