{"title":"Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes","authors":"Xiaoyu Wang, Michael J. Servis","doi":"arxiv-2309.06400","DOIUrl":null,"url":null,"abstract":"Understanding the impact of extractant functionalization in solvent\nextraction is essential to guide the development of better separations\nprocesses. Traditionally, computational extractant design uses electronic\nstructure calculations to determine the metal binding energy of the lowest\nenergy state. Although highly accurate, this approach does not account for all\nthe relevant physics encountered under experimental conditions, such as\ntemperature effects and ligand flexibility, in addition to approximating\nsolvent-extractant interactions with implicit solvent models. In this study, we\nuse classical MD simulations with an advanced sampling method, metadynamics, to\nmap out extractant molecule conformational free energies in the condensed\nphase. We generate the complete conformational landscape in solution for a\nfamily of bidentate malonamide-based extractants with different\nfunctionalizations of the head group and the side chains. In particular, we\nshow how such alkyl functionalization reshapes the free energy landscape,\naffecting the free energy penalty of organizing the extractant into the\ncis-like metal binding conformation from the trans-like conformation of the\nfree extractant in solution. Specifically, functionalizing alkyl tails to the\ncenter of the head group has a greater influence on increasing molecular\nrigidity and disfavoring the binding conformation than functionalizing side\nchains. These findings are consistent with trends in metal binding energetics\nbased on experimentally reported distribution ratios. This study demonstrates\nthe feasibility of using molecular dynamics simulations with advance sampling\ntechniques to investigate extractant conformational energetics in solution,\nwhich, more broadly, will enable extractant design that accounts for entropic\neffects and explicit solvation.","PeriodicalId":501259,"journal":{"name":"arXiv - PHYS - Atomic and Molecular Clusters","volume":"19 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2023-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - PHYS - Atomic and Molecular Clusters","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2309.06400","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Understanding the impact of extractant functionalization in solvent
extraction is essential to guide the development of better separations
processes. Traditionally, computational extractant design uses electronic
structure calculations to determine the metal binding energy of the lowest
energy state. Although highly accurate, this approach does not account for all
the relevant physics encountered under experimental conditions, such as
temperature effects and ligand flexibility, in addition to approximating
solvent-extractant interactions with implicit solvent models. In this study, we
use classical MD simulations with an advanced sampling method, metadynamics, to
map out extractant molecule conformational free energies in the condensed
phase. We generate the complete conformational landscape in solution for a
family of bidentate malonamide-based extractants with different
functionalizations of the head group and the side chains. In particular, we
show how such alkyl functionalization reshapes the free energy landscape,
affecting the free energy penalty of organizing the extractant into the
cis-like metal binding conformation from the trans-like conformation of the
free extractant in solution. Specifically, functionalizing alkyl tails to the
center of the head group has a greater influence on increasing molecular
rigidity and disfavoring the binding conformation than functionalizing side
chains. These findings are consistent with trends in metal binding energetics
based on experimentally reported distribution ratios. This study demonstrates
the feasibility of using molecular dynamics simulations with advance sampling
techniques to investigate extractant conformational energetics in solution,
which, more broadly, will enable extractant design that accounts for entropic
effects and explicit solvation.