The LIGO/Virgo collaborations recently announced the detection of a likely binary neutron star merger, GW190425. The total mass of GW190425 is significantly larger than the masses of Galactic double neutron stars known through radio astronomy. This suggests that GW190425 formed differently from Galactic double neutron stars. We hypothesize that GW190425 formed via unstable "case BB" mass transfer. According to this hypothesis, the progenitor of GW190425 was a binary consisting of a neutron star and a ${\sim} 4-5 {M_\odot}$ helium star, which underwent a common-envelope process. Following the supernovae of the helium star core, a tight, eccentric, double neutron star was formed, which merged in ${\lesssim}$ 10 Myr. The helium star progenitor may explain the unusually large mass of GW190425, while the short time to merger may explain why we do not see similar systems in radio. In order to test this hypothesis, we measure the eccentricity of GW190425 using publicly available LIGO/Virgo data. We constrain the eccentricity at 10 Hz to be $e \leq 0.007$ with $90\%$ confidence. This result provides no evidence for or against the unstable mass transfer scenario because the binary is likely to have circularized to $e\lesssim10^{-4}$ by the time it entered the LIGO/Virgo band. Future detectors operating in lower frequency bands will enable us to discern the formation channel of mergers similar to GW190425 using eccentricity measurements.

The conventional wisdom, dating back to 2012, is that the mass distribution of Galactic double neutron stars is well-fit by a Gaussian distribution with a mean of $1.33 M_\odot$ and a width of $0.09 M_\odot$. With the recent discovery of new Galactic double neutron stars and GW170817, the first neutron star merger event to be observed with gravitational waves, it is timely to revisit this model. In order to constrain the mass distribution of double neutron stars, we perform Bayesian inference using a sample of 17 Galactic double neutron stars effectively doubling the sample used in previous studies. We expand the space of models so that the recycled neutron star need not be drawn from the same distribution as the non-recycled companion. Moreover, we consider different functional forms including uniform, single-Gaussian, and two-Gaussian distributions. While there is insufficient data to draw firm conclusions, we find positive support (a Bayes factor of 9) for the hypothesis that recycled and non-recycled neutron stars have distinct mass distributions. The most probable model---preferred with a Bayes factor of 29 over the conventional model---is one in which the recycled neutron star mass is distributed according to a two-Gaussian distribution and the non-recycled neutron star mass is distributed uniformly. We show that precise component mass measurements of $\approx 20$ double neutron stars are required in order to determine with high confidence (a Bayes factor of 150) if recycled and non-recycled neutron stars come from a common distribution. Approximately $60$ are needed in order to establish the detailed shape of the distributions.