Friday, August 5, 2016

Space Physics: Gradient-B Drift

Now that you understand why a proton or electron from the Sun that bombards the Earth's magnetic field lines undergoes cyclotron motion, let's make things slightly more complicated.  Now, the magnetic fields that we'll deal with have a gradient.  Simply put, this means that the magnetic field line is stronger on one end than on the other sections of the magnetic field line.  In mathematical terms, the magnetic field vector \( \textbf{B} \) is no longer constant with respect to position; thanks to the gradient, changes in \( \textbf{B} \) only depend on location.  (For now, this vector is only a function of position--it remains constant over time.)   Just to clarify, I'll add that the Lorentz force--discussed in the previous post--applied by the \( \textbf{B} \) field onto the particle is dependent on the direction, or angle with respect to the \( \textbf{B} \) field, of the particle's velocity, but this force is not dependent on the location of the particle while it is inside of that magnetic field as long as that magnetic field has no gradient.
Pheww!! ... tough explanation ... it can be hard sometimes to explain the concept of vectors in contrast to vectors that are functions of position to people that might not utilize vector mathematics in their everyday lives.

The simulation I am about to show has a proton traveling in the presence of a magnetic field with a gradient.  For the sake of laziness, I'll just call this magnetic field a Gradient-B field ( \( \nabla B \) ) from hereon.  It is important to note that because we are now dealing with \( \nabla B \) field--a magnetic field where its strength changes depending on your location within the field--the Lorentz force applied to the proton will change as the proton travels from one location to another within this field.  When the proton enters a \( \nabla B \) field, it drifts in a direction other than the direction of the magnetic field vector \( \textbf{B} \).  In my previous post, the particle only drifted parallel to \( \textbf{B} \) but as it drifted in that direction it also traveled in a circle--hence cyclotron motion.

The average velocity due to the \( \nabla B \) drift is described as
\[ v_a  =  \frac{\mu}{q} \frac{ \textbf{B}\times \nabla_{\perp}B }{B^2}  \]
where mu \(\mu\) is
\[ \mu = \frac{m v_{\perp}^2}{2B} \]
and \( \nabla_{\perp}B \) is the strength of the \( \nabla B \)  field in the direction perpendicular to the \( \textbf{B} \) field.

Top View
Star is the beginning and circle is end of simulation.
 For this simulation, I used the following equation to describe the motion of a proton:
\begin{equation} \textbf{r}(t)-\textbf{r}(0) = \frac{v_{\perp}}{\Omega_c}[1 - cos(\Omega_c t)]\hat{\textbf{e}}_x + \frac{v_{\perp}}{\Omega_c}sin(\Omega_c t)\hat{\textbf{e}}_y + v_{\|}t\hat{\textbf{e}}_z + \frac{\mu}{q} \frac{ \textbf{B}\times \nabla_{\perp}B }{B^2}t
however!! it is important to know that the initial velocity I set for my proton was \( v_{\perp} = \textbf{v}(0) = 4000\) \( \hat{\textbf{e}}_y \) m/s which means that the velocity of the proton when \(0\) seconds has gone by is 4,000 meters per second only in the +y-direction.  I will also note that I set the \( \textbf{B} \) field to point in the +z-direction but the \( \nabla B \) vector is pointed in the +x-direction.  These details are important for understanding why the particle traveled the path that is shown in the video.

Because of the mathematics, I knew how the proton should have behaved but I needed to confirm my hypothesis with the video.
Here are some questions I asked myself when I first saw the video my code created:

Did the proton really undergo cyclotron motion?
Yes.  That's why it travels in that loopty-loop fashion, as seen in the "Top View" image.

Why did the proton drift in the +y-direction?
This is because the velocity drift due to \( \nabla B \) is in that direction.  This is shown in the \( v_a \) equation because of the cross product \( \textbf{B}\times \nabla_{\perp}B \).  In terms of the directions, this cross product adheres to the right hand rule: \( \hat{\textbf{e}}_z \times \hat{\textbf{e}}_x = \hat{\textbf{e}}_y \) , which is the +y-direction.

Why didn't the proton ever move in the \( \pm \)z-direction?
This is because my initial condition for velocity \( \textbf{v}(0) = 4000\) \(\hat{ \textbf{e}}_y \) m/s is not in the \(\pm \)z-direction nor does the Lorentz force \( \textbf{F} = q(\textbf{v} \times \textbf{B}) \) ever at any point in time push the proton in the \(\pm\)z-direction.

Why does the proton maneuver in both x and y directions if the initial condition of the proton started it off going in only the +y-direction?
This is because the Lorentz force accelerates the proton into +x-direction at the beginning ( \(\textbf{v} \times \textbf{B} == \hat{y} \times \hat{z} = \hat{x} \) ) and as time goes on the Lorentz force continues to accelerate the particle in different directions until eventually the particle has gone a circular path.

Hopefully, I answered all of the questions that may or may not have already popped into your head when watching my code run.  I think that this is enough fun for today.  Next time, we'll see what flabbergasting events transpire in the simulation when the initial velocity changes.

Friday, July 8, 2016

Space Physics: Cyclotron Motion

I'll discuss a basic topic in Space Physics today that is certainly not a basic idea for the average person outside of science to understand.  However! that should not matter nor should it intimidate you!  Any person (such as myself) that attempts to teach ANY scientific concept should be capable of explaining their science in such a way that even children with little scientific background can understand the concept.

The basic topic is the behavior of a proton as it approaches the Earth's magnetic field.  Remember, an atom--those boxes you see on a periodic table--can consist of a positively charged proton and a negatively charged electron.  Because I am discussing a topic in the field of space physics, I will also remind you that the Earth's magnetic field is always being hit by protons and electrons from the Sun.

Today, I will describe one of the ways a proton behaves when encountering an unchanging magnetic field line. 

When a proton approaches a magnetic field line, the proton spirals around a magnetic field line and continues in the direction of that line.  This is known as cyclotron motion.

Top View
The video illustrates a simulation in which a magnetic field line is aligned with the z-axis and is pointing upward in the +z direction (also denoted as the \( +\hat{\textbf{e}}_z \) direction).  This is why the proton travels in the +z direction over time while moving in a circle.  Mathematically, the +z direction of the particle's motion is shown via the \( v_{\|}t\hat{\textbf{e}}_z \) term at the end of my cyclotron motion derivation.

The image to the right provides a top view of the simulation.  This perspective allowed me to verify if the proton truly has a circular trajectory in 2D space.  It is easy to see (thanks to the numbers on the axes) that the proton traveled 0.23 meters in the +y direction as well as the -y direction.  The trajectory looks visually like an ellipse (ugghhhh, thanks MATLAB...) but after looking at the numbers, I found that the proton truly did travel 0.23 meters from the center of the circle in all directions.  The radius of this circular cyclotron motion is known as the Larmor orbital radius:
\[r_L = \Big| \frac{v_{\perp}}{\Omega} \Big| = \Big| \frac{mv_{\perp}}{qB} \Big| \]
where \(\Omega\) is the frequency of the orbit, m is the mass of the particle (in this case a proton), \(v_{\perp}\) is the component of velocity of the particle that goes in a direction that is perpendicular to the magnetic field B, and q is the charge of that particle (in this case +1.602\(\times \)10\(^{-19}\) Coulombs).

Now, let's see how the electron behaves in the presence of a magnetic field line.  Notice anything different?
Top View
Hmmmm... Interesante.  At first glance, it seems to have done the same thing as the proton.  One difference seen is that the electron's cyclotron motion is counterclockwise while the proton's motion is clockwise.  Another distinction is noticed when looking at the simulation from the top view.  The Larmor radius is now 1.26\(\times \) 10\(^{-4}\) meters.

Mathematically speaking, the reason why the electron's Larmor radius is shorter than the proton's Larmor radius is because the mass of an electron is about 1800 times less than the mass of a proton, and because the \(r_L\) is directly proportional to the m, the super tiny electron mass makes its Larmor radius super tiny as well.

The difference in clockwise vs counterclockwise motion between proton and electron is due to the Lorentz force:
\[ \vec{F} = q(\vec{E} + \vec{v} \times \vec{B} ) \]
where \(\vec{E}\) is the electric field.  Basically, this equation means that whenever a charged particle is in the presence of an electric or magnetic field, a force will be applied to that particle which will cause it to either to accelerate or change directions.  The arrows on top of the variables indicate that they are vectors, which means that the electric field strength, the particle's velocity, and the magnetic field strength are all pointed in some direction.  In my simulation, I have not included an electric field.  Therefore, in this case, \(\vec{E}=0\) and \( \vec{F} = q ( \vec{v} \times \vec{B} ) \).  Between the proton and electron simulations I did not change the velocity nor the magnetic field.  The only thing that was different is the charge q.  Since protons and electrons have opposite charge, the magnetic field's force \(\vec{F}\) was pushing them in opposite directions!  Before clicking on my previous videos, you could have predicted if the motion was going to be clockwise or counterclockwise as long as you would have used the right hand rule for the cross product \( \vec{v} \times \vec{B} \).

Later, I will discuss other behaviors that particles from the Sun exhibit when they encounter the Earth's magnetic field.

Deriving Position from Velocity of Pure Cyclotron Motion
\begin{equation} \textbf{v}(t) = v_{\perp}sin(\Omega_c t)\hat{\textbf{e}}_x + v_{\perp}cos(\Omega_c t)\hat{\textbf{e}}_y + v_{\|}\hat{\textbf{e}}_z
\[ \int_{t=0}^{t} \textbf{v}(t) dt = \int_{t=0}^{t} \big( v_{\perp}sin(\Omega_c t)\hat{\textbf{e}}_x + v_{\perp}cos(\Omega_c t)\hat{\textbf{e}}_y + v_{\|}\hat{\textbf{e}}_z \big)dt  \]
\[ \textbf{r}(t)-\textbf{r}(0) = -\frac{v_{\perp}}{\Omega_c}[cos(\Omega_c t) - 1]\hat{\textbf{e}}_x + \frac{v_{\perp}}{\Omega_c}[sin(\Omega_c t)- 0]\hat{\textbf{e}}_y + v_{\|}[t - 0]\hat{\textbf{e}}_z \]
\begin{equation}\label{Eq2: Analy path} \textbf{r}(t)-\textbf{r}(0) = \frac{v_{\perp}}{\Omega_c}[1 - cos(\Omega_c t)]\hat{\textbf{e}}_x + \frac{v_{\perp}}{\Omega_c}sin(\Omega_c t)\hat{\textbf{e}}_y + v_{\|}t\hat{\textbf{e}}_z

Saturday, January 2, 2016

Long Time, No See

Wow! it feels good to be back.

Many of you know, have seen, and have heard about the grand schemes I have been up to in the past few months, which has kept me too busy to write more blog posts here.  I'll take a moment to get you all up to speed on the big moves I have been making.

Some of you have seen the inception of my Research Website and the drastic changes I have made to it since I started writing the HTML/CSS/JavaScript code for that in September.  I even have a cool background image for the site now!  Look familiar?
Seeing this website come to fruition is like a small dream come true for me.  This started as a simple idea and at the start I only knew a couple HTML commands (< br > for break and < a > for links).  The website is still a work in progress but I love how it's turning out. 

Another endeavor (that only my research advisors probably noticed) is the amount of money I have racked up this year for research--as seen on my CV.  Receiving that much $ $ $ $ $ for research funding as an undergraduate, within one year's time, required me to devote a LOT of time to writing about my research.  (And my CV doesn't even include the generous stipend I received from my summer internship! )

Needless to say, all of that writing about my research took a lot of time away from my academic studies during this past fall semester.  The writing combined with the actual obligations of conducting the research certainly made this semester much more difficult.   However, I am happy to inform you that I persevered.  I didn't get a 4.0 GPA this semester, but I am still happy with the fruits of my labor.  Excellent grades, great progress in research, and an abundance of research funding came out of my hard work this summer and semester.

The McNair Scholars Program, formally known as the Post-Baccalaureate Achievement Program, also kept me busy.  Much of my extra-curricular writing this year was for the McNair Program.  Additionally, because of my completed McNair research project I had to present my research at the FIU McNair conference.
Me delivering my talk.  For those of you who've seen my website, the photo should look familiar.
 Hanging out with other McNair students that completed their research projects was a fun experience too.

We had a great time together.  The conference itself was mediocre though.  (With the AAS Meeting being my first science conference, I admit that I'm probably spoiled when it comes to this stuff.)

In this post, I have given you a brief synopsis for a few of my most time consuming (although fun) extra-curricular activities for this past semester.  However, I didn't even tell you about the hard work I put into my NSF Graduate Research Fellowship application.  Nor did I discuss the essays I have written this semester for graduate school applications.  But most importantly, YOU JUST DON'T KNOW how tortuous the Quantum Mechanics class I took this semester was.  To give you a glimpse of the pain, here is a quote from my professor: "I expect the homework to take 15 hours for all of you."
Unfortunately, this was not a joke.  Each homework assignment given once every week to us--undergraduates who are taking 4 or more classes this semester--truly did take a minimum of 15 hours to complete.  The teacher believed that this absurd amount of work would prepare us well for graduate school.  Literally, every student in his class disagreed with this notion at the beginning, middle, and very end of his class this semester because of the 4+ classes we all were taking.

Without a doubt, this was the hardest semester of my college and high school career.  To be frank, I am willing to say that I had a great/fun semester ONLY because of how good things turned out for me in the end.  Ultimately, the, sometimes overwhelming, amount of work and stress I experienced this semester was mostly my fault.  I decided to work on 3 research projects simultaneously this semester (including 1 that is not mentioned on my CV).  I made the choice to apply for research grants and write essays while taking this tortuous Quantum Mechanics class.  I made the choice to do all of these extra-curricular obligations even though I knew that I had graduate school applications coming up.

Perhaps this semester, I bit a little more than I could chew... Regardless of whether or not that is true, great things were--and still are--happening for me on a weekly basis thanks to the immense amount of time and effort I put into all that I do.  I'm looking forward to seeing more fruits of this semester's labor manifest once the results of my NSF and graduate school applications come out.