$B^0_s - \bar{B^0_s}$ Mixing Parameters from Unquenched Lattice QCD

We determine hadronic matrix elements relevant for the mass and width differences, $\Delta M_s$ & $\Delta \Gamma_s$ in the $B^0_s - \bar{B^0_s}$ meson system using fully unquenched lattice QCD. We employ the MILC collaboration gauge configurations that include $u$, $d$ and $s$ sea quarks using the improved staggered quark (AsqTad) action and a highly improved gluon action. We implement the valence $s$ quark also with the AsqTad action and use Nonrelativistic QCD for the valence $b$ quark. For the nonperturbative QCD input into the Standard Model expression for $\Delta M_s$ we find $f_{B_s} \sqrt{\hat{B}_{B_s}} = 0.281(21)$GeV. Results for four-fermion operator matrix elements entering Standard Model formulas for $\Delta \Gamma_s$ are also presented.