Can the product of a Beta and some other distribution give an Exponential? Let $X$ follow a Beta distribution in $[0,1]$ with some parameters, and $Y$ be a non-negative random variable, independent of $X$.  
Is it possible for $Y$ to have a distribution such that $Z = XY$ follows an Exponential distribution? 
A bit of elaboration
By the formula for the density of the product distribution we can write
$$f_z(z) = \int_0^1\frac {1}{x}f_x(x)f_y(z/x)dx = \int_0^1\frac {dF_x(x)}{dx}\left(\frac {1}{x}f_y(z/x)\right)dx$$
$$=F_x(x)\left(\frac {1}{x}f_y(z/x)\right)|^1_0 - \int_0^1F_x(x)\left(-\frac {1}{x^2}f_y(z/x)-\frac {z}{x^3}f'_y(z/x)\right)dx$$
$$=f_y(z) + \int_0^1\frac {1}{x^2}F_x(x)\left(f_y(z/x)+\frac {z}{x}f'_y(z/x)\right)dx \tag{1}$$
To obtain the last step, I think that we have already assumed that $f_y(\cdot)$ goes to zero "fast enough".
If this is to be a density for an r.v. in $[0.\infty]$, it must be the case that 
$$\int_0^{\infty}f_y(z)dz + \int_0^{\infty}\int_0^1\frac {1}{x^2}F_x(x)\left(f_y(z/x)+\frac {z}{x}f'_y(z/x)\right)dx dz = 1 \tag{2}$$
But if $f_y(\cdot)$ is a densirty of a non-negative randomm variable, it follows that the first intergal should equal unity on its own (here $z$ is just the dummy variable of intergation). So a -necessary- condition is that 
$$\int_0^{\infty}\int_0^1\frac {1}{x^2}F_x(x)\left(f_y(z/x)+\frac {z}{x}f'_y(z/x)\right)dx dz  =0$$
where $F_x(x)$ is the cdf of the Beta distribution.
So the question is: is there a density $f_y(\cdot)$ such that it satisfies $(2)$ and makes $(1)$ equal to $\lambda e^{-\lambda z}$ for $\lambda>0$?
I am concerned mainly with existence, I don't expect to obtain $f_y(\cdot)$.
 A: Based on multiple answers to Distribution of $XY$ if $X \sim$ Beta$(1,K-1)$ and $Y \sim$ chi-squared with $2K$ degrees , let $X \sim \mathrm{Beta} \left(1,K-1 \right) $ and $Y$ be chi-squared with $2K$ degrees of freedom. 
Then the product is exponential with parameter $\frac{1}{2}$ .
A: A thought:
Since we are prescribing the distributions of $Z, X$, then $Z = XY \implies Y = Z/X$, so why not try to calculate $Y$ as a ratio distribution and see what we get?
We have 
$$f_y(y) = \int_0^1 x f_z(yx)f_x(x)dx$$
We assume that $X$ is a $Beta(a,b)$ and $Z$ is an $Exp(\lambda)$ (rate parameter). So 
$$f_y(y) = \frac {\lambda}{B(a,b)} \int_0^1 x \exp \{(-\lambda y)x\} x^{a-1}(1-x)^{b-1}dx$$
$$=\frac {\lambda}{B(a,b)} \int_0^1 x^a (1-x)^{b-1}\exp \{(-\lambda y)x\} dx$$
Looking up Gradshteyn & Ryzhik (2007, 7th ed.) p 347, eq. 3.383(1) and mapping coefficients, we get 
$$f_y(y) = \lambda \frac {B(b,1+a)}{B(a,b)} \cdot  {}_1F_{1}(1+a;1+a+b;-\lambda y)$$
where the function is Kummer's confluent hypergeometric function.
Is this a probability density over $[0,\infty]$?
