Unified formalism for non-autonomous mechanical systems

We present a unified geometric framework for describing both the Lagrangian and Hamiltonian formalisms of regular and non-regular time-dependent mechanical systems, which is based on the approach of Skinner and Rusk (1983). The dynamical equations of motion and their compatibility and consistency are carefully studied, making clear that all the characteristics of the Lagrangian and the Hamiltonian formalisms are recovered in this formulation. As an example, it is studied a semidiscretization of the nonlinear wave equation proving the applicability of the proposed formalism.

A fully adaptive numerical approximation for a two-dimensional epidemic model with nonlinear cross-diffusion

An epidemic model is formulated by a reactionâeuro"diffusion system where the spatial pattern formation is driven by cross-diffusion. The reaction terms describe the local dynamics of susceptible and infected species, whereas the diffusion terms account for the spatial distribution dynamics. For both self-diffusion and cross-diffusion, nonlinear constitutive assumptions are suggested. To simulate the pattern formation two finite volume formulations are proposed, which employ a conservative and a non-conservative discretization, respectively. An efficient simulation is obtained by a fully adaptive multiresolution strategy. Numerical examples illustrate the impact of the cross-diffusion on the pattern formation.

Dielectric constant obtained from TDR and volumetric moisture of soils in southern Brazil

Soil moisture is the property which most greatly influences the soil dielectric constant, which is also influenced by soil mineralogy. The aim of this study was to determine mathematical models for soil moisture and the dielectric constant (Ka) for a Hapludalf, two clayey Hapludox and a very clayey Hapludox and test the reliability of universal models, such as those proposed by Topp and Ledieu and their co-workers in the 80's, and specific models to estimate soil moisture with a TDR. Soil samples were collected from the 0 to 0.30 m layer, sieved through a mesh of 0.002 m diameter and packed in PVC cylinders with a 0.1 m diameter and 0.3 m height. Seven samples of each soil class were saturated by capillarity and a probe composed of two rods was inserted in each one of them. Moisture readings began with the saturated soil and concluded when the soil was near permanent wilting point. In each step, the samples were weighed on a precision scale to calculate volumetric moisture. Linear and polynomial models were adjusted for each soil class and for all soils together between soil moisture and the dielectric constant. Accuracy of the models was evaluated by the coefficient of determination, the standard error of estimate and the 1:1 line. The models proposed by Topp and Ledieu and their co-workers were not adequate for estimating the moisture in the soil classes studied. The adjusted linear and polynomial models for the entire set of data of the four soil classes did not have sufficient accuracy for estimating soil moisture. The greater the soil clay and Fe oxide content, the greater the dielectric constant of the medium for a given volumetric moisture. The specific models, θ = 0.40283 - 0.04231 Ka + 0.00194 Ka² - 0.000022 Ka³ (Hapludox) θ = 0.01971 + 0.02902 Ka - 0.00086 Ka² + 0.000012 Ka³ (Hapludox -PF), θ = 0.01692 - 0.00507 Ka (Hapludalf) and θ = 0.08471 + 0.01145 Ka (Hapludox-CA), show greater accuracy and reliability for estimating soil moisture in the soil classes studied.

High b-value diffusion-weighted imaging: A sensitive method to reveal white matter differences in schizophrenia.

Over the last 10 years, diffusion-weighted imaging (DWI) has become an important tool to investigate white matter (WM) anomalies in schizophrenia. Despite technological improvement and the exponential use of this technique, discrepancies remain and little is known about optimal parameters to apply for diffusion weighting during image acquisition. Specifically, high b-value diffusion-weighted imaging known to be more sensitive to slow diffusion is not widely used, even though subtle myelin alterations as thought to happen in schizophrenia are likely to affect slow-diffusing protons. Schizophrenia patients and healthy controls were scanned with a high b-value (4000s/mm(2)) protocol. Apparent diffusion coefficient (ADC) measures turned out to be very sensitive in detecting differences between schizophrenia patients and healthy volunteers even in a relatively small sample. We speculate that this is related to the sensitivity of high b-value imaging to the slow-diffusing compartment believed to reflect mainly the intra-axonal and myelin bound water pool. We also compared these results to a low b-value imaging experiment performed on the same population in the same scanning session. Even though the acquisition protocols are not strictly comparable, we noticed important differences in sensitivities in the favor of high b-value imaging, warranting further exploration.

Relaxation times of non-Markovian processes

We consider a general class of non-Markovian processes defined by stochastic differential equations with Ornstein-Uhlenbeck noise. We present a general formalism to evaluate relaxation times associated with correlation functions in the steady state. This formalism is a generalization of a previous approach for Markovian processes. The theoretical results are shown to be in satisfactory agreement both with experimental data for a cubic bistable system and also with a computer simulation of the Stratonovich model. We comment on the dynamical role of the non-Markovianicity in different situations.

Fluctuations in domain growth: Ginzburg-Landau equations with multiplicative noise

Ginzburg-Landau equations with multiplicative noise are considered, to study the effects of fluctuations in domain growth. The equations are derived from a coarse-grained methodology and expressions for the resulting concentration-dependent diffusion coefficients are proposed. The multiplicative noise gives contributions to the Cahn-Hilliard linear-stability analysis. In particular, it introduces a delay in the domain-growth dynamics.

l0-deconvolution for compressive diffusion MRI

Diffusion MRI is a well established imaging modality providing a powerful way to non-invasively probe the structure of the white matter. Despite the potential of the technique, the intrinsic long scan times of these sequences have hampered their use in clinical practice. For this reason, a wide variety of methods have been proposed to shorten acquisition times. [...] We here review a recent work where we propose to further exploit the versatility of compressed sensing and convex optimization with the aim to characterize the fiber orientation distribution sparsity more optimally. We re-formulate the spherical deconvolution problem as a constrained l0 minimization.

Classical solutions of the leading-logarithm approximation with non-trivial topology

Exact solutions of the classical equations corresponding to the leading-logarithm approximation are obtained. They are classified by an (integer) topological number.

Diffusion on a solid surface: anomalous is normal

We present a numerical study of classical particles diffusing on a solid surface. The particles motion is modeled by an underdamped Langevin equation with ordinary thermal noise. The particle-surface interaction is described by a periodic or a random two-dimensional potential. The model leads to a rich variety of different transport regimes, some of which correspond to anomalous diffusion such as has recently been observed in experiments and Monte Carlo simulations. We show that this anomalous behavior is controlled by the friction coefficient and stress that it emerges naturally in a system described by ordinary canonical Maxwell-Boltzmann statistics.

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.

Can one or two high doses of oral vitamin D3 correct insufficiency in a non-supplemented rheumatologic population?

We evaluated the effectiveness of supplementation with high dose of oral vitamin D3 to correct vitamin D insufficiency. We have shown that one or two oral bolus of 300,000 IU of vitamin D3 can correct vitamin D insufficiency in 50% of patients and that the patients who benefited more from supplementation were those with the lowest baseline levels. INTRODUCTION: Adherence with daily oral supplements of vitamin D3 is suboptimal. We evaluated the effectiveness of a single high dose of oral vitamin D3 (300,000 IU) to correct vitamin D insufficiency in a rheumatologic population. METHODS: Over 1 month, 292 patients had levels of 25-OH vitamin D determined. Results were classified as: deficiency <10 ng/ml, insufficiency ≥10 to 30 ng/ml, and normal ≥30 ng/ml. We added a category using the IOM recommended cut-off of 20 ng/ml. Patients with deficient or normal levels were excluded, as well as patients already supplemented with vitamin D3. Selected patients (141) with vitamin D insufficiency (18.5 ng/ml (10.2-29.1) received a prescription for 300,000 IU of oral vitamin D3 and were asked to return after 3 (M3) and 6 months (M6). Patients still insufficient at M3 received a second prescription for 300,000 IU of oral vitamin D3. Relation between changes in 25-OH vitamin D between M3 and M0 and baseline values were assessed. RESULTS: Patients (124) had a blood test at M3. Two (2%) had deficiency (8.1 ng/ml (7.5-8.7)) and 50 (40%) normal results (36.7 ng/ml (30.5-5.5)). Seventy-two (58%) were insufficient (23.6 ng/ml (13.8-29.8)) and received a second prescription for 300,000 IU of oral vitamin D3. Of the 50/124 patients who had normal results at M3 and did not receive a second prescription, 36 (72%) had a test at M6. Seventeen (47%) had normal results (34.8 ng/ml (30.3-42.8)) and 19 (53%) were insufficient (25.6 ng/ml (15.2-29.9)). Of the 72/124 patients who receive a second prescription, 54 (75%) had a test at M6. Twenty-eight (52%) had insufficiency (23.2 ng/ml (12.8-28.7)) and 26 (48%) had normal results (33.8 ng/ml (30.0-43.7)). At M3, 84% patients achieved a 25-OH vitamin D level >20 ng/ml. The lowest the baseline value, the highest the change after 3 months (negative relation with a correlation coefficient r = -0.3, p = 0.0007). CONCLUSIONS: We have shown that one or two oral bolus of 300,000 IU of vitamin D3 can correct vitamin D insufficiency in 50% of patients.

Dynamical properties of nonmarkovian stochastic differential equations

We study nonstationary non-Markovian processes defined by Langevin-type stochastic differential equations with an OrnsteinUhlenbeck driving force. We concentrate on the long time limit of the dynamical evolution. We derive an approximate equation for the correlation function of a nonlinear nonstationary non-Markovian process, and we discuss its consequences. Non-Markovicity can introduce a dependence on noise parameters in the dynamics of the correlation function in cases in which it becomes independent of these parameters in the Markovian limit. Several examples are discussed in which the relaxation time increases with respect to the Markovian limit. For a Brownian harmonic oscillator with fluctuating frequency, the non-Markovicity of the process decreases the domain of stability of the system, and it can change an infradamped evolution into an overdamped one.

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.

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.