Robustness of entanglement

In the quest to completely describe entanglement in the general case of a finite number of parties sharing a physical system of finite-dimensional Hilbert space an entanglement magnitude is introduced for its pure and mixed states: robustness. It corresponds to the minimal amount of mixing with locally prepared states which washes out all entanglement. It quantifies in a sense the endurance of entanglement against noise and jamming. Its properties are studied comprehensively. Analytical expressions for the robustness are given for pure states of two-party systems, and analytical bounds for mixed states of two-party systems. Specific results are obtained mainly for the qubit-qubit system (qubit denotes quantum bit). As by-products local pseudomixtures are generalized, a lower bound for the relative volume of separable states is deduced, and arguments for considering convexity a necessary condition of any entanglement measure are put forward.

Collineations of a symmetric 2-covariant tensor: Ricci collineations

The infinitesimal transformations that leave invariant a two-covariant symmetric tensor are studied. The interest of these symmetry transformations lays in the fact that this class of tensors includes the energy-momentum and Ricci tensors. We find that in most cases the class of infinitesimal generators of these transformations is a finite dimensional Lie algebra, but in some cases exhibiting a higher degree of degeneracy, this class is infinite dimensional and may fail to be a Lie algebra. As an application, we study the Ricci collineations of a type B warped spacetime.

Notas del primer seminario de integrabilidad de la Universitat Politècnica de Catalunya

Estas notas corresponden a las exposiciones presentadas en el \emph{Primer Seminario de Integrabilidad}, dentro de lo que se denomina \emph{Aula de Sistemas Din\'amicos}. Durante este evento se realizaron seis conferencias, todas presentadas por miembros del grupo de Sistemas Din\'amicos de la UPC. El programa desarrollado fue el siguiente:\\\begin{center}AULA DE SISTEMAS DIN\'AMICOS\end{center}\begin{center}\texttt{http://www.ma1.upc.es/recerca/seminaris/aulasd-cat.html}\end{center}\begin{center}SEMINARIO DE INTEGRABILIDAD\end{center}\begin{center}Martes 29 y Mi\'ercoles 30 de marzo de 2005\\Facultad de Matem\'aticas y Estad\'{\i}stica, UPC\\Aula: Seminario 1\end{center}\bigskip\begin{center}PROGRAMA Y RES\'UMENES\end{center}{\bf Martes 29 de marzo}\begin{itemize}\item15:30. Juan J. Morales-Ruiz. \emph{El problema de laintegrabilidad en Sistemas Din\'amicos}\medskip {\bf Resumen.} En esta presentaci\'on se pretende dar unaidea de conjunto, pero sin entrar en detalles, sobre las diversasnociones de integrabilidad, asociadas a nombres de matem\'aticostan ilustres como Liouville, Galois-Picard-Vessiot, Lie, Darboux,Kowalevskaya, Painlev\'e, Poincar\'e, Kolchin, Lax, etc. Adem\'astambi\'en mencionaremos la revoluci\'on que supuso en los a\~nossesenta del siglo pasado el descubrimiento de Gardner, Green,Kruskal y Miura sobre un nuevo m\'etodo para resolver en algunoscasos determinadas ecuaciones en derivadas parciales. \medskip\item16:00. David G\'omez-Ullate. \emph{Superintegrabilidad, pares deLax y modelos de $N-$cuerpos en el plano}\medskip{\bf Resumen.} Introduciremos algunas t\'ecnicas cl\'asicas paraconstruir modelos de N-cuerpos integrables, como los pares de Laxo la din\'amica de los ceros de un polinomio. Revisaremos lanoci\'on de integrabilidad Liouville y superintegrabilidad, ydiscutiremos un nuevo m\'etodo debido a F. Calogero para contruirmodelos de N-cuerpos en el plano con muchas \'orbitasperi\'odicas. La exposici\'on se acompa\~nar\'a de animaciones delmovimiento de los cuerpos, y se plantear\'an algunos problemasabiertos.\medskip\item17:00. Pausa\medskip\item17:30. Yuri Fedorov. \emph{An\'alisis de Kovalevskaya--Painlev\'ey Sistemas Algebraicamente Integrables}\medskip{\bf Resumen.} Muchos sistemas integrables poseen una propiedadremarcable: todas sus soluciones son funciones meromorfas deltiempo como una variable compleja. Tal comportamiento, que serefiere como propiedad de Kovalevskaya-Painleve (KP) y que se usafrecuentemente como una ensayo de integrabilidad, no es accidentaly tiene unas ra\'{\i}ces geom\'etricas profundas. En esta charladescribiremos una clase de tales sistemas (conocidos como lossistemas algebraicamente integrables) y subrayaremos suspropiedades geom\'etricas principales que permiten predecir laestructura de las soluciones complejas y adem\'as encontrarlasexpl\'{\i}citamente. Eso lo ilustraremos con algunos sistemas dela mec\'anica cl\'asica. Tambi\'en mencionaremos unasgeneralizaciones \'utiles de la noci\'on de integrabilidadalgebraica y de la propiedad KP.\end{itemize}\medskip{\bf Mi\'ercoles 30 de marzo}\begin{itemize}\item 15:30. Rafael Ram\'{\i}rez-Ros. \emph{El m\'etodo de Poincar\'e}\medskip{\bf Resumen.} Dado un sistema Hamiltoniano aut\'onomo cercano acompletamente integrable Poincar\'e prob\'o que, en general, noexiste ninguna integral primera adicional uniforme en elpar\'ametro de perturbaci\'on salvo el propio Hamiltoniano.Esbozaremos las ideas principales del m\'etodo de prueba ycomentaremos algunas extensiones y generalizaciones.\newpage\item16:30. Chara Pantazi. \emph{El M\'etodo de Darboux}\medskip{\bf Resumen.} Darboux, en 1878, present\'o su m\'etodo paraconstruir integrales primeras de campos vectoriales polinomialesutilizando sus curvas invariantes algebraicas. En estaexposici\'on presentaremos algunas extensiones del m\'etodocl\'asico de Darboux y tambi\'en algunas aplicaciones.\medskip\item17:30. Pausa\medskip\item18:00. Juan J. Morales-Ruiz. \emph{M\'etodos recientes paradetectar la no integrabilidad}\medskip{\bf Resumen.} En 1982 Ziglin utiliza la estructura de laecuaci\'on en variaciones de Poincar\'e (sobre una curva integralparticular) como una herramienta fundamental para detectar la nointegrabilidad de un sistema Hamiltoniano. En esta charla sepretende dar una idea de esta aproximaci\'on a la nointegrabilidad, junto con t\'ecnicas m\'as recientes queinvolucran la teor\'{\i}a de Galois de ecuaciones diferencialeslineales, haciendo \'enfasis en los ejemplos m\'as que en lateor\'{\i}a general. Ilustraremos estos m\'etodos con resultadossobre la no integrabilidad de algunos problemas de $N$ cuerpos enMec\'anica Celeste.\end{itemize}

Indirect likelihood inference (revised)

Standard indirect Inference (II) estimators take a given finite-dimensional statistic, Z_{n} , and then estimate the parameters by matching the sample statistic with the model-implied population moment. We here propose a novel estimation method that utilizes all available information contained in the distribution of Z_{n} , not just its first moment. This is done by computing the likelihood of Z_{n}, and then estimating the parameters by either maximizing the likelihood or computing the posterior mean for a given prior of the parameters. These are referred to as the maximum indirect likelihood (MIL) and Bayesian Indirect Likelihood (BIL) estimators, respectively. We show that the IL estimators are first-order equivalent to the corresponding moment-based II estimator that employs the optimal weighting matrix. However, due to higher-order features of Z_{n} , the IL estimators are higher order efficient relative to the standard II estimator. The likelihood of Z_{n} will in general be unknown and so simulated versions of IL estimators are developed. Monte Carlo results for a structural auction model and a DSGE model show that the proposed estimators indeed have attractive finite sample properties.

On equilibria in duopolies with finite strategy spaces

We will call a game a reachable (pure strategy) equilibria game if startingfrom any strategy by any player, by a sequence of best-response moves weare able to reach a (pure strategy) equilibrium. We give a characterizationof all finite strategy space duopolies with reachable equilibria. Wedescribe some applications of the sufficient conditions of the characterization.

Variational inequalities in Hilbert spaces with measures and optimal stopping problems

We study the existence theory for parabolic variational inequalities in weighted L2 spaces with respect to excessive measures associated with a transition semigroup. We characterize the value function of optimal stopping problems for finite and infinite dimensional diffusions as a generalized solution of such a variational inequality. The weighted L2 setting allows us to cover some singular cases, such as optimal stopping for stochastic equations with degenerate diffusion coeficient. As an application of the theory, we consider the pricing of American-style contingent claims. Among others, we treat the cases of assets with stochastic volatility and with path-dependent payoffs.

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.

Nash equilibrium strategies in discrete-time finite-horizon dynamic games with risk-and effort-averse players

The objective of this paper is to re-examine the risk-and effort attitude in the context of strategic dynamic interactions stated as a discrete-time finite-horizon Nash game. The analysis is based on the assumption that players are endogenously risk-and effort-averse. Each player is characterized by distinct risk-and effort-aversion types that are unknown to his opponent. The goal of the game is the optimal risk-and effort-sharing between the players. It generally depends on the individual strategies adopted and, implicitly, on the the players' types or characteristics.

Discontinuous Galerkin methods for the Multi-dimensional Vlasov-Poisson problem

We introduce and analyze two new semi-discrete numerical methods for the multi-dimensional Vlasov-Poisson system. The schemes are constructed by combing a discontinuous Galerkin approximation to the Vlasov equation together with a mixed finite element method for the Poisson problem. We show optimal error estimates in the case of smooth compactly supported initial data. We propose a scheme that preserves the total energy of the system.

Contact melting of a three-dimensional phase change material on a flat substrate

In this paper a model is developed to describe the three dimensional contact melting process of a cuboid on a heated surface. The mathematical description involves two heat equations (one in the solid and one in the melt), the Navier-Stokes equations for the flow in the melt, a Stefan condition at the phase change interface and a force balance between the weight of the solid and the countering pressure in the melt. In the solid an optimised heat balance integral method is used to approximate the temperature. In the liquid the small aspect ratio allows the Navier-Stokes and heat equations to be simplified considerably so that the liquid pressure may be determined using an igenfunction expansion and finally the problem is reduced to solving three first order ordinary differential equations. Results are presented showing the evolution of the melting process. Further reductions to the system are made to provide simple guidelines concerning the process. Comparison of the solutions with experimental data on the melting of n-octadecane shows excellent agreement.

Fronts from complex two-dimensional dispersal kernels: theory and application to Reid's paradox

Bimodal dispersal probability distributions with characteristic distances differing by several orders of magnitude have been derived and favorably compared to observations by Nathan [Nature (London) 418, 409 (2002)]. For such bimodal kernels, we show that two-dimensional molecular dynamics computer simulations are unable to yield accurate front speeds. Analytically, the usual continuous-space random walks (CSRWs) are applied to two dimensions. We also introduce discrete-space random walks and use them to check the CSRW results (because of the inefficiency of the numerical simulations). The physical results reported are shown to predict front speeds high enough to possibly explain Reid's paradox of rapid tree migration. We also show that, for a time-ordered evolution equation, fronts are always slower in two dimensions than in one dimension and that this difference is important both for unimodal and for bimodal kernels

Monte Carlo study of the relation between vacancy diffusion and domain growth in two-dimensional binary alloys

Domain growth in a two-dimensional binary alloy is studied by means of Monte Carlo simulation of an ABV model. The dynamics consists of exchanges of particles with a small concentration of vacancies. The influence of changing the vacancy concentration and finite-size effects has been analyzed. Features of the vacancy diffusion during domain growth are also studied. The anomalous character of the diffusion due to its correlation with local order is responsible for the obtained fast-growth behavior.

Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.

Spanning avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics: Field dependence and geometrical properties

Spanning avalanches in the 3D Gaussian Random Field Ising Model (3D-GRFIM) with metastable dynamics at T=0 have been studied. Statistical analysis of the field values for which avalanches occur has enabled a Finite-Size Scaling (FSS) study of the avalanche density to be performed. Furthermore, a direct measurement of the geometrical properties of the avalanches has confirmed an earlier hypothesis that several types of spanning avalanches with two different fractal dimensions coexist at the critical point. We finally compare the phase diagram of the 3D-GRFIM with metastable dynamics with the same model in equilibrium at T=0.

Diluted three-dimensional random field Ising model at zero temperature with metastable dynamics.

The influence of vacancy concentration on the behavior of the three-dimensional random field Ising model with metastable dynamics is studied. We have focused our analysis on the number of spanning avalanches which allows us a clear determination of the critical line where the hysteresis loops change from continuous to discontinuous. By a detailed finite-size scaling analysis we determine the phase diagram and numerically estimate the critical exponents along the whole critical line. Finally, we discuss the origin of the curvature of the critical line at high vacancy concentration.