Numerical precision calculations for LHC physics

In dieser Arbeit stelle ich Aspekte zu QCD Berechnungen vor, welche eng verknüpft sind mit der numerischen Auswertung von NLO QCD Amplituden, speziell der entsprechenden Einschleifenbeiträge, und der effizienten Berechnung von damit verbundenen Beschleunigerobservablen. Zwei Themen haben sich in der vorliegenden Arbeit dabei herauskristallisiert, welche den Hauptteil der Arbeit konstituieren. Ein großer Teil konzentriert sich dabei auf das gruppentheoretische Verhalten von Einschleifenamplituden in QCD, um einen Weg zu finden die assoziierten Farbfreiheitsgrade korrekt und effizient zu behandeln. Zu diesem Zweck wird eine neue Herangehensweise eingeführt welche benutzt werden kann, um farbgeordnete Einschleifenpartialamplituden mit mehreren Quark-Antiquark Paaren durch Shufflesummation über zyklisch geordnete primitive Einschleifenamplituden auszudrücken. Ein zweiter großer Teil konzentriert sich auf die lokale Subtraktion von zu Divergenzen führenden Poltermen in primitiven Einschleifenamplituden. Hierbei wurde im Speziellen eine Methode entwickelt, um die primitiven Einchleifenamplituden lokal zu renormieren, welche lokale UV Counterterme und effiziente rekursive Routinen benutzt. Zusammen mit geeigneten lokalen soften und kollinearen Subtraktionstermen wird die Subtraktionsmethode dadurch auf den virtuellen Teil in der Berechnung von NLO Observablen erweitert, was die voll numerische Auswertung der Einschleifenintegrale in den virtuellen Beiträgen der NLO Observablen ermöglicht. Die Methode wurde schließlich erfolgreich auf die Berechnung von NLO Jetraten in Elektron-Positron Annihilation im farbführenden Limes angewandt.

Study of the helicity dependence of pi 0 photoproduction on proton at MAMI

The excitation spectrum is one of the fundamental properties of every spatially extended system. The excitations of the building blocks of normal matter, i.e., protons and neutrons (nucleons), play an important role in our understanding of the low energy regime of the strong interaction. Due to the large coupling, perturbative solutions of quantum chromodynamics (QCD) are not appropriate to calculate long-range phenomena of hadrons. For many years, constituent quark models were used to understand the excitation spectra. Recently, calculations in lattice QCD make first connections between excited nucleons and the fundamental field quanta (quarks and gluons). Due to their short lifetime and large decay width, excited nucleons appear as resonances in scattering processes like pion nucleon scattering or meson photoproduction. In order to disentangle individual resonances with definite spin and parity in experimental data, partial wave analyses are necessary. Unique solutions in these analyses can only be expected if sufficient empirical information about spin degrees of freedom is available. The measurement of spin observables in pion photoproduction is the focus of this thesis. The polarized electron beam of the Mainz Microtron (MAMI) was used to produce high-intensity, polarized photon beams with tagged energies up to 1.47 GeV. A "frozen-spin" Butanol target in combination with an almost 4π detector setup consisting of the Crystal Ball and the TAPS calorimeters allowed the precise determination of the helicity dependence of the γp → π0p reaction. In this thesis, as an improvement of the target setup, an internal polarizing solenoid has been constructed and tested. A magnetic field of 2.32 T and homogeneity of 1.22×10−3 in the target volume have been achieved. The helicity asymmetry E, i.e., the difference of events with total helicity 1/2 and 3/2 divided by the sum, was determined from data taken in the years 2013-14. The subtraction of background events arising from nucleons bound in Carbon and Oxygen was an important part of the analysis. The results for the asymmetry E are compared to existing data and predictions from various models. The results show a reasonable agreement to the models in the energy region of the ∆(1232)-resonance but large discrepancies are observed for energy above 600 MeV. The expansion of the present data in terms of Legendre polynomials, shows the sensitivity of the data to partial wave amplitudes up to F-waves. Additionally, a first, preliminary multipole analysis of the present data together with other results from the Crystal Ball experiment has been as been performed.

Hard thermal loop benchmark for the extraction of the nonperturbative QQ[over ¯] potential

The extraction of the finite temperature heavy quark potential from lattice QCD relies on a spectral analysis of the Wilson loop. General arguments tell us that the lowest lying spectral peak encodes, through its position and shape, the real and imaginary parts of this complex potential. Here we benchmark this extraction strategy using leading order hard-thermal loop (HTL) calculations. In other words, we analytically calculate the Wilson loop and determine the corresponding spectrum. By fitting its lowest lying peak we obtain the real and imaginary parts and confirm that the knowledge of the lowest peak alone is sufficient for obtaining the potential. Access to the full spectrum allows an investigation of spectral features that do not contribute to the potential but can pose a challenge to numerical attempts of an analytic continuation from imaginary time data. Differences in these contributions between the Wilson loop and gauge fixed Wilson line correlators are discussed. To better understand the difficulties in a numerical extraction we deploy the maximum entropy method with extended search space to HTL correlators in Euclidean time and observe how well the known spectral function and values for the real and imaginary parts are reproduced. Possible venues for improvement of the extraction strategy are discussed.

Benchmarking the Bayesian reconstruction of the non-perturbative heavy ǪǬ potential

The extraction of the finite temperature heavy quark potential from lattice QCD relies on a spectral analysis of the real-time Wilson loop. Through its position and shape, the lowest lying spectral peak encodes the real and imaginary part of this complex potential. We benchmark this extraction strategy using leading order hard-thermal loop (HTL) calculations. I.e. we analytically calculate the Wilson loop and determine the corresponding spectrum. By fitting its lowest lying peak we obtain the real- and imaginary part and confirm that the knowledge of the lowest peak alone is sufficient for obtaining the potential. We deploy a novel Bayesian approach to the reconstruction of spectral functions to HTL correlators in Euclidean time and observe how well the known spectral function and values for the real and imaginary part are reproduced. Finally we apply the method to quenched lattice QCD data and perform an improved estimate of both real and imaginary part of the non-perturbative heavy ǪǬ potential.

Optimization of multidimensional cross-section tables for few-group core calculations

Multigroup diffusion codes for three dimensional LWR core analysis use as input data pre-generated homogenized few group cross sections and discontinuity factors for certain combinations of state variables, such as temperatures or densities. The simplest way of compiling those data are tabulated libraries, where a grid covering the domain of state variables is defined and the homogenized cross sections are computed at the grid points. Then, during the core calculation, an interpolation algorithm is used to compute the cross sections from the table values. Since interpolation errors depend on the distance between the grid points, a determined refinement of the mesh is required to reach a target accuracy, which could lead to large data storage volume and a large number of lattice transport calculations. In this paper, a simple and effective procedure to optimize the distribution of grid points for tabulated libraries is presented. Optimality is considered in the sense of building a non-uniform point distribution with the minimum number of grid points for each state variable satisfying a given target accuracy in k-effective. The procedure consists of determining the sensitivity coefficients of k-effective to cross sections using perturbation theory; and estimating the interpolation errors committed with different mesh steps for each state variable. These results allow evaluating the influence of interpolation errors of each cross section on k-effective for any combination of state variables, and estimating the optimal distance between grid points.

Topochemical dimerization of non-parallel double bonds: 7-methoxycoumarin

X-ray crystal structure analysis of 7-methoxycoumarin reveals that the reactive double bonds are rotated by about 65° with respect to each other, the centre-to-centre distance between the double bonds being 3.83 Å. In spite of this unfavourable arrangement, photodimerization occurs in the crystalline state yielding the syn-head-tail dimer as the only product. Lattice energy calculations on ground-state molecules in crystals throw light on the mechanism of the reaction.

Solid-state photobehaviour and crystal packing of o-chlorobenzylidene-DL-piperitone: influence of molecular topology on photobehaviour

C17H19ClO, M(r) = 274.7, triclinic, P1BAR, a = 11.154 (3), b = 12.685 (2), c = 12.713 (2) angstrom, alpha = 100.68 (1), beta = 113.58 (1), gamma = 104.50 (2)-degrees, V = 1511.1 (6) angstrom3, Z = 4, D(m) = 1.22, D(x) = 1.215 Mg m-3, Cu K-alpha, lambda = 1.5418 angstrom, mu = 2.16 mm-1, F(000) = 584, T = 293 K, R = 0.057 for 3481 observed reflections. The title compound is photostable in the crystalline state and lattice-energy calculations have been employed to rationalize the photobehaviour. The well-known beta-steering ability of the chloro group is not operative in this system as there are no Cl...Cl interactions in the crystal lattice. All five benzylidene-DL-piperitone structures so far studied are alpha-packed and the molecular topology appears to be a deciding factor even in the presence of steering groups.

Conformational Polymorphism in a Non-steroidal Anti-inflammatory Drug, Mefenamic Acid

For several years, the non-steroidal anti-inflammatory drug mefenamic acid, MA, has been known to exist as dimorphs (I and II). We report a new metastable polymorph (III) of MA obtained during attempted co-crystallization experiments and establish its stability relationship with existing forms. At elevated temperatures I and III convert to II, as evident from DSC experiments. On the basis of the lattice energy calculations in conjunction with thermal analysis, the stability order is proposed to be I > II > III at ambient conditions, whereas at elevated temperature the order is II > I > III. In either condition III is a metastable form and hence transforms to I at ambient conditions and to II at higher temperatures. Also we report the structural studies of a DMF solvate and a cytosine complex.

gamma cross section in p plus p collisions at root s=200 GeV

We report on a measurement of the gamma(1S + 2S + 3S) -> e(+)e(-) cross section at midrapidity in p + p collisions at root s = 200 GeV. We find the cross section to be 114 +/- 38(stat + fit)(-24)(+23)(syst) pb. Perturbative QCD calculations at next-to-leading order in the color evaporation model are in agreement with our measurement, while calculations in the color singlet model underestimate it by 2 sigma. Our result is consistent with the trend seen in world data as a function of the center-of-mass energy of the collision and extends the availability of gamma data to RHIC energies. The dielectron continuum in the invariant-mass range near the gamma is also studied to obtain a combined yield of e(+)e(-) pairs from the sum of the Drell-Yan process and b-(b) over bar production.

Inclusive pi(0), eta, and direct photon production at high transverse momentum in p plus p and d plus Au collisions at root s(NN)=200 GeV

We report a measurement of high-p(T) inclusive pi(0), eta, and direct photon production in p + p and d + Au collisions at root s(NN) = 200 GeV at midrapidity (0 < eta < 1). Photons from the decay pi(0) -> gamma gamma were detected in the barrel electromagnetic calorimeter of the STAR experiment at the Relativistic Heavy Ion Collider. The eta -> gamma gamma decay was also observed and constituted the first eta measurement by STAR. The first direct photon cross-section measurement by STAR is also presented; the signal was extracted statistically by subtracting the pi(0), eta, and omega(782) decay background from the inclusive photon distribution observed in the calorimeter. The analysis is described in detail, and the results are found to be in good agreement with earlier measurements and with next-to-leading-order perturbative QCD calculations.

Role of five-quark components in radiative and strong decays of the Lambda(1405) resonance

Within an extended chiral constituent quark model, the three- and five-quark structure of the S-01 resonance Lambda(1405) is investigated. Helicity amplitudes for electromagnetic decays [Lambda(1405)->Lambda(1116)gamma, Sigma(1194)gamma] and transition amplitudes for strong decays [Lambda(1405)->Sigma(1194)pi, K- p] are derived, as well as the relevant decay widths. The experimental value for the strong decay width, Gamma(Lambda(1405)->(Sigma pi)degrees) = 50 +/- 2MeV, is well reproduced with about 50% of a five-quark admixture in the Lambda(1405). Important effects owing to the configuration mixing among Lambda P-2(1)A, Lambda P-2(8)M, and Lambda P-4(8)M are found. In addition, transitions between the three- and the five-quark components in the baryons turn out to be significant in both radiative and strong decays of the Lambda(1405) resonance.

Mesure du rapport d'embranchement des désintégrations B^0--> pi^- l^+v et B^+--> êta^(') l^+v, du spectre des facteurs de forme des désintégrations B^0--> pi^- l^+v et B^+--> êta l^+v, et extraction de |V_ub| dans l'expérience BABAR

Nous rapportons les résultats d'une étude des désintégrations semileptoniques non-charmées B^+--> êta^(') l^+v et B^0--> pi^- l^+v, mesurés par le détecteur BABAR avec une production d'environ 464 millions de paires de mésons BBbar issues des collisions e^+e^- à la résonance Upsilon(4S). L'analyse reconstruit les événements avec une technique relâchée des neutrinos. Nous obtenons les rapports d'embranchement partiels pour les désintégrations B^+--> êta l^+v et B^0--> pi^- l^+v en trois et douze intervalles de q^2, respectivement, à partir desquels nous extrayons les facteurs de forme f_+(q^2) et les rapports d'embranchement totaux B(B^+--> êta l^+v) = (3.39 +/- 0.46_stat +/- 0.47_syst) x 10^-5 et B(B^0--> pi^- l^+v) = (1.42 +/- 0.05_stat +/- 0.08_syst) x 10^-4. Nous mesurons aussi B(B^+--> êta' l^+v) = (2.43 +/- 0.80_stat +/- 0.34_syst) x 10^-5. Nous obtenons les valeurs de la norme de l'élément |V_ub| de la matrice CKM en utilisant trois calculs différents de la CDQ.

Thermal photon production in Au plus Au collisions: Viscous corrections in two different hydrodynamic formalisms

We calculate the spectra of produced thermal photons in Au + Au collisions taking into account the nonequilibrium contribution to photon production due to finite shear viscosity. The evolution of the fireball is modeled by second-order as well as by divergence-type 2 + 1 dissipative hydrodynamics, both with an ideal equation of state and with one based on Lattice QCD that includes an analytical crossover. The spectrum calculated in the divergence-type theory is considerably enhanced with respect to the one calculated in the second-order theory, the difference being entirely due to differences in the viscous corrections to photon production. Our results show that the differences in hydrodynamic formalisms are an important source of uncertainty in the extraction of the value of eta/s from measured photon spectra. The uncertainty in the value of eta/s associated with different hydrodynamic models used to compute thermal photon spectra is larger than the one occurring in matching hadron elliptic flow to RHIC data. (C) 2010 Elsevier B.V. All rights reserved.

Dilepton mass spectra in p plus p collisions at root s=200 GeV

PHENIX has measured the electron-positron pair mass spectrum from 0 to 8 GeV/c(2) in p + p collisions at root s = 200 GeV. The contributions from light meson decays to e(+)e(-) pairs have been determined based on measurements of hadron production cross sections by PHENIX. Within the systematic uncertainty of similar to 20% they account for all e(+)e(-) pairs in the mass region below similar to 1 GeV/c(2). The e(+)e(-) pair yield remaining after subtracting these contributions is dominated by semileptonic decays of charmed hadrons correlated through flavor conservation. Using the spectral shape predicted by PYTHIA, we estimate the charm production cross section to be 544 +/- 39(stat) +/- 142(syst) +/- 200(model) pb. which is consistent with QCD calculations and measurements of single leptons by PHENIX. (C) 2008 Elsevier BV. All rights reserved.

NeXSPheRIO results on elliptic-flow fluctuations at RHIC

By using the NeXSPheRIO code, we study the elliptic-flow fluctuations in Au + Au collisions at 200 A GeV. It is shown that, by fixing the parameters of the model to correctly reproduce the charged pseudorapidity and the transverse-momentum distributions, reasonable agreement of < v(2)> with data is obtained, both as function of pseudorapidity as well as of transverse momentum, for charged particles. Our results on elliptic-flow fluctuations are in good agreement with the recently measured data on experiments.