STOCHASTIC & DETERMINISTIC SOLVERS
|
|
- Derick Terry
- 6 years ago
- Views:
Transcription
1 STOCHASTIC & DETERMINISTIC SOLVERS
2 Outline Spatial Scales of Optical Technologies & Mathematical Models Prototype RTE Problems 1-D Transport in Slab Geometry: Exact Solution Stochastic Models: Monte Carlo Simulation Fidelity of Stochastic Model: Model Equivalence P N & Standard Diffusion Approximation (SDA) Fluence Comparisons: RTE vs SDA
3 Spatial Scales of Current Imaging Techniques
4 Need for Model-Based Computation Optical technologies rely upon the precise delivery, detection, and analysis of light signals for (a) accurate dosimetry planning (forward RTE problems) (b) extraction of optical and physiological properties from measurements (inverse RTE problems) (c) guidance in instrument measurement and design (involves both forward and inverse RTE solvers However, model-based computation has lagged behind development of optical technologies Virtual Photonics Core Technology in LAMMP
5 Spatial Scales of Mathematical Models
6 Resolution vs. Depth linked to Physics
7 With Different Optical Regimes
8 Core LAMMP Optical Technologies & Mathematical Models SDA RTE Maxwell s eqs
9 Exact RTE Solutions A small number of RTE equations can be solved in closed form. An important family of such problems describes transport in one dimension These provide interesting and useful benchmarks for approximate RTE methods Surprisingly perhaps, they capture certain features of 3D RTE solutions quite well We introduce this 1D family next
10 RTE in 1-D Slab Geometry A.J. Welch and M.J.C. van Gemert (eds.), Optical-Thermal Response of Laser-Irradiated Tissue, Plenum, 1995 Input: μ a,μ s,d (slab thickness), source Q 0 = 1 p(μ ) = p f δ(μ-1) + p b δ(μ+1) p f + p b = 1, p f p b = g d L x c x c x ( ) 1exp( ) 2exp( ) L x d x d x ( ) 1exp( ) 2exp( ) where c 1, c 2, d 1, d 2 depend on input data and a( a ' s) atr This characteristic decay rate for RTE is different than the characteristic decay rate for SDA, as we will see
11 RTE Solutions We focus on the steady state RTE L( r, ) ( r) L( r, ) ( r) L( r, ') p( r, ' ) d' Q( r, ) s S 2 t Solution L(r, Ω) determines light field everywhere Add detection f(r,ω) to indicate location(s) and angular span(s) for measurement(s) Add target(s): region(s) in phase space needing maximal interrogation These ideas are schematically described on next slide
12 Depictions of RTE Problem Classes Radiance (blue) produced by collimated point source (left) Photon banana (red) indicates paths that reach detector from source (middle) Information density function (red) reveals whether targeted region is receiving adequate light that reaches detector (right)
13 Practical Utility of RTE Problem Types In Lecture 1 it was noted that measurements relevant to therapy & fluorescence excitation rely on virtual detection within the tissue, while diagnostics and imaging are detected by measurements at the boundaries of the tissue. In fact, capturing light signals anywhere can be expressed as weighted integrals of the radiance f ( r,, t) L( r,, t) drddt where the function f ( r,, t), the detector function, restricts the light collected to specific locations r, orientations Ω and times t, either inside the tissue (therapeutics) or on the surface (diagnostics)
14 Forward/Inverse Problems Given problem input: I= {μ t, μ a,p, n, Q, Γ, BC, IC, f₁, f₂,, f D } Determine output O= {I₁, I₂,, I D }, D=# of measurements (detectors) The forward problem map: F: I O The inverse problem map: F -1 :O I In general, F <1 (F is contractive) while O >1 (inverse problem is ill-conditioned). F depends on choice of model and implementation strategy
15 Forward Light Propagation Models Problem Definition layer of VTS establishes the problem input Solver choice (EM, RTE, SDA, P n ) determines mathematical model that produces the RTE solution and (weighted) integrals of describe the problem output Lr (, ) Lr (, ) that VTS platform In this lecture we discuss stochastic (Monte Carlo) & diffusionbased RTE solvers (including SDA)
16 Stochastic Solver Ideas Assumption: RTE describes average behavior of photons; i.e., RTE describes the average behavior of a stochastically varying ensemble of particles RTE coefficients μ t, μ s, p and source Q determine the radiance, L generate photon biographies by A. initiating biography from Q, applying Beer s Law (using μ t, μ s ) to determine intercollision distances, B. using p(μ 0 ) to determine changes in direction; μ 0 = Ω Ω C. alternating sampling between distance & direction until photon leaves tissue or is absorbed
17 Analog Monte Carlo Simulation Each photon represents a unit of radiant energy that is unvarying throughout the biography Energy is conserved: each photon is either absorbed, specularly or diffusely reflected or transmitted Recall Beer s Law: In homogeneous media: L(r, Ω) = L 0 exp(- µ t s) Distances between successive collisions ~ µ t exp(- µ t s) (produces correct total mean free path, 1/ µ t ) Changes in direction ~ p(ω Ω ) = p(μ 0 = cos θ) Photon biographies each contribute tally of 1 either to absorbed energy, reflectance or transmittance
18 Weighted Monte Carlo Simulation Each photon represents a beam of continuously (or discretely) varying radiant intensity, with unit initial intensity Energy is conserved, as before Launch each photon with weight W = 1 Distances between successive scattering collisions ~ µ s exp(- µ s s) (produces correct scattering mean free path, 1/ µ s ) Distribute a weight W [1- exp(- µ a s)] of absorbed energy density along track of length s, and scatter the complementary weight W exp(- µ a s) making use of the phase function p; i.e., changes in direction determined from p(ω Ω ) = p(cos θ) Photon biographies deposit correct weights for absorption, reflectance, transmittance
19 Illustration: Continuous Absorption Weighting This method (CAW) performs weighted MC simulation in tissue in which only scattering is modeled, but absorption is nevertheless accounted for using weights For any µ a > 0, absorbed and scattered energy densities are estimated as on previous slide, using weighted estimators Method applies much more generally, with appropriate modifications A single set of photon biographies can be assigned a set of weights to mimic tissue with arbitrary optical properties, layer thicknesses, etc. This idea is crucial to the use of perturbation Monte Carlo (pmc) in studying sensitivities and solving inverse problems
20 Flow Charts: Weighted (CAW) & Analog MC Weighted MC Analog MC
21 Fidelity of Stochastic Model: Model Equivalence There are (infinitely) many equivalent stochastic solution models for every RTE. E.g, Discrete Absorption Weighting (DAW) de-weights) only at discrete collision events by μ s /μ t. How can we be certain the one we ve chosen models RTE exactly? Each model produces AN N If one can show that N as that is, N 0 A N that, then AN is an unbiased estimator. Intuitively plausible estimators may not possess this property. Conversely, many estimators that are developed strictly from mathematical considerations are not only unbiased but give rise to larger FOMs than more intuitive choices. A
22 Fidelity of Stochastic Model: Model Equivalence Physical/Analytic Model: {Γ, RTE, BCs, f} Stochastic/Probabilistic Model: {B, ξ, M} B = sample space of all possible photon biographies ξ(b) = weight associated with biography b The RTE induces an analog probability measure M on B E[ξ] = (M(B) = 1) B ξ( b) dm ( b) = exact mean value of ξ with respect to M In what sense are these models equivalent? It has been rigorously established J. Spanier and E.M. Gelbard, Monte Carlo Principles and Neutron Transport Problems, Dover, 2008 E[ξ] = B ( b) dm ( b) f ( r, ) L( r, ) drd S 2
23 Spatially Resolved Diffuse Reflectance R( ) L( x, y, z 0, ) d where (ρ 2 = x 2 + y 2 ) and S 2- designates the full hemisphere of exiting directions at the tissue surface, z = 0 (matched n s) The variables/parameters of interest to vary: Distance S 2 between source(s) & detector(s), Number and locations of detectors In practice, light is collected over reduced numerical apertures and small surface areas that include the desired values of ρ; that is, R(, ) L( x, y, z 0, ) dd
24 Time-Resolved Diffuse Reflectance R(,, t) L( x, y, z 0,, t) ddt t Time-resolved measurements depend on introduction of a pulsed light source and collection of light that arrives at each detector at different times (equivalently, distances traversed) For both spatially- and time-resolved measurements, avoiding redundancy requires intelligent placement of detectors
25 Approximating RTE Solutions: Recapturing Some Angular Variation (SDA & P N ) We revisit the steady state RTE L( r, ) ( r) L( r, ) ( r) L( r, ') p( r, ' ) d' Q( r, ) s S 2 t to sketch the ideas that lead to the standard diffusion approximation (SDA) and higher order (P N ) methods. For simplicity we specialize to an infinite plane geometry, reducing L to a function of the spatial variable, x, and the cosine, μ, of arccos (Ω Ω x ). With μ 0 = cos (Ω Ω ), the RTE reduces to dl L ( x, ) t ( x ) L ( x, ) dx 1 s ( x) (, ') (, 0 L x p x ) d ' Q( x, ) 1
26 Derivation of P N Equations The general approach is to expand all functions in Legendre polynomials, P j : 2j 1 L( x, ) Lj( x) Pj( ) 2j 1 Q( x, ) j0 2 Qj( x) Pj( ) j0 2 2j 1 p( 0) pjpj( 0) j0 2 Because the P j are orthogonal 1 and normalized, Lj( x) (, ) j L x P ( x) d and since P 0 =1 and P 1 = μ we find that 1 L0( x) L( x, ) d ( x) is the fluence and N = L1 x L x d J x 1 ( ) (, ) ( ) 1 3 L x L x L x 2 2 (, ) 0( ) 1( ) is the net energy flux
27 The Standard Diffusion Approximation (SDA) Substituting into the RTE we find dl1( x) dl0( x) ( t sp0) L0( x) Q0( x) dx dx 1 3 L( x, ) L0( x) L1( x) K1(1 3 D)exp( x) K 2(1 3 D)3 L1( x)exp( x) 2 where 1 D and SDA a / D 3 a( a s ') 3atr 3( t sg) which compares with RTE a tr 3( t sp1) L1( x) 3 Q1( x) found for RTE
28 The Significance of μ s in SDA Modeling From the SDA equations dl1 x ( ) dx ( t sp0) L0( x) Q0( x) dl0 x ( ) dx 3( t sp1) L1( x) 3 Q1( x) Since p 0 = 1 and p 1 = g, we note that L 0 depends only on μ a + μ s (1-g) = μ a + μ s = μ tr and L 1 only depends on μ a and L 0 For fixed μ a, all (μ s, g) pairs with identical μ s values give rise to identical SDA solutions. this feature can be used advantageously in both SDA and RTE modeling through choices of pairs (μ s, g) that offer computational advantages however, it also means that one can expect difficulty in predicting the optical property distribution (inverse problem solution) from measurements of the RTE solution
29 Model Transport Problem: RTE & SDA S.A. Dupree & S.K. Fraley, A Monte Carlo Primer, Plenum, 2002 It is instructive to examine fluence distributions produced by an isotropic point source in an infinite medium: We compare RTE and SDA solutions in this spherical geometry for different ratios of scattering to absorption and different anisotropies g ( r) ( r) ( r) ( r r ') 2 2 eft 1 ( r) exp( eft r r ' ), eft 3atr 4 r r'
30 Fluences with Varying μ s /μ a, g l* = 1, g = 0, μ s /μ a = 1, 10, 100
31 Fluences with Varying μ s /μ a, g l* = 1, g = 0.8, μ s /μ a = 1, 10, 100
32 Summary and Take Home Messages SDA and RTE are basic and important VTS solvers Monte Carlo simulation is a stochastic model for solving the RTE that offers great generality: FOM measures computational efficiency of Monte Carlo simulations and is easy to estimate There are many stochastically equivalent MC simulations based on weighted Monte Carlo with different FOMs Exact SDA solvers are formula-based and fast (for simple geometries and boundary conditions) Choosing the right solver requires understanding the range and limitations of each
33 The VTS Computational Platform After lunch today, Lisa will present an overview of the VTS computational platform This platform provides access to both command line and GUI-enabled learning & research capabilities We ll use VTS throughout the Short Course and in the hands-on computational labs each afternoon
Monte Carlo Methods:
Short Course on Computational Monte Carlo Methods: Fundamentals Shuang Zhao Assistant Professor Computer Science Department University of California, Irvine Shuang Zhao 1 Teaching Objective Introducing
More informationFrequency domain photon migration in the -P 1 approximation: Analysis of ballistic, transport, and diffuse regimes
Frequency domain photon migration in the -P 1 approximation: Analysis of ballistic, transport, and diffuse regimes J. S. You, 1,2 C. K. Hayakawa, 2 and V. Venugopalan 1,2,3, * 1 Department of Biomedical
More informationMCRT L10: Scattering and clarification of astronomy/medical terminology
MCRT L10: Scattering and clarification of astronomy/medical terminology What does the scattering? Shape of scattering Sampling from scattering phase functions Co-ordinate frames Refractive index changes
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 informationThis is the third of three lectures on cavity theory.
This is the third of three lectures on cavity theory. 1 In this lecture, we are going to go over what is meant by charged particle equilibrium and look at the dose and kerma when you have charged particle
More informationChapter V: Cavity theories
Chapter V: Cavity theories 1 Introduction Goal of radiation dosimetry: measure of the dose absorbed inside a medium (often assimilated to water in calculations) A detector (dosimeter) never measures directly
More informationLight Electromagnetic Radiation
www.atomic.physics.lu.se/biophotonics Biomedicinsk Optik Ljusutbredning i vävnad STEFAN ANDERSSON-ENGELS 1 Light Electromagnetic Radiation Soto Thompson, PhD Thesis, 2004 2 1 Ionizing versus optical radiation
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 informationMultiple-source optical diffusion approximation for a multilayer scattering medium
Multiple-source optical diffusion approximation for a multilayer scattering medium Joseph. Hollmann 1 and ihong V. Wang 1,2, * 1 Optical Imaging aboratory, Department of Biomedical Engineering, Texas A&M
More informationLecture 06. Fundamentals of Lidar Remote Sensing (4) Physical Processes in Lidar
Lecture 06. Fundamentals of Lidar Remote Sensing (4) Physical Processes in Lidar Physical processes in lidar (continued) Doppler effect (Doppler shift and broadening) Boltzmann distribution Reflection
More informationMCRT: L4 A Monte Carlo Scattering Code
MCRT: L4 A Monte Carlo Scattering Code Plane parallel scattering slab Optical depths & physical distances Emergent flux & intensity Internal intensity moments Constant density slab, vertical optical depth
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 informationMean Intensity. Same units as I ν : J/m 2 /s/hz/sr (ergs/cm 2 /s/hz/sr) Function of position (and time), but not direction
MCRT: L5 MCRT estimators for mean intensity, absorbed radiation, radiation pressure, etc Path length sampling: mean intensity, fluence rate, number of absorptions Random number generators J Mean Intensity
More informationLichtausbreitung in streuenden Medien: Prinzip und Anwendungsbeispiele
Lichtausbreitung in streuenden Medien: Prinzip und Anwendungsbeispiele Alwin Kienle 06.12.2013 Institut für Lasertechnologien in der Medizin und Meßtechnik an der Universität Ulm Overview 1) Theory of
More informationPC4262 Remote Sensing Scattering and Absorption
PC46 Remote Sensing Scattering and Absorption Dr. S. C. Liew, Jan 003 crslsc@nus.edu.sg Scattering by a single particle I(θ, φ) dφ dω F γ A parallel beam of light with a flux density F along the incident
More informationDiscrete-ordinates solution of short-pulsed laser transport in two-dimensional turbid media
Discrete-ordinates solution of short-pulsed laser transport in two-dimensional turbid media Zhixiong Guo and Sunil Kumar The discrete-ordinates method is formulated to solve transient radiative transfer
More informationMonte Carlo Radiation Transfer I
Monte Carlo Radiation Transfer I Monte Carlo Photons and interactions Sampling from probability distributions Optical depths, isotropic emission, scattering Monte Carlo Basics Emit energy packet, hereafter
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 informationMonte Carlo Sampling
Monte Carlo Sampling Sampling from PDFs: given F(x) in analytic or tabulated form, generate a random number ξ in the range (0,1) and solve the equation to get the randomly sampled value X X ξ = F(x)dx
More informationMonte Carlo method projective estimators for angular and temporal characteristics evaluation of polarized radiation
Monte Carlo method projective estimators for angular and temporal characteristics evaluation of polarized radiation Natalya Tracheva Sergey Ukhinov Institute of Computational Mathematics and Mathematical
More informationAnalysis of Scattering of Radiation in a Plane-Parallel Atmosphere. Stephanie M. Carney ES 299r May 23, 2007
Analysis of Scattering of Radiation in a Plane-Parallel Atmosphere Stephanie M. Carney ES 299r May 23, 27 TABLE OF CONTENTS. INTRODUCTION... 2. DEFINITION OF PHYSICAL QUANTITIES... 3. DERIVATION OF EQUATION
More informationDigital Holographic Measurement of Nanometric Optical Excitation on Soft Matter by Optical Pressure and Photothermal Interactions
Ph.D. Dissertation Defense September 5, 2012 Digital Holographic Measurement of Nanometric Optical Excitation on Soft Matter by Optical Pressure and Photothermal Interactions David C. Clark Digital Holography
More informationLecture 4* Inherent optical properties, IOP Theory. Loss due to absorption. IOP Theory 12/2/2008
Lecture 4* Inherent optical properties, part I IOP Theory What properties of a medium affect the radiance field as it propagate through it? 1. Sinks of photons (absorbers) 2. Sources of photons (internal
More informationTransport under magnetic fields with the EGSnrc simulation toolkit
Transport under magnetic fields with the EGSnrc simulation toolkit Ernesto Mainegra-Hing, Frédéric Tessier, Blake Walters Measurement Science and Standards, National Research Council Canada Hugo Bouchard
More informationMonte Carlo Methods & Virtual Photonics
Monte Carlo Methods & Virtual Photonics Tuesday-Thursday 9 10:30 Room 3201, Nat. Sci. II Jerry Spanier P220 BLI jspanier@uci.edu Please help yourselves to refreshments at the back of the room and leave
More informationSupplementary Figure S1 Definition of the wave vector components: Parallel and perpendicular wave vector of the exciton and of the emitted photons.
Supplementary Figure S1 Definition of the wave vector components: Parallel and perpendicular wave vector of the exciton and of the emitted photons. Supplementary Figure S2 The calculated temperature dependence
More informationScattering in Cold- Cathode Discharges
Simulating Electron Scattering in Cold- Cathode Discharges Alexander Khrabrov, Igor Kaganovich*, Vladimir I. Demidov**, George Petrov*** *Princeton Plasma Physics Laboratory ** Wright-Patterson Air Force
More informationULTRAFAST LASER PULSE TRAIN RADIATION TRANSFER IN A SCATTERING-ABSORBING 3D MEDIUM WITH AN INHOMOGENEITY
Heat Transfer Research 46(9), 861 879 (2015) ULTRAFAST LASER PULSE TRAIN RADIATION TRANSFER IN A SCATTERING-ABSORBING 3D MEDIUM WITH AN INHOMOGENEITY Masato Akamatsu 1,* & Zhixiong Guo 2 1 Graduate School
More information11. Point Detector Tally
and suppose that we want half the source to be in the range 005 < E < 1 and the other half to be in the range 1 < E < 20 Then the input is SPn -5 a SIn 005 1 20 SBn C 0 5 1 MCNP breaks up the function
More informationReview Article Mathematical Methods in Biomedical Optics
ISRN Biomedical Engineering Volume 2013, Article ID 464293, 8 pages http://dx.doi.org/10.1155/2013/464293 Review Article Mathematical Methods in Biomedical Optics Macaveiu Gabriela Academic Center of Optical
More informationFluorescence Workshop UMN Physics June 8-10, 2006 Quantum Yield and Polarization (1) Joachim Mueller
Fluorescence Workshop UMN Physics June 8-10, 2006 Quantum Yield and Polarization (1) Joachim Mueller Quantum yield, polarized light, dipole moment, photoselection, dipole radiation, polarization and anisotropy
More information8: Source-Sink Problems in 1 Energy Group
8: Source-Sink Problems in 1 Energy Group B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 015 Sept.-Dec. 015 September 1 Contents Solving the 1-group diffusion
More informationFundametals of Rendering - Radiometry / Photometry
Fundametals of Rendering - Radiometry / Photometry Physically Based Rendering by Pharr & Humphreys Chapter 5: Color and Radiometry Chapter 6: Camera Models - we won t cover this in class Realistic Rendering
More informationRadiation Quantities and Units
Radiation Quantities and Units George Starkschall, Ph.D. Lecture Objectives Define and identify units for the following: Exposure Kerma Absorbed dose Dose equivalent Relative biological effectiveness Activity
More informationLesson 6: Diffusion Theory (cf. Transport), Applications
Lesson 6: Diffusion Theory (cf. Transport), Applications Transport Equation Diffusion Theory as Special Case Multi-zone Problems (Passive Media) Self-shielding Effects Diffusion Kernels Typical Values
More informationTERRESTRIAL REDSHIFTS FROM A DIFFUSE LIGHT SOURCE. FOM Institute for Atomic and Molecular Physics, Kruislaan 407, 1098 SJ Amsterdam, The Netherlands
March 25, 1999 TERRESTRIAL REDSHIFTS FROM A DIFFUSE LIGHT SOURCE Ad Lagendijk a FOM Institute for Atomic and Molecular Physics, Kruislaan 407, 1098 SJ Amsterdam, The Netherlands (Received ) Abstract Enhanced
More informationMODERN OPTICS. P47 Optics: Unit 9
MODERN OPTICS P47 Optics: Unit 9 Course Outline Unit 1: Electromagnetic Waves Unit 2: Interaction with Matter Unit 3: Geometric Optics Unit 4: Superposition of Waves Unit 5: Polarization Unit 6: Interference
More informationFast simulation of finite-beam optical coherence tomography of inhomogeneous turbid media
Fast simulation of finite-beam optical coherence tomography of inhomogeneous turbid media by Micheal Sobhy A thesis submitted to the Faculty of Graduate Studies of the University of Manitoba, in partial
More informationClassical and Improved Diffusion Theory for Subsurface Scattering
Classical and Improved Diffusion Theory for Subsurface Scattering Ralf Habel 1 Per H. Christensen 2 Wojciech Jarosz 1 1 Disney Research Zürich 2 Pixar Animation Studios Technical report, June 2013 Abstract
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationEnrico Fermi and the FERMIAC. Mechanical device that plots 2D random walks of slow and fast neutrons in fissile material
Monte Carlo History Statistical sampling Buffon s needles and estimates of π 1940s: neutron transport in fissile material Origin of name People: Ulam, von Neuman, Metropolis, Teller Other areas of use:
More informationEquivalent isotropic scattering formulation for transient
Equivalent isotropic scattering formulation for transient short-pulse radiative transfer in anisotropic scattering planar media Zhixiong Guo and Sunil Kumar An isotropic scaling formulation is evaluated
More informationLecture 20 Reactor Theory-V
Objectives In this lecture you will learn the following We will discuss the criticality condition and then introduce the concept of k eff.. We then will introduce the four factor formula and two group
More informationPhysics of Radiotherapy. Lecture II: Interaction of Ionizing Radiation With Matter
Physics of Radiotherapy Lecture II: Interaction of Ionizing Radiation With Matter Charge Particle Interaction Energetic charged particles interact with matter by electrical forces and lose kinetic energy
More informationgives rise to multitude of four-wave-mixing phenomena which are of great
Module 4 : Third order nonlinear optical processes Lecture 26 : Third-order nonlinearity measurement techniques: Z-Scan Objectives In this lecture you will learn the following Theory of Z-scan technique
More informationTheoretical Analysis of Thermal Damage in Biological Tissues Caused by Laser Irradiation
Copyright c 007 Tech Science Press MCB, vol.4, no., pp.7-39, 007 Theoretical Analysis of Thermal Damage in Biological Tissues Caused by Laser Irradiation Jianhua Zhou,J.K.Chen and Yuwen Zhang Abstract:
More informationThe Scattering of Light by Small Particles. Advanced Laboratory, Physics 407 University of Wisconsin Madison, Wisconsin 53706
(4/28/09) The Scattering of Light by Small Particles Advanced Laboratory, Physics 407 University of Wisconsin Madison, Wisconsin 53706 Abstract In this experiment we study the scattering of light from
More informationWhat is it good for? RT is a key part of remote sensing and climate modeling.
Read Bohren and Clothiaux Ch.; Ch 4.-4. Thomas and Stamnes, Ch..-.6; 4.3.-4.3. Radiative Transfer Applications What is it good for? RT is a key part of remote sensing and climate modeling. Remote sensing:
More informationSound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur
Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-17 HARMONIC PLANE WAVES Introduction In this lecture, we discuss propagation of 1-D planar waves as they travel
More informationDetermining how uncertainties in optical properties affect light dose calculations
Determining how uncertainties in optical properties affect light dose calculations Julia Sandell 1, Jarod Finlay, Timothy Zhu 1 Department of Physics, University of Pennsylvania Radiation Oncology, Hospital
More informationModelling the Directionality of Light Scattered in Translucent Materials
Modelling the Directionality of Light Scattered in Translucent Materials Jeppe Revall Frisvad Joint work with Toshiya Hachisuka, Aarhus University Thomas Kjeldsen, The Alexandra Institute March 2014 Materials
More informationImproved solutions of the steady-state and the time-resolved diffusion equations for reflectance from a semi-infinite turbid medium
246 J. Opt. Soc. Am. A/Vol. 14, No. 1/January 1997 A. Kienle and M. S. Patterson Improved solutions of the steady-state and the time-resolved diffusion equations for reflectance from a semi-infinite turbid
More informationApparent optical properties and radiative transfer theory*
Apparent optical properties and radiative transfer theory* Apparent optical properties. The RTE and Gershun s equation The Secchi disk (and depth, an AOP). *based in part on lectures by Roesler, Mobley,
More information1 CHAPTER 5 ABSORPTION, SCATTERING, EXTINCTION AND THE EQUATION OF TRANSFER
1 CHAPTER 5 ABSORPTION, SCATTERING, EXTINCTION AND THE EQUATION OF TRANSFER 5.1 Introduction As radiation struggles to make its way upwards through a stellar atmosphere, it may be weakened by absorption
More informationA high accuracy approximation for half - space problems with anisotropic scattering
IOSR Journal of Applied Physics (IOSR-JAP e-iss: 78-486. Volume 4, Issue (May. - Jun., PP 7-46 A high accuracy approximation for half - space problems with anisotropic scattering M. H. Haggag,, H.F. Mashali
More informationBeer-Lambert (cont.)
The Beer-Lambert Law: Optical Depth Consider the following process: F(x) Absorbed flux df abs F(x + dx) Scattered flux df scat x x + dx The absorption or scattering of radiation by an optically active
More informationNeutron Diffusion Theory: One Velocity Model
22.05 Reactor Physics - Part Ten 1. Background: Neutron Diffusion Theory: One Velocity Model We now have sufficient tools to begin a study of the second method for the determination of neutron flux as
More informationMCSHAPE: A MONTE CARLO CODE FOR SIMULATION OF POLARIZED PHOTON TRANSPORT
Copyright JCPDS - International Centre for Diffraction Data 2003, Advances in X-ray Analysis, Volume 46. 363 MCSHAPE: A MONTE CARLO CODE FOR SIMULATION OF POLARIZED PHOTON TRANSPORT J.E. Fernández, V.
More informationDetectors in Nuclear Physics: Monte Carlo Methods. Dr. Andrea Mairani. Lectures I-II
Detectors in Nuclear Physics: Monte Carlo Methods Dr. Andrea Mairani Lectures I-II INTRODUCTION Sampling from a probability distribution Sampling from a probability distribution X λ Sampling from a probability
More informationVacuum Science and Technology in Accelerators
Vacuum Science and Technology in Accelerators Ron Reid Consultant ASTeC Vacuum Science Group (ron.reid@stfc.ac.uk) Cockcroft Institute Lectures 2010 Session 9 Calculations to Support the Design Aims To
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 1. Title: Neutron Life Cycle
Lectures on Nuclear Power Safety Lecture No 1 Title: Neutron Life Cycle Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture Infinite Multiplication Factor, k Four Factor Formula
More informationRadiative transport in the delta-p 1 approximation: accuracy of fluence rate and optical penetration depth predictions in turbid semi-infinite media
Journal of Biomedical Optics 9(3), 632 647 (May/June 24) Radiative transport in the delta-p 1 approximation: accuracy of fluence rate and optical penetration depth predictions in turbid semi-infinite media
More informationATM 507 Lecture 4. Text reading Chapters 3 and 4 Today s topics Chemistry, Radiation and Photochemistry review. Problem Set 1: due Sept.
ATM 507 Lecture 4 Text reading Chapters 3 and 4 Today s topics Chemistry, Radiation and Photochemistry review Problem Set 1: due Sept. 11 Temperature Dependence of Rate Constants Reaction rates change
More informationRadiation processes and mechanisms in astrophysics I. R Subrahmanyan Notes on ATA lectures at UWA, Perth 18 May 2009
Radiation processes and mechanisms in astrophysics I R Subrahmanyan Notes on ATA lectures at UWA, Perth 18 May 009 Light of the night sky We learn of the universe around us from EM radiation, neutrinos,
More informationRadiation in the atmosphere
Radiation in the atmosphere Flux and intensity Blackbody radiation in a nutshell Solar constant Interaction of radiation with matter Absorption of solar radiation Scattering Radiative transfer Irradiance
More informationPhD Qualifying Exam Nuclear Engineering Program. Part 1 Core Courses
PhD Qualifying Exam Nuclear Engineering Program Part 1 Core Courses 9:00 am 12:00 noon, November 19, 2016 (1) Nuclear Reactor Analysis During the startup of a one-region, homogeneous slab reactor of size
More informationEEE4101F / EEE4103F Radiation Interactions & Detection
EEE4101F / EEE4103F Radiation Interactions & Detection 1. Interaction of Radiation with Matter Dr. Steve Peterson 5.14 RW James Department of Physics University of Cape Town steve.peterson@uct.ac.za March
More informationModeling Focused Beam Propagation in scattering media. Janaka Ranasinghesagara, Ph.D.
Modeling Focused Beam Propagation in scattering media Janaka Ranasinghesagara, Ph.D. Teaching Objectives The need for computational models of focused beam propagation in scattering media Introduction to
More informationBLUE PRINT FOR MODEL QUESTION PAPER 4
Unit Chapter Teaching Hours Marks allotted (VSA) mark (SA) 3 mark (SA) 5 mark (LA) 5 mark (NP) BLUE PRINT FOR MODEL QUESTION PAPER 4 SUBJECT : PHYSICS (33) CLASS : II PUC Topic Electric Charges and Fields
More informationAT622 Section 2 Elementary Concepts of Radiometry
AT6 Section Elementary Concepts of Radiometry The object of this section is to introduce the student to two radiometric concepts intensity (radiance) and flux (irradiance). These concepts are largely geometrical
More informationOutline. Absorbed Dose in Radioactive Media. Introduction. Radiation equilibrium. Charged-particle equilibrium
Absorbed Dose in Radioactive Media Chapter F.A. Attix, Introduction to Radiological Physics and Radiation Dosimetry Outline General dose calculation considerations, absorbed fraction Radioactive disintegration
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 1 (2/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications
.54 Neutron Interactions and Applications (Spring 004) Chapter 1 (/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications There are many references in the vast literature on nuclear
More informationForward and inverse wave problems:quantum billiards and brain imaging p.1
Forward and inverse wave problems: Quantum billiards and brain imaging Alex Barnett barnett@cims.nyu.edu. Courant Institute Forward and inverse wave problems:quantum billiards and brain imaging p.1 The
More informationWave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces
Lecture 5: Crystal Optics Outline 1 Homogeneous, Anisotropic Media 2 Crystals 3 Plane Waves in Anisotropic Media 4 Wave Propagation in Uniaxial Media 5 Reflection and Transmission at Interfaces Christoph
More informationChapter VIII: Photographic films
Chapter VIII: Photographic films 1 Photographic films First (integrating) dosimeters to be developed based on photographic techniques dosifilms Dosifilms were more and more replaced by dosimeters using
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 informationStatistical regimes of random laser fluctuations
Statistical regimes of random laser fluctuations Stefano Lepri Istituto dei Sistemi Complessi ISC-CNR Firenze S. Cavalieri, G.-L. Oppo, D.S. Wiersma Stefano Lepri (ISC-CNR) Random laser fluctuations 1
More informationNonclassical Particle Transport in Heterogeneous Materials
Nonclassical Particle Transport in Heterogeneous Materials Thomas Camminady, Martin Frank and Edward W. Larsen Center for Computational Engineering Science, RWTH Aachen University, Schinkelstrasse 2, 5262
More informationLecture 2: principles of electromagnetic radiation
Remote sensing for agricultural applications: principles and methods Lecture 2: principles of electromagnetic radiation Instructed by Prof. Tao Cheng Nanjing Agricultural University March Crop 11, Circles
More informationDiffuse Optical Spectroscopy: Analysis of BiophotonicsSignals Short Course in Computational Biophotonics Albert Cerussi Day 4/Lecture #2
Diffuse Optical Spectroscopy: Analysis of BiophotonicsSignals 2013 Short Course in Computational Biophotonics Albert Cerussi Day 4/Lecture #2 Disclosures No financial conflicts to disclose Learning Objectives
More informationUltrafast X-Ray-Matter Interaction and Damage of Inorganic Solids October 10, 2008
Ultrafast X-Ray-Matter Interaction and Damage of Inorganic Solids October 10, 2008 Richard London rlondon@llnl.gov Workshop on Interaction of Free Electron Laser Radiation with Matter Hamburg This work
More informationGoal: The theory behind the electromagnetic radiation in remote sensing. 2.1 Maxwell Equations and Electromagnetic Waves
Chapter 2 Electromagnetic Radiation Goal: The theory behind the electromagnetic radiation in remote sensing. 2.1 Maxwell Equations and Electromagnetic Waves Electromagnetic waves do not need a medium to
More informationTRANSFER OF RADIATION
TRANSFER OF RADIATION Under LTE Local Thermodynamic Equilibrium) condition radiation has a Planck black body) distribution. Radiation energy density is given as U r,ν = 8πh c 3 ν 3, LTE), tr.1) e hν/kt
More informationNIH Public Access Author Manuscript J Comput Phys. Author manuscript; available in PMC 2009 September 22.
NIH Public Access Author Manuscript Published in final edited form as: J Comput Phys. 2008 December 1; 227(23): 9762 9777. doi:10.1016/j.jcp.2008.07.016. A new proof of geometric convergence for general
More informationEmphasis on what happens to emitted particle (if no nuclear reaction and MEDIUM (i.e., atomic effects)
LECTURE 5: INTERACTION OF RADIATION WITH MATTER All radiation is detected through its interaction with matter! INTRODUCTION: What happens when radiation passes through matter? Emphasis on what happens
More informationBrief Introduction to: Transport Theory/Monte Carlo Techniques/MCNP
22:54 Neutron Interactions and Applications Brief Introduction to: Transport Theory/Monte Carlo Techniques/MCNP The behavior of individual neutrons and nuclei cannot be predicted. However, the average
More informationShielded Scintillator for Neutron Characterization
Shielded Scintillator for Neutron Characterization A Thesis Submitted in Partial Fulfillment of the Requirements for Graduation with Research Distinction in Engineering Physics By Patrick X. Belancourt
More informationEnergy Dependence of Neutron Flux
Energy Dependence of Neutron Flux B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec. 2015 September 1 Contents We start the discussion of the energy
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 information2. NOTES ON RADIATIVE TRANSFER The specific intensity I ν
1 2. NOTES ON RADIATIVE TRANSFER 2.1. The specific intensity I ν Let f(x, p) be the photon distribution function in phase space, summed over the two polarization states. Then fdxdp is the number of photons
More informationToday, I will present the first of two lectures on neutron interactions.
Today, I will present the first of two lectures on neutron interactions. I first need to acknowledge that these two lectures were based on lectures presented previously in Med Phys I by Dr Howell. 1 Before
More informationAnswers to questions on exam in laser-based combustion diagnostics on March 10, 2006
Answers to questions on exam in laser-based combustion diagnostics on March 10, 2006 1. Examples of advantages and disadvantages with laser-based combustion diagnostic techniques: + Nonintrusive + High
More informationSecond-Order Linear ODEs (Textbook, Chap 2)
Second-Order Linear ODEs (Textbook, Chap ) Motivation Recall from notes, pp. 58-59, the second example of a DE that we introduced there. d φ 1 1 φ = φ 0 dx λ λ Q w ' (a1) This equation represents conservation
More informationHigh-Order Similarity Relations in Radiative Transfer: Supplementary Document
High-Order Similarity Relations in Radiative Transfer: Supplementary Document Shuang hao Ravi Ramamoorthi Kavita Bala Cornell University University of California, Berkeley Cornell University Overview This
More informationLET! (de / dx) 1 Gy= 1 J/kG 1Gy=100 rad. m(kg) dose rate
Basics of Radiation Dosimetry for the Physicist http://en.wikipedia.org/wiki/ionizing_radiation I. Ionizing radiation consists of subatomic particles or electromagnetic waves that ionize electrons along
More informationCRaTER Pre-Ship Review (PSR) Instrument Calibration Science Requirements Compliance
CRaTER Pre-Ship Review (PSR) Instrument Calibration Science Requirements Compliance Justin C Kasper Smithsonian Astrophysical Observatory January 3, 2008 Outline Calibration Relate the ADU of the Pulse
More informationThermal Radiation Studies for an Electron-Positron Annihilation Propulsion System
Thermal Radiation Studies for an Electron-Positron Annihilation Propulsion System Jonathan A. Webb Embry Riddle Aeronautical University Prescott, AZ 8631 Recent studies have shown the potential of antimatter
More informationAbsorption length for photon propagation in highly dense colloidal dispersions
Journal of MATERIALS RESEARCH Welcome Comments Help Absorption length for photon propagation in highly dense colloidal dispersions Rajeev Garg, Robert K. Prud homme, and Ilhan A. Aksay Department of Chemical
More informationLecture 9: Introduction to Diffraction of Light
Lecture 9: Introduction to Diffraction of Light Lecture aims to explain: 1. Diffraction of waves in everyday life and applications 2. Interference of two one dimensional electromagnetic waves 3. Typical
More informationFundamental Concepts of Radiometry p. 1 Electromagnetic Radiation p. 1 Terminology Conventions p. 3 Wavelength Notations and Solid Angle p.
Preface p. xiii Fundamental Concepts of Radiometry p. 1 Electromagnetic Radiation p. 1 Terminology Conventions p. 3 Wavelength Notations and Solid Angle p. 4 Fundamental Definitions p. 7 Lambertian Radiators
More information