To convert your calculations for $\exp(-\rho\,\sigma_z/2)$ to corresponding calculations for a boost in a general direction, first recall that in $\mathrm{SO}(1,\,3)$ a boost in a general direction has the matrix:
$$\Lambda = R\,\Lambda\,R^{-1}\tag{1}$$
where $R$ is a rotation that rotates a frame whose $z$-axis corresponds with the boost direction to the frame you wish to compute in. Since the Adjoint representation, which maps the double cover $\mathrm{SL}(2,\,\mathbb{C})$ to $\mathrm{SO}(1,\,3)$, is a homomorphism, you can use the same equation as (1) to find the preimage (or one of the two preimages) in $\mathrm{SL}(2,\,\mathbb{C})$ that maps to your general boost in $\mathrm{SO}(1,\,3)$.
So now you need to know the what a preimage of a rotation operator looks in $\mathrm{SU}(2)\subset\mathrm{SL}(2,\,\mathbb{C})$; unsurprisingly, a rotation in $\mathrm{SU}(2)$ is given by:
$$R(\theta, \hat{n}) = \exp\left(-i\,\frac{\theta}{2}\,\hat{n}\cdot\sigma\right) = \cos\left(\frac{\theta}{2}\right)-\sin\left(\frac{\theta}{2}\right)\,\hat{n}\cdot\sigma\tag{2}$$
So that your calculation yields:
$$\Lambda = \left(\cos\left(\frac{\theta}{2}\right)-\sin\left(\frac{\theta}{2}\right)\,\hat{n}\cdot\sigma\right)\,\left(\cosh\left(\frac{\rho}{2}\right)-\sinh\left(\frac{\rho}{2}\right)\,\hat{n}\cdot\sigma_z\right)\,\left(\cos\left(\frac{\theta}{2}\right)+\sin\left(\frac{\theta}{2}\right)\,\hat{n}\cdot\sigma\right)\tag{3}$$
which I'll leave you to multiply out.
But take heed that the calculation's result naturally wipes out the rotational degree of freedom one subtly introduces - any frame that is a rotation about the $z$ axis of any other frame wherein the boost is along the $z$ axis is itself a frame wherein the boost is along the $z$-axis. This is because if we replace the rotation in (2) by any rotation of the form $R(\theta, \hat{n})\,\exp(-\phi\,\sigma_z/2)$ for $\phi\in\mathbb{R}$, the $\exp(-\phi\,\sigma_z/2)$ and $\exp(+\phi\,\sigma_z/2)$ that are introduced into (3) both commute with the central $z$-direction boost in (3) and the calculation's result is unchanged. Thus you can apply these ideas with any rotation that aligns the $z$ axis as you need, and the calculation's result will be independent of this choice.