**Subject:**

gamma (second derivative of option price wrt S). Does it have a value at t = 0?

Can we model gamma as a PDE?

This is for t == 0 -> t = 1.0e-3

When t = T for initial conditions gives a better approximation..

- Cuchulainn
**Posts:**64700**Joined:****Location:**Drosophila melanogaster-
**Contact:**

gamma (second derivative of option price wrt S). Does it have a value at t = 0?

Can we model gamma as a PDE?

This is for t == 0 -> t = 1.0e-3

When t = T for initial conditions gives a better approximation..

"Compatibility means deliberately repeating other people's mistakes."

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

I’m pretty sure Peter Carr derives a PDE for gamma in his “Deriving derivatives of derivative securities” paper. Since delta at time zero is given by the Heaviside function, gamma at time zero is given by the Dirac delta function.

- Cuchulainn
**Posts:**64700**Joined:****Location:**Drosophila melanogaster-
**Contact:**

The reason for this red shift is that I forgot a convection term when differentiating the PDE for gamma. Now it's fine.Subject:

gamma (second derivative of option price wrt S). Does it have a value at t = 0?

Can we model gamma as a PDE?

This is for t == 0 -> t = 1.0e-3

When t = T for initial conditions gives a better approximation..

"Compatibility means deliberately repeating other people's mistakes."

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

Your gamma Dopplered to the left?The reason for this red shift is that I forgot a convection term when differentiating the PDE for gamma. Now it's fine.Subject:

gamma (second derivative of option price wrt S). Does it have a value at t = 0?

Can we model gamma as a PDE?

This is for t == 0 -> t = 1.0e-3

When t = T for initial conditions gives a better approximation..

- Cuchulainn
**Posts:**64700**Joined:****Location:**Drosophila melanogaster-
**Contact:**

Found the article by Peter C. Will read it. Actually, I attended a talk by Peter around 2001 in Broad Street (51?) when he actually introduced the topic!I’m pretty sure Peter Carr derives a PDE for gamma in his “Deriving derivatives of derivative securities” paper. Since delta at time zero is given by the Heaviside function, gamma at time zero is given by the Dirac delta function.

Now, I am fairly confident the maths is solid and the numerics are more robust than other stuff. I tried both Gaussian and box deltas and they are fine. Moving to multiple stat variables with this approach should be easy as well.

An open question is how to compute [$] \frac{\partial C}{\partial K} [$]; the strike [$]K[$] is 'hidden' in the payoff and is nowhere to be seen in the PDE.

Last edited by Cuchulainn on April 23rd, 2019, 4:18 pm, edited 2 times in total.

"Compatibility means deliberately repeating other people's mistakes."

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

Given that standard options are essentially symmetric in S and K, you could try a change of numeraire. Then the “strike delta” becomes the regular delta. This is usually easier to think about in an FX context.

- Cuchulainn
**Posts:**64700**Joined:****Location:**Drosophila melanogaster-
**Contact:**

Don't suppose you have a useful link? (BTW PC also excludes this case from his article).Given that standard options are essentially symmetric in S and K, you could try a change of numeraire. Then the “strike delta” becomes the regular delta. This is usually easier to think about in an FX context.

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

I had taken it for granted that everyone here knows how to differentiate. Looks like I was wrong.

- Cuchulainn
**Posts:**64700**Joined:****Location:**Drosophila melanogaster-
**Contact:**

Mathematica and Maple have a lot to answer for. Just like no one can do log tables not long division any more.I had taken it for granted that everyone here knows how to differentiate. Looks like I was wrong.

Last edited by Cuchulainn on April 23rd, 2019, 8:23 pm, edited 1 time in total.

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

- Cuchulainn
**Posts:**64700**Joined:****Location:**Drosophila melanogaster-
**Contact:**

http://citeseerx.ist.psu.edu/viewdoc/do ... 1&type=pdf

This is the article. I use the approach as in eq (1) but then for all greeks in my case, but not in the article..

Theorem 2 shows a formula for the derivatives wrt S written as series of Hermite polynomials. These are computationally difficult in a number of ways. The alternative here is that I only need to embed the greeks in a PDE (qualitative solution) and them embed them in a parabolic PDE system, much easier.

Eqs. (28)-(30); look like PDEs.

This is the article. I use the approach as in eq (1) but then for all greeks in my case, but not in the article..

Theorem 2 shows a formula for the derivatives wrt S written as series of Hermite polynomials. These are computationally difficult in a number of ways. The alternative here is that I only need to embed the greeks in a PDE (qualitative solution) and them embed them in a parabolic PDE system, much easier.

Eqs. (28)-(30); look like PDEs.

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

- FaridMoussaoui
**Posts:**507**Joined:****Location:**Genève, Genf, Ginevra, Geneva

You have to look into Dupire work.

An open question is how to compute ∂C/∂K; the strike is 'hidden' in the payoff and is nowhere to be seen in the PDE.

That makes it easy. You can still differentiate the PDE => [$]u = C_K[$] follows the same PDE with the differentiated payoff.

- Cuchulainn
**Posts:**64700**Joined:****Location:**Drosophila melanogaster-
**Contact:**

Excellent! I was over-engineering it! When I came back from the dojo it came in a flash.An open question is how to compute ∂C/∂K; the strike is 'hidden' in the payoff and is nowhere to be seen in the PDE.

That makes it easy. You can still differentiate the PDE => [$]u = C_K[$] follows the same PDE with the differentiated payoff.

Here is [$] \frac{\partial P}{\partial K} [$] as a function of the underlying [$]S[$], Is it correct(*)? It is the same as the BS PDE and the same BC/IC as detta.It looks a bit like a call delta.

(*) compared against [$]e^{-rT}N(-d2)[$]

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

- Cuchulainn
**Posts:**64700**Joined:****Location:**Drosophila melanogaster-
**Contact:**

Strike gamma (output as expected I suppose)

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl

- Cuchulainn
**Posts:**64700**Joined:****Location:**Drosophila melanogaster-
**Contact:**

This counterexample was salvaged from the ashes and saved for posterity.

Some related question to follow..

Good remarks by @fogsnow.

Some related question to follow..

Good remarks by @fogsnow.

David Wheeler

http://www.datasimfinancial.com

http://www.datasim.nl