Kinematic reduction of reaction-diffusion fronts with multiplicative noise. Derivation of stochastic sharp-interface equations

We study the dynamics of generic reaction-diffusion fronts, including pulses and chemical waves, in the presence of multiplicative noise. We discuss the connection between the reaction-diffusion Langevin-like field equations and the kinematic (eikonal) description in terms of a stochastic moving-boundary or sharp-interface approximation. We find that the effective noise is additive and we relate its strength to the noise parameters in the original field equations, to first order in noise strength, but including a partial resummation to all orders which captures the singular dependence on the microscopic cutoff associated with the spatial correlation of the noise. This dependence is essential for a quantitative and qualitative understanding of fluctuating fronts, affecting both scaling properties and nonuniversal quantities. Our results predict phenomena such as the shift of the transition point between the pushed and pulled regimes of front propagation, in terms of the noise parameters, and the corresponding transition to a non-Kardar-Parisi-Zhang universality class. We assess the quantitative validity of the results in several examples including equilibrium fluctuations and kinetic roughening. We also predict and observe a noise-induced pushed-pulled transition. The analytical predictions are successfully tested against rigorous results and show excellent agreement with numerical simulations of reaction-diffusion field equations with multiplicative noise.

Vortices in atomic Bose-Einstein condensates in the large-gas-parameter region

In this work we compare the results of the Gross-Pitaevskii and modified Gross-Pitaevskii equations with ab initio variational Monte Carlo calculations for Bose-Einstein condensates of atoms in axially symmetric traps. We examine both the ground state and excited states having a vortex line along the z axis at high values of the gas parameter and demonstrate an excellent agreement between the modified Gross-Pitaevskii and ab initio Monte Carlo methods, both for the ground and vortex states.

Pairing in a two-component ultracold Fermi gas: Phases with broken-space symmetries

We explore the phase diagram of a two-component ultracold atomic Fermi gas interacting with zero-range forces in the limit of weak coupling. We focus on the dependence of the pairing gap and the free energy on the variations in the number densities of the two species while the total density of the system is held fixed. As the density asymmetry is increased, the system exhibits a transition from a homogenous Bardeen-Cooper-Schrieffer (BCS) phase to phases with spontaneously broken global space symmetries. One such realization is the deformed Fermi surface superfluidity (DFS) which exploits the possibility of deforming the Fermi surfaces of the species into ellipsoidal form at zero total momentum of Cooper pairs. The critical asymmetries at which the transition from DFS to the unpaired state occurs are larger than those for the BCS phase. In this precritical region the DFS phase lowers the pairing energy of the asymmetric BCS state. We compare quantitatively the DFS phase to another realization of superconducting phases with broken translational symmetry: the single-plane-wave Larkin-Ovchinnikov-Fulde-Ferrell phase, which is characterized by a nonvanishing center-of-mass momentum of the Cooper pairs. The possibility of the detection of the DFS phase in the time-of-flight experiments is discussed and quantified for the case of 6Li atoms trapped in two different hyperfine states.

Ground-state properties of a dilute homogeneous Bose gas of hard disks in two dimensions

The energy and structure of a dilute hard-disks Bose gas are studied in the framework of a variational many-body approach based on a Jastrow correlated ground-state wave function. The asymptotic behaviors of the radial distribution function and the one-body density matrix are analyzed after solving the Euler equation obtained by a free minimization of the hypernetted chain energy functional. Our results show important deviations from those of the available low density expansions, already at gas parameter values x~0.001 . The condensate fraction in 2D is also computed and found generally lower than the 3D one at the same x.

Modelization of surface diffusion of a molecular dimer

A simple model for a dimer molecular diffusion on a crystalline surface, as a function of temperature, is presented. The dimer is formed by two particles coupled by a quadratic potential. The dimer diffusion is modeled by an overdamped Langevin equation in the presence of a two-dimensional periodic potential. Numerical simulation¿s results exhibit some dynamical properties observed, for example, in Si2 diffusion on a silicon [100] surface. They can be used to predict the value of the effective friction parameter. Comparison between our model and experimental measurements is presented.

Reaction-diffusion fronts under stochastic advection

We study front propagation in stirred media using a simplified modelization of the turbulent flow. Computer simulations reveal the existence of the two limiting propagation modes observed in recent experiments with liquid phase isothermal reactions. These two modes respectively correspond to a wrinkled although sharp propagating interface and to a broadened one. Specific laws relative to the enhancement of the front velocity in each regime are confirmed by our simulations.

Mismatch of arterial and central venous blood gas analysis during haemorrhage.

BACKGROUND AND OBJECTIVE: Arterial base excess and lactate levels are key parameters in the assessment of critically ill patients. The use of venous blood gas analysis may be of clinical interest when no arterial blood is available initially. METHODS: Twenty-four pigs underwent progressive normovolaemic haemodilution and subsequent progressive haemorrhage until the death of the animal. Base excess and lactate levels were determined from arterial and central venous blood after each step. In addition, base excess was calculated by the Van Slyke equation modified by Zander (BE(z)). Continuous variables were summarized as mean +/- SD and represent all measurements (n = 195). RESULTS: Base excess according to National Committee for Clinical Laboratory Standards for arterial blood was 2.27 +/- 4.12 versus 2.48 +/- 4.33 mmol(-l) for central venous blood (P = 0.099) with a strong correlation (r(2) = 0.960, P < 0.001). Standard deviation of the differences between these parameters (SD-DIFBE) did not increase (P = 0.355) during haemorrhage as compared with haemodilution. Arterial lactate was 2.66 +/- 3.23 versus 2.71 +/- 2.80 mmol(-l) in central venous blood (P = 0.330) with a strong correlation (r(2) = 0.983, P < 0.001). SD-DIFLAC increased (P < 0.001) during haemorrhage. BE(z) for central venous blood was 2.22 +/- 4.62 mmol(-l) (P = 0.006 versus arterial base excess according to National Committee for Clinical Laboratory Standards) with strong correlation (r(2) = 0.942, P < 0.001). SD-DIFBE(z)/base excess increased (P < 0.024) during haemorrhage. CONCLUSION: Central venous blood gas analysis is a good predictor for base excess and lactate in arterial blood in steady-state conditions. However, the variation between arterial and central venous lactate increases during haemorrhage. The modification of the Van Slyke equation by Zander did not improve the agreement between central venous and arterial base excess.

Diffusion of passive scalars under stochastic convection

The diffusion of passive scalars convected by turbulent flows is addressed here. A practical procedure to obtain stochastic velocity fields with well¿defined energy spectrum functions is also presented. Analytical results are derived, based on the use of stochastic differential equations, where the basic hypothesis involved refers to a rapidly decaying turbulence. These predictions are favorable compared with direct computer simulations of stochastic differential equations containing multiplicative space¿time correlated noise.

Lattice-gas model of particles with orientational and positional degrees of freedom: Mean-field treatment

We consider a lattice-gas model of particles with internal orientational degrees of freedom. In addition to antiferromagnetic nearest-neighbor (NN) and next-nearest-neighbor (NNN) positional interactions we also consider NN and NNN interactions arising from the internal state of the particles. The system then shows positional and orientational ordering modes with associated phase transitions at Tp and To temperatures at which long-range positional and orientational ordering are, respectively, lost. We use mean-field techniques to obtain a general approach to the study of these systems. By considering particular forms of the orientational interaction function we study coupling effects between both phase transitions arising from the interplay between orientational and positional degrees of freedom. In mean-field approximation coupling effects appear only for the phase transition taking place at lower temperatures. The strength of the coupling depends on the value of the long-range order parameter that remains finite at that temperature.

Lattice-gas model of orientable molecules: Application to liquid crystals

We have analyzed a two-dimensional lattice-gas model of cylindrical molecules which can exhibit four possible orientations. The Hamiltonian of the model contains positional and orientational energy interaction terms. The ground state of the model has been investigated on the basis of Karl¿s theorem. Monte Carlo simulation results have confirmed the predicted ground state. The model is able to reproduce, with appropriate values of the Hamiltonian parameters, both, a smectic-nematic-like transition and a nematic-isotropic-like transition. We have also analyzed the phase diagram of the system by mean-field techniques and Monte Carlo simulations. Mean-field calculations agree well qualitatively with Monte Carlo results but overestimate transition temperatures.

RNA interference against aquaporin-4 decreases the apparent diffusion coefficient in the normal brain

Aim: Diffusion weighted magnetic resonance imaging (MRI) is now widely used in human brain diagnosis.1 To date molecular mechanisms underlying changes in Apparent Diffusion Coefficient (ADC) signals remain poorly understood. AQP4, localized to astrocytes, is one of the most highly expressed cerebral AQPs.2 AQP4 is involved in water movement within the cell membrane of cultured astrocytes.3 We hypothesize that AQP4 contributes to water diffusion and underlying ADC values in normal brain. Methods: We used an RNA interference (RNAi) protocol in vivo, to acutely knockdown expression of AQP4 in rat brain and to determine whether this was associated with changes in brain ADC values using MRI protocols as previously described.4 RNAi was performed using specific small interference RNA (siRNA) against AQP4 (siAQP4) and a non-targeted-siRNA (siGLO) as a control. The specificity and efficiency of the siAQP4 were first tested in vitro in astrocyte and hippocampal slice cultures. In vivo, siRNAs were injected into the rat cortex 3d prior to MRI acquisition and AQP4 was assessed by western blot (n=4) and immunohistochemistry (n=6). Histology was performed on adjacent slices. Results: siAQP4 application on primary astrocyte cultures induced a 76% decrease in AQP4 expression after 4 days. In hippocampal slice cultures; we also found a significant decrease in AQP4 expression in astrocytes after siAQP4. In vivo, injection of non-targeted siRNA (siGLO) tagged with CY3 allowed us to show that GFAP positive cells (astrocytes) were positively stained with CY3-siGLO, showing efficient transfection. Western blot and immunohistochemical analysis showed that siAQP4 induced a ~30% decrease in AQP4 expression without modification of tissue properties or cell death. After siAQP4 treatment, a significant decrease in ADC values (~50%) were observed without altered of T2 values. Conclusions: Together these results suggest that AQP4 reduces water diffusion through the astrocytic plasma membrane and decreases ADC values. Our findings demonstrate for the first time that astrocytic AQP4 contributes significantly to brain water diffusion and ADC values in normal brain. These results open new avenues to interpretation of ADC values under normal physiological conditions and in acute and chronic brain injuries.

Interpretation of primary care physicians' attitude regarding rotavirus immunisation using diffusion of innovation theories.

To evaluate primary care physicians' attitude towards implementation of rotavirus (RV) immunisation into the Swiss immunisation schedule, an eight-question internet-based questionnaire was sent to the 3799 subscribers of InfoVac, a nationwide web-based expert network on immunisation issues, which reaches >95% of paediatricians and smaller proportions of other primary care physicians. Five demographic variables were also inquired. Descriptive statistics and multivariate analyses for the main outcome "acceptance of routine RV immunisation" and other variables were performed. Diffusion of innovation theory was used for data assessment. Nine-hundred seventy-seven questionnaires were returned (26%). Fifty percent of participants were paediatricians. Routine RV immunisation was supported by 146 participants (15%; so called early adopters), dismissed by 620 (64%), leaving 211 (21%) undecided. However, when asked whether they would recommend RV vaccination to parents if it were officially recommended by the federal authorities and reimbursed, 467 (48.5%; so called early majority) agreed to recommend RV immunisation. Multivariate analysis revealed that physicians who would immunise their own child (OR: 5.1; 95% CI: 4.1-6.3), hospital-based physicians (OR: 1.6; 95% CI: 1.1-2.3) and physicians from the French (OR: 1.6; 95% CI: 1.2-2.3) and Italian speaking areas of Switzerland (OR: 2.5; 95% CI: 1.1-5.8) were more likely to support RV immunisation. Diffusion of innovation theory predicts a >80% implementation if approximately 50% of a given population support an innovation. Introduction of RV immunisation in Switzerland is likely to be successful, if (i) the federal authorities issue an official recommendation and (ii) costs are covered by basic health care insurance.

Accelerated Microstructure Imaging via Convex Optimization (AMICO) from diffusion MRI data.

Microstructure imaging from diffusion magnetic resonance (MR) data represents an invaluable tool to study non-invasively the morphology of tissues and to provide a biological insight into their microstructural organization. In recent years, a variety of biophysical models have been proposed to associate particular patterns observed in the measured signal with specific microstructural properties of the neuronal tissue, such as axon diameter and fiber density. Despite very appealing results showing that the estimated microstructure indices agree very well with histological examinations, existing techniques require computationally very expensive non-linear procedures to fit the models to the data which, in practice, demand the use of powerful computer clusters for large-scale applications. In this work, we present a general framework for Accelerated Microstructure Imaging via Convex Optimization (AMICO) and show how to re-formulate this class of techniques as convenient linear systems which, then, can be efficiently solved using very fast algorithms. We demonstrate this linearization of the fitting problem for two specific models, i.e. ActiveAx and NODDI, providing a very attractive alternative for parameter estimation in those techniques; however, the AMICO framework is general and flexible enough to work also for the wider space of microstructure imaging methods. Results demonstrate that AMICO represents an effective means to accelerate the fit of existing techniques drastically (up to four orders of magnitude faster) while preserving accuracy and precision in the estimated model parameters (correlation above 0.9). We believe that the availability of such ultrafast algorithms will help to accelerate the spread of microstructure imaging to larger cohorts of patients and to study a wider spectrum of neurological disorders.

A multi-center study: intra-scan and inter-scan variability of diffusion spectrum imaging.

The objective of this study was to investigate whether it is possible to pool together diffusion spectrum imaging data from four different scanners, located at three different sites. Two of the scanners had identical configuration whereas two did not. To measure the variability, we extracted three scalar maps (ADC, FA and GFA) from the DSI and utilized a region and a tract-based analysis. Additionally, a phantom study was performed to rule out some potential factors arising from the scanner performance in case some systematic bias occurred in the subject study. This work was split into three experiments: intra-scanner reproducibility, reproducibility with twin-scanner settings and reproducibility with other configurations. Overall for the intra-scanner and twin-scanner experiments, the region-based analysis coefficient of variation (CV) was in a range of 1%-4.2% and below 3% for almost every bundle for the tract-based analysis. The uncinate fasciculus showed the worst reproducibility, especially for FA and GFA values (CV 3.7-6%). For the GFA and FA maps, an ICC value of 0.7 and above is observed in almost all the regions/tracts. Looking at the last experiment, it was found that there is a very high similarity of the outcomes from the two scanners with identical setting. However, this was not the case for the two other imagers. Given the fact that the overall variation in our study is low for the imagers with identical settings, our findings support the feasibility of cross-site pooling of DSI data from identical scanners.

Perfusion and diffusion MRI of glioblastoma progression in a four-year prospective temozolomide clinical trial

RESUME BUT Cette étude a été menée sur le suivi de patients traités pour un glioblastome nouvellement diagnostiqué. Son objectif a été de déterminer l'impact des séquences de perfusion et de diffusion en imagerie par résonance magnétique (IRM). Un intérêt particulier a été porté au potentiel de ces nouvelles techniques d'imagerie dans l'anticipation de la progression de la maladie. En effet, l'intervalle de temps libre de progression est une mesure alternative de pronostic fréquemment utilisée. MATERIEL ET METHODE L'étude a porté sur 41 patients participant à un essai clinique de phase II de traitement par temozolomide. Leur suivi radiologique a comporté un examen IRM dans les 21 à 28 jours après radiochimiothérapie et tous les 2 mois par la suite. L'évaluation des images s'est faite sur la base de l'évaluation de l'effet de masse ainsi que de la mesure de la taille de la lésion sur les images suivantes : T1 avec produit de contraste, T2, diffusion, perfusion. Afin de déterminer la date de progression de la maladie, les critères classiques de variation de taille adjoints aux critères cliniques habituels ont été utilisés. RESULAT 311 examens IRM ont été revus. Au moment de la progression (32 patients), une régression multivariée selon Cox a permis de déterminer deux paramètres de survie : diamètre maximal en T1 (p>0.02) et variation de taille en T2 (p<0.05). L'impact de la perfusion et de la diffusion n'a pas été démontré de manière statistiquement significative. CONCLUSION Les techniques de perfusion et de diffusion ne peuvent pas être utilisées pour anticiper la progression tumorale. Alors que la prise de décision au niveau thérapeutique est critique au moment de la progression de la maladie, l'IRM classique en T1 et en T2 reste la méthode d'imagerie de choix. De manière plus spécifique, une prise de contraste en T1 supérieure à 3 cm dans son plus grand diamètre associée à un hypersignal T2 en augmentation forment un marqueur de mauvais pronostic.