Light–matter interactions simulated on the world’s fastest supercomputer

In a study published this month in The International Journal of High Performance Computing Applications, a research team led by the University of Tsukuba describes a highly efficient method for simulating light-matter interactions at the atomic scale.

What makes these interactions so difficult to simulate? One reason is that phenomena associated with the interactions encompass many areas of physics, involving both the propagation of light waves and the dynamics of electrons and ions in matter. Another reason is that such phenomena can cover a wide range of length and time scales.

Given the multiphysics and multiscale nature of the problem, light-matter interactions are typically modeled using two separate computational methods. The first is electromagnetic analysis, whereby the electromagnetic fields of the light are studied; the second is a quantum-mechanical calculation of the optical properties of the matter. But these methods assume that the electromagnetic fields are weak and that there is a difference in the length scale.

"Our approach provides a unified and improved way to simulate light-matter interactions," says senior author of the study Professor Kazuhiro Yabana. "We achieve this feat by simultaneously solving three key physics equations: the Maxwell equation for the electromagnetic fields, the time-dependent Kohn-Sham equation for the electrons, and the Newton equation for the ions."

The researchers implemented the method in their in-house software SALMON (Scalable Ab initio Light-Matter simulator for Optics and Nanoscience), and they thoroughly optimized the simulation computer code to maximize its performance. They then tested the code by modeling light-matter interactions in a thin film of amorphous silicon dioxide, composed of more than 10,000 atoms. This simulation was carried out using almost 28,000 nodes of the fastest supercomputer in the world, Fugaku, at the RIKEN Center for Computational Science in Kobe, Japan.

"We found that our code is extremely efficient, achieving the goal of one second per time step of the calculation that is needed for practical applications," says Professor Yabana. "The performance is close to its maximum possible value, set by the bandwidth of the computer memory, and the code has the desirable property of excellent weak scalability."

Although the team simulated light-matter interactions in a thin film in this work, their approach could be used to explore many phenomena in nanoscale optics and photonics.