We present nonperturbative fragmentation functions (FFs) for bottom-flavored (B) hadrons both at next-to-leading (NLO) and, for the first time, at next-to-next-to-leading order in the ¯¯¯¯¯¯MS factorization scheme with five massless quark flavors. They are determined by fitting all available experimental data of inclusive single B-hadron production in e+e− annihilation, from the ALEPH, DELPHI, and OPAL collaborations at CERN LEP1 and the SLD collaboration at SLAC SLC. The uncertainties in these FFs as well as in the corresponding observables are estimated using the Hessian approach. We perform comparisons with available NLO sets of B-hadron FFs. We apply our new FFs to generate theoretical predictions for the energy distribution of B hadrons produced through the decay of unpolarized or polarized top quarks, to be measured at the CERN LHC.