How does the velocity anisotropy of halo stars, dark matter and satellite galaxies depend on host halo properties?

We investigate the mass ($M_{200}$) and concentration ($c_{200}$) dependencies of the velocity anisotropy ($β$) profiles for different components in the dark matter halo, including halo stars, dark matter and subhalos, using systems from the IllustrisTNG simulations. Beyond a critical radius, $β$ becomes more radial with the increase of $M_{200}$, reflecting more prominent radial accretion around massive halos. The critical radius is $r\sim r_s$, $0.3~r_s$ and $r_s$ for halo stars, dark matter and subhalos, with $r_s$ the scale radius of host halos. This dependence on $M_{200}$ is the strongest for subhalos, and the weakest for halo stars. In central regions, $β$ of halo stars and dark matter particles gets more isotropic with the increase of $M_{200}$ in TNG300 due to baryons. By contrast, $β$ of dark matter from the dark matter only TNG300-Dark run shows much weaker dependence on $M_{200}$ within $r_s$. Dark matter in TNG300 is slightly more isotropic than in TNG300-Dark at $0.2~r_sr_s$ as $β_\mathrm{star}-β_\mathrm{DM}=(-0.028 \pm 0.008)\log_{10}M_{200}/M_\odot + (0.690\pm0.010)$.