A magneto-viscoelasticity problem with a singular memory kernel
Carillo Sandra
;
Chipot Michel
;
Valente Vanda
;
Caffarelli Giorgio Vergara
The existence of solutions to a one-dimensional problem arising in magneto-viscoelasticity is here considered. Specifically, a non-linear system of integro-differential equations is analysed; it is obtained coupling an integro-differential equation modelling the viscoelastic behaviour, in which the kernel represents the relaxation function, with the non-linear partial differential equations modelling the presence of a magnetic field. The case under investigation generalizes a previous study since the relaxation function is allowed to be unbounded at the origin, provided it belongs to L-1; the magnetic model equation adopted, as in the previous results (Garillo et al., 2011, 2012; Chipot et al. 2008, 2009) is the penalized Ginzburg-Landau magnetic evolution equation. (C) 2016 Elsevier Ltd. All rights reserved.
Magneto-viscoelastic materials
Nonlinear integro-differential problem
Materials with memory
Singular kernel
For kernels zi which are positive and integrable we show that the operator g bar right arrow J(v)g = integral(x)(0) v(x-s)g(s)ds on a finite time interval enjoys a regularizing effect when applied to Holder continuous and Lebesgue functions and a "contractive" effect when applied to Sobolev functions. For Holder continuous functions, we establish that the improvement of the regularity of the modulus of continuity is given by the integral of the kernel, namely by the factor N(x) = integral(x)(0) v(s)ds. For functions in Lebesgue spaces, we prove that an improvement always exists, and it can be expressed in terms of Orlicz integrability. Finally, for functions in Sobolev spaces, we show that the operator J. "shrinks" the norm of the argument by a factor that, as in the Holder case, depends on the function N (whereas no regularization result can be obtained).
These results can be applied, for instance, to Abel kernels and to the Volterra function Z(x) = mu(x,0, -1) = integral(infinity)(0)x(s-1)/Gamma(s)ds, the latter being relevant for instance in the analysis of the Schrodinger equation with concentrated nonlinearities in R-2.
Volterra functions
Singular kernels
Volterra integral equations
Sonine kernels
The dynamics of thermally fluctuating conserved order parameters are described by stochastic conservation laws. Thermal equilibrium in such systems requires the dissipative and stochastic components of the flux to be related by detailed balance. Preserving this relation in spatial and temporal discretization is necessary to obtain solutions that have fidelity to the continuum. Here, we propose a finite-difference discretization that preserves the detailed balance on the lattice, has a spatial error that is isotropic to leading order in lattice spacing, and can be integrated accurately in time using a delayed difference method. We benchmark the method for model B dynamics with a phi(4) Landau free energy and obtain excellent agreement with the analytical results.
We present a novel application of the Lattice Boltzmann Method to the study of pulsed reactive flows in transitional Knudsen number regimes, namely 0.1 < Kn < 1.
We characterize the conversion efficiency of catalytic particles for different geometries and configurations, including single catalytic particle and nanoporous gold (npAu) spheres, within pulsed-flow reactors.
For all the explored configurations, the reactivity is found to increase with the Knudsen number of the flow, consistently with previous theoretical models and in reasonable agreement with experimental results in the literature.
Pulsed Reactive Flow
Lattice Boltzmann
Nanoporous Gold Catalyst
TAP Experiments
The surface structure and composition of a multi-component catalyst are critical factors in determining its catalytic performance. The surface composition can depend on the local pressure of the reacting species, leading to the possibility that the flow through a nanoporous catalyst can affect its structure and reactivity. Here, we explore this possibility for oxidation reactions on nanoporous gold, an AgAu bimetallic catalyst. We use microscopy and digital reconstruction to obtain the morphology of a two-dimensional slice of a nanoporous gold sample. Using lattice Boltzmann fluid dynamics simulations along with thermodynamic models based on first-principles total-energy calculations, we show that some sections of this sample have low local O-2 partial pressures when exposed to reaction conditions, which leads to a pure Au surface in these regions, instead of the active bimetallic AgAu phase. We also explore the effect of temperature on the surface structure and find that moderate temperatures (approximate to 300-450 K) should result in the highest intrinsic catalytic performance, in apparent agreement with experimental results. Published by AIP Publishing.
LATTICE-BOLTZMANN METHOD; AU ALLOY CATALYSTS; PEM FUEL-CELL; CO OXIDATION; BIMETALLIC NANOPARTICLES; HETEROGENEOUS CATALYSIS; AU(321) SURFACE; LOW-TEMPERATURE; GOLD CATALYSTS; SILVER
This paper presents a methodology to generate maps of atmosphere's precipitable water vapor (PWV) over large areas with a length of hundreds of kilometers and a width of about 250 km, based on the use of interferometric Sentinel-1A/BC-band synthetic aperture radar (SAR) data with a high spatial resolution of 5 x 20 m(2) and the revisiting time of six days. An algorithm to calibrate and merge PWV maps from different swaths of Sentinel-1 acquired along the same track, using global navigation satellite system (GNSS) measurements, is described. The proposed methodology is tested on Sentinel-1A SAR images acquired over the Iberian Peninsula, along both descending and ascending tracks. The assessment with an independent set of GNSS measurements shows a mean difference of a fraction of millimeter and a dispersion lower than 2 mm. Both the use of Sentinel-1A/B SAR images and the proposed methodology open new perspectives on the application of SAR meteorology for the high-resolution mapping of PWV over large region-spanning areas and the assimilation of interferometric SAR data into numerical weather models.
Global navigation satellite system (GNSS)
moisture content
precipitable water vapor (PWV)
synthetic aperture radar (SAR)
SAR interferometry (InSAR)
Sentinel-1
Objective: Neuroimaging studies provide evidence of disturbed resting-state brain networks in Schizophrenia (SZ). However, untangling the neuronal mechanisms that subserve these baseline alterations requires measurement of their electrophysiological underpinnings. This systematic review specifically investigates the contributions of resting-state Magnetoencephalography (MEG) in elucidating abnormal neural organization in SZ patients.& para;& para;Method: A systematic literature review of resting-state MEG studies in SZ was conducted. This literature is discussed in relation to findings from resting-state fMRI and EEG, as well as to task-based MEG research in SZ population. Importantly, methodological limitations are considered and recommendations to overcome current limitations are proposed.& para;& para;Results: Resting-state MEG literature in SZ points towards altered local and long-range oscillatory network dynamics in various frequency bands. Critical methodological challenges with respect to experiment design, and data collection and analysis need to be taken into consideration.& para;& para;Conclusion: Spontaneous MEG data show that local and global neural organization is altered in SZ patients. MEG is a highly promising tool to fill in knowledge gaps about the neurophysiology of SZ. However, to reach its fullest potential, basic methodological challenges need to be overcome.& para;& para;Significance: MEG-based resting-state power and connectivity findings could be great assets to clinical and translational research in psychiatry, and SZ in particular. (C) 2017 International Federation of Clinical Neurophysiology. Published by Elsevier Ireland Ltd. All rights reserved.
2017Presentazione / Comunicazione non pubblicata (convegno, evento, webinar...)metadata only access
NeuroPycon: A Python-based package for advanced MEG, EEG and fMRI connectivity analyses
David Meunier
;
Annalisa Pascarella
;
Daphné BertrandDubois
;
Jordan Alves
;
Fanny Barlaam
;
Arthur Dehgan
;
Tarek Lajnef
;
Etienne Combrisson
;
Dmitrii Altukhov
;
Karim Jerbi
NeuroPycon is an open-source multi-modal brain data analysis kit which provides Python-based pipelines for advanced multi-thread processing of fMRI, MEG and EEG data, with a focus on connectivity and graph analyses [1].
NeuroPycon is based on NiPype framework [2] which facilitates data analyses by wrapping many commonly-used neuroimaging software into a common python framework. Therefore, a major strength of NeuroPycon is that it relies on (and interfaces with) several freely available Python packages developed for efficient and fast parallel processing and that it seamlessly connects with existing open-science neuroimaging and signal processing toolboxes.
The flexible design allows users to configure analysis pipelines defined by connecting different nodes, where each node may be a user-defined function or a well-established tool or python-wrapped module (e.g. MNE-python for MEG analysis [3], etc.).
The current implementation of NeuroPycon contains three different packages:
- ephypype includes pipelines for electrophysiology analysis; current implementations allow for MEG/EEG data import, data pre-processing and cleaning by an automatic removal of eyes and heart related artefacts, sensor or source-level connectivity analyses
- graphpype allows to study functional connectivity exploiting graph-theoretical metrics including also modular partitions
- clipype is a command line interface for ephypype package.
NeuroPycon will shortly be available for download via github (installation via Docker) and is currently being documented. Future developments include fusion of multi-modal data (ex. MEG and fMRI or iEEG and fMRI).
References
1. Bullmore, Sporns (2009), Nat Rev Neurosci
2. Gorgolewski et al. (2011) Front. Neuroinform
3. Gramfort et al. (2013), Front. Neurosci
2017Poster in Atti di convegnometadata only access
The impact of MEG source reconstruction method on source-space connectivity estimation: A comparison between minimum-norm solution and beamforming
Ana Sofia Hincapie
;
Jan Kujala
;
Jérémie Mattout
;
Annalisa Pascarella
;
Sebastien Daligault
;
Claude Delpuech
;
Domingo Mery
;
Diego Cosmelli
;
Karim Jerbi
The effect of the choice of the inverse method on the cortico-cortical coupling analysis has been largely overlooked in the literature. Here, we set out to investigate the impact of three inverse methods on source coherence detection using simulated MEG data. To this end, we created thousands of randomly located pairs of sources and varied their inter- and intra-source correlation strength, source size and spatial configuration. Then, we used the simulated pairs of sources to generate sensor-level MEG measurements at varying signal-to-noise ratios (SNR). Next, we reconstructed the sources using L2-Minimum-Norm Estimate (MNE), Linearly Constrained Minimum Variance (LCMV) beamforming, and Dynamic Imaging of Coherent Sources (DICS) beamforming; and calculated source level power and coherence maps. We evaluated the performance of the methods using the Receiver Operating Characteristic (ROC) curves. The results indicate that beamformers perform better than MNE for coherence reconstructions of interacting point-like sources; but MNE provides better connectivity estimation than beamformers of interacting extended cortical patches, if each patch consists of dipoles with identical time series (high intra-patch coherence). However, the performance of the beamformers for interacting patches improves substantially if each cortical patch is simulated with partly coherent time series (partial intra-patch coherence). These results demonstrate that the choice of the inverse method impacts the results of MEG source-space coherence analysis, and that the optimal choice of the inverse solution depends on the spatial and synchronization profile of the interacting cortical sources. Our conclusions can guide method selection and help improve data interpretation regarding MEG connectivity estimation.
2017Poster in Atti di convegnometadata only access
Large-scale brain integration patterns differ in focused-attention and open-monitoring meditation
Daphné BertrandDubois
;
David Meunier
;
Annalisa Pascarella
;
Vittorio Pizzella
;
Laura Marzetti
;
Karim Jerbi
An important process underlying meditation and its benefits involves the regulation of attention. Although the two main meditation categories - open-monitoring meditation (OMM) and focused-attention meditation (FAM) - are associated with different benefits and attentional processes, direct comparisons between the attentional neural mechanism of FAM and OMM are rare. This study uses magnetoencephalography (MEG) recordings in 12 expert meditators to compare FAM and OMM by assessing (i) source spectral power, (ii) seed-based functional connectivity of key regions in attention, (including anterior cingulate cortex, dorsolateral prefrontal cortex and the thalamus) and (iii) graph theory metrics that describe brain-wide efficiency of information processing. We reconstructed the source space using minimum norm estimate and computed spectral power and functional connectivity in multiple frequency bands (delta, theta, alpha, beta, gamma) using a custom-designed python-based MEG analysis pipeline (NeuroPycon). The results reveal unique patterns of neural processes specific to FAM or OMM. Among other things, compared to FAM, OMM appears to be characterized by enhanced small-world network properties. By contrast, FAM exhibits greater functional connectivity between the anterior cingulate cortex and frontal regions. These findings shed light onto the mechanisms that potentially mediate the different behavioral and attentional capacities associated with each of the two meditation techniques. Our results are discussed in the context of previous behavioral and fMRI studies on meditation and attention.
Under the current deluge of omics, module networks distinctively emerge as methods capable of not only identifying inherently coherent groups (modules), thus reducing dimensionality, but also hypothesizing cause-effect relationships between modules and their regulators. Module networks were first designed in the transcriptomic era and further exploited in the multi-omic context to assess (for example) miRNA regulation of gene expression. Despite a number of available implementations, expansion of module networks to other omics is constrained by a limited characterization of the solutions' (modules plus regulators) accuracy and stability-an immediate need for the better characterization of molecular biology complexity in silico. We hence carefully assessed for LemonTree-a popular and open source module network implementation-the dependency of the software performances (sensitivity, specificity, false discovery rate, solutions' stability) on the input parameters and on the data quality (sample size, expression noise) based on synthetic and real data. In the process, we uncovered and fixed an issue in the code for the regulator assignment procedure. We concluded this evaluation with a table of recommended parameter settings. Finally, we applied these recommended settings to gut-intestinal metagenomic data from rheumatoid arthritis patients, to characterize the evolution of the gut-intestinal microbiome under different pharmaceutical regimens (methotrexate and prednisone) and we inferred innovative clinical recommendations with therapeutic potential, based on the computed module network.
Motivation: Cells derived by cellular engineering, i.e. differentiation of induced pluripotent stem cells and direct lineage reprogramming, carry a tremendous potential for medical applications and in particular for regenerative therapies. These approaches consist in the definition of lineage-specific experimental protocols that, by manipulation of a limited number of biological cues-niche mimicking factors, (in) activation of transcription factors, to name a few-enforce the final expression of cell-specific (marker) molecules. To date, given the intricate complexity of biological pathways, these approaches still present imperfect reprogramming fidelity, with uncertain consequences on the functional properties of the resulting cells. Results: We propose a novel tool eegc to evaluate cellular engineering processes, in a systemic rather than marker-based fashion, by integrating transcriptome profiling and functional analysis. Our method clusters genes into categories representing different states of (trans) differentiation and further performs functional and gene regulatory network analyses for each of the categories of the engineered cells, thus offering practical indications on the potential lack of the reprogramming protocol.
DEFINED FACTORS; TRANSCRIPTION FACTORS; ENDOTHELIAL-CELLS; DIRECT CONVERSION; STEM-CELLS; FIBROBLASTS; DIFFERENTIATION; IDENTITY; NETWORKS; NEURONS
In this paper the Wiener estimator for signal-denoising is generalized to
finite frame operators. In particular, a two-stage procedure which results
in a non-linear and non-diagonal estimator is proposed. Advantages and
disadvantages with respect to the classical Wiener estimator used with
orthonormal basis operator are discussed showing results on standard and
real test signals.
The numerical study presented in Part I (Dubbioso et al., 2017) on the bearing loads developed by the propellers
of a twin screw model during quasi-steady conditions is extended to transient maneuvers. In the previous study,
numerical simulations highlighted that the hydrodynamic loads might experience significant peak at moderate
turning rates due to complex interaction of the propeller with the wake. In the present paper, the complete turning
circle maneuver at ? 1/4 35 ? at Fr 1/4 0:265 is numerically simulated in order to analyze the character of the blade
loads during the transient phases after the actuation of the rudder (start and pull-out). The analysis shows that the
overall degradation of the propeller performance may occur also at kinematic conditions weaker than those
usually considered as the most critical ones (in general, tight maneuvers); therefore, these conditions should be
accounted for also in the early design phases.
Infection by Leishmania protozoan parasites can cause a variety of disease outcomes in humans and other mammals, from single self-healing cutaneous lesions to a visceral dissemination of the parasite. The correlation between chronic lesions and ecto-nucleotidase enzymes activity on the surface of the parasite is addressed here using damage caused in epithelial cells by nitric oxide. In order to explore the role of purinergic metabolism in lesion formation and the outcome of the infection, we implemented a cellular automata/lattice gas model involving major immune characters (Th1 and Th2 cells, IFN-gamma, IL-4, IL-12, adenosine-Ado-, NO) and parasite players for the dynamic analysis of the disease progress. The model were analyzed using partial ranking correlation coefficient (PRCC) to indicate the components that most influence the disease progression. Results show that low Ado inhibition rate over Th-cells is shared by L. major and L. braziliensis, while in L. amazonensis infection the Ado inhibition rate over Th-cells reaches 30%. IL-4 inhibition rate over Th-cell priming to Th1 independent of IL-12 are exclusive of L. major. The lesion size and progression showed agreement with published biological data and the model was able to simulate cutaneous leishmaniasis outcomes. The sensitivity analysis suggested that Ado inhibition rate over Th-cells followed by Leishmania survival probability were the most important characteristics of the process, with PRCC of 0.89 and 0.77 respectively. The simulations also showed a non-linear relationship between Ado inhibition rate over Th-cells and lesion size measured as number of dead epithelial cells. In conclusion, this model can be a useful tool for the quantitative understanding of the immune response in leishmaniasis.
leishmaniasis
cutaneous
adenosine (Ado)
model
lattice-gas
inflammation
A comprehensive approach to Sobolev-type embeddings, involving arbitrary rearrangement-
invariant norms on the entire Euclidean space R^n, is offered. In particular, the optimal target space in
any such embedding is exhibited. Crucial in our analysis is a new reduction principle for the relevant
embeddings, showing their equivalence to a couple of considerably simpler one-dimensional inequalities.
Applications to the classes of the Orlicz-Sobolev and the Lorentz-Sobolev spaces are also presented.
These contributions fill in a gap in the existing literature, where sharp results in such a general setting
are only available for domains of finite measure.
A comprehensive approach to Sobolev-type embeddings, involving arbitrary rearrangement-
invariant norms on the entire Euclidean space R^n, is offered. In particular, the optimal target space in
any such embedding is exhibited. Crucial in our analysis is a new reduction principle for the relevant
embeddings, showing their equivalence to a couple of considerably simpler one-dimensional inequalities.
Applications to the classes of the Orlicz-Sobolev and the Lorentz-Sobolev spaces are also presented.
These contributions fill in a gap in the existing literature, where sharp results in such a general setting
are only available for domains of finite measure.