Heavy quark and quarkonium production at CERN LEP2: <Emphasis Type="Italic">k</Emphasis><Subscript>T</Subscript>-factorization versus data

We present calculations of heavy quark and quarkonium production at CERN LEP2 in the k _T-factorization QCD approach. Both direct and resolved photon contributions are taken into account. A conservative error analysis is performed. The unintegrated gluon distribution in the photon is obtained from the full CCFM evolution equation. The traditional color-singlet mechanism to describe the non-perturbative transition of a -pair into a final quarkonium is used. Our analysis covers the polarization properties of heavy quarkonia at moderate and large transverse momenta. We find that the total and differential open charm production cross sections are consistent with the recent experimental data taken by the L3, OPAL and ALEPH collaborations. At the same time the DELPHI data for the inclusive production exceed our predictions, but experimental uncertainties are too large to claim a significant inconsistency. The bottom production in photon-photon collisions at CERN LEP2 is hard to explain within the k _T-factorization formalism.Received: 22 December 2004, Revised: 15 February 2005, Published online: 12 April 2005     

Shadowing correction to the gluon distribution behavior at small x

We determined the saturation exponent of the gluon distribution using the solution of the QCD nonlinear Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (NLDGLAP) evolution equation at small x . The very small-x behavior of the gluon distribution is obtained by solving the Gribov, Levin, Ryskin, Mueller and Qiu (GLR-MQ) evolution equation with the nonlinear shadowing term incorporated. The form of the initial condition for the equation is determined. We find, with decreasing x , the emergence of a singular behavior and the eventual taming (at R = 5 GeV^-1) and the essential taming (at R = 2 GeV^-1) of this singular behavior by the shadowing term. The nonlinear gluon density functions are calculated and compared with the results for the integrated gluon density from the Balitsky-Kovchegov (BK) equation for the different values of Q². It is shown that the results for the gluon density function are comparable with the results obtained from the BK equation solution. Also we show that for each x , the Q²-dependence of the data is well described by gluon shadowing corrections to the GLR-MQ equation. The resulting analytic expression allows us to predict the logarithmic derivative \( {\frac{{\partial F^{s}_{2} (x,Q^{2})}}{{\partial \ln Q^{2}}}}\) and to compare the results with H1 data and a QCD analysis fit.     

$\mbox{Log}(1/x)$ gluon distribution and structure functionsin the loop-loop correlation model

We consider the interaction of the partonic fluctuation of a scalar photon with an external color field to calculate the leading and next-to-leading order gluon distribution of the proton following the work done by Dosch-Hebecker-Metz-Pirner. We relate these gluon distributions to the short and long distance behavior of the cross section of an adjoint dipole scattering off a proton. The leading order result is a constant, while the next-to-leading order result shows a enhancement at small x. To get numerical results for the gluon distributions at the initial scale Q²₀ = 1.8 GeV², we compute the adjoint dipole-proton cross section in the loop-loop correlation model. Quark distributions at the same initial scale are parameterized according to Regge theory. We evolve quark and gluon distributions to higher Q² values using the DGLAP equation and compute charm and proton structure functions in the small-x region for different Q² values.Received: 13 September 2003, Revised: 22 November 2003, Published online: 15 January 2003     

Estimating nonlinear effects in forward dijet production in ultra-peripheral heavy ion collisions at the LHC

Using the framework that interpolates between the leading power limit of the color glass condensate and the high energy (or \(k_{T}\)) factorization we calculate the direct component of the forward dijet production in ultra-peripheral \(\mathrm {Pb}\)–\(\mathrm {Pb}\) collisions at CM energy \(5.1\,\mathrm {TeV}\) per nucleon pair. The formalism is applicable when the average transverse momentum of the dijet system \(P_{T}\) is much bigger than the saturation scale \(Q_{s}\), \(P_{T}\gg Q_{s}\), while the imbalance of the dijet system can be arbitrary. The cross section is uniquely sensitive to the Weizsäcker–Williams (WW) unintegrated gluon distribution, which is far less known from experimental data than the most common dipole gluon distribution appearing in inclusive small-x processes. We have calculated cross sections and nuclear modification ratios using WW gluon distribution obtained from the dipole gluon density through the Gaussian approximation. The dipole gluon distribution used to get WW was fitted to the inclusive HERA data with the nonlinear extension of unified BFKL + DGLAP evolution equation. The saturation effects are visible but rather weak for realistic \(p_{T}\) cut on the dijet system, reaching about 20% with the cut as low as \(6\,\mathrm {GeV}\). We find that the LO collinear factorization with nuclear leading-twist shadowing predicts quite similar effects.     

A new approach to calculate the gluon polarization

We derive the Leading-Order (LO) master equation to extract the polarized gluon distribution G(x,Q ²)=xδg(x,Q ²) from polarized proton structure function, g^p₁(x,Q²)g^{p}_{1}(x,Q^{2}). By using a Laplace-transform technique, we solve the master equation and derive the polarized gluon distribution inside the proton. The test of accuracy which is based in our calculations on two different methods, confirms that we achieve to the correct solution for the polarized gluon distribution. To determine the polarized gluon distribution xδg(x,Q ²) more accurately, we only need to have more experimental data on the polarized structure functions, g₁^p(x,Q²)g_{1}^{p}(x,Q^{2}). Our result for polarized gluon distribution is in good agreement with some phenomenological models.     

Higgs boson production at hadron collidersin the kT-factorization approach

We consider the Higgs boson production at high energy hadron colliders in the framework of the k_T-factorization approach. The attention is focused on the dominant gluon-gluon fusion subprocess. We calculate the total cross section and transverse momentum distributions of the inclusive Higgs production using unintegrated gluon distributions in a proton obtained from the full CCFM evolution equation. We show that k_T-factorization gives a possibility to investigate the associated Higgs boson and jets production. We calculate the transverse momentum distributions and study the Higgs-jet and jet-jet azimuthal correlations in the Higgs + one or two jet production processes. We demonstrate the importance of the higher-order corrections within the k_T-factorization approach. These corrections should be developed and taken into account in the future applications. Received: 26 January 2005, Revised: 8 July 2005, Published online: 6 October 2005     

A solution of the DGLAP equation for gluon at low x

We obtain a solution of the DGLAP equation for the gluon at low x first by expanding the gluon in a Taylor series and then using the method of characteristics. We test its validity by comparing it with that of Glück, Reya and Vogt. The convergence criteria of the approximation used are also discussed. We also calculate εF ₂(x,Q)²/ε In Q ² using its approximate relations with the gluon distribution at low x. The predictions are then compared with the HERA data.     

Influence of gravity on surface melting

A chemical network system for use in the study of reaction cascades is described by the nonlinear rate equation , in which (x) is derived from the –G/x of Taylor's expansion of the order parameterx of the thermodynamic potential, Gibbs' functionG(T,P,x), at about the critical pointC(T _C ,P _C ) of the control variables (parameters)T andP. The behavior of the system around the stable pointx=0 is analyzed only with the sign ofk ₁(T, P), becausek ₂(T, P) is always positive. The system is not closed and is affected by physical and chemical changes in the neighbor systems. WhenT andP fluctuate ( ) through the changes andk ₁ passes zero, the system in the steady state becomes instable andk ₁ jumps at the reaction threshold. Then, the products are formed at the number of moleculesx=(|k ₁|/k ₂)^1/2. To describe such a transition,k ₁ giving theG curvature atx=0 is represented phenomenologically by an approximate function, tanh(G/RT), for a metastable state with a relatively smallk ₁(>0). The reaction takes place in a cascade in accordance with a cubic state equation obtained from tanh(G/RT), which describes a jump of the reaction energy G at the reaction threshold.     

Correlator of heavy- quark currents at small q2 in the large-$\beta_0$ limit

The correlator of vector heavy-quark currents at small q ² is considered in the large- limit. The leading IR renormalon ambiguity of the sum of the perturbative series is canceled by the UV renormalon ambiguity of the gluon condensate. The asymptotic behavior of the perturbative series is obtained in a model-independent way, up to a single unknown normalization factor. Gluon-virtuality distribution functions for the perturbative correction are calculated.Received: 3 December 2004, Published online: 21 February 2005     

Analysis of the logarithmic slope of <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> from the Regge gluon density behavior at small <Emphasis Type="Italic">x</Emphasis>

We study the accuracy of the Regge behavior of the gluon distribution function for an approximate relation that is frequently used to extract the logarithmic slopes of the structure function from the gluon distribution at small x. We show that the Regge behavior analysis results are comparable with HERA data and are also better than other methods that expand the gluon density at distinct points of expansion. We also show that for Q ² = 22.4 GeV², the x dependence of the data is well described by gluon shadowing corrections to the GLR-MQ equation. The resulting analytic expression allows us to predict the logarithmic derivative ∂F ₂(x, Q ²)/∂lnQ ² and to compare the results with the H1 data and a QCD analysis fit with the MRST parameterization input.     

Jet production in <Emphasis Type="Italic">pA</Emphasis> and <Emphasis Type="Italic">AA</Emphasis> collisions in the perturbative QCD pomeron model

Inclusive cross-sections for gluon jet production are studied numerically in the perturbative QCD pomeron model for pA and central AA collisions at high energies. Two forms for the inclusive cross-sections, with and without emission from the triple pomeron vertex, are compared. The difference was found to reduce to a numerical factor for momenta below the saturation momentum Q _s. Above Q _s no difference was found at all. For pA collisions the gluon spectrum was found to be at momenta k below Q _s and above it. For central AA collisions it was found to be at momenta k below Q _s and above it. At large k the spectrum goes like , flattening with energy. The multiplicities turned out to be proportional to A ^0.7 for pA collisions and A for central AA collisions with a good precision. In the latter case they are becoming more peaked at the center with the growth of energy. Their absolute values are high and grow rapidly with energy in accordance with the high value of the BFKL intercept.Received: 12 October 2004, Revised: 22 November 2004, Published online: 21 January 2005     

The role of screening corrections in smallx-behaviour of structure functions

We consider the influence of shadowing effects on the proton structure function in the small-x interval. The gluon-gluon shadowing are noticeable in the HERA kinematical region while the screening of the quark component of the structure function effects negligibly the gluon distribution. The only noticeable effect is the decreasing of sea quark densities at small-x. The explicit form of the gluon distribution in the proton depends significantly on the form of used boundary condition atQ ²=Q ₀ ² . We consider also difference of results obtained with the help of Kuraev-Lipatov-Fadin and Altarelli-Parisi evolution equations.     

Nonlinear correction to the longitudinal structure function at small x

We computed the longitudinal proton structure function F_L, using the nonlinear Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (NLDGLAP) evolution equation approach at small x . For the gluon distribution, the nonlinear effects are related to the longitudinal structure function. As the very small-x behavior of the gluon distribution is obtained by solving the Gribov, Levin, Ryskin, Mueller and Qiu (GLR-MQ) evolution equation with the nonlinear shadowing term incorporated, we show that the strong rise that corresponds to the linear QCD evolution equations can be tamed by screening effects. Consequently, the obtained longitudinal structure function shows a tamed growth at small x . We computed the predictions for all details of the nonlinear longitudinal structure function in the kinematic range where it has been measured by the H1 Collaboration and made comparisons with the computation by Moch, Vermaseren and Vogt at the second order with input data from the MRST QCD fit.     

Properties of a pseudo-Hermitian Hamiltonian for harmonic oscillator decorated with Dirac delta interactions

We have investigated bound state solutions of the Schrodinger equation for one-dimensional harmonic oscillator potential together with even number of Dirac delta functions. These point interactions are located at symmetric points x = x _i and x = −x _i (i = 1, 2,..., N) and they have complex conjugate strengths and , respectively. We present explicit forms of eigenfunctions and an algebraic eigenvalue equation and numerical solutions for this -symmetric Hamiltonian. Presented at the 3rd International Workshop “Pseudo-Hermitian Hamiltonians in Quantum Physics”, Istanbul, Turkey, June 20–22, 2005.     

Differential equations in the spectral parameter

We determine all the potentialsV(x) for the Schrödinger equation (– _x ² +V(x))=k² such that some family of eigenfunctions satisfies a differential equation in the spectral parameterk of the formB(k, _k )ø=(x)ø. For each suchV(x) we determine the algebra of all possible operatorsB and the corresponding functions (x)This research was partially supported by NSF grant DMS 84-03232 and ONR contract NOOO14-84-C-0159     

A phenomenological analysis of the longitudinal structure function at small <Emphasis Type="Italic">x</Emphasis> and low <Emphasis Type="Italic">Q</Emphasis><Superscript>2</Superscript>

The longitudinal structure function in deep-inelastic scattering is one of the observables from which the gluon distribution can be unfolded. Consequently, this observable can be used to constrain the QCD dynamics at small x. In this work we compare the predictions of distinct QCD models with the recent experimental results for F _L(x,Q ²) at small x and low Q ² obtained by the H1 Collaboration. We focus mainly on the color dipole approach, selecting those models which include saturation effects. Such models are suitable at this kinematical region and also resum a wide class of higher-twist contributions to the observables. Therefore, we investigate the influence of these corrections to F _L in the present region of interest.Received: 23 June 2004, Revised: 13 July 2004, Published online: 14 September 2004     

On the polynomial τ-functions for the KP hierarchy and the bispectral property

We show that the -functions obtained from Schur polynomials lead to wave functions w(x ₁, x ₂, ... ; k) that possess the following bispectral property: There exists a differential operator B{k,_k}, independent of x ₁ , such that B{k,_k}w = {x ₁}w, where {x ₁} is independent of k. This extends for the KP hierarchy some earlier results of J. J. Duistermaat and F. A. Grünbaum for the rational solutions of KdV and of P. Wright for certain rational solutions of the generalized KdV equations.     

The role of gluon-gluon shadowing in smallx-behaviour of structure functions

We consider influence of gluon shadowing effects on the proton structure function in the smallx interval. These effects are noticeable in the HERA kinematical region. The explicit form of the gluon distribution in the proton depend significantly on the form of used boundary condition atQ ²=Q ₀ ² . We consider also difference of results obtained with the help of Kuraev-Lipatov-Fadin and Altarelli-Parisi evolution equations.     

Dependency of control parameters on a rate coefficient giving the potential curvature in a reaction cascade model

Many chemical reactions in vivo are self-controlled by fluxes of chemical energy and matter through biological systems, so the induction of such reactions can be governed by changes in the control parameters of the rate equation. A potential of a system is assumed to be given by Gibbs' functionG(T, P, x), which is continuously differentiable, and the rate equation can be derived from the differential (–G/x) of Taylor's expansion ofG _(T,P)(x) for the order parameterx, which corresponds to the product number, at around the critical pointC(T _C, P_C). The equation is described bydx/dt=(x)–k₁x–k₂x³, andk ₂>0. In this equation,k ₁ andk ₂ are functions of the control parameters, temperatureT and pressureP, andk ₁ is allowed to have a positive or negative values as (T, P). Thenk ₁ is an important factor that decides the induction conditions of the reactions with a phase transition in the steady statex=0. Because bothk ₁ (the transition parameter) andG are the quantity of state, they are given by the total differential, and functions that decideG andk ₁ are related to a mutual inverse function. From the above relation, the rate of change ink ₁ by G, which corresponds to the reaction energy of the system, is uniquely determined by a function ofk ₁, [f(k ₁ ^± )] andf(k ₁ ^± ) is described approximately by ±₁ k ₁ ^± in the transient process thatk ₁ approaches zero, where ₁ implies 1/RT. These results indicate that internal driving forces caused by a stimulus in a system are proportional tok ₁ ^± and that the system is regulated by competition of the forces. an approximate function fork ₁ in the transient process is described by tanh (G/RT) and Arrhenius' law is elucidated from this theory.Decreased January 19, 1992