It looks like you're new here. If you want to get involved, click one of these buttons!

- 140.5K All Categories
- 103.4K Programming Languages
- 6.4K Assembler Developer
- 1.9K Basic
- 39.9K C and C++
- 2.9K C#
- 7.9K Delphi and Kylix
- 4 Haskell
- 9.6K Java
- 4.1K Pascal
- 1.3K Perl
- 2K PHP
- 527 Python
- 37 Ruby
- 4.4K VB.NET
- 1.6K VBA
- 20.8K Visual Basic
- 2.6K Game programming
- 313 Console programming
- 90 DirectX Game dev
- 1 Minecraft
- 110 Newbie Game Programmers
- 2 Oculus Rift
- 9K Applications
- 1.8K Computer Graphics
- 733 Computer Hardware
- 3.5K Database & SQL
- 535 Electronics development
- 1.6K Matlab
- 628 Sound & Music
- 257 XML Development
- 3.3K Classifieds
- 198 Co-operative Projects
- 192 For sale
- 190 FreeLance Software City
- 1.9K Jobs Available
- 602 Jobs Wanted
- 203 Wanted
- 2.9K Microsoft .NET
- 1.7K ASP.NET
- 1.1K .NET General
- 3.4K Miscellaneous
- 6 Join the Team
- 0 User Profiles
- 354 Comments on this site
- 68 Computer Emulators
- 2.1K General programming
- 187 New programming languages
- 615 Off topic board
- 179 Mobile & Wireless
- 53 Android
- 124 Palm Pilot
- 335 Multimedia
- 151 Demo programming
- 184 MP3 programming
- 0 Bash scripts
- 22 Cloud Computing
- 53 FreeBSD
- 1.7K LINUX programming
- 370 MS-DOS
- 0 Shell scripting
- 320 Windows CE & Pocket PC
- 4.1K Windows programming
- 919 Software Development
- 411 Algorithms
- 68 Object Orientation
- 89 Project Management
- 93 Quality & Testing
- 257 Security
- 7.6K WEB-Development
- 1.8K Active Server Pages
- 61 AJAX
- 2 Bootstrap Themes
- 55 CGI Development
- 19 ColdFusion
- 224 Flash development
- 1.4K HTML & WEB-Design
- 1.4K Internet Development
- 2.2K JavaScript
- 35 JQuery
- 295 WEB Servers
- 133 WEB-Services / SOAP

cHenriques
Member Posts: **1**

in Matlab

hi everyone, I'm making a research work to my master's degree and I have a problem I will appreciate your help.

I want simulate classically (not quantum mechanics) a electron orbiting a proton (H atom) immerse in a random distribution electromagnetic field. This theory is known as SED http://www.bu.edu/simulation/publications/dcole/PD...

For now we will forget the random EM. Since I'm trying first simulate the electron orbit without influence of EM.

As you now classically the electron feel the Coulomb force and the reaction term. The reaction term make the electron spin down and shock with the nucleus, the proton.

The non-relativistic Lorentz/Dirac equation to the movement is:

m*z'' = (-e^2)*z/|z|^3 + REAC - e*(EM)

The EM are the E and B random equations but for now we will forget this term.

REAC is the reaction term =>

2/3 * e^2 /c^3 * z''' (so third order differential eq)

but REAC can be approximated to

2/3 * e^2 /c^3 * d/dt( (-e^2)/m * z / z^3 ) ---> first order equation

the problem is in 2D x-y, so the vector z have two components, the z is the radius = sqrt(x^2 + y^2)

I can plot the electron orbit only with the Coulomb term:

((-e^2)*z/|z|^3)

but I cannot plot correctly with the REAC term. I already try use the ode45, the ode23 and the ode23t

I also tried use the REAC in the first form (third order ode) and in the second form

I also tried use polar coordinates, I'm lost I need help the following code is the my successful implementation of ode without the REAC

function mov = semZPF2(t,y)

r=(sqrt(y(1)^2+y(2)^2));

e = 1.602177e-19;

m=9.10938291*10^-31;

mov(1)=y(3);

mov(2)=y(4);

mov(3)=(-e^2*y(1)/r^3)/m;

mov(4)=(-e^2*y(2)/r^3)/m;

return;

------------------------------------------------------

yi=[0 5.3*10^-11]; % initial position

vi=[23.0583421525753 0]; % initial velocity

[t,y]= ode45(@(t,y) semZPF2(t,y,), [0 1E-12] , [yi vi] );

plot (y(i,1), y(i,2))

Please I need your help to implement the REAC term

m*z'' = (-e^2)*z/|z|^3 + REAC

thanks,

Carlos Henriques

I want simulate classically (not quantum mechanics) a electron orbiting a proton (H atom) immerse in a random distribution electromagnetic field. This theory is known as SED http://www.bu.edu/simulation/publications/dcole/PD...

For now we will forget the random EM. Since I'm trying first simulate the electron orbit without influence of EM.

As you now classically the electron feel the Coulomb force and the reaction term. The reaction term make the electron spin down and shock with the nucleus, the proton.

The non-relativistic Lorentz/Dirac equation to the movement is:

m*z'' = (-e^2)*z/|z|^3 + REAC - e*(EM)

The EM are the E and B random equations but for now we will forget this term.

REAC is the reaction term =>

2/3 * e^2 /c^3 * z''' (so third order differential eq)

but REAC can be approximated to

2/3 * e^2 /c^3 * d/dt( (-e^2)/m * z / z^3 ) ---> first order equation

the problem is in 2D x-y, so the vector z have two components, the z is the radius = sqrt(x^2 + y^2)

I can plot the electron orbit only with the Coulomb term:

((-e^2)*z/|z|^3)

but I cannot plot correctly with the REAC term. I already try use the ode45, the ode23 and the ode23t

I also tried use the REAC in the first form (third order ode) and in the second form

I also tried use polar coordinates, I'm lost I need help the following code is the my successful implementation of ode without the REAC

function mov = semZPF2(t,y)

r=(sqrt(y(1)^2+y(2)^2));

e = 1.602177e-19;

m=9.10938291*10^-31;

mov(1)=y(3);

mov(2)=y(4);

mov(3)=(-e^2*y(1)/r^3)/m;

mov(4)=(-e^2*y(2)/r^3)/m;

return;

------------------------------------------------------

yi=[0 5.3*10^-11]; % initial position

vi=[23.0583421525753 0]; % initial velocity

[t,y]= ode45(@(t,y) semZPF2(t,y,), [0 1E-12] , [yi vi] );

plot (y(i,1), y(i,2))

Please I need your help to implement the REAC term

m*z'' = (-e^2)*z/|z|^3 + REAC

thanks,

Carlos Henriques

Terms of use / Privacy statement / Publisher: Lars Hagelin

Programmers Heaven articles / Programmers Heaven files / Programmers Heaven uploaded content / Programmers Heaven C Sharp ebook / Operated by CommunityHeaven

© 1997-2015 Programmersheaven.com - All rights reserved.