Since the end of 2018, the Transiting Exoplanet Survey Satellite (TESS) provides high-quality space data on stellar photometry to the astronomical community. We present the results of an analysis of TESS photometric data for known slowly rotating magnetic chemically peculiar (mCP) stars. In general, mCP stars show an inhomogeneous distribution of elements in their stellar atmospheres that leads to spectroscopic (line profile) and photometric (light curve) variations over the rotation period. In the frame of the oblique magnetic rotator (OMR) model, patches of enhanced chemical abundance on the stellar surface reveal the frequency of stellar rotation. Using this approach, we have compiled a list of slowly rotating mCP stars with rotation periods longer than two days from the analysis of the photometric data provided by TESS for the first eight sectors of observations. Slowly rotating mCP stars usually possess a hydrodynamically stable stellar atmosphere where a magnetic field can amplify the process of atomic diffusion and leads to the horizontal and vertical stratification of chemical abundances.

A portion of the upper main-sequence stars, called chemically peculiar (CP) stars, show significant abundance anomalies mainly due to atomic diffusion of chemical elements within the stellar atmospheres of these stars. Slowly rotating CP stars may have hydrodynamically stable atmospheres where the competition between the gravitational and radiative forces leads to atomic diffusion that can be responsible for the abundance peculiarities observed in CP stars. Recently, project VeSElkA (Vertical Stratification of Element Abundances) was initiated with the aim to detect and study the vertical stratification of element abundances in the atmospheres of CP stars. Some results from abundance analysis of several slowly rotating (V$\sin{i}$< 40 km/s) CP stars observed recently with ESPaDOnS are presented here. Signatures of vertical abundance stratification for several chemical elements have been found in the stellar atmospheres of HD~22920, HD~41076, HD~95608, HD~116235, HD~148330 and HD~157087.

HD~63401 is known magnetic chemically peculiar (mCP) star that shows a slow rotation and probably possesses a hydrodynamically stable stellar atmosphere. In the latter case the atomic diffusion mechanism enforced by the magnetic field can lead to stratification of elemental abundance with optical depth. HD~63401 was recently observed with the space telescope \textit{TESS}, and its light curve shows distinct variability, which is usually detected in the $\alpha^2$ CVn type stars. Based on the analysis of the light curve we derived its rotational period $P$ = 2.414 $\pm$ 0.018~d and studied variability of the effective temperature with this period. The estimates of the effective temperature and surface gravity were obtained from the best fit of Balmer line profiles observed in seven high-resolution spectra acquired with the spectropolarimeter ESPaDOnS for HD~63401. The same spectra were used to perform an abundance analysis employing the modified ZEEMAN2 code. We have found that He, C, P, V, Y, and Dy are in deficit in stellar atmosphere of HD~63401, while Na, Al, Si, Fe, Zn, and Sr are significantly overabundant.

One of the key questions in the study of relativistic jets is how magnetic reconnection occurs and whether it can effectively accelerate electrons in the jet. We performed 3D particle-in-cell (PIC) simulations of a relativistic electron-proton jet of relatively large radius that carries a helical magnetic field. We focussed our investigation on the interaction between the jet and the ambient plasma and explore how the helical magnetic field affects the excitation of kinetic instabilities such as the Weibel instability (WI), the kinetic Kelvin-Helmholtz instability (kKHI), and the mushroom instability (MI). In our simulations these kinetic instabilities are indeed excited, and particles are accelerated. At the linear stage we observe recollimation shocks near the center of the jet. As the electron-proton jet evolves into the deep nonlinear stage, the helical magnetic field becomes untangled due to reconnection-like phenomena, and electrons are repeatedly accelerated as they encounter magnetic-reconnection events in the turbulent magnetic field.

Magnetic massive and intermediate-mass stars constitute a separate population whose properties are still not fully understood. Increasing the sample of known objects of this type would help answer fundamental questions regarding the origins and characteristics of their magnetic fields. The MOBSTER Collaboration seeks to identify candidate magnetic A, B and O stars and explore the incidence and origins of photometric rotational modulation using high-precision photometry from the Transiting Exoplanet Survey Satellite (\textit{TESS}) mission. In this contribution, we present an overview of our methods and planned targeted spectropolarimetric follow-up surveys.

In this contribution, we present the MOBSTER Collaboration, a large community effort to leverage high-precision photometry from the Transiting Exoplanet Survey Satellite (\textit{TESS}) in order to characterize the variability of magnetic massive and intermediate-mass stars. These data can be used to probe the varying column density of magnetospheric plasma along the line of sight for OB stars, thus improving our understanding of the interaction between surface magnetic fields and massive star winds. They can also be used to map out the brightness inhomogeneities present on the surfaces of Ap/Bp stars, informing present models of atomic diffusion in their atmospheres. Finally, we review our current and ongoing studies, which lead to new insights on this topic.

The new photometric data on pulsating Ap star HD~27463 obtained recently with the Transiting Exoplanet Survey Satellite (\textit{TESS}) are analysed to search for variability. Our analysis shows that HD~27463 exhibits two types of photometric variability. The low frequency variability with the period $P$ =~2.834274 $\pm$ 0.000008 d can be explained in terms of axial stellar rotation assuming the oblique magnetic rotator model and presence of surface abundance/brightness spots, while the detected high-frequency variations are characteristics of $\delta$~Scuti pulsations. From the analysis of Balmer line profiles visible in two FEROS spectra of HD~27463 we have derived its effective temperature and surface gravity, finding values that are close to those published for this star in the \textit{TESS} Input Catalogue (TIC). Knowing the rotation period and the v$\sin{i}$ value estimated from the fitting of Balmer line profiles we found that the rotational axis is inclined to the line of sight with an angle of $i=33\pm8\deg$. Our best-fitting model of the observed pulsation modes results in an overshoot parameter value $f_{ov} = 0.014$ and values of global stellar parameters that are in good agreement with the data reported in the TIC and with the data derived from fitting Balmer line profiles. This model indicates an age of 5.0 $\pm$~0.4 $\times 10^8$~yrs, which corresponds to a core hydrogen fraction of 0.33.

We investigated electromagnetic precursor wave emission in relativistic shocks by using two-dimensional particle-in-cell simulations. We found that the wave amplitude is significantly enhanced by a positive feedback process associated with ion-electron coupling through the wakefields for high magnetization. The wakefields collapse during the nonlinear process of the parametric decay instability in the near-upstream region, where nonthermal electrons and ions are generated. The intense coherent emission and the particle acceleration may opperate in high-energy astrophysical objects.

Recently, high-precision optical 2~min cadence light curves obtained with \emph{TESS} for targets located in the mission's defined first four sectors have been released. The majority of these high-cadence and high-precision measurements currently span $\sim28\,{\rm d}$, thereby allowing periodic variability occurring on timescales $\lesssim14\,{\rm d}$ to potentially be detected. Magnetic chemically peculiar (mCP) A-type stars are well known to exhibit rotationally modulated photometric variability that is produced by inhomogeneous chemical abundance distributions in their atmospheres. While mCP stars typically exhibit rotation periods that are significantly longer than those of non-mCP stars, both populations exhibit typical periods $\lesssim10\,{\rm d}$; therefore, the early \emph{TESS} releases are suitable for searching for rotational modulation of the light curves of both mCP and non-mCP stars. We present the results of our search for A-type stars that exhibit variability in their \emph{TESS} light curves that is consistent with rotational modulation based on the first two data releases obtained from sectors 1 to 4. Our search yielded $134$ high-probability candidate rotational variables -- $60$ of which have not been previously reported. Approximately half of these stars are identified in the literature as Ap (mCP) stars. Comparisons between the subsample of high-probability candidate rotationally variable Ap stars and the subsample of stars that are not identified as Ap reveal that the latter subsample exhibits statistically (i) shorter rotation periods and (ii) significantly lower photometric amplitudes.

In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previously published values. We also discuss the overall photometric phenomenology of the known magnetic massive and intermediate-mass stars and propose an observational strategy to augment this population by taking advantage of the high-quality observations produced by TESS.

Active galactic nuclei (AGNs) are considered to be one of the most appropriate sources of ultra high energy cosmic rays (UHECRs, $E \gtrsim 10^{18} \mathrm{~eV}$). Radiogalaxy Virgo A (M87) in the centre of a cluster of galaxies Virgo Cluster (VC) can be a prominent source of UHECRs. We investigate the possible contribution of Virgo A and the VC to the flux of events with trans-GZK energies - extremely high energy cosmic rays (EHECRs) - from the recent Auger and Telescope Array (TA) data sets ($E > 52 \mathrm{~EeV}$ and $E > 57 \mathrm{~EeV}$, respectively). We simulate EHECR propagation from Virgo A and the VC taking into account their deflections in galactic (GMF) and extragalactic (EGMF) magnetic fields and show that there is no excess of EHECR arrival directions from images of Virgo A/VC at different EHECR rigidities. By means of event-by-event analysis we recover the extragalactic arrival directions of EHECR events detected by Auger and TA for representative set of nuclei H(p), He, N, Si, Fe, and find evidences of enhanced fluxes of N-Si-Fe EHECRs from the Local Filament and Hot/Cold Spot regions. The Local Filament with its enhanced magnetic field is an expected contributor to the UHECR flux as the closest to the Earth last scattering centre, whereas Hot/Cold Spot region is a part of a larger arc-like spot, possibly created by diffusively spreading jet of UHECRs, accelerated in the relativistic jet of Virgo A during a prominent nuclear outburst about 10 - 12 Myr ago.

In this study we investigate jet interaction at a microscopic level in a cosmological environment, which responds to a key open question in the study of relativistic jets. Using small simulation systems during prior research, we initially studied the evolution of both electron-proton and electron-positron relativistic jets containing helical magnetic fields, by focusing on their interactions with an ambient plasma. Here, using larger jet radii, we have performed simulations of global jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). We found that the evolution of global jets strongly depends on the size of the jet radius. For example, phase bunching of jet electrons, in particular in the electron-proton jet, is mixed with larger jet radius due to the more complicated structures of magnetic fields with excited kinetic instabilities. In our simulation study these kinetic instabilities lead to new types of instabilities in global jets. In the electron-proton jet simulation a modified recollimation occurs and jet electrons are strongly perturbed. In the electron-positron jet simulation mixed kinetic instabilities occur at early times followed by a turbulence-like structure. Simulations using much larger (and longer) systems are further required in order to thoroughly investigate the evolution of global jets containing helical magnetic fields.

We perform particle-in-cell simulations of perpendicular nonrelativistic collisionless shocks to study electron heating and pre-acceleration for parameters that permit extrapolation to the conditions at young supernova remnants. Our high-resolution large-scale numerical experiments sample a representative portion of the shock surface and demonstrate that the efficiency of electron injection is strongly modulated with the phase of the shock reformation. For plasmas with low and moderate temperature (plasma beta $\beta_{\rm p}=5\cdot 10^{-4}$ and $\beta_{\rm p}=0.5$), we explore the nonlinear shock structure and electron pre-acceleration for various orientations of the large-scale magnetic field with respect to the simulation plane while keeping it at $90^\circ$ to the shock normal. Ion reflection off the shock leads to the formation of magnetic filaments in the shock ramp, resulting from Weibel-type instabilities, and electrostatic Buneman modes in the shock foot. In all cases under study, the latter provides first-stage electron energization through the shock-surfing acceleration (SSA) mechanism. The subsequent energization strongly depends on the field orientation and proceeds through adiabatic or second-order Fermi acceleration processes for configurations with the out-of-plane and in-plane field components, respectively. For strictly out-of-plane field the fraction of supra-thermal electrons is much higher than for other configurations, because only in this case the Buneman modes are fully captured by the 2D simulation grid. Shocks in plasma with moderate $\beta_{\rm p}$ provide more efficient pre-acceleration. The relevance of our results to the physics of fully three-dimensional systems is discussed.

A nonresonant cosmic-ray-current-driven instability may operate in the shock precursors of young supernova remnants and be responsible for magnetic-field amplification, plasma heating and turbulence. Earlier simulations demonstrated magnetic-field amplification, and in kinetic studies a reduction of the relative drift between cosmic rays and thermal plasma was observed as backreaction. However, all published simulations used periodic boundary conditions, which do not account for mass conservation in decelerating flows and only allow the temporal development to be studied. Here we report results of fully kinetic Particle-In-Cell simulations with open boundaries that permit inflow of plasma on one side of the simulation box and outflow at the other end, hence allowing an investigation of both the temporal and the spatial development of the instability. Magnetic-field amplification proceeds as in studies with periodic boundaries and, observed here for the first time, the reduction of relative drifts causes the formation of a shock-like compression structure at which a fraction of the plasma ions are reflected. Turbulent electric field generated by the nonresonant instability inelastically scatters cosmic rays, modifying and anisotropizing their energy distribution. Spatial CR scattering is compatible with Bohm diffusion. Electromagnetic turbulence leads to significant nonadiabatic heating of the background plasma maintaining bulk equipartition between ions and electrons. The highest temperatures are reached at sites of large-amplitude electrostatic fields. Ion spectra show supra-thermal tails resulting from stochastic scattering in the turbulent electric field. Together, these modifications in the plasma flow will affect the properties of the shock and particle acceleration there.

We study the interaction of relativistic jets with their environment, using 3-dimensional relativistic particle-in-cell simulations for two cases of jet composition: (i) electron-proton ($e^{-}-p^{+}$) and (ii) electron-positron ($e^{\pm}$) plasmas containing helical magnetic fields. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability and the Mushroom instability. We have found that these kinetic instabilities are suppressed and new types of instabilities can grow. For the $e^{-}-p^{+}$ jet, a recollimation-like instability occurs and jet electrons are strongly perturbed, whereas for the $e^{\pm}$ jet, a recollimation-like instability occurs at early times followed by kinetic instability and the general structure is similar to a simulation without a helical magnetic field. We plan to perform further simulations using much larger systems to confirm these new findings.

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

Active galactic nuclei (AGN) are considered as one of the most appropriate sources of cosmic rays with energy exceeding $\sim 10^{18}$ eV. Virgo A (M87 or NGC 4486) is the second closest to the Milky Way active galaxy. According to existing estimations it can be a prominent source of ultra high energy cosmic rays (UHECR). However not many events have been registered in the sky region near Virgo A, possibly due to magnetic field influence. In the present work we check UHECR events from the recent sets of data (AUGER, Telescope Array etc.) for possibility of their origination in this AGN. We carried out the simulation of UHECR motion from Virgo A taking into account their deflections in galactic (GMF) as well as extragalactic (EGMF) magnetic fields according to the several latest models. The maps of expected UHECR arrival directions were obtained as a result. It has been found the following: 1) UHECR deflection caused by EGMF is comparable with GMF one, moreover the influence of EGMF sometimes is dominating; 2) effect of EGMF demonstrates obvious asymmetry in the final distribution of expected UHECR arrival directions; 3) the results of simulations depend on chosen GMF model and are still open for the discussion.

The Workshop on results of the Project Kosmomikrofizyka-2 (Astroparticle Physics) of the National Academy of Sciences (NAS) of Ukraine "Astrophysical and cosmological problems of invisible mass and dark energy in the Universe" was held on November 21-22, 2012 in the Institute for Nuclear Research, Kyiv, Ukraine (http://lpd.kinr.kiev.ua/kmf12). This Project was carried out during three years (2010-2012) by scientists from various universities and institutes of the National Academy of Sciences of Ukraine; it was a logical continuation of the previous scientific program of the NAS of Ukraine "Researches of structure and composition of the Universe, hidden mass and dark energy (Kosmomikrofizyka)" in 2007-2009. These programs were devoted to theoretical and experimental investigations in astronomy, astrophysics, cosmology, physics of atomic nuclei and particle physics, which are related with the problems of dark matter and dark energy in the Universe.

The propagation of ultra high energy cosmic rays in Galactic and extragalactic magnetic fields is investigated in the present paper. The motion of charged particles of different energies and chemical composition is simulated using different Galactic magnetic field models. Positions for the real sources of events registered at the Auger observatory are calculated taking into account the influence of both Galactic and extragalactic turbulent fields. The possibility of their correlation with the Centaurus A radio galaxy is analyzed.