Abstract At long times, a fractional Brownian particle in a confining external potential reaches a non-equilibrium (non-Boltzmann)&#xD;steady state. Here we consider scale-invariant power-law potentials $V(x)\sim |x|^m$, where $m>0$, and employ the optimal fluctuation method (OFM) to determine the large-$|x|$ tails of the steady-state probability distribution $\mathcal{P}(x)$ of the particle position. The calculations involve finding the optimal (that is, the most likely) path of the particle, which determines these tails, via a minimization of the exact action functional for this system, which has recently become available. Exploiting dynamical scale invariance of the model in conjunction with the OFM ansatz, we establish the large-$|x|$ tails of $\ln \mathcal{P}(x)$ up to a dimensionless factor $\alpha(H,m)$, where $0<H<1$ is the Hurst exponent. We determine $\alpha(H,m)$ analytically (i) in the limits of $H\to 0$ and $H\to 1$, and (ii) for $m=2$ and arbitrary $H$, corresponding to the fractional Ornstein-Uhlenbeck (fOU) process. Our results for the fOU process are in agreement with the previously known exact $\mathcal{P}(x)$ and autocovariance. The form of the tails of $\mathcal{P}(x)$ yields exact conditions, in terms of $H$ and $m$, for the particle confinement in the potential.&#xD;For $H\neq 1/2$, the tails encode the non-equilibrium character of the steady state distribution, and we observe violation of time reversibility of the system except for $m=2$. To compute the optimal paths and the factor $\alpha(H,m)$ for arbitrary permissible $H$ and $m$, one needs to solve an (in general nonlinear) integro-differential equation. To this end we develop a specialized numerical iteration algorithm which&#xD;accounts analytically for an intrinsic cusp singularity of the optimal paths for $H<1/2$.&#xD;&#xD;