  The chart of the nuclides is limited by particle drip lines beyond which
nuclear stability to proton or neutron emission is lost. Predicting the range
of particle-bound isotopes poses an appreciable challenge for nuclear theory as
it involves extreme extrapolations of nuclear masses beyond the regions where
experimental information is available. Still, quantified extrapolations are
crucial for a variety of applications, including the modeling of stellar
nucleosynthesis. We use microscopic nuclear mass models and Bayesian
methodology to provide quantified predictions of proton and neutron separation
energies as well as Bayesian probabilities of existence throughout the nuclear
landscape all the way to the particle drip lines. We apply nuclear density
functional theory with several energy density functionals. To account for
uncertainties, Bayesian Gaussian processes are trained on the separation-energy
residuals for each individual model, and the resulting predictions are combined
via Bayesian model averaging. This framework allows to account for systematic
and statistical uncertainties and propagate them to extrapolative predictions.
We characterize the drip-line regions where the probability that the nucleus is
particle-bound decreases from $1$ to $0$. In these regions, we provide
quantified predictions for one- and two-nucleon separation energies. According
to our Bayesian model averaging analysis, 7759 nuclei with $Z\leq 119$ have a
probability of existence $\geq 0.5$. The extrapolations obtained in this study
will be put through stringent tests when new experimental information on exotic
nuclei becomes available. In this respect, the quantified landscape of nuclear
existence obtained in this study should be viewed as a dynamical prediction
that will be fine-tuned when new experimental information and improved global
mass models become available.
