Study interaction between intensive circularly polarized laser and hydrogen atom using a matrix method

Similar documents
Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari

A Note on Estimability in Linear Models

Lecture 23 APPLICATIONS OF FINITE ELEMENT METHOD TO SCALAR TRANSPORT PROBLEMS

External Equivalent. EE 521 Analysis of Power Systems. Chen-Ching Liu, Boeing Distinguished Professor Washington State University

COMPLEX NUMBER PAIRWISE COMPARISON AND COMPLEX NUMBER AHP

18th European Signal Processing Conference (EUSIPCO-2010) Aalborg, Denmark, August 23-27, 2010

Jones vector & matrices

The Hyperelastic material is examined in this section.

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

Radial Cataphoresis in Hg-Ar Fluorescent Lamp Discharges at High Power Density

Lucas Test is based on Euler s theorem which states that if n is any integer and a is coprime to n, then a φ(n) 1modn.

Fine Structure Calculation of Energy Levels, Oscillator Strengths and Lifetimes in Se-Like Ions

HORIZONTAL IMPEDANCE FUNCTION OF SINGLE PILE IN SOIL LAYER WITH VARIABLE PROPERTIES

CHAPTER 7d. DIFFERENTIATION AND INTEGRATION

From Structural Analysis to FEM. Dhiman Basu

Optimal Topology Design for Replaceable of Reticulated Shell Based on Sensitivity Analysis

10/7/14. Mixture Models. Comp 135 Introduction to Machine Learning and Data Mining. Maximum likelihood estimation. Mixture of Normals in 1D

ON THE COMPLEXITY OF K-STEP AND K-HOP DOMINATING SETS IN GRAPHS

CHAPTER 33: PARTICLE PHYSICS

Soft k-means Clustering. Comp 135 Machine Learning Computer Science Tufts University. Mixture Models. Mixture of Normals in 1D

VISUALIZATION OF DIFFERENTIAL GEOMETRY UDC 514.7(045) : : Eberhard Malkowsky 1, Vesna Veličković 2

Physics of Very High Frequency (VHF) Capacitively Coupled Plasma Discharges

Classical Magnetic Dipole

Rayleigh-Schrödinger Perturbation Theory

Journal of Theoretical and Applied Information Technology 10 th January Vol. 47 No JATIT & LLS. All rights reserved.

Lecture 14. Relic neutrinos Temperature at neutrino decoupling and today Effective degeneracy factor Neutrino mass limits Saha equation

Physics 256: Lecture 2. Physics

Folding of Regular CW-Complexes

Hydrogen Atom and One Electron Ions

Group Codes Define Over Dihedral Groups of Small Order

A Probabilistic Characterization of Simulation Model Uncertainties

Decentralized Adaptive Control and the Possibility of Utilization of Networked Control System

Stress-Based Finite Element Methods for Dynamics Analysis of Euler-Bernoulli Beams with Various Boundary Conditions

Grand Canonical Ensemble

PARTIAL DISTRIBUTION FUNCTION AND RADIAL STATISTICAL COEFFICIENTS FOR Be-ATOM

FEFF and Related Codes

A Novel Finite Volume Scheme with Geometric Average Method for Radiative Heat Transfer Problems *

Two-Temperature Plasma Modeling of Argon Gas Tungsten Arcs.

Relate p and T at equilibrium between two phases. An open system where a new phase may form or a new component can be added

Einstein Equations for Tetrad Fields

Lecture 3: Phasor notation, Transfer Functions. Context

Addition of angular momentum

1- Summary of Kinetic Theory of Gases

Electrochemical Equilibrium Electromotive Force. Relation between chemical and electric driving forces

Discrete Shells Simulation

Numerical Study of Two-fluid Flowing Equilibria of Helicity-driven Spherical Torus Plasmas

Phys 774: Nonlinear Spectroscopy: SHG and Raman Scattering

A NEW GENERALISATION OF SAM-SOLAI S MULTIVARIATE ADDITIVE GAMMA DISTRIBUTION*

Representation and Reasoning with Uncertain Temporal Relations

Addition of angular momentum

Outlier-tolerant parameter estimation

Derangements and Applications

Economics 600: August, 2007 Dynamic Part: Problem Set 5. Problems on Differential Equations and Continuous Time Optimization

Chapter 6 Student Lecture Notes 6-1

Naresuan University Journal: Science and Technology 2018; (26)1

Quasi-Classical States of the Simple Harmonic Oscillator

Electrostatic Surface Waves on Semi-Bounded Quantum Electron-Hole Semiconductor Plasmas

A general N-dimensional vector consists of N values. They can be arranged as a column or a row and can be real or complex.

167 T componnt oftforc on atom B can b drvd as: F B =, E =,K (, ) (.2) wr w av usd 2 = ( ) =2 (.3) T scond drvatv: 2 E = K (, ) = K (1, ) + 3 (.4).2.2

Three-Node Euler-Bernoulli Beam Element Based on Positional FEM

Geometric algebra, qubits, geometric evolution, and all that

Modeling and Simulation Analysis of Power Frequency Electric Field of UHV AC Transmission Line

PH300 Modern Physics SP11 Final Essay. Up Next: Periodic Table Molecular Bonding

SCITECH Volume 5, Issue 1 RESEARCH ORGANISATION November 17, 2015

Linear Algebra Provides a Basis for Elasticity without Stress or Strain

INTERFACE CORNERS IN ANISOTROPIC/PIEZOELECTRIC/ VISCOELASTIC MATERIALS

Review - Probabilistic Classification

Analyzing Frequencies

Homotopy perturbation technique

PREDICTION OF STRESS CONCENTRATION FACTORS IN UNLAPPED SQUARE HOLLOW "K" JOINTS BY THE FINITE ELEMENT METHOD

JEE-2017 : Advanced Paper 2 Answers and Explanations

Modelling of new generation plasma optical devices

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding...

Polytropic Process. A polytropic process is a quasiequilibrium process described by

EDGE PEDESTAL STRUCTURE AND TRANSPORT INTERPRETATION (In the absence of or in between ELMs)

As the matrix of operator B is Hermitian so its eigenvalues must be real. It only remains to diagonalize the minor M 11 of matrix B.

September 27, Introduction to Ordinary Differential Equations. ME 501A Seminar in Engineering Analysis Page 1. Outline

Davisson Germer experiment Announcements:

An Effective Technique for Enhancing Anti-Interference Performance of Adaptive Virtual Antenna Array

Α complete processing methodology for 3D monitoring using GNSS receivers

The Fourier Transform

ANALYTICITY THEOREM FOR FRACTIONAL LAPLACE TRANSFORM

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

ST 524 NCSU - Fall 2008 One way Analysis of variance Variances not homogeneous

Poisson-Boltzmann Theory with Non-Linear Ion Correlations

te Finance (4th Edition), July 2017.

Study of Dynamic Aperture for PETRA III Ring K. Balewski, W. Brefeld, W. Decking, Y. Li DESY

Schrodinger Equation in 3-d

On the Hamiltonian of a Multi-Electron Atom

Observer Bias and Reliability By Xunchi Pu

Determination of Optimum Length of Frames and Number of Vectors to Compress ECG Signals

Elements of Statistical Thermodynamics

Phy213: General Physics III 4/10/2008 Chapter 22 Worksheet 1. d = 0.1 m

Application of Local Influence Diagnostics to the Linear Logistic Regression Models

PHYSICS 489/1489 LECTURE 7: QUANTUM ELECTRODYNAMICS

VLSI Implementation and Performance Evaluation of Low Pass Cascade & Linear Phase FIR Filter

THREE DIMENSIONAL GEOMETRY MAINTENANCE FOR FORMATION FLYING ON ELLIPTIC ORBITS

NON-SYMMETRY POWER IN THREE-PHASE SYSTEMS

Fundamentals of Continuum Mechanics. Seoul National University Graphics & Media Lab

Transcription:

ISBN 978-1-84626-020-9 Procdngs of 3 rd Intrnatonal Workshop on Matrx Analyss angzhou,p.r.chna.july 9-13, 2009, pp. 199-202 ( Wll st y th pulshr ) Study ntracton twn ntnsv crcularly polarzd lasr and hydrogn atom usng a matrx mthod L.F. Yang 1, 3, W. Xu 1+, Q. R. Zhang 2, X. L. Ca 1, 3, G. T. Fan 1, 3, G. W. Fan 1,3, W. Luo 1,3, Y. J. L 1,3, Q. Y. Pan 1, B. J. Xu 1, Y. Xu 1,3, Z. Yan 1 1 Shangha Insttut of Appld Physcs, Chns Acadmy of Scncs, 201800, Shangha, Chna 2 Dpartmnts of Tchncal Physcs, Pkng Unvrsty, 100871, Bjng, Chna 3 Graduat School of th Chns Acadmy of Scncs, 100039, Bjng, Chna Astract. Rcntly th study of ntnsv lasr and atom ntracton s on of th most ntrstng topcs n lasr and atom ntrdscplnary physcs. Study of ntracton of hydrogn atom and ntns lasr usng a matrx mthod s rportd n ths papr. Ths mthod ncluds Gvns rducton for ral symmtrc trdagonal matrx form of amltonan matrx, Thomas algorthm for gnvctors, and scton algorthm wth Sturm thorm for gnvalus. Som prlmnary rsults and progrss for ths study ar prsntd. Kywords: ntnsv lasr, hydrogn atom, matrx mthod. 1. Introducton In rcnt dcads, wth dvlopmnt of hgh ntnsty lasr, th ntnsv lasr and mattr ntracton hav com a hot topc n oth xprmntal and thortcal physcs [1-2]. Bcaus th hydrogn s th smplst systm for studyng th ntracton of an atom systm wth a lasr fld, wthout ntrfrnc from multlctron ffcts, th study of ntracton twn hydrogn atom and a lasr provds a clan tst cas for thortcal calculaton.thrfor, w dvlop an approach to solv th tm dpndnt Schrödngr quaton (TDSE) of hydrogn atom n th prsnc of crcularly polarzd lasr fld [3]. Th rmanng of ths papr s organzd as follows: Scton 2 dscrs our thortcal approach. In Scton 3 th corrspondng gnvalus and gnvctors ar calculatd and chckd. Scton 4 concluds th papr and dscusss possl futur rsarch drctons. 2. Thortcal approach Consdrng th hydrogn atom n th prsnc of a rlatv long puls and crcularly polarzd lasr fld, Ths work s supportd y th On undrd Talnts Projct of Chns Acadmy of Scncs (2006) (26010701), y Knowldg Innovaton Projct of Chns Acadmy of Scncs (KJCX2-SW-N13), y Pujang Talnt Projct of th Shangha Scnc and Tchnology Commtt (06PJ14114), and y th Natonal Natural Scnc Foundaton of Chna (10675156). + E-mal addrss: xuwang@snap.ac.cn

th TDSE of th lctron can changd nto a non tm dpndnt Schrödngr quaton (NTDSE) y som rotaton and untary transformaton [4] = E (1), whr E s Egn nrgy of gnvctor and 2 2 2 2 2 = p + V() r k pzl k A A z + lz + + px + kcl. (2) z 2μ μ 2μ 2μ 2μ r, Vr () s coulom potntal, c s lght vlocty, μ s th mass of lctron, k s wav numr, p s momntum, px and p y ar x and y drcton componnts of p rspctvly, l z = m s th z-componnt of th angular momntum, and m s th magntc quantum numr. At frst, w only consdr th part of ound statonary stat of hydrogn atom ψ ( knlm ( ) or to ass to comn th gnvctor of, and n th coordnaton t can wrttn as = C = C ψ = C R r Y θϕ, ( n l) ( ) ( ) ( ) nl lm l l+ 3/2 r +! 2 l 2 n C r 1F1 l 1 n,2l 2, r Ylm ( θϕ, ) 2l 1! 2n n l 1! n n = + + ( + ), whr C s th coffcnt of th orthogonal ass vctor ψ. Thn usng th matrx form of (NTDSE), scular Eq. 1 can wrttn as 11 12 1 f c 1 c 1 21 22 2 f c2 c 2 = ε (4) c f 1 f 2 ff f c f, n whch, = k ( nlm), (5) k k k C s arvaton of C, ε s th gnvalu of ffctv amltonan for gnvctor, and can spcfcally rwrttn as 2 2 2 2 2 2 p k m A A k k km k k = () k k + V r + + + kcm δ + Ψ kk px Ψ Ψ pz Ψ. (6) 2μ 2μ 2μ 2μ μ Fg. 2.1 shows th procss of our matrx mthod. Aftr w calculat th lmnt of ffctv amltonan th ral symmtrc trdagonal matrx form of amltonan matrx s otand y Gvns rducton [5]. Th rlatd gnvalus ar thn calculatd y scton algorthm wth Sturm thorm [6], and gnvctors ar computd y Thomas algorthm [7] wth wll-sparatd gnvalu. Fnally, w can otan som physcs quantts, such as onzaton cross scton and lasr ntnsty dpndnc [8]. In ths papr, w wll only prsnt som rsults of gnvalu and gnvctor of. k Ψ ) (3)

0.71235 0.33973 0.28615 0.30388 0.29401 0.33973 1.18585-0.21846 0.06685-0.37360 0.28615-0.21846 0.18159 0.27955 0.38898 0.30388 0.06685 0.27955 0.23195 0.20496 0.29401-0.37360 0.38898 0.20496 0.46004 Fg. 2.1. Flowchart of our approach. Tal 3.1. A fv ordr squar matrx M0. 3. Som rsults Bfor gvng som rsults of gnvalu and gnvctor of, w would lk to chck th rlalty and accuracy for Gvns rducton and Bscton algorthm usd n our approach. W compar th gnvalus and gnvctors from our approach to a known rsult. Th matrx usd for comparson s a fv ordr squar matrx M0 [5] (shown n Tal 3.1). Tal 3.2 Egnvalu. Egnvalu No No.1 No.2 No.3 No.4 No.5 Matrx calculator -0.141269V 0.050511V 0.092036V 0.092036V 0.092036V Gvns+Bscton* -0.141269V 0.050511V 0.092036V 0.092036V 0.092036V Tal 3.3 No1. Egnvctor. Mthod Gvns+ Matrx Thomas calculator ndx No.1 No.1 1-0.000000-3.316883E-07 2 0.004647 0.004647 3 0.826421 0.826421 4-0.401293-0.401293 5-0.394932-0.394932 Tal 3.4 No3. Egnvctor. Mthod Gvns+ Matrx Thomas Calculator ndx No.3 No.3 1-0.000001-8.184097E-07 2 0.161946 0.161946 3-0.057399-0.057399 4-0.747648-0.747648 5 0.641484 0.641484 Excpt No.1 and No.3 gnvctors, othr gnvctors ar compltly sam. Only No.1 and No3 gnvctors wth som slght dffrncs ar shown n Tal 3.3 and Tal 3.4, rspctvly. Thrfor, w can asly conclud that th rsult of our modl has a 5-dgt prcson for gnvctor and has mor than 5- dgt prcson for gnvalu. 5 7 6 As w st lasr paramtrs A = 1.8 10 V s/ m, and k = 1.03 10 / m (wav lngthy λ = 1.64 10 m ), th gnvalu E of amltonan, and corrspondng gnvctor n th prsnc of ths lasr fld ar calculatd asd on our approach. Tal 3.5 shows th frst tn tms of gnvctor ψ = 4.65 of E V

amltonan (whch has a largst componnt of ground stat of atomc hydrogn) asd on th ass of ound stats of hydrogn atom, and ts corrspondng gn nrgy E= 4.65V. Tal 3.5. Data for th frst tn tms of gnvctorψ = 4.65. Indx of c Valu of c Indx of c Valu of c 1 0.7570579446137287 6 0.00208548190289008 2 0.09483098226226934 7-0.17854082543374933 3-0.3767373556802446 8-4.97999186909243E-7 4-0.0000012743403275 9 0.1615973380163324 5 0.33683326841968736 10 0.11358396205104035 ψ = Usng data n Tal 3.5, on can plot E 4.65 mathmatc softwar. Z 7, Ψ E 4.65 V 2 y 10 V E, an gnvctor of ffctv amltonan 10 V Z 7, Ψ E 13.6 V 2 y wth 5 5 0 x 0 x 5 5 10 10 5 0 5 10 10 10 5 0 5 10 Fg. 3.6. Proalty forψ E = 4.65V, an gnvctor of ffctv amltonan (lft pctur), and th proalty for ground stat of atomc hydrogn wth lasr off (rght pctur). Wht color rprsnts a larg proalty and lu rprsnts a low proalty n oth pcturs. ψ = Bcaus E 4.65V has a largst componnt of ground stat of atomc hydrogn wthout lasr among all gnvctors of, th ground stat ψ E = 13.6V of atomc hydrogn wll most lkly transt to gnvctor 5 7 ψ E = 4.65V of ffctv amltonan as a lasr fld ( A = 1.8 10 V s/ m, and k = 1.03 10 / m) s mployd. Th transton proalty s calculatd as 0.57 y our approach. owvr, from oth plots n Fg. 3.6, on can clarly s that th sphrcal symmtry of th ground stat of atomc hydrogn s rokn as th crcularly polarzd lasr fld s mployd.

4. Concluson and outlook In ths papr, w hav rfly dscrd th procss of how to stalsh th approach, and chckd ts rlalty and accuracy wth an xampl whn th gnvalus ar wll-sparatd. On of th gnvctor of wth a largst componnt of ground stat of hydrogn atom has n studd. In th futur, w wll study th gnvctors whn th corrspondng gnvalus ar not wll-sparatd, onzaton cross scton, transton rat twn two dscrt stats, and spcally th ffct of contnuum part of hydrogn atom wav functon on our approach n th prsnc of th crcularly polarzd lasr fld. 5. Rfrncs [1] M. Protopapas, C.. Ktl and P.L. Knght, Atomc physcs wth supr-hgh ntnsty lasrs, Rp. Prog. Phys, 60(1997), pp 389 486. [2] C.J. Joachan, M. Dörr, and N.J. Kylstra, gh-ntnsty lasr-atom physcs, Adv. At, Mol, and Opt. Phy, pp.42-225, 2000. [3] Q.R. Zhang, Intnsty dpndnt of th th photolctrc ffct nducd y crcularly polarzd lasr am, Phys. Ltt. A, pp.216-125, 1996. [4] D. Zhang, Q.R. Zhang, Th Photoonzaton of Atomc ydrogn n an Intns Crcularly Polarzd Lasr Fld, Comm.Thor.Phys, pp. 685-690,2001. [5] J.. Wlknson, Th Algrac Egnvalu Prolm, Oxford Unvrsty Prss,1965. [6] W. Barth, R. S. Martn, and J.. Wlknson, Calculaton of th Egnvalus of a Symmtrc Trdagonal Matrx y th Mthod of Bscton, Numrsch Mathmatk 9 (1967), pp. 386-393. [7] S.D. Cont, C.W. D Boor, Elmntary Numrcal Analyss. McGraw-ll Educaton, 1972. [8] Q.R. Zhang, Non-ntgr Quantum Transton, a Tru Non-prturaton Ect n Lasr-Atom Intracton Commun. Thor. Phys, 47 (2007), pp. 1017 1023.