#### Abstract

We present a review of the state of the art of our understanding of the intrinsic charm and bottom content of the nucleon. We discuss theoretical calculations, constraints from global analyses, and collider observables sensitive to the intrinsic heavy quark distributions. A particular emphasis is put on the potential of a high energy and high luminosity fixed target experiment using the LHC beams (AFTER@LHC) to search for intrinsic charm.

#### 1. Introduction

The existence of a nonperturbative intrinsic heavy quark component in the nucleon is a rigorous prediction of Quantum Chromodynamics (QCD). An unambiguous experimental confirmation is still missing and would represent a major discovery. The goal of this paper is to summarize our current understanding of this subject with a particular focus on the potential of a high energy and high luminosity fixed target experiment using the LHC beams (AFTER@LHC) [1–4] to search for intrinsic charm.

Production processes sensitive to the intrinsic heavy quark distributions of protons and nuclei are among the most interesting hadronic physics topics that can be investigated with AFTER@LHC. In contrast to the familiar extrinsic contributions which arise from gluon splitting in perturbative QCD, the intrinsic heavy quarks have multiple connections to the valence quarks of the proton and thus are sensitive to its nonperturbative structure. For example, if the gluon-gluon scattering box diagram, (the analog of QED light-by-light scattering), is inserted into the proton self-energy, the cut of this amplitude generates five-quark Fock states of the proton ; see Figure 1.

Intrinsic strange, charm, and bottom quarks are thus a fundamental property of the wavefunctions of hadronic bound states [5–8]. While the extrinsic contributions to the heavy quark parton distribution functions (PDFs) are most important at low and depend logarithmically on the heavy quark mass , the intrinsic heavy quark contributions are dominant at high and depend on . Because the extrinsic heavy quarks are generated by gluon splitting, their PDFs are always softer than those of the parent gluon by a factor of . In contrast, the high intrinsic heavy quark contributions are kinematically dominated by the regime where the state is minimally off shell, corresponding to equal rapidities of the constituent quarks. The resulting momentum and spin distributions of the intrinsic and can be distinct, for example, , since the comoving quarks are sensitive to the global quantum numbers of the protons.

A finite intrinsic charm contribution to the nucleon has been extracted from lattice QCD. An analysis by the MILC collaboration [9] yields a probability for the charm matrix element in the range of 5-6%, consistent with a four-loop perturbative QCD calculation [10].

While the first experimental evidence of intrinsic heavy quarks came from the EMC measurement of the large charm structure function [11], a variety of other charm hadrons and charmonium measurements are consistent with the existence of intrinsic charm. Open charm observables in hadroproduction include forward production at the ISR [12]. Similarly, the coalescence of comoving , , and quarks from the intrinsic bottom Fock state in the proton can explain the high production of the baryon, as observed at the ISR [12]. and asymmetries between leading and nonleading charms ( mesons which share valence quarks with the projectile and mesons which do not, resp.) measured as functions of and in fixed target experiments, WA89 and WA82 at CERN; E791 and SELEX at Fermilab; see [13–15] and references therein. Previous fixed target measurements also give indications of important intrinsic charm contributions, particularly from the nuclear mass, or , dependence, as measured by NA3 at CERN as well as E772 and, later, E866 at Fermilab; see, for example, [16]. Indeed, the dependence, proportional to , is quite different than the expected from extrinsic-type production [17]. At large , there are indications of dependence, consistent with a nuclear surface-type interaction instead of the volume dependence of pQCD. In addition, the NA3 collaboration measured double production at forward in interactions, difficult to explain without an intrinsic charm mechanism [18]. All of these observables can be studied with higher energies and luminosities at AFTER@LHC, making precision measurements possible for the first time.

In addition to the typical observables for intrinsic heavy quarks, these intrinsic heavy quarks also contribute to a number of more exotic observables and inclusive and diffractive Higgs production , in which the Higgs boson carries a significant fraction of the projectile proton momentum [19, 20]. There are also important implications for intrinsic charm and bottom quarks in Standard Model physics, as in the weak decays of the -meson [21] and a novel solution to the problem [22]. AFTER@LHC could also shed light on these topics.

The rest of this paper is organized as follows. In Section 2, we give an overview of the theoretical models predicting the -shape (but not the normalization) of the intrinsic charm and bottom parton distribution functions. In Section 3, we discuss the constraints on the normalization of the intrinsic charm (IC) obtained in global analyses of PDFs. Section 4 is devoted to the intrinsic bottom (IB) content of the nucleon, for which there are currently no quantitative constraints. In Section 5 we review collider observables sensitive to an intrinsic charm or bottom PDF. Finally, in Section 6 we present our conclusions.

#### 2. Theoretical Models

The QCD wavefunction of a hadron can be represented as a superposition of quark and gluon Fock states. For example, at fixed light-front time, a hadron wavefunction can be expanded as a sum over the complete basis of free quark and gluon states: where the color-singlet states, , represent the fluctuations in the hadron wavefunction with the Fock components , , , and so forth. The boost-invariant light-front wavefunctions, , are functions of the relative momentum coordinates and where denotes the parton momenta and the hadron momentum. Momentum conservation demands and where is the number of partons in state . For example, as predicted by Brodsky and collaborators, in the BHPS model intrinsic charm fluctuations [5, 23] can be liberated by a soft interaction which breaks the coherence of the Fock state [24] provided the system is probed during the characteristic time that such fluctuations exist.

Microscopically, the intrinsic heavy quark Fock component in the proton wavefunction, , is generated by virtual interactions such as where the gluons couple to two or more valence quarks. The probability for fluctuations to exist in a hadron is higher twist since it scales as relative to the extrinsic, EC, leading-twist production by photon-gluon fusion [18].

The dominant Fock state configurations are not far off shell and thus have minimal invariant mass, , where is the square of the average transverse mass of parton . The general form of the Fock state wavefunction for a hadron with mass appropriate to any frame at fixed light-front time is where is a vertex function, expected to be a slowly varying, decreasing function of . The particle distributions are then controlled by the light-front energy denominator and phase space. This form for the higher Fock components is applicable to an arbitrary number of light and heavy partons. Intrinsic Fock components with minimum invariant mass correspond to configurations with equal rapidity constituents. Thus, unlike extrinsic heavy quarks generated from a single parton, intrinsic heavy quarks carry a larger fraction of the parent momentum than the light quarks in the state [5, 23].

The parton distributions reflect the underlying shape of the Fock state wavefunction. Assuming that it is sufficient to use for the transverse momentum, the probability distribution as a function of in a general -particle intrinsic Fock state is where normalizes the -particle Fock state probability.

At LO in the heavy quark limit, , , ,leading to

There are many applications of intrinsic charm in charm hadron production. See, for example, [13–16, 18] for more details.

Paiva et al. have also calculated an intrinsic charm component of the nucleon sea within the context of the meson cloud model [25]. They assumed that the nucleon can fluctuate into . The distribution in the nucleon is thenwherewith being a form factor at the vertex and . In this case they chose a monopole form factor with GeV. The coupling constant was assumed to be . From heavy quark effective theories [26], the distribution in the is expected to be hard because, in the bound state, the exchanges momenta much less than . They make the extreme assumption that the entire momentum is carried by the charm quark, .

Next, Steffens et al. investigated all the charm structure function data with two variants of intrinsic charm [27]. The first was that of (4), called IC1 in their paper, while the second was a meson cloud model, IC2. In the second approach, the distribution is obtained from the light-front distribution of mesons in the nucleon:A hard charm momentum distribution was assumed in the , similar to that of [25]. The vertex function is parameterized as where is the square of the center-of-mass energy of the system and the coupling constant at . For an intrinsic charm probability of 1%, GeV. The charm distribution is then where the charm distribution in the is assumed to be and .

Pumplin [28] considered a model where a point scalar particle of mass couples with strength to scalar particles with mass . The probability density is thenwhere . The form factor suppresses higher mass state contributions. If the quark transverse momenta are neglected, with much greater than all other mass scales, and , then the BHPS model is recovered. Two types of form factors were studied, an exponential and a power law where the cutoff is varied between 2 and 10 GeV.

Hobbs et al. employed a meson cloud type approach but specified the spin and parity of all lowest mass charm meson-baryon combinations from the 5-particle Fock states of the proton [29]. They pointed out that treating quarks as scalar point-like particles, as in, for example, [28], does not conserve spin and parity. They calculated the appropriate meson-baryon splitting functions for the meson-baryon combinations and found that the production of charm mesons would be almost entirely through mesons. To study the phenomenological distributions of charm mesons and baryons in this approach, they studied exponential and confining vertex functions, and , respectively. They used these results to compare to the distribution from the ISR [30] and the asymmetry from SELEX [31]. See [29] for details.

#### 3. Global Analyses of PDFs with Intrinsic Charm

In the standard approach employed by almost all global analyses of PDFs, the heavy quark distributions are generated* radiatively*, according to DGLAP evolution equations [32–34], starting with a perturbatively calculable boundary condition [35, 36] at a scale of the order of the heavy quark mass. In other words, there are no free fit parameters associated with the heavy quark distribution and it is entirely related to the gluon distribution function at the scale of the boundary condition. As a consequence, also the PDF uncertainties for the heavy quark and the gluon PDFs are strongly correlated as has been discussed in the context of inclusive Higgs production at the Tevatron and the LHC [37]. However, a purely perturbative treatment might not be adequate in particular for the charm quark with a mass GeV which is not much bigger than typical hadronic scales but also for the bottom quark with a mass GeV. Indeed, as discussed above, light-front models predict a nonperturbative (“intrinsic”) heavy quark component in the proton wavefunction [5, 23]. Motivated by the theoretical predictions of the BHPS light-front model, analyses of the charm distribution in the proton going beyond the common assumption of purely radiatively generated charm date back almost as far as the BHPS predictions themselves. For definiteness, in the following we refer to the radiatively generated charm by and to the intrinsic charm by . The full charm parton distribution is then given by the sum . Strictly speaking, this decomposition is defined at the initial scale of the DGLAP evolution but holds to a good approximation at any scale since the intrinsic component is governed (to a very good approximation) by a standalone nonsinglet evolution equation [38]. A similar decomposition is understood for the bottom quark which will be discussed in Section 4.

The BHPS model of the Fock state predicts a simple form for :If there is a 1% intrinsic charm contribution to the proton PDF, .

Hoffmann and Moore incorporated mass effects and introduced next-to-leading order corrections as well as scale evolution [39]. They compared their result to the EMC data from muon scattering on iron at high and with the intrinsic charm contribution added to the leading order calculation of by photon-gluon fusion.

A complete next-to-leading order analysis of both the “extrinsic” radiatively generated charm component and the intrinsic component was later carried out by Harris et al. [7]. The EMC data with , 95, and 168 GeV were fit by a sum of the extrinsic and intrinsic components [7]. The normalization of the two components was left as free parameters:with the scale . The parameter , typically larger than unity, was considered to be an estimate of the NNLO contribution to the extrinsic contribution. Since a 1% normalization of the IC component was assumed in (11), the fitted value of is the fraction of this normalization. Given the quality of the data, no statement could be made about the intrinsic charm content of the proton when and 95 GeV. However, with GeV an intrinsic charm contribution of ()% was indicated. These results were consistent with those of the original analysis by Hoffmann and Moore [39].

The BHPS light-front model assumes that . Meson cloud models, introduced later, treat the 5-particle Fock state as a combination of (predominantly) . In this case, of course, with the quark in the carrying more momentum than the quark in the charm baryon. An analysis by Steffens et al. in the context of the meson cloud model and using a hybrid scheme to interpolate between massless evolution at high and “extrinsic” production at low found a limit of ~0.4% [27].

Regardless of whether or not the models predict , intrinsic charm should provide the dominant contribution to the charm density in the proton at large [28].

For some time, no other analyses of the charm structure function were made. The EMC data remain the only measurement of the charm structure function in the relevant regime and are the only DIS data cited as evidence for intrinsic charm. The HERA data on were at too low to address the issue.

The first global analyses of the proton PDFs with an intrinsic charm contribution included were performed by members of the CTEQ collaboration [40, 41]. In addition to the BHPS and meson cloud approaches, they also allowed for a “sea-like” contribution with the same shape as the radiatively generated charm distribution. They characterized the magnitude of the intrinsic charm component () by the first moment of the charm distribution at the input scale GeV (Note that at the radiatively generated charm component () vanishes at NLO in the scheme so that .):which translates into a momentum fractionThey found that the global analyses of hard-scattering data provided no evidence for or against the existence of intrinsic charm up to ; that is, the quality of the fit is insensitive to in this interval. They also found that the allowed range was greatest for the sea-like IC expected since this shape is rather easily interchangeable with other sea quark components while the other, harder, charm distributions are not [40]. In addition, they concluded that the enhancement due to IC relative to analyses without it persisted up to scales of ~100 GeV and could have an influence on charm-initiated processes at the LHC, as is discussed later. The CTEQ6.6C proton PDFs were generated as a result of this analysis [41].

There are two recent updates to the global analyses, reaching different conclusions about the importance of intrinsic charm. The first, by Dulat et al. [42], follows the previous work in the context of the CTEQ collaboration [40, 41]. The second, by Jimenez-Delgado et al. [43], included more lower energy data than the previous global analyses.

The result of Dulat et al. [42] was based on the CT10 NNLO parton densities. Here the strong coupling, , the evolution equations, and the matrix elements are calculated at NNLO. Only the inclusive jet data still required NLO expressions. Their analysis included DIS data from BCDMS, NMC, CDHSW, and CCFR; SIDIS data from NuTeV and CCFR; the combined DIS and data from HERA; Drell-Yan production; the charge asymmetry and rapidity from CDF and D0; and the inclusive jet measurements from CDF and D0; see [42] for a complete list.

Two models of IC were considered: the BHPS light-front model and the sea-like IC introduced in [40]. They found a broader possible probability range for IC in this analysis, for BHPS and for the sea-like IC, summarized in Figure 2. This finding differs from the previous work which found a larger upper limit on IC for the sea-like model. They believe that the difference is caused by the improved treatment of the charm quark mass in the later study [42].

In addition to the global fit, they also tested the sensitivity of their result to individual experiments by introducing a penalty factor, , for each experiment . This penalty factor is designed to increase more rapidly than the for that experiment when goes beyond the 90% confidence level. The penalty factor employs an equivalent Gaussian variable which measures the goodness of fit for each individual data set. Values of are considered good fits, is considered to be a poor fit, and values of are better fits than expected from usual statistical analyses. Using the dependence on , they determined which of the data sets used in the global analyses are most sensitive to intrinsic charm. The upper limit on the BHPS value of comes from the CCFR structure function data while the HERA combined charm data sets the upper limit on IC from the sea-like model [42].

They also studied the sensitivity of their sea-like result to the charm quark mass and found that if the charm quark mass was raised from 1.3 GeV, as in the CT10 fits, to 1.67 GeV, then the minimum for the global analyses would support rather than 0 although the global is worse for the larger charm mass [42]. Finally, they showed how and production at the LHC might be affected by a nonzero IC contribution.

In the most recent study, Jimenez-Delgado et al. [43] included the full range of high energy scattering data by using looser kinematic cuts and . In particular, they included the lower energy SLAC fixed target data which did not pass the more stringent standard DIS cuts on the plane applied in the previous work [40–42]. The EMC data, cited as the strongest evidence for intrinsic charm in DIS, are used as a consistency check. The low energy, high , fixed target data lie precisely in the region where IC is expected to be most important. Thus including these data could enhance the sensitivity of the global fit to IC. Note, however, that some of these newly added data are on heavier targets than the deuteron and thus target mass corrections, nuclear corrections for , and higher-twist effects need to be taken into account [43].

They followed the framework of the JR14 [44] global fit which decomposed into light and heavy components. The charm component is itself separated into the “extrinsic” and intrinsic charm components. The fixed-flavor number scheme is used to compute the extrinsic contribution. In this scheme, the charm quark mass enters the PDF evolution only indirectly through the running of [43]. They employed a charm quark mass of 1.3 GeV, as did Dulat et al. [42]. They used all three intrinsic charm models previously considered: BHPS, the meson cloud model (this time including pseudoscalar and vector mesons as well as spin and spin charm baryons—the CTEQ analyses only included the scalar fluctuation), and the sea-like component [43]. The IC contribution was evolved up to NLO.

They found that the total is minimized for with % at the level. When a hadron suppression factor to suppress charm contributions near threshold is applied, they find a minimum at % for the full data set. The SLAC (large ), NMC cross sections (medium ), and HERA (low ) display the greatest sensitivity to IC; see Figure 3 for details. However, fits without the SLAC data still give a low IC contribution [43]. The difference between their results and previous results is in part due to the very different tolerance criteria, for their fit and for Dulat et al. [42]. Increasing the tolerance to would also accommodate % at the level [43]. (For a critical discussion of the analysis in [44] and in particular of the tolerance criterion see [45, 46].)

When checked against the EMC data, a clear preference for IC is found, as expected, for the highest data. (The EMC - data have not been included in some global analyses because of a possible conflict with HERA - data at very low ; however, the relative suppression of the low and low EMC data could be accounted for by nuclear shadowing, suppressed evolution, higher twist, and other effects.)

Given that the two most recent analyses set significantly different limits on IC, it is important to collect further large data, particularly on to try and place greater confidence on the limit of IC in the nucleon. This would be an important measurement at the future electron-ion collider.

#### 4. Predictions for Intrinsic Bottom

In contrast to the case of intrinsic charm, there is currently no global analysis available that investigates the possibility of an intrinsic bottom (IB) content of the nucleon. The main reason for this is the lack of experimental data that could constrain it. The BHPS light-front model [5] predicts the existence of IB with an -shape very similar to the one of IC given in (10) but with a normalization which is parametrically suppressed by the ratio . This fact, together with the observation that the IB PDF is governed (to an excellent approximation) by an independent nonsinglet evolution equation [38], can be used to investigate IB in a flexible way without the need of a dedicated global analysis. Such a study has been done in [38] where a set of decoupled IB (and IC) PDFs has been provided and used together with the CTEQ6.6 PDFs [41] to estimate the impact of the IB on new physics searches at the LHC. The advantage of this approach is that the provided IB (IC) PDF can be used with any standard set of PDFs and the normalization of the intrinsic component can be freely adjusted. This is especially useful for studies of possible IB effects, as, in that case, there are no experimental limits on what amount of IB is allowed.

In the following we show some of the results found in [38]. In this work, the boundary condition for the IB distribution was modeled using the IC distributions in the CTEQ analyses [40, 41] scaled down by the mass factor . The result of such an intrinsic bottom distribution , with normalization , is shown in Figure 4, where the ratio of the intrinsic () and the radiatively generated () component of the bottom PDF is plotted. As always in the light-front models the intrinsic component is mostly present at large values. We can see that for low scales GeV the modification of the bottom PDF, , can reach . However, it decreases rapidly with the rising scale. Since evolves independently of the other PDFs the change in the normalization of the IB component in Figure 4 can be done by simply rescaling the curves in the figure. If we allowed for a normalization of the IB the modification of the bottom PDF would be given by , which for high and GeV would result in an enhancement of the bottom PDF by a factor ~6.25. However, at a scale of around 100 GeV and below 0.2–0.3, even with the higher IB normalization, the effect is becoming negligible.

In Figure 5 we show the sum of the intrinsic bottom PDF and the dynamically generated PDF from CTEQ6.6 for different normalizations of the IB component, namely, 0.01 and . We compare this sum to the asymmetric uncertainties (The asymmetric errors are computed following [47, 48].) of the CTEQ6.6 PDF set (a). In the same figure is also shown the ratio of the same PDFs to the central value of CTEQ6.6 (b). As can be seen, the IB curve with the normalization clearly lies outside the uncertainty band whereas the one with the smaller normalization is marginally outside the band (up to ).

**(a)**

**(b)**

If we are looking for new physics with couplings proportional to the mass, the suppression of IB compared to the IC would be partly compensated by the square of the coupling. For a more detailed study of the relevant parton-parton luminosities please see [38].

#### 5. Collider Observables

Several collider observables receive large contributions from heavy quark initiated subprocesses and are hence potentially sensitive to an intrinsic charm content in the nucleon. In order to expect optimal effects the heavy quark PDF should be probed at large (for light-front models) and not too large factorization scales. This kinematic region is best accessible at lower energies in the center-of-mass system (cms) and/or large rapidities. Therefore, a fixed target experiment like AFTER@LHC [1–4] operating at cms energy GeV with a high luminosity is ideally suited for searches of IC effects. In the following we review some of the collider processes which have been studied in the literature in this respect.

##### 5.1. Open Heavy Flavor Production

Inclusive charm hadron () production in hadronic collisions was advocated in [49] as a laboratory to probe IC inside the colliding hadrons. In this analysis, predictions for the differential cross section in dependence on the transverse momentum were obtained in the general-mass variable-flavor-number scheme (GM-VFNS) [50–52] at next-to-leading order (NLO). In this scheme, the charm quark is an active parton and the differential cross sections of inclusive charm meson production depend heavily on the PDF of the charm quark. The sensitivity of these cross sections to IC was studied for the Tevatron at cms energy of 1960 GeV and the Relativistic Heavy Ion Collider (RHIC) at cms energies of 200 GeV (RHIC200) and 500 GeV (RHIC500). The different IC models from the CTEQ6.5c global analysis [40] were employed together with the fragmentation functions for charm mesons from [53]. While the effects at the Tevatron were found to be very moderate and likely not testable, large enhancements were found at RHIC200 reaching values of ~3 at GeV. Unfortunately, the measurements at RHIC200 are limited by the luminosity. At RHIC500 the cross section is increased by about a factor 3.6. However, the sensitivity to IC for the light-front models is greatly reduced.

More recently, the GM-VFNS was applied to obtain predictions for the production of inclusive mesons at the LHC for cms energy of 7 TeV (LHC7) [54]. It was found that the production cross sections at large rapidities are sensitive to an IC component. These predictions can be tested by measurements at forward rapidities with the LHCb detector.

The ideal experiment to search for the effects of IC would be a high luminosity fixed target experiment such as AFTER@LHC operating at cms energy of 115 GeV. In Figure 6 we show results for inclusive meson production as a function of the transverse momentum of the meson integrated over the rapidity range (in the laboratory frame) in essentially the same setup as in [49] to which we refer for details. The only difference is that, following [55], the default choice for the renormalization and factorization scales is , , where is the transverse mass. The theoretical predictions are shown on an absolute scale in Figure 6(a) and as a ratio with respect to the default results in Figure 6(b). In both figures, the black dotted lines have been obtained by varying the renormalization scale around the central choice to (upper line) and (lower line). In Figure 6(b) we repeat the calculation of the central prediction in turn with PDF sets CTEQ6.5Cn for and normalize the outcome to the default prediction with zero IC of Figure 6(a). We observe that the ratios for corresponding to the BHPS () or meson cloud () models become very large at large . Indeed, the default cross section can be increased by more than a factor 5 at GeV in scenarios with maximally allowed intrinsic charm (). Even for the IC sets with smaller normalization () corresponding to and the cross section would be enhanced by a factor larger than 2 (red solid line) or 3 (blue dashed line) at GeV. It is also interesting to note that the phenomenological models for a sea-like IC () lead to a significant enhancement of the cross section at small which would be probed at AFTER@LHC as well.

**(a)**

**(b)**

##### 5.2. Production of a Photon in Association with a Charm Quark

Another process with a wide range of phenomenological applications in , , and collisions [56–58] which is very sensitive to the heavy quark PDF is the associated production of a photon with a heavy quark. A dedicated study of this process at the LHC operating at TeV (LHC8) was performed in [59, 60] where it was demonstrated that the existence of IC in the proton can be visible at large transverse momenta of the photons and heavy quark jets at rapidities , . Indeed, for the BHPS model the cross section can be enhanced by a factor of 2-3 for GeV (see Figure 5 in [60]). This comes with the penalty that the cross section falls rapidly with increasing transverse momentum so that this measurement will be limited by statistics.

Again, as for open heavy flavor production, the lower cms energy together with the high luminosity makes a fixed target experiment like AFTER@LHC the ideal place to discover IC using production. This can be seen in Figure 7, where the differential cross section is enhanced by a factor 5 at GeV (b) with a not too small cross section (a).

**(a)**

**(b)**

##### 5.3. Vector Boson Production

Dulat et al. [42] studied the sensitivity of and production to the presence of IC. Vector boson production at the LHC is an interesting testing ground for IC because they are produced at relatively large and is a rather clean final state. They did a NNLO calculation of and production including IC based on their global fits at and 14 TeV. They also studied the ratio relative to the result with no IC. Neither of these calculations showed an effect larger than the uncertainties due to the CT10 sets themselves. However, when the distribution with IC was compared to that without, they saw a factor of two enhancement at GeV for TeV in the range . The corresponding enhancement at 14 TeV was smaller at the same because the value reached is reduced at the higher energy [42].

We show a simple test case here for and production to NLO at TeV. We use only the BHPS IC parameterization for the five-particle Fock state, shown in (10). We assume a 1% normalization and no evolution to maximize the possible effect at forward rapidity. The -integrated rapidity distribution is shown in Figure 8, as is the ratio of the result with IC to that without as a function of rapidity. The rapidity distributions without IC are given by the solid curves while the dashed curves are the calculations with the BHPS IC contribution to the charm parton density. With BHPS IC, one expects enhancement only at forward rapidity. The enhancement from IC appears for . Note that if the sea-like IC would be used instead, the enhancement would be small but finite over all rapidity.

**(a)**

**(b)**

The cross section is largest and most forward peaked, because of the contribution. The contribution from the part is a very small addition since the valence contribution is large and peaks at large , making the distribution larger at than at . Indeed, it gives the smallest IC contribution. The distribution should have the largest possible contribution from IC because both and peak at low and because the valence distribution peaks at lower so that the rapidity distribution has a maximum at . At , the IC enhancement is ~40%. Finally, the distribution, with a plateau over , also has a very small IC contribution because the charm enhancement only comes through .

Such IC enhancements are only visible outside the midrapidity acceptance of the collider detector coverage of CMS and ATLAS. However, LHCb or ALICE covers this forward rapidity range with muons and could detect forward . They could also look at the lepton rapidity asymmetry at forward rapidity. The statistical accuracy of the measurement would need to be high to distinguish an IC enhancement from the no IC result, especially since the 1% BHPS IC is likely an upper limit on this enhancement. Note that the higher energy of LHC Run 2 will reduce the potential enhancement even though it would increase the rates.

#### 6. Conclusions

The existence of nonperturbative intrinsic charm and bottom components is a fundamental prediction of QCD. In this paper, we have reviewed the current status of our understanding of this intrinsic heavy quark content of the nucleon which yet remains to be confirmed experimentally. In particular, after introducing theoretical models predicting the intrinsic heavy quark distributions we have turned to a summary of the available information on intrinsic charm coming from global analyses of parton distribution functions. There are no global analyses of intrinsic bottom available and we have described how IB can be modeled in order to explore its impact on collider observables keeping in mind that bottom quark initiated subprocesses play an important role in certain electroweak observables and in models for physics beyond the Standard Model. We then have turned to a discussion of collider processes where IC could be discovered. Generally, the effects of IC are larger at colliders with a lower center-of-mass energy and for hard processes with moderate factorization scales. Therefore, a high luminosity fixed target experiment like AFTER@LHC operating at a center-of-mass energy GeV would be ideally suited to discover or constrain IC.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors are grateful to T. Stavreva for providing Figure 7. The work of S. J. Brodsky was supported by the Department of Energy Contract no. DE-AC02-76SF00515. The work of R. Vogt was performed under the auspices of the US Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.