CRITICAL BEHAVIOR AND THRESHOLD OF COEXISTENCE OF A PREDATOR-PREY STOCHASTIC MODEL IN A 2D LATTICE

We investigate the critical behavior of a stochastic lattice model describing a predator-prey system. By means of Monte Carlo procedure we simulate the model defined on a regular square lattice and determine the threshold of species coexistence, that is, the critical phase boundaries related to the transition between an active state, where both species coexist and an absorbing state where one of the species is extinct. A finite size scaling analysis is employed to determine the order parameter, order parameter fluctuations, correlation length and the critical exponents. Our numerical results for the critical exponents agree with those of the directed percolation universality class. We also check the validity of the hyperscaling relation and present the data collapse curves.

Ground-state energy of a classical artificial molecule

We study the ground-state energy of a classical artificial molecule formed by two-dimensional clusters (artificial atoms) of N/2 charged particles separated by a distance d. For the small molecules of N = 2 and 4, we obtain analytical expressions for this energy. For the larger ones, we calculate the ground-state energy using molecular dynamics simulation for N up to 128. From our numerical results, we are able to find out a function to approximate the ground-state energy of the molecules covering the range from atoms to molecules for any inter-atom distance d and for particle number from N = 8 to 128 within a difference less than one percent from the MD data.

On the generalized eigenvalue method for energies and matrix elements in lattice field theory

We discuss the generalized eigenvalue problem for computing energies and matrix elements in lattice gauge theory, including effective theories such as HQET. It is analyzed how the extracted effective energies and matrix elements converge when the time separations are made large. This suggests a particularly efficient application of the method for which we can prove that corrections vanish asymptotically as exp(-(E(N+1) - E(n))t). The gap E(N+1) - E(n) can be made large by increasing the number N of interpolating fields in the correlation matrix. We also show how excited state matrix elements can be extracted such that contaminations from all other states disappear exponentially in time. As a demonstration we present numerical results for the extraction of ground state and excited B-meson masses and decay constants in static approximation and to order 1/m(b) in HQET.

Orthogonal packing of identical rectangles within isotropic convex regions

A mixed integer continuous nonlinear model and a solution method for the problem of orthogonally packing identical rectangles within an arbitrary convex region are introduced in the present work. The convex region is assumed to be made of an isotropic material in such a way that arbitrary rotations of the items, preserving the orthogonality constraint, are allowed. The solution method is based on a combination of branch and bound and active-set strategies for bound-constrained minimization of smooth functions. Numerical results show the reliability of the presented approach. (C) 2010 Elsevier Ltd. All rights reserved.

Improved maximum likelihood estimators in a heteroskedastic errors-in-variables model

This paper develops a bias correction scheme for a multivariate heteroskedastic errors-in-variables model. The applicability of this model is justified in areas such as astrophysics, epidemiology and analytical chemistry, where the variables are subject to measurement errors and the variances vary with the observations. We conduct Monte Carlo simulations to investigate the performance of the corrected estimators. The numerical results show that the bias correction scheme yields nearly unbiased estimates. We also give an application to a real data set.

Improved likelihood inference in Birnbaum-Saunders regressions

The Birnbaum-Saunders regression model is commonly used in reliability studies. We address the issue of performing inference in this class of models when the number of observations is small. Our simulation results suggest that the likelihood ratio test tends to be liberal when the sample size is small. We obtain a correction factor which reduces the size distortion of the test. Also, we consider a parametric bootstrap scheme to obtain improved critical values and improved p-values for the likelihood ratio test. The numerical results show that the modified tests are more reliable in finite samples than the usual likelihood ratio test. We also present an empirical application. (C) 2009 Elsevier B.V. All rights reserved.

Testing hypotheses in the Birnbaum-Saunders distribution under type-II censored samples

The two-parameter Birnbaum-Saunders distribution has been used successfully to model fatigue failure times. Although censoring is typical in reliability and survival studies, little work has been published on the analysis of censored data for this distribution. In this paper, we address the issue of performing testing inference on the two parameters of the Birnbaum-Saunders distribution under type-II right censored samples. The likelihood ratio statistic and a recently proposed statistic, the gradient statistic, provide a convenient framework for statistical inference in such a case, since they do not require to obtain, estimate or invert an information matrix, which is an advantage in problems involving censored data. An extensive Monte Carlo simulation study is carried out in order to investigate and compare the finite sample performance of the likelihood ratio and the gradient tests. Our numerical results show evidence that the gradient test should be preferred. Further, we also consider the generalized Birnbaum-Saunders distribution under type-II right censored samples and present some Monte Carlo simulations for testing the parameters in this class of models using the likelihood ratio and gradient tests. Three empirical applications are presented. (C) 2011 Elsevier B.V. All rights reserved.

Improved maximum-likelihood estimation in a regression model with general parametrization

We analyse the finite-sample behaviour of two second-order bias-corrected alternatives to the maximum-likelihood estimator of the parameters in a multivariate normal regression model with general parametrization proposed by Patriota and Lemonte [A. G. Patriota and A. J. Lemonte, Bias correction in a multivariate regression model with genereal parameterization, Stat. Prob. Lett. 79 (2009), pp. 1655-1662]. The two finite-sample corrections we consider are the conventional second-order bias-corrected estimator and the bootstrap bias correction. We present the numerical results comparing the performance of these estimators. Our results reveal that analytical bias correction outperforms numerical bias corrections obtained from bootstrapping schemes.

A note on an application of Arrow's theorem: sufficient conditions for Lucas' inflation and welfare

Bellman's methods for dynamic optimization constitute the present mainstream in economics. However, some results associated with optimal controI can be particularly usefuI in certain problems. The purpose of this note is presenting such an example. The value function derived in Lucas' (2000) shopping-time economy in Infiation and Welfare need not be concave, leading this author to develop numerical analyses to determine if consumer utility is in fact maximized along the balanced path constructed from the first order conditions. We use Arrow's generalization of Mangasarian's results in optimal control theory and develop sufficient conditions for the problem. The analytical conclusions and the previous numerical results are compatible .

Modelos cinéticos da equação linearizada de Boltzmann e um problema de transferência de calor em microescala

Neste trabalho, um problema de transferência de calor da dinâmica de gases rarefeitos, causado pela diferença de temperaturas nas superfícies de um canal, é abordado. O problema é formulado através dos modelos cinéticos BGK, S e Gross-Jackson da equação linearizada de Boltzmann e resolvido, de forma unificada, pelo método analítico de ordenadas discretas (método ADO). Resultados numéricos para as perturbações de densidade e temperatura e também para o fluxo de calor são apresentados e comparados, mostrando que não se pode dizer que algum dos três modelos seja uma melhor aproximação da solução aos resultados da equação linearizada de Boltzmann.

Derivativos de volatilidade no mercado brasileiro de câmbio : viabilidade e impactos de sua utilização

A volatilidade possui um papel central na gestão de risco tanto de portfólios de derivativos como de portfólios de ativos não alavancados. Este risco é gerenciado nos mercados financeiros através de diversos instrumentos, incluindo o uso de derivativos de volatilidade. No entanto, um mercado de derivativos de volatilidade no Brasil ainda é uma lacuna a ser preenchida, talvez pela baixa liquidez em determinadas opções ou mesmo pela falta de todos os ativos necessários para se estabelecer o portfólio replicante para os mesmos. O objetivo deste trabalho é apresentar um modelo simples de se apreçar swaps de volatilidade sob o BRL, estimulando um diálogo entre a comunidade acadêmica e os praticantes do mercado que permita o desenvolvimento de derivativos de volatilidade ao considerar o melhor de cada grupo. Para se apreçar este instrumento, a modelagem e os ativos utilizados são apresentados em detalhes como sendo os ingredientes básicos de um produto financeiro de sucesso. Os resultados numéricos demonstram que o modelo proposto pode ser considerado um poderoso instrumento para se realizar o hedge do risco de volatilidade. Um benefício adicional deste trabalho é apresentar os riscos e benefícios de se utilizar este instrumento com o BRL.

Share bidding auctions, sliding scale royalty rates and the new brazilian regulatory framework for the pre-salt areas

O governo brasileiro recentemente aprovou uma legislação instituindo um novo marco regulatório para as reservas petrolíferas do pré-sal. Segundo as novas regras, estas áreas deverão ser licitadas mediante um leilão de partilha de lucro. Motivado por esta mudança, apresentamos um modelo de leilão de partilha sob afiliação, demonstrando a existência de um equilíbrio monótono em estratégias puras e caracterizando a solução. Alem disso, provamos que este mecanismo gera receita esperada maior ou igual a um leilão de primeiro preço usual. Em seguida, introduzimos no modelo uma função representando taxas de royalties que dependem do valor do objeto. Este instrumento permite uma elevação na receita esperada de ambos os modelos, fazendo com que a diferença entre eles encolha. Finalmente, analisando o novo marco regulatório sob o ponto de vista dos resultados obtidos, concluímos que o antigo modelo de concessão utilizado pelo governo brasileiro é mais adequado e lucrativo.

Legal enforcement, collateral and heterogeneity of project financing contracts

This paper employs mechanism design to study the effects of imperfect legal enforcement on optimal scale of projects, borrowing interest rates and the probability of default. The analysis departs from an environment that combines asymmetric information about cash flows and limited commitment by borrowers. Incentive for repayment comes from the possibility of liquidation of projects by a court, but courts are costly and may fail to liquidate. The value of liquidated assets can be used as collateral: it is transferred to the lender when courts liquidate. Examples reveal that costly use of courts may be optimal, which contrasts with results from most limited commitment models, where punishments are just threats, never applied in optimal arrangements. I show that when voluntary liquidation is allowed, both asymmetric information and uncertainty about courts are necessary conditions for legal punishments ever to be applied. Numerical solutions for several parametric specifications are presented, allowing for heterogeneity on initial wealth and variability of project returns. In all such solutions, wealthier individuals borrow with lower interest rates and run higher scale enterprises, which is consistent with stylized facts. The reliability of courts has a consistently positive effect on the scale of projects. However its effect on interest rates is subtler and depends essentially on the degree of curvature of the production function. Numerical results also show that the possibility of collateral seizing allows comovements of the interest rates and the probability of repayment.

Um modelo matemático baseado em wavelets para análise do método térmico de recuperação de óleo pesado aplicando irradiação eletromagnética

This work proposes a model to investigate the use of a cylindrical antenna used in the thermal method of recovering through electromagnetic radiation of high-viscosity oil. The antenna has a simple geometry, adapted dipole type, and it can be modelled by using Maxwell s equation. The wavelet transforms are used as basis functions and applied in conjunction with the method of moments to obtain the current distribution in the antenna. The electric field, power and temperature distribution are carefully calculated for the analysis of the antenna as electromagnetic heating. The energy performance is analyzed based on thermo-fluid dynamic simulations at field scale, and through the adaptation in the Steam Thermal and Advanced Processes Reservoir Simulator (STARS) by Computer Modelling Group (CMG). The model proposed and the numerical results obtained are stable and presented good agreement with the results reported in the specialized literature

Simulação de fluxo de fluidos em meios porosos desordenados uma análise de efeito de escala na estimativa da permeabilidade e do coeficiente de arrasto

The present study provides a methodology that gives a predictive character the computer simulations based on detailed models of the geometry of a porous medium. We using the software FLUENT to investigate the flow of a viscous Newtonian fluid through a random fractal medium which simplifies a two-dimensional disordered porous medium representing a petroleum reservoir. This fractal model is formed by obstacles of various sizes, whose size distribution function follows a power law where exponent is defined as the fractal dimension of fractionation Dff of the model characterizing the process of fragmentation these obstacles. They are randomly disposed in a rectangular channel. The modeling process incorporates modern concepts, scaling laws, to analyze the influence of heterogeneity found in the fields of the porosity and of the permeability in such a way as to characterize the medium in terms of their fractal properties. This procedure allows numerically analyze the measurements of permeability k and the drag coefficient Cd proposed relationships, like power law, for these properties on various modeling schemes. The purpose of this research is to study the variability provided by these heterogeneities where the velocity field and other details of viscous fluid dynamics are obtained by solving numerically the continuity and Navier-Stokes equations at pore level and observe how the fractal dimension of fractionation of the model can affect their hydrodynamic properties. This study were considered two classes of models, models with constant porosity, MPC, and models with varying porosity, MPV. The results have allowed us to find numerical relationship between the permeability, drag coefficient and the fractal dimension of fractionation of the medium. Based on these numerical results we have proposed scaling relations and algebraic expressions involving the relevant parameters of the phenomenon. In this study analytical equations were determined for Dff depending on the geometrical parameters of the models. We also found a relation between the permeability and the drag coefficient which is inversely proportional to one another. As for the difference in behavior it is most striking in the classes of models MPV. That is, the fact that the porosity vary in these models is an additional factor that plays a significant role in flow analysis. Finally, the results proved satisfactory and consistent, which demonstrates the effectiveness of the referred methodology for all applications analyzed in this study.