Termination of Atrial Fibrillation by Catheter Ablation can be Successfully Predicted from Baseline ECG

Purpose Multiple organization indices (OIs) have been used to predict the outcome of stepwise catheter ablation (step-CA) in long-standing persistent AF (pers-AF), however with limited success. Our study aims at developing innovative OIs from baseline (BL, before ablation) ECG in order to predict the outcome of step-CA. Methods Fourteen consecutive patients (pts) (60±5 y, AF duration 21±9 m) underwent a step-CA consisting in pulmonary veins isolation, left atrial (LA) defragmentation and linear ablations, and right atrial (RA) ablations if non terminated. Chest lead V6 was placed in the back (V6b). After QRST cancellation from chest leads V1 to V6b, two OIs were computed to quantify the harmonic components of ECG atrial activity: 1) phase difference variance (PD) between the AF harmonic components as a measure of AF regularity, and 2) adaptive OI (AOI) evaluating the time evolution of the AF harmonic components. Both indices were compared to classical ones: a spectrum-based OI (SOI) and ECG AF cycle length (AFCL). Results Pers-AF was terminated into sinus rhythm or atrial tachycardia in 10/14 pts during step-CA, 8 during LA (LT), 2 during RA (RT) ablation, and 4 were non terminated (NT). The figure shows that LT was best separated from RT/NT before ablation by AOI computed on lead V1 (A) and PD from lead V6b (B) as compared to SOI and AFCL respectively. Conclusion Our results suggest that adaptive OIs computed before ablation perform better than classical OIs for separating LT from RT/NT pts. These findings are indicative of a higher baseline organization in LT pts that could be used to select candidates for the restoration of sinus rhythm by step-CA.

Monte Carlo simulation of interface alloying

An Ising-like model, with interactions ranging up to next-nearest-neighbor pairs, is used to simulate the process of interface alloying. Interactions are chosen to stabilize an intermediate "antiferromagnetic" ordered structure. The dynamics proceeds exclusively by atom-vacancy exchanges. In order to characterize the process, the time evolution of the width of the intermediate ordered region and the diffusion length is studied. Both lengths are found to follow a power-law evolution with exponents depending on the characteristic features of the model.

Monte Carlo Study of the Growth of L12 ordered domains in fcc A3B binary alloys

A Monte Carlo study of the late time growth of L12-ordered domains in a fcc A3B binary alloy is presented. The energy of the alloy has been modeled by a nearest-neighbor interaction Ising Hamiltonian. The system exhibits a fourfold degenerated ground state and two kinds of interfaces separating ordered domains: flat and curved antiphase boundaries. Two different dynamics are used in the simulations: the standard atom-atom exchange mechanism and the more realistic vacancy-atom exchange mechanism. The results obtained by both methods are compared. In particular we study the time evolution of the excess energy, the structure factor and the mean distance between walls. In the case of atom-atom exchange mechanism anisotropic growth has been found: two characteristic lengths are needed in order to describe the evolution. Contrarily, with the vacancyatom exchange mechanism scaling with a single length holds. Results are contrasted with existing experiments in Cu3Au and theories for anisotropic growth.

Decay of unstable states in the presence of colored noise and random initial conditions. II. Analog experiments and digital simulations

The decay of an unstable state under the influence of external colored noise has been studied by means of analog experiments and digital simulations. For both fixed and random initial conditions, the time evolution of the second moment ¿x2(t)¿ of the system variable was determined and then used to evaluate the nonlinear relaxation time. The results obtained are found to be in excellent agreement with the theoretical predictions of the immediately preceding paper [Casademunt, Jiménez-Aquino, and Sancho, Phys. Rev. A 40, 5905 (1989)].

Effects of small surface tension in Hele-Shaw multifinger dynamics: an analytical and numerical study

We study the singular effects of vanishingly small surface tension on the dynamics of finger competition in the Saffman-Taylor problem, using the asymptotic techniques described by Tanveer [Philos. Trans. R. Soc. London, Ser. A 343, 155 (1993)] and Siegel and Tanveer [Phys. Rev. Lett. 76, 419 (1996)], as well as direct numerical computation, following the numerical scheme of Hou, Lowengrub, and Shelley [J. Comput. Phys. 114, 312 (1994)]. We demonstrate the dramatic effects of small surface tension on the late time evolution of two-finger configurations with respect to exact (nonsingular) zero-surface-tension solutions. The effect is present even when the relevant zero-surface-tension solution has asymptotic behavior consistent with selection theory. Such singular effects, therefore, cannot be traced back to steady state selection theory, and imply a drastic global change in the structure of phase-space flow. They can be interpreted in the framework of a recently introduced dynamical solvability scenario according to which surface tension unfolds the structurally unstable flow, restoring the hyperbolicity of multifinger fixed points.

Competitive evaporation in arrays of droplets

We consider the evaporation of periodic arrays of initially equal droplets in two-dimensional systems with open (absorbing) boundaries. Our study is based on the numerical solution of the Cahn-Hilliard equation. We show that due to cooperative effects the droplets which are further from the boundary may evaporate earlier than those in the boundary¿s vicinity. The time evolution of the overall amount of matter in the system is also studied.

Particle agglomeration study in rf silane plasmas: In situ study by polarization-sensitive laser light scattering

To determine self‐consistently the time evolution of particle size and their number density in situ multi‐angle polarization‐sensitive laser light scattering was used. Cross‐polarization intensities (incident and scattered light intensities with opposite polarization) measured at 135° and ex situ transmission electronic microscopy analysis demonstrate the existence of nonspherical agglomerates during the early phase of agglomeration. Later in the particle time development both techniques reveal spherical particles again. The presence of strong cross‐polarization intensities is accompanied by low‐frequency instabilities detected on the scattered light intensities and plasma emission. It is found that the particle radius and particle number density during the agglomeration phase can be well described by the Brownian free molecule coagulation model. Application of this neutral particle coagulation model is justified by calculation of the particle charge whereby it is shown that particles of a few tens of nanometer can be considered as neutral under our experimental conditions. The measured particle dispersion can be well described by a Brownian free molecule coagulation model including a log‐normal particle size distribution.

Glutamine synthesis rate in the hyperammonaemic rat brain using simultaneous localized in vivo H-1 and N-15 MRS

Objectives: Glutamine synthetase is a critical step in the glutamate-glutamine cycle, the major mechanism of glutamate neurotransmission and is implicated in the mechanism of ammonia toxicity. 15N MRS is an alternative approach to 13C MRS in studying glutamate- glutamine metabolism. 15N MRS studies allow to measure an apparent glutamine synthesis rate (Vsyn) which reflects a combination of the glutamate- glutamine cycle activity (Vnt) and net glutamine accumulation. The net glutamine synthesis (Vsyn-Vnt) can be directly measured from 1H NMR. Therefore, the aim of this study was to perform in vivo localized 1H MRS interleaved with 15N MRS to directly measure the net glutamine synthesis rate and the apparent glutamine synthesis rate under 15N labeled ammonia infusion in the rat brain, respectively. Methods: 1H and 15N MRS data were acquired interleaved on a 9.4T system (Varian/Magnex Scientific) using 5 rats. 15NH4Cl solution was infused continuously into the femoral vein for up to 10 h (4.5 mmol/h/kg).1 The plasma ammonia concentration was increased to 0.95±0.08 mmol/L (Analox GM7 analyzer). 1H spectra were acquired and quantified as described previously.2 15N unlocalized and localized spectra were acquired using the sequence;3 and quantified using AMARES and an external reference method.4 The metabolic model used to analyze the total Gln and 5-15N labeled Gln time courses is shown on Figure 1A. Results: Glutamine concentration increased from 2.5±0.3 to 15±3.3 mmol/kg whereas the total glutamate concentrations remained unchanged (Figure 1B). The linear fit of the time-evolution of the total Gln from the 1H spectra gave the net synthesis flux (Vsyn-Vnt), which was 0.021± 0.006 mmol/min per g (Figure 1D). The 5-15N Gln peak (_271 ppm) was visible in the first and all subsequent scans, whereas the 2-15N Gln/Glu peak (_342 ppm) appeared after B1.5 h (Figure 1C). From the in vivo 5-15N Gln time course, Vsyn = 0.29±0.1 mmol/min per g and a plasma NH3 fractional enrichment of 71%±6% were calculated. Vnt was 0.26±0.1 mmol/min/g, obtained assuming a negligible Gln efflux.5 Vsyn and Vnt were within the range of 13C NMR measurements.6 Conclusion: The combination of 1H and 15N NMR allowed for the first time a direct and localized measurement of Vnt and apparent glutamine synthesis rate. Vnt is approximately one order of magnitude faster than the net glutamine accumulation.

Competitive evaporation in arrays of droplets

We consider the evaporation of periodic arrays of initially equal droplets in two-dimensional systems with open (absorbing) boundaries. Our study is based on the numerical solution of the Cahn-Hilliard equation. We show that due to cooperative effects the droplets which are further from the boundary may evaporate earlier than those in the boundary¿s vicinity. The time evolution of the overall amount of matter in the system is also studied.

Formation and distribution of point defects on a disclination line near a free nematic interface

Point defects of opposite signs can alternately nucleate on the -1/2 disclination line that forms near the free surface of a confined nematic liquid crystal. We show the existence of metastable configurations consisting of periodic repetitions of such defects. These configurations are characterized by a minimal interdefect spacing that is seen to depend on sample thickness and on an applied electric field. The time evolution of the defect distribution suggests that the defects attract at small distances and repel at large distances.

Measurements of glial metabolic fluxes with C-11-acetate using positron emission and an adapted NMR-based metabolic modeling approach

Objectives: Acetate brain metabolism has the particularity to occur specifically in glial cells. Labeling studies, using acetate labeled either with 13C (NMR) or 11C (PET), are governed by the same biochemical reactions and thus follow the same mathematical principles. In this study, the objective was to adapt an NMR acetate brain metabolism model to analyse [1-11C]acetate infusion in rats. Methods: Brain acetate infusion experiments were modeled using a two-compartment model approach used in NMR.1-3 The [1-11C]acetate labeling study was done using a beta scintillator.4 The measured radioactive signal represents the time evolution of the sum of all labeled metabolites in the brain. Using a coincidence counter in parallel, an arterial input curve was measured. The 11C at position C-1 of acetate is metabolized in the first turn of the TCA cycle to the position 5 of glutamate (Figure 1A). Through the neurotransmission process, it is further transported to the position 5 of glutamine and the position 5 of neuronal glutamate. After the second turn of the TCA cycle, tracer from [1-11C]acetate (and also a part from glial [5-11C]glutamate) is transferred to glial [1-11C]glutamate and further to [1-11C]glutamine and neuronal glutamate through the neurotransmission cycle. Brain poster session: oxidative mechanisms S460 Journal of Cerebral Blood Flow & Metabolism (2009) 29, S455-S466 Results: The standard acetate two-pool PET model describes the system by a plasma pool and a tissue pool linked by rate constants. Experimental data are not fully described with only one tissue compartment (Figure 1B). The modified NMR model was fitted successfully to tissue time-activity curves from 6 single animals, by varying the glial mitochondrial fluxes and the neurotransmission flux Vnt. A glial composite rate constant Kgtg=Vgtg/[Ace]plasma was extracted. Considering an average acetate concentration in plasma of 1 mmol/g5 and the negligible additional amount injected, we found an average Vgtg = 0.08±0.02 (n = 6), in agreement with previous NMR measurements.1 The tissue time-activity curve is dominated by glial glutamate and later by glutamine (Figure 1B). Labeling of neuronal pools has a low influence, at least for the 20 mins of beta-probe acquisition. Based on the high diffusivity of CO2 across the blood-brain barrier; 11CO2 is not predominant in the total tissue curve, even if the brain CO2 pool is big compared with other metabolites, due to its strong dilution through unlabeled CO2 from neuronal metabolism and diffusion from plasma. Conclusion: The two-compartment model presented here is also able to fit data of positron emission experiments and to extract specific glial metabolic fluxes. 11C-labeled acetate presents an alternative for faster measurements of glial oxidative metabolism compared to NMR, potentially applicable to human PET imaging. However, to quantify the relative value of the TCA cycle flux compared to the transmitochondrial flux, the chemical sensitivity of NMR is required. PET and NMR are thus complementary.

Aggregation of liposomes induced by calcium: A structural and kinetic study

In this work, the calcium-induced aggregation of phosphatidylserine liposomes is probed by means of the analysis of the kinetics of such process as well as the aggregate morphology. This novel characterization of liposome aggregation involves the use of static and dynamic light-scattering techniques to obtain kinetic exponents and fractal dimensions. For salt concentrations larger than 5 mM, a diffusion-limited aggregation regime is observed and the Brownian kernel properly describes the time evolution of the diffusion coefficient. For slow kinetics, a slightly modified multiple contact kernel is required. In any case, a time evolution model based on the numerical resolution of Smoluchowski's equation is proposed in order to establish a theoretical description for the aggregating system. Such a model provides an alternative procedure to determine the dimerization constant, which might supply valuable information about interaction mechanisms between phospholipid vesicles.

Optimization of Snow Drifting Mitigation and Control Methods for Iowa Conditions Part II, TR-626, 2015

Blowing and drifting of snow is a major concern for transportation efficiency and road safety in regions where their development is common. One common way to mitigate snow drift on roadways is to install plastic snow fences. Correct design of snow fences is critical for road safety and maintaining the roads open during winter in the US Midwest and other states affected by large snow events during the winter season and to maintain costs related to accumulation of snow on the roads and repair of roads to minimum levels. Of critical importance for road safety is the protection against snow drifting in regions with narrow rights of way, where standard fences cannot be deployed at the recommended distance from the road. Designing snow fences requires sound engineering judgment and a thorough evaluation of the potential for snow blowing and drifting at the construction site. The evaluation includes site-specific design parameters typically obtained with semi-empirical relations characterizing the local transport conditions. Among the critical parameters involved in fence design and assessment of their post-construction efficiency is the quantification of the snow accumulation at fence sites. The present study proposes a joint experimental and numerical approach to monitor snow deposits around snow fences, quantitatively estimate snow deposits in the field, asses the efficiency and improve the design of snow fences. Snow deposit profiles were mapped using GPS based real-time kinematic surveys (RTK) conducted at the monitored field site during and after snow storms. The monitored site allowed testing different snow fence designs under close to identical conditions over four winter seasons. The study also discusses the detailed monitoring system and analysis of weather forecast and meteorological conditions at the monitored sites. A main goal of the present study was to assess the performance of lightweight plastic snow fences with a lower porosity than the typical 50% porosity used in standard designs of such fences. The field data collected during the first winter was used to identify the best design for snow fences with a porosity of 50%. Flow fields obtained from numerical simulations showed that the fence design that worked the best during the first winter induced the formation of an elongated area of small velocity magnitude close to the ground. This information was used to identify other candidates for optimum design of fences with a lower porosity. Two of the designs with a fence porosity of 30% that were found to perform well based on results of numerical simulations were tested in the field during the second winter along with the best performing design for fences with a porosity of 50%. Field data showed that the length of the snow deposit away from the fence was reduced by about 30% for the two proposed lower-porosity (30%) fence designs compared to the best design identified for fences with a porosity of 50%. Moreover, one of the lower-porosity designs tested in the field showed no significant snow deposition within the bottom gap region beneath the fence. Thus, a major outcome of this study is to recommend using plastic snow fences with a porosity of 30%. It is expected that this lower-porosity design will continue to work well for even more severe snow events or for successive snow events occurring during the same winter. The approach advocated in the present study allowed making general recommendations for optimizing the design of lower-porosity plastic snow fences. This approach can be extended to improve the design of other types of snow fences. Some preliminary work for living snow fences is also discussed. Another major contribution of this study is to propose, develop protocols and test a novel technique based on close range photogrammetry (CRP) to quantify the snow deposits trapped snow fences. As image data can be acquired continuously, the time evolution of the volume of snow retained by a snow fence during a storm or during a whole winter season can, in principle, be obtained. Moreover, CRP is a non-intrusive method that eliminates the need to perform man-made measurements during the storms, which are difficult and sometimes dangerous to perform. Presently, there is lots of empiricism in the design of snow fences due to lack of data on fence storage capacity on how snow deposits change with the fence design and snow storm characteristics and in the estimation of the main parameters used by the state DOTs to design snow fences at a given site. The availability of such information from CRP measurements should provide critical data for the evaluation of the performance of a certain snow fence design that is tested by the IDOT. As part of the present study, the novel CRP method is tested at several sites. The present study also discusses some attempts and preliminary work to determine the snow relocation coefficient which is one of the main variables that has to be estimated by IDOT engineers when using the standard snow fence design software (Snow Drift Profiler, Tabler, 2006). Our analysis showed that standard empirical formulas did not produce reasonable values when applied at the Iowa test sites monitored as part of the present study and that simple methods to estimate this variable are not reliable. The present study makes recommendations for the development of a new methodology based on Large Scale Particle Image Velocimetry that can directly measure the snow drift fluxes and the amount of snow relocated by the fence.

Toxicokinetic modeling of folpet fungicide and its ring-biomarkers of exposure in humans.

A human in vivo toxicokinetic model was built to allow a better understanding of the toxicokinetics of folpet fungicide and its key ring biomarkers of exposure: phthalimide (PI), phthalamic acid (PAA) and phthalic acid (PA). Both PI and the sum of ring metabolites, expressed as PA equivalents (PAeq), may be used as biomarkers of exposure. The conceptual representation of the model was based on the analysis of the time course of these biomarkers in volunteers orally and dermally exposed to folpet. In the model, compartments were also used to represent the body burden of folpet and experimentally relevant PI, PAA and PA ring metabolites in blood and in key tissues as well as in excreta, hence urinary and feces. The time evolution of these biomarkers in each compartment of the model was then mathematically described by a system of coupled differential equations. The mathematical parameters of the model were then determined from best fits to the time courses of PI and PAeq in blood and urine of five volunteers administered orally 1 mg kg(-1) and dermally 10 mg kg(-1) of folpet. In the case of oral administration, the mean elimination half-life of PI from blood (through feces, urine or metabolism) was found to be 39.9 h as compared with 28.0 h for PAeq. In the case of a dermal application, mean elimination half-life of PI and PAeq was estimated to be 34.3 and 29.3 h, respectively. The average final fractions of administered dose recovered in urine as PI over the 0-96 h period were 0.030 and 0.002%, for oral and dermal exposure, respectively. Corresponding values for PAeq were 24.5 and 1.83%, respectively. Finally, the average clearance rate of PI from blood calculated from the oral and dermal data was 0.09 ± 0.03 and 0.13 ± 0.05 ml h(-1) while the volume of distribution was 4.30 ± 1.12 and 6.05 ± 2.22 l, respectively. It was not possible to obtain the corresponding values from PAeq data owing to the lack of blood time course data.

Modeling morphogen gradient formation from arbitrary realistically shaped sources.

Much of the analytical modeling of morphogen profiles is based on simplistic scenarios, where the source is abstracted to be point-like and fixed in time, and where only the steady state solution of the morphogen gradient in one dimension is considered. Here we develop a general formalism allowing to model diffusive gradient formation from an arbitrary source. This mathematical framework, based on the Green's function method, applies to various diffusion problems. In this paper, we illustrate our theory with the explicit example of the Bicoid gradient establishment in Drosophila embryos. The gradient formation arises by protein translation from a mRNA distribution followed by morphogen diffusion with linear degradation. We investigate quantitatively the influence of spatial extension and time evolution of the source on the morphogen profile. For different biologically meaningful cases, we obtain explicit analytical expressions for both the steady state and time-dependent 1D problems. We show that extended sources, whether of finite size or normally distributed, give rise to more realistic gradients compared to a single point-source at the origin. Furthermore, the steady state solutions are fully compatible with a decreasing exponential behavior of the profile. We also consider the case of a dynamic source (e.g. bicoid mRNA diffusion) for which a protein profile similar to the ones obtained from static sources can be achieved.