We present a systematic numerical relativity study of the mass ejection and the associated electromagnetic transients and nucleosynthesis from binary neutron star (NS) mergers. We find that a few $10^{-3}\, M_\odot$ of material are ejected dynamically during the mergers. The amount and the properties of these outflow depend on binary parameters and on the NS equation of state (EOS). A small fraction of these ejecta, typically ${\sim}10^{-6}\, M_\odot$, is accelerated by shocks formed shortly after merger to velocities larger than $0.6\, {\rm c}$ and produces bright radio flares on timescales of weeks, months, or years after merger. Their observation could constrain the strength with which the NSs bounce after merger and, consequently, the EOS of matter at extreme densities. The dynamical ejecta robustly produce second and third $r$-process peak nuclei with relative isotopic abundances close to solar. The production of light $r$-process elements is instead sensitive to the binary mass ratio and the neutrino radiation treatment. Accretion disks of up to ${\sim}0.2\, M_\odot$ are formed after merger, depending on the lifetime of the remnant. In most cases, neutrino- and viscously-driven winds from these disks dominate the overall outflow. Finally, we generate synthetic kilonova light curves and find that kilonovae depend on the merger outcome and could be used to constrain the NS EOS.

General-relativistic simulations of binary neutron star mergers with viscosity reveal a new outflow mechanism operating in unequal mass binaries on dynamical timescales and enabled by turbulent viscosity. These "viscous-dynamical" ejecta are launched during merger due to the thermalization of mass exchange streams between the secondary and the primary neutron star. They are characterized by asymptotic velocities extending up to ${\sim} 0.8\, c$, and have masses that depend on the efficiency of the viscous mechanism. Depending on the unknown strength of the effective viscosity arising from magnetohydrodynamics instabilities operating during merger, the overall mass of the dynamical ejecta could be enhanced by a factor of a few and the mass of the fast tail of the ejecta having asymptotic velocities $\geq 0.6\, c$ by up to four orders of magnitude. The radioactive decay of the expanding viscous-dynamical ejecta could produce bright kilonova transients with signatures of free neutron decay in the first hour and enhanced near infrared flux on a timescale of a few days. The synchrotron remnant produced by the interaction between the ejecta and the interstellar medium could also be significantly enhanced by viscosity. Such remnant could be detected in the case of GW170817 as a rebrightening of the radio signal in the next months to years.