March 20th, 2006, 7:11 pm
zeta:solitons are not the natural symmetry for BS. You would use solitons if you had travelling wave solutions, whereas log S is drift-diffusive, with drift rt + convexity adjustment, and diffusivity vol^2.The FX convention is to work with BS call delta, D = N(d1) where d1 = (log F/K + 1/2 vol^2 t ) / vol sqrt(t) - that rescales the domain to [0, 1] and the S <-> D inversion is straightforward given that you know the BS vols. Advantage, already scaled to optionality so the domain is fixed. Disadvantage, structure in strike space manifests itself through thinning boundary layers in the delta space so you need adaptive discretization.Cuchulainn:You need to ensure that the implied pdf does not give you nonzero probability over negative strikes, otherwise your prices are off. This does happen with interest rates - even though nobody has ever written an option struck at negative rate, the fact we give them nonzero probability for hedging imposes a lot of structure on implied vol and pdf for positive strikes.Also, rather than truncating, maybe you can extrapolate? Something like domain rescaling with moving end points carrying the extrapolation information, where the extrapolation is BS mod your extrapolation of implied volatilities? Just a thought.