Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRadiative MHD Simulations of Photon Bubbles in Radiation-Supported Magnetized Atmospheres of Neutron Stars with Isotropic Thomson Scattering
76
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A major uncertainty in the structure and dynamics of magnetized, radiation pressure dominated neutron star accretion columns in X-ray pulsars and pulsating ultraluminous X-ray sources is that they are thought to be subject to the photon bubble instability. We present the results of two dimensional radiation relativistic magnetohydrodynamic simulations of a non-accreting, static atmosphere to study the development of this instability assuming isotropic Thomson scattering in the slow diffusion regime that is relevant to neutron star accretion columns. Photon bubbles generally grow faster toward shorter wavelengths, until a maximum growth rate is achieved at the radiation viscosity length scale, which is generally quite small and requires high numerical resolution to simulate. We confirm the consistency between our simulation results and linear theory in detail, and show that the nonlinear evolution inevitably leads to collapse of the atmosphere with the higher resolution simulation collapsing faster due to the presence of shorter length scale nonlinear structures. At least in static atmospheres with horizontally periodic boundary conditions, this resolution dependence may make simulations of the nonlinear dynamics of photon bubble instability in neutron star accretion columns challenging. It remains to be seen whether these difficulties will persist upon inclusion of an accretion flow through the top and magnetically-confined horizontal boundaries through which photons can escape. Our results here provide a foundation for such future work.
rate research
Read More
We present the implementation of a radiative transfer solver with coherent scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD) simulations of stellar surface convection. The code is fully parallelized using MPI domain decomposition, which allows for large grid sizes and improved resolution of hydrodynamical structures. We apply the code to simulate the surface granulation in a solar-type star, ignoring magnetic fields, and investigate the importance of coherent scattering for the atmospheric structure. A scattering term is added to the radiative transfer equation, requiring an iterative computation of the radiation field. We use a short-characteristics-based Gauss-Seidel acceleration scheme to compute radiative flux divergences for the energy equation. The effects of coherent scattering are tested by comparing the temperature stratification of three 3D time-dependent hydrodynamical atmosphere models of a solar-type star: without scattering, with continuum scattering only, and with both continuum and line scattering. We show that continuum scattering does not have a significant impact on the photospheric temperature structure for a star like the Sun. Including scattering in line-blanketing, however, leads to a decrease of temperatures by about 350,K below log tau < -4. The effect is opposite to that of 1D hydrostatic models in radiative equilibrium, where scattering reduces the cooling effect of strong LTE lines in the higher layers of the photosphere. Coherent line scattering also changes the temperature distribution in the high atmosphere, where we observe stronger fluctuations compared to a treatment of lines as true absorbers.
Neutron stars are among the most fascinating astrophysical sources, being characterized by strong gravity, densities about the nuclear one or even above, and huge magnetic fields. Their observational signatures can be extremely diverse across the electromagnetic spectrum, ranging from the periodic and low-frequency signals of radio pulsars, up to the abrupt high-energy gamma-ray flares of magnetars, where energies of ~10^46 erg are released in a few seconds. Fast-rotating and highly magnetized neutron stars are expected to launch powerful relativistic winds, whose interaction with the supernova remnants gives rise to the non-thermal emission of pulsar wind nebulae, which are known cosmic accelerators of electrons and positrons up to PeV energies. In the extreme cases of proto-magnetars (magnetic fields of ~10^15 G and millisecond periods), a similar mechanism is likely to provide a viable engine for the still mysterious gamma-ray bursts. The key ingredient in all these spectacular manifestations of neutron stars is the presence of strong magnetic fields in their constituent plasma. Here we will present recent updates of a couple of state-of-the-art numerical investigations by the high-energy astrophysics group in Arcetri: a comprehensive modeling of the steady-state axisymmetric structure of rotating magnetized neutron stars in general relativity, and dynamical 3-D MHD simulations of relativistic pulsar winds and their associated nebulae.
We present a detailed investigation of atmospheres around accreting neutron stars with high magnetic field ($Bgtrsim 10^{12}$ G) and low luminosity ($Llesssim 10^{33}$ erg/s). We compute the atmospheric structure, intensity and emergent spectrum for a plane-parallel, pure hydrogen medium by solving the transfer equations for the normal modes coupled to the hydrostatic and energy balance equations. The hard tail found in previous investigations for accreting, non-magnetic neutron stars with comparable luminosity is suppressed and the X-ray spectrum, although still harder than a blackbody at the star effective temperature, is nearly planckian in shape. Spectra from accreting atmospheres, both with high and low fields, are found to exhibit a significant excess at optical wavelengths above the Rayleigh-Jeans tail of the X-ray continuum.
All the neutron star (NS) atmosphere models published so far have been calculated in the cold plasma approximation, which neglects the relativistic effects in the radiative processes, such as cyclotron emission/absorption at harmonics of cyclotron frequency. Here we present new NS atmosphere models which include such effects. We calculate a set of models for effective temperatures T_eff =1-3 MK and magnetic fields B sim 10^{10}-10^{11} G, typical for the so-called central compact objects (CCOs) in supernova remnants, for which the electron cyclotron energy E_{c,e} and its first harmonics are in the observable soft X-ray range. Although the relativistic parameters, such as kT_eff /(m_e c^2) and E_{c,e} /(m_e c^2), are very small for CCOs, the relativistic effects substantially change the emergent spectra at the cyclotron resonances, E approx sE_{c,e} (s=1, 2,...). Although the cyclotron absorption features can form in a cold plasma due to the quantum oscillations of the free-free opacity, the shape and depth of these features change substantially if the relativistic effects are included. In particular, the features acquire deep Doppler cores, in which the angular distribution of the emergent intensity is quite different from that in the cold plasma approximation. The relative contributions of the Doppler cores to the equivalent widths of the features grow with increasing the quantization parameter b_eff = E_{c,e}/kT_eff and harmonic number s. The total equivalent widths of the features can reach sim 150-250 eV; they increase with growing b_eff and are smaller for higher harmonics.
Pulsars are highly magnetized and rapidly rotating neutron stars. The magnetic field can reach the critical magnetic field from which quantum effects of the vacuum becomes relevant, giving rise to magnetooptic properties of vacuum characterized as an effective non linear medium. One spectacular consequence of this prediction is a macroscopic friction that leads to an additional contribution in the spindown of pulsars. In this paper, we highlight some observational consequences and in particular derive new constraints on the parameters of the Crab pulsar and J0540-6919.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
