• Refine Query
  • Source
  • Publication year
  • to
  • Language
  • 1
  • Tagged with
  • 2
  • 2
  • 1
  • 1
  • 1
  • 1
  • 1
  • 1
  • 1
  • 1
  • 1
  • 1
  • 1
  • 1
  • 1
  • About
  • The Global ETD Search service is a free service for researchers to find electronic theses and dissertations. This service is provided by the Networked Digital Library of Theses and Dissertations.
    Our metadata is collected from universities around the world. If you manage a university/consortium/country archive and want to be added, details can be found on the NDLTD website.
1

Path Integral Approach to Levy Flights and Hindered Rotations

Janakiraman, Deepika January 2013 (has links) (PDF)
Path integral approaches have been widely used for long in both quantum mechanics as well as statistical mechanics. In addition to being a tool for obtaining the probability distributions of interest(wave functions in the case of quantum mechanics),these methods are very instructive and offer great insights into the problem. In this thesis, path integrals are extensively employed to study some very interesting problems in both equilibrium and non-equilibrium statistical mechanics. In the non-equilibrium regime, we have studied, using a path integral approach, a very interesting class of anomalous diffusion, viz. the L´evy flights. In equilibrium statistical mechanics, we have evaluated the partition function for a class of molecules referred to as the hindered rotors which have a barrier for internal rotation. Also, we have evaluated the exact quantum statistical mechanical propagator for a harmonic potential with a time-dependent force constant, valid under certain conditions. Diffusion processes have attracted a great amount of scientific attention because of their presence in a wide range of phenomena. Brownian motion is the most widely known class of diffusion which is usually driven by thermal noise. However ,there are other classes of diffusion which cannot be classified as Brownian motion and therefore, fall under the category of Anomalous diffusion. As the name suggests, the properties of this class of diffusion are very different from those for usual Brownian motion. We are interested in a particular class of anomalous diffusion referred to as L´evy flights in which the step sizes taken by the particle during the random walk are obtained from what is known as a L´evy distribution. The diverging mean square displacement is a very typical feature for L´evy flights as opposed to a finite mean square displacement with a linear dependence on time in the case of Brownian motion. L´evy distributions are characterized by an index α where 0 <α ≤ 2. When α =2, the distribution becomes a Gaussian and when α=1, it reduces to a Cauchy/Lorentzian distribution. In the overdamped limit of friction, the probability density or the propagator associated with L´evy flights can be described by a position space fractional Fokker-Planck equation(FFPE)[1–3]. Jespersen et al. [4]have solved the FFPE in the Fourier domain to obtain the propagator for free L´evy flight(absence of an external potential) and L´evy flights in linear and harmonic potentials. We use a path integral technique to study L´evy flights. L´evy distributions rarely have a compact analytical expression in the position space. However, their Fourier transformations are rather simple and are given by e−D │p│α where D determines the width of the distribution. Due to the absence of a simple analytical expression, attempts in the past to study L´evy flights using path integrals in the position space [5, 6] have not been very successful. In our approach, we have tried to make use of the elegant representation of the L´evy distribution in the Fourier space and therefore, we write the propagator in terms of a two-dimensional path integral –one over paths in the position space(x)and the other over paths in the Fourier space(p). We shall refer to this space as the ‘phase space’. Such a representation is similar to the Hamiltonian path integral of quantum mechanics which was introduced by Garrod[7]. If we try to perform the path integral over Fourier variables first, then what remains is the usual position space path integral for L´evy flights which is rather difficult to solve. Instead, we perform the position space path integral first which results in expressions which are rather simple to handle. Using this approach, we have obtained the propagators for free L´evy flight and L´evy flights in linear and harmonic potentials in the over damped limit [8]. The results obtained by this method are in complete agreement with those obtained by Jesepersen et al. [4]. In addition to these results, we were also able to obtain the exact propagator for L´evy flights in a harmonic potential with a time-dependent force constant which has not been reported in the literature. Another interesting problem that we have considered in the over damped limit is to obtain the probability distribution for the area under the trajectory of a L´evy particle. The distributions, again, were obtained for free L´evy flight and for L´evy flights subjected to linear and harmonic potentials. In the harmonic potential, we have considered situations where the force constant is time-dependent as well as time-independent. Like in the case of the over damped limit, the probability distribution for L´evy flights in the under damped limit of friction can also be described using a fractional Fokker-Planck equation, although in the full phase space. However, this has not yet been solved for any general value of α to obtain the complete propagator in terms of both position and velocity. Using our path integral approach, the exact full phase space propagators have been obtained for all values of α for free L´evy flights as well as in the presence of linear and harmonic potentials[8]. The results that we obtain are all exact when the potential is at the most harmonic. If the potential is higher than harmonic, like the cubic potential, we have used a semi classical evaluation where, we extremize the action using an optimal path and further, account for fluctuations around this optimal path. Such potentials are very useful in describing the problem of escape of a particle over a barrier. The barrier crossing problem is very extensively studied for Brownian motion (Kramers problem) and the associated rate constant has been calculated in a variety of methods, including the path integral approach. We are interested in its L´evy analogue where we consider the escape of a particle driven by a L´evy noise over a barrier. On extremizing the action which depends both on phase space variables, we arrived at optimal paths in both the position space as well as the space of the conjugate variable, p. The paths form an infinite hierarchy of instant on paths, all of which have to be accounted for in order to obtain the correct rate constant. Care has to be taken while accounting for fluctuations around the optimal path since these fluctuations should be independent of the time-translational mode of the instant on paths. We arrived at an ‘orthogonalization’ scheme to perform the same. Our procedure is valid in the limit when the barrier height is large(or when the diffusion constant is very small), which would ensure that there is small but a steady flux of particles over the barrier even at very large times. Unlike the traditional Kramers rate expression, the rate constant for barrier crossing assisted by L´evy noise does not have an exponential dependence on the barrier height. The rate constant for wide range of α, other than for those very close to α = 2, are proportional to Dμ where, µ ≈ 1 and D is the diffusion constant. These observations are consistent with the simulation results obtained by Chechkin et al. [9]. In addition, our approach when applied to Brownian motion, gives the correct dependence on D. In equilibrium statistical mechanics we have considered two problems. In the first one, we have evaluated the imaginary time propagator for a harmonic oscillator with a time-dependent force constant(ω2(t))exactly, when ω2(t) is of the form λ2(t) - λ˙(t)where λ(t) is any arbitrary function of t. We have made use of Hamiltonian path integrals for this. The second problem that we considered was the evaluation of the partition function for hindered rotors. Hindered rotors are molecules which have a barrier for internal rotation. The molecule behaves like free rotor when the barrier is very small in comparison with the thermal energy, and when the barrier is very high compared to thermal energy, it behaves like a harmonic oscillator. Many methods have been developed in order to obtain the partition function for a hindered rotor. However, most of them are some what ad-hoc since they interpolate between free-rotor and the harmonic oscillator limits. We have obtained the approximate partition function by writing it as the trace of the density matrix and performing a harmonic approximation around each point of the potential[10]. The density matrix for a harmonic potential is in turn obtained from a path integral approach[11]. The results that we obtain using this method are very close to the exact results for the problem obtained numerically. Also, we have devised a proper method to take the indistinguishability of particles into account in internal rotation which becomes very crucial while calculating the partition function at low temperatures.
2

Stochastic modelling of financial time series with memory and multifractal scaling

Snguanyat, Ongorn January 2009 (has links)
Financial processes may possess long memory and their probability densities may display heavy tails. Many models have been developed to deal with this tail behaviour, which reflects the jumps in the sample paths. On the other hand, the presence of long memory, which contradicts the efficient market hypothesis, is still an issue for further debates. These difficulties present challenges with the problems of memory detection and modelling the co-presence of long memory and heavy tails. This PhD project aims to respond to these challenges. The first part aims to detect memory in a large number of financial time series on stock prices and exchange rates using their scaling properties. Since financial time series often exhibit stochastic trends, a common form of nonstationarity, strong trends in the data can lead to false detection of memory. We will take advantage of a technique known as multifractal detrended fluctuation analysis (MF-DFA) that can systematically eliminate trends of different orders. This method is based on the identification of scaling of the q-th-order moments and is a generalisation of the standard detrended fluctuation analysis (DFA) which uses only the second moment; that is, q = 2. We also consider the rescaled range R/S analysis and the periodogram method to detect memory in financial time series and compare their results with the MF-DFA. An interesting finding is that short memory is detected for stock prices of the American Stock Exchange (AMEX) and long memory is found present in the time series of two exchange rates, namely the French franc and the Deutsche mark. Electricity price series of the five states of Australia are also found to possess long memory. For these electricity price series, heavy tails are also pronounced in their probability densities. The second part of the thesis develops models to represent short-memory and longmemory financial processes as detected in Part I. These models take the form of continuous-time AR(∞) -type equations whose kernel is the Laplace transform of a finite Borel measure. By imposing appropriate conditions on this measure, short memory or long memory in the dynamics of the solution will result. A specific form of the models, which has a good MA(∞) -type representation, is presented for the short memory case. Parameter estimation of this type of models is performed via least squares, and the models are applied to the stock prices in the AMEX, which have been established in Part I to possess short memory. By selecting the kernel in the continuous-time AR(∞) -type equations to have the form of Riemann-Liouville fractional derivative, we obtain a fractional stochastic differential equation driven by Brownian motion. This type of equations is used to represent financial processes with long memory, whose dynamics is described by the fractional derivative in the equation. These models are estimated via quasi-likelihood, namely via a continuoustime version of the Gauss-Whittle method. The models are applied to the exchange rates and the electricity prices of Part I with the aim of confirming their possible long-range dependence established by MF-DFA. The third part of the thesis provides an application of the results established in Parts I and II to characterise and classify financial markets. We will pay attention to the New York Stock Exchange (NYSE), the American Stock Exchange (AMEX), the NASDAQ Stock Exchange (NASDAQ) and the Toronto Stock Exchange (TSX). The parameters from MF-DFA and those of the short-memory AR(∞) -type models will be employed in this classification. We propose the Fisher discriminant algorithm to find a classifier in the two and three-dimensional spaces of data sets and then provide cross-validation to verify discriminant accuracies. This classification is useful for understanding and predicting the behaviour of different processes within the same market. The fourth part of the thesis investigates the heavy-tailed behaviour of financial processes which may also possess long memory. We consider fractional stochastic differential equations driven by stable noise to model financial processes such as electricity prices. The long memory of electricity prices is represented by a fractional derivative, while the stable noise input models their non-Gaussianity via the tails of their probability density. A method using the empirical densities and MF-DFA will be provided to estimate all the parameters of the model and simulate sample paths of the equation. The method is then applied to analyse daily spot prices for five states of Australia. Comparison with the results obtained from the R/S analysis, periodogram method and MF-DFA are provided. The results from fractional SDEs agree with those from MF-DFA, which are based on multifractal scaling, while those from the periodograms, which are based on the second order, seem to underestimate the long memory dynamics of the process. This highlights the need and usefulness of fractal methods in modelling non-Gaussian financial processes with long memory.

Page generated in 0.0478 seconds