The role of Compton and Raman scattering in the quasar continuum

There are three ways in which an electromagnetic wave can undergo scattering in a plasma: (i) when the scattering of radiation occurs by a single electron, it is called Compton Scattering (CS); (ii) if it occurs by a longitudinal electron plasma mode, it is called Stimulated Raman Scattering (SRS), and (iii) if it occurs by a highly damped electron plasma mode, it is called Stimulated Compton Scattering (SCS). The non-thermal continuum of quasars is believed to be produced through the combined action of synchrotron and inverse Compton processes, which are essentially single-particle processes. Here, we investigate the role of SRS and SCS in the generation of continuum radiation from these compact objects. It is shown as an example that the complete spectrum of 3C 273 can be reproduced by suitably combining SCS and SRS. The differential contributions of SCS and SRS under different values of the plasma parameters are also calculated.

Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. It has been proposed that some of these outflows are driven off gas infalling towards the central supermassive black hole. We perform radiative transfer calculations to compute the gas ionization state and the emergent X-ray spectra for both two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of relatively dense (10 -10 g cm ) outflowing gas. This outflow causes strong absorption for observer orientation angles of ?35°. Particularly in 3D, the strength of this absorption varies significantly for different lines of sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features in the X-ray spectra since the ionization state is predicted to be very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a Compton-scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow are predicted to give rise to several emission features including a multicomponent Fe Ka emission complex for all observer orientations. For the hydrodynamical simulations considered, we predict both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. Higher resolution hydrodynamic simulations are needed to determine whether the outflows fragment on scales unresolved in our current study, which may yield the denser lower ionization material that could reconcile the models and the observations. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Variability of the high-velocity outflow in the quasar PDS 456

We present a comparison of two Suzaku X-ray observations of the nearby (z = 0.184), luminous (L ∼ 10 erg s) type I quasar, PDS 456. A new 125 ks Suzaku observation in 2011 caught the quasar during a period of low X-ray flux and with a hard X-ray spectrum, in contrast with a previous 190 ks Suzaku observation in 2007 when the quasar appeared brighter and had a steep (Γ > 2) X-ray spectrum. The 2011 X-ray spectrum contains a pronounced trough near 9 keV in the quasar rest frame, which can be modeled with blueshifted iron K-shell absorption, most likely from the He- and H-like transitions of iron. The absorption trough is observed at a similar rest-frame energy as in the earlier 2007 observation, which appears to confirm the existence of a persistent high-velocity wind in PDS 456, at an outflow velocity of 0.25-0.30c. The spectral variability between 2007 and 2011 can be accounted for by variations in a partial covering absorber, increasing in covering fraction from the brighter 2007 observation to the hard and faint 2011 observation. Overall, the low-flux 2011 observation can be explained if PDS 456 is observed at relatively low inclination angles through a Compton-thick wind, originating from the accretion disk, which significantly attenuates the X-ray flux from the quasar. © 2014. The American Astronomical Society. All rights reserved.

Constraints from the old quasar APM 08279+5255 on two classes of Lambda(t)-cosmologies

The viability of two different classes of Lambda(t)CDM cosmologies is tested by using the APM 08279+5255, an old quasar at redshift z = 3.91. In the first class of models, the cosmological term scales as Lambda(t) similar to R(-n). The particular case n = 0 describes the standard Lambda CDM model whereas n = 2 stands for the Chen and Wu model. For an estimated age of 2 Gyr, it is found that the power index has a lower limit n > 0.21, whereas for 3 Gyr the limit is n > 0.6. Since n can not be so large as similar to 0.81, the Lambda CDM and Chen and Wu models are also ruled out by this analysis. The second class of models is the one recently proposed by Wang and Meng which describes several Lambda(t)CDM cosmologies discussed in the literature. By assuming that the true age is 2 Gyr it is found that the epsilon parameter satisfies the lower bound epsilon > 0.11 while for 3 Gyr, a lower limit of epsilon > 0.52 is obtained. Such limits are slightly modified when the baryonic component is included.

A faraway quasar in the direction of the highest energy Auger event

The highest energy cosmic ray event reported by the Auger Observatory has an energy of 148 EeV. It does not correlate with any nearby (z<0.024) object capable of originating such a high energy event. Intrigued by the fact that the highest energy event ever recorded (by the Fly`s Eye collaboration) points to a faraway quasar with very high radio luminosity and large Faraday rotation measurement, we have searched for a similar source for the Auger event. We find that the Auger highest energy event points to a quasar with similar characteristics to the one correlated to the Fly`s Eye event. We also find the same kind of correlation for one of the highest energy AGASA events. We conclude that so far these types of quasars are the best source candidates for both Auger and Fly`s Eye highest energy events. We discuss a few exotic candidates that could reach us from gigaparsec distances.

Quasar - variabilità e casualità

Analisi della variabilità temporale dello spettro di AGN e in particolare di Quasar, con definizione di "quasi-periodicità" e "casualità" su lungo e breve periodo, per un suo uso astrofisico mirato alla creazione di un cifrario crittografico.

The 2dF-SDSS LRG and QSO (2SLAQ) survey: the z < 2.1 quasar luminosity function from 5645 quasars to g=21.85

We have used the Two-Degree Field (2dF) instrument on the Anglo-Australian Telescope (AAT) to obtain redshifts of a sample of z < 3 and 18.0 < g < 21.85 quasars selected from Sloan Digital Sky Survey (SDSS) imaging. These data are part of a larger joint programme between the SDSS and 2dF communities to obtain spectra of faint quasars and luminous red galaxies, namely the 2dF-SDSS LRG and QSO (2SLAQ) Survey. We describe the quasar selection algorithm and present the resulting number counts and luminosity function of 5645 quasars in 105.7 deg(2). The bright-end number counts and luminosity functions agree well with determinations from the 2dF QSO Redshift Survey (2QZ) data to g similar to 20.2. However, at the faint end, the 2SLAQ number counts and luminosity functions are steeper (i.e. require more faint quasars) than the final 2QZ results from Croom et al., but are consistent with the preliminary 2QZ results from Boyle et al. Using the functional form adopted for the 2QZ analysis ( a double power law with pure luminosity evolution characterized by a second-order polynomial in redshift), we find a faint-end slope of beta =-1.78 +/- 0.03 if we allow all of the parameters to vary, and beta =-1.45 +/- 0.03 if we allow only the faint-end slope and normalization to vary (holding all other parameters equal to the final 2QZ values). Over the magnitude range covered by the 2SLAQ survey, our maximum-likelihood fit to the data yields 32 per cent more quasars than the final 2QZ parametrization, but is not inconsistent with other g > 21 deep surveys for quasars. The 2SLAQ data exhibit no well-defined 'break' in the number counts or luminosity function, but do clearly flatten with increasing magnitude. Finally, we find that the shape of the quasar luminosity function derived from 2SLAQ is in good agreement with that derived from Type I quasars found in hard X-ray surveys.

Análise wavelet da variabilidade do quasar 3C 273

Discovered in 1963, 3C 273 was the second quasar identified and cataloged in the Third Cambridge Catalog for radio sources, and the first one for which emission lines were identified with a hydrogen sequence redshifted. It is the brightest quasar of the celestial sphere, the most studied, analyzed, and with a resulting abundance of data available in a vast literature. The accurate analysis of the deviations of the spectral lines of quasars provides enough information to put in evidence the variation of fundamental constants of nature and similarly the universe expansion rate. The analysis of the variability of the light curves of these bodies, and the consequent accuracy of their periodicity, is of utmost importance as it provides an efficiency of their observations, enables a greater understanding of the physical phenomena, and makes it possible to conduct spectral observations on more accurate dates (when their light curves show pronounced peaks and therefore richer spectra information). In this master’s thesis twenty eight light curves from the quasar 3C 273 are studied, covering all the electromagnetic spectrum wavebands (radio emission to gamma rays), totaling in the analysis of four light curves for each waveband. We have applied the method of Continuous Wavelet Transform using the sixth-order (!0 = 6) Morlet wavelet function, and obtained excellent results in accordance with the literature.

Análise wavelet da variabilidade do quasar 3C 273

Discovered in 1963, 3C 273 was the second quasar identified and cataloged in the Third Cambridge Catalog for radio sources, and the first one for which emission lines were identified with a hydrogen sequence redshifted. It is the brightest quasar of the celestial sphere, the most studied, analyzed, and with a resulting abundance of data available in a vast literature. The accurate analysis of the deviations of the spectral lines of quasars provides enough information to put in evidence the variation of fundamental constants of nature and similarly the universe expansion rate. The analysis of the variability of the light curves of these bodies, and the consequent accuracy of their periodicity, is of utmost importance as it provides an efficiency of their observations, enables a greater understanding of the physical phenomena, and makes it possible to conduct spectral observations on more accurate dates (when their light curves show pronounced peaks and therefore richer spectra information). In this master’s thesis twenty eight light curves from the quasar 3C 273 are studied, covering all the electromagnetic spectrum wavebands (radio emission to gamma rays), totaling in the analysis of four light curves for each waveband. We have applied the method of Continuous Wavelet Transform using the sixth-order (!0 = 6) Morlet wavelet function, and obtained excellent results in accordance with the literature.

Testing quasar unification: radiative transfer in clumpy winds

Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to 1043 erg s-1. Here, we introduce a simple treatment of clumping, and find that a filling factor of ˜0.01 moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Lyα and C IV 1550 Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question.

The High-Redshift Quasar Luminosity Function from Multi-Epoch Imaging Surveys

Thesis (Ph.D.)--University of Washington, 2016-08

The development of Monte Carlo techniques for the verification of radiotherapy treatments

Introduction: Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo (MC) methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the ‘gold standard’ for predicting dose deposition in the patient [1]. This project has three main aims: 1. To develop tools that enable the transfer of treatment plan information from the treatment planning system (TPS) to a MC dose calculation engine. 2. To develop tools for comparing the 3D dose distributions calculated by the TPS and the MC dose engine. 3. To investigate the radiobiological significance of any errors between the TPS patient dose distribution and the MC dose distribution in terms of Tumour Control Probability (TCP) and Normal Tissue Complication Probabilities (NTCP). The work presented here addresses the first two aims. Methods: (1a) Plan Importing: A database of commissioned accelerator models (Elekta Precise and Varian 2100CD) has been developed for treatment simulations in the MC system (EGSnrc/BEAMnrc). Beam descriptions can be exported from the TPS using the widespread DICOM framework, and the resultant files are parsed with the assistance of a software library (PixelMed Java DICOM Toolkit). The information in these files (such as the monitor units, the jaw positions and gantry orientation) is used to construct a plan-specific accelerator model which allows an accurate simulation of the patient treatment field. (1b) Dose Simulation: The calculation of a dose distribution requires patient CT images which are prepared for the MC simulation using a tool (CTCREATE) packaged with the system. Beam simulation results are converted to absolute dose per- MU using calibration factors recorded during the commissioning process and treatment simulation. These distributions are combined according to the MU meter settings stored in the exported plan to produce an accurate description of the prescribed dose to the patient. (2) Dose Comparison: TPS dose calculations can be obtained using either a DICOM export or by direct retrieval of binary dose files from the file system. Dose difference, gamma evaluation and normalised dose difference algorithms [2] were employed for the comparison of the TPS dose distribution and the MC dose distribution. These implementations are spatial resolution independent and able to interpolate for comparisons. Results and Discussion: The tools successfully produced Monte Carlo input files for a variety of plans exported from the Eclipse (Varian Medical Systems) and Pinnacle (Philips Medical Systems) planning systems: ranging in complexity from a single uniform square field to a five-field step and shoot IMRT treatment. The simulation of collimated beams has been verified geometrically, and validation of dose distributions in a simple body phantom (QUASAR) will follow. The developed dose comparison algorithms have also been tested with controlled dose distribution changes. Conclusion: The capability of the developed code to independently process treatment plans has been demonstrated. A number of limitations exist: only static fields are currently supported (dynamic wedges and dynamic IMRT will require further development), and the process has not been tested for planning systems other than Eclipse and Pinnacle. The tools will be used to independently assess the accuracy of the current treatment planning system dose calculation algorithms for complex treatment deliveries such as IMRT in treatment sites where patient inhomogeneities are expected to be significant. Acknowledgements: Computational resources and services used in this work were provided by the HPC and Research Support Group, Queensland University of Technology, Brisbane, Australia. Pinnacle dose parsing made possible with the help of Paul Reich, North Coast Cancer Institute, North Coast, New South Wales.

The effect of organ motion of radiotherapy plans

Organ motion as a result of respiration is an important field of research for medical physics. Knowledge of magnitude and direction of this motion is necessary to allow for more accurate radiotherapy treatment planning. This will result in higher doses to the tumour whilst sparing healthy tissue. This project involved human trials, where the radiation therapy patient's kidneys were CT scanned under three different conditions; whilst free breathing (FB), breath-hold at normal tidal inspiration (BHIN), and breath-hold at normal tidal expiration (BHEX). The magnitude of motion was measured by recording the outline of the kidney from a Beam's Eye View (BEV). The centre of mass of this 2D shape was calculated for each set using "ImageJ" software and the magnitude of movement determined from the change in the centroid's coordinates between the BHIN and BHEX scans. The movement ranged from, for the left and right kidneys, 4-46mm and 2-44mm in the superior/inferior (axial) plane, 1-21mm and 2- 16mm in the anterior/posterior (coronal) plane, and 0-6mm and 0-8mm in the lateral/medial (sagittal) plane. From exhale to inhale, the kidneys tended to move inferiorly, anteriorly and laterally. A standard radiotherapy plan, designed to treat the para-aortics with opposed lateral fields was performed on the free breathing (planning) CT set. The field size and arrangement was set up using the same parameters for each subject. The prescription was to deliver 45 Gray in 25 fractions. This field arrangement and prescription was then copied over to the breath hold CT sets, and the dosimetric differences were compared using Dose Volume Histograms (DVH). The point of comparison for the three sets was recorded as the percentage volume of kidney receiving less than or equal to 10 Gray. The QUASAR respiratory motion phantom was used with the range of motion determined from the human study. The phantom was imaged, planned and treated with a linear accelerator with dose determined by film. The effect of the motion was measured by the change in the penumbra of the film and compared to the penumbra from the treatment planning system.