Coarse-grid higher order finite-difference time-domain algorithm with low dispersion errors

Higher order (2,4) FDTD schemes used for numerical solutions of Maxwell`s equations are focused on diminishing the truncation errors caused by the Taylor series expansion of the spatial derivatives. These schemes use a larger computational stencil, which generally makes use of the two constant coefficients, C-1 and C-2, for the four-point central-difference operators. In this paper we propose a novel way to diminish these truncation errors, in order to obtain more accurate numerical solutions of Maxwell`s equations. For such purpose, we present a method to individually optimize the pair of coefficients, C-1 and C-2, based on any desired grid size resolution and size of time step. Particularly, we are interested in using coarser grid discretizations to be able to simulate electrically large domains. The results of our optimization algorithm show a significant reduction in dispersion error and numerical anisotropy for all modeled grid size resolutions. Numerical simulations of free-space propagation verifies the very promising theoretical results. The model is also shown to perform well in more complex, realistic scenarios.

Comparison of Fourier-Domain and Time-Domain Optical Coherence Tomography in the Detection of Band Atrophy of the Optic Nerve

PURPOSE: To compare the ability of Fourier-domain (FD) optical coherence tomography (3D OCT-1000; Top, con, Tokyo, Japan) and time domain (TD) OCT (Stratus; Carl Zeiss Meditec Inc, Dublin, California, USA) to detect axonal loss in eyes with band atrophy (BA) of the optic nerve. DESIGN: Cross-sectional study. METHODS: Thirty-six eyes from 36 patients with BA and temporal visual field (VF) defect from chiasmal compression and 36 normal eyes were studied. Subjects were submitted to standard automated perimetry and macular and retinal nerve fiber layer (RNFL) measurements were taken using 3D OCT-1000 and Stratus OCT. Receiver operating characteristic (ROC) curves were calculated for each parameter. Spearman correlation coefficients were obtained to evaluate the relationship between RNFL and macular thickness parameters and severity of VF loss. Measurements from the two devices were compared. RESULTS: Regardless of OCT device, all RNFL and macular thickness parameters were significantly lower in eyes with BA compared with normal eyes, but no statistically significant difference was found with regard to the area under the ROC curve. Structure-function relationships were also similar for the two devices. In both groups, RNFL and macular thickness measurements were generally and in some cases significantly smaller with 3D OCT-1000 than with Stratus OCT. CONCLUSIONS: The introduction of FD technology did not lead to better discrimination ability for detecting BA of the optic nerve compared with TD technology when using the software currently provided by the manufacturer. 3D OCT-1000 FD OCT RNFL and macular measurements were generally smaller than TD Stratus OCT measurements. Investigators should be aware of this fact when comparing measurements obtained with these two devices. (Am J Oplathalmol 2009;147: 56-63. (c) 2009 by Elsevier Inc. All rights reserved.)

A Study on Unconditionally Stable FDTD Methods for the Modeling of Metamaterials

We assess the performance of three unconditionally stable finite-difference time-domain (FDTD) methods for the modeling of doubly dispersive metamaterials: 1) locally one-dimensional FDTD; 2) locally one-dimensional FDTD with Strang splitting; and (3) alternating direction implicit FDTD. We use both double-negative media and zero-index media as benchmarks.

Dynamic Imaging in Electrical Impedance Tomography of the Human Chest With Online Transition Matrix Identification

One of the electrical impedance tomography objectives is to estimate the electrical resistivity distribution in a domain based only on electrical potential measurements at its boundary generated by an imposed electrical current distribution into the boundary. One of the methods used in dynamic estimation is the Kalman filter. In biomedical applications, the random walk model is frequently used as evolution model and, under this conditions, poor tracking ability of the extended Kalman filter (EKF) is achieved. An analytically developed evolution model is not feasible at this moment. The paper investigates the identification of the evolution model in parallel to the EKF and updating the evolution model with certain periodicity. The evolution model transition matrix is identified using the history of the estimated resistivity distribution obtained by a sensitivity matrix based algorithm and a Newton-Raphson algorithm. To numerically identify the linear evolution model, the Ibrahim time-domain method is used. The investigation is performed by numerical simulations of a domain with time-varying resistivity and by experimental data collected from the boundary of a human chest during normal breathing. The obtained dynamic resistivity values lie within the expected values for the tissues of a human chest. The EKF results suggest that the tracking ability is significantly improved with this approach.

Noise reduction in biomedical speech signal processing based on time and frequency Kalman filtering combined with spectral subtraction

The goal of this paper is to study and propose a new technique for noise reduction used during the reconstruction of speech signals, particularly for biomedical applications. The proposed method is based on Kalman filtering in the time domain combined with spectral subtraction. Comparison with discrete Kalman filter in the frequency domain shows better performance of the proposed technique. The performance is evaluated by using the segmental signal-to-noise ratio and the Itakura-Saito`s distance. Results have shown that Kalman`s filter in time combined with spectral subtraction is more robust and efficient, improving the Itakura-Saito`s distance by up to four times. (C) 2007 Elsevier Ltd. All rights reserved.

Wavelet correlation between subjects: A time-scale data driven analysis for brain mapping using fMRI

Functional magnetic resonance imaging (fMRI) based on BOLD signal has been used to indirectly measure the local neural activity induced by cognitive tasks or stimulation. Most fMRI data analysis is carried out using the general linear model (GLM), a statistical approach which predicts the changes in the observed BOLD response based on an expected hemodynamic response function (HRF). In cases when the task is cognitively complex or in cases of diseases, variations in shape and/or delay may reduce the reliability of results. A novel exploratory method using fMRI data, which attempts to discriminate between neurophysiological signals induced by the stimulation protocol from artifacts or other confounding factors, is introduced in this paper. This new method is based on the fusion between correlation analysis and the discrete wavelet transform, to identify similarities in the time course of the BOLD signal in a group of volunteers. We illustrate the usefulness of this approach by analyzing fMRI data from normal subjects presented with standardized human face pictures expressing different degrees of sadness. The results show that the proposed wavelet correlation analysis has greater statistical power than conventional GLM or time domain intersubject correlation analysis. (C) 2010 Elsevier B.V. All rights reserved.

Development of Hepatitis C Virus Genotyping by Real-Time PCR Based on the NS5B Region

Background: Hepatitis C virus (HCV) genotyping is the most significant predictor of the response to antiviral therapy. The aim of this study was to develop and evaluate a novel real-time PCR method for HCV genotyping based on the NS5B region. Methodology/Principal Findings: Two triplex reaction sets were designed, one to detect genotypes 1a, 1b and 3a; and another to detect genotypes 2a, 2b, and 2c. This approach had an overall sensitivity of 97.0%, detecting 295 of the 304 tested samples. All samples genotyped by real-time PCR had the same type that was assigned using LiPA version 1 (Line in Probe Assay). Although LiPA v. 1 was not able to subtype 68 of the 295 samples (23.0%) and rendered different subtype results from those assigned by real-time PCR for 12/295 samples (4.0%), NS5B sequencing and real-time PCR results agreed in all 146 tested cases. Analytical sensitivity of the real-time PCR assay was determined by end-point dilution of the 5000 IU/ml member of the OptiQuant HCV RNA panel. The lower limit of detection was estimated to be 125 IU/ml for genotype 3a, 250 IU/ml for genotypes 1b and 2b, and 500 IU/ml for genotype 1a. Conclusions/Significance: The total time required for performing this assay was two hours, compared to four hours required for LiPA v. 1 after PCR-amplification. Furthermore, the estimated reaction cost was nine times lower than that of available commercial methods in Brazil. Thus, we have developed an efficient, feasible, and affordable method for HCV genotype identification.

Performance of a multi-level TDR system exposed to tropical field conditions - A time-space comparison with tensiometry

Time-domain reflectometry (TDR) is an important technique to obtain series of soil water content measurements in the field. Diode-segmented probes represent an improvement in TDR applicability, allowing measurements of the soil water content profile with a single probe. In this paper we explore an extensive soil water content dataset obtained by tensiometry and TDR from internal drainage experiments in two consecutive years in a tropical soil in Brazil. Comparisons between the variation patterns of the water content estimated by both methods exhibited evidences of deterioration of the TDR system during this two year period at field conditions. The results showed consistency in the variation pattern for the tensiometry data, whereas TDR estimates were inconsistent, with sensitivity decreasing over time. This suggests that difficulties may arise for the long-term use of this TDR system under tropical field conditions. (c) 2008 Elsevier B.V. All rights reserved.

Coda wave attenuation in the Parecis Basin, Amazon Craton, Brazil: sensitivity to basement depth

Small local earthquakes from two aftershock sequences in Porto dos GaA(0)chos, Amazon craton-Brazil, were used to estimate the coda wave attenuation in the frequency band of 1 to 24 Hz. The time-domain coda-decay method of a single backscattering model is employed to estimate frequency dependence of the quality factor (Q (c)) of coda waves modeled usingwhere Q (0) is the coda quality factor at frequency of 1 Hz and eta is the frequency parameter. We also used the independent frequency model approach (Morozov, Geophys J Int, 175:239-252, 2008), based in the temporal attenuation coefficient, chi(f) instead of Q(f), given by the equation for the calculation of the geometrical attenuation (gamma) and effective attenuation Q (c) values have been computed at central frequencies (and band) of 1.5 (1-2), 3.0 (2-4), 6.0 (4-8), 9.0 (6-12), 12 (8-16), and 18 (12-24) Hz for five different datasets selected according to the geotectonic environment as well as the ability to sample shallow or deeper structures, particularly the sediments of the Parecis basin and the crystalline basement of the Amazon craton. For the Parecis basin for the surrounding shield and for the whole region of Porto dos GaA(0)chos Using the independent frequency model, we found: for the cratonic zone, gamma = 0.014 s (-aEuro parts per thousand 1), nu a parts per thousand 1.12; for the basin zone with sediments of similar to 500 m, gamma = 0.031 s (-aEuro parts per thousand 1), nu a parts per thousand 1.27; and for the Parecis basin with sediments of similar to 1,000 m, gamma = 0.047 s (-aEuro parts per thousand 1), nu a parts per thousand 1.42. Analysis of the attenuation factor (Q (c)) for different values of the geometrical spreading parameter (nu) indicated that an increase of nu generally causes an increase in Q (c), both in the basin as well as in the craton. But the differences in the attenuation between different geological environments are maintained for different models of geometrical spreading. It was shown that the energy of coda waves is attenuated more strongly in the sediments, (in the deepest part of the basin), than in the basement, (in the craton). Thus, the coda wave analysis can contribute to studies of geological structures in the upper crust, as the average coda quality factor is dependent on the thickness of sedimentary layer.

A discontinuous-Galerkin-based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user-defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous-Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements. boundary locking is avoided and optimal-order convergence is achieved. This is shown through numerical experiments in reaction-diffusion problems. Copyright (c) 2008 John Wiley & Sons, Ltd.

Fermion perturbations in string theory black holes

In this paper we study fermion perturbations in four-dimensional black holes of string theory, obtained either from a non-extreme configuration of three intersecting five-branes with a boost along the common string or from a non-extreme intersecting system of two two-branes and two five-branes. The Dirac equation for the massless neutrino field, after conformal re-scaling of the metric, is written as a wave equation suitable to study the time evolution of the perturbation. We perform a numerical integration of the evolution equation, and with the aid of Prony fitting of the time-domain profile, we calculate the complex frequencies that dominate the quasinormal ringing stage, and also determine these quantities by the semi-analytical sixth-order WKB method. We also find numerically the decay factor of fermion fields at very late times, and show that the falloff is identical to those showing for massless fields in other four-dimensional black hole spacetimes.

Estimativas do teor de umidade empregando o método GPR: uma avaliação comparativa em experimentos de laboratório e campo

Este trabalho faz uma análise das estimativas de teores de umidade obtidas com o método Ground Penetrating Radar (GPR) comparativamente às determinadas com os métodos Time Domain Reflectometry (TDR) e gravimétrico. Os dados foram obtidos em dois experimentos diferentes: um experimento controlado em laboratório buscando reproduzir um meio homogêneo onde foram obtidas as medidas de umidade com GPR (antenas de 1 GHz) e TDR, e outro experimento de campo onde foram obtidos dados com GPR (antenas de 200 MHz) e de amostras de solos do local. Para a obtenção das estimativas a partir do método GPR foram analisados os eventos relativos à onda de transmissão direta entre as antenas, onda refratada criticamente e onda refletida em interfaces com diferentes propriedades elétricas.O GPR mostrou-se sensível às variações de umidades presentes nos dois experimentos e apresentou boa correlação com os dados obtidos com TDR (REQM de0,007 m³m-3) e das amostras (REQM de 0,039 m³m-3).

Quasinormal modes of black holes in anti-de Sitter space: A numerical study of the eikonal limit

Using series solutions and time-domain evolutions, we probe the eikonal limit of the gravitational and scalar-field quasinormal modes of large black holes and black branes in anti-de Sitter backgrounds. These results are particularly relevant for the AdS/CFT correspondence, since the eikonal regime is characterized by the existence of long-lived modes which (presumably) dominate the decay time scale of the perturbations. We confirm all the main qualitative features of these slowly damped modes as predicted by Festuccia and Liu [G. Festuccia and H. Liu, arXiv:0811.1033.] for the scalar-field (tensor-type gravitational) fluctuations. However, quantitatively we find dimensional-dependent correction factors. We also investigate the dependence of the quasinormal mode frequencies on the horizon radius of the black hole (brane) and the angular momentum (wave number) of vector- and scalar-type gravitational perturbations.

Relationship of autonomic imbalance and circadian disruption with obesity and type 2 diabetes in resistant hypertensive patients

Background: Hypertension, diabetes and obesity are not isolated findings, but a series of interacting interactive physiologic derangements. Taking into account genetic background and lifestyle behavior, AI (autonomic imbalance) could be a common root for RHTN (resistant hypertension) or RHTN plus type 2 diabetes (T2D) comorbidity development. Moreover, circadian disruption can lead to metabolic and vasomotor impairments such as obesity, insulin resistance and resistant hypertension. In order to better understand the triggered emergence of obesity and T2D comorbidity in resistant hypertension, we investigated the pattern of autonomic activity in the circadian rhythm in RHTN with and without type 2 diabetes (T2D), and its relationship with serum adiponectin concentration. Methods: Twenty five RHTN patients (15 non-T2D and 10 T2D, 15 males, 10 females; age range 34 to 70 years) were evaluated using the following parameters: BMI (body mass index), biochemical analysis, serum adiponectinemia, echocardiogram and ambulatory electrocardiograph heart rate variability (HRV) in time and frequency domains stratified into three periods: 24 hour, day time and night time. Results: Both groups demonstrated similar characteristics despite of the laboratory analysis concerning T2D like fasting glucose, HbA1c levels and hypertriglyceridemia. Both groups also revealed disruption of the circadian rhythm: inverted sympathetic and parasympathetic tones during day (parasympathetic > sympathetic tone) and night periods (sympathetic > parasympathetic tone). T2D group had increased BMI and serum triglyceride levels (mean 33.7 +/- 4.0 vs 26.6 +/- 3.7 kg/m(2) - p = 0.00; 254.8 +/- 226.4 vs 108.6 +/- 48.7 mg/dL - p = 0.04), lower levels of adiponectin (6729.7 +/- 3381.5 vs 10911.5 +/- 5554.0 ng/mL - p = 0.04) and greater autonomic imbalance evaluated by HRV parameters in time domain compared to non-T2D RHTN patients. Total patients had HRV correlated positively with serum adiponectin (r = 0.37 [95% CI - 0.04 - 1.00] p = 0.03), negatively with HbA1c levels (r = -0.58 [95% CI -1.00 - -0.3] p = 0.00) and also adiponectin correlated negatively with HbA1c levels (r = -0.40 [95% CI -1.00 - -0.07] p = 0.02). Conclusion: Type 2 diabetes comorbidity is associated with greater autonomic imbalance, lower adiponectin levels and greater BMI in RHTN patients. Similar circadian disruption was also found in both groups indicating the importance of lifestyle behavior in the genesis of RHTN.

Sudan Black B treatment reduces autofluorescence and improves resolution of in situ hybridization specific fluorescent signals of brain sections

Interference by autofluorescence is one of the major concerns of immunofluorescence analysis of in situ hybridization-based diagnostic assays. We present a useful technique that reduces autofluorescent background without affecting the tissue integrity or direct immunofluorescence signals in brain sections. Using six different protocols, such as ammonia/ethanol, Sudan Black B (SBB) in 70% ethanol, photobleaching with UV light and different combinations of them in both formalin-fixed paraffin-embedded and frozen human brain tissue sections, we have found that tissue treatment of SBB in a concentration of 0.1% in 70% ethanol is the best approach to reduce/eliminate tissue autofluorescence and background, while preserving the specific fluorescence hybridization signals. This strategy is a feasible, non-time consuming method that provides a reasonable compromise between total reduction of the tissue autofluorescence and maintenance of specific fluorescent labels.