Advances in characterisation, calibration and data processing speed of optical coherence tomography systems

This thesis describes advances in the characterisation, calibration and data processing of optical coherence tomography (OCT) systems. Femtosecond (fs) laser inscription was used for producing OCT-phantoms. Transparent materials are generally inert to infra-red radiations, but with fs lasers material modification occurs via non-linear processes when the highly focused light source interacts with the materials. This modification is confined to the focal volume and is highly reproducible. In order to select the best inscription parameters, combination of different inscription parameters were tested, using three fs laser systems, with different operating properties, on a variety of materials. This facilitated the understanding of the key characteristics of the produced structures with the aim of producing viable OCT-phantoms. Finally, OCT-phantoms were successfully designed and fabricated in fused silica. The use of these phantoms to characterise many properties (resolution, distortion, sensitivity decay, scan linearity) of an OCT system was demonstrated. Quantitative methods were developed to support the characterisation of an OCT system collecting images from phantoms and also to improve the quality of the OCT images. Characterisation methods include the measurement of the spatially variant resolution (point spread function (PSF) and modulation transfer function (MTF)), sensitivity and distortion. Processing of OCT data is a computer intensive process. Standard central processing unit (CPU) based processing might take several minutes to a few hours to process acquired data, thus data processing is a significant bottleneck. An alternative choice is to use expensive hardware-based processing such as field programmable gate arrays (FPGAs). However, recently graphics processing unit (GPU) based data processing methods have been developed to minimize this data processing and rendering time. These processing techniques include standard-processing methods which includes a set of algorithms to process the raw data (interference) obtained by the detector and generate A-scans. The work presented here describes accelerated data processing and post processing techniques for OCT systems. The GPU based processing developed, during the PhD, was later implemented into a custom built Fourier domain optical coherence tomography (FD-OCT) system. This system currently processes and renders data in real time. Processing throughput of this system is currently limited by the camera capture rate. OCTphantoms have been heavily used for the qualitative characterization and adjustment/ fine tuning of the operating conditions of OCT system. Currently, investigations are under way to characterize OCT systems using our phantoms. The work presented in this thesis demonstrate several novel techniques of fabricating OCT-phantoms and accelerating OCT data processing using GPUs. In the process of developing phantoms and quantitative methods, a thorough understanding and practical knowledge of OCT and fs laser processing systems was developed. This understanding leads to several novel pieces of research that are not only relevant to OCT but have broader importance. For example, extensive understanding of the properties of fs inscribed structures will be useful in other photonic application such as making of phase mask, wave guides and microfluidic channels. Acceleration of data processing with GPUs is also useful in other fields.

Stokes waves revisited:exact solutions in the asymptotic limit

The Stokes perturbative solution of the nonlinear (boundary value dependent) surface gravity wave problem is known to provide results of reasonable accuracy to engineers in estimating the phase speed and amplitudes of such nonlinear waves. The weakling in this structure though is the presence of aperiodic “secular variation” in the solution that does not agree with the known periodic propagation of surface waves. This has historically necessitated increasingly higher-ordered (perturbative) approximations in the representation of the velocity profile. The present article ameliorates this long-standing theoretical insufficiency by invoking a compact exact n-ordered solution in the asymptotic infinite depth limit, primarily based on a representation structured around the third-ordered perturbative solution, that leads to a seamless extension to higher-order (e.g., fifth-order) forms existing in the literature. The result from this study is expected to improve phenomenological engineering estimates, now that any desired higher-ordered expansion may be compacted within the same representation, but without any aperiodicity in the spectral pattern of the wave guides.

Wave trapping in a two-dimensional sound-soft or sound-hard acoustic waveguide of slowly-varying width

In this paper we derive novel approximations to trapped waves in a two-dimensional acoustic waveguide whose walls vary slowly along the guide, and at which either Dirichlet (sound-soft) or Neumann (sound-hard) conditions are imposed. The guide contains a single smoothly bulging region of arbitrary amplitude, but is otherwise straight, and the modes are trapped within this localised increase in width. Using a similar approach to that in Rienstra (2003), a WKBJ-type expansion yields an approximate expression for the modes which can be present, which display either propagating or evanescent behaviour; matched asymptotic expansions are then used to derive connection formulae which bridge the gap across the cut-off between propagating and evanescent solutions in a tapering waveguide. A uniform expansion is then determined, and it is shown that appropriate zeros of this expansion correspond to trapped mode wavenumbers; the trapped modes themselves are then approximated by the uniform expansion. Numerical results determined via a standard iterative method are then compared to results of the full linear problem calculated using a spectral method, and the two are shown to be in excellent agreement, even when $\epsilon$, the parameter characterising the slow variations of the guide’s walls, is relatively large.

T-wave alternans: desempenho clínico, limitações, metodologias de análise

Reconhecer com precisão indivíduos com maior risco imediato de morte súbita cardíaca (MSC) ainda é uma questão em aberto. A natureza fortuita dos eventos cardiovasculares agudos não parece se adequar ao conhecido modelo de indução de taquicardia/fibrilação ventricular por um gatilho em sincronia a um substrato arritmogênico estático. Quanto ao mecanismo da MSC, uma instabilidade elétrica dinâmica explicaria melhor a raridade da associação simultânea de um gatilho certo a um substrato cardíaco apropriado. Diversos estudos tentaram medir essa instabilidade elétrica cardíaca (ou um equivalente válido) em uma sequência de batimentos cardíacos no ECG. Dentre os mecanismos possíveis podemos citar o prolongamento do QT, dispersão do QT, potenciais tardios, alternância de onda T ou T-wave alternans (TWA), e turbulência da frequência cardíaca. Este artigo se atém em particular ao papel da TWA no panorama atual da estratificação de risco cardíaco. Os achados sobre TWA ainda são heterogêneos, variando de um desempenho prognóstico muito bom até um quase nulo, dependendo da população clínica observada e protocolo clínico usado. Para preencher as atuais lacunas no conhecimento sobre TWA, profissionais médicos e pesquisadores devem explorar melhor as características técnicas das diversas tecnologias disponíveis para a avaliação de TWA e atentar ao fato de que os valores de TWA respondem a diversos outros fatores, além de medicamentos. Informações sobre mecanismos celulares e subcelulares da TWA estão fora do escopo deste artigo, mas são referenciados alguns dos principais trabalhos sobre este tópico, com o intuito de auxiliar no entendimento dos conceitos e fatos cobertos neste artigo.

Determination of 5-aminosalicylic acid in pharmaceutical formulations by square wave voltammetry at pencil graphite electrodes

An analytical method for the determination of the anti-inflammatory drug 5-aminosalicylic acid (5-ASA) in pharmaceutical formulations using square wave voltammetry at pencil graphite electrodes was developed. After the optimization of the experimental conditions, calibration curves were obtained in the linear concentration range from 9.78 × 10-7 to 7.25 × 10-5 mol L-1 resulting in a limit of detection of 2.12 ± 0.05 x 10-8 mol L-1. Statistical tests showed that the concentrations of 5-ASA in commercial tablets and enemas obtained with the proposed voltammetric method agreed with HPLC values at a 95% confidence level.

Brazilian offshore wave climate based on NWW3 reanalysis

This paper provides a description of the wave climate off the Brazilian coast based on an eleven-year time series (Jan/1997-Dec/2007) obtained from the NWW3 operational model hindcast reanalysis. Information about wave climate in Brazilian waters is very scarce and mainly based on occasional short-term observations, the present analysis being the first covering such temporal and spatial scales. To define the wave climate, six sectors were defined and analyzed along the Brazilian shelf-break: South (W1), Southeast (W2), Central (W3), East (W4), Northeast (W5) and North (W6). W1, W2 and W3 wave regimes are determined by the South Atlantic High (SAH) and the passage of synoptic cold fronts; W4, W5 and W6 are controlled by the Intertropical Convergence Zone (ITCZ) and its meridional oscillation. The most energetic waves are from the S, generated by the strong winds associated to the passage of cold fronts, which mainly affect the southern region. Wave power presents a decrease in energy levels from south to north, with its annual variation showing that the winter months are the most energetic in W1 to W4, while in W5 and W6 the most energetic conditions occur during the austral summer. The information presented here provides boundary conditions for studies related to coastal processes, fundamental for a better understanding of the Brazilian coastal zone.

Evaluation of left ventricular diastolic echocardiographic parameters in healthy dogs by pulsed-wave Doppler

Left ventricular diastolic dysfunction plays an important role on heart failure progression. In order to obtain additional reference values of left ventricular diastolic parameters and investigate influence of common variables, peak E wave (peak E), peak A wave (peak A), E/A ratio (E/A), E wave deceleration time (EDT) and isovolumic relaxation time (IRVT) were studied in 40 clinically healthy dogs, by pulsed wave Doppler. The following values were obtained: peak E = 0.747 ± 0.117 m/s, peak A = 0.487 ± 0.062 m/s, E/A = 1.533 ± 0.198, EDT = 88.7 ± 9.2 ms and IRVT = 0.080 ± 0.009 s. Some parameters were influenced by heart rate (peak E, peak A and IRVT), by age (peak A and E/A) and by body weight (TRIV). Gender influence was absent. Values obtained can be used as reference for canine specimens but its interpretation should consider on the influence of related variables.

Matter-wave two-dimensional solitons in crossed linear and nonlinear optical lattices

The existence of multidimensional matter-wave solitons in a crossed optical lattice (OL) with a linear optical lattice (LOL) in the x direction and a nonlinear optical lattice (NOL) in the y direction, where the NOL can be generated by a periodic spatial modulation of the scattering length using an optically induced Feshbach resonance is demonstrated. In particular, we show that such crossed LOLs and NOLs allow for stabilizing two-dimensional solitons against decay or collapse for both attractive and repulsive interactions. The solutions for the soliton stability are investigated analytically, by using a multi-Gaussian variational approach, with the Vakhitov-Kolokolov necessary criterion for stability; and numerically, by using the relaxation method and direct numerical time integrations of the Gross-Pitaevskii equation. Very good agreement of the results corresponding to both treatments is observed.

Spin-wave fluctuations in ferrimagnetic Mg(x)Fe(3-x)O(4) nanoparticles

We have performed a systematic study of the magnetic properties of a series of ferrimagnetic nanoparticles of Mg(x)Fe(3-x)O(4) (0.8 <= x <= 1.5) prepared by the combustion reaction method. The magnetization data can be well fitted by Bloch's law with T(3/2). Bloch's constant B determined from the fitting procedure was found to increase with Mg content x from similar to 3.09 X 10(-5) K(-3/2) for x = 0.8 to 6.27 X 10(-5) K(-3/2) for x=1.5. The exchange integral J(AB) and the spin-wave stiffness constant D of Mg(x)Fe(3-x)O(4) nanoparticles were also determined as similar to 0.842 and 0.574 meV and 296 and 202 meV angstrom(2) for specimens with x=0.8 and 1.5, respectively. These results are discussed in terms of cation redistribution among A and B sites on these nanostructured spinel ferrites. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3359709]

Nonmonotonic thickness dependence of spin wave energy in ultrathin Fe films: Experiment and theory

High wave-vector spin waves in ultrathin Fe/W(110) films up to 20 monolayers (MLs) thick have been studied using spin-polarized electron energy-loss spectroscopy. An unusual nonmonotonous dependence of the spin wave energies on the film thickness is observed, featuring a pronounced maximum at 2 ML coverage. First-principles theoretical study reveals the origin of this behavior to be in the localization of the spin waves at the surface of the film, as well as in the properties of the interlayer exchange coupling influenced by the hybridization of the electron states of the film and substrate and by the strain.

Electromagnetic Wave Scattering by Schwarzschild Black Holes

We analyze the scattering of a planar monochromatic electromagnetic wave incident upon a Schwarzschild black hole. We obtain accurate numerical results from the partial wave method for the electromagnetic scattering cross section and show that they are in excellent agreement with analytical approximations. The scattering of electromagnetic waves is compared with the scattering of scalar, spinor, and gravitational waves. We present a unified picture of the scattering of all massless fields for the first time.

Scalar resonances in a unitary pi pi S-wave model for D(+)->pi(+)pi(-)pi(+)

We propose a model for D(+)->pi(+)pi(-)pi(+) decays following experimental results which indicate that the two-pion interaction in the S wave is dominated by the scalar resonances f(0)(600)/sigma and f(0)(980). The weak decay amplitude for D(+)-> R pi(+), where R is a resonance that subsequently decays into pi(+)pi(-), is constructed in a factorization approach. In the S wave, we implement the strong decay R ->pi(+)pi(-) by means of a scalar form factor. This provides a unitary description of the pion-pion interaction in the entire kinematically allowed mass range m(pi pi)(2) from threshold to about 3 GeV(2). In order to reproduce the experimental Dalitz plot for D(+)->pi(+)pi(-)pi(+), we include contributions beyond the S wave. For the P wave, dominated by the rho(770)(0), we use a Breit-Wigner description. Higher waves are accounted for by using the usual isobar prescription for the f(2)(1270) and rho(1450)(0). The major achievement is a good reproduction of the experimental m(pi pi)(2) distribution, and of the partial as well as the total D(+)->pi(+)pi(-)pi(+) branching ratios. Our values are generally smaller than the experimental ones. We discuss this shortcoming and, as a by-product, we predict a value for the poorly known D ->sigma transition form factor at q(2)=m pi(2).

Gravitational wave recoil in Robinson-Trautman spacetimes

We consider the gravitational recoil due to nonreflection-symmetric gravitational wave emission in the context of axisymmetric Robinson-Trautman spacetimes. We show that regular initial data evolve generically into a final configuration corresponding to a Schwarzschild black hole moving with constant speed. For the case of (reflection-)symmetric initial configurations, the mass of the remnant black hole and the total energy radiated away are completely determined by the initial data, allowing us to obtain analytical expressions for some recent numerical results that have appeared in the literature. Moreover, by using the Galerkin spectral method to analyze the nonlinear regime of the Robinson-Trautman equations, we show that the recoil velocity can be estimated with good accuracy from some asymmetry measures (namely the first odd moments) of the initial data. The extension for the nonaxisymmetric case and the implications of our results for realistic situations involving head-on collision of two black holes are also discussed.

Radio-frequency spectroscopy of (6)Li p-wave molecules: Towards photoemission spectroscopy of a p-wave superfluid

We study rf spectroscopy of a lithium gas with the goal to explore the possibilities for photoemission spectroscopy of a strongly interacting p-wave Fermi gas. Radio-frequency spectra of quasibound p-wave molecules and of free atoms in the vicinity of the p-wave Feshbach resonance located at 159.15G are presented. The spectra are free of detrimental final-state effects. The observed relative magnetic-field shifts of the molecular and atomic resonances confirm earlier measurements realized with direct rf association. Furthermore, evidence of molecule production by adiabatically ramping the magnetic field is observed. Finally, we propose the use of a one-dimensional optical lattice to study anisotropic superfluid gaps as most direct proof of p-wave superfluidity.

Quantitative determination of optical transmission through subwavelength slit arrays in Ag films: Role of surface wave interference and local coupling between adjacent slits

Measurement of the transmitted intensity from a coherent monomode light source through a series of subwavelength slit arrays in Ag films, with varying array pitch and number of slits, demonstrates enhancement (suppression) by factors of as much as 6 (9) when normalized to the transmission efficiency of an isolated slit. Pronounced minima in the transmitted intensity are observed at array pitches corresponding to lambda(SPP), 2 lambda(SPP), and 3 lambda(SPP), where lambda(SPP) is the wavelength of the surface plasmon polariton (SPP). The position of these minima arises from destructive interference between incident propagating waves and pi-phase-shifted SPP waves. Increasing the number of slits to four or more does not increase appreciably the per-slit transmission intensity. A simple interference model fits well the measured transmitted intensity profile.