We discuss the two-channel Kondo problem with a pseudogap density of states $\ensuremath{\rho}(\ensuremath{\omega})\ensuremath{\propto}{|\ensuremath{\omega}|}^{r}$ of the bath fermions. Combining both analytical and numerical renormalization group techniques, we characterize the impurity phases and quantum phase transitions of the relevant Kondo and Anderson models. The line of stable points, corresponding to the overscreened non-Fermi-liquid behavior of the metallic
r=0 case, is replaced by a stable particle-hole-symmetric intermediate-coupling fixed point for $0<r<{r}_{\mathrm{max}}\ensuremath{\approx}0.23$. For
r>rmax, this non-Fermi-liquid phase disappears, and instead a critical fixed point with an emergent spin-channel symmetry appears, controlling the quantum phase transition between two phases with stable spin and channel moments, respectively. We propose low-energy field theories to describe the quantum phase transitions, all being formulated in fermionic variables. We employ $\ensuremath{\epsilon}$ expansion techniques to calculate critical properties near the critical dimensions
r=0 and 1, the latter being potentially relevant for two-channel Kondo impurities in neutral graphene. We find the analytical results to be in excellent agreement with those obtained from applying Wilson's numerical renormalization group technique.
We address the phenomenon of critical Kondo destruction in pseudogap Bose-Fermi Anderson and Kondo quantum impurity models. These models describe a localized level coupled both to a fermionic bath having a density of states that vanishes like ${|\ensuremath{\epsilon}|}^{r}$ at the Fermi energy ($\ensuremath{\epsilon}=0$) and, via one component of the impurity spin, to a bosonic bath having a sub-Ohmic spectral density proportional to ${|\ensuremath{\omega}|}^{s}$. Each bath is capable by itself of suppressing the Kondo effect at a continuous quantum phase transition. We study the interplay between these two mechanisms for Kondo destruction using continuous-time quantum Monte Carlo for the pseudogap Bose-Fermi Anderson model with
0<r<21 and $\frac{1}{2}\ensuremath{\le}s<1$, and applying the numerical renormalization group to the corresponding Kondo model. At particle-hole symmetry, the models exhibit a quantum-critical point between a Kondo (fermionic strong-coupling) phase and a localized (Kondo-destroyed) phase. The two solution methods, which are in good agreement in their domain of overlap, provide access to the many-body spectrum, as well as to correlation functions including, in particular, the single-particle Green's function and the static and dynamical local spin susceptibilities. The quantum-critical regime exhibits the hyperscaling of critical exponents and $\ensuremath{\omega}/T$ scaling in the dynamics that characterize an interacting critical point. The
(r,s) plane can be divided into three regions: one each in which the calculated critical properties are dominated by the bosonic bath alone or by the fermionic bath alone, and between these two regions, a third in which the bosonic bath governs the critical spin response but both baths influence the renormalization-group flow near the quantum-critical point.