Pairing a Global Optimization Algorithm with EXAFS to Accelerate Prediction of Lanthanide Structures in Solution

Thomas Summers,Difan Zhang,Josiane Sobrinho,Ana de Bettencourt-Dias,Roger Rousseau,Vassiliki-Alexandra Glezakou,David Cantu
DOI: https://doi.org/10.26434/chemrxiv-2024-lgkjh
2024-09-16
Abstract:Ensemble-average sampling of structures sampled from ab initio molecular dynamics (AIMD) simulations can be used to predict theoretical extended X-ray absorption fine structure (EXAFS) signals that closely match experimental spectra. However, AIMD simulations are time-consuming and resource-intensive, particularly for solvated lanthanide ions, which often form multiple non-rigid geometries with high coordination numbers. To accelerate the characterization of lanthanide structures in solution, we employed the Northwest Potential Energy Surface Search Engine (NWPEsSe), an adaptive-learning global optimization algorithm, to efficiently screen first-shell structures. As case studies, we examine two systems: Eu(NO3)3 dissolved in acetonitrile with a terpyridine ligand (terpyNO2), and Nd(NO3)3 dissolved in acetonitrile. The theoretical spectra for structures identified by NWPEsSe were compared to both experimental and AIMD-derived EXAFS spectra. The NWPEsSe algorithm successfully identified the proper solvation structure for both Eu(NO3)3(terpyNO2) and Nd(NO3)(acetonitrile)3, with the calculated EXAFS signals closely matching the experimental spectra for the Eu-ligand complex and showing good similarity for the Nd salt; the better agreement with the ligand-containing structure is attributed to a less dynamic coordination environment due to the rigid ligand. The key advantage of the global optimization algorithm lies in its ability to sample the coordination environment across the potential energy surface. The structure identification process can be done with a variety of methods and generally has reduced the time required from a month to within a week.
Chemistry
What problem does this paper attempt to address?
The paper aims to address the issue of structural characterization of lanthanide ions in solution. Specifically, the paper attempts to accelerate the prediction of the structure of lanthanide ions in solution by combining a global optimization algorithm (Northwest Potential Energy Surface Search Engine, NWPEsSe) with extended X-ray absorption fine structure (EXAFS) spectroscopy. The paper discusses two systems as case studies: 1. A complex formed by Eu(NO₃)₃ dissolved in acetonitrile with the terpyNO₂ ligand (Eu(terpyNO₂)(NO₃)₃). 2. Nd(NO₃)₃ salt dissolved in acetonitrile. The theoretical spectra screened by the NWPEsSe algorithm were compared with the experimental spectra and the EXAFS spectra obtained from AIMD simulations. The results show that the NWPEsSe algorithm successfully identified the correct solvation structure, and the calculated EXAFS signals were highly consistent with the experimental data. Moreover, this method significantly improved computational efficiency compared to traditional AIMD simulation methods, reducing the time required from one month to within one week.