Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Simulation of the fission as a stochastic process

  1. Jun 10, 2005 #1

    I'm a french student, I'm actually not sure this is the good place to ask my question but as it deals with the nuclear fission I try here... don't hesitate to tell me if there is a better forum... thx..

    well, I'm trying to solve numerically the Langevin equation, initially for the brownian motion, it is also used in nuclear physics as you probably know, to describe the fission of nuclei...

    the Langevin equation is :

    [tex]\ddot{x}+\beta \dot{x} + \frac{1}{m}\frac{\partial U}{\partial x} = \Gamma(t)[/tex]

    where [tex]\beta=\frac{\gamma}{m}[/tex] is a friction term (from a stocke's law) divided by the mass of the brownian particle. (so [tex]\beta[/tex] is [tex]s^{-1}[/tex])

    [tex]U(x)[/tex] is an extern potential

    and [tex]\Gamma(t)[/tex] is a stochastic force divided my the mass.

    well, now I decide, as Kramers did, to tell that my potential U is represented by two parabolas as you can see on http://nicolas.aunai.free.fr/courb.htm [Broken]

    I call B the height of the maximum, and [tex]x_b[/tex] the correspondant X-coordonnate (assuming that the minimum will be at x=0)

    I decide that B will be my caracteristic energy for the problem, and [tex]x_b[/tex] my caracteristc length.

    I want, in order to solve my equation numerically, to write it without dimension, so I introduce the following new variables :

    [tex]Q = \frac{x}{x_b}[/tex]

    [tex]\Pi = \frac{U}{B}[/tex]

    [tex]t'=\omega t[/tex] where [tex]\omega[/tex] is the caracteristic pulsation of my parabolas

    i.e. the equations of the parabolas are :

    [tex]U_1(x)=\frac{1}{2}m\omega^2 x^2[/tex] for [tex]x<x_0[/tex]

    [tex]U_2(x)=-\frac{1}{2}m\omega^2(x-x_b)^2 + B[/tex] for [tex]x>x_0[/tex]

    assuming the continuity of the potential, the functions and the derivatives shoud be equal between them at [tex]x_0[/tex] the junction point. these conditions give us the height B which is equals to :

    [tex]\frac{1}{4}m\omega^2 x_b^2[/tex]

    well now we have to re-write the langevin equation, replacing by the new variables... and here began my problems....

    [tex]\ddot{x}+\beta \dot{x} + \frac{1}{m}\frac{\partial U}{\partial x} = \Gamma(t)[/tex]

    [tex]\dot{x}[/tex] means [tex]\frac{dx}{dt}[/tex], so since we have [tex]t'=\omega t[/tex], the derivation by t' gives us an [tex]\omega[/tex]

    [tex]\omega^2 x_b\ddot{Q}+\beta \omega x_b \dot{Q} + \frac{B}{m}\frac{\partial \Pi}{\partial Q} \frac{\partial Q}{\partial x} = \Gamma(t)[/tex]

    we know that [tex]\frac{\partial Q}{\partial x} = \frac{1}{x_b}[/tex]

    so we have :

    [tex]\omega^2 x_b\ddot{Q}+\beta \omega x_b \dot{Q} + \frac{B}{mx_b}\frac{\partial \Pi}{\partial Q}= \Gamma(t)[/tex]

    and with the B value we have :

    [tex]\omega^2 x_b\ddot{Q}+\beta \omega x_b \dot{Q} + \frac{\omega^2 x_b}{4}\frac{\partial \Pi}{\partial Q}= \Gamma(t)[/tex]

    finally, dividing by [tex]\omega^2 x_b[/tex] we obtain :

    [tex]\ddot{Q}+\frac{\beta}{\omega}\dot{Q} + \frac{1}{4}\frac{\partial \Pi}{\partial Q}= \frac{1}{\omega^2 x_b}\Gamma(t)[/tex]

    which I think is almost the result that I'm looking for.. since the first term is without dimention, [tex]\frac{\beta}{\omega}[/tex] is without dimension too, and the third term too..

    my problem is the stochastic term...

    now I don't know how to do with it.... if you have an idea...

    I now, that the autocorrelation function of the Langevin Force is :

    [tex]<\Gamma(t)\Gamma(t')> = 2\beta T \delta(t_1-t_2)[/tex]

    since [tex]t_i = \frac{t_i'}{\omega}[/tex], can I write the following equation :

    [tex]<\Gamma(t)\Gamma(t')> = 2\beta T \delta(\frac{t_1'}{\omega}-\frac{t_2'}{\omega})[/tex]

    and with [tex]\delta(a x)=\frac{1}{a}\delta(x)[/tex]

    I find :

    [tex]<\Gamma(t)\Gamma(t')> = 2\omega \beta T \delta(t_1'-t_2')[/tex]

    so [tex]\Gamma(t)[/tex] should be proportionnal to [tex]\sqrt{(2\omega\beta T \delta(t_1'-t_2'))}[/tex]

    in the nuclear system, k (the Boltzmann constant) is equal to 1, so we can also introduce the new variable :


    so we can put it in our relation :

    [tex]\Gamma(t)\ \ \ \alpha\ \ \ \sqrt{(2\omega\beta B\Theta\delta(t_1'-t_2'))}[/tex]

    well... I'm not sure where I'm going with this... can someone help me ?

    thanks !
    Last edited by a moderator: May 2, 2017
  2. jcsd
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Can you offer guidance or do you also need help?