Recently the AMS-02 reported the precise measurements of the energy spectra of medium-mass compositions (Neon, Magnesium, Silicon) of primary cosmic rays, which reveal different properties from those of light compositions (Helium, Carbon, Oxygen). Here we propose a nearby source scenario, together with the background source contribution, to explain the newly measured spectra of cosmic ray Ne, Mg, Si, and particularly their differences from that of He, C, O. Their differences at high energies can be naturally accounted for by the element abundance of the nearby source. Specifically, the abundance ratio of the nearby source to the background of the Ne, Mg, Si elements is lower by a factor of $\sim1.7$ than that of the He, C, O elements. Such a difference could be due to the abundance difference of the stellar evolution of the progenitor star or the acceleration process/environment, of the nearby source. This scenario can simultaneously explain the high-energy spectral softening features of cosmic ray spectra revealed recently by CREAM/NUCLEON/DAMPE, as well as the energy-dependent behaviors of the large-scale anisotropies. It is predicted that the dipole anisotropy amplitudes below PeV energies of the Ne, Mg, Si group are smaller than that of the He, C, O group, which can be tested with future measurements.

The spectral hardenings of cosmic ray nuclei above $\sim 200$ GV followed by softenings around 10 TV, the knee of the all-particle spectrum around PeV energies, as well as the pattern change of the amplitude and phase of the large-scale anisotropies around 100 TeV indicate the complexities of the origin and transportation of Galactic cosmic rays. It has been shown that nearby source(s) are most likely to be the cause of such spectral features of both the spectra and the anisotropies. In this work, we study the anisotropy features of different mass composition (or mass groups) of cosmic rays in this nearby source model. We show that even if the spectral features from the nearby source component is less distinctive compared with the background component from e.g., the population of distant sources, the anisotropy features are more remarkable to be identified. Measurements of the anisotropies of each mass composition (group) of cosmic rays by the space experiments such as DAMPE and HERD and the ground-based experiments such as LHAASO in the near future are expected to be able to critically test this scenario.

Imposing an external magnetic field in short-pulse intense laser-plasma interaction is of broad scientific interest in related plasma research areas. We propose a simple method using a virtual current layer by introducing an extra current density term to simulate the external magnetic field, and demonstrate it with three-dimensional particle-in-cell simulations. The field distribution and its evolution in sub-picosecond time scale are obtained. The magnetization process takes a much longer time than that of laser-plasma interaction due to plasma diamagnetism arising from collective response. The long-time evolution of magnetic diffusion and diamagnetic current can be predicted based on a simplified analytic model in combination with simulations.

The composition of the astrophysical relativistic jets remains uncertain. By kinetic particle-in-cell simulations, we show that the baryon component in the jet, or the so-called baryon loading effect (BLE), heavily affects relativistic jets transport dynamics in the interstellar medium. On the one hand, with the BLE, relativistic jets can transport in a much longer distance, because jet electrons draw a significant amount of energy from jet baryons via the Buneman-induced electrostatic waves and the Weibel-mediated collisionless shock; on the other hand, the jet electron phase space distribution may transform from a bottom-wide-single-peak structure to a center-wide-multiple-peak one by increasing the BLE, which largely influences the observed jet morphology. Implications for related astrophysical studies are also discussed.

With the standard deviation for the logarithm of the re-scaled range $\langle |F(t+\tau)-F(t)|\rangle$ of simulated fractal Brownian motions $F(t)$ given in a previous paper \cite{q14}, the method of least squares is adopted to determine the slope, $S$, and intercept, $I$, of the log$(\langle |F(t+\tau)-F(t)|\rangle)$ vs $\rm{log}(\tau)$ plot to investigate the limitation of this procedure. It is found that the reduced $\chi^2$ of the fitting decreases with the increase of the Hurst index, $H$ (the expectation value of $S$), which may be attributed to the correlation among the re-scaled ranges. Similarly, it is found that the errors of the fitting parameters $S$ and $I$ are usually smaller than their corresponding standard deviations. These results show the limitation of using the simple least square method to determine the dimension of a fractal time series. Nevertheless, they may be used to reinterpret the fitting results of the least square method to determine the dimension of fractal Brownian motions more self-consistently. The currency exchange rate between Euro and Dollar is used as an example to demonstrate this procedure and a fractal dimension of 1.511 is obtained for spans greater than 30 transactions.

To model a given time series $F(t)$ with fractal Brownian motions (fBms), it is necessary to have appropriate error assessment for related quantities. Usually the fractal dimension $D$ is derived from the Hurst exponent $H$ via the relation $D=2-H$, and the Hurst exponent can be evaluated by analyzing the dependence of the rescaled range $\langle|F(t+\tau)-F(t)|\rangle$ on the time span $\tau$. For fBms, the error of the rescaled range not only depends on data sampling but also varies with $H$ due to the presence of long term memory. This error for a given time series then can not be assessed without knowing the fractal dimension. We carry out extensive numerical simulations to explore the error of rescaled range of fBms and find that for $0<H<0.5$, $|F(t+\tau)-F(t)|$ can be treated as independent for time spans without overlap; for $0.5<H<1$, the long term memory makes $|F(t+\tau)-F(t)|$ correlated and an approximate method is given to evaluate the error of $\langle|F(t+\tau)-F(t)|\rangle$. The error and fractal dimension can then be determined self-consistently in the modeling of a time series with fBms.