A mixed finite element method for nonlinear diffusion equations

We propose a mixed finite element method for a class of nonlinear diffusion equations, which is based on their interpretation as gradient flows in optimal transportation metrics. We introduce an appropriate linearization of the optimal transport problem, which leads to a mixed symmetric formulation. This formulation preserves the maximum principle in case of the semi-discrete scheme as well as the fully discrete scheme for a certain class of problems. In addition solutions of the mixed formulation maintain exponential convergence in the relative entropy towards the steady state in case of a nonlinear Fokker-Planck equation with uniformly convex potential. We demonstrate the behavior of the proposed scheme with 2D simulations of the porous medium equations and blow-up questions in the Patlak-Keller-Segel model.

Oscillations in a molecular structured cell population model

We consider a nonlinear cyclin content structured model of a cell population divided into proliferative and quiescent cells. We show, for particular values of the parameters, existence of solutions that do not depend on the cyclin content. We make numerical simulations for the general case obtaining, for some values of the parameters convergence to the steady state but also oscillations of the population for others.

A decision-making Fokker-Planck model in computational neuroscience

Minimal models for the explanation of decision-making in computational neuroscience are based on the analysis of the evolution for the average firing rates of two interacting neuron populations. While these models typically lead to multi-stable scenario for the basic derived dynamical systems, noise is an important feature of the model taking into account finite-size effects and robustness of the decisions. These stochastic dynamical systems can be analyzed by studying carefully their associated Fokker-Planck partial differential equation. In particular, we discuss the existence, positivity and uniqueness for the solution of the stationary equation, as well as for the time evolving problem. Moreover, we prove convergence of the solution to the the stationary state representing the probability distribution of finding the neuron families in each of the decision states characterized by their average firing rates. Finally, we propose a numerical scheme allowing for simulations performed on the Fokker-Planck equation which are in agreement with those obtained recently by a moment method applied to the stochastic differential system. Our approach leads to a more detailed analytical and numerical study of this decision-making model in computational neuroscience.

Outcomes and reoperations after total correction of complete atrio-ventricular septal defect

BACKGROUND: Surgical correction of complete atrio-ventricular septal defect (AVSD) achieves satisfactory results with low morbidity and mortality, but may require reoperation. Our recent operative results at mid-term were followed-up. METHODS: From June 2000 to December 2007, 81 patients (Down syndrome; n=60), median age 4.0 months (range 0.7-118.6) and weight 4.7kg (range 2.2-33), underwent complete AVSD correction. Patch closure for the ventricular septal defect (VSD; n=69) and atrial septal defect (ASD; n=42) was performed with left atrio-ventricular valve (LAVV) cleft closure (n=76) and right atrio-ventricular valve (RAVV) repair (n=57). Mortality, morbidity, and indications for reoperation were retrospectively studied; the end point 'time to reoperation' was analyzed using Kaplan-Meier curves. Follow-up was complete except in two patients and spanned a median of 28 months (range 0.4-6.1 years). RESULTS: In-hospital mortality was 3.7% (n=3) and one late death occurred. Reoperation was required in 7/79 patients (8.9%) for LAVV insufficiency (n=4), for a residual ASD (n=1), for right atrio-ventricular valve insufficiency (n=1), and for subaortic stenosis (n=1). At last follow-up, no or only mild LAVV and RAVV insufficiency was present in 81.3% and 92.1% of patients, respectively, and 2/3 of patients were medication-free. Risk factors for reoperation were younger age (<3 months; p=0.001) and lower weight (<4kg; p=0.003), and a trend towards less and later reoperations in Down syndrome (p<0.2). CONCLUSIONS: Surgical correction of AVSD can be achieved with low mortality and need for reoperation, regardless of Down syndrome or not. Immediate postoperative moderate or more residual atrio-ventricular valve insufficiency will eventually require a reoperation, and could be anticipated in patients younger than 3 months and weighing <4kg.

In utero exposure to immunosuppressive drugs.

The number of pregnant women receiving immunosuppressants for anti-rejection therapy or autoimmune diseases is increasing. All immunosuppressive drugs cross the placenta, raising questions about the long-term outcome of the children exposed in utero. There is no higher risk of congenital anomalies. However, an increased incidence of prematurity, intrauterine growth retardation (IUGR) and generally low birth weight has been reported, as well as maternal hypertension and preeclampsia. The most frequent neonatal complications are those associated with prematurity and IUGR, as well as adrenal insufficiency with corticosteroids, immunological disturbances with azathioprine and cyclosporine, and hyperkalemia with tacrolimus. The long-term follow-up of infants exposed to immunosuppressants in utero is still limited and experimental studies raise the question whether there could be an increased incidence at adult age of some pathologies including renal insufficiency, hypertension and diabetes.

A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics

To describe the collective behavior of large ensembles of neurons in neuronal network, a kinetic theory description was developed in [13, 12], where a macroscopic representation of the network dynamics was directly derived from the microscopic dynamics of individual neurons, which are modeled by conductance-based, linear, integrate-and-fire point neurons. A diffusion approximation then led to a nonlinear Fokker-Planck equation for the probability density function of neuronal membrane potentials and synaptic conductances. In this work, we propose a deterministic numerical scheme for a Fokker-Planck model of an excitatory-only network. Our numerical solver allows us to obtain the time evolution of probability distribution functions, and thus, the evolution of all possible macroscopic quantities that are given by suitable moments of the probability density function. We show that this deterministic scheme is capable of capturing the bistability of stationary states observed in Monte Carlo simulations. Moreover, the transient behavior of the firing rates computed from the Fokker-Planck equation is analyzed in this bistable situation, where a bifurcation scenario, of asynchronous convergence towards stationary states, periodic synchronous solutions or damped oscillatory convergence towards stationary states, can be uncovered by increasing the strength of the excitatory coupling. Finally, the computation of moments of the probability distribution allows us to validate the applicability of a moment closure assumption used in [13] to further simplify the kinetic theory.

A one-dimensional Keller-Segel equation with a drift issued from the boundary

We investigate in this note the dynamics of a one-dimensional Keller-Segel type model on the half-line. On the contrary to the classical configuration, the chemical production term is located on the boundary. We prove, under suitable assumptions, the following dichotomy which is reminiscent of the two-dimensional Keller-Segel system. Solutions are global if the mass is below the critical mass, they blow-up in finite time above the critical mass, and they converge to some equilibrium at the critical mass. Entropy techniques are presented which aim at providing quantitative convergence results for the subcritical case. This note is completed with a brief introduction to a more realistic model (still one-dimensional).

Immigrants'assimilation process in a segmented labor market

While much of the literature on immigrants' assimilation has focused on countries with a large tradition of receiving immigrants and with flexible labor markets, very little is known on how immigrants adjust to other types of host economies. With its severe dual labor market, and an unprecedented immigration boom, Spain presents a quite unique experience to analyze immigrations' assimilation process. Using data from the 2000 to 2008 Labor Force Survey, we find that immigrants are more occupationally mobile than natives, and that much of this greater flexibility is explained by immigrants' assimilation process soon after arrival. However, we find little evidence of convergence, especially among women and high skilled immigrants. This suggests that instead of integrating, immigrants occupationally segregate, providing evidence consistent with both imperfect substitutability and immigrants' human capital being under-valued. Additional evidence on the assimilation of earnings and the incidence of permanent employment by different skill levels also supports the hypothesis of segmented labor markets.

The NLRP3 inflammasome promotes renal inflammation and contributes to CKD.

Inflammation significantly contributes to the progression of chronic kidney disease (CKD). Inflammasome-dependent cytokines, such as IL-1β and IL-18, play a role in CKD, but their regulation during renal injury is unknown. Here, we analyzed the processing of caspase-1, IL-1β, and IL-18 after unilateral ureteral obstruction (UUO) in mice, which suggested activation of the Nlrp3 inflammasome during renal injury. Compared with wild-type mice, Nlrp3(-/-) mice had less tubular injury, inflammation, and fibrosis after UUO, associated with a reduction in caspase-1 activation and maturation of IL-1β and IL-18; these data confirm that the Nlrp3 inflammasome upregulates these cytokines in the kidney during injury. Bone marrow chimeras revealed that Nlrp3 mediates the injurious/inflammatory processes in both hematopoietic and nonhematopoietic cellular compartments. In tissue from human renal biopsies, a wide variety of nondiabetic kidney diseases exhibited increased expression of NLRP3 mRNA, which correlated with renal function. Taken together, these results strongly support a role for NLRP3 in renal injury and identify the inflammasome as a possible therapeutic target in the treatment of patients with progressive CKD.