Researchers from Japan successfully estimate the bending energy of disiloxane molecule with ultra-fine molecular simulations run on a supercomputer

Discovering novel materials with desirable properties is the ultimate dream in materials science. This requires an accurate prediction of molecular structures and properties, often accomplished by simulations. However, some apparently simple problems have remained unsolvable. In a new study, a research team tackles one such problem, the bending energy for disiloxane molecule, by using a supercomputer and an ultra-fine simulation method, which correctly predicts the bending energy value, outperforming conventional techniques.

Understanding the structural properties of molecules found in nature or synthesized in the laboratory has always been the bread and butter of materials scientists. But, with advancements in science and technology, the endeavor has become even more ambitious: discovering new materials with highly desirable properties. To accomplish such a feat systematically, materials scientists rely upon sophisticated simulation techniques that incorporate the rules of quantum mechanics, the same rules which govern the molecules themselves.

The simulation-based approach has been remarkably successful, so much so that an entire field of study called materials informatics has been dedicated to it. But there have also been instances of failure. A notable example comes from disiloxane, a silicon (Si)-containing compound consisting of a Si-O-Si bridge with three hydrogen atoms at each end. The structure is simple enough, and yet it has been notoriously difficult to estimate just how much energy is needed to bend the Si-O-Si bridge. Experimental results have been inconsistent and theoretical calculations have yielded widely different values due to the sensitivity of the calculated properties to parameter choices and level of theory.

Fortunately, an international research team led by Dr. Kenta Hongo, Associate Professor at Japan Advanced Institute of Science and Technology, has now managed to solve this problem. In their study published in Physical Chemistry Chemical Physics, the team achieved this feat by using a state-of-the-art simulation technique called "first-principles quantum Monte Carlo method" that finally overcame the difficulties other standard techniques could not surmount.

But does it all just boil down to better simulations? Not quite. "Getting an answer that does not agree with the experimentally known value is, in itself, not surprising. The agreement can improve with more careful, and more expensive, simulations. But with disiloxane, the agreement actually becomes worse with more careful simulations," explains Dr. Hongo. "What our method has achieved, rather, is good results without much dependence on the adjustment parameters, so that we don't need to worry about whether the adjusted values are sufficient."

The team compared the first-principles quantum Monte Carlo approach with other standard techniques, such as "density functional theory" (DFT) calculations and "coupled cluster method with single and double substitutions and noniterative triples" (CCSD(T)), along with empirical measurements from previous studies. The three methods differed mainly in their sensitivity towards the "completeness" of basis sets (a set of functions used to define the quantum wavefunctions).

It turned out that for DFT and CCSD(T), the choice of basis set affected the amplitude as well as positions of zero amplitude for the wavefunctions, while for quantum Monte Carlo, it only affected the zero amplitude positions. This allowed one to tweak the amplitude such that the wavefunction shape approached that of an exact solution. "This self-healing property of the amplitude works well to reduce the basis-set dependence and lower the bias arising from an incomplete basis set in calculating the bending energy barrier," elaborates Dr. Hongo.

While this is an exciting development in itself, Prof. Hongo points out the bigger picture. "Molecular simulations are widely used to design new medicines and catalysts. Getting rid of the fundamental difficulties in using them greatly contributes to the design of such materials. With our powerful supercomputers, the method used in our study could be a standard strategy for overcoming such difficulties," he says.

This is indeed deserving of being called a "quantum leap"!

Figure 1. Disiloxane molecule structure in linear form (left) and bent form (right).
A global team of scientists considered the bending energy of a silicate molecule, disiloxane. Although the system appears to be simple and tractable, the bending energy calculation is actually a difficult problem to solve by the conventional simulation methods.
Image credit: Dr. Kenta Hongo from JAIST.

Reference

Funding information
The study was funded by the US Department of Energy, Office of Science, Basic Energy Sciences, under contract DEAC05-00OR22725, MEXT-KAKENHI (JP16KK0097, JP19H04692, and JP21K03400), FLAGSHIP2020 (project nos. hp190169 and hp190167 at Kcomputer), the Air Force Office of Scientific Research (AFOSRAOARD/FA2386-17-1-4049;FA2386-19-1-4015), JSPS Bilateral Joint Projects (with India DST), the HPCI System Research Project (Project ID: hp210019, hp210131, and jh210045), MEXT-KAKENHI (JP16H06439, JP17K17762, JP19K05029, JP19H05169, and JP21K03400), the Air Force Office of Scientific Research (Award Numbers: FA2386-20-1-4036), the JSPS Overseas Research Fellowships, a Grant-in-Aid for Early-Career Scientists (Grant Number JP21K17752), and a Grant-in-Aid for Scientific Research(C) (Grant Number JP21K03400).

February 4, 2022