A statistical mechanical approach for the computation of the climatic response to general forcings

The climate belongs to the class of non-equilibrium forced and dissipative systems, for which most results of quasi-equilibrium statistical mechanics, including the fluctuation-dissipation theorem, do not apply. In this paper we show for the first time how the Ruelle linear response theory, developed for studying rigorously the impact of perturbations on general observables of non-equilibrium statistical mechanical systems, can be applied with great success to analyze the climatic response to general forcings. The crucial value of the Ruelle theory lies in the fact that it allows to compute the response of the system in terms of expectation values of explicit and computable functions of the phase space averaged over the invariant measure of the unperturbed state. We choose as test bed a classical version of the Lorenz 96 model, which, in spite of its simplicity, has a well-recognized prototypical value as it is a spatially extended one-dimensional model and presents the basic ingredients, such as dissipation, advection and the presence of an external forcing, of the actual atmosphere. We recapitulate the main aspects of the general response theory and propose some new general results. We then analyze the frequency dependence of the response of both local and global observables to perturbations having localized as well as global spatial patterns. We derive analytically several properties of the corresponding susceptibilities, such as asymptotic behavior, validity of Kramers-Kronig relations, and sum rules, whose main ingredient is the causality principle. We show that all the coefficients of the leading asymptotic expansions as well as the integral constraints can be written as linear function of parameters that describe the unperturbed properties of the system, such as its average energy. Some newly obtained empirical closure equations for such parameters allow to define such properties as an explicit function of the unperturbed forcing parameter alone for a general class of chaotic Lorenz 96 models. We then verify the theoretical predictions from the outputs of the simulations up to a high degree of precision. The theory is used to explain differences in the response of local and global observables, to define the intensive properties of the system, which do not depend on the spatial resolution of the Lorenz 96 model, and to generalize the concept of climate sensitivity to all time scales. We also show how to reconstruct the linear Green function, which maps perturbations of general time patterns into changes in the expectation value of the considered observable for finite as well as infinite time. Finally, we propose a simple yet general methodology to study general Climate Change problems on virtually any time scale by resorting to only well selected simulations, and by taking full advantage of ensemble methods. The specific case of globally averaged surface temperature response to a general pattern of change of the CO2 concentration is discussed. We believe that the proposed approach may constitute a mathematically rigorous and practically very effective way to approach the problem of climate sensitivity, climate prediction, and climate change from a radically new perspective.

Response theory for equilibrium and non-equilibrium statistical mechanics: causality and generalized Kramers-Kronig relations

We consider the general response theory recently proposed by Ruelle for describing the impact of small perturbations to the non-equilibrium steady states resulting from Axiom A dynamical systems. We show that the causality of the response functions entails the possibility of writing a set of Kramers-Kronig (K-K) relations for the corresponding susceptibilities at all orders of nonlinearity. Nonetheless, only a special class of directly observable susceptibilities obey K-K relations. Specific results are provided for the case of arbitrary order harmonic response, which allows for a very comprehensive K-K analysis and the establishment of sum rules connecting the asymptotic behavior of the harmonic generation susceptibility to the short-time response of the perturbed system. These results set in a more general theoretical framework previous findings obtained for optical systems and simple mechanical models, and shed light on the very general impact of considering the principle of causality for testing self-consistency: the described dispersion relations constitute unavoidable benchmarks that any experimental and model generated dataset must obey. The theory exposed in the present paper is dual to the time-dependent theory of perturbations to equilibrium states and to non-equilibrium steady states, and has in principle similar range of applicability and limitations. In order to connect the equilibrium and the non equilibrium steady state case, we show how to rewrite the classical response theory by Kubo so that response functions formally identical to those proposed by Ruelle, apart from the measure involved in the phase space integration, are obtained. These results, taking into account the chaotic hypothesis by Gallavotti and Cohen, might be relevant in several fields, including climate research. In particular, whereas the fluctuation-dissipation theorem does not work for non-equilibrium systems, because of the non-equivalence between internal and external fluctuations, K-K relations might be robust tools for the definition of a self-consistent theory of climate change.

Evidence of dispersion relations for the nonlinear response of the Lorenz 63 system

Along the lines of the nonlinear response theory developed by Ruelle, in a previous paper we have proved under rather general conditions that Kramers-Kronig dispersion relations and sum rules apply for a class of susceptibilities describing at any order of perturbation the response of Axiom A non equilibrium steady state systems to weak monochromatic forcings. We present here the first evidence of the validity of these integral relations for the linear and the second harmonic response for the perturbed Lorenz 63 system, by showing that numerical simulations agree up to high degree of accuracy with the theoretical predictions. Some new theoretical results, showing how to derive asymptotic behaviors and how to obtain recursively harmonic generation susceptibilities for general observables, are also presented. Our findings confirm the conceptual validity of the nonlinear response theory, suggest that the theory can be extended for more general non equilibrium steady state systems, and shed new light on the applicability of very general tools, based only upon the principle of causality, for diagnosing the behavior of perturbed chaotic systems and reconstructing their output signals, in situations where the fluctuation-dissipation relation is not of great help.

New exotic charmonium states

In the last years many states in the charmonium mass region were discovery by BABAR, Belle and CDF collaborations I discuss some of these discoveries, and how the QCD Sum Rule approach can be used to understand the structure of these states

Can the X(4350) narrow structure be a 1(-+) exotic state?

Using the QCD sum rules we test if the new narrow structure, the X(4350) recently observed by the Belle Collaboration, can be described as a J(PC) = 1(-+) exotic D(s)(*)D(s0)(*) molecular state. We consider the contributions of condensates up to dimension eight, we work at leading order in alpha(s) and we keep terms which are linear in the strange quark mass Ins. The mass obtained for such state is m(Ds*Ds0*) = (5.05 +/- 0.19) GeV. We also consider a molecular 1(-+), D(s)(*)D(s0)(*); current and we obtain m(D*D0*) = (4.92 +/- 0.08) GeV. We conclude that it is not possible to describe the X(4350) structure as a 1(-+) D(s)(*)D(s0)(*) molecular state. (C) 2010 Elsevier B.V. All rights reserved.

D(s)D* Molecule as an Axial Meson

We use QCD sum rules to study the possible existence of a D(s)(D) over bar* - D(s)*(D) over bar molecule with the quantum number J(PC) = 1(++). We consider the contributions of condensates up to dimension eight and work at leading order in alpha(s). We obtain m(DsD*) = (3.96 +/- 0.10) GeVaround 100 MeV above the mass of the meson X(3872). The proposed state is a natural generalized state to the strangeness sector of the X(3872), which was also found to be consistent with a multiquark state from a previous QCD sum rule analysis.

The meson Z(+) (4430) as a tetraquark state

We test the validity of the QCD sum rules applied to the meson Z(+)(4430). by considering a diquark-antidiquark type of current with J(P) = 0(-) and with J(P) = 1(-). We find that, with the studied currents, it is possible to find an acceptable Borel window. In such a Borel window we have simultaneously a good OPE convergence and a pole contribution which is bigger than the continuum contribution. We get m(z) = (4.52 +/- 0.09) GeV and m(Z) = (4.84 +/- 0.14) GeV for the currents with J(P) = 0(-) and J(P) = 1(-), respectively. We conclude that the QCD sum rules results favors J(P) = 0(-) quantum numbers for the Z(+) (4430) meson. (C) 2008 Elsevier B.V. All rights reserved.

Can the pi(+) chi(c1) resonance structures be D*(D)over-bar* and D(1)(D)over-bar molecules?

We use QCD sum rules to study the recently observed resonance-like structures in the pi(+)chi(c1) mass distribution, Z(1)(+) (4050) and Z(2)(+) (4250), considered as D*(+) (D) over bar*(0) and D(1)(+) (D) over bar (0) + D(+) (D) over bar (0)(1) molecules with the quantum number J(P) = 0(+) and J(P) = 1-, respectively. We consider the contributions of condensates up to dimension eight and work at leading order in alpha(s). We obtain m(D*D*) = (4.15 +/- 0.12) GeV, around 100 MeV above the D*D* threshold, and m(D1D) = (4.19 +/- 0.22) GeV, around 100 MeV below the D(1)D threshold. We conclude that the D*(+)(D) over bar*(0) state is probably a virtual state that is not related with the Z(1)(+) (4050) resonance-like structure. In the case of the D(1)D molecular state, considering the errors, its mass is consistent with both Z(1)(+)(4050) and Z(2)(+)(4250) resonance-like structures. Therefore, we conclude that no definite conclusion can be drawn for this state from the present analysis. (C) 2008 Elsevier B.V All rights reserved.

Five-dimensional f(R) braneworld models

After incorporating f(R) gravity into the general braneworld sum rules scope, it is shown that some particular class of warped five-dimensional nonlinear braneworld models, which may be interesting for the hierarchy problem solution, still require a negative tension brane. For other classes of warp factors (suitable and not suitable for approaching the hierarchy problem) any negative brane tension in the compactification scheme is not necessary. In this vein, it is argued that, in the bulk f(R) gravity context, some types of warp factors may be useful for approaching the hierarchy problem and for evading the necessity of a negative brane tension in the compactification scheme.

Hadronic form factors and the J/psi secondary production cross section: An update

Improving previous calculations, we compute the D + (D) over bar J/psi + pi cross section using the most complete effective Lagrangians available. The new crucial ingredients are the form factors on the charm meson vertices, which are determined from QCD sum rules calculations. Some of them became available only very recently, and the last one, needed for our present purpose, is calculated in this work.

Dressed coordinates: the path-integral approach

The recently introduced dressed coordinates are studied in the path-integral approach. These coordinates are defined in the context of a harmonic oscillator linearly coupled to massless scalar field and it is shown that in this model the dressed coordinates appear as a coordinate transformation preserving the path-integral functional measure. The analysis also generalizes the sum rules established in a previous work. (c) 2006 Elsevier B.V. All rights reserved.

Topological charge screening and pseudoscalar glueballs

Topological charge screening in the QCD vacuum is found to provide crucial nonperturbative contributions to the short-distance expansion of the pseudoscalar (0-+) glueball correlator. The screening contributions enter the Wilson coefficients and are an indispensable complement to the direct instanton contributions. They restore consistency with the anomalous axial Ward identity and remedy several flaws in the 0-+ glueball sum rules caused by direct instantons in the absence of screening (lack of resonance signals, violation of the positivity bound and of the underlying low-energy theorem). The impact of the finite width of the instanton size distribution and the (gauge-invariant) renormalization of the instanton contributions are also discussed. New predictions for the 0-+ glueball mass and decay constant are presented.

Notes on the two-brane model with variable tension

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Quark-antiquark correlation in the pion

The electromagnetic tensor for inclusive electron scattering off the pion Wμν for momentum transfers such that q+ = 0, (q+ = q0 + q3) is shown to obey a sum-rule for the component W++. From this sum-rule, one can define the quark-antiquark correlation function in the pion, which characterizes the transverse distance distribution between the quark and antiquark in the light-front pion wave-function. Within the realistic models of the relativistic pion wave function (including instanton vacuum inspired wave function) it is shown that the value of the two-quark correlation radius (rqq̄) is near twice the pion electromagnetic radius (rπ), where rπ ≈ 2/3 fm. We also define the correlation length lcorr where the two-particle correlation have an extremum. The estimation of lcorr ≈ 0.3-0,5 fm is very close to estimations from instanton models of QCD vacuum. It is also shown that the above correlation is very sensitive to the pion light-front wave-function models. © 1997 Elsevier Science B.V.

Uses of QCD sumrules (QCDSR) from different points of view

We support the idea that the baryon, B with mass MB, couples to its current with a coupling λ2 B ∼ 0.71 M6 B from an analysis of magnetic moment sum rules. And we find a sum rule among the experimental magnetic moments which is independent of the parameters of QCDSR. © 1998 Elsevier Science B.V. All rights reserved.