\section{Introduction}\label{sec:intro}
\noindent Over 60 years ago, Burbidge {\it et al.} suggested that nuclear fission was responsible for the behavior of supernova light curves \cite{Burbidge56}. We now know the process that synthesizes fissioning nuclei---rapid neutron capture, or $r$-process, nucleosynthesis---is unrelated to supernova light curves and unlikely to occur robustly, if at all, in ordinary supernovae \cite{Ji+2016,Roederer+16}. The most attractive astrophysical site for an $r$ process that reaches fissioning nuclei is within a neutron star merger \cite{Lattimer+74,Meyer89,Freiburghaus+99,Goriely+11,Korobkin+12,Wanajo+14,Just+15}. The electromagnetic counterpart to the GW170817 neutron star merger \cite{Cowperthwaite2017,AbbottGW170817,Kasen} indicated some thousandths of a solar mass of lanthanides were produced in the event, possibly enough for mergers to account for all of the $r$-process lanthanides in the galaxy if the event was typical \cite{Cote}. If neutron star mergers are indeed the source of all $r$-process elements, including the actinides, we can look to such events as opportunities to probe fission properties.

Fission processes can play an important role in determining $r$-process observables such as abundance patterns and light curves. For example, lanthanide abundances can be influenced by late-time deposition of fission products \cite{Cote} and neutrons from fission can affect the amount of late-time neutron capture that sets the overall abundance pattern \cite{Eichler15}. Nuclear heating by fission can shape kilonova light curves \cite{Barnes+16}, with the late-time heating possibly dominated by the spontaneous fission of $^{254}$Cf \cite{Cfpaper}.  Understanding these effects requires knowledge of fission properties for hundreds of nuclei on the neutron-rich side of stability, about which little is experimentally known. Calculations of the $r$ process instead rely almost entirely on theoretical descriptions that vary widely. Here we examine the influence of two key fission inputs in $r$-process calculations: the fission fragment distributions and fission barrier heights (i.e.\ the maximum energy along the optimum path toward scission).

The importance of the fission fragment treatment in $r$-process calculations is well established \cite{Goriely+13,GorielyGEF,GorielyGMPGEF,Panov,Kodama,Eichler15,Eichler16}. Parameterized, semi-empirical formulae based on systematics, such as those in Refs. \cite{Panov,Kodama,ShibagakiMathews}, are an improvement over symmetric splits which assume the nucleus to divide in half, but are still a simplification of complex fission dynamics. Phenomenological descriptions, such as ABLA \cite{ABLA91,ABLA07}, Wahl \cite{Wahl}, and GEF (a GEneral description of Fission observables) \cite{GEF}, take into account the influence of shell structure, fission barriers, angular momentum, and neutron emission from the excited fragments. These fission yield descriptions have also been applied in $r$-process simulations \cite{Mendoza-Temis+15,Mendoza-Temis+16,RobertsWahl,GorielyGEF,GorielyGMPGEF}.

An aspect of the phenomenological descriptions of fission that has so far remained relatively unexplored in the $r$ process is the dependence of the fission yields on the excitation energy of the compound nuclei. This effect is generally assumed to be small, as the $r$-process temperatures are low ($\sim 0.1$ MeV), and is often ignored. The possibility for the yields of $r$-process nuclei to explicitly depend on the excitation energy was considered in Ref. \cite{PanovNIF}, but yield distributions were found to vary smoothly with initial energy, thus the yields at vanishing neutron bombarding energy were taken to be appropriate at all energies relevant for the $r$ process. Here we revisit an examination of the role of excitation energy on the fission yields of $r$-process nuclei by treating the three main processes occurring in a fission cycling $r$ process, namely neutron-induced, $\beta$-delayed, and spontaneous fission, with distinct initial excitation energies when determining the fission yields to apply to each process. To examine the impact of such an energy dependence on the $r$ process, we use fission yields from the publicly available GEF code (version GEF-2016-V1-2 \cite{GEFweb}).

We explore the sensitivity of the $r$ process to the assumptions made for the evolution of excited fission fragments by comparing results with GEF 2016 and results given an alternate treatment for the energy sharing and de-excitation. To do so we modified the published FREYA (Fission Reaction Event Yield Algorithm) code \cite{CPC,CPC_NVA} to use the GEF 2016 fission fragment yields (pre prompt neutron emission) as a function of post-scission fission fragment mass, charge and total kinetic energy, as input. The fission product yields (post prompt neutron emission) and average neutron multiplicities from the two codes are compared. We explore the effect of the additional neutron emission predicted by FREYA in the $r$ process for merger dynamical ejecta conditions.

Perhaps the most important characteristic of a heavy nucleus for predicting fission properties within a given model is the height of the fission barrier. It is well known that the fission barriers assumed for heavy, neutron-rich nuclei can have a great influence on the outcome of $r$-process calculations \cite{Petermann,GorielyGEF,GorielyGMPGEF,Eichler15,ShibagakiMathews,Samuel18}. Fission barriers often determine whether the nuclear flow ($\lambda_i Y_i$ where $\lambda_i$ is the rate of the reaction or decay and $Y_i$ is the abundance) will permit the population of nuclei of interest, such as the predicted superheavy island of stability \cite{Thielemann+83,PanovSF,GorielyGEF,GorielyGMPGEF,Petermann,BDFrp}. The influence of fission barriers can also lead to dramatically different conclusions regarding the origin of the second $r$-process peak as discussed in Ref. \cite{ShibagakiMathews}, where high fission barriers near $N=184$ coupled with very broad fission fragment distributions caused a disappearance of this main $r$-process feature in neutron star merger conditions.

Here we examine results using the Finite Range Droplet Model (FRDM2012) \cite{FRDM2012}, Thomas-Fermi (TF) \cite{TFmass}, Hartree-Fock-Bogoliubov (HFB-17) \cite{HFB17PRL}, and Extended Thomas-Fermi with Strutinsky Integral (ETFSI) \cite{ETFSI1992,ETFSI1995} model masses and corresponding barriers along with the GEF+FREYA yields described above. We explore how the different termination points for the $r$ process predicted by these models influence the final abundance pattern, and identify fissioning nuclei most accessed under a range of neutron star merger conditions. We also consider the population of $^{254}$Cf and find it to be sensitive to the fission barrier treatment, as suggested recently in Ref. \cite{WuBGMP}.

The paper is organized as follows: our application of GEF 2016 and implementation of FREYA is described in Section~\ref{sec:GEF}. In Section~\ref{sec:rpimpact} we explore the $r$-process impact of the fission yields and neutron multiplicities that result from this approach. Variations of fission barrier heights and other nuclear physics inputs are considered in Section~\ref{sec:rpMMastro}. In Section~\ref{sec:maxflows} we conclude by identifying the key fissioning nuclei important in a variety of neutron star merger conditions and common to all fission barrier models considered, some of which we find to be potentially within reach of future experimental facilities.

\section{GEF and FREYA fission treatments}\label{sec:GEF}

Before describing the fission properties of neutron-rich nuclei we obtain from the GEF code, we note this and other codes presently applied in the $r$ process, such as ALBA and Wahl, are phenomenological descriptions which extrapolate into unmeasured regions using systematics. Therefore, reminiscent of some nuclear mass models, predictions for the properties of neutron-rich nuclei are subject to large uncertainties. We look forward to progress in theoretical campaigns to calculate the fission yields of neutron-rich nuclei, such as density functional theory and microscopic-macroscopic approaches, in the future. Here we choose GEF 2016 as a tool to explore the sensitivity of the $r$ process to various yield properties due to its public availability and documented physics inputs.

We now describe the GEF and FREYA fission models, the differences between them, and how those differences affect the resulting fission product yields and average neutron multiplicities for nuclei relevant for the $r$ process.  Both codes use Monte Carlo techniques to produce fission events that provide complete kinematic information, including angular momentum, for all fission fragments; see Refs.~\cite{GEF,CPC_NVA} for more details.  (We note that two versions of GEF are available, the stand-alone Monte Carlo version used here and a deterministic subroutine for use in codes like TALYS \cite{TALYSRef,TALYScode}.)

While other codes such as FIFRELIN \cite{FIFRELIN} and CGMF \cite{Talou:2017qlc} are also available, GEF and FREYA are the fastest and thus most suitable for studies that require calculations for hundreds of nuclei, as is the case in this work.  Here we first describe the fission fragment yields, as well as their excitation energy dependence, followed by a discussion of particle emission in GEF and FREYA and how the two models can lead to different average neutron multiplicities, relevant for these studies.  All the calculations in this section, for both GEF and FREYA, are based on one million fission events generated for each fissioning isotope.

While GEF and FREYA can achieve similar end results as far as the output of complete fission events, the general approaches are rather different and worth some discussion.  The published version of FREYA is generally more limited in the number of isotopes available because of its approach. Like FIFRELIN and CGMF, FREYA requires the fission fragment yields and total kinetic energies (TKE) of the fragments in some form as inputs.  Thus the current published version of FREYA, FREYA 2.0.2 \cite{CPC_NVA}, is limited to certain isotopes: spontaneous fission of $^{244}$Cm, $^{252}$Cf, $^{238}$U and $^{238,240,242}$Pu as well as neutron-induced fission of $^{233,235,238}$U and $^{239,241}$Pu.

Although FREYA, FIFRELIN and CGMF differ in detail, the basic numerical approaches are similar.  The yields are sampled to choose one of the fragments with the partner chosen to conserve mass number $A$ and charge $Z$, followed by sampling of the TKE for the event.  Given the masses and TKE, the total excitation energy (TXE) is obtained and shared between the two fragments according to a model-dependent prescription.  Once this excitation energy sharing is complete, neutron evaporation follows. The excitation energy dependence of the input yields is modeled based on limited energy-dependent data.  Some progress has recently been made using models based on potential energy surfaces calculated in the macroscopic-microscopic approach \cite{RandrupMoller,PatrickJetal}.  Other yields calculated from many-body approaches based on density functional theory may become available in the future but these are so far quite limited \cite{schunck2016}.

The philosophy of the GEF code is quite different.  The only user inputs required by GEF 2016 are the charge and mass number of the fissioning nucleus, incident neutron energy or initial excitation energy for a compound nucleus, and the number of simulated events. The GEF yields and fragment excitation energies are then based on incident energy, nuclear mass and charge systematics.  This allows a broader range of isotopes to be studied and extrapolations made to regions of the nuclear chart where no data exist.  Thus GEF requires no input data for the yields and total kinetic energies, unlike all other available complete event fission codes. The fission fragment yields in GEF depend on three things: the fission barriers, the fission modes, and the excitation energy sharing.  We will touch upon all three here.

GEF 2016 employs the Thomas-Fermi (TF) macroscopic fission barriers \cite{TFBH} and nuclear masses \cite{TFmass}, along with a microscopic correction to the ground state mass and a pairing correction to the binding energy at the barrier.  Additionally, GEF makes use of experimentally-inferred fission barriers to derive further parameterized corrections to the TF barrier heights (see \cite{GEF}). The mass and charge systematics determined from these nuclear inputs are then extrapolated to unmeasured regions.

  There are four fission modes used in GEF 2016 to describe the fission fragment mass yields immediately after scission (called pre-neutron in GEF).
  Three of these were introduced by Brosa {\it et al.} \cite{Brosa}: standard 1 (S1, related to near-spherical nuclei near the doubly-magic closed proton and neutron shells at $Z = 50$ and $N=82$), standard 2 (S2, associated with the deformed neutron shell closure at $N=88$), and super long (SL, a symmetric mode).  The fourth is a super asymmetric mode (SA or S3) introduced by Mulgin {\it et al.} \cite{Mulgin}.

  The proximity to a closed neutron or proton shell determines whether the yields are symmetric or antisymmetric and also governs the widths of the yield distributions.  The effects of nuclear deformation, fragment angular momentum and charge polarization are taken into account.  Empirical $A$ and $Z$ systematics are used to determine the yields where no data are available.

To employ GEF 2016 to calculate the fission fragment yields for $r$-process nuclei, we extended the ``range of validity", defined by $A/Z > 172/80$ and $A/Z < 250/90$ for $76<Z<120$, implemented in the default GEF code. We do not, however, go beyond the nuclei included in the default mass and shell correction tables, $1\le Z \le136$ and $1\le N \le203$, to ensure that our output fission data is consistent with standard user outputs. We will later show that, given astrophysical conditions leading to fission, there is minimal $r$-process fission flow beyond $N=203$ for most of the nuclear mass models considered.  Therefore, implementing yields for nuclei with $N\le 203$ is sufficient.

Finally, the energy sharing mechanism in GEF determines how the intrinsic (statistical) excitation energy is divided between the heavy and light fragments.  Before scission, the fragments are coupled so that nucleons unpaired before scission are preferentially transferred to the heavy fragment.  Thus, at higher excitation energies, the transfer of unpaired nucleons to the heavy fragment gives the additional neutron multiplicity to the heavy fragment while the light fragment neutron multiplicity remains relatively constant.  The charge polarization, determining the charge yields, is assumed to remain essentially unchanged by an increase in excitation energy.  While all fission codes assume binary fission, FREYA assumes that neutrons are only emitted from the fully-accelerated fragments.  On the other hand, GEF has some probability for neutron emission prior to scission.  These ``scission neutrons'' are postulated to be emitted even in spontaneous fission.

Because the excitation energy is modeled at scission, no input TKE is required by GEF; instead, the TKE is fixed by energy conservation.  Recall that the TKE is an input in FREYA and other fission models. To then fix the fragment excitation energies at scission, GEF again applies $A$ and $Z$ systematics determined from data which are extrapolated to unmeasured isotopes.

The fission product yields (termed post-neutron emission in GEF) depend on the excitation energy and the identity of the compound nucleus.  In the case of spontaneous and $\beta$-delayed fission, the compound nucleus is simply the mass number $A$ while, for neutron-induced fission, it is $A+1$.  The excitation energy is zero for spontaneous fission. In the case of neutron-induced fission, it is the sum of the incident neutron energy, $E_n$, and the neutron separation energy, $S_n$, with $E_n + S_n$ on the order of a few MeV, except for the most neutron-rich nuclei near the neutron dripline. For $\beta$-delayed fission, the excitation energy (denoted here as $\langle E \rangle_\beta$) can be as large as $9$~MeV.

Prior evaluations of cumulative neutron-induced fission yields (after both prompt neutron emission and delayed emission from the $\beta$-decay of fission products) indicated that the yields become more symmetric with increasing incident neutron energy \cite{england1994endf}. However, a more recent experiment saw the cumulative fission yields for $^{239}$Pu(n,f) in particular to have a non-monotonic energy dependence on the energy of the incident neutron \cite{gooden2016}. Thus the energy dependence of the fission yields may be more complicated than previously assumed.

Neutron-induced fission in the $r$ process occurs late in time at low temperatures, corresponding to incident neutron energies between $\sim 0.01 - 0.2$~MeV, with most fission taking place at $\sim 0.1$~MeV. Although the $r$-process indeed occurs over a range of temperatures, and the neutrons are characterized by a thermal distribution rather than a single energy, the variation in fission yields over the energy range important for the $r$ process is small. Figure~\ref{fig:Yenergy94279} shows that for the incident neutron energy range of relevance to the $r$ process, even a very neutron-rich nucleus such as $^{279}$Pu with a low separation energy will not exhibit significant differences in its fission fragment yields. Percent-level differences begin to appear when comparing yields at  $\sim 0.1$~MeV and  $\sim 1.0$~MeV, however $1.0$~MeV ($\approx 10$~GK) is not a temperature at which fission participates in the $r$-process since here the environment is governed by nuclear statistical equilibrium (NSE) and has yet to synthesize fissioning nuclei.
The differences in the yields of a less neutron-rich nucleus with a separation energy on the order of MeV, such as $^{258}$Pu, are even less relevant. In this case, even a neutron incident energy of $\sim1.0$~MeV produces differences on the sub-percent level relative to $E_n = 0.1$~MeV. Therefore, unless the $r$-process proceeds through conditions in which nuclear reheating produces a rise in late time temperatures which reaches NSE values, the variance in the fission yields over the temperature evolution of the $r$ process can be safely ignored. Thus we apply a constant incident neutron energy of  $\sim 0.1$~MeV.

\begin{figure}
\begin{center}
\includegraphics[scale=0.485]{GEF_YNeutEnergy_94_279_zoom.pdf}
\end{center}
\caption{(Color online) The fission yield for $^{278}$Pu(n,f) with an incoming neutron energy of 0.1 MeV (black) as compared to 0.2 MeV (red), 0.5 MeV (blue), and 1.0 MeV (orange).}
\label{fig:Yenergy94279}
\end{figure}

\begin{figure*}
\begin{center}
\includegraphics[scale=0.43]{NZ_GEF_spf_nif_lblue.pdf}
\end{center}
\caption{(Color online) The GEF 2016 product yields for neutron-induced fission from a $0.1$~MeV neutron (red) compared to the product yields from spontaneous fission of the same compound nucleus (blue).}
\label{fig:yieldsNIF}
\end{figure*}

Even though the incident neutron energy of relevance to the $r$ process is low, the excitation energy of neutron-induced fission is typically a few MeV, considerably higher than the zero excitation energy of spontaneous fission. One could thus expect a difference in the spontaneous and neutron-induced fission yields. The resultant fission product yields are shown in Fig.~\ref{fig:yieldsNIF}. Some increased asymmetry is seen in the tails of the neutron-induced yields.
For the neutron-rich nuclei of interest here, the GEF 2016 systematics suggest a global trend of transition from asymmetric toward symmetric yields along most isotopic chains, with a region of primarily symmetric yields near the shell closure at $N=184$. Along an isotonic chain in this neutron-rich region, on average the yields become increasingly symmetric with increasing $Z$. This behavior can be further demonstrated by examining the trend in the width of the yields, represented in Fig.~\ref{fig:yieldwidthNIF} by the most probable mass number difference between the light and heavy fragment. To obtain this value, we find the maximum, $A'$, of the yield distributions shown in Fig.~\ref{fig:yieldsNIF}, which gives a most probable width of $|A' - (A_f - A' -\bar{\nu})|$ where $A_f$ and $\bar{\nu}$ are the mass number and average neutron multiplicity of the fissioning nucleus, respectively. Since, as can be seen in Fig.~\ref{fig:yieldsNIF}, many of the GEF 2016 yields in the neutron-rich regions contain both a symmetric and asymmetric component, a metric based on the maximum of the yields will not fully capture their complex behavior. However Fig.~\ref{fig:yieldwidthNIF} is still representative of the dominant yield trends predicted for $r$-process isotopes. In Section III, we will examine the impact that such asymmetric-to-symmetric yield trends have on the $r$-process abundance pattern and will show that the enhanced asymmetric yield contributions for finite excitation energies appear in key regions for a fission cycling $r$ process. We note that for some nuclei the fission yields predicted by other versions of the GEF code can be significantly different from those of GEF 2016. Nevertheless, the general arguments we lay out in this work remain the same.

\begin{figure}
\begin{center}
\includegraphics[scale=0.5]{GEF_width_nif_c.pdf}
\end{center}
\caption{(Color online) The most probable mass difference between the light ($A_L$) and heavy ($A_H$) fragments given GEF 2016 product yields for neutron-induced fission from a $0.1$~MeV neutron. The grey region shows the TF dripline.}
\label{fig:yieldwidthNIF}
\end{figure}

\begin{figure*}
\begin{center}
\includegraphics[scale=0.53]{GEF_FREYA_nubar_4pan_nc.pdf}
\end{center}
\caption{(Color online) The average prompt neutron multiplicity for fission induced by the capture of a 0.1 MeV neutron (right panels) compared to the multiplicity from spontaneous fission (left panels) using the energy sharing and de-excitation treatments of GEF 2016 (top panels) and FREYA (bottom panels).}
\label{fig:FREYAnu}
\end{figure*}

  We now turn to neutron emission from the fragments and the resulting fission product yields.  First we discuss how we modify FREYA to make use of the GEF 2016 fission fragment yields.  We then discuss how neutron emission differs in the two codes and how, even though we may start with identical yields in both GEF and FREYA, we may end up with different fission product yields and average neutron multiplicities, $\overline \nu$.  The difference in neutron emission is important for the $r$ process because prompt fission neutrons can be a substantial fraction of the late-time neutrons available for capture \cite{GorielyGEF}, thereby influencing the movement of lighter nuclei near the second and third $r$-process abundance peaks \cite{GorielyGMPGEF}.

For FREYA to be used in calculating the fission product yields and neutron emission relevant to the $r$ process, it needed to be adapted to use the GEF 2016 fission fragment yields as input.  FREYA was thus modified to take the fission fragment yields as function of mass and charge, $Y(A,Z)$, as well as the yields as a function of mass and TKE, $Y(A)$ and $Y({\rm TKE})$, given by GEF.  Thus, in the calculations that follow, the primary fragment yields employed in the two codes are identical and the differences are due to how neutron evaporation is treated.

Recall that, because GEF models the excitation energy partition between the fragments at scission, the energy available for neutron emission is set with no additional parameters required.  Neutron emission from the fully accelerated fragments in GEF then proceeds through a statistical model, using the relative neutron emission width from \cite{GEFGammaNref}.  The width depends on fragment mass, excitation energy, angular momentum, fragment temperature, and neutron separation energy, which are all modeled internally in GEF.

On the other hand, FREYA, like CGMF and FIFRELIN, starts with the TKE as input instead of the total excitation energy.  Thus the excitation energy partition between the fragments in these codes is done empirically, which requires a number of parameters.  FREYA 2.0.2 has five physics-based parameters: $d$TKE, which shifts TKE as a function of heavy fragment mass by some amount to ensure agreement with the average neutron multiplicity; $x$, representing additional excitation energy given to the light fragment, $x > 1$; $c$, setting the level of thermal fluctuations in the fragments; $e_0$, the value of the asymptotic level density parameter and $c_S$, governing the fragment spin magnitude, see Ref.~\cite{RVJR_gamma2}.
Neutron emission in FREYA proceeds assuming a Weisskopf-Ewing energy spectrum with a maximum temperature in the daughter fragment decreased appropriately \cite{CPC}.

Because FREYA does not have any systematic way of setting the parameters for unknown isotopes, we adjust $d$TKE in FREYA to match the known $\overline \nu$ in two specific cases, $^{252}$Cf(sf), for spontaneous fission, and $^{239}$Pu(n,f), for neutron-induced fission, while leaving the other FREYA parameters ($x$, $c$, $c_S$, and $e_0$) fixed at their default values for these isotopes.  In these two cases, it was found that the $d$TKE required in FREYA when using the GEF 2016 yields was negligibly small, and the calculated average neutron multiplicity in both models agreed with each other and with the data.  We use the $^{252}$Cf(sf) parameters for all spontaneous fission and the $^{239}$Pu(n,f) parameters for all neutron-induced and $\beta$-delayed fissions.

Because FREYA and GEF employ different methods of fragment de-excitation and the FREYA parameters are not fixed for each isotope, it is clear that we will not obtain the same neutron multiplicity from FREYA and GEF aside from our two matching points.  Indeed, there are no data for us to test the parameter values in either model.  However, we can take any differences between the fission product yields and neutron multiplicities in FREYA and GEF as an indication of the fission uncertainties affecting the $r$ process.

\begin{figure*}
\begin{center}
\includegraphics[scale=0.42]{NZ_GEF_FREYA_nif_pre.pdf}
\end{center}
\caption{(Color online) The GEF 2016 fragment yields for fission induced by a 0.1 MeV neutron (blue) along with the resultant product yields after applying the de-excitation treatments of GEF (red) and FREYA (black). The increase in the predicted prompt neutron emission shown in Fig.~\ref{fig:FREYAnu} shifts the FREYA yields to lower mass number.}
\label{fig:FREYAy}
\end{figure*}

  Figure~\ref{fig:FREYAnu} shows the resulting neutron multiplicities for spontaneous (left panels) and neutron-induced (right panels) fission from GEF 2016 (top panels) and FREYA (bottom panels).  In both cases the multiplicities are generally higher for neutron-induced fission, with $1-2$ more neutrons emitted than for spontaneous fission.  In the very neutron-rich region, in many cases the neutron multiplicities can be quite high, with some nuclei emitting $8-10$ neutrons, presumably because the outer neutrons are not strongly bound and thus emission is more probable.  It is also clear that more neutrons are emitted through the de-excitation process in FREYA than in GEF, even starting from the same initial fission fragment yields.  Both GEF and FREYA see a diagonal region of systematically higher neutron emission near $N=184$ which comes from the ability of the asymmetric yields seen in this region (recall Fig.~\ref{fig:yieldwidthNIF}) to emit more neutrons from a neutron-rich light fragment. We now describe where the differences in multiplicity may come from and discuss possible consequences of this difference later.

Recall that the model of excitation energy sharing in GEF is replaced in FREYA by the constant parameter $x$.  We note that in FIFRELIN and CGMF, the excitation energy sharing is also parameterized but, in those cases, mass dependent ratios are derived from data on the measured neutron multiplicity for specific isotopes as a function of fragment mass, $\overline \nu(A)$. Both FREYA's $x$ and the parameterized ratios in FIFRELIN and CGMF can reproduce the ``sawtooth" shape of $\overline \nu(A)$ in gross or fine detail, depending on which approach is employed.  The sawtooth shape is thought to arise from closed nuclear shells, in particular at $A = 132$ where there is a doubly closed shell and $\overline \nu(A)$ is at a minimum because closed-shell nuclei are harder to excite than highly deformed nuclei, resulting in fewer neutrons emitted. The $x$ parameter in FREYA gives more excitation energy to fragments with $A < 132$, resulting in a higher neutron multiplicity for the light fragment in spontaneous and thermal-neutron-induced fission, as suggested by data on well-studied nuclei such as $^{252}$Cf(sf), $^{235}$U(n,f) and $^{239}$Pu(n,f) \cite{Lemaire}.

In the case of very neutron rich nuclei, the structure of $\overline \nu(A)$ is wholly unknown.
While FREYA will minimize $\overline \nu(A)$ in the proximity of a closed shell, forcing $x > 1$ for very neutron-rich nuclei, where shell effects may be less important, would artificially increase neutron emission, as seen in Fig.~\ref{fig:FREYAnu}.  As suggested earlier, one can take this difference in neutron emission as a theoretical uncertainty.

The enhanced neutron multiplicities in FREYA relative to GEF 2016, due to the fixed energy sharing in FREYA, as seen in Fig.~\ref{fig:FREYAnu}, also has an effect on the shape of the fission product yields, as shown in Fig.~\ref{fig:FREYAy}.  The fission product yields are shifted toward lower mass numbers for FREYA.  The larger product yields on the low mass side of symmetry, especially for fragment masses $A<132$, is the result of larger neutron emission from the light fragments in FREYA.  We will discuss the sensitivity of the $r$ process to the energy sharing and de-excitation treatment in Sec. III. The sensitivity of the $r$ process depends on how much isotopic material reaches the most neutron-rich region with $N>184$ where FREYA predicts as many as $\sim 3$ more neutrons emitted per fission than GEF.

\begin{figure}
\begin{center}
\includegraphics[scale=0.48]{GEF_nubar_bdf_nc.pdf}
\end{center}
\caption{(Color online) The average excitation energy for the daughter nucleus populated by $\beta$-decay calculated as in \cite{BDFrp} with $\beta$-strength functions from \cite{MollerSd0} (upper panel) and the corresponding average prompt neutron emission as predicted by GEF 2016 when this excited daughter fissions (lower panel).}
\label{fig:enubdaught}
\end{figure}

$\beta$-delayed fission is a good test of the excitation energy dependence, as seen in Fig.~\ref{fig:enubdaught}. The excitation energy of this process can range from near zero, as is the case for spontaneous fission, up to $8-9$~MeV, as seen in the upper panel of Fig.~\ref{fig:enubdaught}, typically somewhat higher than the few MeV excitation energies of neutron-induced fission relevant for the $r$ process (the average excitation energies used to calculate $\beta$-delayed fission yields from GEF 2016, shown in Fig.~\ref{fig:enubdaught}, are tabulated in Supplemental Materials).  The average neutron multiplicity in $\beta$-delayed fission, shown in the lower panel of Fig.~\ref{fig:enubdaught}, correspondingly tends to be higher than that of neutron-induced fission shown in the upper right panel of Fig.~\ref{fig:FREYAnu}.

\begin{figure*}
\begin{center}
\includegraphics[scale=0.44]{NZ_GEFbdf_Kod.pdf}
\end{center}
\caption{(Color online) The product yields predicted for the fission of a daughter nucleus populated by $\beta$-decay after using the excitation energies shown in Fig.~\ref{fig:enubdaught} in GEF 2016 (purple) as compared to the yields from Kodama and Takahasi \cite{Kodama} (green).}
\label{fig:yieldcomp}
\end{figure*}

The $\beta$-delayed fission product yields predicted by GEF 2016, given the average excitation energy for the $\beta$ daughter shown in the upper panel of Fig.~\ref{fig:enubdaught}, also suggest that higher excitation energies will increase the asymmetry of the yields of neutron-rich nuclei.
Compared to the case of neutron-induced fission shown in Fig.~\ref{fig:yieldsNIF}, even larger deviations from symmetry can be observed in the GEF 2016 $\beta$-delayed yields shown in purple in Fig.~\ref{fig:yieldcomp}.  For reference, the parameterized double-Gaussian fission yields of Kodama and Takahashi \cite{Kodama}, widely used in $r$-process calculations, are also shown in Fig.~\ref{fig:yieldcomp} in green.  These yields also show a transition from asymmetric to symmetric fission near neutron number $N=184$, similar to GEF 2016.  However, in regions of the nuclear chart of interest for the $r$ process, the broad distributions predicted by Kodama and Takahashi are centered near $A\sim144$ for both asymmetric and symmetric fission while the asymmetric GEF 2016 yields tend to prefer $A\sim150$ daughters. Although both models predict asymmetric yield contributions for similar nuclei, note that the GEF 2016 yields for such nuclei often contain a symmetric component as well.

\section{Impact of Fission Yield Trends on $\bm{r}$-process abundances}\label{sec:rpimpact}

For nucleosynthesis calculations, we use the network Portable Routines for Integrated nucleoSynthesis Modeling (PRISM) developed jointly at the University of Notre Dame and Los Alamos National Laboratory \cite{BDFrp,Cote,Cfpaper}. PRISM permits a straightforward implementation of mass model-dependent nucleosynthesis rates due to its flexibility with nuclear data inputs. For the masses of neutron-rich nuclei, we first apply the Thomas-Fermi (TF) model in order to employ masses and fission barriers consistent with the GEF 2016 inputs used to determine the yields. For this, we explicitly use the barriers assumed in GEF 2016 which include corrections to TF barriers, as discussed in Sec. II. Where available we use experimental masses \cite{AME2016} as well as experimentally established half-lives and branching ratios from NUBASE \cite{NUBASE2016}. For theoretical $\alpha$-decay rates we use the well-established Viola-Seaborg formula \cite{VS1966}:

\begin{equation}
\text{log}_{10}T_{1/2}^{\alpha}(s) = \frac{aZ+b}{\sqrt{Q_{\alpha}(\text{MeV})}}+cZ+d+h_{\text{log}}
\end{equation}

\noindent where we apply a least-squares fit to NUBASE2016 half-life data that takes into account the reported experimental uncertainties when optimizing coefficients. Using this procedure we find values of $a=1.6606 \pm 0.0007$, $b=-9.2990 \pm 0.0656$, $c=-0.2121 \pm 0.0003$, and $d=-32.5432 \pm 0.0267$. Since the fit is performed using even-even nuclei \cite{VS1966}, species with unpaired nucleons are accounted for by finding the average difference between the fitted and experimental half-life values. These hindrance factors were found to be:
\begin{equation}
  h_{\text{log}}=
  \begin{cases}
    0.5325, & \text{$Z$ odd, $N$ even} \\
    0.5253, & \text{$Z$ even, $N$ odd} \\
    0.9222, & \text{$Z$ odd, $N$ odd.}
  \end{cases}
\end{equation}

\noindent We note, however, that the uncertainties in the treatment of theoretical $\alpha$-decay rates are relatively unimportant in the $r$ process since this decay is most influential at late times when the nuclei populated are within the experimentally probed regions. We use neutron capture, $\beta$-decay, neutron-induced fission and $\beta$-delayed fission rates as in Refs. \cite{Kawano+17,Mumpower+14,Mumpower+16,BDFrp,MollerQRPA}, with all rates determined from the same model masses as in Ref. \cite{Mumpower+15} and updated to be self-consistent with the fission barrier heights of a given model. For spontaneous fission we apply a parameterized prescription with a simple dependence on barrier height as in Refs. \cite{Karpov,Zagrebaev}. Therefore with the same fission barriers used to determine the fission yields and rates of all fission reaction and decay channels, our calculations which apply TF inputs represent the most fully self-consistent fission cycling $r$-process calculations in this work. Note that in this section, we focus on the influence of the fission yields in the $r$ process and therefore keep the fission rates identical in all comparisons shown.

\begin{figure*}
\begin{center}
\includegraphics[scale=0.48]{abA_Sneden_Ross1TFGEFvKodv5050_2panC_nc.pdf}
\end{center}
\caption{(Color online) The $r$-process abundances at 1 Gyr as a function of mass number (upper panel) and charge number (lower panel) using GEF 2016 yields for spontaneous fission, neutron-induced fission and $\beta$-delayed fission (red) as compared to using the fission yields of Kodama and Takahasi (orange) and simple symmetric splits (blue). The solar data is that of Sneden {\it et al.} \cite{Sneden}.}
\label{fig:abFFD}
\end{figure*}

\begin{figure*}
\begin{center}
\includegraphics[scale=0.585]{IntFlow_Ross1412traj1_tfGEF_ncf.pdf}
\end{center}
\caption{(Color online) The fission flow ($\lambda_i Y_i$ where $\lambda_i$ is the rate of the reaction or decay and $Y_i$ is the abundance) for nuclei undergoing neutron-induced (middle panel) and $\beta$-delayed (lower panel) fission integrated over time along with the resultant $r$-process abundance pattern (upper panel). The numbers labeling nuclei in middle and lower panels denote the mass number at the location of the maximum of the daughter fission yield distribution (here from GEF 2016).}
\label{fig:rtraj1}
\end{figure*}

For astrophysical conditions, we consider dynamical ejecta from a 1.2--1.4 $M_{\odot}$ neutron star merger simulation \cite{Rosswog,PiranRoss} for which we calculate nuclear reheating self-consistently with the chosen nuclear inputs when extrapolating beyond the reported simulation trajectory (with an assumed $50\%$ heating efficiency as in \cite{Oleg}; note heating efficiencies have been argued to lie between $25-100\%$ \cite{Metzger2010} and often a variety of efficiency values are considered \cite{Jonas}). We extrapolate the density by assuming free expansion as in Korobkin {\it et al.} \cite{Oleg} and assume NSE at 10 GK to obtain seed nuclei abundances using the SFHo equation of state \cite{SFHoSteiner}. Since these simulation trajectories are publicly available, we refer to the original number labeling (1--30) in order to permit direct comparisons with the results presented here. We first consider a ``cold", very neutron-rich ($Y_e=0.01957$) tidal tail trajectory in this set (traj.\ 1) which permits significant fission cycling due to its extreme neutron richness. Here the term ``cold" when used in reference to an astrophysical trajectory implies photodissociation drops out of equilibrium early, leaving $\beta$-decay to compete with neutron capture. In contrast, the term ``hot" when applied to trajectories implies conditions which support an extended (n,$\gamma$)$\rightleftharpoons$($\gamma$,n) equilibrium.

We show the impact of the asymmetric-to-symmetric yield trends predicted by GEF 2016 as compared to more simplistic descriptions in such fission rich environments in Fig.~\ref{fig:abFFD}. The underproduction at $A\sim144$ seen with a simple symmetric split (50/50) is a consequence of their narrow distribution which exclusively deposits material near $A\sim130$ when $r$-process material encounters the region with $Z\gtrsim90$ at $N<184$. In contrast, the fission yields of Kodama and Takahashi, whose $r$-process impact has been previously explored in Refs. \cite{Eichler15,Eichler16,BDFrp}, are exclusively asymmetric at $N<184$ and transition to symmetric distributions at higher neutron number (see Fig. \ref{fig:yieldcomp}). However the width of these fission yields places material over a fairly broad range around $A\sim144$ at both early and late times. The overproduction of the light lanthanides with the Kodama and Takahashi yields is a direct consequence of its global preference for $A\sim144$ fission daughters. The results obtained with the GEF 2016 yields here predict a more gradual fall off for the right edge of the second $r$-process peak, and the lanthanide abundances in this region, as well as near $A\sim150$, follow solar data trends. This lanthanide abundance behavior given the GEF 2016 yields is due to the transition from asymmetric to symmetric yields discussed in Sec. II.

\begin{figure*}
\begin{center}
\includegraphics[scale=0.51]{ab_GEFvFREYA_Ross22TF_pre_sfall_nc.pdf}
\end{center}
\caption{(Color online) An example of the $r$-process abundance pattern sensitivity to the treatment of excited fission fragments from examining results using GEF 2016 fragment yields (pre-neutron emission, blue), GEF product yields (post-neutron emission, red), FREYA product yields (black), and GEF product yields with zero excitation energy ($E^*=0$) (purple). The upper panel shows abundances at 1 Gyr while the lower panel shows the percent difference relative to results with the GEF 2016 product yields which account for excitation energy ($E^*\neq0$).}
\label{fig:abFREYA}
\end{figure*}

To explicitly demonstrate the role that GEF 2016 yield trends play in determining the $r$-process abundance pattern, the integrated fission flow can be cross checked with the fissioning daughter yields for these nuclei as shown in Fig.~\ref{fig:rtraj1}.  The fission flows are calculated with raw $r$-process abundances which satisfy $\sum Y_iA_i =1$ and no other scalings are applied so that the flows of different fission processes, such as neutron-induced and $\beta$-delayed fission, are directly comparable. In Fig.~\ref{fig:rtraj1}, we label each reactant species at $Z$, $N$ with the mass number which locates the maximum of the fission yield for the fissioning nucleus $Z_f$, $N_f$ (where $Z_f = Z+1$, $N_f = N-1$ for $\beta$df and $Z_f = Z$, $N_f = N+1$ for (n,f)). In the case of symmetric yields, the mass number, $A$, in the boxes of Fig. \ref{fig:rtraj1} will be $A\sim (N_f+Z_f)/2$, where $N_f$ and $Z_f$ correspond to the fissioning nucleus populated from the reactant species at $Z,N$. If the fission is asymmetric, the mass number shown is for one fragment peak so material will be deposited at $A'\sim (N_f+Z_f)-A$ as well. Of course the yields can contain both symmetric and asymmetric contributions that are similarly probable. However, it is still instructive to see explicitly where each fissioning nucleus will preferably deposit material. With the GEF 2016 yields, the $r$-process material in ``cold" dynamical ejecta trajectories, as in the very neutron-rich tidal tail ejecta conditions of Fig.~\ref{fig:abFFD}, first encounters symmetric yields centered at $A\sim144$ beyond the $N=184$ shell closure. Significant fission flow passes through this region as neutron-induced fission acts to terminate the $r$ process. As material decays back to stability, nuclei at neutron numbers below $N=184$ with substantial asymmetric yield contributions are encountered, placing material primarily near $A\sim110,130,$ and $150$. The contributions from fission products at early times near $A\sim144$, followed by late-time contributions near $A\sim130,150$, work to smooth the right edge of the second $r$-process peak. It is interesting to note that with the GEF 2016 yields there are only a few nuclei with $A\sim160$ fission daughters which place material directly into the rare-earth peak region. Therefore, with these yields it is not possible to explain rare-earth peak formation by a late-time deposition of fission material at $A\sim164$, as can occur with the four hump fragment distributions predicted by the SPY yield model \cite{GorielyGEF}.

We next turn to examine the influence of prompt fission neutrons in the $r$-process. Here we consider merger dynamical ejecta trajectory 22 with $Y_e\sim0.054$ which starts similarly cold to traj.\ 1 but later reaches higher temperatures. The most straight-forward approach to evaluate the impact of fission neutrons is to compare the $r$-process abundances using the GEF 2016 fragment yields (pre-neutron emission) and GEF 2016 product yields (post-neutron emission) as is shown in Fig.~\ref{fig:abFREYA}. The effect of prompt neutron emission on the $r$-process abundance pattern is two fold. The widening of the yield distributions toward lower mass number that accompanies neutron emission places more material to the left of the $N=82$ shell closure. This can lead to an increase on the order of $10\%$ in the height of the second peak. The increase of material held at the $N=82$ shell closure effectively decreases the number of available isotopes that can neutron capture back up past the $N=126$ shell closure which therefore reduces overall actinide abundances. This tendency to decrease the ability of material to access the heaviest nuclei also implies less fission activity in the region of mostly asymmetric yields between $N=126$ and $N=184$ which leads to less lanthanide material at $A\sim150$. The influence of the extra neutrons alone can be seen from the overall shift in the pattern to the right of the rare-earth peak toward higher mass number. Such a narrowing of the third $r$-process peak due to late-time neutron capture from $\beta$-delayed neutron emission and prompt neutron emission from excited fission fragments has been noted in previous work \cite{Goriely+13,Oleg,Eichler15,GorielyGMPGEF}.

The influence of the fission fragment energy sharing and de-excitation treatment, which determines prompt neutron emission, is also demonstrated in Fig.~\ref{fig:abFREYA} by the $r$-process abundances using the FREYA fission product yields as compared to those from GEF 2016. The widening of the yields toward lower mass number (recall Fig.~\ref{fig:FREYAy}) is primarily responsible for the differences seen between GEF and FREYA, as can be seen from the increase in height of the second peak produced by FREYA yields which place more material below $N=82$ relative to GEF 2016 yields.
Although the treatment of prompt neutron emission can lead to differences in the main $r$-process peaks, these effects are modest relative to the influence of the global yield trend demonstrated in Fig.~\ref{fig:abFFD}, which points to the fragment yields (prior to neutron emission) as being of primary importance for $r$-process calculations.

The finding that the location of fission fragment deposition is of dominant importance over prompt neutron emission is further demonstrated in Fig.~\ref{fig:abFREYA} by the effect of using excitation-energy dependent sets of fission yields for neutron-induced and $\beta$-delayed fission as compared to applying the zero excitation energy yields of spontaneous fission for all fission processes. The enhanced asymmetries in the fission yields discussed in Sec. II (and shown in Fig. \ref{fig:yieldsNIF}) due to an excited compound parent nucleus places more material at $A\sim110$ and $150$ as compared to the case when all fission processes make use of spontaneous fission yields. These enhanced asymmetric contributions produce a clear signature in the final abundance pattern with the fission yields which account for excitation energy giving a result more consistent with observed solar data. Meanwhile, the differing prompt neutron emission between these cases (recall from Figs.~\ref{fig:FREYAnu} and \ref{fig:enubdaught} that neutron-induced and $\beta$-delayed yields will contribute slightly more prompt neutrons than spontaneous fission) shows no clear effect on the abundances in the third peak region where late-time neutrons shape the final results.

We next consider the effect of more realistic fission yields and prompt neutron emission, as calculated by GEF 2016 and FREYA, on the nuclear heating rates needed to calculate kilonova light curves. For such heating rate calculations, we use the ``cold" dynamical ejecta conditions applied in Figs.~\ref{fig:abFFD} and \ref{fig:rtraj1}. We calculate the heating contributions from $\beta$-decay, neutron capture, $\beta$-delayed fission, and neutron-induced fission as flow $\times$ $Q$-value for each channel and compare the results using GEF and FREYA to simpler treatments in Fig.~\ref{fig:heatyield}.
We first focus on the heating before $\sim 1$ day when $\beta$-delayed and neutron-induced fission are most active.  Both the $\beta$-decay and $\beta$-delayed fission heating rates show some dependence on the yield distribution, however it is the neutron capture reaction channels that show a rather pronounced sensitivity to the fission yields. Applying simple symmetric splits produces early time heating contributions from neutron capture and neutron-induced fission which are about three orders of magnitude lower than those predicted using the GEF 2016 yields. This difference comes from the tendency of symmetric splits to deposit material in a concentrated region near the $N=82$ shell closure as compared to the wider distribution of daughter products seen with GEF 2016. The late-time addition of nuclei to the right of the $N=82$ shell closure permits further neutron capture and therefore more material is driven up in mass number toward the fissioning regions near $N=184$, correspondingly increasing the heating contribution from neutron capture and neutron-induced fission. A comparison of heating rates when using the yield model of Kodama and Takahashi, which does not contain a prescription for neutron emission from excited daughter fragments, confirms that it is the narrow placement of daughter nuclei near $N=82$ that is most responsible for the lower heating rate for neutron-induced processes in the symmetric yield case. The inclusion of prompt neutron emission from the excited fragments as in GEF can further increase the heating from neutron capture reaction channels by a factor of around three as compared to Kodama and Takahashi. We note that applying the yields obtained using the FREYA energy sharing and de-excitation treatment, while not shown in Fig.~\ref{fig:heatyield}, can further increase the heating from neutron capture processes by roughly $10\%$ at these early times.

 \begin{figure}
\begin{center}
\includegraphics[scale=0.49]{HeatingGEFKod5050_Ross1TFC_1pan.pdf}
\end{center}
\caption{(Color online) The nuclear heating rates for $\beta$-decay (blue), $\beta$-delayed fission (light blue), neutron capture (red) and neutron-induced fission (orange) as a function of time using GEF 2016 fission yields (solid lines) as compared to results using simple symmetric splits (dashed lines) and the yields of Kodama and Takahashi (dotted lines).  The masses and fission barriers applied are those of TF and GEF 2016 respectively.}
\label{fig:heatyield}
\end{figure}

When timescales on the order of days or longer are considered, the influence of the fission treatment on nuclear heating becomes more pronounced due to the late-time dominance of fission, specifically the spontaneous fission of $^{254}$Cf \cite{Cfpaper}. Fig.~\ref{fig:spfheat} shows the heating rate for $\beta$-decay, spontaneous fission, and $\alpha$-decay when the yield distributions from GEF 2016 and FREYA are compared with those obtained using simple symmetric splits. The 50/50 splits result in the largest predicted late-time dominance of the spontaneous fission heating with the greatest deviation between $\beta$-decay and spontaneous fission heating curves. The spread of fission recycled material produced by the GEF yields can produce nuclei closer to the neutron dripline where $\beta$-decay rates are faster, thereby increasing the total effective $\beta$ heating. When considering the extra late-time neutrons and increased yield widths predicted by FREYA, Fig.~\ref{fig:spfheat} shows that $\alpha$-decay and spontaneous fission heatings are decreased relative to the GEF case with slightly more material getting stuck near $N=82$ and therefore less material populating the highest mass number regions (as seen in Fig.~\ref{fig:abFREYA}). Figures~\ref{fig:heatyield} and \ref{fig:spfheat} explicitly demonstrate that in fission cycling conditions the exact details of the heating rates applied to calculations of kilonova light curves depend on the nuclear physics assumptions for the heaviest, fissioning $r$-process nuclei.

 \begin{figure}
\begin{center}
\includegraphics[scale=0.49]{LATEheatingGEFFREYA5050_Ross1TF_1pan.pdf}
\end{center}
\caption{(Color online) The nuclear heating rate for $\beta$-decay (blue), alpha decay (black) and spontaneous fission (green) as a function of time using GEF 2016 fission yields (solid lines) as compared to results using simple symmetric splits (dashed lines) and the yields with the FREYA energy sharing and de-excitation treatment (dotted lines). The masses and fission barriers applied are those of TF and GEF 2016 respectively.}
\label{fig:spfheat}
\end{figure}

As described above, exactly how fission yields shape abundance patterns and heating curves depends sensitively on which nuclei are fissioning. This is in part set by the $r$-process astrophysical conditions. To examine this sensitivity, we consider astrophysical conditions from the same 1.2--1.4 M$_{\odot}$ neutron star merger simulation \cite{Rosswog} used throughout this section. While all thirty of these dynamical ejecta trajectories are similarly neutron-rich, with $Y_e$ ranging from $\sim0.015-0.055$, they exhibit a variety of density and temperature profiles. A comparison of the final abundances (upper panel), abundance wighted mass number (middle panel), and temperature profile (lower panel) for the two trajectories found to represent the extremes in temperature evolution can be seen in Fig.~\ref{fig:rtraj117}. In contrast with the ``cold" conditions of trajectory 1, in the ``hot" dynamical ejecta conditions of trajectory 17, the $r$-process path does not significantly populate the mostly symmetric yield region past $N=184$ and therefore underproduces near $A\sim140$ relative to results with cold conditions. This is partially due to the ability of photodissociation to prevent material from reaching the most neutron-rich nuclei past $N=184$ but also due to the slightly higher $Y_e$ ($\sim0.049$) producing somewhat lower fission flow (total integrated fission flow for traj.\ 17 of 0.00488 as compared to 0.00567 for traj.\ 1). The region below $N=184$ with mostly asymmetric yields, however, is still accessed in such hot conditions resulting in increases to the abundances near $A\sim100,150$ (to show this behavior explicitly the version of Fig. \ref{fig:rtraj1} given the astrophysical conditions of traj.\ 17 used in Fig. \ref{fig:rtraj117} has been included in Supplemental Materials). Therefore which fissioning nuclei are accessed is influenced by the impact of the astrophysical conditions on the location and termination point of the $r$-process path.

\begin{figure}
\begin{center}
\includegraphics[scale=0.545]{abA_trajr1r17_tfC_ncsb.pdf}
\end{center}
\caption{(Color online) The $r$-process abundance pattern (upper panel) given two dynamical ejecta trajectories from a 1.2--1.4 M$_{\odot}$ neutron star merger simulation \cite{Rosswog} (lower panel) along with the abundance weighted mass number as a function of time for each set of conditions (middle panel).}
\label{fig:rtraj117}
\end{figure}

\section{Variations in other nuclear inputs influencing fission deposition}\label{sec:rpMMastro}

\begin{figure*}
\begin{center}
\includegraphics[scale=0.53]{BHrp_4pan_Ross1C_nc.pdf}
\end{center}
\caption{(Color online) Fission barrier heights for GEF 2016 (TF) (upper left), FRLDM (upper right), ETFSI (lower left) and HFB-14 (lower right) models along with a snapshot of the $r$-process material with an abundance $\ge 10^{-10}$ (black outline) just before the $r$-process path begins to move back toward stability in the cold, dynamical ejecta conditions of trajectory 1.}
\label{fig:bhrp}
\end{figure*}

Having established that the location and termination point of the $r$-process path influences which fissioning nuclei will most impact the $r$ process, we next consider how other nuclear inputs, such as the nuclear masses, fission barriers, and $\beta$-decay rates, can affect the nuclear flow. The connection between the path termination point and the fission yields most influencing $r$-process abundances has been made previously. For example, in Ref. \cite{ShibagakiMathews} the authors find that termination near $N=184$ along with symmetric yields deposits daughter nuclei near the second $r$-process peak while higher mass number termination points see more influence from their broad yield distributions for high mass nuclei. Here we explore such considerations by examining the termination behavior predicted by different mass models. We show that it is not only how far the $r$-process path proceeds, but also the structure of predicted fission barriers near $N=184$, which determine the fission yields of most relevance.

To study these dependencies, we compare results when employing the Finite Range Droplet Model (FRDM2012), Thomas-Fermi (TF), Hartree-Fock-Bogoliubov (HFB-17), and Extended Thomas-Fermi with Strutinsky Integral (ETFSI) mass models. We chose these four models since the data for both the masses and fission barriers consistent with model masses are publicly available and commonly used in $r$-process calculations. All reaction and decay rates are consistently calculated with model masses as in Ref. \cite{Mumpower+15} and all fission rates are updated based on the fission barrier heights of a given model. The fission barriers are illustrated in Fig.~\ref{fig:bhrp}. The barriers for the Finite Range Liquid Droplet Model (FRLDM) are from \cite{MollerBH1,MollerBH2}, Thomas-Fermi barriers are again considered to be those applied in GEF 2016, and ETFSI barriers taken from \cite{Mamdouh98,Mamdouh01}. The HFB-14 fission paths used to determine the barriers are available for $Z\geq 90$ in the BRUSLIB database \cite{BRUSLIB} as well as TALYS. When the data truncates at a mass or charge number, as for HFB barriers at $Z<90$ and ETFSI masses at $A>300$, we apply rates based on FRDM2012 masses and FRLDM barriers for the absent nuclei. For the fission yields in this section, we apply the default GEF 2016 distributions presented in Sec. II and do not update these to reflect the barriers of each model. With a fixed fission yield model, we can study how the abundance pattern is shaped by these yields when different sets of fissioning nuclei, as determined by the fission barriers, are accessed.

\begin{figure*}
\begin{center}
\includegraphics[scale=0.53]{abA_RossAvg_AllModC_nc.pdf}
\end{center}
\caption{(Color online) Range (grey band) and average (colored line) for the final $r$-process abundances with TF (upper left), FRDM (upper right), ETFSI (lower left) and HFB (lower right) models given thirty 1.2--1.4 M$_{\odot}$ neutron star merger simulation trajectories \cite{Rosswog}.}
\label{fig:rpRossAll}
\end{figure*}

We first consider ``cold" dynamical ejecta conditions of trajectory 1 \cite{Rosswog} in order to examine the case in which the nuclear flow can reach the highest possible mass numbers because it is less impeded by photodissociation. The black outline in Fig. \ref{fig:bhrp} shows the location of the most populated nuclei (abundance $\ge 10^{-10}$) just before the $r$-process path begins to move back toward stability (taken to be the time when the abundance weighted mass number reaches its last maximum). With FRDM and TF, material encounters relatively low barrier heights $\sim 4-5$ MeV as it pushes past the $N=184$ shell closure, permitting fission to occur quickly and thus preventing the nuclear flow from continuing much higher in mass number. As noted in Ref. \cite{Eichler15}, we find HFB permits the synthesis of heavier nuclei than FRDM and TF models due to higher barriers near the predicted $N=184$ shell closure. The same is true for the ETFSI model \cite{Petermann,Eichler16}, which was previously found to have its nuclear flow terminated by $\beta$-delayed fission \cite{Petermann} instead of neutron-induced fission which terminates the path using TF and FRLDM barriers. We also find $\beta$-delayed fission to be more active in the ETFSI case, with a $\sim60\%$ enhancement in the total integrated $\beta$-delayed fission flow relative to the flow found with TF or FRDM. We note that the ability to synthesize nuclei with higher mass number makes $r$-process calculations using ETFSI and HFB models more sensitive than results with FRDM and TF to the increase in prompt neutrons and widening of yields predicted by FREYA which become more significant at higher neutron and proton numbers (as shown in Figs.~\ref{fig:FREYAnu} and \ref{fig:FREYAy}) (to show this explicitly the version of Fig. \ref{fig:abFREYA} given the astrophysical conditions of traj.\ 22 with HFB model inputs has been included in Supplemental Materials).

We next turn to recalculate the thirty trajectories of the 1.2--1.4 M$_{\odot}$ neutron star merger simulation \cite{Rosswog} considered in Sec.\ III with the four sets of mass and barrier models described above. The results appear in Fig.~\ref{fig:rpRossAll}. Although these dynamical ejecta trajectories are all very neutron-rich, the variations in their density and temperature profiles cause a spread in the range of predicted abundances along the right edge of the second peak. As previously discussed in the context of Fig.~\ref{fig:rtraj117}, this is largely due to whether material reaches the region past $N=184$ where GEF 2016 yields are mostly symmetric and deposit material near $A\sim144$.
The ability of ETFSI and HFB-17 to reach nuclei higher in mass number due to higher fission barriers near $N=184$ means these models access more of this region found by GEF 2016 to have mostly symmetric yields as compared to the nuclei accessed by FRDM and TF.  We find that the higher nuclear flow through the symmetric GEF yield region with ETFSI and HFB, seen in Fig.~\ref{fig:bhrp}, contributes to their overproduction of the right edge of the second peak as compared to solar data. The model dependent shell closure predictions also play a role with FRDM having a stronger $N=82$ shell closure than ETFSI and HFB which keeps fission daughter products closer to the $A\sim130$ region at late times. It is therefore the interplay between the barrier height landscape around $N=184$ and the structure of the $N=82$ shell closure which determines the shape of the second $r$-process peak in fission cycling conditions.
It was previously found in Ref. \cite{Eichler15} that given HFB-14 barriers, most of the fissioning nuclei accessed had fission fragments lying between $A = 125$ and $A = 155$, similar to their results with the FRDM case.
We find that the HFB-14 barriers produce fission flow which dominantly accesses fragment distributions centered at $132\leq A\leq 150$, even after material is primarily located at $N<184$. This contributes to the overall shift in the second $r$-process peak with HFB in Fig.~\ref{fig:rpRossAll}  (to show this explicitly the version of Fig. \ref{fig:rtraj1} given the astrophysical conditions of traj.\ 1 and HFB model inputs has been included in Supplemental Materials).

\begin{figure}
\begin{center}
\includegraphics[scale=0.527]{abA_Sneden_Ross1FRDM12GEFMollvMarkC_ncd.pdf}
\end{center}
\caption{(Color online) The $r$-process abundances at 1 Gyr using the cold tidal tail conditions of trajectory 1 \cite{Rosswog} and the FRDM2012 mass model with $\beta$-decay rates determined from M$\ddot{\mathrm{o}}$ller {\it et al.} QRPA calculations \cite{MollerQRPA} (orange) as compared to those from Marketin {\it et al.} \cite{Marketin} (grey).}
\label{fig:rpMark}
\end{figure}

Having considered the role of masses and fission barriers, we next turn to how $\beta$-decay rates influence which fissioning nuclei are accessed in an $r$-process calculation. We repeat the FRLDM/FRDM simulations described above with the M$\ddot{\mathrm{o}}$ller {\it et al.} \cite{MollerQRPA} $\beta$-decay rates replaced by those from Marketin {\it et al.} \cite{Marketin} and show an example abundance pattern comparison in Fig.~\ref{fig:rpMark}. The two sets of rates are generally similar except for nuclei above the $N=126$ shell closure, where the Marketin {\it et al.} rates are faster. Simulations with the faster Marketin {\it et al.} rates tend to show less material hung up in the higher mass regions than with M$\ddot{\mathrm{o}}$ller {\it et al.} rates which reduces the extra post-freeze-out neutrons produced via fission and $\beta$-delayed neutron emission. This in turn thwarts the shifting and narrowing of the third peak from late-time neutron capture (as noted in Refs. \cite{Eichler15,Eichler16}). This reduction of nuclei present near $N=184$ hinders the opportunity for fission to build and shape the second $r$-process peak at late times. Specifically, with nuclei near $N=184$ less populated at freeze-out, the region of the nuclear chart where GEF 2016 predicts fission yields to be centered mostly near $A\sim130$ is not accessed very heavily. This is responsible for the differences in the height of the second $r$-process peak when results using the M$\ddot{\mathrm{o}}$ller {\it et al.} \cite{MollerQRPA} and Marketin {\it et al.} rates are compared (to show this explicitly the versions of Fig. \ref{fig:rtraj1} given the Marketin {\it et al.} and M$\ddot{\mathrm{o}}$ller {\it et al.} conditions applied in Fig. \ref{fig:rpMark} have been included in Supplemental Materials).

 \begin{figure}
\begin{center}
\includegraphics[scale=0.49]{HeatingGEFKod5050_Ross1FRDM12C_1pan.pdf}
\end{center}
\caption{(Color online) The nuclear heating rates for $\beta$-decay (blue), $\beta$-delayed fission (light blue), neutron capture (red) and neutron-induced fission (orange) as a function of time using GEF 2016 fission yields (solid lines) as compared to results using simple symmetric splits (dashed lines) and the yields of Kodama and Takahashi (dotted lines).  The masses and fission barriers applied are those of FRDM2012 and FRLDM respectively.}
\label{fig:heatyield2}
\end{figure}

 \begin{figure}
\begin{center}
\includegraphics[scale=0.49]{LATEheatingGEFFREYA5050_Ross1FRDM12_1pan.pdf}
\end{center}
\caption{(Color online) The nuclear heating rate for $\beta$-decay (blue), alpha decay (black) and spontaneous fission (green) as a function of time using GEF 2016 fission yields (solid lines) as compared to results using simple symmetric splits (dashed lines) and the yields with the FREYA energy sharing and de-excitation treatment (dotted lines). The masses and fission barriers applied are those of FRDM2012 and FRLDM respectively.}
\label{fig:spfheat2}
\end{figure}

We lastly consider how the masses and fission barriers influence the sensitivity of the nuclear heating rates to fission yield treatments. To do so we repeat the calculations presented in Sec. III with nuclear rates determined from FRDM2012 masses and FRLDM fission barriers, shown in Figs.~\ref{fig:heatyield2} and  \ref{fig:spfheat2}. Figure~\ref{fig:heatyield2} confirms that the three order of magnitude discrepancy in the heating for neutron capture channels when 50/50 and GEF 2016 yield results are compared is not isolated to the TF model. This increase in neutron-induced fission in the GEF case is even stronger given FRDM inputs (compare Figs.~\ref{fig:heatyield} and  \ref{fig:heatyield2}) since here the rise in heating rate of neutron capture channels reaches a value $\sim$ one order of magnitude lower than the dominant $\beta$-decay heating channel.

When late time heating is considered, a comparison of Figs.~\ref{fig:spfheat} and  \ref{fig:spfheat2} shows that the spontaneous fission contribution to the late-time nuclear heating is two orders of magnitude larger with the FRDM masses and FRLDM barriers as compared to the result using TF masses and GEF 2016 barriers. In Fig.~\ref{fig:Cfbar} we show that it is the difference between the barriers of these two models along the $A=254$ isobaric chain that is primarily responsible for the difference in the late-time contributions from the long-lived californium nucleus $^{254}$Cf. We find that in this region of the nuclear chart, fission flows are highest in the presence of barrier heights around $4-5$ MeV. In the TF model, this is precisely the height of the barrier which $A=254$ nuclei must overcome in order to eventually populate californium, and we find fission of these nuclei transfers material out of this isobaric chain. In contrast, with the FRLDM model, the nuclei set to populate $^{254}$Cf pass just to the left of the region with barrier heights of $4-5$ MeV and the fission encountered along the path to $^{254}$Cf is insignificant. Note from the discussion in Sec. II that the GEF 2016 barriers include some systematic corrections in this region of the nuclear chart. However, along the $A=254$ isobaric line, the TF barriers without GEF corrections are also $4-5$ MeV or lower. Specifically we find that the lower barriers of TF/GEF roughly translate into an order of magnitude higher neutron-induced fission rate for two key nuclei ($^{254}$Np and $^{255}$U) along the path feeding $^{254}$Cf. Though these two key nuclei are populated late in the $r$ process (on the order of $\sim 1-10$ seconds) when the neutron abundance has dropped significantly, their neutron-induced fission rates can be sufficiently high to interrupt the $\beta$ feeding of $^{254}$Cf. We do not find the differences in $\beta$-delayed fission with TF and FRDM models to significantly affect the population of $^{254}$Cf.

\begin{figure}
\begin{center}
\includegraphics[scale=0.56]{254CfpathBarrierC_nc.pdf}
\end{center}
\caption{(Color online) Fission barrier heights for FRLDM (upper panel) and GEF 2016 (TF) (lower panel) along with the $\beta$-decay path for $^{254}$Cf in each case (black dashed line). The blue outline shows the region of experimentally established decay rates.}
\label{fig:Cfbar}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[scale=0.53]{Cf254_Ross30_AllModYeC_1pan_band_nc.pdf}
\end{center}
\caption{(Color online) The $^{254}$Cf abundance as a function of time with TF (green), FRDM (orange), ETFSI (blue) and HFB (purple) models. The bands represent the range in $^{254}$Cf production given thirty 1.2--1.4 M$_{\odot}$ neutron star merger simulation trajectories \cite{Rosswog}.}
\label{fig:254Cfallmod}
\end{figure}

Finally, we comment on the population of $^{254}$Cf given all the nuclear mass/fission barrier models considered in this section. The $^{254}$Cf abundance as a function of time for all cases is shown in Fig. \ref{fig:254Cfallmod}. The two order of magnitude reduction in the abundance of this nucleus with the TF model relative to FRDM is expected given the differences in barrier heights and fission rates previously discussed. The ETFSI model also populates $^{254}$Cf less than FRDM. In this case, the barrier height landscape near $N=184$ leaves less material available to populate  $^{254}$Cf when it pushes significantly past this predicted shell closure. After fission terminates the flow of material near $N=210$, the significant amount of material piled up here encounters the right edge of a large region with low barriers and fissions quickly. In addition, the ETFSI model predicts the $^{254}$Cf path to encounter a small region of barriers with heights $4-5$ MeV which reduce the population $\beta$ feeding this nucleus. Much like ETFSI, the HFB model barriers allow a significant flow of material past $N=184$ which will then not be available to populate the $A=254$ isobaric chain. However the material that remains piled up near $N=184$ at freeze-out which is capable of populating $^{254}$Cf is not inhibited by low $4-5$ MeV barrier heights during its decay back to stability. This produces a higher predicted $^{254}$Cf abundance with HFB as compared to ETFSI and TF models. We find that it is the FRDM model which populates $^{254}$Cf most strongly through the coupled effects of low barriers just past $N=184$ preventing material from moving to higher mass number as well as sufficiently high barriers along the $^{254}$Cf path that prevent a significant depopulation from fission. It is clear that the fission barrier heights of heavy unstable nuclei are key to assessing the influence of $^{254}$Cf on kilonova light curves. Therefore, potential observations of the increased light curve luminosity associated with the heating ``bump" from $^{254}$Cf at late times would not only be able to confirm the synthesis of long-lived actinides, but actually inform nuclear physics calculations of fission properties in this heavy, unstable region.

\section{Fission Hot Spots}\label{sec:maxflows}

With fission in the $r$ process occurring over a large range of mass numbers, and with much of the initial fission concentrated as far out as the neutron drip line, it is not obvious how practical it is for experimental and theoretical efforts to refine our knowledge of fissioning $r$-process nuclei. To push experiments even a few neutron numbers out from presently-studied nuclei could require tremendous efforts (for a brief review of experimental campaigns to study the neutron-rich isotopes of heavy $r$-process nuclei see \cite{HorowitzRIB2018} and references therein). Theoretical campaigns to calculate fission yields starting from the nuclear properties assumed for a given model, such as density functional theory approaches, can assist with predictions of fission yields for nuclei far from experimental reach. However such methods can be computationally expensive. Therefore given the impracticality for experimental, and even theoretical, studies to provide information on the fission yields of all nuclei of interest to the $r$ process, a guide as to which nuclei participate most during fission in the $r$ process is needed. Here we provide such information by finding the ``hot spots" of nuclei with the highest fission flow in our $r$-process calculations. To do so, we average over $r$-process results given thirty trajectories from a 1.2--1.4 M$_{\odot}$ neutron star merger simulation \cite{Rosswog} (as in Fig.~\ref{fig:rpRossAll}). For each mass model considered, we report the average integrated fission flow for the neutron-induced and $\beta$-delayed fission processes which we find determine the final $r$-process abundances.

We confine our discussion to neutron-induced and $\beta$-delayed fission since here spontaneous fission flows are comparatively much lower. In this work we apply spontaneous fission rates determined by the phenomenological equation of Karpov {\it et al.} \cite{Karpov} which depends on fissility ($Z^2/A$) and barrier height. These spontaneous fission rates are very low until $Z>100$. Therefore, for all of the simulations considered here, $r$-process material tends to only encounter spontaneous fission at late times when the main abundance pattern is nearly finalized. This is consistent with previous studies which surveyed the influence of several descriptions of spontaneous fission rates \cite{PanovSF} and observed this process to weakly influence the abundances of the second to third $r$-process peaks. However, given the sparsity of experimental fission data, there exist phenomenological fits to spontaneous fission half-lives which, when extrapolated into neutron-rich regions, give very high spontaneous fission rates starting at $Z>94$ \cite{XuRen}. Such rates effectively cut-off the influence of neutron-induced and $\beta$-delayed fission at higher proton numbers. Therefore, the ``hotspots" reported here are meant to represent a case where neutron-induced and $\beta$-delayed fission alone shape the final $r$-process abundance pattern.

\begin{figure*}
\begin{center}
\includegraphics[scale=0.54]{AvgNIFIntFlow_Rosswog30_allMMC_nc.pdf}
\end{center}
\caption{(Color online) The integrated neutron-induced fission flow averaged over thirty astrophysical trajectories from a 1.2--1.4 M$_{\odot}$ neutron star merger simulation \cite{Rosswog} assuming GEF (TF) (top left), FRLDM (top right), ETFSI (bottom left), and HFB-14 (bottom right) barrier heights.}
\label{fig:hsnif}
\end{figure*}

\begin{figure*}
\begin{center}
\includegraphics[scale=0.54]{AvgBDFIntFlow_Rosswog30_allMMC_nc.pdf}
\end{center}
\caption{(Color online) The integrated $\beta$-delayed fission flow averaged over thirty astrophysical trajectories from a 1.2--1.4 M$_{\odot}$ neutron star merger simulation \cite{Rosswog} assuming GEF (TF) (top left), FRLDM (top right), ETFSI (bottom left), and HFB-14 (bottom right) barrier heights. The black outline here shows where the probability for nuclei to undergo multi-chance $\beta$-delayed fission (relative to ordinary $\beta$-decay) exceeds $10\%$.}
\label{fig:hsbdf}
\end{figure*}

The averaged integrated neutron-induced and $\beta$-delayed fission flows are shown in Figs.~\ref{fig:hsnif} and \ref{fig:hsbdf}. Each panel shows high fission flows are obtained in two general regions of the nuclear chart: along isotopic chains $90<Z<95$ where fission terminates the $r$-process path in $A$, and along the $\beta$-decay pathways of nuclei initially hung up near the $N=184$ shell closure. The exact nuclei in each general region which will have high fission flow is strongly dependent on nuclear masses and barrier heights.  The Thomas-Fermi model used in GEF 2016 as well as the FRLDM model predict a large range of nuclei with relatively low barrier heights $\sim4-6$ MeV near the $N=184$ shell closure. These models therefore show significant fission flow to the left of $N=184$ for the neutron-rich actinide region of $89<Z<95$ which is where GEF 2016 yields show asymmetric yield contributions. Less of this asymmetric region is accessed by HFB-17 and ETFSI models which have higher barrier heights $\sim6-8$ MeV near $N=184$ which hinders fission of material until it moves higher in mass number during the decay back to stability.

We next identify the nuclei whose fission yields are primarily responsible for the shape of the second $r$-process peak. To do so we consider the threshold of integrated fission fission flow that contributes to finalizing the $r$-process abundance pattern. For each mass model, we apply GEF 2016 neutron-induced fission yields for the nuclei with an average integrated neutron-induced fission above a set threshold with symmetric splits assumed for all remaining nuclei. We apply the same threshold criterion for $\beta$-delayed fission. We find that applying GEF yields to only nuclei with an average integrated fission flow larger than $10^{-5}$ reproduces the final abundance trend almost exactly for all mass models. In the case of TF, HFB, and ETFSI models, implementing the GEF yields for only nuclei having an average integrated fission flow larger than $10^{-4}$ was found to be sufficient to reproduce the relative abundances.

We note Figs.~\ref{fig:hsnif} and \ref{fig:hsbdf} show many of the nuclei with fission flows higher than the $10^{-5}$ threshold have odd neutron number. Are these flows primarily responsible for shaping the final abundance pattern, or do even-$N$ nuclei play a larger role? It is not immediately obvious that high odd-$N$ flow implies a great influence on the abundance pattern since, roughly speaking, the even-$N$ abundances are greater than those of odd-$N$ nuclei throughout the $r$ process. However fission, particularly in the case of neutron-induced, can have rates for an odd-$N$ initial species which are $\sim8-10$ orders of magnitude larger than their even-$N$ isotopic neighbor. We find this disparity in the rates to dominate over the abundance preference for even-$N$ nuclei. Previous work noted the expectation for neutron-induced fission to mainly occur with an even-$N$ fissioning species due to lower values in the difference between the fission barrier and neutron separation energy (since fission barriers do not suffer from strong odd-even effects as is the case for the neutron separation energies) \cite{Petermann}. Here our conclusions are not drawn solely from the odd-even dependence of the fission rates, but also account for the population of fissioning nuclei in the $r$ process by considering flow. We also find that the dominance of odd-$N$ flow, which suggests the importance of the fission yields of even-$N$ species, generally extends to $\beta$-delayed fission as well. In fact, we find fission reactions which have an odd-$N$ initial species are so influential that the $r$-process abundance pattern can nearly be reproduced when GEF 2016 yields are applied to the daughters of just odd-$N$ nuclei, as shown in Fig.~\ref{fig:abthres}. Therefore it is the fission rates of odd-$N$ target nuclei and the fission yields of the corresponding even-$N$ compound nuclei which most impact our $r$-process calculations.

To report the exact nuclei whose daughter yields are of most consequence when assuming a particular mass model, in Supplemental Materials we tabulate the nuclei that satisfy the $10^{-5}$ threshold criterion for each of the mass models considered. We find that all mass models predict the fission outcomes of $200$ nuclei or less to be of relevance in setting the $r$-process abundance pattern. If only the odd-$N$ nuclei are considered, then all mass models predict $120$ nuclei or less to be relevant.  Although the nuclei which most impact the shape of the abundance pattern are dependent on the mass model and fission barriers, in Fig.~\ref{fig:hotspots} we highlight nuclei commonly found to have a high average integrated fission flow ($\ge10^{-5}$) given the range of nuclear inputs considered. The earlier onset of neutron-induced fission, as compared to $\beta$-delayed fission, is reflected in the reach of this fission channel beyond $N=184$ where abundances are high only at early times. Note that Fig.~\ref{fig:hotspots} essentially reports the overlap of nuclei found to have high fission flow in Figs.~\ref{fig:hsnif} and \ref{fig:hsbdf}, which was mostly determined by a given model's fission barriers. Therefore, the shape of these hotspots is highly influenced by where models agree nuclei have low $4-6$ MeV barrier heights.

Our calculations with the nuclear inputs applied here see all four models agree upon 15 nuclei to be of importance for the neutron-induced fission channel, with over half of these found in the $Z=93$ and $94$ isotopic chains (all of which have odd-$N$). The $\beta$-delayed fission channel sees all models agree upon only 7 nuclei of importance, with nearly all of these nuclei in $Z=93$ and $Z=97$ isotopic chains and most (but not all) having odd-$N$. We note that since the fission yields do affect the flow of material in a fission cycling $r$-process, the ``hotspots" found using 50/50 splits differ for a small handful of nuclei, but the overall region of importance remains the same. In previous work, the mass region with $93 \leq Z \leq 95$ and $180 \leq N \leq 186$ was identified as the dominant region for neutron-induced and $\beta$-delayed fission flow for the FRDM case considered in Ref.~\cite{Eichler15}. In Refs.~\cite{GorielyGEF,GorielyGMPGEF}, nuclei with $A\simeq278$ were singled out as the isobars whose fission products determine the abundance of nuclei in the $110\lesssim A \lesssim 170$ region. Although we also see high fission flow in this region near the $N=184$ shell closure, we see much of the high flow concentrated at $N<184$ since fission is still very active during the decay back to stability.

\begin{figure}
\begin{center}
\includegraphics[scale=0.53]{abA_FRDM12_Ross1EvenOddThres105KodC.pdf}
\end{center}
\caption{(Color online) The $r$-process abundance pattern using the cold tidal tail conditions of trajectory 1 \cite{Rosswog} and the FRDM2012 mass model when GEF 2016 yields are applied to all nuclei (dashed red) as compared to GEF yields for only the odd-$N$ fission reactions of ($Z$,$N$)(n,f) and ($Z$,$N$)$\beta$df (solid purple) with all other fissioning nuclei assuming the fission yields of \cite{Kodama} (K$\&$T). For comparison the abundances with GEF yields applied to only the even-$N$ fission reactions of ($Z$,$N$)(n,f) and ($Z$,$N$)$\beta$df (solid blue) as well as K$\&$T applied to all nuclei (dot-dashed dark blue) are also shown.}
\label{fig:abthres}
\end{figure}

\begin{figure*}
\begin{center}
\includegraphics[scale=0.54]{Hotspots2pan_Rosswog30_CountMMThres105C_nc.pdf}
\end{center}
\caption{(Color online) The nuclei found to have an average integrated neutron-induced (upper) /$\beta$-delayed (lower) fission flow above the $10^{-5}$ threshold found to control the $r$-process abundance pattern in one (yellow/green), two (orange/light blue), three (red/dark blue), or four (dark red/purple) of the mass models considered in Figure \ref{fig:hsnif}/\ref{fig:hsbdf}. The black boxes show stable nuclei while the black outline shows the location of known experimental decay rates.}
\label{fig:hotspots}
\end{figure*}

\section{Conclusions}

The influence of the astrophysical conditions on whether a fission cycling $r$ process is achieved, as well as the dependence on unknown nuclear physics properties for the heaviest neutron-rich nuclei, imply astrophysical environments which host fission to involve the greatest number of $r$-process uncertainties. If ``red'" kilonova components originate from the very neutron-rich conditions seen in neutron star merger dynamical ejecta, fission most certainly plays a role in setting the lanthanide mass fraction which determines wavelength and peak luminosity. If the ``red" kilonova is instead generated primarily from an accretion disk wind, astrophysical uncertainties in the exact neutron richness of such conditions make it difficult to know how much the treatment of fission can influence observation. Since neutron star mergers permit many possible routes to a fission cycling $r$ process, it is important to understand the potential impact on $r$-process observables from variations in the fission treatment.

We have shown that taking into account the fission yield dependence on the initial excitation energy can influence the final $r$-process abundances by using the publicly available 2016 GEF code to obtain sets of fission yields for neutron-induced, $\beta$-delayed, and spontaneous fission with an appropriate energy applied for each respective fission process. We considered the sensitivity of our results to the treatment of the excitation-energy sharing and de-excitation of the fission fragments using FREYA and showed that such considerations can change the average prompt neutron multiplicity by $\sim1-3$ neutrons in the most neutron-rich regions. However, the sensitivity of the $r$ process to the energy sharing and de-excitation treatment was found to be secondary compared to the effect of the fission yield dependence on initial excitation energy, demonstrating the treatment of the primary fragment yields to be of dominant influence in the $r$ process. Thus careful theoretical calculations of fission fragment yields (prior to neutron emission) that are consistently derived within the framework of a given mass model are crucial to understanding lanthanide production in a fission cycling $r$ process.

We showed that the trend in GEF 2016 yields, which transition from asymmetric to mostly symmetric yields along an isotopic chain, can reproduce the trend of the right edge of the second $r$-process peak seen in solar data given fission cycling conditions that reach the most neutron-rich regions beyond $N=184$. We considered the influence of nuclear mass models by applying fission rates that self-consistently account for the dependence on nuclear masses and fission barriers. We found that the fission flow explicitly follows the regions of low $4-6$ MeV fission barrier heights making the fissioning nuclei which most impact $r$-process calculations model dependent. The population of key fissioning actinides, such as $^{254}$Cf, was also shown to sensitively depend on the fission barrier assumptions of a given model.

For each of the four sets of masses and barrier heights considered, we reported the integrated fission flow averaged over 30 simulation trajectories for a 1.2--1.4 M$_{\odot}$ neutron star merger to find the fissioning nuclei of most importance. The odd-even behavior observed in these fission flows lead us to identify odd-$N$ species undergoing ($Z$,$N$)(n,f) and ($Z$,$N$)$\beta$df to substantially dominate over even-$N$ species in setting the $r$-process abundance pattern. The ``hot spots" of key fissioning nuclei given all models considered show that nuclei of importance are often found to have $N<184$ due to pile-up of material at this predicted shell closure. The proximity of these ``hot spots" to the region of experimentally established decay data shows the potential for experiments to access some of the fissioning nuclei found to play a key role in setting the $r$-process abundance pattern. Efforts by experimental and theoretical nuclear physics to further the knowledge of fission properties in neutron-rich regions are necessary to develop a more complete picture of heavy element production in neutron star mergers.

\begin{acknowledgments}
N.V. would like to thank Nicolas Schunck, Karl-Heinz Schmidt, Stephan Goriely, Sean Liddick, Samuel Giuliani and Marius Eichler for useful discussions. The work of N.V., R.S., M.M., P.J. and R.V. was partly supported by the Fission In R-process Elements (FIRE) topical collaboration in nuclear theory, funded by the U.S. Department of Energy.  The work of R.V. was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344.  The work of J.R. was performed under the auspices of the U.S. Department of Energy by Lawrence Berkeley National Laboratory under Contract DE-AC02-05CH11231.  The work of M.M. and P.J. was supported in part under the auspices of the National Nuclear Security Administration of the U.S. Department of Energy at Los Alamos National Laboratory under Contract No. DE-AC52-06NA25396.  The work of R.S., T.M.S., and E.M.H. was also supported in part by the U.S. Department of Energy under grant numbers DE-SC0013039 and was enabled by the National Science Foundation under Grant No. PHY-1430152 (JINA Center for the Evolution of the Elements). This manuscript is locatable via Los Alamos and Lawerence Livermore National Laboratory report numbers LA-UR-19-22109 and LLNL-JRNL-760148, respectively.
\end{acknowledgments}

$\,\,$\newline
$\,\,$\newline
$\,\,$\newline
$\,\,$\newline
$\,\,$\newline