基于Kubelka-Munk理论的半色调印刷品光谱反射预测模型

首先介绍了经典的Kubelka-Munk理论及其Saunderson校正,并在此基础上推广到半色调印刷品中,考虑到光在印刷品中的横向散射,引入了散射概率函数描述光子从一个界面元进入而从另一个界面元射出的概率,建立了一个统一的半色调印刷品光谱反射预测模型.     

Photon migration in Raman spectroscopy

Monte Carlo simulation has been applied to study time-resolved Raman and Tyndall photon migration in opaque samples under isotropic and forward scattering conditions. For isotropic scattering, Raman and Tyndall intensities are predicted to decay according to t((1-n)) and t(-n), respectively, where the value of n depends on the ratio of the optical collection aperture to the mean scattering length. The simulation correctly reproduced the analytical results of n = 3/2 and n = 5/2 for a point source in infinite and semi-infinite media, respectively. In addition the model can be used to relate the time at which a Raman photon exits the sample to the mean depth at which it was generated. This could provide a useful tool for depth profiling the chemical composition of turbid systems, and hence be a useful addition to the established array of photon-migration techniques. The model was applied to analyze experimentally observed Raman and Tyndall decay profiles from powdered trans-stilbene. The transport mean free path (l(t)) was calculated to be approximately 400 microm, which was significantly larger than the particle sizes present in the sample (approximately 10-100 microm). This implies that the particles were highly forward scattering, as would be expected for this size range. When highly anisotropic scattering was introduced into the model a much more reasonable scattering length (l(s) approximately 40 microm) was obtained. Finally, a simple analytical model was developed that gives the correct relationship between the Raman and Tyndall decay exponents for isotropic scattering. To the best of our knowledge this work represents the first detailed study of Raman photon migration under time-resolved conditions.     

Power law analysis estimates of analyte concentration and particle size in highly scattering granular samples from photon time-of-flight measurements

Optical measurements of particle size and composition in granular samples are difficult to make due to complex light scattering from particles. These multiple scattering events bias absorption estimates and complicate the calculation of scattering and absorption coefficients used to estimate sample properties. Time series data, such as chromatograms and photon time-of-flight (TOF) profiles, contain self-repeating (fractal) characteristics. Power law analysis of photon TOF profiles allows the determination of absorption coefficients and particle sizes in a single experiment. A correlation dimension algorithm was used on photon TOF data from scattering samples. MLR models were then obtained from correlation dimension plots for the estimation of sample properties. Estimates of particle sizes and absorption coefficients were shown to agree well with theoretical values when compared using independent validation sets. Results show close to a 3-fold and up to a 5-fold decrease in the errors of estimation of dye concentration and particle size, respectively, as compared to steady-state measurements. The power law approach provides a useful means of determining sample properties in highly scattering media.     

Effective Mie scattering of a spherical fractal aggregate and its application in turbid media

An effective Mie-scattering model is developed to deal with the scattering property of a spherical fractal aggregate consisting of scattering particles. In this model the scattered field of a scattering particle is given by the classical Mie-scattering theory. On the basis of the Monte Carlo simulation method, we determine the physical parameters of a scattering aggregate, the scattering efficiency Q, and the anisotropy value g, as well as their dependence on the size and the effective mean-free-path length of a scattering aggregate. Accordingly, photon migration through a microscope objective focused into a turbid medium including scattering aggregates is simulated to understand the effect of complex tissue on image quality.     

Application of frequency domain photon migration to particle size analysis and monitoring of pharmaceutical powders

The frequency domain photon migration (FDPM) technique was employed to determine mean particle size of pharmaceutical powders. Results show that the FDPM-measured scattering coefficient increases linearly with reciprocal mean particle size of powdered samples. In contrast to near-infrared spectroscopy techniques, FDPM technique enables determination of scattering and absorption separately so that it does not require data pretreatment and chemometric calibration models. In addition, this unique advantage provides more detailed information about powder samples, which can be used as a potential tool for on-line monitoring of not only variation of active pharmaceutical ingredient concentrations from changes in the absorption coefficient but also variation of particle sizes from changes in the scattering coefficient.     

Determination of the droplet effective size and optical depth of cloudy media from polarimetric measurements: theory

We present the development of a semi-analytical algorithm for optical particle sizing in disperse media. The algorithm is applied to the specific case of water clouds. However, it can be extended with minor modifications to other types of light-scattering medium. It is assumed that the optical thickness tau of the medium is large and the probability of photon absorption beta is small. Thus the optical particle-sizing problem is studied in the regime of highly developed multiple light scattering. It was found that the degree of polarization in visible and near-infrared channels provides us with information both on the effective size of droplets and on the optical thickness tau.     

Monte Carlo study of coherent diffuse photon transport in a homogeneous turbid medium: a degree-of-coherence based approach

We employ a Monte Carlo (MC) algorithm to investigate the decoherence of diffuse photons in turbid media. For the MC simulation of coherent photons, the degree of coherence, defined as a random variable for a photon packet, is associated with a decoherence function that depends on the scattering angle and is updated as a photon interacts with a medium via scattering. Using a slab model, the effects of medium scattering properties were studied, which reveals that a linear random variable model for the degree of coherence is in better agreement with experimental results than a sinusoidal model and that decoherence is quick for the initial few scattering events followed by a slow and gradual decrease of coherence.     

Light scattering of ultrafine silica particles by VUV synchrotron radiation

Vacuum ultraviolet (VUV) light scattering from ultrafine silica particles is studied with an aerosol instrument recently established at the Advanced Light Source (ALS) in Berkeley. Silica particles, size-selected by a differential mobility analyzer, are introduced into vacuum through a set of aerodynamic lenses to form a particle beam. The scattered photons from the crossing area of the VUV synchrotron beam and particle beam are detected with a rotatable VUV photon detector. The angular distributions of scattered photons (ADSP) originating from 70, 100, 200 nm diameter silica particles are measured with 145.9 and 118.1 nm synchrotron radiation. These angular distributions show strong forward scattering. The measured ADSPs are consistent with simulation of Mie scattering. The refractive indices of silica particles, 2.6 + 1.1i and 1.6 + 0.0001i for 118.1 and 145.9 nm, respectively, are obtained by fitting the measured ADSPs; the least average percentage deviations are 18% and 6%, respectively. The scattered fluxes at widely different wavelengths (visible versus VUV) also exhibit clear size sensitivity. Under comparable experimental conditions of photon fluxes and detection efficiencies, limits of particle size detection of 70 and 250 nm are obtained, respectively, when using 118.1 and 532 nm illumination. As anticipated, VUV scattering is a more sensitive probe for ultrafine particles, which will find application in detection of these ubiquitous species beyond the confines of a laboratory.     

QUADOS intercomparison: a summary of photon and charged particle problems

QUADOS, a Concerted Action of the European Commission, has promoted an intercomparison aimed at evaluating the use of computational codes for dosimetry in radiation protection and medical physics. This intercomparison was open to all users of radiation transport codes. Eight problems were selected for their relevance to the radiation dosimetry community, five of which involved photon and proton transport. This paper focuses on the analysis of the photon and charged particle problems. The neutron problems were presented in a paper at the NEUDOS9 conference.     

Excitation with a focused, pulsed optical beam in scattering media: diffraction effects

To gain a better understanding of the spatiotemporal problems that are encountered in two-photon excitation fluorescence imaging through highly scattering media, we investigate how diffraction affects the three-dimensional intensity distribution of a focused, pulsed optical beam propagating inside a scattering medium. In practice, the full potential of the two-photon excitation fluorescence imaging is unrealized at long scattering depths, owing to the unwanted temporal and spatial broadening of the femtosecond excitation light pulse that reduces the energy density at the geometric focus while it increases the excitation energy density in the out-of-focus regions. To analyze the excitation intensity distribution, we modify the Monte Carlo-based photon-transport model to a semi-quantum-mechanical representation that combines the wave properties of light with the particle behavior of the propagating photons. In our model the propagating photon is represented by a plane wave with its propagation direction in the scattering medium determined by the Monte Carlo technique. The intensity distribution in the focal region is given by the square of the linear superposition of the various plane waves that arrive at different incident angles and optical path lengths. In the absence of scattering, the propagation model yields the intensity distribution that is predicted by the Huygens-Fresnel principle. We quantify the decrease of the energy density delivered at the geometric focus as a function of the optical depth to the mean-free-path ratio that yields the average number of scattering events that a photon encounters as it propagates toward the focus. Both isotropic and anisotropic scattering media are considered. Three values for the numerical aperture (NA) of the focusing lens are considered: NA = 0.25, 0.5, 0.75.     

Monte Carlo simulations of the diffuse backscattering mueller matrix for highly scattering media

We have developed a Monte Carlo algorithm that computes all two-dimensional elements of the diffuse backscattering Mueller matrix for highly scattering media. Using the Stokes-Mueller formalism and scattering amplitudes calculated with Mie theory, we are able to consider polarization-dependent photon propagation in highly scattering media, including linearly and circularly polarized light. The numerically determined matrix elements are compared with experimental data for different particle sizes and show good agreement in both azimuthal and radial direction.     

Pitfalls and modelling inconsistencies in computational radiation dosimetry: lessons learnt from the QUADOS intercomparison. Part II: Photons, electrons and protons

'QUADOS', a concerted action of the European Commission, has promoted an intercomparison aimed at evaluating the use of computational codes for dosimetry in radiation protection and medical physics. This intercomparison was open to all users of radiation transport codes. Eight problems were selected for their relevance to the radiation dosimetry community, five of which involved photon and proton transport. This paper focuses on a discussion of lessons learned from the participation in solving the photon and charged particle problems. The lessons learned from the participation in solving the neutron problems are presented in a companion paper (in this issue).     

Estimation of scattering error in spectrophotometric measurements of light absorption by aquatic particles from three-dimensional radiative transfer simulations

We present an approach based on three-dimensional Monte Carlo radiative transfer simulations for estimating scattering error in measurements of light absorption by aquatic particles with a typical laboratory double-beam spectrophotometer. The scattering error is calculated by combining the weighting function describing the angular distribution of photon losses that are due to scattering on suspended particles with the volume scattering function of particles. We applied this method to absorption measurements made on marine phytoplankton, a diatom Thalassiosira pseudonana and a cyanobacterium Synechococcus. Assuming that the scattering phase function is described by the Henyey-Greenstein formula, we determined the backscatter probability of phytoplankton, which yields the best correction for scattering error at a light wavelength of 750 nm, where true absorption is null. The backscattering ratio estimated for both phytoplankton species is significantly higher than previously reported data based on Mie-scattering calculations for homogeneous spheres. Depending on the type of particles, the corrected absorption spectra obtained with our method may be similar or significantly different from spectra obtained with the null-point correction based on wavelength-independent scattering error.     

Validating the assumption to the interference approximation by use of measurements of absorption efficiency and hindered scattering in dense suspensions

Frequency domain photon migration (FDPM) measurements were employed to accurately quantify optical properties of both the suspending fluid and particles within dense polystyrene suspensions of 143- or 226-nm mean diameter at varying concentrations (5-30% by volume). The measured absorption coefficients varied linearly with particle volume fraction whereas the isotropic scattering coefficients varied nonlinearly in agreement with the prediction that utilizes the hard-sphere structure factor model. These results validate the interference approximation of light scattering to describe light propagation accurately within dense suspensions. Furthermore, owing to the accuracy of FDPM absorption measurements, the imaginary refractive indices for both particles and their suspending fluid were determined and were found to compare favorably with literature values.     

Probabilistic model of multiple light scattering based on rigorous computation of the first and the second moments of photon coordinates

We consider a concise method based on recurrent relations that permit rigorous computing of the first and the second moments of the components of the vector locating a randomly walking photon in an infinite homogeneous light-scattering medium. On assumption that the components obey a three-dimensional Gaussian distribution a probability density for the photon locations at the Nth scattering event can readily be written down and the light-intensity distribution in the medium may be calculated. The results from theoretical analyses are compared with the solution of a light-diffusion equation and with results of Monte Carlo simulations and show a better fit with simulated data than the diffusion approximation.     

Optimal probe geometry for near-infrared spectroscopy of biological tissue

The results of this study clarify the influence of probe geometry on spectroscopic measurements obtained from the surface of a turbid biological tissue. We show that the transition between the measurement of the predominantly backward-propagating and the predominantly forward-propagating photon fluxes is marked by the separation between the source probe and the detector probes at which the dependence of the fluence on small changes in scattering coefficient vanishes. This is the probe separation at which a variable scattering background has the least influence on the measurement of optical absorption in turbid materials. Estimates of the optimum probe spacing for typical values of absorption and scattering coefficients of soft tissue in the near-infrared spectral region (800-2500 nm) are derived from an analytical solution of the diffusion equation. The estimates were verified by Monte Carlo simulations and experiments on particle suspensions with optical properties similar to those of skin tissue.     

基于FPGA的多通道动态光散射自相关器的研制

为提高多角度动态光散射(MDLS)实验的测试效率与测量精度,开发了一种基于现场可编程门阵列(FPGA)的多通道动态光散射自相关器。该自相关器利用FPGA同时采集4个散射角度的光子脉冲信号,实现高频脉冲计数和自相关计算;设计双计数器模块保证高频无丢失计数,结合片内环形寄存区与片外DDR3芯片实现动态存储数据;提出multi-tau相关器结构设计自相关器的相关计算模块,通过USB通信模块传输自相关计算结果。该自相关器利用有限的FPGA内部硬件资源,实现了多通道自相关运算,解决了大量计数数据的存储问题,提高了MDLS实验的测试效率。利用自相关器对50～200 nm的聚苯乙烯颗粒进行MDLS实验,实验结果与样品的标称粒径表现出良好的一致性,表明该设计能有效地表征颗粒粒径信息。     

Spectral model of a fluorescent ink halftone

A model is presented of a fluorescent ink halftone. Unlike a nonfluorescent ink, which only absorbs, a fluorescent ink absorbs higher-energy photons and emits lower-energy photons. The amount of fluorescent light produced depends on the percent absorption of the incident light. For fluorescent ink printed on paper, both photon scattering within the paper substrate and multiple internal reflections between the ink layer and the paper substrate significantly increase the percent absorption, so a realistic model must include these effects. The model presented here utilizes the generalized Clapper-Yule theory, which accounts for photon diffusion that is due to both scatter and internal reflection. It is shown that while multiple internal reflections alone only marginally increase the percent absorption, when there are both scattering and internal reflection, the percent absorption is increased significantly. The current study is a theoretical model and does not present experimental results.     

Monte Carlo Modeling of Phonons at Crystal Interfaces

One of the strategies in the effort to directly detect dark matter particles is to measure the phonon and charge signals produced in an incidental collision of a weakly interacting massive particle with a nucleon in a crystalline detector. Proper calibration of the detected phonon energy relies on detailed models of phonon propagation through the crystal to the instrumented surface. Previous Detector Monte Carlo incorporate probabalistic anharmonic decay and mass defect scattering but neglect the mode dependent scattering at crystal boundaries. We calculate mode specific reflection and transmission at detector interfaces with a simple acoustic mismatch model. We find that mode preserving transmission is the most probable outcome at germanium-aluminum detector interfaces, but the probability of reflection is not negligible. The average phonon reflection probability at near normal angles of incidence at a Ge/Al interface is near 20 %, but grows dramatically for oblique incidence. We develop a code using Geant4, which should allow modeling extensions to all phonon mediated dark matter detection schemes. Our models are adaptable to other crystal materials and are generally useful in any phonon interface problem.     

TURBIDIMETRIC SPECTRUM METHOD FOR FINE PARTICLE SIZE ANALYSIS

Based upon light scattering techniques, a turbidimetric spectrum method and corresponding inversion algorithm for particle size characterization are presented. Comprehensive theoretical analysis and computer simulation show that the presented method is capable of determining not only average particle diameter, but also particle size distribution in a very broad size range. In this paper the emphasis is on and limited to this method's applicability to fine particle size analysis, ranging from tens of nanometers to several micrometers. Experimental studies are carried out with submicron monodispersed lattices. The results of comparative measurements using the photon correlation spectrum (PCS) and presented turbidimetric spectrum method are also given. It shows that the agreement is good.