Real-time evaluation of diffusion of the local anesthetic solution during peribulbar block using ultrasound imaging and clinical correlates of diffusion

The aims of this prospective observational study were to assess the incidence of intraconal spread during peribulbar (extraconal) anesthesia by real-time ultrasound imaging of the retro-orbital compartment and to determine whether a complete sensory and motor block (with akinesia) of the eye is directly related to the intraconal spread.

Real-time simulation of nonequilibrium transport of magnetization in large open quantum spin systems driven by dissipation

Using quantum Monte Carlo, we study the nonequilibrium transport of magnetization in large open strongly correlated quantum spin-12 systems driven by purely dissipative processes that conserve the uniform or staggered magnetization, disregarding unitary Hamiltonian dynamics. We prepare both a low-temperature Heisenberg ferromagnet and an antiferromagnet in two parts of the system that are initially isolated from each other. We then bring the two subsystems in contact and study their real-time dissipative dynamics for different geometries. The flow of the uniform or staggered magnetization from one part of the system to the other is described by a diffusion equation that can be derived analytically.

High-field 1H T1 and T2 NMR Relaxation Time Measurements of H2O in Homeopathic Preparations of Quartz, Sulfur and Copper Sulfate

Quantitative meta-analyses of randomized clinical trials investigating the specific therapeutic efficacy of homeopathic remedies yielded statistically significant differences compared to placebo. Since the remedies used contained mostly only very low concentrations of pharmacologically active compounds, these effects cannot be accounted for within the framework of current pharmacology. Theories to explain clinical effects of homeopathic remedies are partially based upon changes in diluent structure. To investigate the latter, we measured for the first time high-field (600/500 MHz) 1H T1 and T2 nuclear magnetic resonance relaxation times of H2O in homeopathic preparations with concurrent contamination control by inductively coupled plasma mass spectrometry (ICP-MS). Homeopathic preparations of quartz (10c–30c, n = 21, corresponding to iterative dilutions of 100−10–100−30), sulfur (13x–30x, n = 18, 10−13–10−30), and copper sulfate (11c–30c, n = 20, 100−11–100−30) were compared to n = 10 independent controls each (analogously agitated dilution medium) in randomized and blinded experiments. In none of the samples, the concentration of any element analyzed by ICP-MS exceeded 10 ppb. In the first measurement series (600 MHz), there was a significant increase in T1 for all samples as a function of time, and there were no significant differences between homeopathic potencies and controls. In the second measurement series (500 MHz) 1 year after preparation, we observed statistically significant increased T1 relaxation times for homeopathic sulfur preparations compared to controls. Fifteen out of 18 correlations between sample triplicates were higher for controls than for homeopathic preparations. No conclusive explanation for these phenomena can be given at present. Possible hypotheses involve differential leaching from the measurement vessel walls or a change in water molecule dynamics, i.e., in rotational correlation time and/or diffusion. Homeopathic preparations thus may exhibit specific physicochemical properties that need to be determined in detail in future investigations.

A numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.

Geochemical synthesis for the Effingen Member in boreholes at Oftringen, Gösgen and Küttigen

The Effingen Member is a low-permeability rock unit of Oxfordian age (ca. 160 Ma) that occurs across northern Switzerland. It comprises sandy calcareous marls and (argillaceous) limestones. This report describes the hydrogeochemistry, mineralogy and supporting physical properties of the Effingen Member in three boreholes in the Jura-Südfuss area: Oftringen, Gösgen and Küttigen, where it is 220–240 m thick. The top of the Effingen Member is at 420, 66 and 32 m depths at the three sites. Core materials are available from Oftringen and Gösgen, whereas information from Küttigen is limited to cuttings, in-situ hydrogeological testing and geophysical logging. Hydrogeological boundaries of the Effingen Member vary between locations. Ground-water flows were identified during drilling at the top (Geissberg Member), but not at the base, of the Effingen Member at Oftringen, at the base (Hauptrogenstein Formation) of the Effingen Member at Gösgen, and in a limestone layer (Gerstenhübel unit) within the Effingen Member at Küttigen. The marls and limestones of the Effingen Member have carbonate contents of 46–91 wt.-% and clay-mineral contents of 5–37 wt.-%. Pyrite contents are up to 1.6 wt.-%, but no sulphate minerals were detected by routine analyses. Clay minerals are predominantly mixed-layer illite-smectite, illite and kaolinite, with sporadic traces of chlorite and smectite. Veins filled with calcite ± celestite occur through the Effingen Member at Oftringen but not at Gösgen or Küttigen. They formed at 50–70 ºC from externally derived fluids, probably of Miocene age. Water contents are 0.7–4.2 wt.-%, corresponding to a water-loss porosity range of 1.9–10.8 vol.-%. Specific surface areas, measured by the BET method, are 2–30 m2/g, correlating with clay-mineral contents. Water activity has been measured and yielded surprisingly low values down to 0.8. These cannot be explained by pore-water salinity alone and include other effects, such as changes in the fabric due to stress release or partial saturation. Observed variations in measurements are not fully understood. Cation exchange capacity (CEC) and exchangeable cation populations have been studied by the Ni-en method. CEC, derived from the consumption of the index cation Ni, is 9–99 meq/kgrock at a solid:liquid ratio of 1, correlating with the clay-mineral content. Cation concentrations in Ni-en extract solutions are in the order Na+≥Ca2+>Mg2+>K+>Sr2+. However, the analytical results from the Ni-en extractions have additional contributions from cations originating from pore water and from mineral dissolution reactions that occurred during extraction, and it was not possible to reliably quantify these contributions. Therefore, in-situ cation populations and selectivity coefficients could not be derived. A suite of methods have been used for characterising the chemical compositions of pore waters in the Effingen Member. Advective displacement was used on one sample from each Oftringen and Gösgen and is the only method that produces results that approach complete hydrochemical compositions. Aqueous extraction was used on core samples from these two boreholes and gives data only for Cl- and, in some cases, Br-. Out-diffusion was used on core samples from Oftringen and similarly gives data for Cl- and Br- only. For both aqueous extraction and out-diffusion, reaction of the experimental water with rock affected concentrations of cations, SO42 and alkalinity in experimental solutions. Another method, centrifugation, failed to extract pore water. Stable isotope ratios (δ18O and δ2H) of pore waters in core samples from Oftringen were analysed by the diffusive exchange method and helium contents of pore water in Oftringen samples were extracted for mass spectrometric analysis by quantitative outgassing of preserved core samples. Several lines of evidence indicate that drillcore samples might not have been fully saturated when opened and subsampled in the laboratory. These include comparisons of water-loss porosities with physical porosities, water-activity measurements, and high contents of dissolved gas as inferred from ground-water samples. There is no clear proof of partial saturation and it is unclear whether this might represent in-situ conditions or is due to exsolution of gas due to the pressure release since drilling. Partial saturation would have no impact on the recalculation of pore-water compositions from aqueous extraction experiments using water-loss porosity data. The largest uncertainty in the pore-water Cl- concentrations recalculated from aqueous extraction and out-diffusion experiments is the magnitude of the anion-accessible fraction of water-loss porosity. General experience of clay-mineral rich formations suggests that the anion-accessible porosity fraction is very often about 0.5 and generally in a range of 0.3 to 0.6 and tends to be inversely correlated with clay-mineral contents. Comparisons of the Cl- concentration in pore water obtained by advective displacement with that recalculated from aqueous extraction of an adjacent core sample suggests a fraction of 0.27 for an Oftringen sample, whereas the same procedure for a Gösgen sample suggests a value of 0.64. The former value for anion-accessible porosity fraction is presumed to be unrepresentative given the local mineralogical heterogeneity at that depth. Through-diffusion experiments with HTO and 36Cl- suggest that the anion-accessible porosity fraction in the Effingen Member at Oftringen and Gösgen is around 0.5. This value is proposed as a typical average for rocks of the Effingen Member, bearing in mind that it varies on a local scale in response to the heterogeneity of lithology and pore-space architecture. The substantial uncertainties associated with the approaches to estimating anion-accessible porosity propagate into the calculated values of in-situ pore-water Cl- concentrations. On the basis of aqueous extraction experiments, and using an anion-accessible porosity fraction of 0.5, Cl- concentrations in the Effingen Member at Oftringen reach a maximum of about 14 g/L in the centre. Cl- decreases upwards and downwards from that, forming a curved depth profile. Cl- contents in the Effingen Member at Gösgen increase with depth from about 3.5 g/L to about 14 g/L at the base of the cored profile (which corresponds to the centre of the formation). Out-diffusion experiments were carried out on four samples from Oftringen, distributed through the Effingen Member. Recalculated Cl- concentrations are similar to those from aqueous extraction for 3 out of the 4 samples, and somewhat lower for one sample. Concentrations of other components, i.e. Na+, K+, Ca2+, Mg2+, Sr2+, SO42- and HCO3- cannot be obtained from the aqueous extraction and out-diffusion experimental data because of mineral dissolution and cation exchange reactions during the experiments. Pore-water pH also is not constrained by those extraction experiments. The only experimental approach to obtain complete pore-water compositions for samples from Oftringen and Gösgen is advective displacement of pore water. The sample from Oftringen used for this experiment is from 445 m depth in the upper part of the Effingen Member and gave eluate with 16.5 g/L Cl- whereas aqueous extraction from a nearby sample indicated about 9 g/L Cl-. The sample from Gösgen used for advective displacement is from 123 m depth in the centre of the Effingen Member sequence and gave eluate with about 9 g/L Cl- whereas aqueous extraction gave 11.5 g/L Cl-. In both cases the pore waters have Na-(Ca)-Cl compositions and SO42- concentrations of about 1.1 g/L. The Gösgen sample has a Br/Cl ratio similar to that of sea water, whereas this ratio is lower for the Oftringen sample. Taking account of uncertainties in the applied experimental approaches, it is reasonable to place an upper limit of ca. 20 g/L on Cl- concentration for pore water in the Effingen Member in this area. There are major discrepancies between pore-water SO42- concentrations inferred from aqueous extraction or out-diffusion experiments and those obtained from advective displacement in both the Oftringen and Gösgen cases. A general conclusion is that all or at least part of the discrepancies are attributable to perturbation of the sulphur system and enhancement of SO42- by sulphate mineral dissolution and possibly minor pyrite oxidation during aqueous extraction and out-diffusion. Therefore, data for SO42- calculated from those pore-water sampling methods are considered not to be representative of in-situ conditions. A reference pore-water composition was defined for the Effingen Member in the Jura Südfuss area. It represents the probable upper limits of Cl- contents and corresponding anion and cation concentrations that are reasonably constrained by experimental data. Except for Cl- and possibly Na+ concentrations, this composition is poorly constrained especially with respect to SO42- and Ca2+ concentrations, and pH and alkalinity. Stable isotope compositions, δ18O and δ2H, of pore waters in the Effingen Member at Oftringen plot to the right of the meteoric water line, suggesting that 18O has been enriched by water-rock exchange, which indicates that the pore waters have a long residence time. A long residence time of pore water is supported by the level of dissolved 4He that has accumulated in pore water of the Effingen Member at Oftringen. This is comparable with, or slightly higher than, the amounts of 4He in the Opalinus Clay at Benken. Ground waters were sampled from flowing zones intersected by boreholes at the three locations. The general interpretation is that pore waters and ground-water solutes may have similar origins in Mesozoic and Cenozoic brackish-marine formations waters, but ground-water solutes have been diluted rather more than pore waters by ingress of Tertiary and Quaternary meteoric waters. The available hydrochemical data for pore waters from the Effingen Member at these three locations in the Jura-Südfuss area suggest that the geochemical system evolved slowly over geological periods of time, in which diffusion was an important mechanism of solute transport. The irregularity of Cl- and δ18O profiles and spatial variability of advective ground-water flows in the Malm-Dogger system suggests that palaeohydrogeological and hydrochemical responses to changing tectonic and surface environmental conditions were complex.

Linking brain connectivity across different time scales with electroencephalogram, functional magnetic resonance imaging, and diffusion tensor imaging

Structural and functional connectivity are intrinsic properties of the human brain and represent the amount of cognitive capacities of individual subjects. These connections are modulated due to development, learning, and disease. Momentary adaptations in functional connectivity alter the structural connections, which in turn affect the functional connectivity. Thus, structural and functional connectivity interact on a broad timescale. In this study, we aimed to explore distinct measures of connectivity assessed by functional magnetic resonance imaging and diffusion tensor imaging and their association to the dominant electroencephalogram oscillatory property at rest: the individual alpha frequency (IAF). We found that in 21 healthy young subjects, small intraindividual temporal IAF fluctuations were correlated to increased blood oxygenation level-dependent signal in brain areas associated to working memory functions and to the modulation of attention. These areas colocalized with functionally connected networks supporting the respective functions. Furthermore, subjects with higher IAF show increased fractional anisotropy values in fascicles connecting the above-mentioned areas and networks. Hence, due to a multimodal approach a consistent functionally and structurally connected network related to IAF was observed.

Spinal Cord Diffusion-Tensor Imaging and Motor-evoked Potentials in Multiple Sclerosis Patients: Microstructural and Functional Asymmetry

Purpose: To assess possible association between intrinsic structural damage and clinical disability by correlating spinal cord diffusion-tensor (DT) imaging data with electrophysiological parameters in patients with a diagnosis of multiple sclerosis (MS). Materials and Methods: This study was approved by the local ethical committee according to the declaration of Helsinki and written informed consent was obtained. DT images and T1- and T2-weighted images of the spinal cord were acquired in 28 healthy volunteers and 41 MS patients. Fractional anisotropy (FA) and apparent diffusion coefficients were evaluated in normal-appearing white matter (NAWM) at the cervical level and were correlated with motor-evoked potentials (n = 34). Asymmetry index was calculated for FA values with corresponding left and right regions of interest as percentage of the absolute difference between these values relative to the sum of the respective FA values. Statistical analysis included Spearman rank correlations, Mann-Whitney test, and reliability analysis. Results: Healthy volunteers had low asymmetry index (1.5%-2.2%). In MS patients, structural abnormalities were reflected by asymmetric decrease of FA (asymmetry index: 3.6%; P = .15). Frequently asymmetrically affected among MS patients was left and right central motor conduction time (CMCT) to abductor digiti minimi muscle (ADMM) (asymmetry index, 15%-16%) and tibialis anterior muscle (TAM) (asymmetry index, 9.5%-14.1%). Statistically significant correlations of functional (ie, electrophysiological) and structural (ie, DT imaging) asymmetries were found (P = .005 for CMCT to ADMM; P = .007 for CMCT to TAM) for the cervical lateral funiculi, which comprise the crossed pyramidal tract. Interobserver reliability for DT imaging measurements was excellent (78%-87%). Conclusion: DT imaging revealed asymmetric anatomic changes in spinal cord NAWM, which corresponded to asymmetric electrophysiological deficits for both arms and legs, and reflected a specific structure-function relationship in the human spinal cord. © RSNA, 2013.

Non-exponential return time distributions for vorticity extremes explained by fractional Poisson processes

Serial correlation of extreme midlatitude cyclones observed at the storm track exits is explained by deviations from a Poisson process. To model these deviations, we apply fractional Poisson processes (FPPs) to extreme midlatitude cyclones, which are defined by the 850 hPa relative vorticity of the ERA interim reanalysis during boreal winter (DJF) and summer (JJA) seasons. Extremes are defined by a 99% quantile threshold in the grid-point time series. In general, FPPs are based on long-term memory and lead to non-exponential return time distributions. The return times are described by a Weibull distribution to approximate the Mittag–Leffler function in the FPPs. The Weibull shape parameter yields a dispersion parameter that agrees with results found for midlatitude cyclones. The memory of the FPP, which is determined by detrended fluctuation analysis, provides an independent estimate for the shape parameter. Thus, the analysis exhibits a concise framework of the deviation from Poisson statistics (by a dispersion parameter), non-exponential return times and memory (correlation) on the basis of a single parameter. The results have potential implications for the predictability of extreme cyclones.

Resolving diffusion in clay minerals at different time scales: Combination of experimental and modeling approaches

Since no single experimental or modeling technique provides data that allow a description of transport processes in clays and clay minerals at all relevant scales, several complementary approaches have to be combined to understand and explain the interplay between transport relevant phenomena. In this paper molecular dynamics simulations (MD) were used to investigate the mobility of water in the interlayer of montmorillonite (Mt), and to estimate the influence of mineral surfaces and interlayer ions on the water diffusion. Random Walk (RW) simulations based on a simplified representation of pore space in Mt were used to estimate and understand the effect of the arrangement of Mt particles on the meso- to macroscopic diffusivity of water. These theoretical calculations were complemented with quasielastic neutron scattering (QENS) measurements of aqueous diffusion in Mt with two pseudo-layers of water performed at four significantly different energy resolutions (i.e. observation times). The size of the interlayer and the size of Mt particles are two characteristic dimensions which determine the time dependent behavior of water diffusion in Mt. MD simulations show that at very short time scales water dynamics has the characteristic features of an oscillatory motion in the cage formed by neighbors in the first coordination shell. At longer time scales, the interaction of water with the surface determines the water dynamics, and the effect of confinement on the overall water mobility within the interlayer becomes evident. At time scales corresponding to an average water displacement equivalent to the average size of Mt particles, the effects of tortuosity are observed in the meso- to macroscopic pore scale simulations. Consistent with the picture obtained in the simulations, the QENS data can be described using a (local) 3D diffusion at short observation times, whereas at sufficiently long observation times a 2D diffusive motion is clearly observed. The effects of tortuosity measured in macroscopic tracer diffusion experiments are in qualitative agreement with RW simulations. By using experimental data to calibrate molecular and mesoscopic theoretical models, a consistent description of water mobility in clay minerals from the molecular to the macroscopic scale can be achieved. In turn, simulations help in choosing optimal conditions for the experimental measurements and the data interpretation. (C) 2014 Elsevier B.V. All rights reserved.

Real-time dynamics of open quantum spin systems driven by dissipative processes

We study the real-time evolution of large open quantum spin systems in two spatial dimensions, whose dynamics is entirely driven by a dissipative coupling to the environment. We consider different dissipative processes and investigate the real-time evolution from an ordered phase of the Heisenberg or XY model towards a disordered phase at late times, disregarding unitary Hamiltonian dynamics. The corresponding Kossakowski-Lindblad equation is solved via an efficient cluster algorithm. We find that the symmetry of the dissipative process determines the time scales, which govern the approach towards a new equilibrium phase at late times. Most notably, we find a slow equilibration if the dissipative process conserves any of the magnetization Fourier modes. In these cases, the dynamics can be interpreted as a diffusion process of the conserved quantity.

Saddlepoint Approximations to the Probability of Ruin in Finite Time for the Compound Poisson Risk Process Perturbed by Diffusion

A large deviations type approximation to the probability of ruin within a finite time for the compound Poisson risk process perturbed by diffusion is derived. This approximation is based on the saddlepoint method and generalizes the approximation for the non-perturbed risk process by Barndorff-Nielsen and Schmidli (Scand Actuar J 1995(2):169–186, 1995). An importance sampling approximation to this probability of ruin is also provided. Numerical illustrations assess the accuracy of the saddlepoint approximation using importance sampling as a benchmark. The relative deviations between saddlepoint approximation and importance sampling are very small, even for extremely small probabilities of ruin. The saddlepoint approximation is however substantially faster to compute.

Accelerated magnetic resonance diffusion tensor imaging of the median nerve using simultaneous multi-slice echo planar imaging with blipped CAIPIRINHA.

PURPOSE To investigate the feasibility of MR diffusion tensor imaging (DTI) of the median nerve using simultaneous multi-slice echo planar imaging (EPI) with blipped CAIPIRINHA. MATERIALS AND METHODS After federal ethics board approval, MR imaging of the median nerves of eight healthy volunteers (mean age, 29.4 years; range, 25-32) was performed at 3 T using a 16-channel hand/wrist coil. An EPI sequence (b-value, 1,000 s/mm(2); 20 gradient directions) was acquired without acceleration as well as with twofold and threefold slice acceleration. Fractional anisotropy (FA), mean diffusivity (MD) and quality of nerve tractography (number of tracks, average track length, track homogeneity, anatomical accuracy) were compared between the acquisitions using multivariate ANOVA and the Kruskal-Wallis test. RESULTS Acquisition time was 6:08 min for standard DTI, 3:38 min for twofold and 2:31 min for threefold acceleration. No differences were found regarding FA (standard DTI: 0.620 ± 0.058; twofold acceleration: 0.642 ± 0.058; threefold acceleration: 0.644 ± 0.061; p ≥ 0.217) and MD (standard DTI: 1.076 ± 0.080 mm(2)/s; twofold acceleration: 1.016 ± 0.123 mm(2)/s; threefold acceleration: 0.979 ± 0.153 mm(2)/s; p ≥ 0.074). Twofold acceleration yielded similar tractography quality compared to standard DTI (p > 0.05). With threefold acceleration, however, average track length and track homogeneity decreased (p = 0.004-0.021). CONCLUSION Accelerated DTI of the median nerve is feasible. Twofold acceleration yields similar results to standard DTI. KEY POINTS • Standard DTI of the median nerve is limited by its long acquisition time. • Simultaneous multi-slice acquisition is a new technique for accelerated DTI. • Accelerated DTI of the median nerve yields similar results to standard DTI.

Optimized CUDA-Based PDE Solver for Reaction Diffusion Systems on Arbitrary Surfaces

Partial differential equation (PDE) solvers are commonly employed to study and characterize the parameter space for reaction-diffusion (RD) systems while investigating biological pattern formation. Increasingly, biologists wish to perform such studies with arbitrary surfaces representing ‘real’ 3D geometries for better insights. In this paper, we present a highly optimized CUDA-based solver for RD equations on triangulated meshes in 3D. We demonstrate our solver using a chemotactic model that can be used to study snakeskin pigmentation, for example. We employ a finite element based approach to perform explicit Euler time integrations. We compare our approach to a naive GPU implementation and provide an in-depth performance analysis, demonstrating the significant speedup afforded by our optimizations. The optimization strategies that we exploit could be generalized to other mesh based processing applications with PDE simulations.