1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Calculus of variations, Fermat's principle

  1. Apr 6, 2010 #1

    fluidistic

    User Avatar
    Gold Member

    1. The problem statement, all variables and given/known data
    Fermat's principle establishes that the path taken by a light ray between 2 given points is such that the time that the light takes is the smallest possible.
    1)Demonstrate that a light ray which propagates through a medium with a constant refractive index follows a straight line.
    2)Demonstrate that Snell's law holds for a ray of light passing by 2 points separated by a plane inter-phase which separates 2 different mediums (i.e. with different refractive indexes).
    3)Determine the equation of the trajectory of a ray of light in the y-z plane assuming that the refractive index changes with z, i.e. n=n(z).



    2. Relevant equations

    [tex]S=\int _{t_1} ^{t_2} L dt[/tex].
    Optical path: [tex]L=\int _a ^b n(\vec r) ds[/tex]. I'd bet it should be [tex]L=\int _a ^b n(\vec r) d\vec s[/tex], right?
    Euler-Lagrange equations: [tex]\frac{d}{dt} \left ( \frac{\partial L}{\partial \dot q} \right ) - \frac{\partial L}{\partial q}=0[/tex].

    3. The attempt at a solution
    1)[tex]S=n \int _a ^b dr = n (b-a) \square[/tex].
    2)[tex]S= \int _a ^b n(\vec r) d\vec r[/tex]. Here I'm not sure, but I think I should reach something like [tex]n_1 \int a ^c d\vec r + n_2 \int c ^b d\vec r[/tex] but I don't really trust it. I don't think I can even use the fact that [tex]n_1=\frac{c}{v _1}[/tex]. It's a problem arising in Classical Mechanics, not optics. And we're dealing with variational calculus and Euler-Lagrange's equations.
    Do you have any idea how to tackle part 2? I'll try to do part 3 right after I've done part 2.
    Thanks for any help.
     
  2. jcsd
  3. Apr 7, 2010 #2

    Matterwave

    User Avatar
    Science Advisor
    Gold Member

    How do you do Fermat's principle while incorporating the principle of least action?

    What does action have to do with these problems?

    I'm not 100% but my professor once showed 1) and 2) using only maximizing and minimizing tools...without using any concept of action.
     
  4. Apr 7, 2010 #3

    Cyosis

    User Avatar
    Homework Helper

    Not sure what you're doing in 1).

    To determine the shortest path light will travel between two points you will need to minimize the action, which yields the EL equations, which you then have to solve.

    It is given that the refractive index is constant hence v is constant (v=cn).

    The distance between two points on a curve is given by ds=vdt:

    Fermat's principle then says that the path the light actually travels minimizes the time integral:

    [tex]
    \int \frac{ds}{v}=\int \frac{\sqrt{dx^2+dy^2}}{v}=\int \frac{\sqrt{1+y'^2}}{v}dx
    [/tex]

    2) If you've done 1) correctly you should have found that the path is a straight line. Now draw the path plus refraction in an xy diagram. Also add the normal and label the refractive angle and incident angle. Now you can calculate the length of the hypotenuse of the two triangles. Find an expression for the time and minimize this expression then calculate the sin of both angles.

    You can of course also split up the integral in two two integrals going from x0 to x1 and x1 to x2 then minimizing again and solving the equation, however since you already know the path the light will travel from 1) you may as well just use Pythagoras.

    3)See 1) however v is no longer constant.
     
    Last edited: Apr 7, 2010
  5. Apr 8, 2010 #4

    fluidistic

    User Avatar
    Gold Member

    Ok thanks to both, this should get me started. I'm going to start this weekend.
     
  6. Apr 13, 2010 #5

    fluidistic

    User Avatar
    Gold Member

    Ok, I need some further help for 1).
    So the Lagrangian is [tex]\sqrt {1+y' ^2}[/tex] since v is a constant and minimizing the integral is the same as minimizing the integral without the 1/v factor.
    So now I have to put this Lagrangian in Euler-Lagrange's equations.
    When [tex]q_i =x[/tex], there's nothing special, 0=0.
    When [tex]q_i = y[/tex], I get that [tex]\frac{\partial \sqrt {1 + y'^2}}{\partial y}=\frac{y '}{\sqrt {1+y'^2}}[/tex].
    So E-L equations become [tex]\frac{y '}{\sqrt {1+y'^2}}-\frac{d}{dt} \left ( \frac{y'}{\sqrt {1+y'^2}} \right ) =0[/tex]. Is the second term just 0? I'm not 100% sure because y' might depend on ds which depends on t.
    If it is 0 then y'=0 and so y, which is false. So clearly the second term is not null. How should I calculate it?
     
  7. Apr 13, 2010 #6

    gabbagabbahey

    User Avatar
    Homework Helper
    Gold Member

    No,

    [tex]\frac{\partial}{\partial y'}\sqrt{1+(y')^2}=\frac{y'}{\sqrt{1+(y')^2}}[/itex]

    But,

    [tex]\frac{\partial}{\partial y}\sqrt{1+(y')^2}=0[/itex]

    Since there is no explicit y-dependence.
     
  8. Apr 13, 2010 #7

    fluidistic

    User Avatar
    Gold Member

    Oops, you're right. A little question: how do you know that y' does not depend on y?
     
  9. Apr 13, 2010 #8

    fluidistic

    User Avatar
    Gold Member

    Now I get [tex]\frac{d}{dt} \left ( \frac{y'}{\sqrt {1+(y')^2}} \right )=0 \Rightarrow y' = C \Rightarrow y(x)=ax+b[/tex]. Which is the equation of a straight line as required.
    Now, is it safe to assume that y depends on x? I don't understand why.
     
  10. Apr 13, 2010 #9

    gabbagabbahey

    User Avatar
    Homework Helper
    Gold Member

    Why not? Keep in mind that even a function that doesn't depend on [itex]x[/itex], say [itex]y=2[/itex] , Can be thought of as being constant w.r.t [itex]x[/itex]. It is more general to assume that [itex]y[/itex] depends on both [itex]x[/itex] (explicitly) and [itex]t[/itex] (implicitly)...and if you find that [itex]y'(x)=0[/itex]...then you know that it is constant w.r.t. [itex]x[/itex].
     
  11. Apr 13, 2010 #10

    gabbagabbahey

    User Avatar
    Homework Helper
    Gold Member

    Because it's a partial derivative, you are only concerned with explicit dependence and not implicit dependence.
     
  12. Apr 13, 2010 #11

    fluidistic

    User Avatar
    Gold Member

    Thank you, it makes sense.

    I tried 2), but I get into a horrible Lagrangian.
    I took the point (x_1,y_1) as in the medium with refractive index n_1 and (x_2, y_2) the point in the medium with refractive index of n_2.
    [tex]t=\frac{s_1}{v_1}+\frac{s_2}{v_2}=\frac{\sqrt {x_1 ^2 + y_1 ^2}}{v_1}+\frac{\sqrt {x_2^2+ y_2^2}}{v_2}[/tex] where [tex]v_1=cn_1[/tex] and [tex]v_2=cn_2[/tex].
    But [tex]x_i[/tex] are given (they are not variables), so I can rewrite [tex]s_1[/tex] as [tex]\frac{x_1}{\cos \theta _i}[/tex] and [tex]s_2[/tex] as [tex]\frac{x_2}{\cos \theta _r}[/tex].
    Thus [tex]s=\frac{x_1}{\cos \theta _i}+\frac{x_2}{\cos \theta _r}[/tex].
    Since [tex]t=\int \frac{ds}{v}[/tex], I must find the total derivative of s.
    I found it to be ... [tex]d\theta _i \left ( \frac{x_1 \sin \theta _i}{(\cos \theta _i)^2} \right ) + d\theta _r \left ( \frac{x_2 \sin \theta _r}{(\cos \theta _r)^2} \right )[/tex]. I don't really trust this result. How could I find the Lagrangian from here? (factorizing by [tex]d\theta _i[/tex] for instance?)
     
  13. Apr 13, 2010 #12

    Cyosis

    User Avatar
    Homework Helper

    The two light rays are connected with each other. With the way you have set up your problem there is some gap between the two rays (which is really just one ray).

    Call the x coordinate where the light ray gets refracted x. The start of the beam before refraction x1 and the 'end' of the second beam x2. Then we can calculate the length of the path as [itex]s_1^2=(x-x_1)^2+y_1^2[/itex] and [itex] s_2^2=(x_1-x)^2+y_2^2[/itex]. Now you can minimize.
     
  14. Apr 13, 2010 #13

    fluidistic

    User Avatar
    Gold Member

    Thanks for your help.
    Did you mean [tex]s_2^2=(x_2-x)^2+y_2^2[/tex]?
    When you say "the start of the beam before refraction x1", do you mean the orthogonal projection of the point (x_1,y_1) into the x axis, axis of the refraction?
     
  15. Apr 14, 2010 #14

    Cyosis

    User Avatar
    Homework Helper

    Yes it should be x_2 and yes I mean the orthogonal projection. Just make a quick drawing if you haven't already.
     
  16. Apr 14, 2010 #15

    fluidistic

    User Avatar
    Gold Member

    Thank you very much. I had made a drawing and wanted to be sure I wasn't missing something.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Calculus of variations, Fermat's principle
  1. Variational Principle (Replies: 3)

  2. Variational Principle (Replies: 4)

  3. Fermat's Principle (Replies: 15)

  4. Variational Principle (Replies: 7)

Loading...