Given $X_1 \sim \exp(\lambda_1)$ and $X_2 \sim \exp(\lambda_2)$, and that they are independent, how can I calculate the probability density function of $X_1+X_2$?
I tried to define $Z=X_1+X_2$ and then: $f_Z(z)=\int_{-\infty}^\infty f_{Z,X_1}(z,x) \, dx = \int_0^\infty f_{Z,X_1}(z,x) \, dx$.
And I don't know how to continue from this point.
$$f_Z(z)=\int f_{X_1}(x)f_{X_2}(z-x)dx$$
Note that in your case the RHS has integrand $0$ if $z\leq0$ so that $f_Z(z)=0$ if $z\leq0$.
For $z>0$ we have:$$f_Z(z)=\int f_{X_1}(x)f_{X_2}(z-x)dx=\int_0^{z}f_{X_1}(x)f_{X_2}(z-x)dx$$
Work this out yourself.