Diffusion in spatially and temporarily inhomogeneous media

In this paper we consider diffusion of a passive substance C in a temporarily and spatially inhomogeneous two-dimensional medium. As a realization for the latter we choose a phase-separating medium consisting of two substances A and B, whose dynamics is determined by the Cahn-Hilliard equation. Assuming different diffusion coefficients of C in A and B, we find that the variance of the distribution function of the said substance grows less than linearly in time. We derive a simple identity for the variance using a probabilistic ansatz and are then able to identify the interface between A and B as the main cause for this nonlinear dependence. We argue that, finally, for very large times the here temporarily dependent diffusion "constant" goes like t-1/3 to a constant asymptotic value D¿. The latter is calculated approximately by employing the effective-medium approximation and by fitting the simulation data to the said time dependence.

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.

Mean field model for spatially extended systems in the presence of multiplicative noise

We present a mean field model that describes the effect of multiplicative noise in spatially extended systems. The model can be solved analytically. For the case of the phi4 potential it predicts that the phase transition is shifted. This conclusion is supported by numerical simulations of this model in two dimensions.

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.

Coarsening of solid-liquid mixtures in a random acceleration field

The effects of flow induced by a random acceleration field (g-jitter) are considered in two related situations that are of interest for microgravity fluid experiments: the random motion of isolated buoyant particles, and diffusion driven coarsening of a solid-liquid mixture. We start by analyzing in detail actual accelerometer data gathered during a recent microgravity mission, and obtain the values of the parameters defining a previously introduced stochastic model of this acceleration field. The diffusive motion of a single solid particle suspended in an incompressible fluid that is subjected to such random accelerations is considered, and mean squared velocities and effective diffusion coefficients are explicitly given. We next study the flow induced by an ensemble of such particles, and show the existence of a hydrodynamically induced attraction between pairs of particles at distances large compared with their radii, and repulsion at short distances. Finally, a mean field analysis is used to estimate the effect of g-jitter on diffusion controlled coarsening of a solid-liquid mixture. Corrections to classical coarsening rates due to the induced fluid motion are calculated, and estimates are given for coarsening of Sn-rich particles in a Sn-Pb eutectic fluid, an experiment to be conducted in microgravity in the near future.

From lattice-gas models to nonlinear diffusion models: A derivation of macroscopic equations and fluctuations

We derive nonlinear diffusion equations and equations containing corrections due to fluctuations for a coarse-grained concentration field. To deal with diffusion coefficients with an explicit dependence on the concentration values, we generalize the Van Kampen method of expansion of the master equation to field variables. We apply these results to the derivation of equations of phase-separation dynamics and interfacial growth instabilities.

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.

Relativistic mean-field interaction with density-dependent meson-nucleon vertices based on microscopical calculations

Although ab initio calculations of relativistic Brueckner theory lead to large scalar isovector fields in nuclear matter, at present, successful versions of covariant density functional theory neglect the interactions in this channel. A new high-precision density functional DD-MEδ is presented which includes four mesons, σ, ω, δ, and ρ, with density-dependent meson-nucleon couplings. It is based to a large extent on microscopic ab initiocalculations in nuclear matter. Only four of its parameters are determined by adjusting to binding energies and charge radii of finite nuclei. The other parameters, in particular the density dependence of the meson-nucleon vertices, are adjusted to nonrelativistic and relativistic Brueckner calculations of symmetric and asymmetric nuclear matter. The isovector effective mass mp*−mn* derived from relativistic Brueckner theory is used to determine the coupling strength of the δ meson and its density dependence.

Origin of the neutron skin thickness of 208Pb in nuclear mean-field models

We study whether the neutron skin thickness Δrnp of 208Pb originates from the bulk or from the surface of the nucleon density distributions, according to the mean-field models of nuclear structure, and find that it depends on the stiffness of the nuclear symmetry energy. The bulk contribution to Δrnp arises from an extended sharp radius of neutrons, whereas the surface contribution arises from different widths of the neutron and proton surfaces. Nuclear models where the symmetry energy is stiff, as typical of relativistic models, predict a bulk contribution in Δrnp of 208Pb about twice as large as the surface contribution. In contrast, models with a soft symmetry energy like common nonrelativistic models predict that Δrnp of 208Pb is divided similarly into bulk and surface parts. Indeed, if the symmetry energy is supersoft, the surface contribution becomes dominant. We note that the linear correlation of Δrnp of 208Pb with the density derivative of the nuclear symmetry energy arises from the bulk part of Δrnp. We also note that most models predict a mixed-type (between halo and skin) neutron distribution for 208Pb. Although the halo-type limit is actually found in the models with a supersoft symmetry energy, the skin-type limit is not supported by any mean-field model. Finally, we compute parity-violating electron scattering in the conditions of the 208Pb parity radius experiment (PREX) and obtain a pocket formula for the parity-violating asymmetry in terms of the parameters that characterize the shape of the 208Pb nucleon densities.

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.

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.