Jumps are an attempt to solve a math mistake in Modern Portfolio Theory. In the 19502-70s, economists were working on solving the variance-mean tradeoff. Furthermore, they needed to do so with punchcard computing. That radically restricted the set of computable, potential solutions. Both the normal distribution and the log-normal distribution are tractable with punchcard computing.
There are two problems, however. The first is one that most economists are unaware of. In 1958, a mathematician by the name of John White proved that there is no solution to an equation of the form $w_{t+1}=Rw_t+\epsilon_{t+1},R>1$ in Frequentist statistics. Of course, you would not invest money if $R\le{1}$. We will come back to this. The second is that Mandelbrot, beginning in 1963, started publishing articles that returns had heavy tails and could not be from a distribution with a variance. In other words, there is no variance-mean tradeoff because the first central moment does not exist.
Going back to the sixties and seventies with its heavy-tailed discussions and the results of the Fama-MacBeth work excluding the CAPM from empirical science, there was sort of a choice to be made. Embrace distributions without a mean and for which there was no undergirding of math for economists to work in, or decide for reality that there is a mean and just add jumps to try and cover the large shifts. That math was easily tractable.
That was an unfortunate choice. What makes it unfortunate is that the distribution of returns is $$R_{total}=R_G\times{\Pr(G)}+R_M\times{\Pr(M)}+0\times{\Pr(B)}+R_D\times{\Pr(D)}-R_L,$$ where $G$ denotes a going concern, $M$ denotes mergers, $B$ denotes bankruptcy, $D$ denotes dividends and $L$ denotes the lost return from liquidity costs. The distribution of $R_G$ is $$\left[\frac{\pi}{2}+\tan\left(\frac{\mu}{\sigma}\right)\right]^{-1}\frac{\sigma}{\sigma^2+(R_G-\mu)^2}.$$
Going back to White, from earlier, his proof was that the sampling distribution of the slope estimator was the Cauchy distribution. Models like Black-Scholes are built on either Ito or Stratonovich calculus. Both assume that all parameters are known. As such, it is a model built on parameters. If you didn't know them, then you cannot build models on them. You would have built them on sufficient statistics instead. As sufficient statistics are independent of the parameter, you wouldn't reference the parameter.
So models like Black-Scholes are valid if the parameters are known but invalid, as per White and a later generalizing article by Sen, if the parameters are unknown. There cannot exist a Frequentist solution to models like the CAPM or Black-Scholes as it is known to be impossible unless you use non-mean and non-variance based tools.
That opens up the possibility of a Bayesian solution, except the Bayesian solution ends up having no mean or variance because the results do not come out the same. That should serve as a deep warning as well.
All Bayesian estimators are admissible estimators. Frequentist estimators are admissible only in two cases. The first is that the Bayesian and the Frequentist solution are the same for every sample. The second is that the Frequentist solution matches the Bayesian solution at the limit. That is why $\bar{x}$ is an admissible solution to estimate $\mu$ for the normal distribution but $\frac{\sum(\sin(x_i))}{n-33}$ is not.
Although Frequentist estimators do not have to be admissible and sometimes are not, when they are not, there should be an investigation. It could be the model was derived incorrectly.
The skew in the volatility is an artifact of the algorithm. Consider, instead, for returns rather than option prices as it is simpler to discuss, where returns are $R=\frac{FV}{PV}-1,$ an algorithm such as $$\Pr(\sigma|X,\mu)=\int_{-1}^\infty\frac{\prod_{i=1}^I\left[\frac{\pi}{2}+\tan\left(\frac{\mu}{\sigma}\right)\right]^{-1}\frac{\sigma}{\sigma^2+(x_i-\mu)^2}\times{1}}{\int_{-1}^\infty\int_0^\infty\prod_{i=1}^I\left[\frac{\pi}{2}+\tan\left(\frac{\mu}{\sigma}\right)\right]^{-1}\frac{\sigma}{\sigma^2+(x_i-\mu)^2}\times{1}\mathrm{d}\mu\mathrm{d}\sigma}\mathrm{d}\mu,\forall\sigma\in\Re^{++}.$$
There will be a small natural amount of skew because the distribution should converge to the standard deviation ratio distribution but it will be small. It is related to the Snecdor's F distribution.
Note that I multiplied by $1$ and really should not have. Good priors exist for this but I didn't want to impose a prior on it as you should use your own.
The volatility skew is an artifact of the tool used to measure it and the fact that there is a non-existence theorem surrounding Ito models.