Semidiscretization and long-time asymptotics of nonlinear diffusion equations

We review several results concerning the long time asymptotics of nonlinear diffusion models based on entropy and mass transport methods. Semidiscretization of these nonlinear diffusion models are proposed and their numerical properties analysed. We demonstrate the long time asymptotic results by numerical simulation and we discuss several open problems based on these numerical results. We show that for general nonlinear diffusion equations the long-time asymptotics can be characterized in terms of fixed points of certain maps which are contractions for the euclidean Wasserstein distance. In fact, we propose a new scaling for which we can prove that this family of fixed points converges to the Barenblatt solution for perturbations of homogeneous nonlinearities for values close to zero.

Generalized backward doubly stochastic differential equations and SPDEs with nonlinear Neumann boundary conditions

In this paper, a new class of generalized backward doubly stochastic differential equations is investigated. This class involves an integral with respect to an adapted continuous increasing process. A probabilistic representation for viscosity solutions of semi-linear stochastic partial differential equations with a Neumann boundary condition is given.

Determinació de la ubicació estructural dels ions Yb3+ i Nb3+ en Yb:RTP, Nb:RTP i (Nb,Yb):RTP

Estudi elaborat a partir d’una estada al Stony Brook University al juliol del 2006. El RbTiOPO4 (RTP) monocristal•lí és un material d' òptica no lineal molt rellevant i utilitzat en la tecnologia làser actual, químicament molt estable i amb unes propietats físiques molt destacades, entre elles destaquen els alts coeficients electro-òptics i l'alt llindar de dany òptic que presenta. En els últims anys s’està utilitzant tecnològicament en aplicacions d'òptica no lineal en general i electro-òptiques en particular. En alguns casos ja ha substituït, millorant prestacions, a materials tals com el KTP o el LNB(1). Dopant RTP amb ions lantànids (Ln3+) (2-4), el material es converteix en un material làser auto-doblador de freqüència, combinant les seves propietats no lineals amb les de matriu làser. El RTP genera radiació de segon harmònic (SHG) a partir d’un feix fonamental amb longituds d’ona inferiors a 990 nm, que és el límit que presenta el KTP.La determinació de la ubicació estructural i l’estudi de l'entorn local del ions actius làser és de fonamental importància per a la correcta interpretació de les propietats espectroscòpiques d’aquest material. Mesures de difracció de neutrons sobre mostra de pols cristal•lí mostren que els ions Nb5+ i Ln3+ només substitueixin posicions de Ti4+ (8-9). Estudis molt recents d'EPR (electron paramagnetic resonance) semblen indicar que quan la concentració d'ió Ln3+ es baixa, aquest ió presenta la tendència a substituir l'ió alcalí present a l'estructura (10).Després dels resultats obtinguts en el present treball a partir de la tècnica EXAFS a la instal•lació sincrotò del Brookhaven National Laboratory/State University of New York (Stony Brook) es pot concloure definitivament que els ions Nb s’ubiquen en la posició Ti (1) i que els ions Yb3+ es distribueixen paritariament en les dues posicions del Ti (1 i 2). Aquests resultats aporten una valuosa informació per a la correcta interpretació dels espectres, tant d’absorció com d’emissió, del material i per la avaluació dels paràmetres del seu comportament durant l'acció làser.

A nonlinear threshold model for the dependence of extremes of stationary sequences

One of the main implications of the efficient market hypothesis (EMH) is that expected future returns on financial assets are not predictable if investors are risk neutral. In this paper we argue that financial time series offer more information than that this hypothesis seems to supply. In particular we postulate that runs of very large returns can be predictable for small time periods. In order to prove this we propose a TAR(3,1)-GARCH(1,1) model that is able to describe two different types of extreme events: a first type generated by large uncertainty regimes where runs of extremes are not predictable and a second type where extremes come from isolated dread/joy events. This model is new in the literature in nonlinear processes. Its novelty resides on two features of the model that make it different from previous TAR methodologies. The regimes are motivated by the occurrence of extreme values and the threshold variable is defined by the shock affecting the process in the preceding period. In this way this model is able to uncover dependence and clustering of extremes in high as well as in low volatility periods. This model is tested with data from General Motors stocks prices corresponding to two crises that had a substantial impact in financial markets worldwide; the Black Monday of October 1987 and September 11th, 2001. By analyzing the periods around these crises we find evidence of statistical significance of our model and thereby of predictability of extremes for September 11th but not for Black Monday. These findings support the hypotheses of a big negative event producing runs of negative returns in the first case, and of the burst of a worldwide stock market bubble in the second example. JEL classification: C12; C15; C22; C51 Keywords and Phrases: asymmetries, crises, extreme values, hypothesis testing, leverage effect, nonlinearities, threshold models

Is the consumption-income ratio stationary? Evidence from linear and nonlinear panel unit root tests for OECD and non-OECD countries

This paper applies recently developed heterogeneous nonlinear and linear panel unit root tests that account for cross-sectional dependence to 24 OECD and 33 non-OECD countries’ consumption-income ratios over the period 1951–2003. We apply a recently developed methodology that facilitates the use of panel tests to identify which individual cross-sectional units are stationary and which are nonstationary. This extends evidence provided in the recent literature to consider both linear and nonlinear adjustment in panel unit root tests, to address the issue of cross-sectional dependence, and to substantially expand both time-series and cross sectional dimensions of the data analysed. We find that the majority (65%) of the series are nonstationary with slightly fewer OECD countries’ (61%) series exhibiting a unit root than non-OECD countries (68%).

Microstructure Order Flow: Statistical and Economic Evaluation of Nonlinear Forecasts

In this paper we propose a novel empirical extension of the standard market microstructure order flow model. The main idea is that heterogeneity of beliefs in the foreign exchange market can cause model instability and such instability has not been fully accounted for in the existing empirical literature. We investigate this issue using two di¤erent data sets and focusing on out- of-sample forecasts. Forecasting power is measured using standard statistical tests and, additionally, using an alternative approach based on measuring the economic value of forecasts after building a portfolio of assets. We nd there is a substantial economic value on conditioning on the proposed models.

A Nonlinear Panel Unit Root Test under Cross Section Dependence

We propose a nonlinear heterogeneous panel unit root test for testing the null hypothesis of unit-roots processes against the alternative that allows a proportion of units to be generated by globally stationary ESTAR processes and a remaining non-zero proportion to be generated by unit root processes. The proposed test is simple to implement and accommodates cross sectional dependence. We show that the distribution of the test statistic is free of nuisance parameters as (N, T) −! 1. Monte Carlo simulation shows that our test holds correct size and under the hypothesis that data are generated by globally stationary ESTAR processes has a better power than the recent test proposed in Pesaran [2007]. Various applications are provided.

Positive solutions of nonlinear problems involving the square root of the Laplacian

We consider nonlinear elliptic problems involving a nonlocal operator: the square root of the Laplacian in a bounded domain with zero Dirichlet boundary conditions. For positive solutions to problems with power nonlinearities, we establish existence and regularity results, as well as a priori estimates of Gidas-Spruck type. In addition, among other results, we prove a symmetry theorem of Gidas-Ni-Nirenberg type.

Normal kinetics of intestinal glucose absorption in the absence of GLUT2: evidence for a transport pathway requiring glucose phosphorylation and transfer into the endoplasmic reticulum.

Glucose is absorbed through the intestine by a transepithelial transport system initiated at the apical membrane by the cotransporter SGLT-1; intracellular glucose is then assumed to diffuse across the basolateral membrane through GLUT2. Here, we evaluated the impact of GLUT2 gene inactivation on this transepithelial transport process. We report that the kinetics of transepithelial glucose transport, as assessed in oral glucose tolerance tests, was identical in the presence or absence of GLUT2; that the transport was transcellular because it could be inhibited by the SGLT-1 inhibitor phlorizin, and that it could not be explained by overexpression of another known glucose transporter. By using an isolated intestine perfusion system, we demonstrated that the rate of transepithelial transport was similar in control and GLUT2(-/-) intestine and that it was increased to the same extent by cAMP in both situations. However, in the absence, but not in the presence, of GLUT2, the transport was inhibited dose-dependently by the glucose-6-phosphate translocase inhibitor S4048. Furthermore, whereas transport of [(14)C]glucose proceeded with the same kinetics in control and GLUT2(-/-) intestine, [(14)C]3-O-methylglucose was transported in intestine of control but not of mutant mice. Together our data demonstrate the existence of a transepithelial glucose transport system in GLUT2(-/-) intestine that requires glucose phosphorylation and transfer of glucose-6-phosphate into the endoplasmic reticulum. Glucose may then be released out of the cells by a membrane traffic-based pathway similar to the one we previously described in GLUT2-null hepatocytes.

Rigorous derivation of a nonlinear diffusion equation as fast-reaction limit of a continuous coagulation-fragmentation model with diffusion

Weak solutions of the spatially inhomogeneous (diffusive) Aizenmann-Bak model of coagulation-breakup within a bounded domain with homogeneous Neumann boundary conditions are shown to converge, in the fast reaction limit, towards local equilibria determined by their mass. Moreover, this mass is the solution of a nonlinear diffusion equation whose nonlinearity depends on the (size-dependent) diffusion coefficient. Initial data are assumed to have integrable zero order moment and square integrable first order moment in size, and finite entropy. In contrast to our previous result [CDF2], we are able to show the convergence without assuming uniform bounds from above and below on the number density of clusters.

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.

Effect of somatostatin on duodenal glucose absorption in man.

OBJECTIVE: The hyperglycemic hyperinsulinemic clamp technique using intraduodenally infused glucose is an attractive tool for studying postprandial glucose metabolism under strictly controlled conditions. Because it requires the use of somatostatin (SST), we examined, in this study, the effect of SST on intestinal glucose absorption. CONTEXT: Twenty-six normal volunteers were given a constant 3-h intraduodenal infusion of glucose (6 mg.kg(-1).min(-1)) labeled with [2-(3)H]glucose for glucose absorption measurement. During glucose infusion, 19 subjects received iv SST at doses of 10-100 ng.kg(-1).min(-1) plus insulin and glucagon, and seven subjects were studied under control conditions. In the controls, glucose was absorbed at a rate that, after a 20-min lag period, equaled the infusion rate. RESULTS: With all the doses of SST tested, absorption was considerably delayed but equaled the rate of infusion after 3 h. At that time, only 5 +/- 2% of the total amount of infused glucose was unabsorbed in the control subjects vs. 36 +/- 2% (P < 0.001) in the SST-infused subjects. In the latter, the intraluminal residue was almost totally absorbed within 40 min of the cessation of SST infusion. At the lowest dose of SST tested (10 ng.kg(-1).min(-1)), suppression of insulin secretion was incomplete. CONCLUSION: These properties of SST hamper the use of intraduodenal hyperglycemic hyperinsulinemic clamps as a tool for exploring postprandial glucose metabolism.

Intestinal protein absorption in malnourished mice with acute schistosomiasis mansoni

Intestinal protein absorption was studied in undernourished albino Swiss mice with acute schistosomiasis mansoni. Undernutrition was induced by feeding mice with the Regional Basic Diet (RBD) ingested by human populations in Northeast Brazil, an experimental model previously developed in our laboratory. Weaning mice were infected with 40 cercariae and compared to undernourished non-infected mice and/or to infected mice fed a balanced control diet. Apparent and True Protein Absorption Coefficients were determined by nitrogen balance during five consecutive days ending at the 63rd day of the trial (acute phase of murine schistosomiasis). Fecal metabolic nitrogen (FMN) was determined after administration of a non-protein diet and was also calculated through linear regression. Our results showed a reduced protein absorption in non-infected RBD-fed mice as compared to mice fed a casein control diet. Infection with Schistosoma mansoni had apparently no effect on intestinal protein absorption in well-nourished mice. However, infection seemed to interfere with protein absorption in under-nourished animals, since the lowest absorption ratios have been detected among RBD-fed infected mice. A brief discussion is made on the advantages of using the method of linear regression for the determination of FMN.

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.