Pricing of Barrier Options on Underlying Assets with Jump-Diffusion Dynamics: A Mellin Transform Approach

A barrier option is an exotic path-dependent option contract where the right to buy or sell is activated or extinguished when the underlying asset reaches a certain barrier price during the lifetime of the contract. In this article we use a Mellin transform approach to derive exact pricing formulas for barrier options with general payoffs and exponential barriers on underlying assets that have jump-diffusion dynamics. With the same approach we also price barrier options on underlying futures contracts.


Introduction
A European call option is a financial contract that gives the holder the right, but not the obligation, to buy an underlying asset from the writer at an agreed strike price on a predetermined expiry date. A European put option is similar but gives the right to sell instead. Examples of the underlying asset, or simply the underlying, are stocks or futures contracts.
Options are mainly used for speculation and hedging. For example, an investor who believes that the share price for a certain stock is going to rise within the next month may invest by buying a call option on that stock. On the other hand, an investor who already owns shares of a certain stock may insure against a temporary fall in the share price by deciding to buy a put option to minimize the risk of a potential loss.
At the time the option contract is agreed upon, the holder must pay a certain amount, known as the premium or time-zero option price, to the writer. Option valuation, or option pricing, is the fundamental problem of determining a fair price for this premium. For European-style contracts, analytical expressions for the call and put premiums are given by the Nobel Prize-winning Black-Scholes formulas [1].
Call and put options can be characterized by their so-called payoff functions. We denote a payoff function by g : R`Ñ R, where R`" p0, 8q, which is typically piecewise linear. If SpTq is the asset price at the expiry date T and K is the strike price, then the call and put payoffs are gpSpTqq " pSpTq´Kq`and gpSpTqq " pK´SpTqq`, respectively, where pzq`" maxpz, 0q for any z P R. Thus the respective call and put payoff functions are gpxq " px´Kq`and gpxq " pK´xq`.
Options are attractive because they can be used to create a wide range of trading strategies characterised by different payoff functions. Suppose that K 1 ă K 2 ă K 3 ă K 4 . Denote the usual indicator function of a set A by 1 A , i.e., 1 A pxq " 1 if x P A and 1 A pxq " 0 if x R A. Some popular trading strategies with their corresponding payoff functions are given in Table 1 (see Hull [2] and Wilmott et al. [3] for more details): Table 1. Examples of options and their payoff functions.

Type
Payoff Function g(x) put pK´xqb ear spread pK 2´x q`´pK 1´x qc all px´Kqb ull spread px´K 1 q`´px´K 2 qd igital call 1 rK,8q pxq asset-or-nothing call x1 rK,8q pxq butterfly spread px´K 1 q``px´K 3 q`´2px´K 2 q`,´K 1´K3`2 K 2 " 0 iron condor px´K 1 q`´px´K 2 q`´px´K 3 q``px´K 4 q`,´K 1`K2`K3´K4 " 0 straddle px´Kq``pK´xqs trip px´Kq``2pK´xqs trap pK´xq``2px´Kqs trangle px´K 2 q``pK 1´x qC onsider the stochastic differential equation where S " {Sptq : t ě 0} is the underlying asset price process and W " {Wptq : t ě 0} is a Wiener process with respect to the risk-neutral measure. Here, the risk-free rate r, the dividend yield D, and the volatility σ are assumed to be constants with r, σ ą 0 and D ě 0. Denote the generic European option price at time t by Vptq and the corresponding payoff function by g. At expiry we therefore have VpTq " gpSpTqq. It is well known that Vptq " vpSptq, tq, where the option pricing function v " vpx, tq satisfies the Black-Scholes partial differential equation (PDE) Bv Bt`1 2 The premium is obtained by setting Vp0q " vpSp0q, 0q, where Sp0q is today's known asset price. While geometric Brownian motion assumed in the Black-Scholes asset price model (1) is convenient, it cannot capture many of the features of asset price returns, e.g., the skew/smile features of the implied volatility surface. In the absence of dividend payments, Merton [4] considered a jump-diffusion process that allows for the probability of the asset price to change at large magnitudes irrespective of the time interval between successive observations. The jumps in the asset price can be incorporated by introducing an additional source of uncertainty into the asset price dynamics. Empirical studies have revealed that the asset price is better described by a process with a discontinuous sample path (see, for instance, Rosenfeld [5], Jarrow and Rosenfeld [6], Ball and Torous [7], and Brown and Dybvig [8]). Cont and Tankov [9] showed that unlike standard diffusion models such as (1), jump-diffusion models produce rich structures of the distribution of asset returns and implied volatility surfaces.
To account for the possibility of instantaneous jumps in the asset price, Merton [4] proposed the following modification of (1) by assuming that the discontinuous jumps arrive as a Poisson process (here we incorporate dividend payments): where Y is a nonnegative continuous random variable with Y´1 denoting the impulse change in the asset price from Sptq to YSptq as a result of the jump, E is the expectation operator, and N " {Nptq : t ě 0} is a Poisson process with constant intensity λ and such that dNptq " 1 (respectively, dNptq " 0) with probability λdt (respectively, 1´λdt). It is assumed that Wptq, Nptq, and Y are independent for each t, and that the asset price jumps occur independently and identically. Analogous to the Black-Scholes PDE (2), it can be shown [4] that the European option pricing function v satisfies the Black-Scholes partial integro-differential equation (PIDE) where Ipv; f Y q is the integral operator defined by and f Y is the probability density function of Y. Using Mellin transform techniques adapted from Rodrigo and Mamon [10], Li and Rodrigo [11] studied (4) and found exact pricing formulas for European options with general payoffs such as those in Table 1. A PIDE analogue of Dupire's PDE was also derived and used to find an explicit formula for the implied volatility. In the same vein, Rodrigo and Goard [12] considered (4) with a time varying D and obtained exact pricing formulas for European options on discrete dividend-paying assets. A barrier option is an exotic path-dependent option contract where the right to buy or sell is activated (in the case of a knock-in barrier option) or extinguished (in the case of a knock-out barrier option) when the underlying reaches a certain barrier price during the lifetime of the contract. If the option expires inactive or extinguishes, then it may be worthless or there may be a cash rebate paid out. Since payoff opportunities are more limited, a barrier option is cheaper than a similar European option. Barriers are generally fixed but time-dependent barriers can be considered as well. The rationale for a barrier option is to provide a hedge at a lower premium than a conventional option.
There are many results in the academic literature on continuously monitored barrier options. One of the earliest dates back to the work of Merton [13], who gave a closed-form solution for the price of a continuously monitored down-and-out European call. One approach, mainly for fixed barriers, identifies pathwise hedging strategies for European-style derivatives that either uniquely determine or provide an admissible range for the barrier option price (see, for instance, Carr et al. [14] and Brown et al. [15]). A static hedge using calls and puts for a time-dependent single barrier option was given by Andersen et al. [16]. Their result also applies to linear diffusions with compound Poisson jumps but the hedging strategy depends on knowing the values of the barrier contract to be hedged at certain times before expiry. Geman and Yor [17] used a probabilistic approach for constant double barrier options in the Black-Scholes model. Kunitomo and Ikeda [18] introduced a method for pricing time-dependent barrier options in the Black-Scholes model with the help of the joint density of the asset and its maximum and minimum. Lattice methods were utilized by Boyle and Lau [19] and Ritchken [20], while finite difference and finite element methods were employed by Boyle and Tian [21] and Zvan et al. [22], respectively.
Recent results on time-dependent double barrier options include Fourier series expansions (Hui and Lo [23]), Green's functions (Dorleitner et al. [24]) and Laplace transforms (Pelsser [25]). Davydov and Linetsky [26] applied spectral methods to obtain constant double barrier option prices in constant elasticity of variance models. The boundary element method was used to derive an integral representation of the barrier option price in Guardasoni and Sanfelici [27] and Shen and Hsiao [28] under a Black-Scholes framework; in Guardasoni and Sanfelici [29] under stochastic volatility and jumps; and in Ballestra et al. [30] under a mixed fractional Brownian motion. Buchen and Konstandatos [31] proposed a method of images approach to price double barrier options with exponential barriers, extending the results by Buchen [32] for single fixed barriers.
Many options have as the underlying not the cash product but the corresponding futures contract, which is often more liquid and involves lower transaction costs [3]. Recall that a forward contract, made at time t˚, is a contract where the holder pays the writer the deterministic amount Fpt˚; T˚q (known as the forward price) at the delivery date T˚and then receives the stochastic amount SpT˚q at the same time. Under the asset price dynamics (1), it is known [3] that It can be shown [2,33] that when interest rates are deterministic, the price of a futures contract is the same as the price of a forward contract. One can then consider a European option, with payoff function g and expiry date T, on an underlying futures contract on an asset with price process S satisfying (1) and with delivery date T˚such that T ă T˚. Thus the payoff of a European option on a futures contract is gpFpT; T˚qq " gpSpTqe pr´DqpT˚´Tq q.
In the case of a call, where gpxq " px´Kq`, then after some modification to include D the so-called "Black-76" formula [34] is obtained.
The primary objective of this article is to employ a Mellin transform approach to price barrier options with general payoffs when the underlying is modeled by the jump-diffusion dynamics (3). The use of Mellin transforms in option pricing has been developed by the present author in a series of articles [10][11][12]35,36]. In the absence of jumps (i.e., λ " 0 in (4)), we recover the classical Black-Scholes framework (1). In fact, the corresponding barrier option pricing problem was considered by Guardasoni, Rodrigo, and Sanfelici in [36], where the Mellin transform was used to price single and double barrier options. Although the results in [36] are applicable to general time-dependent barriers, it is necessary to solve an associated linear Volterra integral equation of the first kind (or a coupled system of two linear Volterra integral equations of the first kind for double barriers). This integral cannot be solved analytically and one must resort to a numerical approximation. In this article we incorporate jumps as in (3) but consider only exponential barriers akin to that studied in Buchen and Konstandatos [31]. Note that this includes fixed barriers as a special case. We will also use a Mellin transform technique but avoid the introduction of Volterra integral equations of the first kind.
The secondary objective of this article is to obtain, with essentially the same amount of work, exact pricing formulas for barrier options with payoffs that include those in Table 1 but when the underlying is a futures contract, with the possibility of jumps in the corresponding asset. To my knowledge, even in the absence of jumps, barrier options on futures have not been previously considered in this general framework but the proposed Mellin transform approach can handle such exotic options, with or without jumps in the asset price dynamics associated with the futures contract.
The outline of this paper is as follows. In Section 2 we recall and also derive some preliminary results involving the generalized Black-Scholes kernel and jump function with the aim of using these to obtain the European option pricing function when the underlying behaves according to jump-diffusion dynamics. We also find image function solutions in the jump-diffusion case that will be needed to price barrier options. Here the reason for choosing an exponential barrier will be evident. We formulate the barrier option pricing problem in Section 4. Knock-out barrier options are considered in Section 4 while Section 5 deals with knock-in options. Illustrative examples for both types are also given. We discuss how to price options on futures in Section 6 and give brief concluding remarks in Section 7.

Preliminary Results
To ensure that our results are applicable to other problems, here we will consider the PIDE where c 0 , c 1 , c 2 P R with c 2 ą 0, and Ipv; f Y q is the integral operator defined in (5). Note that λ " 0 reduces L 0 to a Black-Scholes-type differential operator. If the underlying asset dynamics is given by (3), then we see from (4) that c 2 " σ 2 {2, c 1 " r´D, and c 0 "´r. On the other hand, for options on futures, if we make the substitution vpx, tq "vpx,tq into (4), wherex " xe pr´DqpT˚´tq andt " t (cf. (6)), then we obtain the PIDE Bv which is of the form (8) with c 2 " σ 2 {2, c 1 " 0, and c 0 "´r. We observe that when λ " 0, (9) simplifies to the PDE for pricing options on futures; see for instance Wilmott et al. [3]. We will return to the pricing of options on futures in Section 6. For the convenience of the reader we summarize here some results pertaining to the Mellin transform (see, for example, ( [37], pp. 362-363) or [38]). The Mellin transform of a function f : R`Ñ R is provided the improper integral converges at the complex number ξ. This transform is useful for Black-Scholes-type equations because of the following properties for the derivatives of f , namely, The convolution of f : R`Ñ R and g : R`Ñ R is defined as Note that the convolution operator is both commutative and associative. It follows that also known as the convolution property.

Generalized Black-Scholes Kernel and Properties
The following lemma motivates an extension (see (18) below) of the Black-Scholes kernel originally introduced in [10,35]: where Φ is the cumulative distribution function on a standard normal random variable. Hence defining it follows that Proof. Let α 1 , α 2 , β, γ 1 , and γ 2 be independent of x. Suppose that Z " Np0, 1q. Using the definition of the Mellin transform (10) and a property of normal random variables, we obtain Choosing and algebraic manipulations yield (14). Equation (15) is justified by continuity of the arguments of the Mellin transforms in (14). The properties in (16) can be easily shown by evaluating the derivatives on the right-hand side and comparing with (15).
Next we set a " c 2 pu´tq, b "´pc 1´c2 qpu´tq, c " c 0 pu´tq in Lemma 1. Define the auxiliary functions and the generalized Black-Scholes kernel Then (16) gives the Black-Scholes kernel identities Finally, from (14) we havê We see that (17)- (20) reduce to analogous results introduced in [10] in the special case when c 2 " σ 2 {2, c 1 " r´D, and c 0 "´r.

Properties of the Jump Function
We refer the reader to [11,12] for further details regarding this subsection. Construct a sequence where f Y is the probability density function of the random variable Y (see (3) and (5)) and δ is the Dirac delta function. It was shown in [11] that Defining the jump function J px, t, uq " e´λ pu´tq we see from (22) thatĴ pξ, t, uq " e´λ r1´EpY´ξ qspu´tq .
Example 1. Let Y be a lognormal random variable with parameters µ Y and σ 2 Y . Then Y " e X , where X " using a property of normal random variables. Hence (22) giveŝ and where we used (14) with a " jσ 2 Y {2, b "´jµ Y , and c " 0. Thus the jump function for a lognormal random variable Y with parameters µ Y and σ 2 Y is J px, t, uq " e´λ pu´tq δpx´1q`e´λ pu´tq

European Option Pricing Function for Underlyings with Jump-Diffusion Dynamics
We are ready to solve the PIDE (8) together with the final condition vpx, Tq " gpxq using the Mellin transform. Let v 0 " v 0 px, tq satisfy the final value problem for the PDE the Mellin transforms of v and v 0 with respect to x, respectively. We assume that the random variable Y is such that EpY´ξq is finite. Taking the Mellin transform of (8) and using the derivative properties in (11), we getL λv pξ, tq " Bv Bt pξ, tq´rppξq`λ´λEpY´ξqsvpξ, tq " 0, where p is defined in (20). Moreover,vpξ, Tq "ĝpξq, whereĝ is the Mellin transform of g. The solution of (27) is thereforev pξ, tq " e´p pξqpT´tq e´λ r1´EpY´ξ qspT´tqĝ pξq.

Example 2.
Continuing with Example 1 and recalling (29) and (26), we see that the Mellin transform of (32) iŝ vpξ, tq " e´λ pT´tq However, (20) gives´p If we define K j " K j px, t, uq such that its Mellin transform iŝ then (15) gives where The analogous Black-Scholes kernel identities (16) are Thus (33) can be expressed aŝ vpξ, tq " e´λ pT´tq Therefore if Y is a lognormal random variable with parameters µ Y and σ 2 Y , then the exact solution of (8) satisfying the final condition vpx, Tq " gpxq is vpx, tq " e´λ pT´tq with K j defined as in (34).

The Black-Scholes kernel identities in
Therefore (37) becomes vpx, tq " Ke pc 0´λ qpT´tq where z j,1 and z j,2 are given by (35). Note that in the special case when λ " 0, c 2 " σ 2 {2, c 1 " r´D, and c 0 "´r, this recovers the Black-Scholes formula for a European put. A general piecewise linear payoff function g, corresponding to any of the option strategies given in Table 1 among others, can similarly be considered since the integrals appearing in (37) can always be evaluated with the help of the Black-Scholes kernel identities (16) (or (36) if Y is lognormal).

Image Function Solutions
Here we construct so-called image function solutions of (8) that will be used to "piece together" the barrier option pricing functions later on (cf. [31] for the special case when λ " 0, c 2 " σ 2 {2, c 1 " r, and c 0 "´r).
Suppose that B " Bptq is an exponential barrier, i.e., where µ P R. Note that (38) includes the fixed barrier case if we set µ " 0, so that Bptq " BpTq for all 0 ď t ď T.
Let v 1 " v 1 px, tq be any solution of (8). We claim that its image function is also a solution of (8) for some α P R to be determined. Rather than use a "brute-force" verification by substituting v 2 into (8), we follow a more instructive Mellin transform route. We know that v " vpx, tq is a solution of (8) if and only if its Mellin transformv "vpξ, tq is a solution of (27). Note that we are not imposing any final condition like vpx, Tq " gpxq here but are considering solutions of the PIDE and the transformed PIDE only.
Our task is to introduce (40) to simplify the right-hand side. Now let us define p λ pξq " ppξq´µξ´λEpY´ξq and assume that there exists α such that p λ pα´ξq " p λ pξq for every ξ.
Recalling (40), we deduce thatv 2 is a solution of (27) and hence v 2 given in (39) is also a solution of (8), as was to be shown. We remark that image functions become solutions only for exponential (or fixed) barriers.
Let us now take a closer look at the assumption (42). This imposes a condition on α and/or the distribution of Y. A pair of sufficient, but not necessary, conditions for (42) to be true is ppα´ξq´µpα´ξq " ppξq´µξ, EpY´ξq " EpY´p α´ξq q for every ξ. (43) Of course, if λ " 0, then the first relation above already implies (42).
In the absence of dividend yields (i.e. D " 0), this is precisely the value of α obtained in [31] to arrive at the mirror function solution (39) of the Black-Scholes PDE (2). Furthermore, for a fixed barrier (i.e., µ " 0), this result is well known [3]. Example 4. Now suppose that we include jumps (i.e., λ ‰ 0). Then α is still given by (44) but the second condition in (43) is more restrictive as it makes further assumptions about the distribution of Y. A particular case is when Y is lognormal with parameters µ Y and σ 2 Y . Then EpY´ξq " EpY´p α´ξq q for every ξ if and only if µ Y and σ 2 Y are such that µ Y " ασ 2 Y {2 (see (26)). If we assume that µ is given, and choose α in (44) and let µ Y " ασ 2 Y {2, then (42) holds. Alternatively, if µ Y and σ 2 Y are arbitrary but given, then we take α " 2µ Y {σ 2 Y and µ " c 1´c2´α c 2 from (44). This would also imply that (42) is also true. As pointed out above, (43) is not a necessary condition so there may exist other α and Y when λ ‰ 0 such that (42) is true.

Formulation of the Barrier Option Pricing Problem for Underlyings with Jump-Diffusion Dynamics
Here we formulate the barrier option pricing problem associated with the PIDE (8). The active domain is defined [32,36]
The active domain for down-and-out and down-and-in options is Aptq " pBptq, 8q. Similarly, the active domain for up-and-out and up-and-in options is Aptq " p0, Bptqq. Let v 0 " v 0 px, tq, v 1 " v 1 px, tq, and v e " v e px, tq be the solutions of the auxiliary problems respectively. We observe that v 1 and v e are solutions of the same PIDE (8) but with different final conditions (we associate v e with the "usual" European option pricing function). Moreover, v 0 is essentially v 1 but with λ " 0. Equation (31) with g replaced by g1 ApTq gives v 0 px, tq " An analogous expression for v e can also be obtained from (31) (with λ " 0) and (32) (with λ ‰ 0) by considering the payoff function g on the entire R`. and v 1 px, tq " e´λ pT´tq The barrier option pricing function v " vpx, tq satisfies the PIDE Note that the PIDE is only considered in the active domain. The final condition is vpx, Tq " gpxq, x P ApTq (52) for knock-out options and vpx, Tq " 0, x P ApTq for knock-in options. For knock-out options the condition at the barrier is vpBptq, tq " 0, 0 ď t ă T, while for knock-in options the barrier condition is vpBptq, tq " v e pBptq, tq, 0 ď t ă T.
Further analysis necessitates that we consider knock-out and knock-in options separately.

Pricing of Knock-out Barrier Options
A knock-out barrier option pricing problem is to solve We claim that is the solution of (56). As v out is a linear combination of a solution and its image function solution, linear superposition implies that v out satisfies the PIDE in (56). It is clear that v out pBptq, tq " 0 for all 0 ď t ă T. Furthermore, by construction v 1 px, Tq " gpxq1 ApTq pxq for x ą 0; hence v 1 px, Tq " gpxq for x P ApTq. If we can show that then v out px, Tq " gpxq for x P ApTq in (57) and this would verify the final condition in (56). We remark that this is the reason for the choice of the final condition of v 1 in (46). Using (48), we have Since all of the terms in the above summation are zero when t " T, it suffices to look at lim tÑT´v 0 T gpyq dy and show that the limit is zero for all x P ApTq. To evaluate the limit, we first investigate the bounds for Our goal here is to show that the integral on the right-hand side tends to zero as t Ñ T´and we would be done. As g is assumed to be piecewise linear, it is either (i) bounded or (ii) unbounded but gpxq " Opxq as x Ñ 8.

g Is Bounded
If g is bounded, then there exists M ą 0 such that |gpxq| ď M for all x ą 0, from which we deduce that For a down-and-out barrier option when Aptq " pBptq, 8q, Equation (19) yields which tends to zero as t Ñ T´since BpTq{x ă 1.
In the up-and-out barrier option case when Aptq " p0, Bptqq, from (19) we obtain which also tends to zero as t Ñ T´since BpTq{x ą 1 this time. Therefore for the case when g is bounded we see that
Consider a down-and-out barrier option, i.e. Aptq " pBptq, 8q. Then On rBpTq, x 8 s we know that g is bounded, so there exists M 1 ą 0 such that |gpyq| ď M 1 for BpTq ď y ď x 8 . Thus (19) implies that which tends to zero as t Ñ T´since BpTq{x ă 1 and BpTq{x 8 ă 1. On px 8 , 8q, we deduce from (19) that which also tends to zero as t Ñ T´because BpTq{x ă 1 and BpTq{x 8 ă 1. Hence lim tÑT´ ApTq for the down-and-out case. Now let us look at an up-and-out barrier option, so that Aptq " p0, Bptqq. Since g is piecewise linear, hence bounded on r0, BpTqs, there exists M 2 ą 0 such that |gpyq| ď M 2 for 0 ď y ď BpTq. Therefore which tends to zero as t Ñ T´since BpTq{x ą 1 in this case. Therefore for the case when g is unbounded but gpxq " Opxq as x Ñ 8, we conclude that Bptq 2 x , t " 0, x P ApTq.
For definiteness we assume that Y is lognormal with parameters µ Y and σ 2 Y . Then (50) gives v 1 px, tq " e´λ pT´tq while the down-and-out barrier option pricing function from (57) is Case (i). BpTq ă K Using (36), we have Case (ii). BpTq ě K Using (36) and noting that y ě BpTq implies y ě K, we have , t, T and therefore v 1 px, tq " xe pc 0`c1´λ qpT´tq , t, T .

Pricing of Knock-in Barrier Options
A knock-in barrier option pricing problem is to solve The previous analysis for knock-out options can be used for knock-in options.
We claim that v in is the knock-in barrier option pricing function. By linear superposition it is clear that v in satisfies the PIDE in (59). At the boundary there holds v in pBptq, tq " v e pBptq, tq´v out pBptq, tq " v e pBptq, tq, 0 ď t ă T, while at the expiry date we see that v in px, Tq " v e px, Tq´v out px, Tq " gpxq´gpxq " 0, x P ApTq.
This proves the claim.

Example 7.
Continuing with Example 6, let us price a down-and-in asset-or-nothing call. All we need is to determine v e , which from (37) and (36) is v e px, tq " e´λ pT´tq The corresponding down-and-out asset-or-nothing call option pricing function v out is given in (58). Note that there will also be two cases for v 1 here: BpTq ă K and BpTq ě K. Therefore the down-and-in asset-or-nothing call option pricing function v in is obtained from (60).

Pricing of Barrier Options on Futures under Jump-Diffusion Dynamics
If the underlying asset pays a constant dividend and has dynamics decribed by (3), then in (57) and (60) we simply take c 2 " σ 2 {2, c 1 " r´D, and c 0 "´r. In the absence of jumps, setting λ " 0 and D " 0 recovers the well-known results for fixed barriers [32] and exponential barriers [31].
On the other hand, if the underlying is a futures contract, as described in Section 2, we take c 2 " σ 2 {2, c 1 " 0, c 0 "´r,x " xe pr´DqpT˚´tq ,t " t in (8), andv "vpx,tq in (32) would give the standard European option pricing function. For example, if Y is lognormal with parameters µ Y and σ 2 Y , then (37) (which is (32) in the lognormal case) gives vpx,tq " e´λ pT´tq where K a j px,t, Tq "x , z a j,2 px,t, Tq " When λ " 0, D " 0 and gpxq " px´Kq`, (61) reduces to the "Black-76" formula [34] for a call option on a futures contract. The determination of the barrier option pricing function can then proceed as given in Sections 4 and 5.
Alternatively, instead of introducing variable transformations in (8), we can price barrier options on futures as follows. Again it suffices to consider the option pricing function of a standard European option with payoff function g on an underlying asset with jump-diffusion dynamics (3). Suppose that Y is lognormal with parameters µ Y and σ 2 Y . Take c 2 " σ 2 {2, c 1 " r´D, and c 0 "´r in (8). The variable 'x' here is the placeholder for the asset price while 'x' above is for the futures price. Recall that for a forward contract, the simple contingent claim Π h pT˚q with contract function h (see [33] for an explanation of the terminology) is now where Fpt˚; T˚q is the forward price to be determined and hpxq " x´Fpt˚; T˚q. Note that the price process for the contingent claim at any time t can be described by Π h ptq " vpSptq, tq, where v " vpx, tq satisfies (8). As it costs nothing to enter into a forward contract at time t˚, we must have Π h pt˚q " 0. Using the pricing formula (37) where we used (36) to evaluate the integrals. The condition Π h pt˚q " vpSpt˚q, t˚q " 0 implies that the forward price is Fpt˚; T˚q " Spt˚qe pr´DqpT˚´t˚q . This is the same as (6) but in the jump-diffusion case. Moreover, the prices of a futures contract and a forward contract are the same if the interest rate is deterministic since this result is model free [2]. As seen in (7), we can therefore price a European option on a futures contract using (37) by substituting gpxe pr´DqpT˚´Tq q, i.e., the European option pricing function v " vpx, tq (with a slight abuse of notation) on a futures contract is vpx, tq " e´λ pT´tq where K b j px, t, Tq " Using (62) and (64), it is straightforward to verify that K a j px,t, Tq " K b j pxe pr´DqpT˚´Tq , t, Tq.
Therefore the European option pricing function on a futures contract can be obtained from either (61) or (63). This is then used to price barrier options in Sections 4 and 5.

Discussion and Concluding Remarks
The knock-out barrier option pricing formula (57) and knock-in barrier option pricing formula (60) are solutions of the general PIDE (8). We can therefore price barrier options with general payoffs and exponential barriers when the underlyings exhibit jump-diffusion dynamics. Some of the more popular trading strategies are given in Table 1, but any barrier option with a piecewise linear payoff can be priced because the resulting integrals in (57) and (60) are evaluated explicitly with the help of the fundamental kernel identities in (19).
For more general barriers we have to use the results in [36] (with λ " 0), which also employs a Mellin transform approach. The tradeoff is that the pricing formulas in [36] are semi-analytic since they involve linear Volterra integral equations of the first kind which have to be evaluated numerically. Although we do not consider them here, double exponential barriers can also be studied. Indeed, in the absence of jumps, Equation (7) in [31] expresses the option pricing formula as a doubly infinite series of image function solutions. In principle, with the addition of jump-diffusion dynamics, one can show that an analogous series expansion is also a solution of (8) by showing that its Mellin transform is a solution of (27), just like what was done in Section 2. 4 The use of the Mellin transform in the pricing of financial derivatives where the underlying dynamics are governed by geometric Brownian motion has proven to be very powerful and provides a useful tool in the quant's toolbox [10][11][12]35,36]. In two articles currently under review, the author has applied this tool to price perpetual American options with general payoffs, as well as a combined Mellin-Laplace transform approach to price American options with general payoffs (the latter article is joint work with Mamon). Other pricing problems that are currently being investigated by the author are lookbacks, compounds, and Parisians, among others, with the assumption of jump-diffusion dynamics. A far-reaching goal is to be able to extend the results to include stochastic volatility and early exercise features to other pricing problems. However, it should be remarked that Mellin transforms are not as useful when considering other asset price dynamics precisely because the derivative properties in (11) are not valid anymore. In these scenarios a different integral transform has to be used, although it is not clear what the appropriate transform should be and therefore the pricing problem has to be handled on a case-to-case basis.