Numerical Exploration of the RI/MOM Scheme Gauge Dependence

The gauge dependence of some fermion bilinear RI/MOM renormalization constants is studied by comparing data which have been gauge-fixed in two different realizations of the Landau gauge and in a generic covariant gauge. The very good agreement between the various sets of results and the theory indicates that the numerical uncertainty induced by the lattice gauge-fixing procedure is below the statistical errors of our data sample which is of the order of (1-1.5)%.Comment: 3 pages, 2 figures, Lattice2002(theoretical

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The system of equations consisting of the full Navier-Stokes equations and two turbulence equations was solved for in the steady state using a multigrid strategy on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time stepping scheme with a stability bound local time step, while the turbulence equations are advanced in a point-implicit scheme with a time step which guarantees stability and positively. Low Reynolds number modifications to the original two equation model are incorporated in a manner which results in well behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved for, initializing all quantities with uniform freestream values, and resulting in rapid and uniform convergence rates for the flow and turbulence equations

A High Statistics Lattice Calculation of Quark Masses with a Non-Perturbative Renormalization Procedure

We present results of a high statistics study (O(2000) configurations) of the quark masses in the MS-bar scheme from Lattice QCD in the quenched approximation at beta=6.0, beta=6.2 and beta=6.4 using both the Wilson and the tree-level improved SW-Clover fermion action. We extract quark masses from the meson spectroscopy and from the axial Ward Identity using non-perturbative values of the renormalization constants. We compare the results obtained with the two methods and we study the O(a) dependence of the quark masses for both actions. Our best results are m_s^(MS-bar)(2 GeV)=(123 +- 4 +- 15) MeV and m_c^(MS-bar)(2 GeV)=(1525 +- 40 +- 100) MeV.Comment: 3 pages, Talk presented at LATTICE97(QCD Spectrum and Quark Masses). Small misprints corrected, no conclusions change

Scheme Independence of the Effective Hamiltonian for b→s γb \rightarrow s \, \gamma and b→s gb \rightarrow s \, g Decays

We present a calculation of the effective weak Hamiltonian which governs $b \rightarrow s\, \gamma$ and $b \rightarrow s \, g$ transitions in two different renormalization schemes (NDR and HV). In the leading logarithmic approximation, we show that the coefficients of the effective Hamiltonian are scheme independent only when one takes correctly into account the scheme dependence of one- and two- loop diagrams. We demonstrate that in NDR there are contributions which were missed in previous calculations. These contributions are necessary to obtain scheme independent coefficients in the final results.Comment: 16 pp + 5 figures not included (available by anonymous ftp at amisan.iss.infn.it (141.108.15.215), directory /ftp/bsgamma), LaTeX, LPTENS 93/28, ROME 93/958, ULB-TH 93/0

RI/MOM Renormalization Window and Goldstone Pole Contamination

We perform a comparative study of the ratio of lattice (Wilson fermion) renormalization constants Z_P/Z_S, obtained non-perturbatively from the RI/MOM renormalization conditions and from Ward Identities of on- and off-shell Green's functions. The off-shell Ward Identity used in this work relies on correlation functions with non-degenerate quark masses. We find that, due to discretization effects, there is a 10-15% discrepancy between the two Ward Identity determinations at current bare couplings (beta values). The RI/MOM result is in the same range and has a similar systematic error of 10-15%. Thus, contrary to a previous claim, the contamination of the RI/MOM result from the presence of a Goldstone pole at scales of about 2 GeV is subdominant, compared to finite cutoff effects.Comment: LATEX, 12 pages final version to appear on Phys. Lett.

Charming Penguins in B decays

Full expressions of the $B^0_d \to \pi^+ \pi^-$ and $B^0_d \to \pi^0 \pi^0$ amplitudes, given in terms of matrix elements of operators of the effective weak Hamiltonian, are used to study the dependence of the relevant branching ratios on the different contributions. The uncertainty in the extraction of the weak phase $\alpha$ from the measurement of the time-dependent asymmetry in $B^0_d \to \pi^+ \pi^-$ decays is also analyzed. We find that, among several effects which may enhance the $B^0_d \to \pi^0 \pi^0$ branching ratio, the most important is due to ``charming penguin" diagrams that have never been studied before. These diagrams easily increase $BR(B^0_d \to \pi^0 \pi^0)$ up to a value of $1-3 \times 10^{-6}$. The same effect produces, however, a large error in the extraction of $\alpha$ from the measurement of the $B^0_d \to \pi^+ \pi^-$ time-dependent asymmetry. We show that it is possible to determine charming-penguin amplitudes from the experimental measurement of many decay rates. Their effect is impressive in $B^+ \to \pi^+ K^0$ and $B^0_d \to K^+ \pi^-$ decays, where charming-penguin contributions easily give values of $BR(B^+ \to \pi^+ K^0)$ and $BR(B^0_d \to K^+ \pi^-)$ of about $1 \times 10^{-5}$. Among other possibilities, we also suggest to use $B^0_d \to K^0 \bar K^0$, the BR of which can be as large as $2-3 \times 10^{-6}$, to determine the size of charming-penguin amplitudes.Comment: LaTeX, 28 pages, 8 figure