In this work, we study the merger of two neutron stars with a gravitational mass of 1.4Mmiddot in circle each, employing the Shen-Horowitz-Teige equation of state. This equation of state is a corner ...case, allowing the formation of a stable neutron star with the given total baryonic mass of 3.03Mmiddot in circle. We investigate in unprecedented detail the structure of the remnant, in particular the mass distribution, the thermal structure, and the rotation profile. We also compute fluid trajectories both inside the remnant and those relevant for the formation of the disk. We find a peanut-shaped fluid flow inside the remnant following a strong m=2 perturbation. Moreover, the flow is locally compressive, causing the appearance of dynamic hot spots. Further, we introduce new diagnostic measures that are easy to implement in numeric simulations and that allow one to quantify mass and compactness of merger remnants in a well-defined way. As in previous studies of supra- and hypermassive stars, we find a remnant with a slowly rotating core and an outer envelope rotating at nearly Keplerian velocity. We compute a Tolman-Oppenheimer-Volkoff star model which agrees well with that of the remnant in the core, while the latter possesses extensive outer layers rotating close to Kepler velocity. Finally, we extract the gravitational wave signal and discuss the detectability with modern observatories. This study has implications for the interpretation of gravitational wave detections from the postmerger phase and is relevant for short gamma-ray burst models.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
After the detection of gravitational waves caused by the coalescence of compact objects in the mass range of neutron stars, GW170817, several studies have searched for an imprint of tidal effects in ...the signal, employing different model assumptions. One important distinction is whether or not to assume that both objects are neutron stars and obey the same equation of state. Some studies assumed independent properties. Others assume a universal equation of state, and in addition that the tidal deformability follows certain phenomenological relations. An important question is whether the gravitational-wave data alone constitute observational evidence for finite tidal effects. All studies provide Bayesian credible intervals, often without sufficiently discussing the impact of prior assumptions, especially in the case of lower limits on the neutron-star tidal deformability or radius. In this article, we scrutinize the implicit and explicit prior assumptions made in those studies. Our findings strongly indicate that existing lower credible bounds are mainly a consequence of prior assumptions combined with information gained about the system's masses. Importantly, those lower bounds are typically not informed by the observation of tidal effects in the gravitational-wave signal. Thus, regarding them as observational evidence might be misleading without a more detailed discussion. Further, we point out technical problems and conceptual inconsistencies in existing studies. We also assess the limitations due to systematic waveform model uncertainties in a novel way, demonstrating that differences between existing models are not guaranteed to be small enough for an unbiased estimation of lower bounds on the tidal deformability. Finally, we propose strategies for gravitational-wave data analysis designed to avoid some of the problems we uncovered.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
We present general relativistic numerical simulations of binary neutron star (BNS) mergers with different initial spin configurations. We focus on models with stars of mass 1.4 M⊙ each, which employ ...the equation of state (EOS) by Shen, Horowitz, and Teige, and which result in stable NSs as merger remnants. For comparison, we consider two irrotational equal mass (M=1.35 M⊙) and unequal mass (M=1.29, 1.42 M⊙) BNS models using the APR4 EOS, which result in a supramassive merger remnant. We present visualizations of the fluid flow and temperature distribution and find a strong impact of the spin on vortex structure and nonaxisymmetric deformation. We compute the radial mass distribution and the rotation profile in the equatorial plane using recently developed measures independent of spatial gauge, revealing slowly rotating cores that can be well approximated by the cores of spherical stars. We also study the influence of the spin on the inspiral phase and the gravitational wave (GW) signal. Using a newly developed analysis method, we further show that gravitational waveforms from BNS mergers can exhibit one or more phase jumps after merger, which occur together with minima of the strain amplitude. We provide a natural explanation in terms of the remnant’s quadrupole moment, and show that cancellation effects due to phase jumps can have a strong impact on the GW power spectrum. Finally, we discuss the impact of the spin on the amount of ejected matter.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
The first multimessenger observation attributed to a merging neutron star binary provided an enormous amount of observational data. Unlocking the full potential of this data requires a better ...understanding of the merger process and the early postmerger phase, which are crucial for the later evolution that eventually leads to observable counterparts. In this work, we perform standard hydrodynamical numerical simulations of a system compatible with GW170817. We focus on a single equation of state and two mass ratios, while neglecting magnetic fields and neutrino radiation. We then apply newly developed postprocessing and visualization techniques to the results obtained for this basic setting. The focus lies on understanding the three-dimensional structure of the remnant, most notably the fluid flow pattern, and its evolution until collapse. We investigate the evolution of mass and angular momentum distribution up to collapse, as well as the differential rotation along and perpendicular to the equatorial plane. For the cases that we studied, the remnant cannot be adequately modeled as a differentially rotating axisymetric neutron star. Further, the dominant aspect leading to collapse is the gravitational wave radiation and not internal redistribution of angular momentum. We relate features of the gravitational wave signal to the evolution of the merger remnant and make the waveforms publicly available. Finally, we find that the three-dimensional vorticity field inside the disk is dominated by medium-scale disturbances and not the orbital velocity, with potential consequences for magnetic field amplification effects.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
Weak reactions are critical for the neutron richness of the matter dynamically ejected after the merger of two neutron stars. The neutron richness, defined by the electron fraction (Ye), determines ...which heavy elements are produced by the r-process and thus directly impacts the kilonova light curve. In this work, we have performed a systematic and detailed post-processing study of the impact of weak reactions on the distribution of the electron fraction and of the entropy on the dynamic ejecta obtained from an equal mass neutron star binary merger simulated in full general relativity and with microscopic equation of state. Previous investigations indicated that shocks increase Ye, however our results show that shocks can also decrease Ye, depending on their thermodynamical conditions. Moreover, we have found that neutrino absorption is key and need to be considered in future simulations. We also demonstrated that the angular dependence of the neutrino luminosity and the spatial distribution of the ejecta can lead to significant difference in the electron fraction distribution. In addition to the detailed study of the Ye evolution and its dependences, we have performed nucleosynthesis calculations. They clearly point to the necessity of improving the neutrino treatment in current simulations to be able to predict the contribution of neutron star mergers to the chemical history of the universe and to reliable calculate their kilonova light curves.
We present new results of fully general relativistic magnetohydrodynamic simulations of binary neutron star (BNS) mergers performed with the Whisky code. All the models use a piecewise polytropic ...approximation of the APR4 equation of state for cold matter, together with a 'hybrid' part to incorporate thermal effects during the evolution. We consider both equal and unequal-mass models, with total masses such that either a supramassive NS or a black hole is formed after merger. Each model is evolved with and without a magnetic field initially confined to the stellar interior. We present the different gravitational wave (GW) signals as well as a detailed description of the matter dynamics (magnetic field evolution, ejected mass, post-merger remnant/disk properties). Our simulations provide new insights into BNS mergers, the associated GW emission and the possible connection with the engine of short gamma-ray bursts (both in the 'standard' and in the 'time-reversal' scenarios) and other electromagnetic counterparts.
First measurements of the Collins and Sivers asymmetries of charged hadrons produced in deep-inelastic scattering of muons on a transversely polarized 6LiD target are presented. The data were taken ...in 2002 with the COMPASS spectrometer using the muon beam of the CERN SPS at 160 GeV/c. The Collins asymmetry turns out to be compatible with zero, as does the measured Sivers asymmetry within the present statistical errors.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
The COMPASS data acquisition system Fischer, H.; Franz, J.; Grunemaier, A. ...
IEEE transactions on nuclear science,
04/2002, Volume:
49, Issue:
2
Journal Article
Peer reviewed
A fully pipelined and massively parallel data acquisition system has been developed for the COMPASS experiment at CERN. The main requirements are to read 250000 detector channels at a trigger rate of ...up to 100 kHz. Such high rates are only possible when using a hit selection mechanism on the front-end combined with dead-time free readout. For this purpose, a time-to-digital converter (TDC) chip has been developed and is used for all time measurement applications in COMPASS. Distributed, field programmable gate array (FPGA)-based readout-driver modules handle parallel front-end initialization, synchronous trigger and control-signal distribution, and local event building at a processing speed of 160 Mbyte/s. Each of the 160 readout-driver modules connects to 16 front-end boards through independent twisted pair cables (CAT 7, 600 MHz) or optical fibers using an industrial (ESCON), self synchronizing link at 40 Mbyte/s. Automatic configuration through unique module and link identification ensures the flexibility and scalability to very large detector systems. The preprocessed data are transmitted through optical fibers at 160 Mbyte/s to the master event building system. Here global event building is realized on high-performance personal computers (PCs) and Gigabit ethernet network components. The complete events are sent to the central data recording at the CERN main site at an average rate of 40 Mbyte/s and stored in an object oriented database. A reduced system was set up for the commissioning run of COMPASS in 2000. Operation of the full system starts in July 2001.
Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations ...offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insufficient to fully capture the small-scale dynamo, initially triggered by the Kelvin-Helmholtz instability, and later sourced by several MHD processes involving differential rotation. Here, we alleviate this limitation by using explicit large-eddy simulations, a technique where the unresolved dynamics occurring at the sub-grid scales (SGS) is modeled by extra terms, which are functions of the resolved fields and their derivatives. The combination of high-order numerical schemes, high resolutions, and the gradient SGS model allow us to capture the small-scale dynamos produced during the binary neutron star mergers. Here we follow the first 50 milliseconds after the merger and, for the first time, we find numerical convergence on the magnetic field amplification, in terms of integrated energy and spectral distribution over spatial scales. We also find that the average intensity of the magnetic field in the remnant saturates at \(\sim 10^{16}\)~G around \(5\)~ms after the merger. After \(20-30\)~ms, both toroidal and poloidal magnetic field components grow continuously, fed by the winding mechanism that provides a slow inverse cascade. We find no clear hints for magneto-rotational instabilities, and no significant impact of the magnetic field on the redistribution of angular momentum in the remnant in our simulations, probably due to the very turbulent and dynamical topology of the magnetic field at all stages, with small-scale components largely dominating over the large-scale ones.
Weak reactions are critical for the neutron richness of the matter dynamically ejected after the merger of two neutron stars. The neutron richness, defined by the electron fraction (Ye), determines ...which heavy elements are produced by the r-process and thus directly impacts the kilonova light curve. In this work, we have performed a systematic and detailed post-processing study of the impact of weak reactions on the distribution of the electron fraction and of the entropy on the dynamic ejecta obtained from an equal mass neutron star binary merger simulated in full general relativity and with microscopic equation of state. Previous investigations indicated that shocks increase Ye, however our results show that shocks can also decrease Ye, depending on their thermodynamical conditions. Moreover, we have found that neutrino absorption are key and need to be considered in future simulations. We also demonstrated that the angular dependence of the neutrino luminosity and the spatial distribution of the ejecta can lead to significant difference in the electron fraction distribution. In addition to the detailed study of the Ye evolution and its dependences, we have performed nucleosynthesis calculations. They clearly point to the necessity of improving the neutrino treatment in current simulations to be able to predict the contribution of neutron star mergers to the chemical history of the universe and to reliable calculate their kilonova light curves.