Refine
Year of publication
Document Type
- Article (18)
- Doctoral Thesis (5)
- Postprint (2)
- Habilitation Thesis (1)
Is part of the Bibliography
- yes (26)
Keywords
- turbulence (26) (remove)
Institute
Magnetic fields play important roles in many astrophysical processes. However, there is no universal diagnostic for the magnetic fields in the interstellar medium (ISM) and each magnetic tracer has its limitation. Any new detection method is thus valuable. Theoretical studies have shown that submillimetre fine-structure lines are polarized due to atomic alignment by ultraviolet photon-excitation, which opens up a new avenue to probe interstellar magnetic fields. We will, for the first time, perform synthetic observations on the simulated three-dimensional ISM to demonstrate the measurability of the polarization of submillimetre atomic lines. The maximum polarization for different absorption and emission lines expected from various sources, including star-forming regions are provided. Our results demonstrate that the polarization of submillimetre atomic lines is a powerful magnetic tracer and add great value to the observational studies of the submilimetre astronomy.
The galactic interstellar medium is magnetized and turbulent. The magnetic field and turbulence play important roles in many astrophysical mechanisms, including cosmic ray transport, star formation, etc. Therefore, measurements of magnetic field and turbulence information are crucial for the proper interpretation of astronomical observations. Nonetheless, the magnetic field observation is quite challenging, especially, there is not universal magnetic tracer for diffuse medium. Moreover, the modelling of turbulence can be oversimplified due to the lack of observational tools to diagnose the plasma properties of the turbulence in the galactic interstellar medium. The studies presented in this thesis have addressed these challenges by bridging the theoretical studies of magnetic field and turbulence with numerical simulations and observations.
The following research are presented in this thesis. The first observational evidence of the novel magnetic tracer, ground state alignment (GSA), is discovered, revealing the three-dimensional magnetic field as well as 2 orders of magnitude higher precision comparing to previous observational study in the stellar atmosphere of the post-AGB 89 Herculis. Moreover, the application of GSA in the sub-millimeter fine-structure lines is comprehensively studied for different elements and with magnetohydrodynamic simulations. Furthermore, the influence of GSA effect on the spectroscopy is analyzed and it is found that measurable variation will be produced on the spectral line intensity and the line ratio without accounting for the optical pumping process or magnetic field.
Additionally, a novel method to measure plasma modes in the interstellar medium, Signatures from Polarization Analysis (SPA), is proposed and applied to real observations. Magneto-sonic modes are discovered in different types of interstellar medium. An explanation is provided for the long-standing mystery, the origin of γ-ray enhanced emission “Cygnus Cocoon”, based on the comparison between the outcome of SPA and multi-waveband observational data. These novel methods have strong potentials for broader observational applications and will play crucial roles in future multi-wavelength astronomy.
We study the damping processes of both incompressible and compressible magnetohydrodynamic (MHD) turbulence in a partially ionized medium. We start from the linear analysis of MHD waves, applying both single-fluid and two-fluid treatments. The damping rates derived from the linear analysis are then used in determining the damping scales of MHD turbulence. The physical connection between the damping scale of MHD turbulence and the cutoff boundary of linear MHD waves is investigated. We find two branches of slow modes propagating in ions and neutrals, respectively, below the damping scale of slow MHD turbulence, and offer a thorough discussion of their propagation and dissipation behavior. Our analytical results are shown to be applicable in a variety of partially ionized interstellar medium (ISM) phases and the solar chromosphere. The importance of neutral viscosity in damping the Alfvenic turbulence in the interstellar warm neutral medium and the solar chromosphere is demonstrated. As a significant astrophysical utility, we introduce damping effects to the propagation of cosmic rays in partially ionized ISM. The important role of turbulence damping in both transit-time damping and gyroresonance is identified.
The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, N(E) proportional to E-p for p greater than or similar to 2. However, in standard Fermi shock acceleration, a substantial fraction of the swept-up particles do not enter the acceleration process at all. Additionally, if acceleration is efficient, then the nonlinear back-reaction of accelerated particles on the shock structure modifies the shape of the nonthermal tail of the particle spectra. Both of these modifications to the standard synchrotron afterglow impact the luminosity, spectra, and temporal variation of the afterglow. To examine the effects of including thermal particles and nonlinear particle acceleration on afterglow emission, we follow a hydrodynamical model for an afterglow jet and simulate acceleration at numerous points during the evolution. When thermal particles are included, we find that the electron population is at no time well fitted by a single power law, though the highest-energy electrons are; if the acceleration is efficient, then the power-law region is even smaller. Our model predicts hard-soft-hard spectral evolution at X-ray energies, as well as an uncoupled X-ray and optical light curve. Additionally, we show that including emission from thermal particles has drastic effects (increases by factors of 100 and 30, respectively) on the observed flux at optical and GeV energies. This enhancement of GeV emission makes afterglow detections by future gamma-ray observatories, such as CTA, very likely.
The general electromagnetic fluctuation theory is a powerful tool to analyze the magnetic fluctuation spectrum of a plasma. Recent works utilizing this theory for a magnetized non-relativistic isotropic Maxwellian electron-proton plasma have demonstrated that the equilibrium ratio of vertical bar delta B vertical bar/B-0 can be as high as 10(-12). This value results from the balance between spontaneous emission of fluctuations and their damping, and it is considerably smaller than the observed value vertical bar delta B vertical bar/B-0 in the solar wind at 1 au, where 10(-3) less than or similar to vertical bar delta B vertical bar/B-0 less than or similar to 10(-1). In the present manuscript, we consider an anisotropic bi-Maxwellian distribution function to investigate the effect of plasma instabilities on the magnetic field fluctuations. We demonstrate that these instabilities strongly amplify the magnetic field fluctuations and provide a sufficient mechanism to explain the observed value of vertical bar delta B vertical bar/B-0 in the solar wind at 1 au.
An interesting feature of solar wind fluctuations is the occasional presence of a well-pronounced peak near the spectral knee. These peaks are well investigated in the context of magnetic field fluctuations in the magnetosheath and they are typically related to kinetic plasma instabilities. Recently, similar peaks were observed in the spectrum of ion flux fluctuations of the solar wind and magnetosheath. In this paper, we propose a simple analytical model to describe such peaks in the ion flux spectrum based on the linear theory of plasma fluctuations. We compare our predictions with a sample observation in the solar wind. For the given observation, the peak requires similar to 10 minutes to grow up to the observed level that agrees with the quasi-linear relaxation time. Moreover, our model well reproduces the form of the measured peak in the ion flux spectrum. The observed lifetime of the peak is about 50 minutes, which is relatively close to the nonlinear Landau damping time of 30-40 minutes. Overall, our model proposes a plausible scenario explaining the observation.
Due to the relatively low rate of particle-particle collisions in the solar wind, kinetic instabilities (e.g., the mirror and firehose) play an important role in regulating large deviations from temperature isotropy. These instabilities operate in the high beta(parallel to) > 1 plasmas, and cannot explain the other limits of the temperature anisotropy reported by observations in the low beta beta(parallel to) < 1 regimes. However, the instability conditions are drastically modified in the presence of streaming (or counterstreaming) components, which are ubiquitous in space plasmas. These effects have been analyzed for the solar wind conditions in a large interval of heliospheric distances, 0.3-2.5 AU. It was found that proton counter-streams are much more crucial for plasma stability than electron ones. Moreover, new instability thresholds can potentially explain all observed bounds on the temperature anisotropy, and also the level of differential streaming in the solar wind.
Stars are born in turbulent molecular clouds that fragment and collapse under the influence of their own gravity, forming a cluster of hundred or more stars. The star formation process is controlled by the interplay between supersonic turbulence and gravity. In this work, the properties of stellar clusters created by numerical simulations of gravoturbulent fragmentation are compared to those from observations. This includes the analysis of properties of individual protostars as well as statistical properties of the entire cluster. It is demonstrated that protostellar mass accretion is a highly dynamical and time-variant process. The peak accretion rate is reached shortly after the formation of the protostellar core. It is about one order of magnitude higher than the constant accretion rate predicted by the collapse of a classical singular isothermal sphere, in agreement with the observations. For a more reasonable comparison, the model accretion rates are converted to the observables bolometric temperature, bolometric luminosity, and envelope mass. The accretion rates from the simulations are used as input for an evolutionary scheme. The resulting distribution in the Tbol-Lbol-Menv parameter space is then compared to observational data by means of a 3D Kolmogorov-Smirnov test. The highest probability found that the distributions of model tracks and observational data points are drawn from the same population is 70%. The ratios of objects belonging to different evolutionary classes in observed star-forming clusters are compared to the temporal evolution of the gravoturbulent models in order to estimate the evolutionary stage of a cluster. While it is difficult to estimate absolute ages, the realtive numbers of young stars reveal the evolutionary status of a cluster with respect to other clusters. The sequence shows Serpens as the youngest and IC 348 as the most evolved of the investigated clusters. Finally the structures of young star clusters are investigated by applying different statistical methods like the normalised mean correlation length and the minimum spanning tree technique and by a newly defined measure for the cluster elongation. The clustering parameters of the model clusters correspond in many cases well to those from observed ones. The temporal evolution of the clustering parameters shows that the star cluster builds up from several subclusters and evolves to a more centrally concentrated cluster, while the cluster expands slower than new stars are formed.
Turbulence in the weakly collisional intracluster medium (ICM) of galaxies is able to generate strong thermal velocity anisotropies in the ions (with respect to the local magnetic field direction), if the magnetic moment of the particles is conserved in the absence of Coulomb collisions. In this scenario, the anisotropic pressure magnetohydrodynamic (AMHD) turbulence shows a very different statistical behaviour from the standard MHD one and is unable to amplify seed magnetic fields. This is in contrast to previous cosmological MHD simulations that are successful in explaining the observed magnetic fields in the ICM. On the other hand, temperature anisotropies can also drive plasma instabilities that can relax the anisotropy. This work aims to compare the relaxation rate with the growth rate of the anisotropies driven by the turbulence. We employ quasi-linear theory to estimate the ion scattering rate resulting from the parallel firehose, mirror and ion-cyclotron instabilities, for a set of plasma parameters resulting from AMHD simulations of the turbulent ICM. We show that the ICM turbulence can sustain only anisotropy levels very close to the instability thresholds. We argue that the AMHD model that bounds the anisotropies at the marginal stability levels can describe the Alfv,nic turbulence cascade in the ICM.
Rotating stellar convection transports angular momentum towards the equator, generating the characteristic equatorial acceleration of the solar rotation while the radial flux of angular momentum is always inwards. New numerical box simulations for the meridional cross-correlation < u(theta)u(phi)>, however, reveal the angular momentum transport towards the poles for slow rotation and towards the equator for fast rotation. The explanation is that for slow rotation a negative radial gradient of the angular velocity always appears, which in combination with a so-far neglected rotation-induced off-diagonal eddy viscosity term nu(perpendicular to) provides "antisolar rotation" laws with a decelerated equator Similarly, the simulations provided positive values for the rotation-induced correlation < u(r)u(theta)>, which is relevant for the resulting latitudinal temperature profiles (cool or warm poles) for slow rotation and negative values for fast rotation. Observations of the differential rotation of slowly rotating stars will therefore lead to a better understanding of the actual stress-strain relation, the heat transport, and the underlying model of the rotating convection.