Periodic Oscillations in LASCO Coronal Mass Ejection Speeds: Space Seismology
Coronal mass ejections (CMEs) are energetic expulsions of organized magnetized plasma from the Sun. They are intensively studied because of their significant impact on Earth. The first CME was discovered in 1971 using the seventh Orbiting Solar Observatory (OSO-7). The longest series of observations have been carried out by the Large Angle and Spectrometric Coronagraph (LASCO) on board the Solar and Heliospheric Observatory (SOHO) mission.
Now we know that the magnetic field in CMEs and in the sheath ahead of shock-driving CMEs, especially the strength of the magnetic field component oriented in the direction opposite to that of Earth’s horizontal magnetic field (Bz), is crucial in causing geomagnetic storms along with the speed (V) with which the CME hits Earth’s magnetosphere. Unfortunately, routine measurements of magnetic fields, due to high temperature, in the outer solar atmosphere, where CMEs originate, are not available. To determine the magnetic field in the solar corona, we must apply indirect methods. It seems that one such method may be based on oscillations in the speed of CME. CMEs have an organized flux-rope magnetic structure. Therefore, the reconnection of the magnetic field during ejection may excite periodic speed oscillations. Periodic speed and intensity oscillations are well established in the case of quiescent prominences, which in many cases are progenitors of CMEs. These oscillations have a range of periodicity from a few minutes up to the order of an hour. We considered all the CMEs with at least 11 height–time measurements included in the CDAW database in the period of 1996–2019. It allowed us to evaluate speed–time profiles with enough accuracy to analyze oscillations in the speed of CMEs. Of the considered CMEs, 22% revealed periodic speed fluctuations during their expansion in the interplanetary medium. This means that speed oscillation is a frequent phenomenon associated with CME propagation.
We have obtained the following new important results. We demonstrated that the average values of the basic attributes of oscillations (amplitude and period) are significantly correlated with cycles of solar activity. In the maxima of solar activity, on average, pulsation periods are ≈240 minutes and in minima as much as ≈300 minutes. This means that periods of oscillations depend on the physical properties of ejections. It is assumed that the flux-rope structure inside CMEs can be approximated as a stretched elastic string of non-uniform density. Using this approximation, we estimated that, on average, the CME internal magnetic field varies from 18 up to 25 mG between the minimum and maximum of solar activity.
The obtained results show that a detailed analysis of speed oscillations can be a very efficient tool for studying not only the physical properties of the ejections themselves but also the condition of the interplanetary medium in which they expand.
Grzegorz Michalek et al., Periodic Oscillations in LASCO Coronal Mass Ejection Speeds: Space Seismology, The Astrophysical Journal Letters (marzec 2022)
The research was conducted at the Department of High Energy Astrophysics of the Jagiellonian University’s Astronomical Observatory (OA UJ).
13 years of monitoring the optical variability of quasars at OAUJ
Active galaxies are some of the most energetic objects in the Universe. Processes of accretion of matter onto supermassive black holes occur in their active nuclei (AGNs), constituting a source of radiation over a wide range of the electromagnetic spectrum. They are also objects exhibiting radiation variability. Studying the nature of this observed variability we can analyse the physical processes responsible for changes in the brightness of the active nucleus.
The physical processes behind the variability in the optical range (on time scales from days to decades) are the subject of numerous studies. Theoretical considerations suggest that several mechanisms can explain both the time scales and the amplitudes of the brightness variations observed in AGNs. These may include instabilities in the accretion disk, processes associated with the evolution of massive stars that explode as supernovae, or variabilities caused by the gravitational microlensing. Each of these mechanisms results in a different character of variability, possible to study using statistical methods.
In 2009, systematic monitoring of the optical variability of the eight quasars has started in our Observatory. The selected objects have extended radio structures, with five of which are classified as so-called giant radio quasars, with sizes exceeding 0.7 Mpc. Initially, observations were performed using a fifty centimetre Cassegrain telescope located at the Astronomical Observatory of the Jagiellonian University. However, after some time, the telescope located at Suhora Observatory (UP) has also been used for the observations, as well as the telescopes belonging to the network of robotic telescopes – SKYNET. Thanks to the many people involved in the observations, it has been possible to obtain the light curves of the quasars, covering 13-year period of their variability. These are unique data not only because of the length of the observations period in question, but also due to the frequent sampling – the observations of each quasar were made on average 3 times a month.
The obtained light curves were used to investigate the characteristic time scales of the variation and the physical process that would explain the observed nature of the variation. For this purpose, two statistical methods was used: the so-called structure function (FS) analysis and power spectral density analysis (PSD). Based on the obtained FS and PSD slopes, it can be estimated that the most likely process that could be responsible for the variability of the observed quasars is the instabilities of the accretion disk, and the nature of this variability can be described as red noise or a damped random-walk process. Our analyses have also shown that the characteristic time scales of the variability of this type of objects are much longer than their monitoring time so far. They can span up to several hundred or even millions years, which is not achievable using conventional single-object observation methods.
An interesting though not fully understood obtained result is a strong anti-correlation between the size of radio structures and the slope of the PSD function (correlation coefficient of 0.86). The anti-correlation may indicate some link between the size of the radio source and the nature of its variability, but it is now based on data for eight objects only, making further study of a larger sample of this kind of objects necessary to confirm it.
Illustration 2: Structure function for [HB89] 1721+343. The plot shows the value of the coefficient of the FS slope (α) and the inflection point (τ) corresponding to the characteristic time scale of the variability. Right: PSD function for [HB89] 1721+343 with the value of PSD slope (β) given. Credit: The Authors.
Illustration 3: The anti-correlation between the slope of the PSD function (β) and the size of the radio structure of quasars. Credit: The Authors.
A Kuźmicz et al., Optical Variability of Eight FRII-type Quasars with 13 yr Photometric Light Curves, ApJS 263 16 (2022).
The research was conducted at the Department of Stellar and Extragalactic Astronomy and Department of Radioastronomy and Space Physics of the Jagiellonian University’s Astronomical Observatory (OA UJ). This research was supported by the Polish National Science Centre grants UMO-2018/29/B/ST9/01793 and UMO-2018/29/B/ST9/02298. The quasar light-curve simulations were performed using the PLGrid Infrastructure.
Crisis in Cosmology: statistics on the rescue
International team of scientists from Italy, Poland, and Japan has released a groundbreaking article with a solution to one of the biggest challenges of modern observational cosmology. A new statistical method was applied in order to remove the impact of selection biases on cosmological computations while avoiding circular reasoning. The article has been accepted for publication in The Astrophysical Journal Supplement Series.
Figure 1. In the center of the image, there is a quasar lensed by the galaxy. Four evenly distributed dots are images of the same quasar, as an effect of gravitational lensing of beams of light by the galaxy in front of the quasar. Credit: ESA/Hubble, NASA, Suyu et al.
Modern cosmology has many open questions: how to solve the discrepancy between the measurements of the Hubble constant? What is the curvature of space-time? What is the nature of dark energy? And many, many more. Astronomers are trying to shed new light on these issues by observing the Universe farther and farther away from the Earth. This is driven by the fact that predictions from different cosmological models have minor differences at a low distance. Still, those differences begin to be significant for objects very far away from us.
In order to probe the high-distance Universe, scientists have to use the most luminous objects that ever existed, quasars and Gamma-Ray Bursts. Quasars are powered by the accretion of gas onto the supermassive black holes, residing in the centers of galaxies. Gamma-Ray Bursts originate when a massive star explodes, or when two neutron stars, or a neutron star and a black hole merge. Since the discovery of these objects, many correlations between their physical quantities have been studied.
These correlations allow us to determine the distance of these objects from us. This is crucial for cosmological analyses, since we have the theoretical formula for the distance in a given cosmological model, which includes the speed of the expansion of the given point in the Universe. This given speed can be measured spectroscopically, thus we can compute the theoretical distance of a particular object for a given cosmological model. The comparison of the theoretical distance and the observed one allows us to test how well the given cosmological model fits the data.
But looking at high distances has shortcomings: the farther away we look, the less luminous events we see. This phenomenon is called the Malmquist bias effect. This can induce an artificial correlation between the observed parameters, so how do we know if the given correlation is a physical feature of the objects, or just induced by a selection effect? This can be tested via a method proposed by Efron & Petrosian (1992) and already successfully developed for the Risaliti-Lusso relation of the quasars by Dainotti et al. 2022 (article by the same authors as the mainly discussed work). This technique eliminates the correlation between the studied parameters and the observational spectroscopic parameter determining the distance (the redshift).
This method creates yet another problem. Many correlations involve luminosity, and in order to calculate the luminosity, one needs to fix the parameters of the cosmological model. Then, with those fixed values, one determines the correction for the selection bias effect. Consequently, this correction has been applied within the assumed cosmological model. Thus, if we would like to use these corrected values to test some cosmological model, our results will be again biased by the model assumed a priori. This represents the so-called circular reasoning problem and it cannot be applied to scientific approaches.
The new article completely overcomes this issue, by creating a new, more general treatment than the one of Efron & Petrosian (1992). This treatment was first discussed by Dainotti et al. (2022), and later has been fully implemented by Lenart & Bargiacchi et. al (2023). Lenart & Bargiacchi have proposed a method which treats the correction for selection bias as a function of the cosmological parameters, allowing to obtain a general formula of correction, valid for any value of cosmological parameters, which can be applied without assuming any values a priori.
The application of this method to the so-called Risaliti-Lusso correlation linking two luminosities for quasars allowed for the first time in the literature to obtain a 2 sigma value of ΩM, the parameter denoting the current amount of matter in the Universe. Moreover, the authors applied the following analysis to study the impact of high-distance measurements on the so-called H0 tension (the incompatibility of the two measurements of the Hubble constant: the one involving supernovae type Ia and the one obtained from the Microwave Cosmic Radiation). Surprisingly, quasars seem to prefer the value of the Hubble constant in between the values obtained for supernovae type Ia and Microwave Cosmic Radiation. This suggests that this discrepancy might be driven by some still unknown physics, and not by the statistical effects. Although in order to confirm such a conclusion, the new correction has yet to be applied also to other samples like supernovae type Ia. The overall results enhance a further significant improvement in cosmological measurements and an extension of our knowledge about the Universe at high distances.
Figure 2. On the left, uncorrected correlation between luminosity measured in an X-ray band and visible light, while on the right, the same correlation, but corrected for evolution and selection bias with the described method. On the left, we can see a clear evolution of the correlation with redshift, while on the right no such evolution can be spotted. Credit: Original publication.
Figure 3. The values of the H0 obtained for different samples and methods. Grey vertical lines represent the values of H0 computed with the Cosmic Microwave Background radiation data (on the left) and supernovae Ia data (on the right). After the addition of quasars data to the supernovae Ia, one can spot that results fall perfectly between the two before mentioned values. Credit: Original publication.
Aleksander Łukasz Lenart, Giada Bargiacchi, Maria Giovanna Dainotti, Shigehiro Nagataki, Salvatore Capozziello, A bias-free cosmological analysis with quasars alleviating H0 tension
The publication under discussion was made possible, amongst others, thanks to financial support from the Astronomical Observatory of Jagiellonian University in Cracow and the Program Council for Mathematical and Natural Sciences Studies at Jagiellonian University. We also acknowledge the National Astronomical Observatory of Japan and RIKEN for their support in realising the work and partially funding the publication.