Tomography of Highly Scattering Media with the Method of Rotated Reference Frames. Vadim Markel Manabu Machida George Panasyuk John Schotland
|
|
- Morris Barrett
- 6 years ago
- Views:
Transcription
1 Tomography of Highly Scattering Media with the Method of Rotated Reference Frames Vadim Markel Manabu Machida George Panasyuk John Schotland Radiology/Bioengineering UPenn, Philadelphia
2 MOTIVATION (The Inverse Problem Perspective)
3 ρ s S zs α(r) zd z Given a data function φ( ρ, ρ ) s d which is measured for multiple pairs ( ρ, ρ ), find the absorption s d D ρ d coefficient α( r) inside the slab 3
4 Linearized Integral Equation φ( ρ, ρ ) = Γ( ρ, ρ ; r) ( r)d 3 s d s d δα r φ ( ρ, ρ ) s d = I ( ρs, zs; ρd, zd) I0 ( ρs, zs; ρd, zd) I ( ρ, z ; ρ, z ) (measurable data-function) 0 α ( r) = α + δα( r) 0 s s d d Γ ( ρ, ρ ) = G ( ρ, z ; r) G ( r; ρ, z ) s d 0 s s 0 d d (first Born approximation)
5 How To Invert? Numerical SVD approach: φ = Γ δα i ij j j j φ i φ =Γδα Γ δα φ = Γ Γδα * * + * 1 * ( ) = ΓΓ Γ reg φ
6 Analytical SVD approach: Making use of the translational invariance φ ( q, q ) = φ ( ρ, ρ ) e d ρ d ρ i ( q s ρ s+ q d ρ d ) s d s d s d qs = q / + p, qd = q / p; Data function: ψ ( q, p) = φ ( q / + p, q / p) L ψ ( q, p) = g ( q / + p; z) g ( q / p; z) δα ( q; z)dz 0 iρ q δα ( q; ) = δα ( ρ, ) d ρ z dz d q G ( ρ, z ; ρ, z) = g s z e ( π ) 0 s s s ( π ) d ( q; z) e d q G0 ( ρ, z; ρd, zd) = g ( ; ) d q z e iq ( ρ ρ ) s iq ( ρ ρ) d 6
7 6cm S.D.Konecky, G.Y.Panasyuk, K.Lee, V.Markel, A.G.Yodh and J.C.Schotland Imaging complex structures with diffuse light Optics Express 16(7), (008)
8 In the case of RTE, we need the plane-wave decomposition of the Greens function, of the form d q G (,; rsr ˆ, sˆ ) = g(;,; qz sˆ z, sˆ ) e ( π ) 0 iq ( ρ ρ ) and the integral kernel Γ ( q, p; z) = g( q/ + p; z, zˆ ; z, sˆ ) g( q/ p; z,; sˆ z, zˆ )ds We then get the integral equations of the form z d ψ( qp, ) = Γ( qp, ; z) δα ( q; z)dz z s s d
9 Motivation Continued (The Forward Model Perspective)
10 Plane Wave Modes The Weyl Expanson The Helmholtz Equation ( + k ) u= 0 0 k = const e ikr k k = k 0 e r ik r 0 i = π Q= k q 0 Q e d q 1 i( q ρ+ Qz) The Diffusion Equation ( D + α) u= 0 D, α = const e k kr k = k = α 0 / D e r kr 0 1 = π Q= k + q 0 Q e d q 1 iq ρ Qz RTE ( sˆ + μ ) I( r, sˆ) = t Arnold Kim and us (the past ~10 years) Arnold Kim and us (the past ~10 years) = μ ( ˆˆ, ) (, ˆ s A ss I rs ) μ, μ = const t s
11 Spectral Method for the RTE? sˆ μ ˆ ˆ ˆ ˆ ˆ ˆ t I r s μs A s s I r s s ε r s RTE: ( + )(,) = (, )(, )d + (,) Where is the "spectral variable"? How can we write this equation in the form ( z+ W) I = ε? μ and μ do not qualify... t s We can try to expand I( rs, ˆ) into a 3D Fourier integral with respect to r and into the basis of ordinary spherical harmonics Y ( θϕ, ) with respect to sˆ... lm...and see if the equation can be cast into the desired form.
12 The Conventional Method of Spherical Harmonics This results in the following system of equations with respect to the vector of the expansion coefficients I( k) ( k - the Fourier variable) : ia k I( k) + ia k I( k) + ia k I( k) + D I( k) = ε ( k) ( x) ( y) ( z) x y z A, A, ( x) ( y) A ( z) are different matrices. ( x) * Alm, l m = sinθcos ϕylm( θϕ, ) Yl m ( θϕ, )sin θdθdϕ, etc.... This rather awe-inspiring set of equations has perhaps only academic interest. K.M.Case, P.F.Zweifel, Linear Transport Theory
13 Rotated Reference Frames The usual sperical harmonics are defined in the laboratory reference frame. Then θ and ϕ are the polar angles of the unit vector sˆ in that frame. THE MAIN IDEA: For each value of the Fourier variable k, use spherical harmonics defined in a reference frame whose z-axis is aligned with the direction of k. We call such frames "rotated". Spherical harmonics defined in the rotaded frame are denoted by Y ( sk ˆ; ˆ ).
14 Rotation of the Laboratory Frame (x,y,z). z y' z' l ˆ ˆ l Y (; ) (,,0) () ˆ lm sk = Dm m ϕk θk Ylm s m = l θ k y Wigner D-functions Euler angles x ϕ θ Spherical functions in the laboratory frame x'
15 ( sˆ + μ ) I = μ AI + ε t (,) ˆ (,) ˆ d ikr 3 I rs = I kse k ikr ε(,) rsˆ = ε( ks,) ˆ e d ( ik sˆ + μ ) I = μ AI + ε t I ( ks, ˆ) = F ( k) Y ( sk ˆ; ˆ ) lm, ε ( ks, ˆ) = E ( k) Y ( sk ˆ; ˆ ) lm, lm lm A(, ss) AY (; skˆ) Y (; skˆ) * ˆˆ = ˆ ˆ l lm lm lm, s s lm lm 3 k
16 lm lm l m l lm lm ik R F ( k) + σ F ( k) = E ( k) lm σ μ μ = + (1 A ) l a s l R = ˆ Y ˆ Y ˆ s= lm sˆ k * (; sk ˆ ) (; sk ˆ )d lm lm l m [ b( m) b ( m) ] = + δ δ δ mm l l = l 1 l+ 1 l = l+ 1 b( m) l = l m 4l 1
17 RTE in the Angular Basis of Rotated Spherical Functions ikr I( k) + D I( k) = E( k) Scalar spectral parameter Block-tridiagonal real symmetric matrix Source term Diagonal matrix Slm, l m = δ ll δmm [ μa + μs(1 Al )] Parameters of the phase function. (For the HG l model, A = g ) l
18 Let D W = SS = S RS 1 1 ikr I( k) + D I( k) = E( k) ikw S I k S E 1 ( + 1) ( ) = ( ) I k S ikw S E ( ) = (1 + ) ( ) k k W ψ = λ ψ μ μ μ I ( k) = μ 1 1 S ψ ψ μ μ S E 1+ ikλ μ ( k)
19 The Spectral Solution ˆ ˆ (,) d σ 1+ ikλ 1 1 lm ψμ ψμ S E( k) Ylm( s; k) ikr 3 I rsˆ = e k lm ε(,) rsˆ = δ( r r) δ( sˆ sˆ ) ikr 0 * ˆ E ˆ lm( k) = e Y ( 3 lm s0; k) ( π ) l μ μ The integral is not easy but doable. Details in J.Phys.A 39, 115 (006)
20 ˆ m * ˆ ˆ ˆ ˆ = lm χll l m 0 m= l, l = m G (,; rsr, s) Y (; sr) ( R) Y ( s; Rˆ) m m ( 1) χll ( R) = ( 1) π σ σ l l l M n 3 M= l n, λ > 0 λn l φ ( M) φ ( M) l ( M ) l l l + j,0 l l + j,0 R ClMl,,, MClml,,, m kl l + j j= 0 λn( M ) l = min( l, l ) μ = ( Mn, ) lm ψ ( ) Mn = δmm l φn M n n r ŝ ŝ 0 R r 0
21 The Plane-Wave Decomposition ˆ ˆ (,) d σ 1+ ikλ 1 1 lm ψμ ψμ S E( k) Ylm( s; k) ikr 3 I rsˆ = e k dq ( 0 ) * ˆ ˆ iq ρ ρ lm ˆ ˆ ˆ lm lm 0 l m 0 lm, l m G (,; rsr, s) = e Y (;) szg (;, qz z ) Y ( s;) zˆ π ( ) dk g q z z ˆ ( ) ˆ 0 = e D q+ zk K q + k D q+ zk π z ikz ( z z0 ) ( ;, ) ( ) ( ) Kw ( ) = lm μ S 1 l ψ μ ψ 1+ iwλ S μ z z z 1 μ μ iϕ J iθ J μ ; D ( kˆ ) = e cos θ = k / q + k k z k z e k y z
22 i( m m ) ϕq l l lm e g ( q; z, z ) = sgn( z z ) d [ iτ ( qλ )] lm Q ( q) = q + 1/ λ μ 0 0 σσ l l μ cos[ τ( )] 1, sin[ τ( )] l+ l + m+ m l [ ] mm1 m = l m = l μλ, > 0 1 Qμ ( q) z z0 e l 1ψ μ ψ μ mm τ λ μ λμqμ( q) lm l m d [ i ( q )] i x = + x i x = ix α D ( αβγ,, ) = e d ( β) e l im l im γ mm mm μ μ
23 Plane waves: Evanescent Waves kr I = e F sˆ k k = k = knˆ k k( ), 1/ λμ, nˆ - real unit vector F ( ˆ) = lm Y ( ˆ; ˆ k s ψ μ lm s k) lm Evanescent waves: k = q± zˆ + q zˆ = k k = i q 1/ λμ, 0, 1/ λμ iq k =± q + 1/ λ z z n
24 kˆ r exp( imϕ k ) l Iˆ = exp Y ( ˆ; ˆ lm ) dmm ( θ ) l φn ( M) Mn sz k k λmn lm σ l (plane wave modes) exp( imϕ ) I = exp i Q ( q) z Y ( sz ˆ; ˆ) ( 1) d ( θ ) l φ ( M) [ q ρ ] ( ± ) q l+ m l qmn Mn lm m, M k n lm σ l (evanescent modes) G V 0 dq ( ) ( ) (,; ˆ 0, ˆ ˆ ˆ 0) I ± rsr s = qmn(,) rs V Mn I q qmn( r0, s0) ( π ) qmn = Q Mn 1 ( q) λ Mn Mn
25 Half-space problem Z>0 Scattering medium Propagationd escribed by the RTE z<0 Nonscattering medium d q I = A ( q) I ( ) Mn π Mn ( + ) qmn
26
27
28
29 RESULTS: SIMULATIONS
30 A set of 5 point absorbers in an L=6l* slab The field of view is 16l* RTE Diffusion approximation
31 A bar target in the center of the same slab RTE Diffusion approximation
32 RESULTS: EXPERIMENT
33 Two thin vertical wires in a 1cm thick slab filled with intralipid solution (looks like milk) RTE Diffusion approximation
34 CONCLUSIONS The method of rotated reference frames can be used in optical tomography of mesoscopic samples The images are of superior quality compared to those obtained by using the diffusion approximation The quality of images can be comparable to that in X-ray tomography because the RTE retains some information about ballistic rays, single-scattered rays, etc.
35 1. V.A.Markel, Modified spherical harmonics method for solving the radiative transport equation, Letter to the Editor, Waves in Random Media 14(1), L13-L19 (004).. G.Y.Panasyuk, J.C.Schotland, and V.A.Markel, Radiative transport equation in rotated reference frames, Journal of Physics A, 39(1), (006). 3. J.C.Schotland and V.A.Markel, Fourier-Laplace structure of the inverse scattering problem for the radiative transport equation, Inverse Problems and Imaging 1(1), (007). 4. M.Machida, J.C.Schotland and V.A.Markel, MRRF with boundary conditios, under consideration in J.Phys.A 5. M.Machida, J.C.Schotland and V.A.Markel, MRRF with boundary conditios, under consideration in J.Phys.A 6. M.Machida, J.C.Schotland and V.A.Markel, Inverse problem with MRRF, to be submitted to Inv.Prob. in the near future..
Method of Rotated Reference Frames
Method of Rotated Reference Frames Vadim Markel in collaboration with George Panasyuk Manabu Machida John Schotland Radiology/Bioengineering UPenn, Philadelphia Outline MOTIVATION (the inverse problems
More informationRadiative transport equation in rotated reference frames
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 39 (26) 115 137 doi:1.188/35-447/39/1/9 Radiative transport equation in rotated reference frames George
More informationGreen's functions for the two-dimensional radiative transfer equation in bounded media
Home Search Collections Journals About Contact us My IOPscience Green's functions for the two-dimensional radiative transfer equation in bounded media This article has been downloaded from IOPscience Please
More informationReconstructing a Function from its V-line Radon Transform in a D
Reconstructing a Function from its V-line Radon Transform in a Disc University of Texas at Arlington Inverse Problems and Application Linköping, Sweden April 5, 2013 Acknowledgements The talk is based
More informationSolutions to exam : 1FA352 Quantum Mechanics 10 hp 1
Solutions to exam 6--6: FA35 Quantum Mechanics hp Problem (4 p): (a) Define the concept of unitary operator and show that the operator e ipa/ is unitary (p is the momentum operator in one dimension) (b)
More informationRadiative transport and optical tomography with large data sets
Radiative transport and optical tomography with large data sets Manabu Machida, George Panasyuk, Zheng-Min Wang, Vadim Markel, John Schotland To cite this version: Manabu Machida, George Panasyuk, Zheng-Min
More informationFINITE SELF MASS OF THE ELECTRON II
FINITE SELF MASS OF THE ELECTRON II ABSTRACT The charge of the electron is postulated to be distributed or extended in space. The differential of the electron charge is set equal to a function of electron
More informationSurface Sensitive X-ray Scattering
Surface Sensitive X-ray Scattering Introduction Concepts of surfaces Scattering (Born approximation) Crystal Truncation Rods The basic idea How to calculate Examples Reflectivity In Born approximation
More informationLecture 4. Diffusing photons and superradiance in cold gases
Lecture 4 Diffusing photons and superradiance in cold gases Model of disorder-elastic mean free path and group velocity. Dicke states- Super- and sub-radiance. Scattering properties of Dicke states. Multiple
More information3. Calculating Electrostatic Potential
3. Calculating Electrostatic Potential 3. Laplace s Equation 3. The Method of Images 3.3 Separation of Variables 3.4 Multipole Expansion 3.. Introduction The primary task of electrostatics is to study
More information6. LIGHT SCATTERING 6.1 The first Born approximation
6. LIGHT SCATTERING 6.1 The first Born approximation In many situations, light interacts with inhomogeneous systems, in which case the generic light-matter interaction process is referred to as scattering
More informationRepresentation theory and quantum mechanics tutorial Spin and the hydrogen atom
Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition
More informationPreliminary Examination - Day 1 Thursday, August 10, 2017
UNL - Department of Physics and Astronomy Preliminary Examination - Day Thursday, August, 7 This test covers the topics of Quantum Mechanics (Topic ) and Electrodynamics (Topic ). Each topic has 4 A questions
More informationSpherical Harmonics on S 2
7 August 00 1 Spherical Harmonics on 1 The Laplace-Beltrami Operator In what follows, we describe points on using the parametrization x = cos ϕ sin θ, y = sin ϕ sin θ, z = cos θ, where θ is the colatitude
More informationAppendix A Spin-Weighted Spherical Harmonic Function
Appendix A Spin-Weighted Spherical Harmonic Function Here, we review the properties of the spin-weighted spherical harmonic function. In the past, this was mainly applied to the analysis of the gravitational
More informationQuantum Mechanics II Lecture 11 (www.sp.phy.cam.ac.uk/~dar11/pdf) David Ritchie
Quantum Mechanics II Lecture (www.sp.phy.cam.ac.u/~dar/pdf) David Ritchie Michaelmas. So far we have found solutions to Section 4:Transitions Ĥ ψ Eψ Solutions stationary states time dependence with time
More informationQuantum Mechanics II
Quantum Mechanics II Prof. Boris Altshuler April, 20 Lecture 2. Scattering Theory Reviewed Remember what we have covered so far for scattering theory. We separated the Hamiltonian similar to the way in
More informationAkira Ishimaru, Sermsak Jaruwatanadilok and Yasuo Kuga
INSTITUTE OF PHYSICS PUBLISHING Waves Random Media 4 (4) 499 5 WAVES IN RANDOMMEDIA PII: S959-774(4)789- Multiple scattering effects on the radar cross section (RCS) of objects in a random medium including
More informationFun With Carbon Monoxide. p. 1/2
Fun With Carbon Monoxide p. 1/2 p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results C V (J/K-mole) 35 30 25
More informationLet b be the distance of closest approach between the trajectory of the center of the moving ball and the center of the stationary one.
Scattering Classical model As a model for the classical approach to collision, consider the case of a billiard ball colliding with a stationary one. The scattering direction quite clearly depends rather
More informationEffects of Entanglement during Inflation on Cosmological Observables
Effects of Entanglement during Inflation on Cosmological Observables Nadia Bolis 1 Andreas Albrecht 1 Rich Holman 2 1 University of California Davis 2 Carnegie Mellon University September 5, 2015 Inflation
More informationOptimized diffusion approximation
356 Vol. 35, No. 2 / February 208 / Journal of the Optical Society of America A Research Article Optimized diffusion approximation UGO TRICOLI, CALLUM M. MACDONALD, ANABELA DA SILVA, AND VADIM A. MARKEL*
More informationOn Impossibility of Negative Refraction
On Impossibility of Negative Refraction Vadim A. Markel Radiology/Bioengeneering, UPenn, Philly REFERENCES: V.A.Markel, Correct definition of the Poynting vector in electrically and magnetically polarizable
More informationGroup Representations
Group Representations Alex Alemi November 5, 2012 Group Theory You ve been using it this whole time. Things I hope to cover And Introduction to Groups Representation theory Crystallagraphic Groups Continuous
More informationdoes not change the dynamics of the system, i.e. that it leaves the Schrödinger equation invariant,
FYST5 Quantum Mechanics II 9..212 1. intermediate eam (1. välikoe): 4 problems, 4 hours 1. As you remember, the Hamilton operator for a charged particle interacting with an electromagentic field can be
More informationThe Radiative Transfer Equation
The Radiative Transfer Equation R. Wordsworth April 11, 215 1 Objectives Derive the general RTE equation Derive the atmospheric 1D horizontally homogenous RTE equation Look at heating/cooling rates in
More informationFORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationNotes 19 Gradient and Laplacian
ECE 3318 Applied Electricity and Magnetism Spring 218 Prof. David R. Jackson Dept. of ECE Notes 19 Gradient and Laplacian 1 Gradient Φ ( x, y, z) =scalar function Φ Φ Φ grad Φ xˆ + yˆ + zˆ x y z We can
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 18 Properties of Sound Sound Waves Requires medium for propagation Mainly
More informationOptical Imaging Chapter 5 Light Scattering
Optical Imaging Chapter 5 Light Scattering Gabriel Popescu University of Illinois at Urbana-Champaign Beckman Institute Quantitative Light Imaging Laboratory http://light.ece.uiuc.edu Principles of Optical
More informationGeneralized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields
PHYSICAL REVIEW E 71, 5661 5 Generalized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields D. R. Lytle II Department of Electrical and Computer Engineering,
More informationCold atoms in the presence of disorder and interactions
Cold atoms in the presence of disorder and interactions Collaboration: A. Minguzzi, S. Skipetrov, B. van-tiggelen (Grenoble), P. Henseler (Bonn), J. Chalker (Oxford), L. Beilin, E. Gurevich (Technion).
More informationPhys 622 Problems Chapter 6
1 Problem 1 Elastic scattering Phys 622 Problems Chapter 6 A heavy scatterer interacts with a fast electron with a potential V (r) = V e r/r. (a) Find the differential cross section dσ dω = f(θ) 2 in the
More informationBroken-Ray and Star Transforms and their Inversion
Broen-Ray and Star Transforms and their Inversion Vadim A Marel, Institut Fresnel Joint wor with Fan Zhao (f. @ UPenn) John Schotland (f. @ Upenn, now @UMich) Lucia Florescu (f. @ UPenn) Interest in using
More informationSummary: angular momentum derivation
Summary: angular momentum derivation L = r p L x = yp z zp y, etc. [x, p y ] = 0, etc. (-) (-) (-3) Angular momentum commutation relations [L x, L y ] = i hl z (-4) [L i, L j ] = i hɛ ijk L k (-5) Levi-Civita
More informationComputation of the scattering amplitude in the spheroidal coordinates
Computation of the scattering amplitude in the spheroidal coordinates Takuya MINE Kyoto Institute of Technology 12 October 2015 Lab Seminar at Kochi University of Technology Takuya MINE (KIT) Spheroidal
More informationLecture 6 Scattering theory Partial Wave Analysis. SS2011: Introduction to Nuclear and Particle Physics, Part 2 2
Lecture 6 Scattering theory Partial Wave Analysis SS2011: Introduction to Nuclear and Particle Physics, Part 2 2 1 The Born approximation for the differential cross section is valid if the interaction
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More information( ) ( ) QM A1. The operator ˆR is defined by R ˆ ψ( x) = Re[ ψ( x)] ). Is ˆR a linear operator? Explain. (it returns the real part of ψ ( x) SOLUTION
QM A The operator ˆR is defined by R ˆ ψ( x) = Re[ ψ( x)] (it returns the real part of ψ ( x) ). Is ˆR a linear operator? Explain. SOLUTION ˆR is not linear. It s easy to find a counterexample against
More informationMath review. Math review
Math review 1 Math review 3 1 series approximations 3 Taylor s Theorem 3 Binomial approximation 3 sin(x), for x in radians and x close to zero 4 cos(x), for x in radians and x close to zero 5 2 some geometry
More informationI. Introduction. A New Method for Inverting Integrals Attenuated Radon Transform (SPECT) D to N map for Moving Boundary Value Problems
I. Introduction A New Method for Inverting Integrals Attenuated Radon Transform (SPECT) D to N map for Moving Boundary Value Problems F(k) = T 0 e k2 t+ikl(t) f (t)dt, k C. Integrable Nonlinear PDEs in
More informationHigher dimensional Kerr-Schild spacetimes 1
Higher dimensional Kerr-Schild spacetimes 1 Marcello Ortaggio Institute of Mathematics Academy of Sciences of the Czech Republic Bremen August 2008 1 Joint work with V. Pravda and A. Pravdová, arxiv:0808.2165
More informationI. Rayleigh Scattering. EE Lecture 4. II. Dipole interpretation
I. Rayleigh Scattering 1. Rayleigh scattering 2. Dipole interpretation 3. Cross sections 4. Other approximations EE 816 - Lecture 4 Rayleigh scattering is an approximation used to predict scattering from
More informationSolution Set of Homework # 2. Friday, September 09, 2017
Temple University Department of Physics Quantum Mechanics II Physics 57 Fall Semester 17 Z. Meziani Quantum Mechanics Textboo Volume II Solution Set of Homewor # Friday, September 9, 17 Problem # 1 In
More informationInverse Problems in Quantum Optics
Inverse Problems in Quantum Optics John C Schotland Department of Mathematics University of Michigan Ann Arbor, MI schotland@umich.edu Motivation Inverse problems of optical imaging are based on classical
More informationIntro to harmonic analysis on groups Risi Kondor
Risi Kondor Any (sufficiently smooth) function f on the unit circle (equivalently, any 2π periodic f ) can be decomposed into a sum of sinusoidal waves f(x) = k= c n e ikx c n = 1 2π f(x) e ikx dx 2π 0
More informationFFTs in Graphics and Vision. Fast Alignment of Spherical Functions
FFTs in Graphics and Vision Fast Alignment of Spherical Functions Outline Math Review Fast Rotational Alignment Review Recall 1: We can represent any rotation R in terms of the triplet of Euler angles
More informationExact Solutions to the Focusing Discrete Nonlinear Schrödinger Equation
Exact Solutions to the Focusing Discrete Nonlinear Schrödinger Equation Francesco Demontis (based on a joint work with C. van der Mee) Università degli Studi di Cagliari Dipartimento di Matematica e Informatica
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationProblem Solving 1: Line Integrals and Surface Integrals
A. Line Integrals MASSACHUSETTS INSTITUTE OF TECHNOLOY Department of Physics Problem Solving 1: Line Integrals and Surface Integrals The line integral of a scalar function f ( xyz),, along a path C is
More informationRadiative transfer equation in spherically symmetric NLTE model stellar atmospheres
Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres Jiří Kubát Astronomický ústav AV ČR Ondřejov Zářivě (magneto)hydrodynamický seminář Ondřejov 20.03.2008 p. Outline 1.
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 11
WiSe 22..23 Prof. Dr. A-S. Smith Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg Problem. Theoretische Physik 2:
More informationMie Theory, Cross-Sections, Efficiency Factors and Phase Function
Appendix C Mie Theory, Cross-Sections, Efficiency Factors and Phase Function C.1 Mie Theory The concept behind the Mie Theory is explained in section.6.1. However, the mathematical derivation of the theory
More informationGROUP THEORY IN PHYSICS
GROUP THEORY IN PHYSICS Wu-Ki Tung World Scientific Philadelphia Singapore CONTENTS CHAPTER 1 CHAPTER 2 CHAPTER 3 CHAPTER 4 PREFACE INTRODUCTION 1.1 Particle on a One-Dimensional Lattice 1.2 Representations
More informationAngular Momentum - set 1
Angular Momentum - set PH0 - QM II August 6, 07 First of all, let us practise evaluating commutators. Consider these as warm up problems. Problem : Show the following commutation relations ˆx, ˆL x = 0,
More informationIntroduction to Seismology
1.510 Introduction to Seismology Lecture 5 Feb., 005 1 Introduction At previous lectures, we derived the equation of motion (λ + µ) ( u(x, t)) µ ( u(x, t)) = ρ u(x, t) (1) t This equation of motion can
More informationQuantization of the open string on exact plane waves and non-commutative wave fronts
Quantization of the open string on exact plane waves and non-commutative wave fronts F. Ruiz Ruiz (UCM Madrid) Miami 2007, December 13-18 arxiv:0711.2991 [hep-th], with G. Horcajada Motivation On-going
More informationEfficient Solution Methods for Inverse Problems with Application to Tomography Radon Transform and Friends
Efficient Solution Methods for Inverse Problems with Application to Tomography Radon Transform and Friends Alfred K. Louis Institut für Angewandte Mathematik Universität des Saarlandes 66041 Saarbrücken
More informationPhysics 506 Winter 2008 Homework Assignment #4 Solutions. Textbook problems: Ch. 9: 9.6, 9.11, 9.16, 9.17
Physics 56 Winter 28 Homework Assignment #4 Solutions Textbook problems: Ch. 9: 9.6, 9., 9.6, 9.7 9.6 a) Starting from the general expression (9.2) for A and the corresponding expression for Φ, expand
More informationScattering of light from quasi-homogeneous sources by quasi-homogeneous media
Visser et al. Vol. 23, No. 7/July 2006/J. Opt. Soc. Am. A 1631 Scattering of light from quasi-homogeneous sources by quasi-homogeneous media Taco D. Visser* Department of Physics and Astronomy, University
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 3 June, 2004 9 to 12 PAPER 67 PHYSICAL COSMOLOGY Attempt THREE questions. There are four questions in total. The questions carry equal weight. You may not start to
More informationTHE NONCOMMUTATIVE TORUS
THE NONCOMMUTATIVE TORUS The noncommutative torus as a twisted convolution An ordinary two-torus T 2 with coordinate functions given by where x 1, x 2 [0, 1]. U 1 = e 2πix 1, U 2 = e 2πix 2, (1) By Fourier
More informationLongitudinal Waves in Scalar, Three-Vector Gravity
Longitudinal Waves in Scalar, Three-Vector Gravity Kenneth Dalton email: kxdalton@yahoo.com Abstract The linear field equations are solved for the metrical component g 00. The solution is applied to the
More informationQuantization of the E-M field
April 6, 20 Lecture XXVI Quantization of the E-M field 2.0. Electric quadrupole transition If E transitions are forbidden by selection rules, then we consider the next term in the expansion of the spatial
More informationReview of Fundamental Equations Supplementary notes on Section 1.2 and 1.3
Review of Fundamental Equations Supplementary notes on Section. and.3 Introduction of the velocity potential: irrotational motion: ω = u = identity in the vector analysis: ϕ u = ϕ Basic conservation principles:
More informationCoupling of Angular Momenta Isospin Nucleon-Nucleon Interaction
Lecture 5 Coupling of Angular Momenta Isospin Nucleon-Nucleon Interaction WS0/3: Introduction to Nuclear and Particle Physics,, Part I I. Angular Momentum Operator Rotation R(θ): in polar coordinates the
More information5. Suggestions for the Formula Sheets
5. uggestions for the Formula heets Below are some suggestions for many more formulae than can be placed easily on both sides of the two standard 8½" " sheets of paper for the final examination. It is
More informationAnouncements. Assignment 3 has been posted!
Anouncements Assignment 3 has been posted! FFTs in Graphics and Vision Correlation of Spherical Functions Outline Math Review Spherical Correlation Review Dimensionality: Given a complex n-dimensional
More informationAngular Momentum - set 1
Angular Momentum - set PH0 - QM II August 6, 07 First of all, let us practise evaluating commutators. Consider these as warm up problems. Problem : Show the following commutation relations ˆx, ˆL x ] =
More informationAstronomy 9620a / Physics 9302a - 1st Problem List and Assignment
Astronomy 9620a / Physics 9302a - st Problem List and Assignment Your solution to problems 4, 7, 3, and 6 has to be handed in on Thursday, Oct. 7, 203.. Prove the following relations (going from the left-
More informationInverse Transport Problems (Generalized) Stability Estimates. Guillaume Bal
Inverse Transport Problems (Generalized) Stability Estimates Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University Joint with Alexandre Jollivet Data Acquisition in CT-scan
More informationMultipole Expansion for Radiation;Vector Spherical Harmonics
Multipole Expansion for Radiation;Vector Spherical Harmonics Michael Dine Department of Physics University of California, Santa Cruz February 2013 We seek a more systematic treatment of the multipole expansion
More informationCourse Outline. Date Lecture Topic Reading
Course Outline Date Lecture Topic Reading Graduate Mathematical Physics Tue 24 Aug Linear Algebra: Theory 744 756 Vectors, bases and components Linear maps and dual vectors Inner products and adjoint operators
More information1 Coherent-Mode Representation of Optical Fields and Sources
1 Coherent-Mode Representation of Optical Fields and Sources 1.1 Introduction In the 1980s, E. Wolf proposed a new theory of partial coherence formulated in the space-frequency domain. 1,2 The fundamental
More information12 Thermal Radiation including Participating Media
Thermal Radiation including Participating Media Nomenclature A area [m ] c speed of light c specific heat [J/(kgK)] p E emissive power [W/m ] F view factor [-] ij G incident radiation [W/m ] h Planck's
More information3 Angular Momentum and Spin
In this chapter we review the notions surrounding the different forms of angular momenta in quantum mechanics, including the spin angular momentum, which is entirely quantum mechanical in nature. Some
More informationPhysics 216 Problem Set 4 Spring 2010 DUE: MAY 25, 2010
Physics 216 Problem Set 4 Spring 2010 DUE: MAY 25, 2010 1. (a) Consider the Born approximation as the first term of the Born series. Show that: (i) the Born approximation for the forward scattering amplitude
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationWKB Approximation in 3D
1 WKB Approximation in 3D We see solutions ψr of the stationary Schrodinger equations for a spinless particle of energy E: 2 2m 2 ψ + V rψ = Eψ At rst, we just rewrite the Schrodinger equation in the following
More informationPolarized light propagation and scattering in random media
Polarized light propagation and scattering in random media Arnold D. Kim a, Sermsak Jaruwatanadilok b, Akira Ishimaru b, and Yasuo Kuga b a Department of Mathematics, Stanford University, Stanford, CA
More informationUniversity of Illinois at Chicago Department of Physics. Quantum Mechanics Qualifying Examination. January 7, 2013 (Tuesday) 9:00 am - 12:00 noon
University of Illinois at Chicago Department of Physics Quantum Mechanics Qualifying Examination January 7, 13 Tuesday 9: am - 1: noon Full credit can be achieved from completely correct answers to 4 questions
More informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More informationLetting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc.
1 Polynomial Matrices 1.1 Polynomials Letting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc., n ws ( ) as a
More informationPARAXIAL AND EXACT HIGHLY LOCALIZED WAVE SOLUTIONS. Aleksei P. Kiselev. Steklov Math Institute St.Petersburg Department
PARAXIAL AND EXACT HIGHLY LOCALIZED WAVE SOLUTIONS Aleksei P. Kiselev Steklov Math Institute St.Petersburg Department The EXACT wave equation: u xx + u yy + u zz c -2 u tt = 0, c = const, The PARAXIAL
More informationKerr black hole and rotating wormhole
Kerr Fest (Christchurch, August 26-28, 2004) Kerr black hole and rotating wormhole Sung-Won Kim(Ewha Womans Univ.) August 27, 2004 INTRODUCTION STATIC WORMHOLE ROTATING WORMHOLE KERR METRIC SUMMARY AND
More informationTo remain at 0 K heat absorbed by the medium must be removed in the amount of. dq dx = q(l) q(0) dx. Q = [1 2E 3 (τ L )] σ(t T 4 2 ).
294 RADIATIVE HEAT TRANSFER 13.7 Two infinite, isothermal plates at temperatures T 1 and T 2 are separated by a cold, gray medium of optical thickness τ L = L (no scattering). (a) Calculate the radiative
More informationMORE NOTES FOR MATH 823, FALL 2007
MORE NOTES FOR MATH 83, FALL 007 Prop 1.1 Prop 1. Lemma 1.3 1. The Siegel upper half space 1.1. The Siegel upper half space and its Bergman kernel. The Siegel upper half space is the domain { U n+1 z C
More informationPHYS 705: Classical Mechanics. Rigid Body Motion Introduction + Math Review
1 PHYS 705: Classical Mechanics Rigid Body Motion Introduction + Math Review 2 How to describe a rigid body? Rigid Body - a system of point particles fixed in space i r ij j subject to a holonomic constraint:
More informationThe Fractional Fourier Transform with Applications in Optics and Signal Processing
* The Fractional Fourier Transform with Applications in Optics and Signal Processing Haldun M. Ozaktas Bilkent University, Ankara, Turkey Zeev Zalevsky Tel Aviv University, Tel Aviv, Israel M. Alper Kutay
More informationLaser Beam Interactions with Solids In absorbing materials photons deposit energy hc λ. h λ. p =
Laser Beam Interactions with Solids In absorbing materials photons deposit energy E = hv = hc λ where h = Plank's constant = 6.63 x 10-34 J s c = speed of light Also photons also transfer momentum p p
More informationINFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS
INFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS Michael Newman Department of Aerospace Engineering Texas A&M University 34 TAMU 744D H.R. Bright Building College Station,
More informationKIPMU Set 1: CMB Statistics. Wayne Hu
KIPMU Set 1: CMB Statistics Wayne Hu CMB Blackbody COBE FIRAS spectral measurement. yellblackbody spectrum. T = 2.725K giving Ω γ h 2 = 2.471 10 5 12 GHz 200 400 600 10 B ν ( 10 5 ) 8 6 4 error 50 2 0
More informationFourier Approach to Wave Propagation
Phys 531 Lecture 15 13 October 005 Fourier Approach to Wave Propagation Last time, reviewed Fourier transform Write any function of space/time = sum of harmonic functions e i(k r ωt) Actual waves: harmonic
More informationPhoton-atom scattering
Photon-atom scattering Aussois 005 Part Ohad Assaf and Aharon Gero Technion Photon-atom scattering Complex system: Spin of the photon: dephasing Internal atomic degrees of freedom (Zeeman sublevels): decoherence
More informationQFT Perturbation Theory
QFT Perturbation Theory Ling-Fong Li Institute) Slide_04 1 / 44 Interaction Theory As an illustration, take electromagnetic interaction. Lagrangian density is The combination is the covariant derivative.
More informationLecture 11 March 1, 2010
Physics 54, Spring 21 Lecture 11 March 1, 21 Last time we mentioned scattering removes power from beam. Today we treat this more generally, to find the optical theorem: the relationship of the index of
More informationMatrices. Background mathematics review
Matrices Background mathematics review David Miller Matrices Matrix notation Background mathematics review David Miller Matrix notation A matrix is, first of all, a rectangular array of numbers An M N
More informationChapter 2 Acoustical Background
Chapter 2 Acoustical Background Abstract The mathematical background for functions defined on the unit sphere was presented in Chap. 1. Spherical harmonics played an important role in presenting and manipulating
More informationScattering Techniques in Soft Condensed Matter
Scattering Techniques in Soft Condensed Matter Arnout Imhof Soft Condensed Matter, Debye Institute, Utrecht University Scattering of waves Scattering occurs whenever a medium is inhomogeneous Waves scatter
More information