A non-perturbative estimate of the heavy quark momentum diffusion coefficient

We estimate the momentum diffusion coefficient of a heavy quark within a pure SU(3) plasma at a temperature of about 1.5Tc. Large-scale Monte Carlo simulations on a series of lattices extending up to 192^3*48 permit us to carry out a continuum extrapolation of the so-called colour-electric imaginary-time correlator. The extrapolated correlator is analyzed with the help of theoretically motivated models for the corresponding spectral function. Evidence for a non-zero transport coefficient is found and, incorporating systematic uncertainties reflecting model assumptions, we obtain kappa = (1.8 - 3.4)T^3. This implies that the "drag coefficient", characterizing the time scale at which heavy quarks adjust to hydrodynamic flow, is (1.8 - 3.4) (Tc/T)^2 (M/1.5GeV) fm/c, where M is the heavy quark kinetic mass. The results apply to bottom and, with somewhat larger systematic uncertainties, to charm quarks.Comment: 18 pages. v2: clarifications adde

Critical point and scale setting in SU(3) plasma: An update

We explore a method developed in statistical physics which has been argued to have exponentially small finite-volume effects, in order to determine the critical temperature Tc of pure SU(3) gauge theory close to the continuum limit. The method allows us to estimate the critical coupling betac of the Wilson action for temporal extents up to Nt ~ 20 with < 0.1% uncertainties. Making use of the scale setting parameters r0 and sqrt{t0} in the same range of beta-values, these results lead to the independent continuum extrapolations Tc r0 = 0.7457(45) and Tc sqrt{t0} = 0.2489(14), with the latter originating from a more convincing fit. Inserting a conversion of r0 from literature (unfortunately with much larger errors) yields Tc / LambdaMSbar = 1.24(10).Comment: 12 pages. v2: clarifications and references added, published versio

Towards the continuum limit in transport coefficient computations

The analytic continuation needed for the extraction of transport coefficients necessitates in principle a continuous function of the Euclidean time variable. We report on progress towards achieving the continuum limit for 2-point correlator measurements in thermal SU(3) gauge theory, with specific attention paid to scale setting. In particular, we improve upon the determination of the critical lattice coupling and the critical temperature of pure SU(3) gauge theory, estimating r0*Tc ~ 0.7470(7) after a continuum extrapolation. As an application the determination of the heavy quark momentum diffusion coefficient from a correlator of colour-electric fields attached to a Polyakov loop is discussed.Comment: 7 pages. To appear in the Proceedings of the 31st International Symposium on Lattice Field Theory, July 29 - August 3, 2013, Mainz, German

The variability of peridotite composition across a mantle shear zone (Lanzo massif, Italy): interplay of melt focusing and deformation

In this paper we present new data on the spatial variability of peridotite composition across a kilometer-scale mantle shear zone within the Lanzo massif (Western Alps, Italy). The shear zone separates the central from the northern part of the massif. Plagioclase peridotite shows gradually increasing deformation towards the shear zone, from porphyroclastic to mylonitic textures in the central body, while the northern body is composed of porphyroclastic rocks. The peridotite displays a large range of compositions, from fertile peridotite to refractory harzburgite and dunite. Deformed peridotites (proto-mylonite and mylonites) tend to be compositionally more homogeneous and fertile than weakly deformed peridotites. The composition of most plagioclase peridotites show rather high and constant (Ce/Yb)N ratios, and YbN that cannot be explained by any simple melting model. Instead, refertilization modeling, consisting of melt increments from spinel peridotite sources, particularly with E-MORB melt, reasonably reproduces the plagioclase peridotite whole rock composition. Combined with constraints from Ce-Nb and Ce-Th systematics, we speculate that peridotites such as those from Lanzo record pervasive refertilization processes in the thermal boundary layer. In this scenario, mantle shear zones might act as important areas of melt focusing in the upper mantle that separates the thermal boundary layer from the conductively cooled mantl

Interplay between chiral and axial symmetries in a SU(2) Nambu--Jona-Lasinio Model with the Polyakov loop

We consider a two flavor Polyakov--Nambu--Jona-Lasinio (PNJL) model where the Lagrangian includes an interaction term that explicitly breaks the U$_A(1)$ anomaly. At finite temperature, the restoration of chiral and axial symmetries, signaled by the behavior of several observables, is investigated. We compare the effects of two regularizations at finite temperature, one of them, that allows high momentum quarks states, leading to the full recovery of chiral symmetry. From the analysis of the behavior of the topological susceptibility and of the mesonic masses of the axial partners, it is found in the SU(2) model that, unlike the SU(3) results, the recovery of the axial symmetry is not a consequence of the full recovery of the chiral symmetry. Thus, one needs to use an additional idea, by means of a temperature dependence of the anomaly coefficient, that simulates instanton suppression effects.Comment: 21 pages, 5 figures; PRD versio

Electric Flux Tube in Magnetic Plasma

In this paper we study a methodical problem related to the magnetic scenario recently suggested and initiated by the authors \cite{Liao_ES_mono} to understand the strongly coupled quark-gluon plasma (sQGP): the electric flux tube in monopole plasma. A macroscopic approach, interpolating between Bose condensed (dual superconductor) and classical gas medium is developed first. Then we work out a microscopic approach based on detailed quantum mechanical calculation of the monopole scattering on electric flux tube, evaluating induced currents for all partial waves. As expected, the flux tube looses its stability when particles can penetrate it: we make this condition precise by calculating the critical value for the product of the flux tube size times the particle momentum, above which the flux tube dissolves. Lattice static potentials indicate that flux tubes seem to dissolve at $T>T_{dissolution} \approx 1.3 T_c$. Using our criterion one gets an estimate of the magnetic density $n\approx 4.4 \sim 6.6 fm^{-3}$ at this temperature.Comment: New version with new referecences added and minor changes. 15 pages, 8 figure

Trace element chemistry and U-Pb dating of zircons from oceanic gabbros and their relationship with whole rock composition (Lanzo, Italian Alps)

The U-Pb ages and the trace element content of zircon U-Pb along with major and trace element whole rock data on gabbroic dikes from the Lanzo lherzolitic massif, N-Italy, have been determined to constrain crustal accretion in ocean-continent transition zones. Three Fe-Ti gabbros were dated from the central and the southern part of the massif providing middle Jurassic ages of 161±2, 158±2 and 163±1Ma, which argue for magmatic activity over few millions of years. Zircon crystals are characterized by high but variable Th/U ratios, rare earth element patterns enriched in heavy rare earths, pronounced positive Ce and negative Eu-anomalies consistent with crystallization after substantial plagioclase fractionation. The zircon trace element composition coupled with whole rock chemistry was used to reconstruct the crystallization history of the gabbros. A number of gabbros crystallized in situ, and zircon precipitated from trapped, intercumulus liquid, while other gabbros represent residual liquids that were extracted from a cumulus pile and crystallized along syn-magmatic shear zones. We propose a model in which the emplacement mechanism of gabbroic rocks in ocean-continent transition zones evolves from in situ crystallization to stratified crystallization with efficient extraction of residual liquid along syn-magmatic shear zones. Such an evolution of the crystallization history is probably related to the thermal evolution of the underlying mantle lithospher

Control of serpentinisation rate by reaction-induced cracking

Serpentinisation of mantle rocks requires the generation and maintenance of transport pathways for water. The solid volume increase during serpentinisation can lead to stress build-up and trigger cracking, which ease fluid penetration into the rock. The quantitative effect of this reaction-induced cracking mechanism on reactive surface generation is poorly constrained, thus hampering our ability to predict serpentinisation rate in geological environments. Here we use a combined approach with numerical modelling and observations in natural samples to provide estimates of serpentinisation rate at mid-ocean ridges. We develop a micromechanical model to quantify the propagation of serpentinisation-induced cracks in olivine. The maximum crystallisation pressure deduced from thermodynamic calculations reaches several hundreds of megapascals but does not necessary lead to crack propagation if the olivine grain is subjected to high compressive stresses. The micromechanical model is then coupled to a simple geometrical model to predict reactive surface area formation during grain splitting, and thus bulk reaction rate. Our model reproduces quantitatively experimental kinetic data and the typical mesh texture formed during serpentinisation. We also compare the model results with olivine grain size distribution data obtained on natural serpentinised peridotites from the Marum ophiolite and the Papuan ultramafic belt (Papua New Guinea). The natural serpentinised peridotites show an increase of the number of olivine grains for a decrease of the mean grain size by one order of magnitude as reaction progresses from 5 to 40%. These results are in agreement with our model predictions, suggesting that reaction-induced cracking controls the serpentinisation rate. We use our model to estimate that, at mid-ocean ridges, serpentinisation occurs up to 12km depth and reaction-induced cracking reduces the characteristic time of serpentinisation by one order of magnitude, down to values comprised between 10 and 1000yr. The increase of effective pressure with depth also prevents cracking, which positions the peak in serpentinisation rate at shallower depths, 4km above previous predictions