Convergence of the mass-transport steepest descent scheme for the sub-critical Patlak-Keller-Segel model

Variational steepest descent approximation schemes for the modified Patlak-Keller-Segel equation with a logarithmic interaction kernel in any dimension are considered. We prove the convergence of the suitably interpolated in time implicit Euler scheme, defined in terms of the Euclidean Wasserstein distance, associated to this equation for sub-critical masses. As a consequence, we recover the recent result about the global in time existence of weak-solutions to the modified Patlak-Keller-Segel equation for the logarithmic interaction kernel in any dimension in the sub-critical case. Moreover, we show how this method performs numerically in one dimension. In this particular case, this numerical scheme corresponds to a standard implicit Euler method for the pseudo-inverse of the cumulative distribution function. We demonstrate its capabilities to reproduce easily without the need of mesh-refinement the blow-up of solutions for super-critical masses.

Diseño de filtros con topología apilada mediante resonadores FBAR

Els continguts de la memòria es divideixen en dues parts fonamentals, precedides d’una breu explicació on s’introdueix al lector en el tema; primer s’exposen les bases d’un ressonador simple FBAR, d’on s’obtenen equacions de control vitals, i seguidament es plantegen les bases per a un ressonador apilat SCR, derivades de l’estudi previ d’un ressonador simple. La primera part fonamental del treball es centra en l’anàlisi d’un ressonador SCR. Aquest anàlisi es recolza sobre el punt de vista teòric dels paràmetres imatge, el nostre punt de sortida. Aquesta primera part és la més teòrica: obtenció i aplicació dels paràmetres imatge i obtenció dels elements discrets que conformen el circuit equivalent SCR, una xarxa de dos ports. Posteriorment, s’analitza què succeeix modificant els valors dels elements discrets sense variar determinats paràmetres imatge i finalment, es proposa una aproximació per a controlar determinades especificacions de disseny, com són l’ample de banda de transmissió de la xarxa i el factor de qualitat, en funció de les modificacions dels elements discrets. En la segona part s’analitzen, de forma qualitativa, xarxes compostes per diferents ressonadors apilats connectats en cascada. Aquest segon estudi es divideix en dues subparts. En la primera connectem N ressonadors idèntics, plantegem algunes equacions de control i analitzem les respostes. En la segona, es planteja la connexió de dos (N=2) ressonadors diferents amb freqüències de ressonància properes. Aquest segon anàlisi es també totalment qualitatiu, però ens aporta informació que amb la unió de N ressonadors iguals no aconseguíem. Finalitzant aquesta segona part, es planteja l’optimització dels resultats obtinguts per a N=2, mitjançant estructures N=3 ressonadors.

Ordering Policy Rules with an Unconditional Welfare Measure

The unconditional expectation of social welfare is often used to assess alternative macroeconomic policy rules in applied quantitative research. It is shown that it is generally possible to derive a linear - quadratic problem that approximates the exact non-linear problem where the unconditional expectation of the objective is maximised and the steady-state is distorted. Thus, the measure of pol icy performance is a linear combinat ion of second moments of economic variables which is relatively easy to compute numerically, and can be used to rank alternative policy rules. The approach is applied to a simple Calvo-type model under various monetary policy rules.

Tractable Consumer Choice

We derive a rational model of separable consumer choice which can also serve as a behavioral model. The central construct is [lambda] , the marginal utility of money, derived from the consumer's rest-of-life problem. We present a robust approximation of [lambda], and show how to incorporate liquidity constraints, indivisibilities and adaptation to a changing environment. We fi nd connections with numerous historical and recent constructs, both behavioral and neoclassical, and draw contrasts with standard partial equilibrium analysis. The result is a better grounded, more flexible and more intuitive description of consumer choice.

The Inflation Bias under Calvo and Rotemberg Pricing.

New Keynesian models rely heavily on two workhorse models of nominal inertia - price contracts of random duration (Calvo, 1983) and price adjustment costs (Rotemberg, 1982) - to generate a meaningful role for monetary policy. These alternative descriptions of price stickiness are often used interchangeably since, to a first order of approximation they imply an isomorphic Phillips curve and, if the steady-state is efficient, identical objectives for the policy maker and as a result in an LQ framework, the same policy conclusions. In this paper we compute time-consistent optimal monetary policy in bench-mark New Keynesian models containing each form of price stickiness. Using global solution techniques we find that the inflation bias problem under Calvo contracts is significantly greater than under Rotemberg pricing, despite the fact that the former typically significant exhibits far greater welfare costs of inflation. The rates of inflation observed under this policy are non-trivial and suggest that the model can comfortably generate the rates of inflation at which the problematic issues highlighted in the trend inflation literature emerge, as well as the movements in trend inflation emphasized in empirical studies of the evolution of inflation. Finally, we consider the response to cost push shocks across both models and find these can also be significantly different. The choice of which form of nominal inertia to adopt is not innocuous.

Bilateral oligopoly and quantity competition

Bilateral oligopoly is a strategic market game with two commodities, allowing strategic behavior on both sides of the market. When the number of buyers is large, such a game approximates a game of quantity competition played by sellers. We present examples which show that this is not typically a Cournot game. Rather, we introduce an alternative game of quantity competition (the market share game) and, appealing to results in the literature on contests, show that this yields the same equilibria as the many-buyer limit of bilateral oligopoly, under standard assumptions on costs and preferences. We also show that the market share and Cournot games have the same equilibria if and only if the price elasticity of the latter is one. These results lead to necessary and sufficient conditions for the Cournot game to be a good approximation to bilateral oligopoly with many buyers and to an ordering of total output when they are not satisfied.

The non-linear Dirichlet problem in Hadamard manifolds

We prove existence theorems for the Dirichlet problem for hypersurfaces of constant special Lagrangian curvature in Hadamard manifolds. The first results are obtained using the continuity method and approximation and then refined using two iterations of the Perron method. The a-priori estimates used in the continuity method are valid in any ambient manifold.

Quantification of the neurochemical profile using simulated macromolecule resonances at 3 T.

The broad resonances underlying the entire (1) H NMR spectrum of the brain, ascribed to macromolecules, can influence metabolite quantification. At the intermediate field strength of 3 T, distinct approaches for the determination of the macromolecule signal, previously used at either 1.5 or 7 T and higher, may become equivalent. The aim of this study was to evaluate, at 3 T for healthy subjects using LCModel, the impact on the metabolite quantification of two different macromolecule approaches: (i) experimentally measured macromolecules; and (ii) mathematically estimated macromolecules. Although small, but significant, differences in metabolite quantification (up to 23% for glutamate) were noted for some metabolites, 10 metabolites were quantified reproducibly with both approaches with a Cramer-Rao lower bound below 20%, and the neurochemical profiles were therefore similar. We conclude that the mathematical approximation can provide sufficiently accurate and reproducible estimation of the macromolecule contribution to the (1) H spectrum at 3 T. Copyright © 2013 John Wiley & Sons, Ltd.

Discontinuous Galerkin methods for the one-dimensional Vlasov-Poisson system

We construct a new family of semi-discrete numerical schemes for the approximation of the one-dimensional periodic Vlasov-Poisson system. The methods are based on the coupling of discontinuous Galerkin approximation to the Vlasov equation and several finite element (conforming, non-conforming and mixed) approximations for the Poisson problem. We show optimal error estimates for the all proposed methods in the case of smooth compactly supported initial data. The issue of energy conservation is also analyzed for some of the methods.

Application of standard and refined heat balance integral methods to one-dimensional Stefan problems

The work in this paper concerns the study of conventional and refined heat balance integral methods for a number of phase change problems. These include standard test problems, both with one and two phase changes, which have exact solutions to enable us to test the accuracy of the approximate solutions. We also consider situations where no analytical solution is available and compare these to numerical solutions. It is popular to use a quadratic profile as an approximation of the temperature, but we show that a cubic profile, seldom considered in the literature, is far more accurate in most circumstances. In addition, the refined integral method can give greater improvement still and we develop a variation on this method which turns out to be optimal in some cases. We assess which integral method is better for various problems, showing that it is largely dependent on the specified boundary conditions.

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.

Generic Nekhoroshev theory without small divisors

In this article, we present a new approach of Nekhoroshev theory for a generic unperturbed Hamiltonian which completely avoids small divisors problems. The proof is an extension of a method introduced by P. Lochak which combines averaging along periodic orbits with simultaneous Diophantine approximation and uses geometric arguments designed by the second author to handle generic integrable Hamiltonians. This method allows to deal with generic non-analytic Hamiltonians and to obtain new results of generic stability around linearly stable tori.

Metropolitan Areas in Spain and Italy

Metropolitan areas concentrate the main share of population, production and consumption in OECD countries. They are likely to be the most important units for economic, social and environmental analysis as well as for the development of policy strategies. However, one of the main problems that occur when adopting metropolitan areas as units of analysis and policy in European countries is the absence of widely accepted standards for identifying them. This severe problem appeared when we tried to perform comparative research between Spain and Italy using metropolitan areas as units of analysis. The aim of this paper is to identify metropolitan areas in Spain and Italy using similar methodologies. The results allow comparing the metropolitan realities of both countries as well as providing the metropolitan units that can be used in subsequent comparative researches. Two methodologies are proposed: the Cheshire-GEMACA methodology (FUR) and an iterative version of the USA-MSA algorithm, particularly adapted to deal with polycentric metropolitan areas (DMA). Both methods show a good approximation to the metropolitan reality and produce very similar results: 75 FUR and 67 DMA in Spain (75% of total population and employment), and 81 FUR and 86 DMA in Italy (70% of total population and employment).

The brazilian leishmaniasis control program

The control methods used in Brazil for both cutaneous and visceral leishmaniasis by the Ministry of Health's National Health Foundation (Fundação Nacional de Saúde) are presented. Data referring to the programme between 1980 and 1991 is presented. There was an increase in the total number of notified cases of cutaneous leishmaniasis until 1987 and from this date onwards the number has remained between 23,000 and 26,000 per year. There were significantly fewer cases of visceral leishmaniasis during the same period and until 1991 the greatest number of cases reported in one year was 2,511 in 1985. The overall picture is that it is difficult to evaluate to what extent the control methods employed have reduced the number of cases of either form of the disease. It is suggested that there is an urgent need for a closer approximation of research and health workers to review the actual control strategies and to define procedures capable of assessing accurately their impact.

On estimation and identifiability issues of sex-linked inheritance with a case study of pigmentation in Swiss barn owl (Tyto alba)

Genetic evaluation using animal models or pedigree-based models generally assume only autosomal inheritance. Bayesian animal models provide a flexible framework for genetic evaluation, and we show how the model readily can accommodate situations where the trait of interest is influenced by both autosomal and sex-linked inheritance. This allows for simultaneous calculation of autosomal and sex-chromosomal additive genetic effects. Inferences were performed using integrated nested Laplace approximations (INLA), a nonsampling-based Bayesian inference methodology. We provide a detailed description of how to calculate the inverse of the X- or Z-chromosomal additive genetic relationship matrix, needed for inference. The case study of eumelanic spot diameter in a Swiss barn owl (Tyto alba) population shows that this trait is substantially influenced by variation in genes on the Z-chromosome (sigma(2)(z) = 0.2719 and sigma(2)(a) = 0.4405). Further, a simulation study for this study system shows that the animal model accounting for both autosomal and sex-chromosome-linked inheritance is identifiable, that is, the two effects can be distinguished, and provides accurate inference on the variance components.