Full Minimal Coupling All-Electron Real-Time TDDFT for X-Ray-Matter Interactions

Light behaves in two distinct ways depending on how closely we look at it. When we shine a beam of light on a molecule, the standard way scientists describe the interaction assumes the light wave is so long and gentle that it looks like a flat, uniform sheet pushing on the entire object at once. This "dipole" view works perfectly for visible light, where the waves are hundreds of times larger than the atoms they touch. However, when scientists use X-rays, the waves become incredibly short, shrinking down to the size of the atoms themselves. In this regime, the light is no longer a flat sheet; it is a rippling wave that hits one side of an atom before the other. For decades, computer models used to predict how matter reacts to these intense X-rays have ignored this rippling effect, treating the light as if it were still uniform. This simplification breaks down when the light is short enough to resolve the tiny details of the atomic structure, leaving researchers unable to accurately simulate how materials absorb or scatter high-energy X-rays.

A team of researchers has now built a new, more precise way to run these simulations. They developed a method that keeps the full, rippling shape of the light wave as it interacts with every single electron in a material, from isolated molecules to solid crystals. By doing this, they can see physical effects that were previously invisible to standard models. Their work shows that when X-rays hit matter, the light's momentum—the push it carries as it travels—plays a critical role in how electrons move. This momentum allows the light to trigger specific reactions that are strictly forbidden when the light is treated as a uniform sheet. The researchers tested their new method across a wide range of energies, from soft X-rays to hard X-rays, and found that the extra detail provided by the rippling wave changes the outcome of the interaction in measurable and significant ways.

The researchers implemented this new approach in a powerful computer program called FHI-aims, which is designed to model the behavior of electrons in atoms and solids. Instead of approximating the light as a simple, uniform push, their code treats the light as a vector potential that varies in space, carrying the specific wave pattern of the X-ray beam. This allows the simulation to capture two types of interactions simultaneously: one where the light pushes electrons into new states, and another where the light's energy shifts the electrons' energy levels locally. Crucially, the method also tracks the flow of electrons not just as a total average, but as a pattern that varies across the crystal lattice. This is essential because many of the interesting effects of X-rays, such as the generation of light at twice the frequency of the incoming beam, happen at specific points in the material's internal structure that standard models miss entirely.

To prove their method works, the team ran four distinct tests, each designed to challenge the limits of the new approach. First, they simulated a benzene molecule hit by an X-ray pulse tuned to the carbon K-edge, an energy level where carbon atoms strongly absorb light. They compared their results with a different, highly detailed simulation that had already accounted for the rippling light. Their new method matched the other simulation almost perfectly, showing a small but real correction to the molecule's response that matched the expected size of the effect. This confirmed that their code correctly captures the subtle ways the light's wave shape influences the electrons, even in a relatively small molecule.

Next, they turned to a molecule called thiophene, which contains sulfur, and scanned it with X-rays ranging from 100 electron volts to 5,000 electron volts. In this test, they looked at how the electrons responded to the light's push versus how they responded to the light's energy density. They found that at most energies, the light's energy density created a steady, flat background response. However, when the X-ray energy matched the specific resonance of the sulfur or carbon atoms, the "push" response surged dramatically, overtaking the background. This confirmed that their method could correctly identify when the light's momentum becomes the dominant force in the interaction, a transition that is critical for understanding how materials behave under intense X-ray beams.

The third test moved from molecules to a solid crystal: diamond. The researchers simulated a scenario where a diamond crystal was hit by hard X-rays to generate a second harmonic, a beam of light with exactly twice the frequency of the incoming beam. In a perfectly symmetrical diamond crystal, standard physics says this forward-moving beam should not exist at all if the light is treated as a uniform sheet. However, because the researchers kept the full spatial shape of the light wave, their simulation showed that the beam did appear. The signal was weak, but it was real and scaled correctly with the strength of the light, proving that the rippling nature of the X-ray wave can break the symmetry rules that usually forbid such reactions. This result is significant because it suggests that experiments looking for these signals in the forward direction are actually seeing a bulk property of the material, not just a surface effect.

Finally, the team studied a catalytic surface where carbon monoxide and oxygen atoms sit on a sheet of ruthenium metal. This is a system relevant to chemical reactions, such as those in fuel cells. They simulated the absorption of X-rays at the oxygen K-edge, an energy where oxygen atoms absorb light. They found that while the standard model could predict the main features of the absorption, it missed a specific type of transition where an electron jumps to a higher energy state in a way that requires the light's momentum. Their new method revealed that this "forbidden" transition, which involves the electron moving from an s-orbital to a d-orbital, contributes significantly to the energy absorbed by the material. This effect was particularly strong when the atoms were in a transition state, suggesting that the momentum of the X-ray light plays a vital role in how catalytic surfaces absorb energy during chemical reactions.

The implications of this work extend beyond just getting the numbers right. By retaining the full spatial structure of the light, the researchers have opened the door to simulating processes that were previously impossible to model from first principles. They can now predict how X-rays will interact with complex materials, from the surfaces of catalysts to the interiors of solids, with a level of detail that includes the momentum transfer of the photons. This capability is essential for interpreting the results of modern X-ray experiments, where the light is often so intense and the wavelengths so short that the old approximations no longer hold. The method provides a foundation for understanding nonlinear X-ray phenomena, such as wave mixing and diffraction, where the signal is carried by the specific wavevector of the light rather than just its intensity.

While the new method is powerful, the researchers acknowledge that it is a step forward rather than a final solution. Their current implementation tracks the primary interaction between the light and the electrons but does not yet include the feedback loop where the electrons' movement generates new light fields that then act back on the electrons. This means that in extremely strong fields, the results might still underestimate the true effect. Additionally, the method currently focuses on the leading order of the momentum interaction, which is sufficient for many cases but may need to be expanded to capture higher-order effects in more complex scenarios. Despite these limitations, the work demonstrates that the spatial phase of the X-ray light is not just a minor detail but a physically consequential factor that shapes how matter responds to high-energy radiation.

The success of these simulations across such a wide energy range, from 100 electron volts to 5,000 electron volts, validates the approach as a robust tool for the scientific community. It bridges the gap between the simple models used for visible light and the complex reality of X-ray interactions. By showing that the momentum of the photon is a key player in these interactions, the researchers have provided a clearer picture of the fundamental physics governing light-matter coupling. This clarity will help experimentalists design better tests and interpret their data more accurately, particularly in the emerging field of nonlinear X-ray science where the rules of the game are being rewritten by the intensity and coherence of modern light sources.

Ultimately, this work represents a shift in how we model the interaction between light and matter. It moves away from the convenient but inaccurate assumption that light is a uniform force and embraces the reality that light is a wave with a specific shape and momentum. This shift allows scientists to see the hidden channels of interaction that open up when the wavelength of light matches the scale of the atoms. Whether it is the subtle correction to a molecule's dipole, the emergence of a forbidden signal in a crystal, or the enhanced absorption on a catalytic surface, the inclusion of the light's spatial structure reveals a richer, more complex world of physical phenomena. As X-ray sources continue to become more powerful and precise, having a tool that can accurately simulate these interactions will be crucial for unlocking the secrets of matter at the atomic scale.