Search for: All records

We report on the effects of strong magnetic fields on neutrino emission in the modified Urca process. We show that the effect of Landau levels on the various Urca pairs affects the neutrino emission spectrum and leads to an angular asymmetry in the neutrino emission. For low magnetic fields, the Landau levels have almost no effect on the cooling. However, as the field strength increases, the electron chemical potential increases resulting in a lower density at which Urca pairs can exist. For intermediate field strength, there is an interesting interference between the Landau level distribution and the Fermi distribution. For high enough field strength, the entire electron energy spectrum is eventually confined to a single Landau level producing dramatic spikes in the emission spectrum.

We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection frommore »the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

The local escape velocity provides valuable inputs to the mass profile of the galaxy, and requires understanding the tail of the stellar speed distribution. Following Leonard & Tremaine, various works have since modeled the tail of the stellar speed distribution as$\propto {(v_{\mathrm{esc}}-v)}^k$, wherev_escis the escape velocity, andkis the slope of the distribution. In such studies, however, these two parameters were found to be largely degenerate and often a narrow prior is imposed onkin order to constrainv_esc. Furthermore, the validity of the power-law form can breakdown in the presence of multiple kinematic substructures or other mis-modeled features in the data. In this paper, we introduce a strategy that for the first time takes into account the presence of kinematic substructure. We model the tail of the velocity distribution as a sum of multiple power laws as a way of introducing a more flexible fitting framework. Using mock data and data from FIRE simulations of Milky Way-like galaxies, we show the robustness of this method in the presence of kinematic structure that is similar to the recently discovered Gaia Sausage. In a companion paper, we present the new measurement of the escape velocity and subsequently the massmore »of the Milky Way using Gaia eDR3 data.

Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless, a full understanding of the interior configuration of the field remains elusive. In this work, we present general relativistic magnetohydrodynamics (MHD) simulations of the magnetic field evolution in neutron stars lasting ${\sim } {880}\,$ms (∼6.5 Alfvén crossing periods) and up to resolutions of $0.1155\,$km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes, and turbulence. We find that the purely poloidal setup generates a toroidal field, which later decays exponentially reaching $1{{\ \rm per\ cent}}$ of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20 per cent of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi-equilibrium state. An analysis of the spectra in our higher resolution setups reveals, however, that in most cases we are not observing turbulence atmore »small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.

Measuring the escape velocity of the Milky Way is critical in obtaining the mass of the Milky Way, understanding the dark matter velocity distribution, and building the dark matter density profile. In Necib & Lin, we introduced a strategy to robustly measure the escape velocity. Our approach takes into account the presence of kinematic substructures by modeling the tail of the stellar distribution with multiple components, including the stellar halo and the debris flow called the Gaia Sausage (Enceladus). In doing so, we can test the robustness of the escape velocity measurement for different definitions of the “tail” of the velocity distribution and the consistency of the data with different underlying models. In this paper, we apply this method to the Gaia eDR3 data release and find that a model with two components is preferred, although results from a single-component fit are also consistent. Based on a fit to retrograde data with two bound components to account for the relaxed halo and the Gaia Sausage, we find the escape velocity of the Milky Way at the solar position to be$v_{\mathrm{esc}}={445}_{-8}^{+25}$km s⁻¹. A fit with a single component to the same data gives$v_{\mathrm{esc}}={472}_{-12}^{+17}$km s⁻¹.more »Assuming a Navarro−Frenck−White dark matter profile, we find a Milky Way concentration of$c_{200}={19}_{-7}^{+11}$and a mass of$M_{200}={4.6}_{-0.8}^{+1.5}\times {10}^{11}{M}_{\odot }$, which is considerably lighter than previous measurements.