Strong Bounds on Sum of Neutrino Masses in a 12 Parameter Extended
  Scenario with Non-Phantom Dynamical Dark Energy ($w(z)\geq -1$)
  We obtained constraints on a 12 parameter extended cosmological scenario
including non-phantom dynamical dark energy (NPDDE) with CPL parametrization.
We also include the six $\Lambda$CDM parameters, number of relativistic
neutrino species ($N_{\textrm{eff}}$) and sum over active neutrino masses
($\sum m_{\nu}$), tensor-to-scalar ratio ($r_{0.05}$), and running of the
spectral index ($n_{run}$). We use CMB Data from Planck 2015; BAO Measurements
from SDSS BOSS DR12, MGS, and 6dFS; SNe Ia Luminosity Distance measurements
from the Pantheon Sample; CMB B-mode polarization data from BICEP2/Keck
collaboration (BK14); Planck lensing data; and a prior on Hubble constant
($73.24\pm1.74$ km/sec/Mpc) from local measurements (HST). We have found strong
bounds on the sum of the active neutrino masses. For instance, a strong bound
of $\sum m_{\nu} <$ 0.123 eV (95\% C.L.) comes from Planck+BK14+BAO. Although
we are in such an extended parameter space, this bound is stronger than a bound
of $\sum m_{\nu} <$ 0.158 eV (95\% C.L.) obtained in $\Lambda \textrm{CDM}+\sum
m_{\nu}$ with Planck+BAO. Varying $A_{\textrm{lens}}$ instead of $r_{0.05}$
however leads to weaker bounds on $\sum m_{\nu}$. Inclusion of the HST leads to
the standard value of $N_{\textrm{eff}} = 3.045$ being discarded at more than
68\% C.L., which increases to 95\% C.L. when we vary $A_{\textrm{lens}}$
instead of $r_{0.05}$, implying a small preference for dark radiation, driven
by the $H_0$ tension.
