We determine the normalization of scalar and pseudoscalar current operators made from nonrelativistic b quarks and highly improved staggered light quarks in lattice quantum chromodynamics (QCD) through Oðα s Þ and Λ QCD =m b . We use matrix elements of these operators to extract B meson decay constants and form factors, and then compare to those obtained using the standard vector and axial-vector operators. This provides a test of systematic errors in the lattice QCD determination of the B meson decay constants and form factors. We provide a new value for the B and B s meson decay constants from lattice QCD calculations on ensembles that include u, d, s, and c quarks in the sea and those that have the u=d quark mass going down to its physical value. Our results are f B ¼ 0.196ð6Þ GeV, f B s ¼ 0.236ð7Þ GeV, and f B s =f B ¼ 1.207ð7Þ, agreeing well with earlier results using the temporal axial current. By combining with these previous results, we provide updated values of f B ¼ 0.190ð4Þ GeV, f B s ¼ 0.229ð5Þ GeV, and f B s =f B ¼ 1.206ð5Þ.