Double Parton Distributions in the nucleon: Lattice parameter dependence

Inside every proton, the building block of atomic nuclei, there is a chaotic, seething sea of smaller particles called quarks and gluons. For decades, physicists have understood how these particles behave when they interact one at a time. However, the universe is rarely so simple. In the high-energy collisions that occur inside particle accelerators, it is possible for two separate pairs of quarks to collide simultaneously within a single proton-proton crash. These events, known as double parton interactions, are crucial for understanding the background noise in experiments searching for new, unknown physics. If scientists cannot precisely predict how these double collisions happen, they might mistake a standard background event for a discovery of new particles. To solve this, researchers need a map of how two quarks are arranged inside a proton at the same time, a map that describes their positions and momenta in a way that goes far beyond our current knowledge.

The challenge lies in the fact that these double interactions are incredibly difficult to isolate in a real-world experiment. The signals are faint, and the data is messy. To get a clearer picture, a team of researchers turned to the most powerful tool available for studying the strong force that binds quarks together: lattice quantum chromodynamics. This method treats space and time not as a smooth continuum, but as a grid of points, allowing scientists to simulate the behavior of protons on a supercomputer. By calculating the probability of finding two quarks at specific distances from each other, they can build the theoretical foundation needed to interpret future data from the Large Hadron Collider and other facilities. The goal is to move from rough approximations to precise predictions, ensuring that when new physics is found, it is not hidden by a misunderstanding of the old kind.

In this specific study, the researchers set out to test the reliability of their simulations by checking how sensitive their results are to the settings of their computer models. In their previous work, they had relied on a single set of simulation parameters, which included a grid size that was not yet fine enough to perfectly mimic reality and a proton mass that was heavier than the one found in nature. To ensure their map of double parton distributions was trustworthy, they needed to see if changing these settings would alter the outcome. They expanded their work to include six different simulation setups, varying the spacing between the grid points and the masses of the quarks inside the proton. This allowed them to distinguish between genuine physical effects and artifacts that might simply be errors introduced by the limitations of the computer grid.

The team performed these calculations using a specific type of mathematical framework known as the Wilson-Clover fermion action, which is a standard way to represent quarks on a lattice. They generated data using ensembles provided by the CLS collaboration, a group dedicated to creating high-quality simulation environments. By running simulations on grids with different spacings, ranging from 0.064 femtometers to 0.085 femtometers, and with pion masses varying between 214 MeV and 422 MeV, they could observe how the results shifted. They focused on the "invariant functions," which are the core mathematical quantities that describe the relationship between the two quarks, stripping away the noise of the simulation to find the underlying physical truth.

One of the primary concerns was whether the size of the simulated universe would distort the results. Because the computer grid is finite, particles can interact with their own "mirror images" across the boundary, creating a false signal. The researchers analyzed the data by looking at the direction of the separation between the two quarks. They found that when they excluded data points where the separation was aligned too closely with the edges of the grid, the mirror effects disappeared. This confirmed that for larger distances, the finite size of the simulation box was no longer a problem, and the results were stable.

They also investigated whether the spacing of the grid itself was introducing errors. In some channels, particularly those involving specific combinations of quark flavors, they observed that the results did change slightly as the grid became finer. However, these changes were small and mostly confined to very short distances. For the most important channels, which involve vector currents, the results from the coarsest and finest grids agreed with each other within the statistical uncertainty. This suggests that the current grid sizes are already sufficient to capture the essential physics, and that the simulations are converging toward the correct answer as the grid becomes finer.

The study also examined how the mass of the quarks influenced the double parton distributions. They compared simulations with heavier quarks to those with lighter quarks, moving closer to the physical mass found in nature. They discovered a clear trend: as the light quark mass decreased, the strength of the signal in certain channels also decreased. This dependence was particularly noticeable in combinations involving up and down quarks. Interestingly, they found that the strange quark mass, which is heavier, had a much smaller effect, though there were hints of a slight influence in specific channels. The fact that the signal remained strong even at the lightest quark masses tested, down to 214 MeV, confirmed that these double-parton effects are robust and can be studied even when the simulation is close to the physical reality of our universe.

The researchers concluded that their method is stable and reliable. The double parton distributions they calculated do not change wildly with the simulation settings, which gives confidence that the results reflect real physics rather than computer artifacts. While some channels showed a mild dependence on the grid spacing or quark mass, the overall picture is consistent. The work provides a solid foundation for future studies, showing that lattice simulations can now be used to constrain theoretical models of double parton interactions with a high degree of precision. This progress is a vital step toward understanding the complex inner workings of the proton, ensuring that when the next generation of particle accelerators comes online, physicists will have the tools they need to separate the known from the unknown.