Engineering PapersSearch

Engineering topics

Kazanas, Demosthenes

Publications and source records attributed to Kazanas, Demosthenes.

At least 19 records

Dust Formation in AGN Winds

Infrared observations of active galactic nuclei (AGNs) reveal emission from the putative dusty circumnuclear "torus" invoked by AGN unification, which is heated up by radiation from the central accreting black hole (BH).The strong 9.7 and 18 m silicate features observed in the AGN spectra, in both emission and absorption, further indicate the presence of such dusty environments. We present detailed calculations of the chemistry of silicate dust formation in AGN accretion disk winds. The winds considered herein are magnetohydrodynamic winds driven off the entire accretion disk domain that extends from the BH vicinity to the radius of BH influence, of order ~1-100 pc depending on the mass of the resident BH. Our results indicate that these winds provide conditions conducive to the formation of significant amounts of dust, especially for objects accreting close to their Eddington limit, making AGNs a significant source of dust in the universe, especially for luminous quasars. Our models justify the importance of an r1 density law in the winds for efficient formation and survival of dust grains. The dust production rate scales linearly with the mass of the central BH and varies as a power law of index between 2 and 2.5 with the dimensionless mass accretion rate. The resultant distribution of the dense dusty gas resembles atoroidal shape, with high column density and optical depths along the equatorial viewing angles, in agreement with the AGN unification picture

Arkaprabha Sarangi

Meta-Analysis of Electron Cyclotron Resonance Absorption Features Detected in High-Mass X-Ray Binaries

Using recent compilations of detailed X-ray observations and spectral models of exceptional quality, we record the electron cyclotron resonance absorption (ECRA) features that have been detected in 45 pulsating high-mass X-ray binaries (HMXBs) and ultraluminous X-ray (ULX) sources harboring neutron stars, although seven of these detections are still questionable and another 21 are single and/or not independently confirmed. From the comprehensive catalogs of Jaisawal & Naik and Staubert et al. and from several additional recent observations, we produce two lists of HMXB ECRA sources: a list of 17 sources in which multiple ECRA lines or single very low-energy lines are seen, in which we can reasonably assume that the lowest energy reveals the fundamental cyclotron level for each source; and a "contaminated" list of 38 sources including the 21 detections of single ECRA lines that may (not) be higher-level harmonics. Both lists confirm a previous result that we have obtained independently by modeling the propeller lines of Magellanic HMXB pulsars: the surface dipolar magnetic fields B * of HMXB neutron stars are segregated around five distinct values with B * = 0.28 ± 0.08, 0.55 ± 0.11, 1.3 ± 0.37, 3.0 ± 0.68 and 7.9 ± 3.1, in units of TG. However, an explanation of this phenomenon is currently lacking. We have found no correlation between these B * values and the corresponding observed spin periods, spin period derivatives, orbital periods, maximum X-ray luminosities, neutron star masses or companion star masses.

Christodoulou, Dimitris M.

A Fundamental Plane for Gamma-Ray Pulsars

We show that the γ-ray pulsar observables, i.e., their total γ-ray luminosity, L(sub γ), spectral cutoff energy, ϵ (sub cut), stellar surface magnetic field, B(sub ⋆), and spin-down power ε, obey a relation of the form L(sub γ) = f (ϵ (sub cut), B(sub ⋆), ε ), which represents a 3D plane in their 4D logspace. Fitting the data of 88 pulsars of the second Fermi pulsar catalog, we show this relation to be L(sub γ) ∝ ϵ (sub cut)(sup 1.18±0.24 B(sub ⋆)(sup 0.17± 0.05) ε(sup 0.41±0.08), a pulsar fundamental plane (FP). We show that the observed FP is remarkably close to the theoretical relation L(sub γ) ∝ ϵ (sub cut)(sup 4/3) B(sub ⋆)(sup 1/6) ε(sup 5/12) obtained assuming that the pulsar γ-ray emission is due to curvature radiation by particles accelerated at the pulsar equatorial current sheet just outside the light cylinder. Interestingly, the FP seems incompatible with emission by synchrotron radiation. The corresponding scatter about the FP is ∼0.35 dex and can only partly be explained by the observational errors while the rest is probably due to the variation of the inclination and observer angles. We predict also that ϵ (sub cut) ∝ ε (sup 7/16) toward low ε for both young and millisecond pulsars implying that the observed death line of γ-ray pulsars is due to ϵ (sub cut) dropping below the Fermi band. Our results provide a comprehensive interpretation of the observations of γ-ray pulsars, setting requirements for successful theoretical modeling.

Kalapotharakos, Constantinos

Universal Expansion with Spatially Varying G

We calculate the expansion of the Universe under the assumptions that G varies in space and the radial size r of the Universe is very large (we call this the MOND regime of varying-G gravity). The inferred asymptotic behaviour turns out to be different from that found by McCrea & Milne in 1934 and our equations bear no resemblance to those of the relativistic case. In this cosmology, the scale factor R(t) increases linearly with time t, the radial velocity is driven by inertia, and gravity is incapable of hindering the expansion. Yet, Hubble’s law is borne out without any additional assumptions. When we include a repulsive acceleration a(sub de) due to dark energy, the resulting universal expansion is then driven totally by this new term and the solutions for a(sub de) → 0 do not reduce to those of the a(sub de) ≡ 0 case. This is a realization of a new Thom catastrophe: The inclusion of the new term alters the conservation of energy and the dark energy solutions are not reducible to those in the case without dark energy.

Christodoulou, Dimitris M.

Gravitational Potential and Non-relativistic Lagrangian in Modified Gravity with Varying G

We have recently shown that the baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations imply that the gravitational "constant" G in the force law vary with acceleration alpha as 1/alpha . Here we derive the converse from first principles. First we obtain the gravitational potential for all accelerations and we formulate the Lagrangian for the central-force problem. Then action minimization implies the BTF/BFJ relations in the deep MOND limit as well as weak-field Weyl gravity in the Newtonian limit. The results show how we can properly formulate a nonrelativistic conformal theory of modified dynamics that reduces to MOND in its low-acceleration limit and to Weyl gravity in the opposite limit. An unavoidable conclusion is that alpha-0 , the transitional acceleration in modified dynamics, does not have a cosmological origin and it may not even be constant among galaxies and galaxy clusters.

Christodoulou, Dimitris M.

Accretion Disc MHD Winds and Blazar Classification

The Fermi Gamma-Ray Space Telescope observations of blazars show a strong correlation between the spectral index of their -ray spectra and their synchrotron peak frequency synpk ;additionally, the rate of Compton dominance of these sources also seems to be a function of synpk .In this work, we adopt the assumption that the non-thermal emission of blazars is primarily due to radiation by a population of Fermi-accelerated electrons in a relativistic outflow (jet) along the symmetry axis of the blazar's accretion disc. Furthermore, we assume that the Compton component is related to an external photon field of photons, which are scattered from particles of the magnetohydrodynamic wind emanating from the accretion disc. Our results reproduce well the aforementioned basic observational trends of blazar classification by varying just oneparameter, namely the mass accretion rate on to the central black hole.

Boula, Stella

Flaring Activity from Quiescent States in Neutron-Star X-Ray Binaries

We examine systematically the observed X-ray luminosity jumps (or flares) from quiescent states in millisecond binary pulsars (MSBPs) and high-mass X-ray binary pulsars (HMXBPs). We rely on the published X-ray light curves of seven pulsars: four HMXBPs, two MSBPs and the ultraluminous X-ray pulsar M82 X-2. We discuss the physics of their flaring activities or lack thereof, paying special attention to their emission properties when they are found on the propeller line, inside the Corbet gap or near the light-cylinder barrier.We provide guiding principles for future interpretations of faint X-ray observations, as well as a method of constraining the propeller lines and the dipolar surfacemagnetic fields of pulsars using a variety of quiescent states. In the process, we clarify some disturbing inaccuracies that have made their way into the published literature.

Christodoulou, Dimitris M.

Magnetized Disk Winds in NGC 3783

We analyze a 900 kilosecond stacked Chandra/HETG (High-Energy Transmission Grating) spectrum of NGC 3783 in the context of magnetically driven accretion-disk wind models in an effort to provide tight constraints on the global conditions of the underlying absorbers. Motivated by the earlier measurements of its absorption measure distribution (AMD) indicating X-ray-absorbing ionic columns that decrease slowly with decreasing ionization parameter, we employ 2-dimension (2-D) magnetohydrodynamic (MHD) disk wind models to describe the global outflow. We compute its photoionization structure along with the wind kinematic properties, allowing us to further calculate in a self-consistent fashion the shapes of the major X-ray absorption lines. With the wind radial density profile determined by the AMD, the profiles of the ensemble of the observed absorption features are determined by the two global parameters of the MHD wind; i.e., disk inclination theta (sub obs) and wind density normalization n (sub o). Considering the most significant absorption features in the approximately 1.8-20 angstrom range, we show that the MHD wind is best described by n(r) approximately equal to 6.9 times 10 (sup 11) (r/r (sub o)) (sup - 1.15) cubic centimeters and theta (sub obs). We argue that winds launched by X-ray heating or radiation pressure, or even MHD winds but with steeper radial density profiles, are strongly disfavored by data. Considering the properties of Fe K-band absorption features (i.e., Fe XXV and Fe XXVI), while typically prominent in the active galactic nucleus X-ray spectra, they appear to be weak in NGC 3783. For the specific parameters of our model obtained by fitting the AMD and the rest of the absorption features, these features are found to be weak, in agreement with observations.

Fukumura, Keigo

A Physical Interpretation of the Titius-Bode Rule and Its Connection to the Closed Orbits of Bertrand's Theorem

We consider the geometric Titius-Bode rule for the semimajor axes of planetary orbits. We derive an equivalent rule for the midpoints of the segments between consecutive orbits along the radial direction and we interpret it physically in terms of the work done in the gravitational field of the Sun by particles whose orbits are perturbed around each planetary orbit. On such energetic grounds, it is not surprising that some exoplanets in multiple-planet extrasolar systems obey the same relation. However, it is surprising that this simple interpretation of the Titius-Bode rule also reveals new properties of the bound closed orbits predicted by Bertrand's theorem, which has been known since 1873.

Titus-Bode rule

Fermi Gamma-Ray Pulsars: Understanding the High-Energy Emission from Dissipative Magnetospheres

Based on the Fermi observational data, we reveal meaningful constraints for the dependence of the macroscopic conductivity (sigma) of dissipative pulsar magnetosphere models on the corresponding spin-down rate, epsilon(sup dot). Our models are refinements of the FIDO (Force-free Inside, Dissipative Outside) models, which have dissipative regions that are restricted on the equatorial current sheet outside of the the light-cylinder. Taking into account the observed cutoff energies of all of the Fermi pulsars and assuming that (a) the corresponding gamma-ray pulsed emission is due to curvature radiation at the radiation-reaction-limit regime, and (b) this emission is produced at the equatorial current sheet near the light cylinder, we show that the Fermi data provide clear indications about the corresponding accelerating electric-field components. A direct comparison between the Fermi cutoff energies and the model ones reveals that if sigma increases with epsilon(sup dot) for high epsilon(sup dot)-values, while it saturates for low ones. This comparison indicates also that the corresponding gap width increases toward low epsilon(sup dot)-values. Assuming the Goldreich-Julian flux for the emitting particles, we calculate the total gamma-ray luminosity (L(sub gamma)). A comparison between the dependence of the Fermi L(sub gamma)-values and the model ones on epsilon(sup dot) indicates an increase of the emitting particle multiplicity with epsilon(sup dot). Our modeling, guided by the Fermi data alone, enhances our understanding of the physical mechanisms behind the high-energy emission in pulsar magnetospheres.

Kalapotharakos, Constantinos

Testing Dissipative Magnetosphere Model Light Curves and Spectra with Fermi Pulsars

We explore the emission properties of a dissipative pulsar magnetosphere model introduced by Kalapotharakos et al. comparing its high-energy light curves and spectra, due to curvature radiation, with data collected by the Fermi LAT. The magnetosphere structure is assumed to be near the force-free solution. The accelerating electric field, inside the light cylinder (LC), is assumed to be negligible, while outside the LC it rescales with a finite conductivity (sigma). In our approach we calculate the corresponding high-energy emission by integrating the trajectories of test particles that originate from the stellar surface, taking into account both the accelerating electric field components and the radiation reaction forces. First, we explore the parameter space assuming different value sets for the stellar magnetic field, stellar period, and conductivity. We show that the general properties of the model are in a good agreement with observed emission characteristics of young gamma-ray pulsars, including features of the phase-resolved spectra. Second, we find model parameters that fit each pulsar belonging to a group of eight bright pulsars that have a published phase-resolved spectrum. The sigma values that best describe each of the pulsars in this group show an increase with the spin-down rate (E˙ ) and a decrease with the pulsar age, expected if pair cascades are providing the magnetospheric conductivity. Finally, we explore the limits of our analysis and suggest future directions for improving such models.

acceleration of particles

Gamma-Ray Emission in Dissipative Pulsar Magnetospheres: from Theory to Fermi Observations

We compute the patterns of gamma-ray emission due to curvature radiation in dissipative pulsar magnetospheres. Our ultimate goal is to construct macrophysical models that are able to reproduce the observed gamma-ray light curve phenomenology recently published in the Second Fermi Pulsar Catalog. We apply specific forms of Ohm's law on the open field lines using a broad range for the macroscopic conductivity values that result in solutions ranging, from near-vacuum to near-force-free. Using these solutions, we generate model gamma-ray light curves by calculating realistic trajectories and Lorentz factors of radiating particles under the influence of both the accelerating electric fields and curvature radiation reaction. We further constrain our models using the observed dependence of the phase lags between the radio and gamma-ray emission on the gamma-ray peak separation. We perform a statistical comparison of our model radio-lag versus peak-separation diagram and the one obtained for the Fermi standard pulsars. We find that for models of uniform conductivity over the entire open magnetic field line region, agreement with observations favors higher values of this parameter. We find, however, significant improvement in fitting the data with models that employ a hybrid form of conductivity, specifically, infinite conductivity interior to the light cylinder and high but finite conductivity on the outside. In these models the gamma-ray emission is produced in regions near the equatorial current sheet but modulated by the local physical properties. These models have radio lags near the observed values and statistically best reproduce the observed light curve phenomenology. Additionally, they also produce GeV photon cut-off energies.

stars

A New Standard Pulsar Magnetosphere

In view of recent efforts to probe the physical conditions in the pulsar current sheet, we revisit the standard solution that describes the main elements of the ideal force-free pulsar magnetosphere. The simple physical requirement that the electric current contained in the current layer consists of the local electric charge moving outward at close to the speed of light yields a new solution for the pulsar magnetosphere everywhere that is ideal force-free except in the current layer. The main elements of the new solution are as follows: (1) the pulsar spindown rate of the aligned rotator is 23% larger than that of the orthogonal vacuum rotator; (2) only 60% of the magnetic flux that crosses the light cylinder opens up to infinity; (3) the electric current closes along the other 40%, which gradually converges to the equator; (4) this transfers 40% of the total pulsar spindown energy flux in the equatorial current sheet, which is then dissipated in the acceleration of particles and in high-energy electromagnetic radiation; and (5) there is no separatrix current layer. Our solution is a minimum free-parameter solution in that the equatorial current layer is electrostatically supported against collapse and thus does not require a thermal particle population. In this respect, it is one more step toward the development of a new standard solution. We discuss the implications for intermittent pulsars and long-duration gamma-ray bursts. We conclude that the physical conditions in the equatorial current layer determine the global structure of the pulsar magnetosphere.

pulsars: general

Stratified Magnetically Driven Accretion-Disk Winds and Their Relations To Jets

We explore the poloidal structure of two-dimensional magnetohydrodynamic (MHD) winds in relation to their potential association with the X-ray warm absorbers (WAs) and the highly ionized ultra-fast outflows (UFOs) in active galactic nuclei (AGNs), in a single unifying approach. We present the density n(r, theta), ionization parameter xi(r, theta), and velocity structure v(r, theta) of such ionized winds for typical values of their fluid-to-magnetic flux ratio, F, and specific angular momentum, H, for which wind solutions become super-Alfvenic. We explore the geometrical shape of winds for different values of these parameters and delineate the values that produce the widest and narrowest opening angles of these winds, quantities necessary in the determination of the statistics of AGN obscuration. We find that winds with smaller H show a poloidal geometry of narrower opening angles with their Alfv´en surface at lower inclination angles and therefore they produce the highest line of sight (LoS) velocities for observers at higher latitudes with the respect to the disk plane. We further note a physical and spatial correlation between the X-ray WAs and UFOs that form along the same LoS to the observer but at different radii, r, and distinct values of n, xi, and v consistent with the latest spectroscopic data of radio-quiet Seyfert galaxies. We also show that, at least in the case of 3C 111, the winds' pressure is sufficient to contain the relativistic plasma responsible for its radio emission. Stratified MHD disk winds could therefore serve as a unique means to understand and unify the diverse AGN outflows.

accretion

The Force-Free Magnetosphere of a Rotating Black Hole

We revisit the Blandford-Znajek process and solve the fundamental equation that governs the structure of the steady-state force-free magnetosphere around a Kerr black hole. The solution depends on the distributions of the magnetic field angular velocity and the poloidal electric current. These are not arbitrary. They are determined self-consistently by requiring that magnetic field lines cross smoothly the two singular surfaces of the problem: the inner "light surface" located inside the ergosphere and the outer "light surface" which is the generalization of the pulsar light cylinder.We find the solution for the simplest possible magnetic field configuration, the split monopole, through a numerical iterative relaxation method analogous to the one that yields the structure of the steady-state axisymmetric force-free pulsar magnetosphere. We obtain the rate of electromagnetic extraction of energy and confirm the results of Blandford and Znajek and of previous time-dependent simulations. Furthermore, we discuss the physical applicability of magnetic field configurations that do not cross both "light surfaces."

Contopoulos, Ioannis

A Method for Localizing Energy Dissipation in Blazars Using Fermi Variability

The distance of the Fermi-detected blazar gamma-ray emission site from the supermassive black hole is a matter of active debate. Here we present a method for testing if the GeV emission of powerful blazars is produced within the sub-pc scale broad line region (BLR) or farther out in the pc-scale molecular torus (MT) environment. If the GeV emission takes place within the BLR, the inverse Compton (IC) scattering of the BLR ultraviolet (UV) seed photons that produces the gamma-rays takes place at the onset of the Klein-Nishina regime. This causes the electron cooling time to become practically energy independent and the variation of the gamma-ray emission to be almost achromatic. If on the other hand the -ray emission is produced farther out in the pc-scale MT, the IC scattering of the infrared (IR) MT seed photons that produces the gamma-rays takes place in the Thomson regime, resulting to energy-dependent electron cooling times, manifested as faster cooling times for higher Fermi energies. We demonstrate these characteristics and discuss the applicability and limitations of our method.

Dotson, Amanda

Toward a Unified AGN Structure

We present a unified model for the structure and appearance of accretion powered sources across their entire luminosity range from galactic X-ray binaries (XRB) to luminous quasars, with emphasis on AG N and their phenomenology. Central to this model is the notion of MHD winds launched by the accretion disks that power these objects. These winds provide the matter that manifests as blueshifted absorption features in the UV and X-ray spectra of a large fraction of these sources; furthermore, their density distribution in the poloidal plane determines their "appearance" (i.e. the column and velocity structure of these absorption features and the obscuration of the continuum source) as a function of the observer inclination angle (a feature to which INTEGRAL has made significant contributions). This work focuses on just the broadest characteristics of these objects; nonetheless, it provides scaling laws that allow one to reproduce within this model the properties of objects extending in luminosity from luminous quasars to XRBs. Our general conclusion is that the AGN phenomenology can be accounted for in terms of three parameters: The wind maSS flux in units of the Eddington value, m(dot), the observers' inclination angle Theta and the logarithmic slope between the 0/UV and X-ray fluxes alpha(sub ox); however because of a correlation between alpha(sub ox) and UV luminosity the number of significant parameters is two. The AGN correlations implied by this model appear to extend to and consistent with the XRB phenomenology, suggesting the presence of a truly unified underlying structure for accretion powered sources.

Kazanas, Demosthenes

Conformal Weyl Gravity and Perihelion Precession

We investigate the perihelion shift of planetary motion in conformal Weyl gravity using the metric of the static, spherically symmetric solution discovered by Mannheim and Kazanas. To this end we employ a procedure similar to that used by Weinberg for the Schwarzschild solution, which has also been used recently to study the solar system effects of the cosmological constant Lambda. We show that besides the general relativistic terms obtained earlier from the Schwarzschild-de Sitter solution, the expression for the perihelion shift includes a negative contribution which arises from the linear term gamma ray in the metric. Using data for perihelion shift observations, we obtain constraints on the value of the constant gammma similar to that obtained earlier using galactic rotational curves.

Sultana, Joseph