Showing 1–50 of 111 results for author: Tassis, K

Systems design, assembly, integration and lab testing of WALOP-South Polarimeter

Authors: Siddharth Maharana, A. N. Ramaprakash, Chaitanya Rajarshi, Pravin Khodade, Bhushan Joshi, Pravin Chordia, Abhay Kohok, Ramya M. Anche, Deepa Modi, John A. Kypriotakis, Amit Deokar, Aditya Kinjawadekar, Stephen B. Potter, Dmitry Blinov, Hans Kristian Eriksen, Myrto Falalaki, Hitesh Gajjar, Tuhin Ghosh, Eirik Gjerløw, Sebastain Kiehlmann, Ioannis Liodakis, Nikolaos Mandarakas, Georgia V. Panopoulou, Vasiliki Pavlidou, Timothy J. Pearson , et al. (6 additional authors not shown)

Abstract: Wide-Area Linear Optical Polarimeter (WALOP)-South is the first wide-field and survey-capacity polarimeter in the optical wavelengths. On schedule for commissioning in 2024, it will be mounted on the 1 m SAAO telescope in Sutherland Observatory, South Africa to undertake the PASIPHAE sky survey. PASIPHAE program will create the first polarimetric sky map in the optical wavelengths, spanning more t… ▽ More Wide-Area Linear Optical Polarimeter (WALOP)-South is the first wide-field and survey-capacity polarimeter in the optical wavelengths. On schedule for commissioning in 2024, it will be mounted on the 1 m SAAO telescope in Sutherland Observatory, South Africa to undertake the PASIPHAE sky survey. PASIPHAE program will create the first polarimetric sky map in the optical wavelengths, spanning more than 2000 square degrees of the southern Galactic region. The innovative design of WALOP-South will enable it to measure the linear polarization (Stokes parameters $q$ and $u$), in a single exposure, of all sources in a field of view (FoV) of $35\times35$ arcminutes-squared in the SDSS-r broadband and narrowband filters between 500-750 nm with 0.1 % polarization accuracy. The unique goals of the instrument place very stringent systems engineering goals, including on the performance of the optical, polarimetric, optomechanical, and electronic subsystems. All the subsystems have been designed carefully to meet the overall instrument performance goals. As of May 2024, all the instrument optical and mechanical subsystems have been assembled and are currently getting tested and integrated. The complete testing and characterization of the instrument in the lab is expected to be completed by August 2024. In this paper, we will present (a) the design and development of the entire instrument and its major subsystems, focusing on the opto-mechanical design which has not been reported before, and (b) assembly and integration of the instrument in the lab and early results from lab characterization of the instrument. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 20 pages; presented at 2024 SPIE Astronomical Telescopes + Instrumentation conference

The LiteBIRD mission to explore cosmic inflation

Authors: T. Ghigna, A. Adler, K. Aizawa, H. Akamatsu, R. Akizawa, E. Allys, A. Anand, J. Aumont, J. Austermann, S. Azzoni, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, A. Basyrov, S. Beckman, M. Bersanelli, M. Bortolami, F. Bouchet, T. Brinckmann, P. Campeti, E. Carinos, A. Carones , et al. (134 additional authors not shown)

Abstract: LiteBIRD, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-… ▽ More LiteBIRD, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-year mission, LiteBIRD will employ three telescopes within 15 unique frequency bands (ranging from 34 through 448 GHz), targeting a sensitivity of 2.2\,$μ$K-arcmin and a resolution of 0.5$^\circ$ at 100\,GHz. Its primary goal is to measure the tensor-to-scalar ratio $r$ with an uncertainty $δr = 0.001$, including systematic errors and margin. If $r \geq 0.01$, LiteBIRD expects to achieve a $>5σ$ detection in the $\ell=$2-10 and $\ell=$11-200 ranges separately, providing crucial insight into the early Universe. We describe LiteBIRD's scientific objectives, the application of systems engineering to mission requirements, the anticipated scientific impact, and the operations and scanning strategies vital to minimizing systematic effects. We will also highlight LiteBIRD's synergies with concurrent CMB projects. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 23 pages, 9 figures, 1 table, SPIE Astronomical Telescopes + Instrumentation 2024

The first degree-scale starlight-polarization-based tomography map of the magnetized interstellar medium

Authors: V. Pelgrims, N. Mandarakas, R. Skalidis, K. Tassis, G. V. Panopoulou, V. Pavlidou, D. Blinov, S. Kiehlmann, S. E. Clark, B. S. Hensley, S. Romanopoulos, A. Basyrov, H. K. Eriksen, M. Falalaki, T. Ghosh, E. Gjerløw, J. A. Kypriotakis, S. Maharana, A. Papadaki, T. J. Pearson, S. B. Potter, A. N. Ramaprakash, A. C. S. Readhead, I. K. Wehus

Abstract: We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of starlight polarization in a region of the sky of 4 square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to i… ▽ More We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of starlight polarization in a region of the sky of 4 square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to invert the 3D volume occupied by the observed stars. We used it to obtain the first 3D map of the dusty magnetized ISM. Specifically, we produced a tomography map of the orientation of the plane-of-sky (POS) component of the magnetic field threading the diffuse, dusty regions responsible for the stellar polarization. For the targeted region centered on Galactic coordinates $(l,b) \approx (103.3^\circ, 22.3^\circ)$, we identified several ISM clouds. Most of the lines of sight intersect more than one cloud. A very nearby component was detected in the foreground of a dominant component from which most of the polarization signal comes. Farther clouds, with a distance of up to 2~kpc, were similarly detected. Some of them likely correspond to intermediate-velocity clouds seen in HI spectra in this region of the sky. We found that the orientation of the POS component of the magnetic field changes along distance for most of the lines of sight. Our study demonstrates that starlight polarization data coupled to distance measures have the power to reveal the great complexity of the dusty magnetized ISM in 3D and, in particular, to provide local measurements of the POS component of the magnetic field. This demonstrates that the inversion of large data volumes, as expected from the PASIPHAE survey, will provide the necessary means to move forward in the modeling of the Galactic magnetic field and of the dusty magnetized ISM as a contaminant in observations of the cosmic microwave background polarization. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Preprint of the accepted version at Astronomy & Astrophysics. The 3D map obtained in this paper can be visualized online at https://pasiphae.science/visualization and a video featuring it is accessible at https://www.youtube.com/watch?v=dB_6J1zhmPI

LiteBIRD Science Goals and Forecasts: Primordial Magnetic Fields

Authors: D. Paoletti, J. Rubino-Martin, M. Shiraishi, D. Molinari, J. Chluba, F. Finelli, C. Baccigalupi, J. Errard, A. Gruppuso, A. I. Lonappan, A. Tartari, E. Allys, A. Anand, J. Aumont, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, M. Bersanelli, M. Bortolami, T. Brinckmann, E. Calabrese, P. Campeti, A. Carones, F. J. Casas , et al. (75 additional authors not shown)

Abstract: We present detailed forecasts for the constraints on primordial magnetic fields (PMFs) that will be obtained with the LiteBIRD satellite. The constraints are driven by the effects of PMFs on the CMB anisotropies: the gravitational effects of magnetically-induced perturbations; the effects on the thermal and ionization history of the Universe; the Faraday rotation imprint on the CMB polarization; a… ▽ More We present detailed forecasts for the constraints on primordial magnetic fields (PMFs) that will be obtained with the LiteBIRD satellite. The constraints are driven by the effects of PMFs on the CMB anisotropies: the gravitational effects of magnetically-induced perturbations; the effects on the thermal and ionization history of the Universe; the Faraday rotation imprint on the CMB polarization; and the non-Gaussianities induced in polarization anisotropies. LiteBIRD represents a sensitive probe for PMFs and by exploiting all the physical effects, it will be able to improve the current limit coming from Planck. In particular, thanks to its accurate $B$-mode polarization measurement, LiteBIRD will improve the constraints on infrared configurations for the gravitational effect, giving $B_{\rm 1\,Mpc}^{n_{\rm B} =-2.9} < 0.8$ nG at 95% C.L., potentially opening the possibility to detect nanogauss fields with high significance. We also observe a significant improvement in the limits when marginalized over the spectral index, $B_{1\,{\rm Mpc}}^{\rm marg}< 2.2$ nG at 95% C.L. From the thermal history effect, which relies mainly on $E$-mode polarization data, we obtain a significant improvement for all PMF configurations, with the marginalized case, $\sqrt{\langle B^2\rangle}^{\rm marg}<0.50$ nG at 95% C.L. Faraday rotation constraints will take advantage of the wide frequency coverage of LiteBIRD and the high sensitivity in $B$ modes, improving the limits by orders of magnitude with respect to current results, $B_{1\,{\rm Mpc}}^{n_{\rm B} =-2.9} < 3.2$ nG at 95% C.L. Finally, non-Gaussianities of the $B$-mode polarization can probe PMFs at the level of 1 nG, again significantly improving the current bounds from Planck. Altogether our forecasts represent a broad collection of complementary probes, providing conservative limits on PMF characteristics that will be achieved with LiteBIRD. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 51 pages, 24 figures, abstract shortened

Time evolution of the galactic $B- ρ$ relation: the impact of the magnetic field morphology

Authors: A. Konstantinou, E. Ntormousi, K. Tassis, A. Pallottini

Abstract: One of the most frequently used indicators to characterize the magnetic field's influence on star formation is the relation between magnetic field strength and gas density ($B-ρ$ relation), usually expressed as $B \propto ρ^κ$. The value of $κ$ is an indication of the dynamical importance of the magnetic field during gas compression. Investigating the global magnetic field's impact on this relatio… ▽ More One of the most frequently used indicators to characterize the magnetic field's influence on star formation is the relation between magnetic field strength and gas density ($B-ρ$ relation), usually expressed as $B \propto ρ^κ$. The value of $κ$ is an indication of the dynamical importance of the magnetic field during gas compression. Investigating the global magnetic field's impact on this relation and its evolution, we conduct MHD simulations of Milky-Way-like galaxies including gravity, star formation, and supernova feedback along with non-equilibrium chemistry up to $H_2$ formation fueling star formation. Two initial magnetic field morphologies are studied: one completely ordered (toroidal) and the other completely random. In these models, we study the dynamical importance of the magnetic field through the plasma $β$ and the $B-ρ$ relation. For both magnetic morphologies, low-density regions are thermally supported, while high-density regions are magnetically dominated. Equipartition is reached earlier and at lower densities in the toroidal model. However, the $B-ρ$ relation is not unique even within the same galaxy, as it consistently includes two different branches for a given density, with $κ$ ranging from about 0.2 to 0.8. The mean value of $κ$ for each model also displays significant variations over time, which supersede the differences between the two models. While our findings suggest that the magnetic field morphology does influence the galactic $B-ρ$ relation, its impact is transient, since time-averaged differences between the models fall within the large temporal scatter. The context and time-dependent nature of the $B-ρ$ relation underscore the need for comprehensive research and observations to understand the intricate role of magnetic fields in star formation processes across diverse galactic environments. △ Less

Submitted 19 February, 2024; v1 submitted 15 February, 2024; originally announced February 2024.

Comments: The paper has been accepted for publication at A&A

Journal ref: A&A 686, A8 (2024)

Impact of beam far side-lobe knowledge in the presence of foregrounds for LiteBIRD

Authors: C. Leloup, G. Patanchon, J. Errard, C. Franceschet, J. E. Gudmundsson, S. Henrot-Versillé, H. Imada, H. Ishino, T. Matsumura, G. Puglisi, W. Wang, A. Adler, J. Aumont, R. Aurlien, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, A. Basyrov, M. Bersanelli, D. Blinov, M. Bortolami, T. Brinckmann, P. Campeti , et al. (86 additional authors not shown)

Abstract: We present a study of the impact of an uncertainty in the beam far side-lobe knowledge on the measurement of the Cosmic Microwave Background $B$-mode signal at large scale. It is expected to be one of the main source of systematic effects in future CMB observations. Because it is crucial for all-sky survey missions to take into account the interplays between beam systematic effects and all the dat… ▽ More We present a study of the impact of an uncertainty in the beam far side-lobe knowledge on the measurement of the Cosmic Microwave Background $B$-mode signal at large scale. It is expected to be one of the main source of systematic effects in future CMB observations. Because it is crucial for all-sky survey missions to take into account the interplays between beam systematic effects and all the data analysis steps, the primary goal of this paper is to provide the methodology to carry out the end-to-end study of their effect for a space-borne CMB polarization experiment, up to the cosmological results in the form of a bias $δr$ on the tensor-to-scalar ratio $r$. LiteBIRD is dedicated to target the measurement of CMB primordial $B$ modes by reaching a sensitivity of $σ\left( r \right) \leq 10^{-3}$ assuming $r=0$. As a demonstration of our framework, we derive the relationship between the knowledge of the beam far side-lobes and the tentatively allocated error budget under given assumptions on design, simulation and component separation method. We assume no mitigation of the far side-lobes effect at any stage of the analysis pipeline. We show that $δr$ is mostly due to the integrated fractional power difference between the estimated beams and the true beams in the far side-lobes region, with little dependence on the actual shape of the beams, for low enough $δr$. Under our set of assumptions, in particular considering the specific foreground cleaning method we used, we find that the integrated fractional power in the far side-lobes should be known at a level as tight as $\sim 10^{-4}$, to achieve the required limit on the bias $δr < 1.9 \times 10^{-5}$. The framework and tools developed for this study can be easily adapted to provide requirements under different design, data analysis frameworks and for other future space-borne experiments beyond LiteBIRD. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Zero-polarization candidate regions for calibration of wide-field optical polarimeters

Authors: N. Mandarakas, G. V. Panopoulou, V. Pelgrims, S. B. Potter, V. Pavlidou, A. Ramaprakash, K. Tassis, D. Blinov, S. Kiehlmann, E. Koutsiona, S. Maharana, S. Romanopoulos, R. Skalidis, A. Vervelaki, S. E. Clark, J. A. Kypriotakis, A. C. S. Readhead

Abstract: Context. Calibration of optical polarimeters relies on the use of stars with negligible polarization (unpolarized standard stars) for determining the instrumental polarization zero-point. For wide-field polarimeters, calibration is often done by imaging the same star over multiple positions in the field of view - a process which is time-consuming. A more effective technique is to target fields con… ▽ More Context. Calibration of optical polarimeters relies on the use of stars with negligible polarization (unpolarized standard stars) for determining the instrumental polarization zero-point. For wide-field polarimeters, calibration is often done by imaging the same star over multiple positions in the field of view - a process which is time-consuming. A more effective technique is to target fields containing multiple standard stars. While this method has been used for fields with highly polarized stars, there are no such sky regions with well-measured unpolarized standard stars. Aims. We aim to identify sky regions with tens of stars exhibiting negligible polarization, which are suitable for zero-point calibration of wide-field polarimeters. Methods. We selected stars in regions with extremely low reddening, located at high Galactic latitudes. We targeted four ~ 400 x 400 fields in the northern, and eight in the southern Equatorial hemisphere. Observations were carried out at the Skinakas Observatory and the South African Astronomical Observatory respectively. Results. We find two fields in the North and seven in the South with mean polarization lower than p < 0.1%. Conclusions. At least nine out of twelve fields can be used for zero-point calibration of wide-field polarimeters. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 15 pages, 25 figures

LiteBIRD Science Goals and Forecasts: Improving Sensitivity to Inflationary Gravitational Waves with Multitracer Delensing

Authors: T. Namikawa, A. I. Lonappan, C. Baccigalupi, N. Bartolo, D. Beck, K. Benabed, A. Challinor, P. Diego-Palazuelos, J. Errard, S. Farrens, A. Gruppuso, N. Krachmalnicoff, M. Migliaccio, E. Martínez-González, V. Pettorino, G. Piccirilli, M. Ruiz-Granda, B. Sherwin, J. Starck, P. Vielva, R. Akizawa, A. Anand, J. Aumont, R. Aurlien, S. Azzoni , et al. (97 additional authors not shown)

Abstract: We estimate the efficiency of mitigating the lensing $B$-mode polarization, the so-called delensing, for the $LiteBIRD$ experiment with multiple external data sets of lensing-mass tracers. The current best bound on the tensor-to-scalar ratio, $r$, is limited by lensing rather than Galactic foregrounds. Delensing will be a critical step to improve sensitivity to $r$ as measurements of $r$ become mo… ▽ More We estimate the efficiency of mitigating the lensing $B$-mode polarization, the so-called delensing, for the $LiteBIRD$ experiment with multiple external data sets of lensing-mass tracers. The current best bound on the tensor-to-scalar ratio, $r$, is limited by lensing rather than Galactic foregrounds. Delensing will be a critical step to improve sensitivity to $r$ as measurements of $r$ become more and more limited by lensing. In this paper, we extend the analysis of the recent $LiteBIRD$ forecast paper to include multiple mass tracers, i.e., the CMB lensing maps from $LiteBIRD$ and CMB-S4-like experiment, cosmic infrared background, and galaxy number density from $Euclid$- and LSST-like survey. We find that multi-tracer delensing will further improve the constraint on $r$ by about $20\%$. In $LiteBIRD$, the residual Galactic foregrounds also significantly contribute to uncertainties of the $B$-modes, and delensing becomes more important if the residual foregrounds are further reduced by an improved component separation method. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: 21 pages, 7 figures

LiteBIRD Science Goals and Forecasts: A full-sky measurement of gravitational lensing of the CMB

Authors: A. I. Lonappan, T. Namikawa, G. Piccirilli, P. Diego-Palazuelos, M. Ruiz-Granda, M. Migliaccio, C. Baccigalupi, N. Bartolo, D. Beck, K. Benabed, A. Challinor, J. Errard, S. Farrens, A. Gruppuso, N. Krachmalnicoff, E. Martínez-González, V. Pettorino, B. Sherwin, J. Starck, P. Vielva, R. Akizawa, A. Anand, J. Aumont, R. Aurlien, S. Azzoni , et al. (97 additional authors not shown)

Abstract: We explore the capability of measuring lensing signals in $LiteBIRD$ full-sky polarization maps. With a $30$ arcmin beam width and an impressively low polarization noise of $2.16\,μ$K-arcmin, $LiteBIRD$ will be able to measure the full-sky polarization of the cosmic microwave background (CMB) very precisely. This unique sensitivity also enables the reconstruction of a nearly full-sky lensing map u… ▽ More We explore the capability of measuring lensing signals in $LiteBIRD$ full-sky polarization maps. With a $30$ arcmin beam width and an impressively low polarization noise of $2.16\,μ$K-arcmin, $LiteBIRD$ will be able to measure the full-sky polarization of the cosmic microwave background (CMB) very precisely. This unique sensitivity also enables the reconstruction of a nearly full-sky lensing map using only polarization data, even considering its limited capability to capture small-scale CMB anisotropies. In this paper, we investigate the ability to construct a full-sky lensing measurement in the presence of Galactic foregrounds, finding that several possible biases from Galactic foregrounds should be negligible after component separation by harmonic-space internal linear combination. We find that the signal-to-noise ratio of the lensing is approximately $40$ using only polarization data measured over $90\%$ of the sky. This achievement is comparable to $Planck$'s recent lensing measurement with both temperature and polarization and represents a four-fold improvement over $Planck$'s polarization-only lensing measurement. The $LiteBIRD$ lensing map will complement the $Planck$ lensing map and provide several opportunities for cross-correlation science, especially in the northern hemisphere. △ Less

Submitted 8 December, 2023; originally announced December 2023.

LiteBIRD Science Goals and Forecasts. A Case Study of the Origin of Primordial Gravitational Waves using Large-Scale CMB Polarization

Authors: P. Campeti, E. Komatsu, C. Baccigalupi, M. Ballardini, N. Bartolo, A. Carones, J. Errard, F. Finelli, R. Flauger, S. Galli, G. Galloni, S. Giardiello, M. Hazumi, S. Henrot-Versillé, L. T. Hergt, K. Kohri, C. Leloup, J. Lesgourgues, J. Macias-Perez, E. Martínez-González, S. Matarrese, T. Matsumura, L. Montier, T. Namikawa, D. Paoletti , et al. (85 additional authors not shown)

Abstract: We study the possibility of using the $LiteBIRD$ satellite $B$-mode survey to constrain models of inflation producing specific features in CMB angular power spectra. We explore a particular model example, i.e. spectator axion-SU(2) gauge field inflation. This model can source parity-violating gravitational waves from the amplification of gauge field fluctuations driven by a pseudoscalar "axionlike… ▽ More We study the possibility of using the $LiteBIRD$ satellite $B$-mode survey to constrain models of inflation producing specific features in CMB angular power spectra. We explore a particular model example, i.e. spectator axion-SU(2) gauge field inflation. This model can source parity-violating gravitational waves from the amplification of gauge field fluctuations driven by a pseudoscalar "axionlike" field, rolling for a few e-folds during inflation. The sourced gravitational waves can exceed the vacuum contribution at reionization bump scales by about an order of magnitude and can be comparable to the vacuum contribution at recombination bump scales. We argue that a satellite mission with full sky coverage and access to the reionization bump scales is necessary to understand the origin of the primordial gravitational wave signal and distinguish among two production mechanisms: quantum vacuum fluctuations of spacetime and matter sources during inflation. We present the expected constraints on model parameters from $LiteBIRD$ satellite simulations, which complement and expand previous studies in the literature. We find that $LiteBIRD$ will be able to exclude with high significance standard single-field slow-roll models, such as the Starobinsky model, if the true model is the axion-SU(2) model with a feature at CMB scales. We further investigate the possibility of using the parity-violating signature of the model, such as the $TB$ and $EB$ angular power spectra, to disentangle it from the standard single-field slow-roll scenario. We find that most of the discriminating power of $LiteBIRD$ will reside in $BB$ angular power spectra rather than in $TB$ and $EB$ correlations. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: 22 pages, 13 figures. Submitted to JCAP

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

Authors: Sergio Martin-Alvarez, Enrique Lopez-Rodriguez, Tara Dacunha, Susan E. Clark, Alejandro S. Borlaff, Rainer Beck, Francisco Rodríguez Montero, S. Lyla Jung, Julien Devriendt, Adrianne Slyz, Julia Roman-Duval, Evangelia Ntormousi, Mehrnoosh Tahani, Kandaswamy Subramanian, Daniel A. Dale, Pamela M. Marcum, Konstantinos Tassis, Ignacio del Moral-Castro, Le Ngoc Tram, Matt J. Jarvis

Abstract: The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy t… ▽ More The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of $h_\text{I radio} \sim h_\text{PI radio} = 0.23 \pm 0.03$ kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane ($h_\text{I FIR} \sim h_\text{PI FIR} = 0.11 \pm 0.01$ kpc). Radio emission traces magnetic fields at scales of $\gtrsim 300$ pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution ($<300$ pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism. △ Less

Submitted 24 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the ApJ. Replaced to match final version. 35 pages, 17 figures, 4 tables

The RoboPol sample of optical polarimetric standards

Authors: D. Blinov, S. Maharana, F. Bouzelou, C. Casadio, E. Gjerløw, J. Jormanainen, S. Kiehlmann, J. A. Kypriotakis, I. Liodakis, N. Mandarakas, L. Markopoulioti, G. V. Panopoulou, V. Pelgrims, A. Pouliasi, S. Romanopoulos, R. Skalidis, R. M. Anche, E. Angelakis, J. Antoniadis, B. J. Medhi, T. Hovatta, A. Kus, N. Kylafis, A. Mahabal, I. Myserlis , et al. (12 additional authors not shown)

Abstract: Optical polarimeters are typically calibrated using measurements of stars with known and stable polarization parameters. However, there is a lack of such stars available across the sky. Many of the currently available standards are not suitable for medium and large telescopes due to their high brightness. Moreover, as we find, some of the used polarimetric standards are in fact variable or have po… ▽ More Optical polarimeters are typically calibrated using measurements of stars with known and stable polarization parameters. However, there is a lack of such stars available across the sky. Many of the currently available standards are not suitable for medium and large telescopes due to their high brightness. Moreover, as we find, some of the used polarimetric standards are in fact variable or have polarization parameters that differ from their cataloged values. Our goal is to establish a sample of stable standards suitable for calibrating linear optical polarimeters with an accuracy down to $10^{-3}$ in fractional polarization. For five years, we have been running a monitoring campaign of a sample of standard candidates comprised of 107 stars distributed across the northern sky. We analyzed the variability of the linear polarization of these stars, taking into account the non-Gaussian nature of fractional polarization measurements. For a subsample of nine stars, we also performed multiband polarization measurements. We created a new catalog of 65 stars (see Table 2) that are stable, have small uncertainties of measured polarimetric parameters, and can be used as calibrators of polarimeters at medium- and large-size telescopes. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Journal ref: A&A 677, A144 (2023)

Bright-Moon Sky as a Wide-Field Linear Polarimetric Flat Source for Calibration

Authors: S. Maharana, S. Kiehlmann, D. Blinov, V. Pelgrims, V. Pavlidou, K. Tassis, J. A. Kypriotakis, A. N. Ramaprakash, R. M. Anche, A. Basyrov, K. Deka, H. K. Eriksen, T. Ghosh, E. Gjerløw, N. Mandarakas, E. Ntormousi, G. V. Panopoulou, A. Papadaki, T. Pearson, S. B. Potter, A. C. S. Readhead, R. Skalidis, I. K. Wehus

Abstract: Next-generation wide-field optical polarimeters like the Wide-Area Linear Optical Polarimeters (WALOPs) have a field of view (FoV) of tens of arcminutes. For efficient and accurate calibration of these instruments, wide-field polarimetric flat sources will be essential. Currently, no established wide-field polarimetric standard or flat sources exist. This paper tests the feasibility of using the p… ▽ More Next-generation wide-field optical polarimeters like the Wide-Area Linear Optical Polarimeters (WALOPs) have a field of view (FoV) of tens of arcminutes. For efficient and accurate calibration of these instruments, wide-field polarimetric flat sources will be essential. Currently, no established wide-field polarimetric standard or flat sources exist. This paper tests the feasibility of using the polarized sky patches of the size of around ten-by-ten arcminutes, at a distance of up to 20 degrees from the Moon, on bright-Moon nights as a wide-field linear polarimetric flat source. We observed 19 patches of the sky adjacent to the bright-Moon with the RoboPol instrument in the SDSS-r broadband filter. These were observed on five nights within two days of the full-Moon across two RoboPol observing seasons. We find that for 18 of the 19 patches, the uniformity in the measured normalized Stokes parameters $q$ and $u$ is within 0.2 %, with 12 patches exhibiting uniformity within 0.07 % or better for both $q$ and $u$ simultaneously, making them reliable and stable wide-field linear polarization flats. We demonstrate that the sky on bright-Moon nights is an excellent wide-field linear polarization flat source. Various combinations of the normalized Stokes parameters $q$ and $u$ can be obtained by choosing suitable locations of the sky patch with respect to the Moon △ Less

Submitted 7 May, 2023; originally announced May 2023.

Comments: 8 pages including appendix, 6 figures and 3 tables. Submitted to Astronomy and Astrophysics for review. Comments are welcome

Journal ref: A&A 679, A68 (2023)

Turnaround density evolution encodes cosmology in simulations

Authors: Giorgos Korkidis, Vasiliki Pavlidou, Konstantinos Tassis

Abstract: The mean matter density within the turnaround radius, which is the boundary that separates a nonexpanding structure from the Hubble flow, was recently proposed as a novel cosmological probe. According to the spherical collapse model, the evolution with cosmic time of this turnaround density, $\rm ρ_{ta}(z)$, can be used to determine both $\rm Ω_m$ and $Ω_Λ$, independently of any other currently us… ▽ More The mean matter density within the turnaround radius, which is the boundary that separates a nonexpanding structure from the Hubble flow, was recently proposed as a novel cosmological probe. According to the spherical collapse model, the evolution with cosmic time of this turnaround density, $\rm ρ_{ta}(z)$, can be used to determine both $\rm Ω_m$ and $Ω_Λ$, independently of any other currently used probe. The properties of $\rm ρ_{ta}$ predicted by the spherical collapse model were also shown to persist in the presence of full three-dimensional effects in $\rm Λ$CDM N-body cosmological simulations when considering galaxy clusters at the present time, $z=0$. However, a small offset was discovered between the spherical-collapse prediction of the value of $ρ_{ta}$ at $z=0$ and its value measured in simulations. In this letter, we explore whether this offset evolves with cosmic time; whether it differs in different cosmologies; whether its origin can be confidently identified; and whether it can be corrected. We found that the offset does evolve slightly with redshift, and that it correlates strongly with the deviation from spherical symmetry of the dark matter halo distribution inside and outside of the turnaround radius. We used an appropriate metric to quantify deviations in the environment of a structure from spherical symmetry. We found that using this metric, we can construct a sphericity-selected sample of halos for which the offset of $ρ_{ta}$ from the spherical collapse prediction is zero, independently of redshift and cosmology. We found that a sphericity-selected halo sample allows us to recover the simulated cosmology, and we conclude that the turnaround density evolution indeed encodes the cosmology in N-body simulations. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: 8 pages, 4 figures, accepted for publication in A&A

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- V: First results on the magnetic field orientation of galaxies

Authors: Alejandro S. Borlaff, Enrique Lopez-Rodriguez, Rainer Beck, Susan E. Clark, Evangelia Ntormousi, Konstantinos Tassis, Sergio Martin-Alvarez, Mehrnoosh Tahani, Daniel A. Dale, Ignacio del Moral Castro, Julia Roman-Duval, Pamela M. Marcum, John E. Beckman, Kandaswamy Subramanian, Sarah Eftekharzadeh, Leslie Proudfit

Abstract: We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from… ▽ More We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from $53$ to $214$ \um. We compute the magnetic pitch angle ($Ψ_{B}$) profiles as a function of the galactrocentric radius. We introduce a new magnetic alignment parameter ($ζ$) to estimate the disordered-to-ordered $B$-field ratio in spiral $B$-fields. We find FIR and radio wavelengths to not generally trace the same $B$-field morphology in galaxies. The $Ψ_{B}$ profiles tend to be more ordered with galactocentric radius in radio ($ζ_{\rm{6cm}} = 0.93\pm0.03$) than in FIR ($ζ_{\rm{154μm}} = 0.84\pm0.14$). For spiral galaxies, FIR $B$-fields are $2-75$\% more turbulent than the radio $B$-fields. For starburst galaxies, we find that FIR polarization is a better tracer of the $B$-fields along the galactic outflows than radio polarization. Our results suggest that the $B$-fields associated with dense, dusty, turbulent star-forming regions, those traced at FIR, are less ordered than warmer, less-dense regions, those traced at radio, of the interstellar medium. The FIR $B$-fields seem to be more sensitive to the activity of the star-forming regions and the morphology of the molecular clouds within a vertical height of few hundred pc in the disk of spiral galaxies than the radio $B$-fields. △ Less

Submitted 13 June, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: 28 pages, 14 figures. Accepted for publication in The Astrophysical Journal

Compact Symmetric Objects -- III Evolution of the High-Luminosity Branch and a Possible Connection with Tidal Disruption Events

Authors: A. C. S. Readhead, V. Ravi, R. D. Blandford, A. G. Sullivan, J. Somalwar, M. C. Begelman, M. Birkinshaw, I. Liodakis, M. L. Lister, T. J. Pearson, G. B. Taylor, P. N. Wilkinson, N. Globus, S. Kiehlmann, C. R. Lawrence, D. Murphy, S. O'Neill, V. Pavlidou, E. Sheldahl, A. Siemiginowska, K. Tassis

Abstract: We use a sample of 54 Compact Symmetric Objects (CSOs) to confirm that there are two unrelated CSO classes: an edge-dimmed, low-luminosity class (CSO~1), and an edge-brightened, high-luminosity class (CSO~2). Using blind tests, we show that CSO~2s consist of three sub-classes: CSO 2.0, having prominent hot-spots at the leading edges of narrow jets and/or narrow lobes; CSO~2.2, without prominent ho… ▽ More We use a sample of 54 Compact Symmetric Objects (CSOs) to confirm that there are two unrelated CSO classes: an edge-dimmed, low-luminosity class (CSO~1), and an edge-brightened, high-luminosity class (CSO~2). Using blind tests, we show that CSO~2s consist of three sub-classes: CSO 2.0, having prominent hot-spots at the leading edges of narrow jets and/or narrow lobes; CSO~2.2, without prominent hot-spots, and with broad jets and/or lobes; and CSO~2.1, which exhibit mixed properties. Most CSO 2s do not evolve into larger jetted-AGN, but spend their whole life-cycle as CSOs of size $\lesssim$500 pc and age $\lesssim$5000 yr. The minimum energies needed to produce the radio luminosity and structure in CSO~2s range from $\sim~10^{-4}\,M_\odot{c}^2$ to $\sim7\,M_\odot{c}^2$. We show that the transient nature of most CSO~2s, and their birthrate, can be explained through ignition in the tidal disruption events of giant stars. We also consider possibilities of tapping the spin energy of the supermassive black hole, and tapping the energy of the accretion disk. Our results demonstrate that CSOs constitute a large family of AGN in which we have thus far studied only the brightest. More comprehensive CSO studies, with higher sensitivity, resolution, and dynamic range, will revolutionize our understanding of AGN and the central engines that power them. △ Less

Submitted 26 November, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: 44 pages, 16 figures, 9 tables, accepted for publication

Compact Symmetric Objects -- II Confirmation of a Distinct Population of High-Luminosity Jetted Active Galaxies

Authors: S. Kiehlmann, A. C. S. Readhead, S. O'Neill, P. N. Wilkinson, M. L. Lister, I. Liodakis, S. Bruzewski, V. Pavlidou, T. J. Pearson, E. Sheldahl, A. Siemiginowska, K. Tassis, G. B. Taylor

Abstract: Compact Symmetric Objects (CSOs) are compact (<1 kpc), jetted Active Galactic Nuclei (AGN), whose jet axes are not aligned close to the line of sight, and whose observed emission is not predominantly relativistically boosted towards us. Two classes of CSOs have previously been identified: approximately one fifth are edge-dimmed and designated as CSO 1s, while the rest are edge brightened and desig… ▽ More Compact Symmetric Objects (CSOs) are compact (<1 kpc), jetted Active Galactic Nuclei (AGN), whose jet axes are not aligned close to the line of sight, and whose observed emission is not predominantly relativistically boosted towards us. Two classes of CSOs have previously been identified: approximately one fifth are edge-dimmed and designated as CSO 1s, while the rest are edge brightened and designated as CSO 2s. This paper focuses almost exclusively on CSO 2s. Using complete samples of CSO 2s we present three independent lines of evidence, based on their relative numbers, redshift distributions, and size distributions, which show conclusively that the vast majority (> 99%) of CSO 2s do not evolve into larger-scale radio sources. These CSO 2s belong to a distinct population of jetted-AGN, which should be characterized as ``short-lived'' compared to the classes of larger jetted-AGN, as opposed to ``young''. We show that there is a sharp upper cutoff in the CSO 2 size distribution at $\approx 500$ pc. The distinct differences between most CSO 2s and other jetted-AGN provides a crucial new time domain window on the formation and evolution of relativistic jets in AGN and the supermassive black holes that drive them. △ Less

Submitted 26 November, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: 29 pages, 10 figures, 7 tables, accepted for publication

Compact Symmetric Objects -- I Towards a Comprehensive Bona Fide Catalog

Authors: S. Kiehlmann, M. L. Lister, A. C. S. Readhead, I. Liodakis, S. O'Neill, T. J. Pearson, E. Sheldahl, A. Siemiginowska, K. Tassis, G. B. Taylor, P. N. Wilkinson

Abstract: Compact Symmetric Objects (CSOs) are jetted Active Galactic Nuclei (AGN) with overall projected size <1 kpc. The classification was introduced to distinguish these objects from the majority of compact jetted-AGN in centimeter wavelength very long baseline interferometry observations, where the observed emission is relativistically boosted towards the observer. The original classification criteria… ▽ More Compact Symmetric Objects (CSOs) are jetted Active Galactic Nuclei (AGN) with overall projected size <1 kpc. The classification was introduced to distinguish these objects from the majority of compact jetted-AGN in centimeter wavelength very long baseline interferometry observations, where the observed emission is relativistically boosted towards the observer. The original classification criteria for CSOs were: (i) evidence of emission on both sides of the center of activity, and (ii) overall size <1 kpc. However some relativistically boosted objects with jet axes close to the line of sight appear symmetric and have been mis-classified as CSOs, thereby undermining the CSO classification. This is because two essential CSO properties, pointed out in the original papers, have been neglected: (iii) low variability, and (iv) low apparent speeds along the jets. As a first step towards creating a comprehensive catalog of ``bona fide'' CSOs, we identify 79 bona fide CSOs, including 15 objects claimed as confirmed CSOs here for the first time, that match the CSO selection criteria. This sample of bona fide CSOs can be used for astrophysical studies of CSOs without contamination by mis-classified CSOs. We show that the fraction of CSOs in complete flux density limited AGN samples with S$_{\rm 5\,GHz}$ >700 mJy is between $(6.8\pm1.6)$% and $(8.5\pm1.8)$%. △ Less

Submitted 26 November, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: 28 pages, 9 figures, 3 tables, accepted for publication

CO enhancement by magnetohydrodynamic waves; Striations in the Polaris Flare

Authors: R. Skalidis, K. Gkimisi, K. Tassis, G. V. Panopoulou, V. Pelgrims, A. Tritsis, P. F. Goldsmith

Abstract: The formation of molecular gas in interstellar clouds is a slow process, but is enhanced by gas compression. Magnetohydrodynamic (MHD) waves create compressed quasiperiodic linear structures, referred to as striations. Striations are observed at column densities where the atomic to molecular gas transition takes place. We explore the role of MHD waves in the CO chemistry in regions with striations… ▽ More The formation of molecular gas in interstellar clouds is a slow process, but is enhanced by gas compression. Magnetohydrodynamic (MHD) waves create compressed quasiperiodic linear structures, referred to as striations. Striations are observed at column densities where the atomic to molecular gas transition takes place. We explore the role of MHD waves in the CO chemistry in regions with striations within molecular clouds. We target a region with striations in the Polaris Flare cloud. We conduct a CO J=2-1 survey in order to probe the molecular gas properties. We use archival starlight polarization data and dust emission maps in order to probe the magnetic field properties and compare against the CO properties. We assess the interaction of compressible MHD wave modes with CO chemistry by comparing their characteristic timescales. The estimated magnetic field is 38 - 76 $μ$G. In the CO integrated intensity map, we observe a dominant quasi-periodic intensity structure, which tends to be parallel to the magnetic field orientation and has a wavelength of one parsec approximately. The periodicity axis is $\sim$ 17 degrees off from the mean magnetic field orientation and is also observed in the dust intensity map. The contrast in the CO integrated intensity map is $\sim 2.4$ times larger than the contrast of the column density map, indicating that CO formation is enhanced locally. We suggest that a dominant slow magnetosonic mode with estimated period $2.1 - 3.4$ Myr, and propagation speed $0.30 - 0.45$ km~s$^{-1}$, is likely to have enhanced the formation of CO, hence created the observed periodic pattern. We also suggest that, within uncertainties, a fast magnetosonic mode with period 0.48 Myr and velocity $2.0$ km~s$^{-1}$ could have played some role in increasing the CO abundance. Quasiperiodic CO structures observed in striation regions may be the imprint of MHD wave modes. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: Accepted for publication in A&A, 16 pages, 13 figures

Journal ref: A&A 673, A76 (2023)

Extragalactic magnetism with SOFIA (SALSA Legacy Program). VI. The magnetic fields in the multi-phase interstellar medium of the Antennae galaxies

Authors: Enrique Lopez-Rodriguez, Alejandro S. Borlaff, Rainer Beck, William T. Reach, Sui Ann Mao, Evangelia Ntormousi, Konstantinos Tassis, Sergio Martin-Alvarez, Susan E. Clark, Daniel A. Dale, Ignacio del Moral-Castro

Abstract: Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multi-phase ISM of the closest merger of two spiral galaxies, the Antennae galaxie… ▽ More Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multi-phase ISM of the closest merger of two spiral galaxies, the Antennae galaxies. We compare the inferred B-fields using $154~μ$m thermal dust and $11$ cm radio synchrotron emission polarimetric observations. We find that the $154~μ$m B-fields are more ordered across the Antennae galaxies than the $11$ cm B-fields. The turbulent-to-ordered $154~μ$m B-field increases at the galaxy cores and star-forming regions. The relic spiral arm has an ordered spiral $154~μ$m B-field, while the $11$ cm B-field is radial. The $154~μ$m B-field may be dominated by turbulent dynamos with high $^{12}$CO(1-0) velocity dispersion driven by star-forming regions, while the $11$ cm B-field is cospatial with high HI velocity dispersion driven by galaxy interaction. This result shows the dissociation between the warm gas mainly disturbed by the merger, and the dense gas still following the dynamics of the relic spiral arm. We find a $\sim8.9$ kpc scale ordered B-field connecting the two galaxies. The base of the tidal tail is cospatial with the HI and $^{12}$CO(1-0) emission and has compressed and/or sheared $154~μ$m and $11$ cm B-fields driven by the merger. We suggest that amplify B-fields, with respect to the rest of the system and other spiral galaxies, may be supporting the gas flow between both galaxies and the tidal tail. △ Less

Submitted 9 December, 2022; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: 11 pages, 5 figures, Accepted for publication in ApJ Letters

Analytic characterization of sub-Alfvénic turbulence energetics

Authors: R. Skalidis, K. Tassis, V. Pavlidou

Abstract: Magnetohydrodynamic (MHD) turbulence is a cross-field process relevant to many systems. A prerequisite for understanding these systems is to constrain the role of MHD turbulence, and in particular the energy exchange between kinetic and magnetic forms. The energetics of strongly magnetized and compressible turbulence has so far resisted attempts to understand them. Numerical simulations reveal tha… ▽ More Magnetohydrodynamic (MHD) turbulence is a cross-field process relevant to many systems. A prerequisite for understanding these systems is to constrain the role of MHD turbulence, and in particular the energy exchange between kinetic and magnetic forms. The energetics of strongly magnetized and compressible turbulence has so far resisted attempts to understand them. Numerical simulations reveal that kinetic energy can be orders of magnitude larger than fluctuating magnetic energy. We solve this lack-of-balance puzzle by calculating the energetics of compressible and sub-Alfvénic turbulence based on the dynamics of coherent cylindrical fluid parcels. Using the MHD Lagrangian, we prove analytically that the bulk of the magnetic energy transferred to kinetic is the energy stored in the coupling between the ordered and fluctuating magnetic field. The analytical relations are in striking agreement with numerical data, up to second order terms. △ Less

Submitted 6 April, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Journal ref: A&A 672, L3 (2023)

Discriminating power of milli-lensing observations for dark matter models

Authors: Nick Loudas, Vasiliki Pavlidou, Carolina Casadio, Kostas Tassis

Abstract: The nature of dark matter (DM) is still under intense debate. Sub-galactic scales are particularly critical, as different, currently viable DM models make diverse predictions on the expected abundance and density profile of DM haloes on these scales. We investigate the ability of sub-galactic DM haloes to act as strong lenses on background compact sources, producing gravitational lensing events on… ▽ More The nature of dark matter (DM) is still under intense debate. Sub-galactic scales are particularly critical, as different, currently viable DM models make diverse predictions on the expected abundance and density profile of DM haloes on these scales. We investigate the ability of sub-galactic DM haloes to act as strong lenses on background compact sources, producing gravitational lensing events on milli-arcsecond scales (milli-lenses), for different DM models. For each DM scenario, we explore whether a sample of $\sim$ 5000 distant sources is sufficient to detect at least one milli-lens. We develop a semi-analytical model to estimate the milli-lensing optical depth as a function of the source's redshift for various DM models. We employ the Press-Schechter formalism, as well as results from recent N-body simulations to compute the halo mass function, taking into account the appropriate spherically averaged density profile of haloes for each DM model. We treat the lensing system as a point-mass lens and invoke the effective surface mass density threshold to calculate the fraction of a halo that acts as a gravitational lens. We study three classes of dark matter models: cold DM, warm DM, and self-interacting DM. We find that haloes consisting of warm DM turn out to be optically thin for strong gravitational milli-lensing (zero expected lensing events). CDM haloes may produce lensing events depending on the steepness of the concentration-mass relation. Self-interacting DM haloes can efficiently act as gravitational milli-lenses only if haloes experience gravothermal collapse, resulting in highly dense central cores. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Comments: 11 pages, 6 figures. Submitted for publication in A&A

Journal ref: A&A 668, A166 (2022)

GRB 210619B optical afterglow polarization

Authors: N. Mandarakas, D. Blinov, D. R. Aguilera-Dena, S. Romanopoulos, V. Pavlidou, K. Tassis, J. Antoniadis, S. Kiehlmann, A. Lychoudis, L. F. Tsemperof Kataivatis

Abstract: We report on the follow-up of the extremely bright long gamma-ray burst GRB~210619B with optical polarimetry. We conducted optopolarimetric observations of the optical afterglow of GRB~210619B in the SDSS-r band in the time window ~ 5967 - 8245 seconds after the burst, using the RoboPol instrument at the Skinakas observatory. We report a $5\,σ$ detection of polarization $P=1.5\pm0.3$ at polarizati… ▽ More We report on the follow-up of the extremely bright long gamma-ray burst GRB~210619B with optical polarimetry. We conducted optopolarimetric observations of the optical afterglow of GRB~210619B in the SDSS-r band in the time window ~ 5967 - 8245 seconds after the burst, using the RoboPol instrument at the Skinakas observatory. We report a $5\,σ$ detection of polarization $P=1.5\pm0.3$ at polarization angle $EVPA=8\pm6^\circ$. We find that during our observations the polarization is likely constant. These values are corrected for polarization induced by the interstellar medium of the Milky Way and host-induced polarization is likely negligible. Thus the polarization we quote is intrinsic to the GRB afterglow. △ Less

Submitted 11 December, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: 6 pages, 4 figures

Journal ref: A&A 670, A144 (2023)

WALOP-South: A Four-Camera One-Shot Imaging Polarimeter for PASIPHAE Survey. Paper II -- Polarimetric Modelling and Calibration

Authors: Siddharth Maharana, Ramya M. Anche, A. N. Ramaprakash, Bhushan Joshi, Artem Basyrov, Dmitry Blinov, Carolina Casadio, Kishan Deka, Hans Kristian Eriksen, Tuhin Ghosh, Eirik Gjerløw, John A. Kypriotakis, Sebastian Kiehlmann, Nikolaos Mandarakas, Georgia V. Panopoulou, Katerina Papadaki, Vasiliki Pavlidou, Timothy J. Pearson, Vincent Pelgrims, Stephen B. Potter, Anthony C. S. Readhead, Raphael Skalidis, Trygve Leithe Svalheim, Konstantinos Tassis, Ingunn K. Wehus

Abstract: The Wide-Area Linear Optical Polarimeter (WALOP)-South instrument is an upcoming wide-field and high-accuracy optical polarimeter to be used as a survey instrument for carrying out the Polar-Areas Stellar Imaging in Polarization High Accuracy Experiment (PASIPHAE) program. Designed to operate as a one-shot four-channel and four-camera imaging polarimeter, it will have a field of view of… ▽ More The Wide-Area Linear Optical Polarimeter (WALOP)-South instrument is an upcoming wide-field and high-accuracy optical polarimeter to be used as a survey instrument for carrying out the Polar-Areas Stellar Imaging in Polarization High Accuracy Experiment (PASIPHAE) program. Designed to operate as a one-shot four-channel and four-camera imaging polarimeter, it will have a field of view of $35\times 35$ arcminutes and will measure the Stokes parameters $I$, $q$, and $u$ in a single exposure in the SDSS-r broadband filter. The design goal for the instrument is to achieve an overall polarimetric measurement accuracy of 0.1 % over the entire field of view. We present here the complete polarimetric modeling of the instrument, characterizing the amount and sources of instrumental polarization. To accurately retrieve the real Stokes parameters of a source from the measured values, we have developed a calibration method for the instrument. Using this calibration method and simulated data, we demonstrate how to correct instrumental polarization and obtain 0.1 % accuracy in the degree of polarization, $p$. Additionally, we tested and validated the calibration method by implementing it on a table-top WALOP-like test-bed polarimeter in the laboratory. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: 34 pages, 24 figures. Accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems

Journal ref: J. of Astronomical Telescopes, Instruments, and Systems, 8(3), 038004 (2022)

Starlight-polarization-based tomography of the magnetized interstellar medium: PASIPHAE's line-of-sight inversion method

Authors: V. Pelgrims, G. V. Panopoulou, K. Tassis, V. Pavlidou, A. Basyrov, D. Blinov, E. Gjerløw, S. Kiehlmann, N. Mandarakas, A. Papadaki, R. Skalidis, A. Tsouros, R. M. Anche, H. K. Eriksen, T. Ghosh, J. A. Kypriotakis, S. Maharana, E. Ntormousi, T. J. Pearson, S. B. Potter, A. N. Ramaprakash, A. C. S. Readhead, I. K. Wehus

Abstract: We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in 3D within dusty regions. We develop a model in which the polarization signal from the m… ▽ More We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in 3D within dusty regions. We develop a model in which the polarization signal from the magnetized and dusty ISM is described by thin layers at various distances. Our modeling makes it possible to infer the mean polarization (amplitude and orientation) induced by individual dusty clouds and to account for the turbulence-induced scatter in a generic way. We present a likelihood function that explicitly accounts for uncertainties in polarization and parallax. We develop a framework for reconstructing the magnetized ISM through the maximization of the log-likelihood using a nested sampling method. We test our Bayesian inversion method on mock data taking into account realistic uncertainties from Gaia and as expected for the optical polarization survey PASIPHAE according to the currently planned observing strategy. We demonstrate that our method is effective at recovering the cloud properties as soon as the polarization induced by a cloud to its background stars is higher than $\sim 0.1\%$ for the adopted survey exposure time and level of systematic uncertainty. Our method makes it possible to recover not only the mean polarization properties but also to characterize the intrinsic scatter, thus creating new ways to characterize ISM turbulence and the magnetic field strength. Finally, we apply our method to an existing data set of starlight polarization with known line-of-sight decomposition, demonstrating agreement with previous results and an improved quantification of uncertainties in cloud properties. △ Less

Submitted 28 February, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: Closely matches the published version: A&A 670, A164, 2023 The accompanying code (BISP-1) is available at https://github.com/vpelgrims/Bisp_1

Journal ref: A&A 670, A164 (2023)

The Musca molecular cloud: The perfect "filament" is still a sheet

Authors: A. Tritsis, F. Bouzelou, R. Skalidis, K. Tassis, T. Enßlin, G. Edenhofer

Abstract: The true 3-dimensional (3D) morphology of the Musca molecular cloud is a topic that has received significant attention lately. Given that Musca does not exhibit intense star-formation activity, unveiling its shape has the potential of also revealing crucial information regarding the physics that dictates the formation of the first generation of stars within molecular clouds. Here, we revisit the s… ▽ More The true 3-dimensional (3D) morphology of the Musca molecular cloud is a topic that has received significant attention lately. Given that Musca does not exhibit intense star-formation activity, unveiling its shape has the potential of also revealing crucial information regarding the physics that dictates the formation of the first generation of stars within molecular clouds. Here, we revisit the shape of Musca and we present a comprehensive array of evidence pointing towards a shape that is extended along the line-of-sight dimension: (a) 3D maps of differential extinction; (b) new non-local thermodynamic equilibrium radiative transfer simulations of CO rotational transitions from a sheet-like, magnetically-dominated simulated cloud; (c) an effective/critical density analysis of available CO observations; (d) indirect consequences that a filamentary structure would have had, from a theoretical star-formation perspective. We conclude that the full collection of observational evidence strongly suggests that Musca has a sheet-like geometry. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: 13 pages, 14 figures. Accepted for publication in MNRAS

Dancing with the stars: Stirring up extraordinary turbulence in Galactic center clouds

Authors: Konstantinos Tassis, Vasiliki Pavlidou

Abstract: Molecular clouds in the central molecular zone (CMZ) have been observed to feature turbulent line widths that are significantly higher, and scale with cloud size more steeply, than in the rest of the Milky Way. In the same Galactic region, the stellar density is also much higher than in the rest of the Milky Way, and the vertical stellar velocity dispersion is large, meaning that even young stars… ▽ More Molecular clouds in the central molecular zone (CMZ) have been observed to feature turbulent line widths that are significantly higher, and scale with cloud size more steeply, than in the rest of the Milky Way. In the same Galactic region, the stellar density is also much higher than in the rest of the Milky Way, and the vertical stellar velocity dispersion is large, meaning that even young stars are likely to cross the entire vertical extent of the CMZ within their lifetimes. Here, we investigate whether interactions of CMZ molecular clouds with crossing stars can account for the extraordinary properties of observed turbulence in this part of the Galaxy. We calculated the rate of energy deposition by stars crossing CMZ clouds due to (a) stellar winds and (b) dynamical friction, and compared it to the rate of turbulence decay. We calculated the predicted scaling of turbulence line width with cloud size in each case. We find that energy deposition by stellar winds of crossing massive stars can account for both the level and the scaling of CMZ cloud turbulence with cloud size. We also find that the mechanism stops being effective at a Galactocentric distance comparable to the CMZ extent. On the other hand, we find that dynamical friction by crossing stars does not constitute a significant driver of turbulence for CMZ clouds. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Comments: 3 pages, comments welcome, accepted for publication in Astronomy and Astrophysics Letters

Journal ref: A&A 662, L1 (2022)

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- IV: Program overview and first results on the polarization fraction

Authors: Enrique Lopez-Rodriguez, Sui Ann Mao, Rainer Beck, Alejandro S. Borlaff, Evangelia Ntormousi, Konstantinos Tassis, Daniel A. Dale, Julia Roman-Duval, Kandaswamy Subramanian, Sergio Martin-Alvarez, Pamela M. Marcum, Susan E. Clark, William T. Reach, Doyal A. Harper, Ellen G. Zweibel

Abstract: We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation an… ▽ More We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154μm} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $μ$m. The falling $53-154$ $μ$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy. △ Less

Submitted 7 July, 2022; v1 submitted 2 May, 2022; originally announced May 2022.

Comments: 43 pages, 28 figures. Accepted for publication to ApJ

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- III: First data release and on-the-fly polarization mapping characterization

Authors: Enrique Lopez-Rodriguez, Melanie Clarke, Sachin Shenoy, William Vacca, Simon Coude, Ryan Arneson, Peter Ashton, Sarah Eftekharzadeh, Rainer Beck, John E. Beckman, Alejandro S. Borlaff, Susan E. Clark, Daniel A. Dale, Sergio Martin-Alvarez, Evangelia Ntormousi, William T. Reach, Julia Roman-Duval, Konstantinos Tassis, Doyal A. Harper, Pamela M. Marcum

Abstract: We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipe… ▽ More We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipeline steps to obtain homogeneously reduced high-level data products of polarimetric maps of galaxies for use in scientific analysis. Our approach has a general design and can be applied to sources smaller than the field-of-view of the HAWC+ array in any given band. We estimate that the OTFMAP polarimetric mode offers a reduction of observing overheads by a factor 2.34, and an improvement in sensitivity by a factor 1.80 when compared to previously obtained polarimetric observations using the chopping and nodding mode. The OTFMAP is a significant optimization of the polarimetric mode of HAWC+ as it ultimately reduces the cost of operations of SOFIA/HAWC+ by increasing the science collected per hour of observation up to an overall factor of 2.49. The OTFMAP polarimetric mode is the standard observing strategy of SALSA. The results and quantitative analysis of this first data release are presented in Papers IV and V of the series. △ Less

Submitted 28 April, 2022; originally announced April 2022.

Comments: 26 pages, 16 figures, 5 tables. Submitted to ApJ (Comments welcome)

Polarization power spectra and dust cloud morphology

Authors: A. Konstantinou, V. Pelgrims, F. Fuchs, K. Tassis

Abstract: In the framework of studies of the CMB polarization and its Galactic foregrounds, the angular power spectra of thermal dust polarization maps have revealed an intriguing E/B asymmetry and a positive TE correlation. In interpretation studies of these observations, magnetized ISM dust clouds have been treated as filamentary structures only; however, sheet-like shapes are also supported by observatio… ▽ More In the framework of studies of the CMB polarization and its Galactic foregrounds, the angular power spectra of thermal dust polarization maps have revealed an intriguing E/B asymmetry and a positive TE correlation. In interpretation studies of these observations, magnetized ISM dust clouds have been treated as filamentary structures only; however, sheet-like shapes are also supported by observational and theoretical evidence. In this work, we study the influence of cloud shape and its connection to the local magnetic field on angular power spectra of thermal dust polarization maps. We simulate realistic filament-like and sheet-like interstellar clouds, and generate synthetic maps of their thermal dust polarized emission using the software $Asterion$. We compute their polarization power spectra in multipole range $\ell \in [100,500]$ and quantify the E/B power asymmetry through the $R_{EB}$ ratio, and the correlation coefficient $r^{TE}$ between T and E modes. We quantify the dependence of $R_{EB}$ and $r^{TE}$ values on the offset angle (between longest cloud axis and magnetic field) and inclination angle (between line-of-sight and magnetic field) for both cloud shapes embedded either in a regular or a turbulent magnetic field. We find that both cloud shapes cover the same regions of the ($R_{EB}$, $r^{TE}$) parameter space. The dependence on inclination and offset angles are similar for both shapes although sheet-like structures generally show larger scatter. In addition to the known dependence on the offset angle, we find a strong dependence of $R_{EB}$ and $r^{TE}$ on the inclination angle. The fact that filament-like and sheet-like structures may lead to polarization power spectra with similar ($R_{EB}$, $r^{TE}$) values complicates their interpretation. In future analyses, this degeneracy should be accounted for as well as the connection to the magnetic field geometry. △ Less

Submitted 27 April, 2022; originally announced April 2022.

Comments: 16 pages, 24 figures. Accepted for publication by A&A

Journal ref: A&A 663, A175 (2022)

The effect of cosmic variance on the characteristics of dust polarization power spectra

Authors: V. Pelgrims, E. Ntormousi, K. Tassis

Abstract: In the context of cosmic microwave background polarization studies and the characterization of the Galactic foregrounds, the power spectrum analysis of the thermal dust polarization sky has led to intriguing evidence of an E/B asymmetry and a positive TE correlation. In this work, we produce synthesized dust polarization maps from a set of global magneto-hydrodynamic (MHD) simulations of Milky-Way… ▽ More In the context of cosmic microwave background polarization studies and the characterization of the Galactic foregrounds, the power spectrum analysis of the thermal dust polarization sky has led to intriguing evidence of an E/B asymmetry and a positive TE correlation. In this work, we produce synthesized dust polarization maps from a set of global magneto-hydrodynamic (MHD) simulations of Milky-Way-sized galaxies, and analyze their power spectra at intermediate angular scales (angular multipoles $\ell \in \left[60 ,\, 140\right]$). We study the role of the initial configuration of the large-scale magnetic field, its strength, and the feedback on the power spectrum characteristics. Using full-galaxy MHD simulations, we were able to estimate the variance induced by the peculiar location of the observer in the galaxy. We find that the polarization power spectra sensitively depend on the observer's location, impeding a distinction between different simulation setups. There is a clear statistical difference between the power spectra measured from within the spiral arms and those measured from the inter-arm regions. Also, power spectra from within supernova-driven bubbles share common characteristics, regardless of the underlying model. However, no correlation was found between the properties of the polarization power spectra and the local (with respect to the observer) mean values of physical quantities such as the density and the strength of the magnetic field. Finally, we find indications that the global strength of the magnetic field may play a role in shaping the power spectrum characteristics; as the global magnetic field strength increases, the E/B asymmetry and the TE correlation increase, whereas the viewpoint-induced variance decreases. However, we find no direct correlation with the strength of the local magnetic field that permeates the mapped region of the interstellar medium. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 20 pages including 3 appendices. Accepted version at A&A

Journal ref: A&A 658, A134 (2022)

Non-ideal MHD simulations of subcritical prestellar cores with non-equilibrium chemistry

Authors: Aris Tritsis, Christoph Federrath, Karen Willacy, Konstantinos Tassis

Abstract: Non-ideal magnetohydrodynamic (MHD) effects are thought to be gravity's closest ally in overcoming the support of magnetic fields and in forming stars. Here, we modify the publicly available version of the adaptive mesh refinement code FLASH (Fryxell et al. 2000; Dubey et al. 2008) to include a detailed treatment of non-ideal MHD and study such effects in collapsing prestellar cores. We implement… ▽ More Non-ideal magnetohydrodynamic (MHD) effects are thought to be gravity's closest ally in overcoming the support of magnetic fields and in forming stars. Here, we modify the publicly available version of the adaptive mesh refinement code FLASH (Fryxell et al. 2000; Dubey et al. 2008) to include a detailed treatment of non-ideal MHD and study such effects in collapsing prestellar cores. We implement two very extended non-equilibrium chemical networks, the largest of which is comprised of $\sim$ 300 species and includes a detailed description of deuterium chemistry. The ambipolar-diffusion, Ohmic and Hall resistivities are then self-consistently calculated from the abundances of charged species. We present a series of 2-dimensional axisymmetric simulations where we vary the chemical model, cosmic-ray ionization rate, and grain distribution. We benchmark our implementation against ideal MHD simulations and previously-published results. We show that, at high densities ($n_{\rm{H_2}}>~10^6~\rm{cm^{-3}}$), the ion that carries most of the perpendicular and parallel conductivities is not $\rm{H_3^+}$ as was previously thought, but is instead $\rm{D_3^+}$. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Comments: 18 pages, 14 figures. Accepted for publication in MNRAS

HI-H$_2$ transition: exploring the role of the magnetic field

Authors: R. Skalidis, K. Tassis, G. V. Panopoulou, J. L. Pineda, Y. Gong, N. Mandarakas, D. Blinov, S. Kiehlmann, J. A. Kypriotakis

Abstract: Atomic gas in the diffuse interstellar medium (ISM) is organized in filamentary structures. These structures usually host cold and dense molecular clumps. The Galactic magnetic field is considered to play an important role in the formation of these clumps. Our goal is to explore the role of the magnetic field in the HI - H$_{2}$ transition process. We targeted a filamentary cloud where gas transit… ▽ More Atomic gas in the diffuse interstellar medium (ISM) is organized in filamentary structures. These structures usually host cold and dense molecular clumps. The Galactic magnetic field is considered to play an important role in the formation of these clumps. Our goal is to explore the role of the magnetic field in the HI - H$_{2}$ transition process. We targeted a filamentary cloud where gas transitions from atomic to molecular. This cloud is located at the edges of an expanding structure, known as the North Celestial Pole Loop (NCPL). We probed the magnetic field properties of the cloud with optical polarization observations. We performed multi-wavelength spectroscopic observations of different species in order to probe the gas phase properties of the cloud. We identified two distinct sub-regions within the cloud. One of the regions hosts purely atomic gas, while the other is dominated by molecular gas although most of it is CO-dark. The estimated plane-of-the-sky magnetic field strength between the two regions remains constant within uncertainties and lies in the range 20 ~ 30$~μ$G. The total magnetic field strength does not scale with density which implies that gas is compressed along the field lines. We also found that turbulence is sub-Alfvénic. The HI velocity gradients are in general perpendicular to the mean magnetic field orientation, except for the region close to the CO clump where they tend to become parallel. The latter is likely related to gas undergoing gravitational infall. The magnetic field morphology of the target cloud is parallel to the HI column density structure of the cloud in the atomic region, while it tends to become perpendicular to the HI structure in the molecular region. If this is verified in more cases it has important consequences for the ISM magnetic field modeling with HI data. △ Less

Submitted 27 May, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 665, A77 (2022)

Why take the square root? An assessment of interstellar magnetic field strength estimation methods

Authors: R. Skalidis, J. Sternberg, J. R. Beattie, V. Pavlidou, K. Tassis

Abstract: The magnetic field strength in interstellar clouds can be estimated indirectly by using the spread of dust polarization angles ($δθ$). The method developed by Davis 1951 and by Chandrasekhar and Fermi 1953 (DCF) assumes that incompressible magnetohydrodynamic (MHD) fluctuations induce the observed dispersion of polarization angles, deriving $B\propto 1/δθ$ (or, $δθ\propto M_{A}$, in terms of the A… ▽ More The magnetic field strength in interstellar clouds can be estimated indirectly by using the spread of dust polarization angles ($δθ$). The method developed by Davis 1951 and by Chandrasekhar and Fermi 1953 (DCF) assumes that incompressible magnetohydrodynamic (MHD) fluctuations induce the observed dispersion of polarization angles, deriving $B\propto 1/δθ$ (or, $δθ\propto M_{A}$, in terms of the Alfvénic Mach number). However, observations show that the interstellar medium (ISM) is highly compressible. Recently, Skalidis & Tassis 2021 (ST) relaxed the incompressibility assumption and derived instead $B\propto 1/\sqrt{δθ}$ ($δθ\propto M_{A}^2$). We explored what the correct scaling is in compressible and magnetized turbulence with numerical simulations. We used 26 magnetized, ideal-MHD numerical simulations with different types of forcing. The range of $M_{A}$ and sonic Mach numbers $M_{s}$ explored are $0.1 \leq M_{A} \leq 2.0$ and $0.5 \leq M_{s} \leq 20$. We created synthetic polarization maps and tested the assumptions and accuracy of the two methods. The synthetic data have a remarkable consistency with the $δθ\propto M_{A}^{2}$ scaling, which is inferred by ST, while the DCF scaling fails to follow the data. The ST method shows an accuracy better than $50\%$ over the entire range of $M_{A}$ explored; DCF performs adequately only in the range of $M_{A}$ for which it has been optimized through the use of a "fudge factor". For low $M_{A}$, DCF is inaccurate by factors of tens. The assumptions of the ST method reflect better the physical reality in clouds with compressible and magnetized turbulence, and for this reason the method provides a much better estimate of the magnetic field strength over the DCF method. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 656, A118 (2021)

Extragalactic magnetism with SOFIA (Legacy Program) -- II: A magnetically-driven flow in the starburst ring of NGC 1097

Authors: Enrique Lopez-Rodriguez, Rainer Beck, Susan E. Clark, Annie Hughes, Alejandro S. Borlaff, Evangelia Ntormousi, Lucas Grosset, Konstantinos Tassis, John E. Beckman, Kandaswamy Subramanian, Daniel Dale, Tanio Diaz-Santos, Legacy Team

Abstract: Galactic bars are frequent in disk galaxies and they may support the transfer of matter towards the central engine of active nuclei. The barred galaxy NGC 1097 has magnetic forces controlling the gas flow at several kpc scales, which suggest that magnetic fields (B-fields) are dynamically important along the bar and nuclear ring. However, the effect of the B-field on the gas flows in the central k… ▽ More Galactic bars are frequent in disk galaxies and they may support the transfer of matter towards the central engine of active nuclei. The barred galaxy NGC 1097 has magnetic forces controlling the gas flow at several kpc scales, which suggest that magnetic fields (B-fields) are dynamically important along the bar and nuclear ring. However, the effect of the B-field on the gas flows in the central kpc scale has not been characterized. Using thermal polarized emission at $89$ $μ$m with HAWC+/SOFIA, here, we measure that the polarized flux is spatially located at the contact regions of the outer-bar with the starburst ring. The linear polarization decomposition analysis shows that the $89$ $μ$m and radio ($3.5$ and $6.2$ cm) polarization traces two different modes, $m$, of the B-field: a constant B-field orientation and dominated by $m=0$ at $89$ $μ$m, and a spiral B-field dominated by $m=2$ at radio. We show that the B-field at 89 $μ$m is concentrated in the warmest region of a shock driven by the galactic-bar dynamics in the contact regions between the outer-bar with the starburst ring. Radio polarization traces a superposition of the spiral B-field outside and within the starburst ring. According to Faraday rotation measures between $3.5$ and $6.2$ cm, the radial component of the B-field along the contact regions points toward the galaxy's center on both sides. We conclude that gas streams outside and within the starburst ring follow the B-field, which feeds the black hole with matter from the host galaxy. △ Less

Submitted 7 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 17 pages, 10 figures, Accepted for publication to ApJ

SMILE: Search for MIlli-LEnses

Authors: C. Casadio, D. Blinov, A. C. S. Readhead, I. W. A. Browne, P. N. Wilkinson, T. Hovatta, N. Mandarakas, V. Pavlidou, K. Tassis, H. K. Vedantham, J. A. Zensus, V. Diamantopoulos, K. E. Dolapsaki, K. Gkimisi, G. Kalaitzidakis, M. Mastorakis, K. Nikolaou, E. Ntormousi, V. Pelgrims, K. Psarras

Abstract: Dark Matter (DM) halos with masses below $\sim10^{8}$ $M_{\odot}$, which would help to discriminate between DM models, may be detected through their gravitational effect on distant sources. The same applies to primordial black holes, considered as an alternative scenario to DM particle models. However, there is still no evidence for the existence of such objects. With the aim of finding compact ob… ▽ More Dark Matter (DM) halos with masses below $\sim10^{8}$ $M_{\odot}$, which would help to discriminate between DM models, may be detected through their gravitational effect on distant sources. The same applies to primordial black holes, considered as an alternative scenario to DM particle models. However, there is still no evidence for the existence of such objects. With the aim of finding compact objects in the mass range $\sim$ 10$^{6}$ -- 10$^{9}$$M_{\odot}$, we search for strong gravitational lenses on milli (mas)-arcseconds scales (< 150 mas). For our search, we used the Astrogeo VLBI FITS image database -- the largest publicly available database, containing multi-frequency VLBI data of 13828 individual sources. We used the citizen science approach to visually inspect all sources in all available frequencies in search for images with multiple compact components on mas-scales. At the final stage, sources were excluded based on the surface brightness preservation criterion. We obtained a sample of 40 sources that passed all steps and therefore are judged to be milli-arcsecond lens candidates. These sources are currently followed-up with on-going European VLBI Network (EVN) observations at 5 and 22 GHz. Based on spectral index measurements, we suggest that two of our candidates have a higher probability to be associated with gravitational lenses. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 6 pages, 1 figure; accepted for publication

Local alignments of parsec-scale AGN radiojets

Authors: N. Mandarakas, D. Blinov, C. Casadio, V. Pelgrims, S. Kiehlmann, V. Pavlidou, K. Tassis

Abstract: Context.Coherence in the characteristics of neighboring sources in 2D and 3D space may suggest the existence of large-scale cosmic structures, which are useful for cosmological studies. Numerous works have been conducted to detect such features in global scalesas well as in confined areas of the sky. However, results are often contradictory and their interpretation remains controversial. Aims.We i… ▽ More Context.Coherence in the characteristics of neighboring sources in 2D and 3D space may suggest the existence of large-scale cosmic structures, which are useful for cosmological studies. Numerous works have been conducted to detect such features in global scalesas well as in confined areas of the sky. However, results are often contradictory and their interpretation remains controversial. Aims.We investigate the potential alignment of parsec-scale radio jets in localized regions of the coordinates-redshift space. Methods.We use data from the Astrogeo VLBI FITS image database to deduce jet directions of radio sources. We perform the search for statistical alignments between nearby sources and explore the impact of instrumental biases. Results.We unveil four regions for which the alignment between jet directions deviates from randomness at a significance level of more than 5 sigma and is unlikely due to instrumental systematics. Intriguingly, their locations coincide with other known large-scale cosmic structures and/or regions of alignments. Conclusions.If the alignments found are the result of physical processes, the discovered regions may designate some of the largest structures known to date. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: 11 pages, 9 figures

Journal ref: A&A 653, A123 (2021)

Extragalactic Magnetism with SOFIA (Legacy Program) -- I: The magnetic field in the multi-phase interstellar medium of M51

Authors: Alejandro S. Borlaff, Enrique Lopez-Rodriguez, Rainer Beck, Rodion Stepanov, Eva Ntormousi, Annie Hughes, Konstantinos Tassis, Pamela M. Marcum, Lucas Grosset, John E. Beckman, Leslie Proudfit, Susan E. Clark, Tanio Díaz-Santos, Sui Ann Mao, William T. Reach, Julia Roman-Duval, Kandaswamy Subramanian, Le Ngoc Tram, Ellen G. Zweibel, SOFIA Legacy Team

Abstract: The recent availability of high-resolution far-infrared (FIR) polarization observations of galaxies using HAWC+/SOFIA has facilitated studies of extragalactic magnetic fields in the cold and dense molecular disks.We investigate if any significant structural differences are detectable in the kpc-scale magnetic field of the grand design face-on spiral galaxy M51 when traced within the diffuse (radio… ▽ More The recent availability of high-resolution far-infrared (FIR) polarization observations of galaxies using HAWC+/SOFIA has facilitated studies of extragalactic magnetic fields in the cold and dense molecular disks.We investigate if any significant structural differences are detectable in the kpc-scale magnetic field of the grand design face-on spiral galaxy M51 when traced within the diffuse (radio) and the dense and cold (FIR) interstellar medium (ISM). Our analysis reveals a complex scenario where radio and FIR polarization observations do not necessarily trace the same magnetic field structure. We find that the magnetic field in the arms is wrapped tighter at 154um than at 3 and 6 cm; statistically significant lower values for the magnetic pitch angle are measured at FIR in the outskirts (R > 7 kpc) of the galaxy. This difference is not detected in the interarm region. We find strong correlations of the polarization fraction and total intensity at FIR and radio with the gas column density and 12CO(1-0) velocity dispersion. We conclude that the arms show a relative increase of small-scale turbulent B-fields at regions with increasing column density and dispersion velocities of the molecular gas. No correlations are found with HI neutral gas. The star formation rate shows a clear correlation with the radio polarized intensity, which is not found in FIR, pointing to a small-scale dynamo-driven B-field amplification scenario. This work shows that multi-wavelength polarization observations are key to disentangling the interlocked relation between star formation, magnetic fields, and gas kinematics in the multi-phase ISM. △ Less

Submitted 8 September, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: 35 pages, 25 figures, 5 tables. Accepted for publication in ApJ

The Time-Dependent Distribution of Optical Polarization Angle Changes in Blazars

Authors: S. Kiehlmann, D. Blinov, I. Liodakis, V. Pavlidou, A. C. S. Readhead, E. Angelakis, C. Casadio, T. Hovatta, N. Kylafis, A. Mahabal, N. Mandarakas, I. Myserlis, G. V. Panopoulou, T. J. Pearson, A. Ramaprakash, P. Reig, R. Skalidis, A. Slowikowska, K. Tassis, J. A. Zensus

Abstract: At optical wavelengths, blazar Electric Vector Position Angle (EVPA) rotations linked with gamma-ray activity have been the subject of intense interest and systematic investigation for over a decade. One difficulty in the interpretation of EVPA rotations is the inherent 180° ambiguity in the measurements. It is therefore essential, when studying EVPA rotations, to ensure that the typical time-inte… ▽ More At optical wavelengths, blazar Electric Vector Position Angle (EVPA) rotations linked with gamma-ray activity have been the subject of intense interest and systematic investigation for over a decade. One difficulty in the interpretation of EVPA rotations is the inherent 180° ambiguity in the measurements. It is therefore essential, when studying EVPA rotations, to ensure that the typical time-interval between successive observations -- i.e. the cadence -- is short enough to ensure that the correct modulo 180° value is selected. This optimal cadence depends on the maximum intrinsic EVPA rotation speed in blazars, which is currently not known. In this paper we address the following questions for the RoboPol sample: What range of rotation speeds for rotations greater than 90° can we expect? What observation cadence is required to detect such rotations? Have rapid rotations been missed in EVPA rotation studies thus far? What fraction of data is affected by the ambiguity? And how likely are detected rotations affected by the ambiguity? We answer these questions with three seasons of optical polarimetric observations of a statistical sample of blazars sampled weekly with the RoboPol instrument and an additional season with daily observations. We model the distribution of EVPA changes on time scales from 1-30 days and estimate the fraction of changes exceeding 90°. We show that at least daily observations are necessary to measure >96% of optical EVPA variability in the RoboPol sample of blazars correctly and that intra-day observations are needed to measure the fastest rotations that have been seen thus far. △ Less

Submitted 14 July, 2021; v1 submitted 6 April, 2021; originally announced April 2021.

Comments: accepted to MNRAS, 20 pages, 15 figures

WALOP-South: A wide-field one-shot linear optical polarimeter for PASIPHAE survey

Authors: Siddharth Maharana, John A. Kypriotakis, A. N. Ramaprakash, Pravin Khodade, Chaitanya Rajarshi, Bhushan S. Joshi, Pravin Chordia, Ramya M. Anche, Shrish Mishra, Dmitry Blinov, Hans Kristian Eriksen, Tuhin Ghosh, Eirik Gjerløw, Nikolaos Mandarakas, Georgia V. Panopoulou, Vasiliki Pavlidou, Timothy J. Pearson, Vincent Pelgrims, Stephen B. Potter, Anthony C. S. Readhead, Raphail Skalidis, Konstantinos Tassis, Ingunn K. Wehus

Abstract: WALOP (Wide-Area Linear Optical Polarimeter)-South, to be mounted on the 1m SAAO telescope in South Africa, is first of the two WALOP instruments currently under development for carrying out the PASIPHAE survey. Scheduled for commissioning in the year 2021, the WALOP instruments will be used to measure the linear polarization of around $10^{6}$ stars in the SDSS-r broadband with $0.1~\%$ polarimet… ▽ More WALOP (Wide-Area Linear Optical Polarimeter)-South, to be mounted on the 1m SAAO telescope in South Africa, is first of the two WALOP instruments currently under development for carrying out the PASIPHAE survey. Scheduled for commissioning in the year 2021, the WALOP instruments will be used to measure the linear polarization of around $10^{6}$ stars in the SDSS-r broadband with $0.1~\%$ polarimetric accuracy, covering 4000 square degrees in the Galactic polar regions. The combined capabilities of one-shot linear polarimetry, high polarimetric accuracy ($< 0.1~\%$) and polarimetric sensitivity ($< 0.05~\%$), and a large field of view (FOV) of $35\times35~arcminutes$ make WALOP-South a unique astronomical instrument. In a single exposure, it is designed to measure the Stokes parameters $I$, $q$ and $u$ in the SDSS-r broadband and narrowband filters between $500-700~nm$. During each measurement, four images of the full field corresponding to the polarization angles of $0^{\circ}$, $45^{\circ}$, $90^{\circ}$ and $135^{\circ}$ will be imaged on four detectors and carrying out differential photometry on these images will yield the Stokes parameters. Major challenges in designing WALOP-South instrument include- (a) in the optical design, correcting for the spectral dispersion introduced by large split angle Wollaston Prisms used as polarization analyzers as well as aberrations from the wide field, and (b) making an optomechanical design adherent to the tolerances required to obtain good imaging and polarimetric performance under all temperature conditions as well as telescope pointing positions. We present the optical and optomechanical design for WALOP-South which overcomes these challenges. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: Presented at SPIE Astronomical Telescopes + Instrumentation 2020. 12 Pages, 6 Figures, 6 Tables

Journal ref: Proceedings Volume 11447, Ground-based and Airborne Instrumentation for Astronomy VIII; 114475E (2020)

WALOP-South: A Four Camera One Shot Imaging Polarimeter for PASIPHAE Survey. Paper I -- Optical Design

Authors: Siddharth Maharana, John A. Kypriotakis, A. N. Ramaprakash, Chaitanya Rajarshi, Ramya M. Anche, Shrish, Dmitry Blinov, Hans Kristian Eriksen, Tuhin Ghosh, Eirik Gjerløw, Nikolaos Mandarakas, Georgia V. Panopoulou, Vasiliki Pavlidou, Timothy J. Pearson, Vincent Pelgrims, Stephen B. Potter, Anthony C. S. Readhead, Raphael Skalidis, Konstantinos Tassis, Ingunn K. Wehus

Abstract: The WALOP-South instrument will be mounted on the 1 m SAAO telescope in South Africa as part of the PASIPHAE program to carry out a linear imaging polarization survey of the Galactic polar regions in the optical band. Designed to achieve polarimetric sensitivity of $0.05~\%$ across a $35\times35$ arcminute field of view, it will be capable of measuring the Stokes parameters I, q and u in a single… ▽ More The WALOP-South instrument will be mounted on the 1 m SAAO telescope in South Africa as part of the PASIPHAE program to carry out a linear imaging polarization survey of the Galactic polar regions in the optical band. Designed to achieve polarimetric sensitivity of $0.05~\%$ across a $35\times35$ arcminute field of view, it will be capable of measuring the Stokes parameters I, q and u in a single exposure in the SDSS-r broadband and narrowband filters between $0.5~μm - 0.7~μm$. For each measurement, four images of the full field corresponding to linear polarization angles of 0 deg, 45 deg, 90 deg and 135 deg in the instrument coordinate system will be created on four detectors from which the Stokes parameters can be found using differential photometry. In designing the optical system, major challenges included correcting for the dispersion introduced by large split angle Wollaston Prisms used as analysers as well as other aberrations from the entire field to obtain imaging quality PSF at the detector. We present the optical design of the WALOP-South instrument which overcomes these challenges and delivers near seeing limited PSFs for the entire field of view. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 31 pages, 18 Figures and 8 Tables. Accepted in the Journal of Astronomical Telescopes, Instruments, and Systems

Evidence for Line-of-Sight Frequency Decorrelation of Polarized Dust Emission in $Planck$ Data

Authors: V. Pelgrims, S. E. Clark, B. S. Hensley, G. V. Panopoulou, V. Pavlidou, K. Tassis, H. K. Eriksen, I. K. Wehus

Abstract: If a single line of sight (LOS) intercepts multiple dust clouds of different spectral energy distributions and magnetic field orientations, the frequency scaling of each of the Stokes $Q$ and $U$ parameters of thermal dust emission may be different ("LOS frequency decorrelation"). We present first evidence for LOS frequency decorrelation in $Planck$ data. We use independent, neutral-hydrogen--meas… ▽ More If a single line of sight (LOS) intercepts multiple dust clouds of different spectral energy distributions and magnetic field orientations, the frequency scaling of each of the Stokes $Q$ and $U$ parameters of thermal dust emission may be different ("LOS frequency decorrelation"). We present first evidence for LOS frequency decorrelation in $Planck$ data. We use independent, neutral-hydrogen--measurements of the number of clouds per LOS and the magnetic field orientation in each cloud to select two sets of sightlines: (i) a target sample (pixels likely to exhibit LOS frequency decorrelation); (ii) a control sample (pixels lacking complex LOS structure). We test the null hypothesis that LOS frequency decorrelation is not detectable in $Planck$ 353 and 217~GHz polarization data at high Galactic latitudes. The data reject this hypothesis at high significance. The detection is robust against choice of CMB map and map-making pipeline. The observed change in polarization angle due to LOS frequency decorrelation is detectable above the $Planck$ noise level. The probability that the detected effect is due to noise alone ranges from $5\times 10^{-2}$ to $4\times 10^{-7}$, depending on the CMB subtraction algorithm and treatment of residual systematics; correcting for residual systematics increases the significance of the effect. The LOS decorrelation effect is stronger for sightlines with more misaligned magnetic fields, as expected. We estimate that an intrinsic variation of $\sim15\%$ in the ratio of 353 to 217~GHz polarized emission between clouds is sufficient to reproduce the measured effect. Our finding underlines the importance of ongoing studies to map the 3D structure of the magnetized dusty ISM that could help component separation methods to account for frequency decorrelation effects in CMB polarization studies. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Comments: 17 pages, 15 figures. Accepted for publication by A&A

Journal ref: A&A 647, A16 (2021)

RoboPol: AGN polarimetric monitoring data

Authors: D. Blinov, S. Kiehlmann, V. Pavlidou, G. V. Panopoulou, R. Skalidis, E. Angelakis, C. Casadio, E. N. Einoder, T. Hovatta, K. Kokolakis, A. Kougentakis, A. Kus, N. Kylafis, E. Kyritsis, A. Lalakos, I. Liodakis, S. Maharana, E. Makrydopoulou, N. Mandarakas, G. M. Maragkakis, I. Myserlis, I. Papadakis, G. Paterakis, T. J. Pearson, A. N. Ramaprakash , et al. (7 additional authors not shown)

Abstract: We present uniformly reprocessed and re-calibrated data from the RoboPol programme of optopolarimetric monitoring of active galactic nuclei (AGN), covering observations between 2013, when the instrument was commissioned, and 2017. In total, the dataset presented in this paper includes 5068 observations of 222 AGN with Dec > -25 deg. We describe the current version of the RoboPol pipeline that was… ▽ More We present uniformly reprocessed and re-calibrated data from the RoboPol programme of optopolarimetric monitoring of active galactic nuclei (AGN), covering observations between 2013, when the instrument was commissioned, and 2017. In total, the dataset presented in this paper includes 5068 observations of 222 AGN with Dec > -25 deg. We describe the current version of the RoboPol pipeline that was used to process and calibrate the entire dataset, and we make the data publicly available for use by the astronomical community. Average quantities summarising optopolarimetric behaviour (average degree of polarization, polarization variability index) are also provided for each source we have observed and for the time interval we have followed it. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Comments: Accepted to MNRAS

High-accuracy estimation of magnetic field strength in the interstellar medium from dust polarization

Authors: Raphael Skalidis, Konstantinos Tassis

Abstract: Dust polarization is a powerful tool for studying the magnetic field properties in the interstellar medium (ISM). However, it does not provide a direct measurement of its strength. Different methods havebeen developed which employ both polarization and spectroscopic data in order to infer the field strength. The most widely applied methods have been developed by Davis (1951), Chandrasekhar & Fermi… ▽ More Dust polarization is a powerful tool for studying the magnetic field properties in the interstellar medium (ISM). However, it does not provide a direct measurement of its strength. Different methods havebeen developed which employ both polarization and spectroscopic data in order to infer the field strength. The most widely applied methods have been developed by Davis (1951), Chandrasekhar & Fermi (1953) (DCF), Hildebrand et al. (2009) and Houde et al.(2009) (HH09). They rely on the assumption that isotropic turbulent motions initiate the propagation of Alvfén waves. Observations,however, indicate that turbulence in the ISM is anisotropic and non-Alfvénic (compressible) modes may be important. Our goal is to develop a new method for estimating the field strength in the ISM, which includes the compressible modes and does not contradict the anisotropic properties of turbulence. We use simple energetics arguments that take into account the compressible modes to estimate the strength of the magnetic field. We derive the following equation: $B_{0}=\sqrt{2 πρ} δv /\sqrt{δθ}$, where $ρ$ is the gas density, $δv$ is the rms velocity as derived from the spread of emission lines, and $δθ$ is the dispersion of polarization angles. We produce synthetic observations from 3D MHD simulationsand we assess the accuracy of our method by comparing the true field strength with the estimates derived from our equation. We find a mean relative deviation of $17 \%$. The accuracy of our method does not depend on the turbulence properties of the simulated model. In contrast DCF and HH09 systematically overestimate the field strength. HH09 produces accurate results only for simulations with high sonic Mach numbers. △ Less

Submitted 4 February, 2021; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: Accepted for publication to Astronomy & Astrophysics

Journal ref: A&A 647, A186 (2021)

HAWC+ Far-Infrared Observations of the Magnetic Field Geometry in M51 and NGC 891

Authors: Terry Jay Jones, Jin-Ah Kim, C. Darren Dowell, Mark R. Morris, Jorge L. Pineda, Dominic J. Benford, Marc Berthoud, David T. Chuss, Daniel A. Dale, L. M. Fissel, Paul F. Goldsmith, Ryan T. Hamilton, Shaul Hanany, Doyal A. Harper, Thomas K. Henning, Alex Lazarian, Leslie W. Looney, Joseph M. Michail, Giles Novak, Fabio P. Santos, Kartik Sheth, Javad Siah, Gordon J. Stacey, Johannes Staguhn, Ian W. Stephens , et al. (7 additional authors not shown)

Abstract: SOFIA HAWC+ polarimetry at $154~\micron$ is reported for the face-on galaxy M51 and the edge-on galaxy NGC 891. For M51, the polarization vectors generally follow the spiral pattern defined by the molecular gas distribution, the far-infrared (FIR) intensity contours, and other tracers of star formation. The fractional polarization is much lower in the FIR-bright central regions than in the outer r… ▽ More SOFIA HAWC+ polarimetry at $154~\micron$ is reported for the face-on galaxy M51 and the edge-on galaxy NGC 891. For M51, the polarization vectors generally follow the spiral pattern defined by the molecular gas distribution, the far-infrared (FIR) intensity contours, and other tracers of star formation. The fractional polarization is much lower in the FIR-bright central regions than in the outer regions, and we rule out loss of grain alignment and variations in magnetic field strength as causes. When compared with existing synchrotron observations, which sample different regions with different weighting, we find the net position angles are strongly correlated, the fractional polarizations are moderately correlated, but the polarized intensities are uncorrelated. We argue that the low fractional polarization in the central regions must be due to significant numbers of highly turbulent segments across the beam and along lines of sight in the beam in the central 3 kpc of M51. For NGC 891, the FIR polarization vectors within an intensity contour of 1500 $\rm{MJy~sr^{-1}}$ are oriented very close to the plane of the galaxy. The FIR polarimetry is probably sampling the magnetic field geometry in NGC 891 much deeper into the disk than is possible with NIR polarimetry and radio synchrotron measurements. In some locations in NGC 891 the FIR polarization is very low, suggesting we are preferentially viewing the magnetic field mostly along the line of sight, down the length of embedded spiral arms. There is tentative evidence for a vertical field in the polarized emission off the plane of the disk. △ Less

Submitted 18 August, 2020; originally announced August 2020.

Comments: AJ, Accepted, In Press

A dynamo amplifies the magnetic field of a Milky-Way-like galaxy

Authors: Evangelia Ntormousi, Konstantinos Tassis, Fabio Del Sordo, Francesca Fragkoudi, Rüdiger Pakmor

Abstract: The magnetic fields of spiral galaxies are so strong that they cannot be primordial. Their typical values are over one billion times higher than any value predicted for the early Universe. Explaining this immense growth and incorporating it in galaxy evolution theories is one of the long-standing challenges in astrophysics. So far, the most successful theory for the sustained growth of the galacti… ▽ More The magnetic fields of spiral galaxies are so strong that they cannot be primordial. Their typical values are over one billion times higher than any value predicted for the early Universe. Explaining this immense growth and incorporating it in galaxy evolution theories is one of the long-standing challenges in astrophysics. So far, the most successful theory for the sustained growth of the galactic magnetic field is the alpha-omega dynamo. This theory predicts a characteristic dipolar or quadrupolar morphology for the galactic magnetic field, which has been observed in external galaxies. However, so far, there has been no direct demonstration of a mean-field dynamo operating in direct, multi-physics simulations of spiral galaxies. We do so in this work. We employ numerical models of isolated, star-forming spiral galaxies that include a magnetized gaseous disk, a dark matter halo, stars, and stellar feedback. Naturally, the resulting magnetic field has a complex morphology that includes a strong random component. Using a smoothing of the magnetic field on small scales, we are able to separate the mean from the turbulent component and analyze them individually. We find that a mean-field dynamo naturally occurs as a result of the dynamical evolution of the galaxy and amplifies the magnetic field by an order of magnitude over half a Gyr. Despite the highly dynamical nature of these models, the morphology of the mean component of the field is identical to analytical predictions. This result underlines the importance of the mean-field dynamo in galactic evolution. Moreover, by demonstrating the natural growth of the magnetic field in a complex galactic environment, it brings us a step closer to understanding the cosmic origin of magnetic fields. △ Less

Submitted 22 June, 2020; originally announced June 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 641, A165 (2020)

Turnaround radius of galaxy clusters in N-body simulations

Authors: Giorgos Korkidis, Vasiliki Pavlidou, Konstantinos Tassis, Evangelia Ntormousi, Theodore N. Tomaras, Konstantinos Kovlakas

Abstract: We use N-body simulations to examine whether a characteristic turnaround radius, as predicted from the spherical collapse model in a $\rm {ΛCDM}$ Universe, can be meaningfully identified for galaxy clusters, in the presence of full three-dimensional effects. We use The Dark Sky Simulations and Illustris-TNG dark-matter--only cosmological runs to calculate radial velocity profiles around collapsed… ▽ More We use N-body simulations to examine whether a characteristic turnaround radius, as predicted from the spherical collapse model in a $\rm {ΛCDM}$ Universe, can be meaningfully identified for galaxy clusters, in the presence of full three-dimensional effects. We use The Dark Sky Simulations and Illustris-TNG dark-matter--only cosmological runs to calculate radial velocity profiles around collapsed structures, extending out to many times the virial radius $R_{200}$. There, the turnaround radius can be unambiguously identified as the largest non-expanding scale around a center of gravity. We find that: (a) Indeed, a single turnaround scale can meaningfully describe strongly non-spherical structures. (b) For halos of masses $M_{200}>10^{13}M_\odot$, the turnaround radius $R_{ta}$ scales with the enclosed mass $M_{ta}$ as $M_{ta}^{1/3}$, as predicted by the spherical collapse model. (c) The deviation of $R_{ta}$ in simulated halos from the spherical collapse model prediction is insensitive to halo asphericity. Rather, it is sensitive to the tidal forces due to massive neighbors when such are present. (d) Halos exhibit a characteristic average density within the turnaround scale. This characteristic density is dependent on cosmology and redshift. For the present cosmic epoch and for concordance cosmological parameters ($Ω_m \sim 0.7$; $Ω_Λ\sim 0.3$) turnaround structures exhibit an average matter density contrast with the background Universe of $δ\sim 11$. Thus $R_{ta}$ is equivalent to $R_{11}$ -- in a way analogous to defining the "virial" radius as $R_{200}$ -- with the advantage that $R_{11}$ is shown in this work to correspond to a kinematically relevant scale in N-body simulations. △ Less

Submitted 17 December, 2019; originally announced December 2019.

Comments: Submitted to A&A, nine pages, 10 figures

Journal ref: A&A 639, A122 (2020)

Eliminating artefacts in Polarimetric Images using Deep Learning

Authors: Dhruv Paranjpye, Ashish Mahabal, A. N. Ramaprakash, Gina Panopoulou, Kieran Cleary, Anthony Readhead, Dmitry Blinov, Kostas Tassis

Abstract: Polarization measurements done using Imaging Polarimeters such as the Robotic Polarimeter are very sensitive to the presence of artefacts in images. Artefacts can range from internal reflections in a telescope to satellite trails that could contaminate an area of interest in the image. With the advent of wide-field polarimetry surveys, it is imperative to develop methods that automatically flag ar… ▽ More Polarization measurements done using Imaging Polarimeters such as the Robotic Polarimeter are very sensitive to the presence of artefacts in images. Artefacts can range from internal reflections in a telescope to satellite trails that could contaminate an area of interest in the image. With the advent of wide-field polarimetry surveys, it is imperative to develop methods that automatically flag artefacts in images. In this paper, we implement a Convolutional Neural Network to identify the most dominant artefacts in the images. We find that our model can successfully classify sources with 98\% true positive and 97\% true negative rates. Such models, combined with transfer learning, will give us a running start in artefact elimination for near-future surveys like WALOP. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: 7 pages, 15 figures

AEDGE: Atomic Experiment for Dark Matter and Gravity Exploration in Space

Authors: Yousef Abou El-Neaj, Cristiano Alpigiani, Sana Amairi-Pyka, Henrique Araujo, Antun Balaz, Angelo Bassi, Lars Bathe-Peters, Baptiste Battelier, Aleksandar Belic, Elliot Bentine, Jose Bernabeu, Andrea Bertoldi, Robert Bingham, Diego Blas, Vasiliki Bolpasi, Kai Bongs, Sougato Bose, Philippe Bouyer, Themis Bowcock, William Bowden, Oliver Buchmueller, Clare Burrage, Xavier Calmet, Benjamin Canuel, Laurentiu-Ioan Caramete , et al. (107 additional authors not shown)

Abstract: We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also compl… ▽ More We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material. △ Less

Submitted 10 October, 2019; v1 submitted 2 August, 2019; originally announced August 2019.

Comments: V2 -- added support authors

Report number: KCL-PH-TH/2019-65, CERN-TH-2019-126

Journal ref: EPJ Quantum Technol. 7, 6 (2020)

SOFIA/HAWC+ traces the magnetic fields in NGC 1068

Authors: E. Lopez-Rodriguez, C. D. Dowell, T. J. Jones, D. A. Harper, M. Berthoud, D. Chuss, D. A. Dale, J. A. Guerra, R. T. Hamilton, L. W. Looney, J. M. Michail, R. Nikutta, G. Novak, F. P. Santos, K. Sheth, J. Siah, J. Staguhn, I. W. Stephens, K. Tassis, C. Q. Trinh, D. Ward-Thompson, M. Werner, E. J. Wollack, E. Zweibel

Abstract: We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $μ$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and… ▽ More We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $μ$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $μ$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um. △ Less

Submitted 15 November, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: 13 pages, 7 figures, Accepted for publication by ApJ