Analysis of Fluctuations in the Early Universe Using Fokker-Planck Formalism

The thesis begins with a review of basic elements of general theory of relativity (GTR) which forms the basis for the theoretical interpretation of the observations in cosmology. The first chapter also discusses the standard model in cosmology, namely the Friedmann model, its predictions and problems. We have also made a brief discussion on fractals and inflation of the early universe in the first chapter. In the second chapter we discuss the formulation of a new approach to cosmology namely a stochastic approach. In this model, the dynam ics of the early universe is described by a set of non-deterministic, Langevin type equations and we derive the solutions using the Fokker—Planck formalism. Here we demonstrate how the problems with the standard model, can be eliminated by introducing the idea of stochastic fluctuations in the early universe. Many recent observations indicate that the present universe may be approximated by a many component fluid and we assume that only the total energy density is conserved. This, in turn, leads to energy transfer between different components of the cosmic fluid and fluctuations in such energy transfer can certainly induce fluctuations in the mean to factor in the equation of state p = wp, resulting in a fluctuating expansion rate for the universe. The third chapter discusses the stochastic evolution of the cosmological parameters in the early universe, using the new approach. The penultimate chapter is about the refinements to be made in the present model, by means of a new deterministic model The concluding chapter presents a discussion on other problems with the conventional cosmology, like fractal correlation of galactic distribution. The author attempts an explanation for this problem using the stochastic approach.

A Novel High Order Space-Time Spectral Method for the Time Fractional Fokker--Planck Equation

The fractional Fokker-Planck equation is an important physical model for simulating anomalous diffusions with external forces. Because of the non-local property of the fractional derivative an interesting problem is to explore high accuracy numerical methods for fractional differential equations. In this paper, a space-time spectral method is presented for the numerical solution of the time fractional Fokker-Planck initial-boundary value problem. The proposed method employs the Jacobi polynomials for the temporal discretization and Fourier-like basis functions for the spatial discretization. Due to the diagonalizable trait of the Fourier-like basis functions, this leads to a reduced representation of the inner product in the Galerkin analysis. We prove that the time fractional Fokker-Planck equation attains the same approximation order as the time fractional diffusion equation developed in [23] by using the present method. That indicates an exponential decay may be achieved if the exact solution is sufficiently smooth. Finally, some numerical results are given to demonstrate the high order accuracy and efficiency of the new numerical scheme. The results show that the errors of the numerical solutions obtained by the space-time spectral method decay exponentially.

Generalizations and extensions of the Fokker-Planck-Kolmogorov equations

The problem of determining probability density functions of general transformations of random processes is considered in this thesis. A method of solution is developed in which partial differential equations satisfied by the unknown density function are derived. These partial differential equations are interpreted as generalized forms of the classical Fokker-Planck-Kolmogorov equations and are shown to imply the classical equations for certain classes of Markov processes. Extensions of the generalized equations which overcome degeneracy occurring in the steady-state case are also obtained.

The equations of Darling and Siegert are derived as special cases of the generalized equations thereby providing unity to two previously existing theories. A technique for treating non-Markov processes by studying closely related Markov processes is proposed and is seen to yield the Darling and Siegert equations directly from the classical Fokker-Planck-Kolmogorov equations.

As illustrations of their applicability, the generalized Fokker-Planck-Kolmogorov equations are presented for certain joint probability density functions associated with the linear filter. These equations are solved for the density of the output of an arbitrary linear filter excited by Markov Gaussian noise and for the density of the output of an RC filter excited by the Poisson square wave. This latter density is also found by using the extensions of the generalized equations mentioned above. Finally, some new approaches for finding the output probability density function of an RC filter-limiter-RC filter system driven by white Gaussian noise are included. The results in this case exhibit the data required for complete solution and clearly illustrate some of the mathematical difficulties inherent to the use of the generalized equations.

Stability and convergence of an effective numerical method for the time-space fractional Fokker-Planck equation with a nonlinear source term

Fractional Fokker-Planck equations (FFPEs) have gained much interest recently for describing transport dynamics in complex systems that are governed by anomalous diffusion and nonexponential relaxation patterns. However, effective numerical methods and analytic techniques for the FFPE are still in their embryonic state. In this paper, we consider a class of time-space fractional Fokker-Planck equations with a nonlinear source term (TSFFPE-NST), which involve the Caputo time fractional derivative (CTFD) of order α ∈ (0, 1) and the symmetric Riesz space fractional derivative (RSFD) of order μ ∈ (1, 2). Approximating the CTFD and RSFD using the L1-algorithm and shifted Grunwald method, respectively, a computationally effective numerical method is presented to solve the TSFFPE-NST. The stability and convergence of the proposed numerical method are investigated. Finally, numerical experiments are carried out to support the theoretical claims.

The subordinated processes controlled by a family of subordinators and corresponding Fokker-Planck type equations

In this work, we consider subordinated processes controlled by a family of subordinators which consist of a power function of a time variable and a negative power function of an α-stable random variable. The effect of parameters in the subordinators on the subordinated process is discussed. By suitable variable substitutions and the Laplace transform technique, the corresponding fractional Fokker–Planck-type equations are derived. We also compute their mean square displacements in a free force field. By choosing suitable ranges of parameters, the resulting subordinated processes may be subdiffusive, normal diffusive or superdiffusive

Fokker-Planck equation for a test-particle weakly coupled to a magnetized plasma

We consider the derivation of a kinetic equation for a charged test particle weakly interacting with an electrostatic plasma in thermal equilibrium, subject to a uniform external magnetic field. The Liouville equation leads to a generalized master equation to second order in the `weak' interaction; a Fokker-Planck-type equation then follows as a `Markovian' approximation. It is shown that such an equation does not preserve the positivity of the distribution function f(x,v;t). By applying techniques developed in the theory of open systems, a correct Fokker-Planck equation is derived. Explicit expressions for the diffusion and drift coefficients, depending on the magnetic field, are obtained.

Uma solução híbrida para a equação de Fokker-Planck dependente da energia usando a técnica da transformada de Laplace e diferenças finitas

Neste trabalho é obtida uma solução híbrida para a equação de Fokker-Planck dependente da energia, muito utilizada em problemas de implantação iônica. A idéia consiste na aplicação da transformada de Laplace na variável de energia e aplicação de um esquema de diferenças finitas nas variáveis espacial e angular desta equação. Tal procedimento gera um problema matricial simbólico para a energia transformada. Para resolver este sistema, procede-se a inversão de Laplace da matriz (sI+A), onde s é um parâmetro complexo, I a matriz identidade e A uma matriz quadrada gerada pela discretização das variáveis espacial e angular. A matriz A não é diagonalizável, desta forma, contorna-se este problema decompondo esta matriz na soma de outras duas, onde uma delas é diagonalizável. É gerado então um método iterativo de inversão, semelhante ao método da fonte fixa associado ao método de diagonalização, do qual o resultado fornecido são os valores para o fluxo de partículas do sistema. A partir disto pode-se determinar a energia depositada no sistema eletrônico e nuclear do alvo. Para validar os resultados obtidos faz-se a simulação de implantação de íons de B em Si numa faixa energética de 1keV a 50MeV, comparam-se os resultados com simulação gerada numericamente pelo software SRIM2003.

Solutions of the Fokker-Planck equation for a Morse isospectral potential

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)