Single-resolution and multiresolution extended-Kalman-filter-based reconstruction approaches to optical refraction tomography

The problem of reconstruction of a refractive-index distribution (RID) in optical refraction tomography (ORT) with optical path-length difference (OPD) data is solved using two adaptive-estimation-based extended-Kalman-filter (EKF) approaches. First, a basic single-resolution EKF (SR-EKF) is applied to a state variable model describing the tomographic process, to estimate the RID of an optically transparent refracting object from noisy OPD data. The initialization of the biases and covariances corresponding to the state and measurement noise is discussed. The state and measurement noise biases and covariances are adaptively estimated. An EKF is then applied to the wavelet-transformed state variable model to yield a wavelet-based multiresolution EKF (MR-EKF) solution approach. To numerically validate the adaptive EKF approaches, we evaluate them with benchmark studies of standard stationary cases, where comparative results with commonly used efficient deterministic approaches can be obtained. Detailed reconstruction studies for the SR-EKF and two versions of the MR-EKF (with Haar and Daubechies-4 wavelets) compare well with those obtained from a typically used variant of the (deterministic) algebraic reconstruction technique, the average correction per projection method, thus establishing the capability of the EKF for ORT. To the best of our knowledge, the present work contains unique reconstruction studies encompassing the use of EKF for ORT in single-resolution and multiresolution formulations, and also in the use of adaptive estimation of the EKF's noise covariances. (C) 2010 Optical Society of America

Quantitative photoacoustic tomography from boundary pressure measurements: noniterative recovery of optical absorption coefficient from the reconstructed absorbed energy map

We describe a noniterative method for recovering optical absorption coefficient distribution from the absorbed energy map reconstructed using simulated and noisy boundary pressure measurements. The source reconstruction problem is first solved for the absorbed energy map corresponding to single- and multiple-source illuminations from the side of the imaging plane. It is shown that the absorbed energy map and the absorption coefficient distribution, recovered from the single-source illumination with a large variation in photon flux distribution, have signal-to-noise ratios comparable to those of the reconstructed parameters from a more uniform photon density distribution corresponding to multiple-source illuminations. The absorbed energy map is input as absorption coefficient times photon flux in the time-independent diffusion equation (DE) governing photon transport to recover the photon flux in a single step. The recovered photon flux is used to compute the optical absorption coefficient distribution from the absorbed energy map. In the absence of experimental data, we obtain the boundary measurements through Monte Carlo simulations, and we attempt to address the possible limitations of the DE model in the overall reconstruction procedure.

Quasiequilibrium optical nonlinearities from spin-polarized carriers in GaAs

Semiconductor Bloch equations, which microscopically describe the dynamics of a Coulomb interacting, spin-unpolarized electron-hole plasma, can be solved in two limits: the coherent and the quasiequilibrium regimes. These equations have been recently extended to include the spin degree of freedom and used to explain spin dynamics in the coherent regime. In the quasiequilibrium limit, one solves the Bethe-Salpeter equation in a two-band model to describe how optical absorption is affected by Coulomb interactions within a spin unpolarized plasma of arbitrary density. In this work, we modified the solution of the Bethe-Salpeter equation to include spin polarization and light holes in a three-band model, which allowed us to account for spin-polarized versions of many-body effects in absorption. The calculated absorption reproduced the spin-dependent, density-dependent, and spectral trends observed in bulk GaAs at room temperature, in a recent pump-probe experiment with circularly polarized light. Hence, our results may be useful in the microscopic modeling of density-dependent optical nonlinearities due to spin-polarized carriers in semiconductors.

Structural, ferroelectric and optical properties of Bi2VO5.5 thin films deposited on platinized silicon {(100) Pt/TiO2/SiO2/Si} substrates

Bismuth vanadate (Bi2VO5.5, BVO) thin films have been deposited by a pulsed laser ablation technique on platinized silicon substrates. The surface morphology of the BVO thin films has been studied by atomic force microscopy (AFM). The optical properties of the BVO thin films were investigated using spectroscopic ellipsometric measurements in the 300–820 nm wavelength range. The refractive index (n), extinction coefficient (k) and thickness of the BVO thin films have been obtained by fitting the ellipsometric experimental data in a four-phase model (air/BVOrough/BVO/Pt). The values of the optical constants n and k that were determined through multilayer analysis at 600 nm were 2.31 and 0.056, respectively. For fitting the ellipsometric data and to interpret the optical constants, the unknown dielectric function of the BVO films was constructed using a Lorentz model. The roughness of the films was modeled in the Brugmann effective medium approximation and the results were compared with the AFM observations.

Optical Studies of the Antarctic Glacio-Oceanic System

The Antarctic system comprises of the continent itself, Antarctica, and the ocean surrounding it, the Southern Ocean. The system has an important part in the global climate due to its size, its high latitude location and the negative radiation balance of its large ice sheets. Antarctica has also been in focus for several decades due to increased ultraviolet (UV) levels caused by stratospheric ozone depletion, and the disintegration of its ice shelves. In this study, measurements were made during three Austral summers to study the optical properties of the Antarctic system and to produce radiation information for additional modeling studies. These are related to specific phenomena found in the system. During the summer of 1997-1998, measurements of beam absorption and beam attenuation coefficients, and downwelling and upwelling irradiance were made in the Southern Ocean along a S-N transect at 6°E. The attenuation of photosynthetically active radiation (PAR) was calculated and used together with hydrographic measurements to judge whether the phytoplankton in the investigated areas of the Southern Ocean are light limited. By using the Kirk formula the diffuse attenuation coefficient was linked to the absorption and scattering coefficients. The diffuse attenuation coefficients (Kpar) for PAR were found to vary between 0.03 and 0.09 1/m. Using the values for KPAR and the definition of the Sverdrup critical depth, the studied Southern Ocean plankton systems were found not to be light limited. Variabilities in the spectral and total albedo of snow were studied in the Queen Maud Land region of Antarctica during the summers of 1999-2000 and 2000-2001. The measurement areas were the vicinity of the South African Antarctic research station SANAE 4, and a traverse near the Finnish Antarctic research station Aboa. The midday mean total albedos for snow were between 0.83, for clear skies, and 0.86, for overcast skies, at Aboa and between 0.81 and 0.83 for SANAE 4. The mean spectral albedo levels at Aboa and SANAE 4 were very close to each other. The variations in the spectral albedos were due more to differences in ambient conditions than variations in snow properties. A Monte-Carlo model was developed to study the spectral albedo and to develop a novel nondestructive method to measure the diffuse attenuation coefficient of snow. The method was based on the decay of upwelling radiation moving horizontally away from a source of downwelling light. This was assumed to have a relation to the diffuse attenuation coefficient. In the model, the attenuation coefficient obtained from the upwelling irradiance was higher than that obtained using vertical profiles of downwelling irradiance. The model results were compared to field measurements made on dry snow in Finnish Lapland and they correlated reasonably well. Low-elevation (below 1000 m) blue-ice areas may experience substantial melt-freeze cycles due to absorbed solar radiation and the small heat conductivity in the ice. A two-dimensional (x-z) model has been developed to simulate the formation and water circulation in the subsurface ponds. The model results show that for a physically reasonable parameter set the formation of liquid water within the ice can be reproduced. The results however are sensitive to the chosen parameter values, and their exact values are not well known. Vertical convection and a weak overturning circulation is generated stratifying the fluid and transporting warmer water downward, thereby causing additional melting at the base of the pond. In a 50-year integration, a global warming scenario mimicked by a decadal scale increase of 3 degrees per 100 years in air temperature, leads to a general increase in subsurface water volume. The ice did not disintegrate due to the air temperature increase after the 50 year integration.

The third post-Newtonian gravitational wave polarizations and associated spherical harmonic modes for inspiralling compact binaries in quasi-circular orbits

The gravitational waveform (GWF) generated by inspiralling compact binaries moving in quasi-circular orbits is computed at the third post-Newtonian (3PN) approximation to general relativity. Our motivation is two-fold: (i) to provide accurate templates for the data analysis of gravitational wave inspiral signals in laser interferometric detectors; (ii) to provide the associated spin-weighted spherical harmonic decomposition to facilitate comparison and match of the high post-Newtonian prediction for the inspiral waveform to the numerically-generated waveforms for the merger and ringdown. This extension of the GWF by half a PN order (with respect to previous work at 2.5PN order) is based on the algorithm of the multipolar post-Minkowskian formalism, and mandates the computation of the relations between the radiative, canonical and source multipole moments for general sources at 3PN order. We also obtain the 3PN extension of the source multipole moments in the case of compact binaries, and compute the contributions of hereditary terms (tails, tails-of-tails and memory integrals) up to 3PN order. The end results are given for both the complete plus and cross polarizations and the separate spin-weighted spherical harmonic modes.

The partition of unity finite element method for elastic wave propagation in Reissner-Mindlin plates

This paper reports a numerical method for modelling the elastic wave propagation in plates. The method is based on the partition of unity approach, in which the approximate spectral properties of the infinite dimensional system are embedded within the space of a conventional finite element method through a consistent technique of waveform enrichment. The technique is general, such that it can be applied to the Lagrangian family of finite elements with specific waveform enrichment schemes, depending on the dominant modes of wave propagation in the physical system. A four-noded element for the Reissner-indlin plate is derived in this paper, which is free of shear locking. Such a locking-free property is achieved by removing the transverse displacement degrees of freedom from the element nodal variables and by recovering the same through a line integral and a weak constraint in the frequency domain. As a result, the frequency-dependent stiffness matrix and the mass matrix are obtained, which capture the higher frequency response with even coarse meshes, accurately. The steps involved in the numerical implementation of such element are discussed in details. Numerical studies on the performance of the proposed element are reported by considering a number of cases, which show very good accuracy and low computational cost. Copyright (C)006 John Wiley & Sons, Ltd.

Formal Verification of Full-Wave Rectifier: A Case Study

We present a case study of formal verification of full-wave rectifier for analog and mixed signal designs. We have used the Checkmate tool from CMU [1], which is a public domain formal verification tool for hybrid systems. Due to the restriction imposed by Checkmate it necessitates to make the changes in the Checkmate implementation to implement the complex and non-linear system. Full-wave rectifier has been implemented by using the Checkmate custom blocks and the Simulink blocks from MATLAB from Math works. After establishing the required changes in the Checkmate implementation we are able to efficiently verify, the safety properties of the full-wave rectifier.

Spectroscopic ellipsometry investigations of the optical properties of manganese doped bismuth vanadate thin films

The optical properties of Bi(2)V(1-x)MnxO(5.5-x) (x=0.05, 0.1, 0.15 and 0.2 at.%) thin films fabricated by pulsed laser deposition on platinized Silicon Substrates were Studied in UV-visible spectral region (1.51-4.17 CV) using spectroscopic ellipsometry. The optical constants and thicknesses of these films have been obtained by fitting the ellipsometric data (Psi and Delta) using a multilayer four-phase model system and a relaxed Lorentz oscillator dispersion relation. The surface roughness and film thickness obtained by spectroscopic ellipsometry were found to be consistent with the results obtained by atomic force and scanning electron microscopy. The refractive index measured at 650 nm does not show any marginal increase with Mn content. Further, the extinction coefficient does not show much decrease with increasing Mn content. An increase in optical band gap energy from 2.52 to 2.77 eV with increasing Mn Content from x = 0.05 to 0.15 was attributed to the increase in oxygen ion vacancy disorder. (C) 2009 Elsevier Ltd. All rights reserved.

Simplified theory of optical nonlinearities in spin-polarized bulk GaAs

Apart from their intrinsic physical interest, spin-polarized many-body effects are expected to be important to the working of spintronic devices. A vast literature exists on the effects of a spin-unpolarized electron-hole plasma on the optical properties of a semiconductor. Here, we include the spin degree of freedom to model optical absorption of circularly polarized light by spin-polarized bulk GaAs. Our model is easy to implement and does not require elaborate numerics, since it is based on the closed-form analytical pair-equation formula that is valid in 3d. The efficacy of our approach is demonstrated by a comparison with recent experimental data.

Structure and optical properties of Cd-substituted ZnO (Zn1-xCdxO) nanostructures synthesized by the high-pressure solution route

We report the synthesis of Cd-substituted ZnO nanostructures (Zn1-xCdxO with x up to approximate to 0.09) by the high-pressure solution growth method. The synthesized nanostructures comprise nanocrystals that are both particles (similar to 10-15 nm) and rods which grow along the [002] direction as established by transmission electron microscope (TEM) and x-ray diffraction (XRD) analysis. Rietveld analysis of the XRD data shows a monotonic increase of the unit cell volume with the increase of Cd concentration. The optical absorption, as well as the photoluminescence (PL), shows a red shift on Cd substitution. The line width of the PL spectrum is related to the strain inhomogeneity and it peaks in the region where the CdO phase separates from the Zn1-xCdxO nanostructures. The time-resolved photoemission showed a long-lived (similar to 10 ns) component. We propose that the PL behaviour of the Zn1-xCdxO is dominated by strain in the sample with the red shift of the PL linked to the expansion of the unit cell volume on Cd substitution.