The mean square displacement and instantaneous diffusion coefficient for different configurations of charged particles in stochastic motion are calculated by numerically solving the associated equations of motion. The method is suitable for obtaining accurate descriptions of diffusion in both intermediate and long time regimes. It is also appropriate for studying a variety of astrophysical configurations since it may incorporate microscopic physics that analytical methods cannot cope with. The results show that, in the intermediary time regime, the diffusion coefficient has an irregular behavior, which can be described in terms of the complex interplay appearing between the physical parameters describing the configuration. The main conclusion is that such an approach may serve at differential diagnosis of different astrophysical configurations.

The study of the Brownian motion of a charged particle in electric and magnetic fields fields has many important applications in plasma and heavy ions physics, as well as in astrophysics. In the present paper we consider the electromagnetic radiation properties of a charged non-relativistic particle in the presence of electric and magnetic fields, of an exterior non-electromagnetic potential, and of a friction and stochastic force, respectively. We describe the motion of the charged particle by a Langevin and generalized Langevin type stochastic differential equation. We investigate in detail the cases of the Brownian motion with or without memory in a constant electric field, in the presence of an external harmonic potential, and of a constant magnetic field. In all cases the corresponding Langevin equations are solved numerically, and a full description of the spectrum of the emitted radiation and of the physical properties of the motion is obtained. The Power Spectral Density (PSD) of the emitted power is also obtained for each case, and, for all considered oscillating systems, it shows the presence of peaks, corresponding to certain intervals of the frequency.

The possibility that dark matter may be in the form of a Bose-Einstein Condensate (BEC) has been extensively explored at galactic scale. In particular, good fits for the galactic rotations curves have been obtained, and upper limits for the dark matter particle mass and scattering length have been estimated. In the present paper we extend the investigation of the properties of the BEC dark matter to the galactic cluster scale, involving dark matter dominated astrophysical systems formed of thousands of galaxies each. By considering that one of the major components of a galactic cluster, the intra-cluster hot gas, is described by King's $\beta$-model, and that both intra-cluster gas and dark matter are in hydrostatic equilibrium, bound by the same total mass profile, we derive the mass and density profiles of the BEC dark matter. In our analysis we consider several theoretical models, corresponding to isothermal hot gas and zero temperature BEC dark matter, non-isothermal gas and zero temperature dark matter, and isothermal gas and finite temperature BEC, respectively. The properties of the finite temperature BEC dark matter cluster are investigated in detail numerically. We compare our theoretical results with the observational data of 106 galactic clusters. Using a least-squares fitting, as well as the observational results for the dark matter self-interaction cross section, we obtain some upper bounds for the mass and scattering length of the dark matter particle. Our results suggest that the mass of the dark matter particle is of the order of $\mu $eV, while the scattering length has values in the range of $10^{-7}$ fm.

We investigate the transition to Self Organized Criticality in a two-dimensional model of a flux tube with a background flow. The magnetic induction equation, represented by a partial differential equation with a stochastic source term, is discretized and implemented on a two dimensional cellular automaton. The energy released by the automaton during one relaxation event is the magnetic energy. As a result of the simulations we obtain the time evolution of the energy release, of the system control parameter, of the event lifetime distribution and of the event size distribution, respectively, and we establish that a Self Organized Critical state is indeed reached by the system. Moreover, energetic initial impulses in the magnetohydrodynamic flow can lead to one dimensional signatures in the magnetic two dimensional system, once the Self Organized Critical regime is established. The applications of the model for the study of Gamma Ray Bursts is briefly considered, and it is shown that some astrophysical parameters of the bursts, like the light curves, the maximum released energy, and the number of peaks in the light curve can be reproduced and explained, at least on a qualitative level, by working in a framework in which the systems settles in a Self Organized Critical state via magnetic reconnection processes in the magnetized Gamma Ray Burst fireball.

A simplified one dimensional grid is used to model the evolution of magnetized plasma flow. We implement diffusion laws similar to those so-far used to model magnetic reconnection with Cellular Automata. As a novelty, we also explicitly superimpose a background flow. The aim is to numerically investigate the possibility that Self-Organized Criticality appears in a one dimensional magnetized flow. The cellular automaton's cells store information about the parameter relevant to the evolution of the system being modelled. Under the assumption that this parameter stands for the magnetic field, the magnetic energy released by one volume during one individual relaxation event is also computed. Our results show that indeed in this system Self-Organized Criticality is established. The possible applications of this model to the study of the X-ray afterglows of GRBs is also briefly considered.

This paper reports on the statistical behaviour of the optical IntraDay Variability of BL Lac S5 0716+714. Available IntraDay Variability data in the optical is tested to see whether or not the magnitude is log-normally distributed. It was consistently found that this is not the case. This is in agreement with a previous discussion for data for the same object but in a different observational period. Simultaneously, the spectral slope of the light curves is calculated. The implications of these findings for models which discuss both the location and the source of IntraDay Variability are presented.

We consider a description of the stochastic oscillations of the general relativistic accretion disks around compact astrophysical objects interacting with their external medium based on a generalized Langevin equation with colored noise, which accounts for the general memory and retarded effects of the frictional force, and on the fluctuation-dissipation theorem. The presence of the memory effects influences the response of the disk to external random interactions, and modifies the dynamical behavior of the disk, as well as the energy dissipation processes. The generalized Langevin equation of the motion of the disk in the vertical direction is studied numerically, and the vertical displacements, velocities and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases. The Power Spectral Distribution (PSD) of the disk luminosity is also obtained. As a possible astrophysical application of the formalism we investigate the possibility that the Intra Day Variability (IDV) of the Active Galactic Nuclei (AGN) may be due to the stochastic disk instabilities. The perturbations due to colored/nontrivially correlated noise induce a complicated disk dynamics, which could explain some astrophysical observational features related to disk variability.

This paper aims to investigate/map the effects that perturbations applied to an accretion disk might produce on the registered Light Curves (LC). The case of accretion disks around supermassive active black holes (AGNs) is studied with the goal to explain some of the statistical properties of the observed IntraDay Variability (IDV). The region producing optical IDV is perturbed by allowing it to develop a mass density of a fractional Brownian Motion-like type. The light curves and spectral slopes are calculated and compared to observational data for different Hurst parameters. The spectral slopes of the simulated light curves vary in the range $(0.4,2.5)$. The agreement with observational data shows that a magnetized disk subjected to stochastic perturbations can produce some of the features observed in the light curves.

We consider a description of the stochastic oscillations of the general relativistic accretion disks around compact astrophysical objects based on the generalized Langevin equation, which accounts for the general retarded effects of the frictional force, and on the fluctuation-dissipation theorems. The vertical displacements, velocities and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases. The Power Spectral Distribution of the luminosity it is also obtained, and it is shown that it has non-standard values. The theoretical predictions of the model are compared with the observational data for the luminosity time variation of the BL Lac S5 0716+714 object.

The possibility that BL Lac S5 0716+714 exhibits a linear root mean square (rms)-flux relation in its IntraDay Variability (IDV) is analysed. The results may be used as an argument in the existing debate regarding the source of optical IDV in Active Galactic Nuclei. 63 time series in different optical bands were used. A linear rms-flux relation at a confidence level higher than 65% was recovered for less than 8% of the cases. We were able to check if the magnitude is log-normally distributed for eight timeseries and found, with a confidence > 95%, that this is not the case.

Boson clouds around black holes exhibit interesting physical phenomena through the Penrose process of superradiance, leading to black hole spin-down. Axionic clouds are of particular interest, since the axion Compton wavelength could be comparable to the Schwarzschild radius, leading to the formation of "gravitational atoms" with a black hole nucleus. These clouds collapse under certain conditions, leading to a "Bosenova". We model the dynamics of such unstable boson clouds by a simple cellular automaton and show that it exhibits self-organized criticality. Our results suggest that the evolution through the black hole Regge plane is due to self-organized criticality.

Once the temperature of a bosonic gas is smaller than the critical, density dependent, transition temperature, a Bose - Einstein Condensation process can take place during the cosmological evolution of the Universe. Bose - Einstein Condensates are very strong candidates for dark matter, since they can solve some major issues in observational astrophysics, like, for example, the galactic core/cusp problem. The presence of the dark matter condensates also drastically affects the cosmic history of the Universe. In the present paper we analyze the effects of the finite dark matter temperature on the cosmological evolution of the Bose-Einstein Condensate dark matter systems. We formulate the basic equations describing the finite temperature condensate, representing a generalized Gross-Pitaevskii equation that takes into account the presence of the thermal cloud in thermodynamic equilibrium with the condensate. The temperature dependent equations of state of the thermal cloud and of the condensate are explicitly obtained in an analytical form. By assuming a flat Friedmann-Robertson-Walker (FRW) geometry, the cosmological evolution of the finite temperature dark matter filled Universe is considered in detail in the framework of a two interacting fluid dark matter model, describing the transition from the initial thermal cloud to the low temperature condensate state. The dynamics of the cosmological parameters during the finite temperature dominated phase of the dark matter evolution are investigated in detail, and it is shown that the presence of the thermal excitations leads to an overall increase in the expansion rate of the Universe.

Observational data in the BVRI bands of the variable BL Lacertae Object S5 0716+714 is discussed from the point of view of its Power Spectral Distribution (PSD). A model of the type $P(f) = \beta f^{-1} [1 + (\frac{f} {\delta}) ^{\alpha -1}]^{-1} + \gamma $ is fitted to the data for four null hypothesis and the Bayesian $p$ parameter for the fits is calculated. Spectral slopes with values ranging from 1.083 to 2.65 are obtained, with medium values for each band of $\bar{\alpha}_B =2.028$, $\bar{\alpha}_V = 1.809$, $\bar{\alpha}_R = 1.932$ and $\bar{\alpha}_I = 1.54$ respectively. These values confirm conclusions of previous studies, namely that the source is turbulent. Two disk models, the standard prescription of the Shakura-Sunyaev disk and magnetized disks exhibiting MagnetoRotational Instability, were discussed. We found that it is unlikely that they explain this set of observational data.

Observational data in the BVRI bands of the variable BL Lacertae Object S5 0716+714 is discussed from the point of view of its Power Spectral Distribution (PSD). Analysis of the noise properties exhibited by the source. Assessment of the validity of a turbulent model. We fitted a model of the type $P(f) = \beta f^{-\alpha} + \gamma$ to the data for two null hypothesis and calculated the Bayesian p parameter for the fits. Different theoretical models are discussed and applied to these data and a special focus is placed on the MagnetoRotational Instability (MRI). The $P(f) = \beta f^{-\alpha}$ fit provides slopes with values ranging from 0.376 to 2.237, with medium values for each band of $\bar{\alpha}_B = 1.756$, $\bar{\alpha}_V = 1.704$, $\bar{\alpha}_R = 1.780$ and $\bar{\alpha}_I = 1.272$ respectively. An interval for the calculated effective $\alpha_{SS}$ parameter is obtained to be $[0.18,1.72]\cdot 10^{-2}$, corresponding to the order of magnitude previously inferred through both observational and numerical investigations. Structure function analysis, fractal dimension analysis and DFT analysis (all performed in other studies) plus our own analysis establish that the source is turbulent in a nontrivial way, i.e. there is "intrinsic" noise superimposed on the deterministic behavior of the source. We propose that embedding a weakly stochastic magnetic field component in the MRI (MagnetoRotational Instability) framework and taking into account stochastic reconnection of magnetic field lines might explain microvariability in AGN objects.

We analyze the general relativistic oscillations of thin accretion disks around compact astrophysical objects interacting with the surrounding medium through non-gravitational forces. The interaction with the external medium (a thermal bath) is modeled via a friction force, and a random force, respectively. The general equations describing the stochastically perturbed disks are derived by considering the perturbations of trajectories of the test particles in equatorial orbits, assumed to move along the geodesic lines. By taking into account the presence of a viscous dissipation and of a stochastic force we show that the dynamics of the stochastically perturbed disks can be formulated in terms of a general relativistic Langevin equation. The stochastic energy transport equation is also obtained. The vertical oscillations of the disks in the Schwarzschild and Kerr geometries are considered in detail, and they are analyzed by numerically integrating the corresponding Langevin equations. The vertical displacements, velocities and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases.