Astrophysics

We discuss the acceleration and escape of secondary cosmic-ray (CR) nuclei, such as lithium, beryllium and boron, produced by spallation of primary CR nuclei like carbon, nitrogen, and oxygen accelerated at the shock in supernova remnants (SNRs) surrounded by the interstellar medium (ISM) or a circumstellar medium (CSM). We take into account the energy-dependent escape of CR particles from the SNR shocks, which is supported by gamma-ray observations of SNRs, to calculate the spectra of primary and secondary CR nuclei running away into the ambient medium. We find that if the SNR is surrounded by a CSM with a wind-like density distribution (i.e., n CSM ??r ?? ), the spectra of the escaping secondary nuclei are harder than those of the escaping primary nuclei, while if the SNR is surrounded by a uniform ISM, the spectra of the escaping secondaries are always softer than those of the escaping primaries. Using this result, we show that if there was a past supernova surrounded by a dense wind-like CSM ( ??.5? 10 ?? M ?? yr ?? ) which happened ??.6? 10 5 yr ago at a distance of ??.6 kpc , we can simultaneously reproduce the spectral hardening of primary and secondary CRs above ??00 GV that have recently been reported by AMS-02.

IceCube has observed a flux of cosmic neutrinos, with a "bump" in the energy range 10?�E/TeV??00 that creates a 3? tension with gamma-ray data from the Fermi satellite. This has been interpreted as evidence for a population of hidden cosmic-ray accelerators. We propose an alternative explanation of this conundrum on the basis of cold dark matter which decays into sterile neutrinos that after oscillations produce the bump in the cosmic neutrino spectrum.

It is widely believed that super-Eddington accretion flow can produce powerful outflow, but where it originates from and how much mass and energy are carried away to which directions? To answer to these questions, we newly perform a large-box, two-dimensional radiation hydrodynamic simulation, paying special attention lest the results should depend on adopted initial and boundary conditions. We could achieve a quasi-steady state in an unprecedentedly large range, r=2 r S - 600 r S (with r S being the Schwarzschild radius) from the black hole. The accretion rate onto the central 10 M ??black hole is M ? BH ??80 L Edd / c 2 , whereas the mass outflow rate is M ? outflow ??4 L Edd / c 2 (where L Edd and c are the Eddington luminosity and the speed of light, respectively). The ratio ( M ? outflow / M ? BH ??.14 ) is much less than those reported previously. By careful inspection we find that most of outflowing gas which reach the outer boundary originates from the region at R??40 r S , while gas at 140 r S - 230 r S forms failed outflow. Therefore, significant outflow occurs inside the trapping radius ??50 r S . The mechanical energy flux (or mass flux) reaches its maximum in the direction of ??15 ??( ??80 ??) from the rotation axis. The total mechanical luminosity is L mec ??.16 L Edd , while the isotropic X-ray luminosity varies from L ISO X ??.9 L Edd , (for a face-on observer) to ??.1 L Edd (for a nearly edge-on observer). The power ratio is L mec / L ISO X ??.05 - 0.08 , in good agreement with the observations of Ultra-Luminous X-ray sources surrounded by optical nebulae.

We report detection of PSR B0656 + 14 with the Gran Telescopio Canarias in narrow optical F657 , F754 , F802 , and F902 and near-infrared JH K s bands. The pulsar detection in the K s band extends its spectrum to 2.2 μ m and confirms its flux increase towards the infrared. We also present a thorough analysis of the optical spectrum obtained by us with the VLT. For a consistency check, we revised the pulsar near-infrared and narrow-band photometry obtained with the \textit{HST}. We find no narrow spectral lines in the optical spectrum. We compile available near-infrared-optical-UV and archival 0.3-20keV X-ray data and perform a self-consistent analysis of the rotation phase-integrated spectrum of the pulsar using unified spectral models. The spectrum is best fitted by the four-component model including two blackbodies, describing the thermal emission from the neutron star surface and its hot polar cap, the broken power-law, originating from the pulsar magnetosphere, and an absorption line near ??0.5 keV detected previously. The fit provides better constraints on the model parameters than using only a single spectral domain. The derived surface temperature is T ??NS =7.9(3)? 10 5 K. The intrinsic radius (7.8-9.9 km) of the emitting region is smaller than a typical neutron star radius (13km) and suggests a nonuniform temperature distribution over the star surface. In contrast, the derived radius of the hot polar cap is about twice as large as the `canonical' one. The spectrum of the nonthermal emission steepens from the optical to X-rays and has a break near 0.1 keV. The X-ray data suggest the presence of another absorption line near 0.3keV.

Keck-telescope spectrophotometry of the companion of PSR J1810+1744 shows a flat, but asymmetric light-curve maximum and a deep, narrow minimum. The maximum indicates strong gravity darkening near the L_1 point, along with a heated pole and surface winds. The minimum indicates a low underlying temperature and substantial limb darkening. The gravity darkening is a consequence of extreme pulsar heating and the near-filling of the Roche lobe. Light-curve modeling gives a binary inclination i=65.7+/-0.4deg. With the Keck-measured radial-velocity amplitude K_c=462.3+/-2.2km/s, this gives an accurate neutron star mass M_NS=2.13+/-0.04M_o, with important implications for the dense-matter equation of state. A classic direct-heating model, ignoring the L_1 gravitational darkening, would predict an unphysical M_NS>3M_o. A few other ``spider" pulsar binaries have similar large heating and fill factor; thus, they should be checked for such effects.

It is an open challenge to estimate systematically the physical parameters of neutron star interiors from pulsar timing data while separating spin wandering intrinsic to the pulsar (achromatic timing noise) from measurement noise and chromatic timing noise (due to propagation effects). In this paper we formulate the classic two-component, crust-superfluid model of neutron star interiors as a noise-driven, linear dynamical system and use a state-space-based expectation-maximization method to estimate the system parameters using gravitational-wave and electromagnetic timing data. Monte Carlo simulations show that we can accurately estimate all six parameters of the two-component model provided that electromagnetic measurements of the crust angular velocity, and gravitational-wave measurements of the core angular velocity, are both available. When only electromagnetic data are available we can recover the overall relaxation time-scale, the ensemble-averaged spin-down rate, and the strength of the white-noise torque on the crust. However, the estimates of the secular torques on the two components and white noise torque on the superfluid are biased significantly.

We present the results the photometric observations of the Type IIP supernova SN 2012aw obtained for the time interval from 7 till 371 days after the explosion. Using the previously published values of the photospheric velocities we've computed the hydrodynamic model which simultaneously reproduced the photometry observations and velocity measurements. The model was calculated with the multi-energy group radiation hydrodynamics code STELLA. We found the parameters of the pre-supernova: radius R=500 R ??, nickel mass M ( 56 Ni ) ??.06 M ??, pre-supernova mass 25 M ??, mass of ejected envelope 23.6 M ??, explosion energy E??? 10 51 erg. The model progenitor mass M=25 M ??significantly exceeds the upper limit mass M=17 M ??, obtained from analysis the pre-SNe observations. This result confirms once more that the 'Red Supergiant Problem' must be resolved by stellar evolution and supernova explosion theories in interaction with observations.

The properties of relativistic jets, their interaction with the ambient environment and particle acceleration due to kinetic instabilities are studied self-consistently with Particle-in-Cell (PIC) simulations. In this work we study how a relativistic electron-positron jet containing a helical magnetic field evolves by focusing on its interaction with the external ambient plasma. Particularly, 3D PIC simulations are performed using a longer simulation system than previous studies with an embedded helical magnetic field. An important key issue in this work is how such a magnetic field affects an electron-positron jet and how this excites kinetic instabilities such as the Weibel instability (WI), the kinetic Kelvin-Helmholtz instability (kKHI) and others by further focusing on how particles accelerate. We do find that kinetic instabilities along with generated magnetic turbulence are present and consequently accelerate particles. At the linear stage we observe recollimation-like features at the center of the simulated jet and later-on as the electron-positron jet evolves, the magnetic fields generated by the instabilities become untangled and reorganized into a new topology near the non-linear phase. We additionally report indications of reconnection near the end of the non-linear stage, before the magnetic-field becomes untangled, as electrons get accelerated by multiple magnetic islands in the jet. In the present study the untangled magnetic field becomes turbulent without any reformation as it happened in our previous study of an electron-proton jet, which we will use to additionally compare the present results, obtaining important insights about the nature of these phenomena applicable to high-energy astrophysical environments such as Active Galactic Nuclei jets and Gamma-ray bursts.

Synchrotron radio emission from non-relativistic jets powered by massive protostars has been reported, indicating the presence of relativistic electrons and magnetic fields of strength ~0.3-5 mG. We study diffusive shock acceleration and magnetic field amplification in protostellar jets with speeds between 300 and 1500 km/s. We show that the magnetic field in the synchrotron emitter can be amplified by the non-resonant hybrid (Bell) instability excited by the cosmic-ray streaming. By combining the synchrotron data with basic theory of Bell instability we estimate the magnetic field in the synchrotron emitter and the maximum energy of protons. Protons can achieve maximum energies in the range 0.04-0.65 TeV and emit gamma rays in their interaction with matter fields. We predict detectable levels of gamma rays in IRAS 16547-5247 and IRAS 16848-4603. The gamma ray flux can be significantly enhanced by the gas mixing due to Rayleigh-Taylor instability. The detection of this radiation by the Fermi satellite in the GeV domain and the forthcoming Cherenkov Telescope Array at higher energies may open a new window to study the formation of massive stars, as well as diffusive acceleration and magnetic field amplification in shocks with velocities of about 1000 km/s.

Strongly magnetized and fast rotating neutron stars are known to be efficient particle accelerators within their magnetosphere and wind. They are suspected to accelerate leptons, protons and maybe ions to extreme relativistic regimes where the radiation reaction significantly feeds back to their motion. In the vicinity of neutron stars, magnetic field strengths are close to the critical value of B c ??.4? 10 9 T and particle Lorentz factors of the order γ??10 9 are expected. In this paper, we investigate the acceleration and radiation reaction feedback in the pulsar wind zone where a large amplitude low frequency electromagnetic wave is launched starting from the light-cylinder. We design a semi-analytical code solving exactly the particle equation of motion including radiation reaction in the Landau-Lifshits approximation for a null-like electromagnetic wave of arbitrary strength parameter and elliptical polarization. Under conventional pulsar conditions, asymptotic Lorentz factor as high as 10 8 ??10 9 are reached at large distances from the neutron star. However, we demonstrate that in the wind zone, within the spherical wave approximation, radiation reaction feedback remains negligible.