Magnetic fields of relativistic jets in active galactic nuclei

The motivation of my work is to understand the role played by the magnetic field in the dynamics and emission of relativistic jets. In order to achieve this objective I have been carrying it out in two parallel ways. In the first of them I have performed numerical magnetohydrodynamic and emission (eRMHD) simulations of relativistic jets. The simulations have been performed in collaboration with the Relativistic Astrophysics Group in the University of Valencia, using a numerical code that solves the RMHD equations in conservative form and cylindrical coordinates with axial symmetry (see Leismann et al., 2005, for more details). I have focused the study on the role played by the magnetic field in the dynamics of the jet, analyzing the balance of the main driving forces which determine the jet evolution. By using these relativistic magnetohydrodynamic (RMHD) simulations as input I have computed the non-thermal (synchrotron) emission which allows to obtain synthetic radio maps that can be directly compared with actual observations (Roca-Sogorb et al., 2008a,b, 2009). The synergy between simulations and observations is proven to be a powerful tool in the under- standing of the physical processes taking place in jets. For this, my second line of work has been the comparison of the eRMHD results with actual sources. I have started the study with a very well known source: the radiogalaxy 3C 120. I have carried out new observations, taken in November 2007 making use of all available observing frequencies (from 1 to 86 GHz) with the VLBA and VLA. These observations, taken also in polarimetric mode, allow to study the jet in 3C120 from pc to Kpc scales with great detail. The comparison of these observations with those from 1999 to 2001 (Gomez et al., o 2008) provides information about the source of Faraday rotation in the jet of the radio galaxy 3C 120. The results indicate that the emitting jet and the source of Faraday rotation are not closely connected physically and have different configurations for the magnetic field and/or kinematical properties, favoring a model in which a significant fraction of the RM originates in foreground clouds (Gomez et al., o 2011). The higher frequency 2007 observations reveal a new component located 80 mas from the core (which corresponds to a deprojected distance of 140 pc), with a brightness temperature about 600 times higher than expected at such distances. I have analyzed the different processes that could be responsible for the enhanced brightness temperature observed, its sudden appearance, and apparent 1 stationary (Roca-Sogorb et al., 2010)

A recollimation shock 80 mas from the core in the jet of the radio galaxy 3C120: Observational evidence and modeling

We present Very Long Baseline Array observations of the radio galaxy 3C120 at 5, 8, 12, and 15 GHz designed to study a peculiar stationary jet feature (hereafter C80) located ~80 mas from the core, which was previously shown to display a brightness temperature ~600 times lager than expected at such distances. The high sensitivity of the images -- obtained between December 2009 and June 2010 -- has revealed that C80 corresponds to the eastern flux density peak of an arc of emission (hereafter A80), downstream of which extends a large (~20 mas in size) bubble-like structure that resembles an inverted bow shock. The linearly polarized emission closely follows that of the total intensity in A80, with the electric vector position angle distributed nearly perpendicular to the arc-shaped structure. Despite the stationary nature of C80/A80, superluminal components with speeds up to ~3 c have been detected downstream from its position, resembling the behavior observed in the HST-1 emission complex in M87. The total and polarized emission of the C80/A80 structure, its lack of motion, and brightness temperature excess are best reproduced by a model based on synchrotron emission from a conical shock with cone opening angle \eta=10 degrees, jet viewing angle \theta=16 degrees, a completely tangled upstream magnetic field, and upstream Lorentz factor \gamma=8.4. The good agreement between our observations and numerical modeling leads us to conclude that the peculiar feature associated with C80/A80 corresponds to a conical recollimation shock in the jet of 3C120 located at a de-projected distance of ~190 pc downstream from the nucleus.Comment: Accepted for publication in Ap

APRICOT: Advanced Platform for Reproducible Infrastructures in the Cloud via Open Tools

[EN] Background Scientific publications are meant to exchange knowledge among researchers but the inability to properly reproduce computational experiments limits the quality of scientific research. Furthermore, bibliography shows that irreproducible preclinical research exceeds 50%, which produces a huge waste of resources on nonprofitable research at Life Sciences field. As a consequence, scientific reproducibility is being fostered to promote Open Science through open databases and software tools that are typically deployed on existing computational resources. However, some computational experiments require complex virtual infrastructures, such as elastic clusters of PCs, that can be dynamically provided from multiple clouds. Obtaining these infrastructures requires not only an infrastructure provider, but also advanced knowledge in the cloud computing field. Objectives The main aim of this paper is to improve reproducibility in life sciences to produce better and more cost-effective research. For that purpose, our intention is to simplify the infrastructure usage and deployment for researchers. Methods This paper introduces Advanced Platform for Reproducible Infrastructures in the Cloud via Open Tools (APRICOT), an open source extension for Jupyter to deploy deterministic virtual infrastructures across multiclouds for reproducible scientific computational experiments. To exemplify its utilization and how APRICOT can improve the reproduction of experiments with complex computation requirements, two examples in the field of life sciences are provided. All requirements to reproduce both experiments are disclosed within APRICOT and, therefore, can be reproduced by the users. Results To show the capabilities of APRICOT, we have processed a real magnetic resonance image to accurately characterize a prostate cancer using a Message Passing Interface cluster deployed automatically with APRICOT. In addition, the second example shows how APRICOT scales the deployed infrastructure, according to the workload, using a batch cluster. This example consists of a multiparametric study of a positron emission tomography image reconstruction. Conclusion APRICOT's benefits are the integration of specific infrastructure deployment, the management and usage for Open Science, making experiments that involve specific computational infrastructures reproducible. All the experiment steps and details can be documented at the same Jupyter notebook which includes infrastructure specifications, data storage, experimentation execution, results gathering, and infrastructure termination. Thus, distributing the experimentation notebook and needed data should be enough to reproduce the experiment.This study was supported by the program "Ayudas para la contratación de personal investigador en formación de carácter predoctoral, programa VALi+d" under grant number ACIF/2018/148 from the Conselleria d'Educació of the Generalitat Valenciana and the "Fondo Social Europeo" (FSE). The authors would like to thank the Spanish "Ministerio de Economía, Industria y Competitividad" for the project "BigCLOE" with reference number TIN2016-79951-R and the European Commission, Horizon 2020 grant agreement No 826494 (PRIMAGE). The MRI prostate study case used in this article has been retrospectively collected from a project of prostate MRI biomarkers validation.Giménez-Alventosa, V.; Segrelles Quilis, JD.; Moltó, G.; Roca-Sogorb, M. (2020). APRICOT: Advanced Platform for Reproducible Infrastructures in the Cloud via Open Tools. Methods of Information in Medicine. 59(S 02):e33-e45. https://doi.org/10.1055/s-0040-1712460Se33e4559S 02Donoho, D. L., Maleki, A., Rahman, I. U., Shahram, M., & Stodden, V. (2009). Reproducible Research in Computational Harmonic Analysis. Computing in Science & Engineering, 11(1), 8-18. doi:10.1109/mcse.2009.15Freedman, L. P., Cockburn, I. M., & Simcoe, T. S. (2015). The Economics of Reproducibility in Preclinical Research. PLOS Biology, 13(6), e1002165. doi:10.1371/journal.pbio.1002165Chillarón, M., Vidal, V., & Verdú, G. (2020). CT image reconstruction with SuiteSparseQR factorization package. Radiation Physics and Chemistry, 167, 108289. doi:10.1016/j.radphyschem.2019.04.039Reader, A. J., Ally, S., Bakatselos, F., Manavaki, R., Walledge, R. J., Jeavons, A. P., … Zweit, J. (2002). One-pass list-mode EM algorithm for high-resolution 3-D PET image reconstruction into large arrays. IEEE Transactions on Nuclear Science, 49(3), 693-699. doi:10.1109/tns.2002.1039550Giménez-Alventosa, V., Antunes, P. C. G., Vijande, J., Ballester, F., Pérez-Calatayud, J., & Andreo, P. (2016). Collision-kerma conversion between dose-to-tissue and dose-to-water by photon energy-fluence corrections in low-energy brachytherapy. Physics in Medicine and Biology, 62(1), 146-164. doi:10.1088/1361-6560/aa4f6aWilkinson, M. D., Dumontier, M., Aalbersberg, Ij. J., Appleton, G., Axton, M., Baak, A., … Bourne, P. E. (2016). The FAIR Guiding Principles for scientific data management and stewardship. Scientific Data, 3(1). doi:10.1038/sdata.2016.18Calatrava, A., Romero, E., Moltó, G., Caballer, M., & Alonso, J. M. (2016). Self-managed cost-efficient virtual elastic clusters on hybrid Cloud infrastructures. Future Generation Computer Systems, 61, 13-25. doi:10.1016/j.future.2016.01.018Caballer, M., Blanquer, I., Moltó, G., & de Alfonso, C. (2014). Dynamic Management of Virtual Infrastructures. Journal of Grid Computing, 13(1), 53-70. doi:10.1007/s10723-014-9296-5Wolstencroft, K., Owen, S., Krebs, O., Nguyen, Q., Stanford, N. J., Golebiewski, M., … Goble, C. (2015). SEEK: a systems biology data and model management platform. BMC Systems Biology, 9(1). doi:10.1186/s12918-015-0174-yDe Alfonso, C., Caballer, M., Calatrava, A., Moltó, G., & Blanquer, I. (2018). Multi-elastic Datacenters: Auto-scaled Virtual Clusters on Energy-Aware Physical Infrastructures. Journal of Grid Computing, 17(1), 191-204. doi:10.1007/s10723-018-9449-zRawla, P. (2019). Epidemiology of Prostate Cancer. World Journal of Oncology, 10(2), 63-89. doi:10.14740/wjon1191Bratan, F., Niaf, E., Melodelima, C., Chesnais, A. L., Souchon, R., Mège-Lechevallier, F., … Rouvière, O. (2013). Influence of imaging and histological factors on prostate cancer detection and localisation on multiparametric MRI: a prospective study. European Radiology, 23(7), 2019-2029. doi:10.1007/s00330-013-2795-0Le, J. D., Tan, N., Shkolyar, E., Lu, D. Y., Kwan, L., Marks, L. S., … Reiter, R. E. (2015). Multifocality and Prostate Cancer Detection by Multiparametric Magnetic Resonance Imaging: Correlation with Whole-mount Histopathology. European Urology, 67(3), 569-576. doi:10.1016/j.eururo.2014.08.079Brix, G., Semmler, W., Port, R., Schad, L. R., Layer, G., & Lorenz, W. J. (1991). Pharmacokinetic Parameters in CNS Gd-DTPA Enhanced MR Imaging. Journal of Computer Assisted Tomography, 15(4), 621-628. doi:10.1097/00004728-199107000-00018Larsson, H. B. W., Stubgaard, M., Frederiksen, J. L., Jensen, M., Henriksen, O., & Paulson, O. B. (1990). Quantitation of blood-brain barrier defect by magnetic resonance imaging and gadolinium-DTPA in patients with multiple sclerosis and brain tumors. Magnetic Resonance in Medicine, 16(1), 117-131. doi:10.1002/mrm.1910160111Tofts, P. S., & Kermode, A. G. (1991). Measurement of the blood-brain barrier permeability and leakage space using dynamic MR imaging. 1. Fundamental concepts. Magnetic Resonance in Medicine, 17(2), 357-367. doi:10.1002/mrm.1910170208Donahue, K. M., Weisskoff, R. M., & Burstein, D. (1997). Water diffusion and exchange as they influence contrast enhancement. Journal of Magnetic Resonance Imaging, 7(1), 102-110. doi:10.1002/jmri.1880070114Flouri, D., Lesnic, D., & Sourbron, S. P. (2015). Fitting the two-compartment model in DCE-MRI by linear inversion. Magnetic Resonance in Medicine, 76(3), 998-1006. doi:10.1002/mrm.25991Brun, R., & Rademakers, F. (1997). ROOT — An object oriented data analysis framework. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 389(1-2), 81-86. doi:10.1016/s0168-9002(97)00048-xXuan Liu, Comtat, C., Michel, C., Kinahan, P., Defrise, M., & Townsend, D. (2001). Comparison of 3-D reconstruction with 3D-OSEM and with FORE+OSEM for PET. IEEE Transactions on Medical Imaging, 20(8), 804-814. doi:10.1109/42.938248Singh, S., Kalra, M. K., Hsieh, J., Licato, P. E., Do, S., Pien, H. H., & Blake, M. A. (2010). Abdominal CT: Comparison of Adaptive Statistical Iterative and Filtered Back Projection Reconstruction Techniques. Radiology, 257(2), 373-383. doi:10.1148/radiol.10092212Shepp, L. A., & Vardi, Y. (1982). Maximum Likelihood Reconstruction for Emission Tomography. IEEE Transactions on Medical Imaging, 1(2), 113-122. doi:10.1109/tmi.1982.4307558Goo, J. M., Tongdee, T., Tongdee, R., Yeo, K., Hildebolt, C. F., & Bae, K. T. (2005). Volumetric Measurement of Synthetic Lung Nodules with Multi–Detector Row CT: Effect of Various Image Reconstruction Parameters and Segmentation Thresholds on Measurement Accuracy. Radiology, 235(3), 850-856. doi:10.1148/radiol.2353040737Ravenel, J. G., Leue, W. M., Nietert, P. J., Miller, J. V., Taylor, K. K., & Silvestri, G. A. (2008). Pulmonary Nodule Volume: Effects of Reconstruction Parameters on Automated Measurements—A Phantom Study. Radiology, 247(2), 400-408. doi:10.1148/radiol.2472070868Hu, Y.-H., Zhao, B., & Zhao, W. (2008). Image artifacts in digital breast tomosynthesis: Investigation of the effects of system geometry and reconstruction parameters using a linear system approach. Medical Physics, 35(12), 5242-5252. doi:10.1118/1.2996110Lyra, M., & Ploussi, A. (2011). Filtering in SPECT Image Reconstruction. International Journal of Biomedical Imaging, 2011, 1-14. doi:10.1155/2011/69379

Faraday rotation and polarization gradients in the jet of 3C~120: Interaction with the external medium and a helical magnetic field?

We present a sequence of 12 monthly polarimetric 15, 22, and 43 GHz VLBA observations of the radio galaxy 3C 120 revealing a systematic presence of gradients in Faraday rotation and degree of polarization across and along the jet. The degree of polarization increases with distance from the core and toward the jet edges, and has an asymmetric profile in which the northern side of the jet is more highly polarized. The Faraday rotation measure is also stratified across the jet width, with larger values for the southern side. We find a localized region of high Faraday rotation measure superposed on this structure between approximately 3 and 4 mas from the core, with a peak of about 6000 rad/m^2. Interaction of the jet with the external medium or a cloud would explain the confined region of enhanced Faraday rotation, as well as the stratification in degree of polarization and the flaring of superluminal knots when crossing this region. The data are also consistent with a helical field in a two-fluid jet model, consisting of an inner, emitting jet and a sheath containing nonrelativistic electrons. However, this helical magnetic field model cannot by itself explain the localized region of enhanced Faraday rotation. The polarization electric vectors, predominantly perpendicular to the jet axis once corrected for Faraday rotation, require a dominant component parallel to the jet axis (in the frame of the emitting plasma) for the magnetic field in the emitting region.Comment: Accepted for publication in ApJ Letters. 4 pages (including 5 figures

On the source of Faraday rotation in the jet of the radio galaxy 3C120

The source of Faraday rotation in the jet of the radio galaxy 3C120 is analyzed through Very Long Baseline Array observations carried out between 1999 and 2007 at 86, 43, 22, 15, 12, 8, 5, 2, and 1.7 GHz. Comparison of observations from 1999 to 2001 reveals uncorrelated changes in the linear polarization of the underlying jet emission and the Faraday rotation screen: while the rotation measure (RM) remains constant between approximately 2 and 5 mas from the core, the RM-corrected electric vector position angles (EVPAs) of two superluminal components are rotated by almost 90 degrees when compared to other components moving through similar jet locations. On the other hand, the innermost 2 mas experiences a significant change in RM -- including a sign reversal -- but without variations in the RM-corrected EVPAs. Similarly, observations in 2007 reveal a double sign reversal in RM along the jet, while the RM-corrected EVPAs remain perpendicular to the jet axis. Although the observed coherent structure and gradient of the RM along the jet supports the idea that the Faraday rotation is produced by a sheath of thermal electrons that surrounds the emitting jet, the uncorrelated changes in the RM and RM-corrected EVPAs indicate that the emitting jet and the source of Faraday rotation are not closely connected physically and have different configurations for the magnetic field and/or kinematical properties. Furthermore, the existence of a region of enhanced RM whose properties remain constant over three years requires a localized source of Faraday rotation, favoring a model in which a significant fraction of the RM originates in foreground clouds.Comment: 12 pages, 11 figures; Accepted for publication in Ap

Probing the Inner Jet of the Quasar PKS 1510-089 with Multi-waveband Monitoring during Strong Gamma-ray Activity

We present results from monitoring the multi-waveband flux, linear polarization, and parsec-scale structure of the quasar PKS 1510-089, concentrating on eight major gamma-ray flares that occurred during the interval 2009.0-2009.5. The gamma-ray peaks were essentially simultaneous with maxima at optical wavelengths, although the flux ratio of the two wavebands varied by an order of magnitude. The optical polarization vector rotated by 720 degrees during a 5-day period encompassing six of these flares. This culminated in a very bright, roughly 1 day, optical and gamma-ray flare as a bright knot of emission passed through the highest-intensity, stationary feature (the "core") seen in 43 GHz Very Long Baseline Array images. The knot continued to propagate down the jet at an apparent speed of 22c and emit strongly at gamma-ray energies as a months-long X-ray/radio outburst intensified. We interpret these events as the result of the knot following a spiral path through a mainly toroidal magnetic field pattern in the acceleration and collimation zone of the jet, after which it passes through a standing shock in the 43 GHz core and then continues downstream. In this picture, the rapid gamma-ray flares result from scattering of infrared seed photons from a relatively slow sheath of the jet as well as from optical synchrotron radiation in the faster spine. The 2006-2009.7 radio and X-ray flux variations are correlated at very high significance; we conclude that the X-rays are mainly from inverse Compton scattering of infrared seed photons by 20-40 MeV electrons.Comment: 10 pages of text + 5 figures, to be published in Astrophysical Journal Letters in 201

On the Location of the Gamma-ray Emission in the 2008 Outburst in the BL Lacertae Object AO 0235+164 through Observations across the Electromagnetic Spectrum

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with 0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7 mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary "core" and in the superluminal knot, both parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes.Comment: Accepted for Publication in the Astrophysical Journal Letters. 7 pages (including 5 figures). Minor corrections with regard to previous version, as proposed by the refere

Flaring Behavior of the Quasar 3C~454.3 across the Electromagnetic Spectrum

We analyze the behavior of the parsec-scale jet of the quasar 3C~454.3 during pronounced flaring activity in 2005-2008. Three major disturbances propagated down the jet along different trajectories with Lorentz factors $\Gamma>$10. The disturbances show a clear connection with millimeter-wave outbursts, in 2005 May/June, 2007 July, and 2007 December. High-amplitude optical events in the $R$-band light curve precede peaks of the millimeter-wave outbursts by 15-50 days. Each optical outburst is accompanied by an increase in X-ray activity. We associate the optical outbursts with propagation of the superluminal knots and derive the location of sites of energy dissipation in the form of radiation. The most prominent and long-lasting of these, in 2005 May, occurred closer to the black hole, while the outbursts with a shorter duration in 2005 Autumn and in 2007 might be connected with the passage of a disturbance through the millimeter-wave core of the jet. The optical outbursts, which coincide with the passage of superluminal radio knots through the core, are accompanied by systematic rotation of the position angle of optical linear polarization. Such rotation appears to be a common feature during the early stages of flares in blazars. We find correlations between optical variations and those at X-ray and $\gamma$-ray energies. We conclude that the emergence of a superluminal knot from the core yields a series of optical and high-energy outbursts, and that the mm-wave core lies at the end of the jet's acceleration and collimation zone.Comment: 57 pages, 23 figures, 8 tables (submitted to ApJ