narrow line seyfert 1 galaxies ( osterbrock & pogge 1985 ) form a special class of active galactic nuclei .
narrowness of the optical lines frequently correlates with the steepness of the soft x - ray spectra ( boller , brandt & fink 1996 ) hinting for a connection of this class with soft state galactic sources .
the overall spectra are sometimes well described by a disk emission superimposed on a broad ir - x - ray power law component ( puchnarewicz et al .
2001 ) but in most cases the soft x - ray slope seems to be too shallow for an extension of a disk black body emission .
such a steep soft x - ray spectrum of a roughly power law shape may form either as a result of comptonization ( e.g. czerny & elvis 1987 ) by a moderately hot material ( in comparison with the very hot plasma responsible for the hard x - ray component in normal seyfert 1 galaxies ) or as a result of reflection by a strongly ionized material ( e.g. czerny & .
zycki 1994 , czerny & dumont 1998 )
. the established presence of warm absorber seen in many of those objects ( collinge et al .
2001 , kastra et al .
2002 , kaspi et al .
2002 , kinkhabwala et al .
2002 , leighly et al .
2002 , yaqoob et al .
2002 ) complicates the analysis .
however , with the extremely good spectral resolution of _ chandra _ and _ xmm _ the issue may be addressed through the careful analysis of the existing spectral features .
ton s180 ( @xmath0 , corresponding to 286 mpc for hubble constant 65 km s@xmath1 mpc@xmath1 ) is an extreme case of a nls1 with its low value of the fwhm of @xmath2 line ( @xmath3 km s@xmath1 ; wisotzki et al .
1995 ) and steep soft x - ray spectrum of photon index @xmath4 ( comastri et al .
1998 ) the galactic absorption toward the object is low ( @xmath5 @xmath6 ; stark et al .
the source was extensively observed in all energy bands . the broad band spectral energy distribution ( sed ) determined by turner et al .
( 2002 ) indicates that most of the energy is emitted in the band @xmath7 ev .
furthermore , _ hst _
data show typical emission lines as expected for a seyfert 1 galaxy , and the lack of absorption features .
high resolution spectra obtained by _ fuse _ reveal uv absorption by ovi and the absence of neutral hydrogen absorption , indicating high ionization state for the absorbing gas .
high resolution x - ray spectra of this object were recently obtained with _ xmm - newton _ and _ chandra _ satellites .
_ xmm - newton _ spectrum was steep and featureless ( vaughan et al .
_ chandra _ data also appeared predominantly featureless and well represented by a single power law with some soft excess ( turner et al .
2001 , turner et al .
the shape of the soft excess was fitted by a diskbb model , however plotted residua hinted for the possible presence of sharp spectral features ( fig . 1 turner et al .
2001 ) .    therefore , in the present paper we reanalyze _ chandra _ data in some detail ( sec .
[ sec : cha ] and [ sec : spe ] ) , using a new response matrix acisleg1d1999-07-22rmfn0002.fits and including any effects related to the contamination layer on acis - s .
we note that the observation of ton s180 was performed on day 143 of the chandra mission when the effects of the contaminant were much smaller than currently observed .
we detected several weak narrow absorption lines due to the material located at the line of sight to the source .
also , several individual disk lines are fitted with equivalent widths in rough agreement with the predictions of the model of illuminated accretion disk spectrum presented by raska et al .
( 2002 a , b ) .    in sec .
[ sec : apl ] we discuss existing models ( kaastra et al .
2002 , kinkhabwala et al .
2002 and others ) for fitting absorption by warm , partially ionized gas present in many agn .
we point out , after krolik ( 2001 ) , that the warm absorber may not be the slab of constant density matter , but it is rather the distribution of regions with different densities and temperatures . such a multi - phase mixture of gas can exist under constant pressure , as a result of thermal instabilities ( krolik , mckee & tarter 1981 ) .
till now , the multi - phase models of warm absorber were not constructed , because popular photoionization codes ( xstar , cloudy ) , compute radiative transfer through the constant density matter . in sec .
[ sec : mod ] , we present new set of models computed by the code titan ( dumont , abrassart & collin 2000 ) , which solves radiative transfer through the matter in non - lte approach .
introducing assumption of constant pressure ( not density ) of warm absorber , we obtain spectra with many absorption lines from the gas with strong temperature and density gradient . our code is suitable for any optical thickness of the absorber , nevertheless , we consider regions with total column density of order of @xmath8 , as suggested by krolik ( 2001 ) to be a typical property of warm absorber .    in sec .
[ sec : bbc ] we construct two global disk corona models trying to predict the slope of observed continuum .
the model with higher black hole mass ( @xmath9 ) is more promising in the sense that emission lines from the disk are expected in the spectra , since hot skin does not cover disk completely .    sec .
[ sec : dis ] is devoted to the discussion , while conclusions are listed in sec .
[ sec : conc ] .
ton s180 was observed with the _ chandra _ satellite in december 14 , 1999 , for about 80 ksec in low - energy transmission grating in acis - s / letg configuration , as described by turner et al .
( 2001 ) and turner et al .
the resolution of letg corresponds to an @xmath10 km s@xmath1 for the range of energies from 0.1 up to 6 kev .
the source luminosity at this epoch ( 210 kev unabsorbed flux equal @xmath11 erg s@xmath1 @xmath6 and 0.310 kev unabsorbed flux equal @xmath12 erg s@xmath1 @xmath6 ) was comparable to that during the _ beppo - sax _ observation ( 210 kev : @xmath13 erg s@xmath1 @xmath6 ; comastri et al . 1998 ) and _ xmm - newton _ observation ( 0.3 - 10 kev : @xmath14 erg s@xmath1 @xmath6 ; vaughan et al .
2002 ) and somewhat fainter than during the 12 days _ asca _ monitoring ( 210 kev : @xmath15 erg s@xmath1 @xmath6 ; romano et al . 2002 ) .
we obtained the data from the chandra archive .
the data were processed with the standard cxc data processing pipeline ascdsver 6.7.0 .
we analyzed the data using ciao version 2.2 . and follow grating analysis thread ( cxc web page ) to extract the spectrum and create appropriate effective area files ( arf files ) .
we combined both sides of the first order spectrum .
we used calibration data base version caldb  v.2.14 , in particular we applied acisd1997-04-17qen0003.fits , quantum efficiency file and acisleg1d1999-07-22rmfn0002.fits as response matrix .
this formally allows us to study the entire energy range between 0.12 till 12.2 kev , although calibration uncertainties are higher below 0.3  kev , while the effective area drops significantly above @xmath16 kev . after removing bad pixels and rebinning the data to achieve at least 20 counts per channel we examined the energy between 0.18 and 6  kev .
all spectral fits have been performed with the xspec 9.0 fitting package .
since the begining of the _ chandra _ mission acis - s detector has been accumulating a layer of a contaminant which affects any measurements at low energies ( below @xmath17  kev ) .
detailed effects of the contaminant are still being investigated by the chandra x - ray center , however , it is clear that the overall degradation of acis quantum efficiency and presence of additional absorption features is time dependent .
ton s180 was observed at the day 143 of the _ chandra _ mission and based on the currently available calibration data we expect that the instrument shows maximum of @xmath1840 - 15 % degradation within 0.28 - 0.6  kev range ( i.e. maximum @xmath19 % in 0.28  kev and @xmath20 % in 0.6  kev ) . in the following we identify the possible absorption features due to the absorber in ton s180 using the entire available energy range .
we note here that some of the features at energies close to the edges of the contaminant ( c , n , o , f k - edge ) may be affected .
we use the available model , acisabs , to include the contamination effects in the analysis .
in the first step , we fit the @xmath21 value of galactic absorption together with a broken power law continuum . our best fit value of @xmath22 @xmath6 is slightly larger than galactic absorption @xmath23 @xmath6 determined by dickey & lockman ( 1990 ) and used by comastri et al .
( 1998 ) or than @xmath24 @xmath6 determined by stark et al .
( 1992 ) and used by turner et al .
( 1998 , 2001 ) .
all fitted parameters for `` plain continuum '' are listed in table  [ tab0 ] .    in the next step
we fit broken power law with value of galactic absorption fixed at @xmath25 @xmath6 ( see table  [ tab0 ] second line ) .
it results in slightly higher @xmath26 , but the power law index below 2 kev is @xmath27 , and agrees with the best fit slope @xmath28 presented by comastri et al .
1998 for beppo sax data .
lower value obtained by turner et al .
( 2001 and 2002 ) , equal to @xmath29 resulted from the use of additional component to model the soft x - ray excess .
the energy break obtained by us in both models ( about @xmath30 kev ) is somewhat lower than the value determined by comastri et al .
1998 ( @xmath31 ) , but the spectral slope above this energy break , @xmath32 agrees with results comastri et al .
1998 .    since both broken power - law indices are very similar to those obtained from other satellites we adopt second model ( fixed galactic absorption ) for further analysis , leaving @xmath33 per 1682 d.o.f .
( reduced @xmath34 ) .
llllll    plain cont . &
@xmath35 & @xmath36 & @xmath37 & @xmath38 & 1185/1681 + fixed @xmath39 & 1.52 & @xmath40 & @xmath41 & @xmath42 & 1197/1682 + with lines & @xmath43 & @xmath44 & @xmath45 & @xmath46 & 932/1624 +    after fitting continuum , there is still some structure in the residuals of the soft x - ray spectrum , suggesting the presence of absorbing matter along the line of sight toward the source .
therefore , we have tried to determine column density of this gas in the simplest way , by fitting ovii , oviii , and cvi edges often seen in seyfert galaxies ( mrk 509 yaqoob et al .
2002 , ngc 3783 kaspi et al .
surprisingly , such edges appeared to be very weak in the data of ton s180 .
we were able only to give upper limits for @xmath47 @xmath48 , and @xmath49 of order of @xmath50 , which practically indicates the lack of those edges .
x - ray absorption lines appear to be common feature in many agn .
the widely acceptable explanation is the existence of warm , ionized gas which is illuminated by the energetic radiation from nucleus . trying to find properties of warm absorber we need proper tools to fit those data .
till now , models of absorption by the partially ionized plasma available in xspec ( absori , photo , ) do assume slab with constant density .    in absori ,
the radiative transfer of lines is not computed at all , only photoionization cross sections are determined and absorption `` on the spot '' is computed @xmath51 ( done et al .
constructed by kinkhabawa et al .
( 2002 ) local photo model ( http://xmm.astro.columbia.edu/research.html ) neglects photo - electric absorption , but allows for individual line absorption , treating each ionic column density separately with an initially unabsorbed power law .
this method produces consistent results only when absorbing matter is optically thin , and transmitted spectra do not display photo - electric edges .
furthermore , they have assumed that plasma is optically thin to reemitted photons .    in the aim to find evidences for warm absorber in ton s180 , we have applied absori to our data . allowing for three parameters to be fitted in this model , @xmath26 decreases by 38 ( now reduced @xmath52 ) .
the fit indicates that the medium has column density @xmath53 @xmath6 and ionization parameter ( @xmath54 ) for the fixed default temperature @xmath55 k. this model , however , leaves considerable residual trends and does not allow for a deeper insight into the ionization state of the absorbing material .
it is not surprising since the model is based on simple computations of photoionization cross sections without accurate transfer through the medium .
we have applied photo to our analysis .
the model again indicates very low column densities of each ion , about @xmath56 @xmath6 for cvi , ovii , and oviii .
we conclude that most probably , the warm absorber in ton s180 is more complicated than optically thin constant density slab . it may be the stratified density gas , where accurate radiative transfer , including all photoionization processes , should be calculate to model transmitted spectrum .
such models are not available for fitting the data yet , so two or three separate absorbers with different constant densities and ionization states are used together to explain spectral features in some agn ( kaastra et al .
2002 , kaspi et al .
this motivates us to present in sec .
[ sec : mod ] the model of the multi - phase gas with stratified density and temperature , where the radiative transfer of incident power law is carefully computed .
the model is not ready for fitting yet , but gives interesting results based on equivalent widths of various lines .
observed line properties in the ton s180 data are determined in the section below .
more advanced modeling was done with the code xstar ( developed by kallman & krolik 1999 ) which allows to estimate the total column densities of absorbed gas .
models are computed assuming constant density slab .
such models are very useful for reconstruction of the flat part of curve of growth ( yaqoob et al .
2002 ) , and for subsequent estimation of turbulent velocities of the absorber .
nevertheless , at the last stage of analysis separate gaussian components are fitted to the data .
we have successfully fitted several absorption lines using standard gaussian model for each line separately ( see figs .
[ fig : cont ] , [ fig2:cont ] , [ fig3:cont ] , [ fig4:cont ] , [ fig5:cont ] , [ fig6:cont ] , and [ fig7:cont ] .
* gaussian fit requires three free parameters , energy of the line centroid , gaussian width and normalizations , which has to be negative in case of absorption lines .
* we include in the fit only those features which are significant .
all absorption lines are presented in table  [ tab1 ] . for each line
we list its central energy in kevs , equivalent width ( ew ) in evs , the decrease of @xmath26 , best identification with observed and rest frame energy , and with possible blends .    some of those lines , especially the softest ones , may be caused by dust presented on chandra mirrors .
this should show up a common feature in all objects even those with different redshifts . to check whether lines are real or instrumental , we compare observations of ton s180 with _
chandra _ data of quasar 3c273 observed two weeks latter with the same configuration as ton s180 .
exemplary residuals ( for soft range between 0.2 - 0.24 kev ) for both objects are shown in fig .
[ por:22 ] . in both cases local power - law
was fitted with fixed value of galactic absorption appropriate for each object . the same features around energy 0.22 kev are present as well in 3c273 with redshift 0.156 as in ton s180 .
we repeated such analysis in the whole energy range .
all possible instrumental lines are labeled in table  [ tab1 ] .
some of the detected features are seen only in ton s180 spectrum and therefore seem to be real , intrinsic to the source .
for example , in fig  [ por:44 ] we show that line cvi with observed energy 0.4428 kev which is detected in ton s180 data is not seen in 3c273 .
lines are very faint .
they usually have equivalent width less or about 2 evs , and several of them are just marginal detections , as clearly seen from the line ew error , given at 90% confidence level ( * each error is computed for the one parameter of interest * ) .
only two of them are stronger , but from theoretical predictions by analysing our modeled spectra from sec  [ sec : mod ] , they are not blends .
all possible blends according of our modeled specta obtained with spectral resolution 100 are listed in table  [ tab1 ] .
there is therefore no major discrepancy between the _ chandra _ results and the _ xmm _ observations ( vaughan et al .
2002 ) where only upper limits of order of 1 ev at energies 0.35 - 0.7 kev were provided .
most lines originate from hydrogen or helium like ions indicating that matter is highly ionized .
there is an absence of lines from mildly ionized spices .
interestingly , we found two different transitions from the same element at the given ionization level , ( cvi@xmath57 and cvi@xmath58 where i  j denotes number of levels from and to which the transition proceeds ) .
this is helpful in the aim to estimate the optical thickness of cvi ions and then total optical thickness of the warm absorber .
if the medium is optically thin in those lines ( i.e. weak line case ) , we are on the linear part of the curve of growth , where equivalent width is proportional to the column density @xmath59 of the absorbed ion : @xmath60 where @xmath61 is the oscillator strength transition and stimulated emission has been set equal to unity .
thus , for the given population of ions at the same ionization level , the ratio of equivalent widths is given by : @xmath62 since the oscillator strengths are strictly atomic features independent from the radiation field and are calculated theoretically we use them to find the ratio @xmath63}{ew[{\rm c vi}_{(1 - 5 ) } ] } =   29.86 .
\label{equ : tra}\ ] ] on the other hand , the same ratio determined for ton s180 ( see table  [ tab1 ] ) equals @xmath64 .
error is computed assuming most pessimistic case , when numerator has its maximal value @xmath65 , and denominator has its minimal value @xmath66 .
although the error is huge , the value of the ratio determined from observation is still far from theoretical predictions for optically thin medium ( eq .  [ equ : tra ] ) .
our theoretical spectra ( see sec .  [
sec : mod ] ) suggest that c@xmath57 line is possibly a blend with mgxii , sixiv and most strongly with cv .
we have computed that real ew of the c@xmath57 line can be even up to four times lower than the ew of the whole blend .
but this fact even lowers the observed ratio of ews of two carbon lines .
the straight conclusion is that the emitting region is optically thick for those lines and the transfer in lines should be taken into account for proper modeling of the warm absorber .
it is also possible that the identification of the lines is not correct or the line is a blend of the carbon line with another transition .
nevertheless , we further explore the possibility that the carbon line ratio is determined correctly .
therefore , the equation  [ equ : column ] gives only a lower limit for the ion column density , equal to @xmath67 @xmath6 from the cvi@xmath57 line and @xmath68 @xmath6 from cvi@xmath58 , which under assumption of all carbon being in the form of cvi ions and with solar abundances translates into the lower limit for the warm absorber column equal @xmath69 @xmath6 and @xmath70 @xmath6 , correspondingly .
exact determination of the column density is complicated . at the flat curve of growth
we need additional information about velocity dispersion @xmath71 to indicate column densities .
in the case of considered data this is very difficult , because @xmath72 error of full width of half maximum ( fwhm ) is too high to estimate @xmath71 parameter .
for instance for sxii line , presented in fig .
[ fig3:cont ] , @xmath73 km / s .
this fact also precludes us to integrate over fitted line profile in the aim of determining exact flux emitted in this line .
so , we are unable to compute exact optical depth of the line and therefore the column density of absorbing ion .
the number of lines is to small to fit full model ( for instance xstar ) .
another lower limit on column density of cvi ion can be derived knowing that the optical thickness of the cvi lines is higher than 1 .
in such a case we can write : @xmath74 where @xmath75 is the temperature of the medium and @xmath76 is atomic mass of the ion . for the cvi@xmath58 line
we obtain : @xmath77 for the temperature approximately equal to @xmath78 k this value is twice higher than the one determined by mason et al.(2002 ) in case of mrk 766 .
this limit is over an order of magnitude weaker than the limit obtained directly from equation  [ equ : column ] .
fwhm estimated from the @xmath79 k is quite narrow and the line width is rather determined by much larger turbulent motion .
the value of reduced @xmath26 was equal 0.6363 after fitting all possible absorption lines .
after absorption lines have been identified , we were able to fit a few emission lines from an accretion disk .
we applied the standard diskline model with fixed power law index in the radial dependence of the irradiating flux @xmath80 .
we assumed @xmath81 and @xmath82 ( where @xmath83 ) .
we have found the best value of inclination angle @xmath84 for the line cv , and we adopted it ( keeping frozen ) for the other lines .
results are listed in table  [ tab2 ] , and shown in figs .  [ fig : cont ] , and [ fig8:cont ] .
the last column of table  [ tab2 ] presents theoretical values of ews of the same emission lines computed for accretion disk illuminated by hard x - ray flux equaled to the local disk flux ( full model presented in raska et al .
2002 a , b ) .
the model takes into account harder incident radiation ( @xmath85 ) , than it is seen in case of ton s180 , and solar metal abundances . for moderate illumination in the model
we have found for instance the equivalent width of nex equal 11 ev , oviii equal 18.2 , mgxii equal 10.5 , and fexvii about 8 ev , or many others at the same order of magnitude in respect to the reflected continuum ( factor of 2 - 3 smaller in respect to the total continuum ) .
models do not predict lines with ews higher than @xmath86 evs , when the disk is in the hydrostatic equilibrium ( raska et al . 2002 a , b ; see also figures in nayakshin , kazanas & kallman 2000 , ballantyne , ross & fabian 2001 ) .
higher values can be only obtained from constant density medium , as discussed by ballantyne , ross & fabian ( 2002 ) , particularly when adopting enhanced metalicity .
emission lines also become by a factor of a few stronger when the dynamical flare model is adopted , i.e. the time - scale of the flare is shorter than the time - scale for a disk surface layer to achieve a hydrostatic equilibrium ( collin et al . 2002 ) .
however , in that case lines should be narrow and not represented well by a diskline shape .
soft disk line detection is therefore a promising tool to constrain the conditions within the cold reprocessing medium .
detection of standard broad disk lines presented in this paper , if real , suggest no strong departure from hydrostatic equilibrium and the extension of the cold disk roughly down to the marginally stable orbit .
after all lines have been fitted , we fit again broken power - law with galactic absorption .
results presented in table  [ tab0 ] last verse indicate that lines do not affect continuum very strongly .
the final value of reduced @xmath26 appeared to be 0.5531 .
llcc    @xmath87 & @xmath88 & 7.42 & _ instumental _
+ @xmath89 & @xmath90 & 6.6 & _ instumental _
+ @xmath91 & @xmath92 & 22.29 & _ instumental _
+ @xmath93 & @xmath94 & 6.34 & _ instumental _
+ @xmath95 & @xmath96 & 15.35 & _ instumental _
+ @xmath97 & @xmath98 & 6.92 & _ instumental _ + @xmath99 & @xmath100 & 14 . & mgxi ( 0.2971 ; 0.3155 ) + @xmath101 & @xmath102 & 5.18 & sixiii
( 0.30878 ; 0.3279 ) + @xmath103 & @xmath104 & 7.55 & sxii ( 0.31899 ; 0.3387 ) + @xmath105 & @xmath106 & 3.47 & cv ( 0.33386 ; 0.3545 ) + @xmath107 & @xmath108 & 4.45 & cvi ( 0.3459 ; 0.3673 ) + & & & blend with mgxii ( 0.368 ) , + & & & sixiv ( 0.3703 ) , cv ( 0.3709 ) + @xmath109 & @xmath110 & 5.16 & cvi ( 0.4428 ; 0.4702 ) + @xmath111 & @xmath112 & 3.47 & fexvii ( 0.6868 ; 0.72932 ) + & & & blend with sxvi ( 0.731 ) + @xmath113 & @xmath114 & 3.29 & nex ( 1.1387 ; 1.209 ) + & & & blend with fexxv ( 1.206 ) , + & & & fexxv ( 1.208 ) , and fexxv ( 1.213 ) + @xmath115 & @xmath116 & 2.95 & mgxii ( 1.3861 ; 1.472 ) + @xmath117 & @xmath118 & 3.7 & mgxii
( 1.7327 ; 1.840 ) + & & & blend with sixiii ( 1.853 ) + @xmath119 & @xmath120 & 6.09 & sixiv ( 1.8831 ; 2.000 ) +    lllcc    @xmath121 & @xmath122 & 72.73 & cv ( 0.2899 ; 0.3079 ) & 8.69 + @xmath123 & @xmath124 & 21.26 & oviii ( 0.6149 ; 0.653 ) & 19.0 + @xmath125 & @xmath126 & 46.26 & nex ( 0.961 ; 1.020 ) & 11.29 +
in this section , we present models of warm absorber , allowing density , temperature and so consequently ionization state to be stratified ( krolik 2002 ) .
the only assumption which we made is the constant pressure .
such an approach is more realistic than the constant density model , since thermal instabilities can produce gas with steep ionization structure keeping all zones in pressure equilibrium .
we have used code titan ( dumont , abrassart & collin 2000 ) to compute temperature and density structure of the warm absorber together with radiative transfer .
the code is based on two - stream approximation and works in plane - parallel geometry .
radiative transfer is computed both in lines and in continuum so the code can be used for very inhomogeneous even optically thick media .
thermal , ionization and statistical equilibrium of ions are computed in complete non - lte for 10 elements and all corresponding ions .
we consider following elemental abundances : h : @xmath127 , he : @xmath128 , c : @xmath129 , n : @xmath130 , o : @xmath131 , ne : @xmath132 , mg : @xmath133 , si : @xmath134 , s : @xmath135 , fe : @xmath136 , .
we transfer totally 721 lines of different ions , taking their oscillator strengths and energies from national institute of standards and technology ( nist , http://physics.nist.gov ) .
photo - electric absorption is fully taken into account .    in our models , we have assumed the shape of incident radiation consistent with this fitted to the data of ton s180 .
we have adopted power law photon index @xmath137 ( sec  [ sec : cont ] ) , and spectral limits from @xmath138 up to @xmath139 kev as we have concluded from broad band spectral energy distribution presented by turner et al .
( 2002 ) .
the assumption of constant pressure refers to the total : radiation and gas pressure , and its value is self consistently determined by the code .
the ionization parameter on the top of the absorber @xmath140 and the total column density are free parameters of our models .
we have computed models for @xmath141 and @xmath142 named a1 , @xmath143 and @xmath144 named a2 , and @xmath145 and @xmath142 named a3 .
such values of column densities for given @xmath140 are upper limits above which optical depth of cloud is large enough that strong absorption edges become present in outgoing spectra .    even for high @xmath140 ,
the maximum compton temperature is low , @xmath146 k , because the incident spectrum is soft ( @xmath147 ) .    in fig .
[ fig : mod ] we show the temperature structure for a1 ( dashed line ) ( higher @xmath140 does not change this structure ) compared to a2 ( solid line ) and in fig .
[ fig : mod1 ] we show the corresponding spectra . in case of higher
@xmath140 ( model a1 ) the temperature range within the medium is much broader and the effective spectrum from such region ( presented in fig .
[ fig : mod1 ] dashed line ) displays slightly stronger lines than for a2 ( solid line spectrum ) .
we have enlarged range of frequencies to illustrate region observed in ton s180 data .
we have also performed comparison of a1 model to the warm absorber with constant density ( fig .
[ fig : mod ] dotted - dashed line ) assuming that both clouds have the same total column densities and they are illuminated by the equivalent fluxes .
thus , choosing the volume density for constant density cloud , the ionization parameter is determined by the value of flux the same as for stratified cloud .
we have chosen volume density equal @xmath148 , which implies @xmath149 .
spectrum for constant density cloud presented in fig .
[ fig : mod1 ] by dotted - dashed line displays strong edges and can not be applicable to explain our observations .
decreasing volume density we increase ionization parameter , and it may remove edges .
we plan to make proper calculations in the future work to show how to distinguish between low density constant cloud and stratified cloud .
we have also calculated theoretical values of the equivalent widths of transfered lines , the same which are found in ton s180 data .
results are presented in table  [ tab3 ] .
ews for high surface ionization ( a1 ) are almost a factor of two higher that for a2 model . in our data
we see slightly larger ews and we can not model them by increasing @xmath140 , since model a3 shows the same ews as a1 .    therefore , we have to introduce model labeled a4 in table  [ tab3 ] , which has exactly the same parameters as a2 i.e. @xmath143 and @xmath144 and additionally we assume turbulent velocity @xmath150 .
turbulent velocity does not change the total temperature structure , it only affects lines , making their ews higher , as shown in table  [ tab3 ] . for model a4 several lines , for instance mgxi , sixiii , sxii cvi , and sixiv )
have ews in quite good agreement with those achieved from observations ( see table  [ tab1 ] ) .    in all presented models
the ratio of @xmath151[cvi@xmath152 $ ] / @xmath151[cvi@xmath153 $ ] equals 0.7 , in very good agreement with observed value of 0.6 .
llllll    sviii & 2s@xmath154 2p@xmath155 @xmath154p@xmath156 - 2s@xmath154 2p@xmath155 3d @xmath154l & 0.522 & 0.293 & 0.506 & 1.431 + fexv & 3s@xmath154 @xmath157s - 3s4p @xmath157p@xmath156 & 0.406 & 0.191 & 0.392 & 1.164 + six & 2s@xmath154 2p @xmath154p@xmath156 - 2s@xmath154 3d @xmath154d & 0.584 & 0.286 & 0.568 & 1.193 + sx & 2s@xmath154 2p@xmath158 @xmath159s@xmath156 -2s@xmath154 2p@xmath154(@xmath158p ) 3l @xmath159p & 0.585 & 0.305 & 0.567 & 1.549 + cv & 1s@xmath154 @xmath157s - 2p @xmath157p@xmath156 & 1.112 & 0.587 & 1.079 & 1.849 + sxi & 2s@xmath154 2p@xmath154 @xmath158p - 2s@xmath154 2p 3d @xmath158l@xmath156 & 0.680 & 0.333 & 0.659 & 1.500 + mgxi & 2p @xmath158p@xmath156 - 4l @xmath158l & 0.834 & 0.468 & 0.809 & 1.995 + sixiii & 2s @xmath157s - 3p @xmath157p@xmath156 & 0.783 & 0.439 & 0.759 & 2.019 + sxii & 2s@xmath154 2p @xmath154p@xmath156 - 2s@xmath154 3d @xmath154d & 0.719 & 0.348 & 0.697 & 1.580 + cv & 1s@xmath154 @xmath157s - 3p @xmath157p@xmath156 & 1.287 & 0.708 & 1.249 & 2.151 + cvi & 1s @xmath154s - 2l @xmath154l & 1.186 & 0.599 & 1.136 & 1.772 + cvi & 1s @xmath154s - 5l @xmath154l & 1.650 & 0.881 & 1.585 & 2.585 + fexvii & 2s@xmath154 2p@xmath160 @xmath157s - 2s@xmath154 2p@xmath1613s @xmath157p@xmath156 & 1.216 & 0.621 & 1.159 & 4.008 + nex & 1s @xmath154s - 3l @xmath154l & 3.151 & 1.856 & 3.063 & 7.485 + mgxii & 1s @xmath154s - 2l @xmath154l & 3.561 & 2.101 & 3.561 & 9.107 + mgxii & 1s @xmath154s - 4l @xmath154l & 5.224 & 2.947 & 5.330 & 12.26 + sixiv & 1s @xmath154s - 2l @xmath154l & 4.594 & 2.426 & 4.598 & 11.56 +
in order to draw attention to the possible constraints resulting from detection of the broad disk lines we model the broad band spectrum of ton s180 assuming three basic spectral components : ( a ) an illuminated accretion disk ( b ) an optically thick comptonizing plasma ( c ) an optically thin ( non - thermal ) comptonizing plasma . for low or moderate accretion rate
we take the simplest keplerian disk model , emitting locally as a black body and parameterized fully by the black hole mass , @xmath162 , and accretion rate , @xmath163 , given in eddington units ( taking efficiency 1/12 appropriate for newtonian model ) .
we consider also highly supercritical accretion disk applying the model described by kawaguchi ( 2002 ) .
this model includes the effect of strong radial advection and the departure from the local black body emission due to electron scattering which are essential at very high accretion rates .
we assume that the energy generated in the inner disk illuminates outer disk , with the irradiating flux @xmath164 , with @xmath165 .
this irradiation effect may come from the scattering of the disk radiation back toward the disk by the warm absorber as well as from the direct interception of the radiation by the flaring outer disk .
we adopted a simple parameterization and did not model this irradiation since the process is quite complex
. it would require assumption about the radial and angular distribution of the density of the scattering medium ( see e.g. kurpiewski et al .
1997 ) and the computations of the viscosity dependent disk vertical structure , as in raska et al .
we did not compute the disk vertical structure in order to verify whether the outer disk is flaring strongly enough to allow for such a strong effect .
we assume that the inner disk for ( @xmath166 ) is fully covered by the optically thick plasma .
thermal plasma parameters are : the optical depth , @xmath167 , and the temperature , @xmath75 .
only soft photons from the inner disk come through this plasma , and we compute the spectrum using the simple code based on sunyaev & titarchuk ( 1980 ) description ( used for example by collin - souffrin et al .
we adopt the plane parallel geometry and the generation of the photons at the bottom of the zone .
this component forms a soft x - ray component of the spectrum .
we next added a contribution from much hotter , or non - thermal , plasma component in order to provide the appropriate hard x - ray tail but the parameters of this component are poorly constraint and we do not address this issue in detail .
our model provides a one - parameter family of solutions for the broad band spectrum , i.e. we can find a good model for any mass of the black hole between @xmath168 and @xmath169 , but for a fixed mass all other parameters ( @xmath163 , @xmath170 , @xmath167 , and @xmath75 ) are uniquely determined .
the appropriate accretion rate is needed to reproduce the normalization of the spectrum in the optical band , where it is mostly given by the product of @xmath171 ( for illumination case ) .
the value of @xmath170 must be chosen in order not to overproduce / underproduce far - uv flux .
the slope of the soft x - ray spectrum gives one constraint on the combination of @xmath167 and @xmath75 , but for low values of @xmath75 which are always derived the normalization of the soft comptonized spectrum also depends significantly on the plasma parameters ( for a number of soft photons already specified by @xmath170 ) .
examples of the solution are therefore as follows : for @xmath172 : we have @xmath173 , @xmath174 , @xmath175 , @xmath176 kev , and for high accretion rate case : @xmath177 : we have @xmath178 , @xmath179 , @xmath180 , @xmath181 kev ; the broad band data fit for the largest mass is shown in fig .
[ fig : bbwidmo_a ] and for the smallest mass  in fig .
[ fig : toshi ] .
those two extreme models differ considerably with the extension of the disk part covered with the hot skin . in the case of small mass
, the disk model would overpredict the far - uv data points considerably if the large portion of the inner disk is not fully covered with the thermal plasma .
plasma parameters are such that they make formation of any emission lines rather difficult .    in the case of high mass
the inner disk does not really have to be covered , as indicated by much smaller value of @xmath170 .
clumpy corona model would be easily applicable to this case , or even the x - ray emission could come from some shocks in outflowing material at the rotation axis
. the fraction of the disk covered may be always further reduced assuming that synchrotron emission of the plasma itself provides additional soft photons . in such a case the lines from the irradiated disk are easily expected .
reanalyzing _ chandra _ data of narrow line seyfert 1 galaxy - ton s180 we were able to identify several narrow absorption lines due to the warm absorber .
such lines were neither discussed by turner et al .
( 2001 ) and turner et al .
( 2002 ) where they presented the data for the first time nor seen in _ xmm - newton _ data ( vaughna et al . 2002 ) because the lines are indeed very faint .
we detected lines of highly ionized carbon and other elements , but we do not seem to detect any oxygen lines while in the fuse data ( turner et al .
2002 ) ovi narrow absorption line components were seen .
five of those absorption lines detected in ton s180 ( cvi(@xmath182 kev ) , nex , mgxii(@xmath183 kev ) , mgxii(@xmath184 kev ) , and sixiv ) are also seen in another seyfert 1 galaxy - ngc 3783 analyzed by kaspi et al .
( 2002 ) , with similar values of ews ( 0.46 , 2.3 , 4.3 , 1.6 and 6.6 ev respectively ) . also kaastra et al .
( 2002 ) in ngc 5548 have found the same cv , cvi(@xmath182 kev ) , nex , mgxii(@xmath183 kev ) , and sixiv , with following ews - 0.26 , 0.99 , 0.47 , 3.09 , and 3.88 ev .
we found a lower limit for the warm absorber column of @xmath70 @xmath6 from the equivalent with of the one of c vi lines ( see sect .  [
sec : abs ] ) and an upper limit of @xmath185 @xmath6 from our model of the warm absorber ( see sect .  [
sec : mod ] ) .
the upper limit is not firm since it is sensitive to the adopted turbulent velocity as well as abundances .
the _ chandra _ data show that the effect of the warm absorber does not modify significantly the x - ray spectrum in the case of ton s180 .
therefore our analysis do not change the overall picture of the ton s180 continuum , as shown by turner et al .
( 2002 ) and discussed by vaughan et al .
in particular , the huge soft x - ray excess is reasonably explained only within the frame of the comptonization model .
variability arguments presented by vaughan et al .
( 2002 ) in favor of hybrid plasma model ( thermal plasma being responsible for the soft x - ray slope and non - thermal plasma component leading to hard x - ray power law ) are convincing .
however , the data do not allow to determine all the hybrid plasma parameters uniquely , as explained by vaughan et al .
( 2002 ) , so in the present paper we have used much simpler model in order to draw attention to the possible constraints resulting from detection of the broad disk lines .
those emission lines can be present in the broad band spectrum rather for higher black hole mass @xmath186 and accretion rate @xmath187 , since with those assumptions hot skin does not cover disk completely .
this conclusion is opposite to the black hole mass estimations done for ton s180 by wang & lu ( 2001 ) , where @xmath188 .
but our result is consistent with studies done by siemiginowska ( 1997 ) , where simple disk corona model was fitted to the _ beppo sax _ data ( those presented by comastri et al .
we argue , contrary to the previous papers by turner et al .
( 2001 ) and vaughan et al .
( 2002 ) that spectral features in absorption and even in emission are present in the _ chandra _ data of ton s180 .
the quality of the data does not allow us to determine column densities of individual metals , we can only put lower limit on column density of cvi .
a few emission lines fitted to the ton s180 _ chandra _ data have lower equivalent widths by factor of 6 than those reported by mason et al .
( 2002 ) for seyfert i galaxy mrk 766 , but of similar order as disk lines found by kaastra et al .
( 2002 ) in ngc 5548 .
estimation of the mass of the black hole in ton s180 based on fwhm of the @xmath2 line and the formula of kaspi et al .
( 2000 ) for the relation between the size of the broad line region and the object luminosity at 5100    rather suggests small black hole mass , @xmath190 ( wang & lu 2001 ) , in contradiction with our conclusion based on the detection of the soft x - ray lines .
_ acknowledgments .
_ we are grateful to suzy collin for discussion regarding spectral models of warm absorber computed by titan .
we thank liz galle for help in processing the archival data and antonella fruscione and fabrizio nicastro for discussion .
this work was supported in part by grant 2p03d00322 of the polish state committee for scientific research and by jumelage / cnrs no .
16 `` astronomie france / pologne '' . as acknowledge support through nasa contract nas8 - 39073(cxc )    ballantyne , d.r . ,
ross , r.r .
, & fabian , a.c . , 2001 , , 327 , 10 ballantyne , d.r . ,
ross , r.r .
, & fabian , a.c . , 2002 , 336 , 867 boller , t. , brandt , w.n . , fink , h. , 1996 , , 305 , 53 collin - souffrin , s. , czerny , b. , dumont , a .- m . , zycki , p. t. , 1996 , , 314 , 393 collin s. , coup , s. , dumont , a .- m . , petrucci , p .- o . ,
raska , a. , 2002 , ( accepted to ) collinge , m.j . , brandt , w.n .
, kaspi , s. , crenshaw , d.m . , elvis , m. , kraemer , s.b . , reynolds , c.s . ,
sambruna , r.m . ,
wills , b.j . , 2001 , , 557 , 2 comastri a. , fiore , f. , guainazzi , m. , matt , g. , stirpe , g. m. , zamorani , g. , brandt , w. n. , leighly , k. m. , piro , l. , molendi , s. , parmar , a. n. , siemiginowska , a. , & puchnarewicz , e. m. , 1998 , 333 , 31 czerny , b. & dumont , a .- m .
, 1998 , , 338 , 386 czerny , b. & elvis , m. , 1987 , , 321 , 305 czerny , b. & .
zycki , p.t . , 1994 , , 431 , l5 dickey , j.m . & lockman , f.j . , 1990 ,
, 28 , 215 done , c. , mulchaey , j.s . ,
mushotzky , r.f . ,
arnaud , k.a .
, 1992 , , 395 , 275 dumont , a .- m . , abrassart a. , & collin s. 2000 357 , 823 kaastra , j.s . ,
steenbrugge , k.c . ,
raassen , a.j.j . , van der meer , r.l.j . ,
brinkman , a.c .
, liedahl , d.a . , behar , e. , de rosa , a. , 2002 , , 386 , 427 kaspi , s. , et al .
2000 , , 533 , 631 kaspi , s. , et al .
, 2002 , , 574 , 643 kawagushi , t. , 2002 , ( submitted to ) kinkhabwala , a. , sako , m. , behar , e. , kahn , s. m. , paerels , f. , brinkman , a.c . ,
kaastra , j.s .
, gu , m.f . , liedahl , d.a . , 2002 , , 575 , 732 krolik , j.h . ,
2002 proceedings of the workshop : `` x - ray spectroscopy of agn with chandra and xmm - newton '' , eds . th .
boller , s. komossa , s. kahn , h. kunieda krolik , j.h .
, kriss , g.a . , 2001 , , 561 , 684 krolik , j.h . ,
mckee , c.f .
, & tarter , c.b .
, 1981 , 249 , 422 kurpiewski a. , kuraszkiewicz j. , czerny b. , 1997 , , 285 , 725 leighly , k.m . ,
zdziarski , a.a .
, kawaguchi , t. , matsumoto , c. , 2002 proceedings of the workshop : `` x - ray spectroscopy of agn with chandra and xmm - newton '' , eds .
boller , s. komossa , s. kahn , h. kunieda mihalas , d. 1978 , stellar atmospheres , w.h .
freeman & co. , san francisco nayakshin , s. , kazanas , d. , & kalman , t. r. , 2000 , , 537 , 833 osterbrock , d.e . ,
pogge , r.w . , 1985 , , 297 , 166o puchnarewicz , e.m . , mason k.o . , siemiginowska a. , fruscione a. , comastri a. , fiore f. & cagoni i. , 2001 , , 550 , 644 romano p. , turner t.j .
, mathur s. , & george i.m . , 2002 , , 564 , 162 raska a. , czerny b. , ycki p.t . , pojmanski g. , 1999 , , 305 , 481 raska , a. , dumont , a .- m . , czerny , b. , & collin , s. , 2002a , , 332 , 799 raska , a. , czerny , b. , dumont , a .- m . ,
collin , s. , siemiginowska , a. `` the radiation spectra of illuminated accretion disks in agn '' , active galactic nuclei : from central engine to host galaxy , 2002b , eds .
: s. collin , f. combes , i. shlosman , asp conference series ( in press ) stark , a.a . , gammie , c.f . ,
wilson , r.w . ,
bally , j. , linke , r.a . , heiles , c . ,
& hurwitz , m. 1992 , , 79 , 77 sunyaev r.a .
, titarchuk l.g . , 1980 , , 86 , 121 turner , t. j. , george , i. m. , & nandra , k , 1998 , , 508 , 648 turner , t. j. , george , i. m. , yaqoob , t. , kriss , g. , crenshaw , d. m. , kraemer , s. , zheng , w. , wang , j. , & nandra , k. , 2001 , , 548 , l13 turner , t. j. , romano , p. , kraemer , s. b. , george , i. m. , yaqoob , t. , crenshaw , d. m. , storm , j. , alloin , d. , lazzaro , d. , da silva , l. , pritchard , j. d. , kriss , g. , zheng , w. , mathur , s. , wang , j. , dobbie , p. , & collins , n. r. , 2002 , , 568 , 120 vaughan , s. , boller , th . ,
fabian , a.c .
, ballantyne , d.r . , brandt , w.n . ,
tr.umper , j. , ( astro - ph/0207650 ) wang & lu 2001 , aa , 377 , 52 .
wisotzki , l. , dreizler , s. , engels , d. , fink , h .- h . , & heber , u , 1995 , , 297 , l55 yaqoob , t. 2002 , ( submitted to )


by now , jump - diffusion processes and more general lvy processes or lvy driven processes are widely used in models of the stock dynamics and term structure models , under both historic and risk neutral measure . however , only discrete time series of empirical data are available , therefore the problem of identification of the jump and diffusion components is a non - trivial task .
the problem of the identification of the diffusion component was considered in detail in ait - sahalia ( 2002 , 2003 ) ; in the present paper , the focus is on the jump component . clearly , the presence of jumps is most prominent in prices of contingent claims near expiry , especially for out - of - the money options , and therefore it is natural to extract information about the jump part of the process by using empirical data for out - of - the - money european options .
carr and wu ( 2003 ) showed that in the presence of jumps , the prices of out - of - the - money european options near expiry are of order @xmath3 , where @xmath3 is the time to expiry , and use this observation to test for the presence of jumps .
however , the results in carr and wu ( 2003 ) are qualitative rather than quantitative .
the principal aim of the present paper is to derive simple analytical formulas which can be used for the identification of the jump part of a lvy process ( under a risk - neutral measure chosen by the market ) ; the straightforward generalizations to mean - reverting processes driven by lvy processes ( the ornstein - uhlenbeck type processes ) , stochastic volatility models with jumps and lvy - driven term structure models will be published elsewhere .
the first group of results of the present paper are formulas for @xmath0 , or time decay , of out - of - the money european call and put options on a stock , at expiry .
contrary to the gaussian case , @xmath0 for the former ( resp .
, the latter ) is shown to be negative if there are positive ( resp . , negative ) jumps . in other words , we prove that near expiry , these prices are linear functions of time to expiry ( up to relatively small errors ) , and calculate the coefficients .
the results are explicit analytical expressions in terms of the parameters which characterize the process . by using these formulas and empirical data for options near expiry as in carr and wu ( 2003 )
, it is possible to identify the jump part of a particular lvy model .
after that , by using the data for options far from expiry , one can identify the first instantaneous moment of the process , and calculate the diffusion coefficient by subtracting the second instantaneous moment of the jump part of the process .
this gives a new identification method for the diffusion component . by using the formulas for the @xmath0 s for the european
put and call , we obtain approximations for the american put price and early exercise boundary near expiry .
the pricing of american options in the gaussian case is well understood
see e.g. musiela and rutkowski ( 1997 ) , karatzas and shreve ( 1998 ) , carr and faguet ( 1994 ) , carr ( 1998 ) , longstaff and schwartz ( 2001 ) , clment et al .
( 2002 ) , and the bibliography therein . in the infinite horizon case ,
explicit analytical formulas are obtained , and in the finite horizon case , a number of efficient numerical methods are developed .
although explicit formulas are not available in the finite horizon case , several general results are proved for this case as well .
consider the american put with the strike price @xmath4 and maturity @xmath5 , on a non - dividend paying stock ; the riskless rate @xmath6 is constant .
let @xmath7 be the optimal exercise price of the american put .
if the stock log - price @xmath8 follows the brownian motion , then the asymptotic behavior of @xmath9 is known both near expiry : @xmath10 as @xmath11 ( see , e.g. , barles _
( 1995 ) and lamberton ( 1995 ) ; an asymptotic formula for the put price is also available ) , and for small values of @xmath12 : for any @xmath13 , @xmath14 ( see , e.g. musiela and rutkowski ( 1997 ) ; ( [ asymprgauss ] ) is formulated as follows : if the riskless interest rate is zero , then the early exercise of the american put is not optimal ) .
notice that ( [ asbgauss ] ) implies that @xmath15 and that ( [ exlim ] ) and ( [ asymprgauss ] ) do not agree well when both the riskless interest rate and time to expiry vanish .
in levendorski ( 2003 ) , a variant of carr s randomization method was developed for wide families of lvy processes , and as a by - product , it was shown that for many families of non - gaussian processes used in empirical studies of financial markets , the analogs of ( [ exlim ] ) and ( [ asymprgauss ] ) agree much better .
namely , @xmath16 where @xmath17 depends on some of parameters which characterize the price process , and is _ smaller _ than the strike price , @xmath4 .
moreover , it was proved that @xmath17 vanishes with @xmath12 , and therefore , @xmath18 notice that ( [ asymprlevy ] ) looks more natural than ( [ asymprgauss ] ) , and may be of some importance in today s almost deflationary world .
further , if the difference @xmath19 is sizable , then any numerical procedure or parameters fitting procedure which does not take this fact into account , may produce significant errors near expiry and strike . hence , the non - standard behavior of the early exercise boundary has practical implications for the parameters fitting and development of accurate numerical methods .    however , the proof in levendorski ( 2003 ) has three disadvantages , which stem from the use of the wiener - hopf method .
first , the non - standard behavior of the early exercise boundary results from long calculations without clear intuition .
next , the approximate formulas for the put price are expressed in terms of the factors in the wiener - hopf factorization formula , hence , they are complicated for practical use , and finally , the method in levendorski ( 2003 ) is not directly applicable to a popular class of variance gamma processes , which was used in carr et al ( 2002 ) , madan and hirsa ( 2003 ) , and carr and hirsa ( 2003 ) for pricing of american put .
the generalization of the method for more general classes of processes seems impossible .
here we prove ( [ exlimlevy ] ) and ( [ asymprlevy ] ) by using a much simpler argument .
we obtain an efficient approximation for the put price near expiry , and formulate the result about the non - standard behavior of the early exercise boundary in a more meaningful form . also , we extend the result in levendorski ( 2003 ) to variance gamma processes , and in addition , we find a condition on a process which ensures that the behavior of the boundary is even more regular than in the gaussian case : there exists a negative constant @xmath20 s.t . @xmath21
we prove that the limit of the boundary at expiry is located where @xmath22 of the european call is equal to the riskless rate .
( if the @xmath22 of the european call at expiry is less than the riskless rate for all spot prices less than the strike , then the limit is the strike . ) when the corresponding equation is solved , one can use the formulas for @xmath0 s of the european call and put and obtain an efficient approximation to the put price near expiry .
the formulas are fairly simple and can be used for parameters fitting purposes near expiry , both for the european options and american put .
they also can be used to improve performance of numerical methods : when using the backward induction , one can start not at expiry , where the value function is most irregular but close to it , and use the analytical approximation as the basis for the backward induction .
we demonstrate how large the deviations from the standard behavior of the early exercise boundary and put price can be by using the empirical data for the american put documented in carr et al .
we found that the non - standard behavior should be observed for reasonable values of the riskless rate : typically , the margin between the strike and boundary is 10% or more , and for the spot price in the 5%interval around the strike , the  jump premium " in the put price over the pay - off is @xmath23 or more , where @xmath3 is the time to expiry ( with the normalization @xmath24 for a business year ) .
the last aim of the paper is to describe the restrictions on the log - price process under an emm , which should be satisfied if the non - standard behavior of the early exercise boundary is observed . as ( [ exlim ] ) , ( [ exlimlevy ] ) and ( [ superregular ] ) imply , various types of the behavior of @xmath9 near expiry are possible
however , for many families of processes used in empirical studies of financial markets , namely , hyperbolic processes , normal inverse gaussian processes ( nig ) , and truncated lvy processes of the extended koponen s family ( a.k.a .
kobol processes , a.k.a .
cgmy model ) of order @xmath25 , the behavior of the early exercise boundary is non - standard for all parameters values .
this means that if in the empirical studies , the early exercise boundary for put options on a particular stock exhibits the standard behavior near expiry , then the use of any of the processes above to model the log - price dynamics under a risk - neutral measure chosen by the market is inconsistent with the data .
normal inverse gaussian processes , hyperbolic processes and kobol processes of order greater than 1 are processes with sample paths of infinite variation .
notice that the non - standard behavior of the early exercise boundary is proved for much wider classes of infinite variation lvy processes as well . for processes with finite variation jump component , the prominent examples being variance gamma processes and kobol processes of order less than 1 , both the non - standard behavior ( [ exlimlevy ] ) and  super - regular " behavior ( [ superregular ] ) of the early exercise boundary are possible , depending on parameters of the process
; there are also borderline cases , when we are able only to prove that ( [ exlim ] ) holds . according to a recent empirical study in carr et al .
( 2002 ) , finite variation processes fit empirical data better .    even for a finite variation process , if we fix all the parameters of a process except for the drift @xmath26 , let @xmath27 and change @xmath26 accordingly so that the stock remains a martingale , then below a certain critical value of @xmath12 ( depending on the parameters of the process ) ,
the non - standard behavior is observed , and as @xmath12 decreases further , the optimal exercise price goes to zero uniformly over the interval @xmath2 .
it is natural to ask how robust the result about the non - standard behavior of the early exercise boundary is . in levendorski ( 2003 ) and
the present paper , it is assumed that dynamics of log - prices @xmath8 follows a non - gaussian lvy process under a risk - neutral measure chosen by the market .
this simplest choice of a non - gaussian risk - neutral measure for pricing of the american put was employed in carr et al ( 2002 ) , madan and hirsa ( 2003 ) , carr and hirsa ( 2003 ) , where parameters fitting and computational issues where addressed .
better fit can be achieved with stochastic volatility models with jumps but the method in levendorski ( 2003 ) can not be generalized for stochastic volatility models .
the method of the present paper can be applied to stochastic volatility models ( and term structure models ) , and under certain parameter restrictions , the margin between the early exercise boundary and strike can be derived .
further , in levendorski ( 2003 ) and the present paper , the formulas for the margin are derived in two different approximate models , and in cases which were considered in the both papers , the formulas for the margin are the same .
the identical results in two approximations allows one to hope that the margin in the continuous time model is the same . as a proven theorem
, we can only claim that the margin in the approximate model is not less than in the continuous time model .
the rest of the paper is organized as follows . in section [ levy ] , we recall the definition of lvy processes , and give examples of several families of lvy processes used in empirical studies of financial markets . in section
[ europe ] , @xmath0 s of out - of - the - money european put and call options are calculated , and in section [ amerput ] , the approximate formulas for the american put price and early exercise boundary near expiry are derived . in section
[ implications ] , we consider several families of lvy processes used ( and suggested for use ) in empirical studies of financial markets , and discuss the implications of non - standard behavior of prices and early exercise boundary for the choice of a family of processes , and parameters fitting in more detail . in section [ empirics ] , we derive properties of the early exercise boundary for kobol processes , for parameters values documented in carr et al .
( 2002 ) , and show that the margin between the strike price and early exercise boundary can be significant indeed : more than 10% . in the appendix , the technical proofs are given .
recall that a lvy process is a process with stationary independent increments ( for general definitions , see e.g. sato ( 1999 ) ) .
a lvy process may have a gaussian component and/or pure jump component .
the latter is characterized by the density of jumps , which is called the lvy density .
we denote it by @xmath28 .
also , a lvy process can be completely specified by its characteristic exponent , @xmath29 , definable from the equality @xmath30=e^{-t\psi(\xi)}$ ] ( we confine ourselves to the one - dimensional case ) .
if @xmath31 is a lvy process with finite variation jump component , then the characteristic exponent is given by @xmath32 where @xmath33 and @xmath26 are the variance and drift coefficient of the gaussian component , and @xmath34 satisfies @xmath35 equation ( [ psi1 ] ) is a special case of the lvy - khintchine formula ; for the general case , see e.g. sato ( 1999 ) .
wide families of jump - diffusion processes used in the theoretical and empirical studies of financial markets are lvy processes with finite variation jump component .
0.1cm_example 2.1 .
_ let @xmath36 be a lvy process with the lvy density @xmath37 where @xmath38 , and @xmath39 .
then @xmath40 where @xmath41 and @xmath42 are the variance and drift of the gaussian component .
the @xmath43 is holomorphic in the strip @xmath44 . in order that @xmath45 $ ] be finite , we need to impose an additional condition @xmath46 . for the sake of brevity
, we will consider below the exponential jump - diffusions with the characteristic exponent of the form ( [ jumpdif ] ) .
the generalization to the case of lvy densities given by general exponential polynomials on positive and negative half - axis is straightforward .
0.1cm_example 2.2 . _
variance gamma processes ( vgp ) have been developed and used by madan and co - authors in a series of papers during 90th ( see madan et al .
( 1998 ) and the bibliography there ) .
the characteristic exponent of a vgp without the diffusion component can be represented in the form @xmath47,\ ] ] where @xmath48 and @xmath49 .
the condition @xmath45<+\infty$ ] imposes an additional restriction @xmath46 .
0.1cm_example 2.3 .
_ truncated lvy processes ( tlp ) constructed by koponen ( 1995 ) were used for modeling in real financial markets in bouchaud and potters ( 1997 ) , cont et al .
( 1997 ) and matacz ( 2000 ) ; a generalization of this family was constructed in boyarchenko and levendorski ( 1999 , 2000 ) and called the extended koponen family of tlp processes .
later , this generalization was used in carr et al .
( 2002 ) under the name cgmy - model . as a.n .
shiryaev remarked , the name tlp was misleading , and so starting with boyarchenko and levendorski ( 2002a-2002c ) , we use the name kobol processes .
the characteristic exponent of a kobol process is of the form @xmath50,\ ] ] where @xmath51 , and @xmath42 .
the condition @xmath45<+\infty$ ] imposes an additional restriction @xmath52 .
( in boyarchenko and levendorski ( 1999 , 2000 , 2002b ) , the reader can also find a formula for the kobol process of order @xmath53 ) .    for @xmath54 , the equation ( [ kbl1d ] )
is obtained from ( [ psi1 ] ) with the lvy density given by @xmath55 and @xmath56 ( that is , there is no gaussian component ) . in the case
@xmath57 , instead of ( [ psi1 ] ) , the general form of the lvy - khintchine formula is needed ( see sato ( 1999 ) and boyarchenko and levendorski ( 2002b ) ) .
kobol processes of order @xmath58 are finite - variation processes , and the ones of order @xmath59 are infinite - variation processes
.    other examples of infinite variation processes are hyperbolic processes ( hp ) , normal inverse gaussian processes ( nig ) and normal tempered stable ( nts ) lvy processes of order @xmath60 ; in the general classification of boyarchenko and levendorski ( 1999 , 2000 , 2002a-2002c ) , hp and nig are processes of order 1 .
hyperbolic processes were constructed and used by eberlein and co - authors ( see e.g. eberlein and keller ( 1995 ) , eberlein et al .
( 1998 ) , eberlein and raible ( 1999 ) ) ; hyperbolic distributions were constructed in barndorff - nielsen ( 1977 ) .
normal inverse gaussian processes ( nig ) were introduced in barndorff - nielsen ( 1998 ) and used to model german stocks in barndorff - nielsen and jiang ( 1998 ) ; generalization of the class nig , namely , the class of normal tempered stable ( nts ) lvy processes , was constructed in barndorff - nielsen and levendorski ( 2001 ) and barndorff - nielsen and shephard ( 2001 ) . as kobol processes , ntslp can be of any order between 0 and 2 .    0.1 cm _ example 2.4 . _
the characteristic exponent of an nts lvy process is of the form @xmath61,\ ] ] where @xmath62 @xmath63 , and @xmath64 .
the condition @xmath45<+\infty$ ] imposes an additional restriction @xmath65 . with @xmath53
, we obtain the characteristic exponent of nig processes .
0.1cm_remark 2.1 .
_ for the explicit calculations in the next two sections , it is important that the characteristic exponent of any of lvy processes considered above admit the analytic continuation into the complex plane with cuts @xmath66 , [ i{{\lambda_+ } } , + i\infty)$ ] , where @xmath52 and @xmath67 , and the analytic continuation is defined by the same formula ( in example 2.4 , @xmath68 and @xmath69 , and in example 2.1 , the characteristic exponent is analytic in the complex plane with two poles at @xmath70 and @xmath71 ) .
this property is unnecessary if we want to obtain only qualitative results .
for instance , the existence of a non - vanishing margin between the strike and early exercise boundary will be proved for more general class of lvy processes with the jump part of infinite variation .
to simplify the presentation , we normalize the strike price to 1 .
let @xmath31 be a lvy process of any class considered above .
let @xmath72 and @xmath73 be the price of the european call and put option , respectively , at time @xmath74 , and the spot price @xmath75 .
denote by @xmath76 , and @xmath77 , the opposite to the @xmath78 of the out - of - the - money european call and put options , respectively , at expiry ( the option s @xmath0 is the derivative of the option price with respect to time ) : for @xmath79 , @xmath80 and for @xmath81 , @xmath82 if the process is gaussian , then the prices of out - of - the - money european call , @xmath83 , and put , @xmath84 , vanish faster than any power of @xmath3 , as @xmath85 , hence the limits ( [ cc ] ) and ( [ cp ] ) equal 0 . the next theorem and lemma
show that in the case of non - gaussian lvy processes which we consider , both @xmath86 and @xmath77 , exist and do not vanish .
the intuition is that when the jumps are present , the value of waiting of a positive ( respectively , negative ) movement in the price is larger than in the gaussian case .
[ margin1 ] let @xmath31 be any of lvy processes considered above .
then 0.cm a ) for @xmath79 , the limit ( [ cc ] ) exists , and it is positive ; 0.cmb ) for @xmath81 , the limit ( [ cp ] ) exists , and it is positive ; 0.cmc ) the time decay at expiry , of out - of - the money european call , @xmath87 @xmath88 and put , @xmath77 , depends only on the positive and negative jump parts of the process , respectively , but not on the drift and gaussian component .
notice that the intuition behind the independence of the time decay at expiry of the gaussian component is simple : the movements caused by the latter are of order @xmath89 , on average , and therefore the main contribution to the price of an out - of - the - money option , ( almost ) at expiry , can come from the jumps in the corresponding direction only ( for a rigorous proof , see the appendix ) .
we start with the simplest case of a pure jump process ( example 2.1 with @xmath56 ) .
let @xmath90 be the infinitesimal generator of @xmath36 .
then the price @xmath91 of a contingent claim with the terminal pay - off @xmath92 can be represented in the form @xmath93g(x+\mu \tau).\ ] ] for the pure jump process , @xmath90 acts as follows : @xmath94 and since @xmath95 for processes in example 1.1 , we conclude that @xmath90 is a bounded operator in @xmath96 . hence , if @xmath97 is bounded ( the case of the put option : @xmath98 ) , we can write ( [ europclaim ] ) in the form @xmath99 if the put option is out - of - the - money , that is , @xmath81 , then for small @xmath3 , @xmath100 , hence the zero - order term in ( [ europclaim2 ] ) is zero , and @xmath101 by dividing by @xmath102 and passing to the limit , we obtain @xmath103 it remains to calculate @xmath104 at @xmath81 .
we have @xmath105 and by calculating the integral , we obtain @xmath106 the case of the call option requires essentially the same calculations but @xmath90 should be regarded as a bounded operator in a space of continuous functions , which grow as @xmath107 as @xmath108 ( for details , see boyarchenko and levendorski ( 2002b ) ) , and the result is @xmath109 similarly , @xmath110 and @xmath111 for more general exponential jump - diffusions can be calculated
. moreover , the proof above gives @xmath112 and @xmath113 provided @xmath114 and @xmath115 ( the latter condition is necessary for the stock to be priced ) .
we conjecture that ( [ ccgen ] ) and ( [ cpgen ] ) are valid for any lvy process satisfying ( [ denspos1 ] ) .
however , we were unable to prove ( [ cc])([cp ] ) in the full generality , and so in the proof of theorem [ margin1 ] for the other families of lvy processes listed above , we calculate the limits @xmath110 and @xmath111 directly , by using the formulas for the characteristic exponents . in the case of nts
lvy processes and nig , set @xmath116 .
now , for any @xmath31 of kobol , nts , and vgp classes , and for @xmath117 and @xmath118 , define @xmath119.\ ] ]    [ gen ] let @xmath31 be a vgp , kobol , or nts lvy process .
then for @xmath79 , @xmath120 and for @xmath81 , @xmath121    in the appendix .    in the gaussian case , @xmath122 , and since @xmath123 depends linearly on @xmath29 , we conclude that @xmath123 depends on the jump component of the process only .
hence , @xmath76 , and @xmath77 , depend on the jump component only .
the fact that the former depends on the positive jumps only , and the latter on the negative ones , follows by inspection of the lvy - khintchine formula : for @xmath118 ( resp . , @xmath117 ) , @xmath124 is independent of the density of negative ( resp .
, positive ) jumps .    in the appendix
, we will calculate @xmath124 for each family of processes , and derive from ( [ genc])([genp ] ) the following formulas :    \1 ) for kobol processes of order @xmath125 : @xmath126 the integrals converge since @xmath127 and @xmath128 ;    \2 ) for vgp : @xmath129 the integrals converge since @xmath130 ;    \3 ) for nts lvy processes @xmath131^{\nu/2}}{z(z+1)}dz,\quad x<0;\\ \label{cpnts } { { \mathcal p}}(x)&=&\frac{{\delta}}{\pi}\sin\frac{\pi\nu}{2}\int^{+\infty}_{{\alpha}+{\beta } } \frac{e^{-xz}[(z-{\beta})^2-{\alpha}^2]^{\nu/2}}{z(z+1)}dz,\quad x>0;\end{aligned}\ ] ] with @xmath53 , formulas for nig processes obtain .
the integrals converge since @xmath132 , @xmath133 , and @xmath128 .    by direct inspection of formulas ( [ cckbl])([cpnts ] ) , we conclude that with @xmath134 , the integrands decay as @xmath135 as @xmath136 for processes of order @xmath137 ; if @xmath31 is a vgp or exponential jump - diffusion , then the integrands decay faster than @xmath138 , for any @xmath139 . therefore , with @xmath134 , the integrals diverge if and only if the jump component of @xmath31 is a process of order @xmath140 .
further , the integrands are positive monotone continuous functions of @xmath141 . by using the monotone convergence theorem
, we deduce    [ margin2 ] a ) for processes of order @xmath58 , vgp , and exponential jump - diffusion processes in example 2.1 , there exist finite limits
@xmath142    \b ) for nig , and kobol and nts lvy processes of order @xmath143 , the limits @xmath144 and @xmath145 exist but are infinite .
\c ) @xmath146 is an increasing continuous function , which maps @xmath147 onto @xmath148 , and @xmath149 is a decreasing continuous function , which maps @xmath150 onto @xmath151 .
_ remark 3.1 . _
a ) it is possible to derive formulas for @xmath86 and @xmath77 , for the case of hyperbolic processes , and that part b ) of theorem [ margin2 ] holds in this case .
b ) can be proved for more general classes of lvy processes .
for instance , @xmath152 can be proved if in a neighborhood of 0 , the density of positive jumps , @xmath153 , admits a lower bound via @xmath154 , where @xmath155 .
then there exists a density @xmath156 such that @xmath157 and @xmath158 where @xmath159
. clearly , @xmath110 can not increase if we replace @xmath153 with @xmath160 .
the argument in the proof of theorem [ margin1 ] for the jump - diffusion case shows that the contribution of @xmath156 to @xmath110 is bounded uniformly in @xmath79 , and the proof for kobol processes shows that @xmath161 gives a contribution which is unbounded as @xmath162 .
0.1 cm    to finish this section , we list approximate formulas for european call and put options near expiry , which follow from theorem [ margin1 ] and the put - call pairity . of course , these approximations are too inaccurate in a very small neighborhood of the strike if the process is of order @xmath163 .    for @xmath81 :
@xmath164 for @xmath79 : @xmath165
assume that @xmath144 and @xmath145 exist .
they are measures of the influence of positive and negative jumps , respectively , and the emm condition is naturally formulated in terms of @xmath166 and @xmath144 and @xmath145 . as it was shown above , the limits ( finite or infinite ) exist for the families of processes in examples 2.12.4 .
[ emm ] let @xmath31 be a lvy process with finite variation jump component .
then @xmath167 is a local martingale if and only if @xmath168    since the gaussian component @xmath169 and jump component @xmath170 of a lvy process are independent , @xmath167 is a local martingale if and only if for any @xmath171 and @xmath172 , @xmath173e\left[e^{y_t}\right]\\   & = & e^{(\mu - r+{\sigma}^2/2)t+x_s}e\left[e^{y_t}\right ] ,   \end{aligned}\ ] ] whereupon @xmath174,\ ] ] and @xmath175 - 1,\ ] ] which is @xmath176 ( here @xmath177 denotes the transition kernel of the semigroup @xmath178 , @xmath179=\int_{-\infty}^{+\infty}p_t(z)f(y+z)dz.\right)\ ] ] divide by @xmath180 : @xmath181 and pass to the limit @xmath182 ; the result is ( [ emm1 ] ) .
we rewrite ( [ emm1 ] ) as @xmath183 for processes without the gaussian component , ( [ emm2 ] ) is simpler : @xmath184 the first theorem below shows that if negative jumps are present but their influence is not very large , and compensated by a positive drift ( which clearly decreases the value of waiting ) , then the early exercise boundary converges to the strike price .
moreover , the convergence is more regular than in the gaussian case in the sense that ( [ superregular ] ) holds .
let @xmath185 be the early exercise boundary .
the continuous time model can be approximated by a discrete time model with small time steps @xmath186 . in the discrete time approximation ,
the option can be exercised at @xmath187 only .
denote by @xmath188 the optimal exercise boundary at the last moment before the expiry , and set @xmath189 .
one should expect that the behavior of @xmath190 for small @xmath191 is a good proxy for the behavior of the optimal exercise boundary near expiry .
in particular , @xmath192 if the limit in the lhs exists .
the proof of a simplified version @xmath193 is straightforward : as we make the class of admissible stopping times smaller , the price of the american put does not increase , therefore the early exercise boundary does not go down .
we were unable to prove the equality ( [ just ] ) but the inequality ( [ just2 ] ) suffices to make the conclusion that if the early exercise boundary in the discrete time approximation exhibits the non - standard behavior in the limit @xmath194 , then the behavior of the early exercise boundary in the continuous time model is also non - standard .    [ extrastan ] let @xmath31 be a vgp or kobol process of finite variation , without gaussian component , or an exponential jump process .
let @xmath195 then @xmath196    in other words , _ if @xmath22 of the out - of - the - money european put , at the expiry and strike , is less than the drift , then the slope of the boundary near expiry equals the @xmath22 .
_    at @xmath197 , a put owner must be indifferent between exercising and holding the put , therefore @xmath190 is the solution to the equation @xmath198 ( the lhs is the value of the put if it is exercised , and the rhs is the expected present value of the put if it is kept alive ) , equivalently , @xmath199 assume that for small @xmath200 , @xmath201 , then @xmath202 , hence @xmath203 . by using the taylor formula and dividing by @xmath191
, we obtain @xmath204 then we pass to the limit @xmath205 @xmath206 condition @xmath207 is necessary lest the assumption @xmath208 lead to a contradiction ; and the condition ( [ cneg ] ) suffices for the inequality @xmath208 to hold for small @xmath191 .
now assume that ( [ cneg ] ) fails .
this happens when either the drift is negative or the influence of negative jumps is too large .
hence , the value of waiting should increase , and the optimal exercise boundary should be lower .
as we will see , apart from the borderline case @xmath209 , the value of waiting becomes so large that @xmath190 remains bounded away from 0 , as @xmath205 .
[ thm : finvar ] let @xmath31 be a vgp or kobol process of finite variation , without gaussian component , or an exponential jump process .
then a ) if @xmath209 , then the limit @xmath210 exists , and it is equal to 0 ;    \b ) if @xmath211 then the limit ( [ barh ] ) exists .
it is negative , and can be found from the equation @xmath212 which has a unique solution .    in other words , _ if @xmath22 of the out - of - the - money european put , at the expiry and strike , is bigger than the drift , then the limit of the early exercise boundary at expiry is less than the strike .
it is located where @xmath22 of the out - of - the - money call , at expiry , is equal to the riskless rate .
_    we start with part b ) .
first of all , we notice that due to ( [ emm3 ] ) , ( [ cneg2 ] ) is equivalent to @xmath213 therefore on the strength of theorem [ margin2 ] , c ) , the solution @xmath214 to ( [ formpos1 ] ) exists and it is unique .
secondly , the reader may be surprised that the answer is formulated in terms of the density of positive jumps , but the proof shows that this is a consequence of the put - call pairity ( for european options ) .
since @xmath215 is a ( local ) martingale under the risk - neutral measure , we have @xmath216 hence , we can rewrite the rhs in ( [ exb1 ] ) as @xmath217 ( essentially , this is the put - call pairity ) , and therefore , ( [ exb1 ] ) can be written as @xmath218 we see that if @xmath219 is negative for small @xmath191 , then @xmath220 depends only on the upper tail of the probability density . by dividing by @xmath191 and passing to the limit , we obtain ( [ formpos1 ] ) .
thus , part b ) has been proved .    in part
a ) , @xmath209 , hence from ( [ emm1 ] ) , @xmath221 , and if we assume that for some sequence @xmath222 , the sequence @xmath223 remains bounded away from 0 , then by dividing by @xmath224 in ( [ exb2 ] ) and passing to the limit , we obtain an inequality , contradiction .    in the following theorem , we allow for a gaussian component .
[ gauss ] let @xmath31 be a vgp or kobol of order @xmath58 , or exponential jump - diffusion process , with a gaussian component .
then a ) if @xmath225 , then the limit @xmath214 exists .
it is negative , and can be found from ( [ formpos1 ] ) ;    \b )
if @xmath226 , then the limit @xmath214 exists , and it is equal to 0 .    by repeating the proof of theorem [ thm : finvar ] , we obtain a ) , and similar arguments show that under condition @xmath226 , an assumption that the behavior of the boundary is non - standard , leads to contradiction . this gives b ) .
we see that in the case @xmath226 , the result is weaker than in the pure jump case , where we can single out the case when the boundary is more regular than in the gaussian case .
according to part b ) of theorem [ margin2 ] , now @xmath227 as @xmath228 . by using this observation , and repeating the proof of part a ) of theorem [ thm : finvar ]
, we obtain    [ thm : infvar ] let @xmath31 be a nig or kobol or ntsl process of infinite variation . then the limit @xmath229 exists .
it is negative , and can be found from the equation ( [ formpos1 ] ) .    _
remark 4.1 .
_ if the density of positive jumps satisfies the conditions in remark 3.1 b ) , then @xmath152 , and repeating the proof of theorem [ thm : finvar ] , we obtain that @xmath214 is negative .
let @xmath31 be a process from any of families in examples 2.12.4 .
let @xmath230 but the lvy density and the diffusion coefficient of the gaussian component remain fixed .
then the drift of the process changes with @xmath12 and converges to a finite value . by the direct inspection of ( [ formpos1 ] )
, we conclude that if the early exercise boundary exhibits the non - gaussian behavior ( [ exlimlevy ] ) near expiry , then the early exercise price near expiry vanishes with @xmath12 as well .
the same argument as in the gaussian case shows that for any @xmath102 , @xmath231 as @xmath230 ( see ( [ asymprgauss ] ) ) , therefore we conclude that the optimal exercise price tends to zero uniformly on the interval @xmath2 .    notice that if @xmath31 is a vgp or a kobol process of order @xmath232 , or exponential jump - diffusion ( as in example 2.1 ) , then for large @xmath12 , ( [ formpos1 ] ) has no solution , but for sufficiently small @xmath12 , equation ( [ formpos1 ] ) has a unique solution , and this solution tends to 0 together with @xmath12 .
this means that for sufficiently low levels of the riskless rate , all the processes considered above should lead to the non - standard behavior of the early exercise boundary .      by using the approximate formulas for the european
put and call near expiry at the end of section 3 , we obtain approximate formulas for the american put price near expiry : @xmath233 notice that this approximation is too inaccurate very close to the strike , especially for infinite variation processes .
the behavior of the early exercise boundary near expiry can be used to determine the type of a lvy process , which can be used to model the log - price dynamics under the risk - neutral measure .
this choice of a risk - neutral measure was used in carr et al .
( 2002 ) , hirsa and madan ( 2003 ) and carr and hirsa ( 2003 ) .
empirical studies of financial markets clearly indicate the presence of both downward and upward jumps ( under historic measure ) , therefore if the boundary exhibits the standard behavior , then one may not use simultaneously a process of order @xmath234 under the historic measure and lvy process under a risk - neutral measure ( chosen by the market ) . indeed , as it is shown in carr et al .
( 2002 ) , p.312 , the difference between the historic and risk - neutral lvy densities must be of finite variation , therefore a process of infinite ( finite ) variation under the historic measure must remain a process of infinite ( finite ) variation under an emm .
processes of order @xmath235 can be used for modelling of log - price dynamics both under historic and risk - neutral measure but additional restrictions on the density of the upward jumps must be imposed .
another possibility is to use one - sided kobol processes with densities of positive jumps of an order different from the one for positive jumps ( see boyarchenko and levendorski ( 1999 , 2000 , 2002b ) ) .
the characteristic exponents of these strongly asymmetric kobol processes are of the form @xmath236\\ & & \hskip1 cm + c_-\gamma(-\nu_-)[(-{{\lambda_-}})^{\nu_-}-(-{{\lambda_-}}-i\xi)^{\nu_-}],\end{aligned}\ ] ] where @xmath237 , @xmath238 , and @xmath239 ( for the formulas in the case @xmath240 , see boyarchenko and levendorski ( 1999 , 2000 , 2002b ) ) .
the analogs of ( [ cckbl ] ) and ( [ cpkbl ] ) are @xmath241 if a finite variation process is chosen for modelling of the american put on a particular stock , then the choice of parameters of the process is constrained by the observed behavior of the early exercise boundary :    a.   if the non - standard behavior ( [ exlimlevy ] ) is observed , then the constraint is @xmath242 ; b.   if the  super - regular " behavior ( [ superregular ] ) is observed , then the constraint is @xmath243 ( and probably , there is no gaussian component ) ; however , as empirical examples in the next section show , this possibility is unlikely to realize ; c.   if @xmath244 but the  super - regular " behavior ( [ superregular ] ) is not documented , then the constraint is @xmath245 , and the gaussian component is admissible .
still another possibility is to assume that under an emm chosen by the market , the rate of decay of the density of positive jumps is very large ; then @xmath244 for processes of order @xmath235 , and @xmath214 is negative but very close to 0 for processes of order @xmath246 . for kobol
, this means that the steepness parameter @xmath247 must be large in modulus .
if under the historic measure , the @xmath247 is not large , this interpretation presumes that the agents in the market are very risk averse .
the formulas for @xmath0 of out - of - the money european call and put options , at expiry , can also be used to infer the type of an appropriate family of processes , and parameters of the process . in particular ,
if @xmath0 at expiry quickly grows as the spot price approaches the strike price , then processes with infinite variation jump component must be used , otherwise the processes with finite variation jump component .
the aim of this section is to demonstrate how large the margin between the strike and early exercise boundary and the  jump premium " in the american option price over the payoff @xmath248 can be for realistic values of parameters .
we use risk - neutral parameters values of the kobol processes with a diffusion component , which were obtained for several stocks in carr et al .
( 2002 ) ( table 3 on p.327 in op .
cit . ) .    0.1 cm _ example 5.1 .
_ only in one case : ibm1111 , the process is of order @xmath249 : the risk - neutral parameters are @xmath250 , and @xmath251 ( in carr et al .
, these parameters are denoted @xmath252 , and @xmath253 ) . since @xmath249 , the non - standard behavior of the early exercise boundary
is guaranteed but the gap between the strike price and boundary is small unless the interest rate is very small indeed .
it is seen from figure 1 , where we plot the graph of @xmath214 as a function of @xmath12 , that for @xmath254 , @xmath214 is about @xmath255 , and for @xmath256 , @xmath214 is close to 0.01 , so the gap is clearly seen . in figure 2 , we plot the graph of @xmath0 at expiry , for out - of - the - money european call and put ( graph of @xmath110 , for small @xmath79 near zero , and put , @xmath111 , for @xmath81 ) .
both are unbounded in a neighborhood of @xmath257 .
since the density of positive jumps decays very fast ( the steepness parameter @xmath258 is very large ) , the time decay of out - of - the - money puts is very small unless the spot price is very close to the strike .
the density of positive jumps decays faster than the one of negative jumps ( @xmath259 ) , therefore the @xmath0 for the out - of - the - money call is smaller than the one for out - of - the - money put , for the same distance from the strike .
0.1 cm _ example 5.2 . _ for amzn1014 , the parameters are @xmath260 the fairly large @xmath261 indicates that the intensity of jumps is large , and relatively small values of @xmath262 and @xmath247 imply that there are many large jumps , especially negative ones ( recall that @xmath262 characterizes the rate of decay of the density of negative jumps ) .
therefore , even not very close to expiry , the owner of an out - of - the - money european option can expect that with a non - negligible probability the option will become in - the - money .
hence , the price of these out - of - the - money options is sizable , especially the one of the put .
figure 3 demonstrates this effect .
we plot graphs of @xmath146 and @xmath149 in a neighborhood of zero .
further , @xmath263 , therefore for reasonable @xmath12 , we have @xmath242 , and the behavior of the early exercise boundary is non - standard . since the intensity of large negative jumps is significant , the value of waiting is very large , and the gap between the strike and optimal exercise boundary is also large . from ( [ amputapprox ] ) and
figure 3 , we can conclude that even at 5% below the strike , the option value of waiting ( the difference between the american put price and the pay - off @xmath248 ) is of order @xmath23 ; @xmath24 is one year .    in the next two figures 4 and 5 , we plot the graph of @xmath146 farther from the strike , and the graph of @xmath214 as the function of @xmath12 .    in figure 6 , we plot both the early exercise boundary , at expiry , @xmath264 , as a function of @xmath12 , and that of the early exercise boundary far from expiry , @xmath265 ( the early exercise boundary in the infinite time horizon case ) .
similar results are valid for the other stocks documented in carr et al .
for several families of lvy processes used in empirical studies of financial markets , we have proved that the time decay of out - of - the - money european options at expiry is not negligible , and derived explicit formulas .
these formulas can be used to make numerical calculations near expiry more accurate  both for european options , and american ones . by using the same formulas
, we proved that the optimal exercise boundary for the american put is separated from the strike price by a non - vanishing margin on the interval @xmath2 , and that as the riskless rate vanishes , the optimal exercise price goes to zero uniformly over the interval @xmath2 .
the last result may be especially interesting since the current levels of interest rates are rather low .
furthermore , we discussed the restrictions on the process under the risk - neutral measure which the non - gaussian behavior of the early exercise price near expiry imposes , and suggested a natural asymmetric extension of the kobol family ( a.k.a . cgmy model ) , which can be used for any type of the behavior of the boundary .
we showed that for the risk - neutral parameters obtained for several stocks in carr et al .
( 2002 ) , the behavior should be non - standard , and the margin between the strike and the early exercise boundary be quite sizable : more than 10% , up to expiry .
the jump premium over the payoff is shown to be sizable as well .
further , we assume that @xmath46 ; the results for the case @xmath266 obtain by passing to the limit @xmath267 in the formulas for the case @xmath46 . for any @xmath268 ,
the price of the european call option is given by @xmath269 } { -i\xi(1-i\xi)}d\xi\ ] ] ( see e.g. boyarchenko and levendorski ( 1999 , 2002b ) ) . clearly , the limit @xmath110 is independent of @xmath12 , therefore we may set @xmath270 .
the characteristic exponent of a process without the gaussian component is of the form @xmath271 where @xmath272 grows as @xmath273 in the complex plane with the cuts @xmath66 $ ] and @xmath274 , for processes of order @xmath137 . in the complex plane outside two poles @xmath275 . for vgp
, @xmath272 grows as @xmath276 , and for exponential jump processes , @xmath272 decays as @xmath277 .
set @xmath278 , and write ( [ callprice ] ) as @xmath279 } { -i\xi(1-i\xi)}d\xi.\ ] ] let @xmath31 be a kobol or nts lvy process of order @xmath280
. then we can find @xmath281 s.t .
since @xmath283 is bounded away from zero on the line @xmath284 , the integrand in ( [ callprice2 ] ) admits an upper bound via @xmath285 it follows that the integral over @xmath286 makes a contribution of order @xmath287 since @xmath288 , we have @xmath289 on the set @xmath290 the @xmath79 is fixed , therefore for sufficiently small @xmath102 , @xmath291 is negative and bounded away from zero , and therefore , @xmath292 is negative and bounded away from zero , uniformly in @xmath293 and @xmath294 in the half - plane @xmath295 .
this means that we can use the cauchy theorem and push the line of the integration in ( [ au1 ] ) down : @xmath296 , and obtain that the integral in ( [ au1 ] ) is zero .
hence , @xmath297 by using essentially the same argument as in boyarchenko and levendorski ( 2002b ) , where a similar integral ( 5.33 ) was transformed to the integral over the banks of the cut @xmath298 $ ] in the complex plane , and after the change of variables @xmath299 and @xmath300 on the left and right banks , respectively , to an integral over @xmath301 ( equation ( 5.37 ) in boyarchenko and levendorski ( 2002b ) ) , we obtain , for @xmath302 , @xmath303 as @xmath304 , @xmath305 , therefore by passing to the limit in ( [ genc2 ] ) , we obtain ( [ genc ] ) .
recall that the calculations above were made for a kobol or nts lvy process of order @xmath306 .
if @xmath31 is a vgp , then we make essentially the same argument by using any @xmath307 .
if @xmath31 is a hyperbolic process , nig , or kobol or nts lvy process of order @xmath59 , then we start with the integration by part in ( [ callprice2 ] ) , by using @xmath308 after that we obtain an absolutely converging integral , and can repeat the constructions above with straightforward changes . in the end , we obtain @xmath309dz+o(1).\ ] ] by integrating by part back , and passing to the limit @xmath293 , we obtain ( [ genc ] ) .
to finish the proof of ( [ genc ] ) , it remains to show that the result does not change when we add a gaussian component .
let @xmath31 be any of the processes above , without a gaussian component , and let @xmath310 be a lvy process with the same  drift " term and jump density as @xmath31 , and with the diffusion coefficient @xmath33 .
let @xmath311 be the probability density of @xmath310 , and @xmath312 be the probability density of the brownian motion with the volatility @xmath33 .
the @xmath312 is the convolution of @xmath177 and @xmath312 , therefore for @xmath81 , @xmath313 take @xmath314 . in the region @xmath315 , the gaussian density is @xmath316 , for any @xmath317 , therefore we may replace the outer integral with the integral over @xmath318 , and add @xmath319 .
but then on the support of the integrand , for a fixed @xmath79 , @xmath320 is bounded away from zero , uniformly in @xmath3 .
define @xmath321 as @xmath110 but assuming that the log - price follows the jump part of the process
. then we can represent the inner integral in the form @xmath322 and therefore , @xmath323 now , define @xmath324 to be @xmath325 for all @xmath326 . then for small @xmath3 , the integral above does not change but now @xmath327 is a continuous bounded function on @xmath328 . by using the super - exponential decay of the gaussian density , we can restore the @xmath329 , and after that pass to the limit and obtain @xmath330 .
this finishes the proof of ( [ genc ] ) .
[ [ proof - of - cckblcpnts - and - cpjumpccjump - for - jump - processes - with - non - zero - drift ] ] proof of ( [ cckbl])([cpnts ] ) , and ( [ cpjump])([ccjump ] ) for jump processes with non - zero drift ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
let @xmath31 be a vgp .
we have @xmath333.\ ] ] if @xmath118 , then the terms with @xmath247 cancel out , and @xmath334 hence , @xmath335 . if @xmath117 , then the terms with @xmath262 cancel out , and @xmath336 thus , @xmath337 .
\2 . let @xmath31 be a kobol process . in the case
@xmath118 , @xmath338\\   & = & ic\gamma(-\nu)[({{\lambda_+}}-z - i0)^\nu-({{\lambda_+}}-z+i0)^\nu]\\   & = & ic\gamma(-\nu)(z-{{\lambda_+}})^\nu i(e^{-i\pi\nu}-e^{i\pi\nu})\\ & = & 2c\gamma(-\nu)\sin(\pi\nu)(z-{{\lambda_+}})^\nu.\end{aligned}\ ] ] similarly , in the case @xmath117 , @xmath339    \3 .
let @xmath31 be an nts lvy process .
if @xmath340 , then @xmath341\\ & = & i{\delta}[({\alpha}^2-({\beta}-z+i0)^2)^{\nu/2}-({\alpha}^2-({\beta}-z - i0)^2)^{\nu/2}]\\ & = & i{\delta}({\alpha}-{\beta}+z)^{\nu/2}[({\alpha}+{\beta}-z+i0)^{\nu/2}-({\alpha}+{\beta}-z - i0)^{\nu/2}]\\ & = & i{\delta}({\alpha}-{\beta}+z)^{\nu/2}(z-{\alpha}-{\beta})^{\nu/2}(e^{i\pi\nu/2}-e^{-i\pi\nu/2})\\ & = & -2\sin\frac{\pi\nu}{2}((z-{\beta})^2-{\alpha}^2)^{\nu/2}.\end{aligned}\ ] ] if @xmath342 , then similarly , we obtain the same expression but without the minus sign .
finally , if @xmath31 is an exponential jump - diffusion process , then the same proof as in lemma [ gen ] shows that for the call , it suffices to consider the pure jump process , and that in this case , ( [ callprice3 ] ) holds . since @xmath343 the integrand in ( [ callprice3 ] ) has the only pole @xmath344 in the plane @xmath345 .
we push down the line of integration , and on crossing the pole , apply the residue theorem .
the result is @xmath346 and by passing to the limit , we obtain ( [ ccjump ] ) .
( [ cpjump ] ) is proved similarly .
cont , r. , potters , m. , and bouchaud , j .-
p . : scaling in stock market data : stable laws and beyond " , in dubrulle , b. , f. graner , and d. sornette ( eds . ) .
scale invariance and beyond ( proceedings of the cnrs workshop on scale invariance , les houches , march 1997 ) .
berlin : springer 1997                \28 .
lamberton , d. : critical price for an american option near maturity . in bolthausen ,
e.  , m.  dozzi , f.  russo ( eds . ) seminar on stochastic analysis , random fields and applications .
boston basel berlin : birkhaser 1995