Spin-orbit coupling as a probe to decipher halogen bonding

The Chemical Bond at the Bottom of the Periodic Table: The Case of the Heavy Astatine and the Super-Heavy Tennessine

In this work, we study the chemical bond in molecules containing heavy and super-heavy elements according to the current state-of-the-art bonding models. An Energy Decomposition Analysis in combination with Natural Orbital for Chemical Valence (EDA-NOCV) within the relativistic four-component Dirac-Kohn-Sham (DKS) framework is employed, which allows to successfully include the spin-orbit coupling (SOC) effects on the chemical bond description. Simple halogen-bonded adducts ClX⋯L (X=At, Ts; L=NH3, Br-, H2O, CO) of astatine and tennessine have been selected to assess a trend on descending along a group, while modulating the ClX⋯L bond features through the different electronic nature of the ligand L. Interesting effects caused by SOC have been revealed: i) a huge increase of the ClTs dipole moment (which is almost twice as that of ClAt), ii) a lowering of the ClX⋯L bonding energy arising from different contributions to the ClX…L interaction energy strongly depending on the nature of L, iii) a quenching of one of the π back-donation components to the bond. In the ClTs(CO) adduct, the back-donation from ClTs to CO becomes the most important component. The analysis of the electronic structure of the ClX dimers allows for a clear interpretation of the SOC effects in these systems.

Estimation of thermodynamic and physicochemical properties of the alkali astatides: On the bond strength of molecular astatine (At2 ) and the hydration enthalpy of astatide (At- )

The recent accurate and precise determination of the electron affinity (EA) of the astatine atom At0 warrants a re-investigation of the estimated thermodynamic properties of At0 and astatine containing molecules as this EA was found to be much lower (by 0.4 eV) than previous estimated values. In this contribution we estimate, from available data sources, the following thermodynamic and physicochemical properties of the alkali astatides (MAt, M = Li, Na, K, Rb, Cs): their solid and gaseous heats of formation, lattice and gas-phase binding enthalpies, sublimation energies and melting temperatures. Gas-phase charge-transfer dissociation energies for the alkali astatides (the energy requirement for M+ At- ➔ M0  + At0 ) have been obtained and are compared with those for the other alkali halides. Use of Born-Haber cycles together with the new AE (At0 ) value allows the re-evaluation of ΔHf (At0 )g (=56 ± 5 kJ/mol); it is concluded that (At2 )g is a weakly bonded species (bond strength <50 kJ/mol), significantly weaker bonded than previously estimated (116 kJ/mol) and much weaker bonded than I2 (148 kJ/mol), but in agreement with the finding from theory that spin-orbit coupling considerably reduces the bond strength in At2 . The hydration enthalpy (ΔHaq ) of At- is estimated to be -230 ± 2 kJ/mol (using ΔHaq [H+ ] = -1150.1 kJ/mol), in good agreement with molecular dynamics calculations. Arguments are presented that the largest alkali halide, CsAt, like the smallest, LiF, will be only sparingly soluble in water, following the generalization from hard/soft acid/base principles that "small likes small" and "large likes large."

In vivo stability of 211At-radiopharmaceuticals: on the impact of halogen bond formation

211At, when coupled to a targeting agent, is one of the most promising radionuclides for therapeutic applications. The main labelling approach consists in the formation of astatoaryl compounds, which often show a lack of in vivo stability. The hypothesis that halogen bond (XB) interactions with protein functional groups initiate a deastatination mechanism is investigated through radiochemical experiments and DFT modelling. Several descriptors agree on the known mechanism of iodoaryl substrates dehalogenation by iodothyronine deiodinases, supporting the higher in vivo dehalogenation of N-succinimidyl 3-[211At]astatobenzoate (SAB) conjugates in comparison with their iodinated counterparts. The guanidinium group in 3-[211At]astato-4-guanidinomethylbenzoate (SAGMB) prevents the formation of At-mediated XBs with the selenocysteine active site in iodothyronine deiodinases. The initial step of At-aryl bond dissociation is inhibited, elucidating the better in vivo stability of SAGMB conjugates compared with those of SAB. The impact of astatine's ability to form XB interactions on radiopharmaceutical degradation may not be limited to the case of aryl radiolabeling.

Solid-State 19F NMR Chemical Shift in Square-Planar Nickel-Fluoride Complexes Linked by Halogen Bonds

The halogen bond (XB) is a highly directional class of noncovalent interactions widely explored by experimental and computational studies. However, the NMR signature of the XB has attracted limited attention. The prediction and analysis of the solid-state NMR (SSNMR) chemical shift tensor provide useful strategies to better understand XB interactions. In this work, we employ a computational protocol for modeling and analyzing the ¹⁹F SSNMR chemical shifts previously measured in a family of square-planar trans Ni^II-L₂-iodoaryl-fluoride (L = PEt₃) complexes capable of forming self-complementary networks held by a NiF···I(C) halogen bond [Thangavadivale, V.; Chem. Sci. 2018, 9, 3767-3781]. To understand how the ¹⁹F NMR resonances of the nickel-bonded fluoride are affected by the XB, we investigate the origin of the shielding in trans-[NiF(2,3,5,6-C₆F₄I)(PEt₃)₂], trans-[NiF(2,3,4,5-C₆F₄I)(PEt₃)₂], and trans-[NiF(C₆F₅)(PEt₃)₂] in the solid state, where a XB is present in the two former systems but not in the last. We perform the ¹⁹F NMR chemical shift calculations both in periodic and molecular models. The results show that the crystal packing has little influence on the NMR signatures of the XB, and the NMR can be modeled successfully with a pair of molecules interacting via the XB. Thus, the observed difference in chemical shift between solid-state and solution NMR can be essentially attributed to the XB interaction. The very high shielding of the fluoride and its driving contributor, the most shielded component of the chemical shift tensor, are well reproduced at the 2c-ZORA level. Analysis of the factors controlling the shielding shows how the highest occupied Ni/F orbitals shield the fluoride in the directions perpendicular to the Ni-F bond and specifically perpendicular to the coordination plane. This shielding arises from the magnetic coupling of the Ni(3d)/F(2p lone pair) orbitals with the vacant σ_Ni-F^* orbital, thereby rationalizing the very highly upfield (shielded) resonance of the component (δ₃₃) along this direction. We show that these features are characteristic of square-planar nickel-fluoride complexes. The deshielding of the fluoride in the halogen-bonded systems is attributed to an increase in the energy gap between the occupied and vacant orbitals that are mostly responsible for the paramagnetic terms, notably along the most shielded direction.

Hydrogen Bonding with Polonium

Hydrogen bonding (H-bonding) with heavier chalcogens such as polonium and tellurium is almost unexplored owing to their lower electronegativities, providing us an opportunity to delve into the uncharted territory of...

An expanded halogen bonding scale using astatine

As a non-covalent interaction, halogen bonding is now acknowledged to be useful in all fields where the control of intermolecular recognition plays a pivotal role. Halogen-bond basicity scales allow quantification of the halogen bonding of referential donors with organic functional groups from a thermodynamic point of view. Herein we present the pK BAtI basicity scale to provide the community an overview of halogen-bond acceptor strength towards astatine, the most potent halogen-bond donor element. This experimental scale is erected on the basis of complexation constants measured between astatine monoiodide (AtI) and sixteen selected Lewis bases. It spans over 6 log units and culminates with a value of 5.69 ± 0.32 for N,N,N',N'-tetramethylthiourea. On this scale, the carbon π-bases are the weakest acceptors, the oxygen derivatives cover almost two-thirds of the scale, and sulphur bases exhibit the highest AtI basicity. Regarding the applications of 211At in targeted radionuclide therapy, stronger labelling of carrier agents could be envisaged on the basis of the pK BAtI scale.

Advances in the Chemistry of Astatine and Implications for the Development of Radiopharmaceuticals

ConspectusAstatine (At) is the rarest on Earth of all naturally occurring elements, situated below iodine in the periodic table. While only short-lived isotopes (t1/2 ≤ 8.1 h) are known, 211At is the object of growing attention due to its emission of high-energy alpha particles. Such radiation is highly efficient to eradicate disseminated tumors, provided that the radionuclide is attached to a cancer-targeting molecule. The interest in applications of 211At in nuclear medicine translates into the increasing number of cyclotrons able to produce it. Yet, many challenges related to the minute amounts of available astatine are to be overcome in order to characterize its physical and chemical properties. This point is of paramount importance to develop synthetic strategies and solve the labeling instability in current approaches that limits the use of 211At-labeled radiopharmaceuticals. Despite its discovery in the 1940s, only the past decade has seen a significant rise in the understanding of astatine's basic chemical and radiochemical properties, thanks to the development of new analytical and computational tools.In this Account, we give a concise summary of recent advances in the determination of the physicochemical properties of astatine, putting in perspective the duality of this element which exhibits the characteristics both of a halogen and of a metal. Striking features were evidenced in the recent determination of its Pourbaix diagram such as the identification of stable cationic species, At+ and AtO+, contrasting with other halogens. Like metals, these species were shown to form complexes with anionic ligands and to exhibit a particular affinity for organic species bearing soft donor atoms. On the other hand, astatine shares many characteristics with other halogen elements. For instance, the At- species exists in water, but with the least range of EH-pH stability in the halogen series. Astatine can form molecular interactions through halogen bonding, and it was only recently identified as the strongest halogen-bond donor. This ability is nonetheless affected by relativistic effects, which translate to other peculiarities for this heavy element. For instance, the spin-orbit coupling boosts astatine's propensity to form charge-shift bonds, catching up with the behavior of the lightest halogens (fluorine, chlorine).All these new data have an impact on the development of radiolabeling strategies to turn 211At into radiopharmaceuticals. Inspired by the chemistry of iodine, the chemical approaches have sparsely evolved over the past decades and have long been limited to electrophilic halodemetalation reactions to form astatoaryl compounds. Conversely, recent developments have favored the use of the more stable At- species including the aromatic nucleophilic substitution with diaryliodonium salts or the copper-catalyzed halodeboronation of arylboron precursors. However, it is clear that new bonding modalities are necessary to improve the in vivo stability of 211At-labeled aryl compounds. The tools and data gathered over the past decade will contribute to instigate original strategies for overcoming the challenges offered by this enigmatic element. Alternatives to the C-At bond such as the B-At and the metal-At bonds are typical examples of exciting new axes of research.

Astatine Facing Janus: Halogen Bonding vs. Charge-Shift Bonding

The nature of halogen-bond interactions was scrutinized from the perspective of astatine, potentially the strongest halogen-bond donor atom. In addition to its remarkable electronic properties (e.g., its higher aromaticity compared to benzene), C₆At₆ can be involved as a halogen-bond donor and acceptor. Two-component relativistic calculations and quantum chemical topology analyses were performed on C₆At₆ and its complexes as well as on their iodinated analogues for comparative purposes. The relativistic spin–orbit interaction was used as a tool to disclose the bonding patterns and the mechanisms that contribute to halogen-bond interactions. Despite the stronger polarizability of astatine, halogen bonds formed by C₆At₆ can be comparable or weaker than those of C₆I₆. This unexpected finding comes from the charge-shift bonding character of the C–At bonds. Because charge-shift bonding is connected to the Pauli repulsion between the bonding σ electrons and the σ lone-pair of astatine, it weakens the astatine electrophilicity at its σ-hole (reducing the charge transfer contribution to halogen bonding). These two antinomic characters, charge-shift bonding and halogen bonding, can result in weaker At-mediated interactions than their iodinated counterparts.

Relativistic four-component potential energy curves for the lowest 23 covalent states of molecular astatine (At2)

The potential energy curves (PECs) of all covalent states of Molecular Astatine (At₂) have been investigated in this work within a four-component relativistic framework using the MOLFDIR program package. The ground state was determined using multireference configuration interaction with all single and double excitations including Davidson size-extensivity correction (MRCISD+Q) whereas the 22 excited states were treated by complete open shell configuration interaction (COSCI). Spectroscopic constants (R_e,ω_e,ω_ex_e,ω_ey_e, D_e,B_e,α_e,β_e,T_e) are presented for all states as well as vertical excitations obtained at COSCI, MRCISD and MRCISD+Q levels. In addition, it is also presented accurate extended Rydberg analytical form for the ground state X: (1)0_g⁺.

Delocalized relativistic effects, from the viewpoint of halogen bonding

The ability of organic and inorganic compounds bearing both iodine and astatine atoms to form halogen-bond interactions is theoretically investigated. Upon inclusion of the relativistic spin-orbit interaction, the I-mediated halogen bonds are more affected than the At-mediated ones in many cases. This unusual outcome is disconnected from the behavior of iodine's electrons. The significant decrease of astatine electronegativity with the spin-orbit coupling triggers a redistribution of the electron density, which propagates relativistic effects toward the distant iodine atom. This mechanism can be controlled by introducing suitable substituents and, in particular, strengthened by taking advantage of electron-withdrawing inductive and mesomeric effects. Noticeable relativistic effects can actually be transferred to light atoms properties, e.g., the halogen-bond basicity of bridgehead carbon atoms doubled in propellane derivatives.

The electron affinity of astatine

One of the most important properties influencing the chemical behavior of an element is the electron affinity (EA). Among the remaining elements with unknown EA is astatine, where one of its isotopes, 211At, is remarkably well suited for targeted radionuclide therapy of cancer. With the At- anion being involved in many aspects of current astatine labeling protocols, the knowledge of the electron affinity of this element is of prime importance. Here we report the measured value of the EA of astatine to be 2.41578(7) eV. This result is compared to state-of-the-art relativistic quantum mechanical calculations that incorporate both the Breit and the quantum electrodynamics (QED) corrections and the electron-electron correlation effects on the highest level that can be currently achieved for many-electron systems. The developed technique of laser-photodetachment spectroscopy of radioisotopes opens the path for future EA measurements of other radioelements such as polonium, and eventually super-heavy elements.

Astatine partitioning between nitric acid and conventional solvents: indication of covalency in ketone complexation of AtO. Chem Commun (Camb) 2020;56:9004-9007. [PMID: 32638758 DOI: 10.1039/d0cc03804k]

Astatine-211 has been produced at Texas A&M University on the K150 cyclotron, with a yield of 890 ± 80 MBq through the ²⁰⁹Bi(α,2n)²¹¹At reaction via an 8 h bombardment with a beam current of 4-8 μA and an α-particle beam energy of 28.8 MeV. The target was then dissolved in HNO₃ and the extraction of ²¹¹At was investigated into a variety of organic solvents in 1-3 M HNO₃. Extraction of ²¹¹At with distribution ratios as high as 11.3 ± 0.6, 12.3 ± 0.8, 42.2 ± 2.2, 69 ± 4, and 95 ± 6 were observed for diisopropyl ether, 1-decanol, 1-octanol, 3-octanone, and methyl isobutyl ketone, respectively, while the distribution ratios for ²⁰⁷Bi were ≤0.05 in all cases. The extraction of ²¹¹At into both methyl isobutyl ketone and 3-octanone showed a strong, linear dependence on the HNO₃ initial aqueous concentration and better extraction than other solvents. DFT calculations show stronger binding between the carbonyl oxygen of the ketone and the At metal center.

Quantum chemical topology at the spin-orbit configuration interaction level: Application to astatine compounds

We report a methodology that allows the investigation of the consequences of the spin-orbit coupling by means of the QTAIM and ELF topological analyses performed on top of relativistic and multiconfigurational wave functions. In practice, it relies on the "state-specific" natural orbitals (NOs; expressed in a Cartesian Gaussian-type orbital basis) and their occupation numbers (ONs) for the quantum state of interest, arising from a spin-orbit configuration interaction calculation. The ground states of astatine diatomic molecules (AtX with X = AtF) and trihalide anions (IAtI^- , BrAtBr^- , and IAtBr^- ) are studied, at exact two-component relativistic coupled cluster geometries, revealing unusual topological properties as well as a significant role of the spin-orbit coupling on these. In essence, the presented methodology can also be applied to the ground and/or excited states of any compound, with controlled validity up to including elements with active 5d, 6p, and/or 5f shells, and potential limitations starting with active 6d, 7p, and/or 6f shells bearing strong spin-orbit couplings.

Towards a Stronger Halogen Bond Involving Astatine: Unexpected Adduct with Bu3 PO Stabilized by Hydrogen Bonding

The halogen bond is a powerful tool for the molecular design and pushing the limits of its strength is of major interest. Bearing the most potent halogen-bond donor atom, astatine monoiodide (AtI) was recently successfully probed [Nat. Chem. 2018, 10, 428-434]. In this work, we continue the exploration of adducts between AtI and Lewis bases with the tributylphosphine oxide (Bu₃ PO) ligand, revealing the unexpected experimental occurrence of two distinct chemical species with 1:1 and 2:1 stoichiometries. The 1:1 Bu₃ PO⋅⋅⋅AtI complex is found to exhibit the strongest astatine-mediated halogen bond so far (with a formation constant of 10^(4.24±0.35) ). Quantum chemical calculations unveil the intriguing nature of the 2:1 2Bu₃ PO⋅⋅⋅AtI adduct, involving a halogen bond between AtI and one Bu₃ PO molecular unit plus CH⋅⋅⋅O hydrogen bonds chelating the second Bu₃ PO unit.

Spin-orbit coupling is the key to unraveling intriguing features of the halogen bond involving astatine

The effect of spin-orbit coupling (SOC) on the halogen bond involving astatine has been investigated using state-of-the-art two- and four-component relativistic calculations. Adducts between Cl-X (X = Cl, Br, I and At) and ammonia have been selected to establish a trend on going down the periodic table. The SOC influence has been explored not only on the geometric and energetic features that can be used to characterize the halogen bond strength but also on the three main contributions to it that are the charge transfer, the "σ-hole" (i.e. the localized region with a net positive electrostatic potential at the halogen site) and the "polar flattening" (which is related to the effective shape of the halogen site). A surprisingly large increase of the Cl-At dipole moment, due to the inclusion of SOC, has been worked out using four-component CCSD(T) reference calculations, indicating that this bond is significantly more ionic than one may predict. Due to the SOC effect, which induces a peculiar charge accumulation on the At side in the Cl-At dimer, a weakening of the astatine-mediated halogen bond occurs arising from the (i) reduced amount of charge transfer, (ii) decrease of the polar flattening and (iii) lowering of the short-range Coulomb potential. The analysis of the electronic structure of the Cl-At moiety allows for a rationalization of the SOC effects on all the considered features of the halogen bond, including an unprecedented unsymmetrical charge back-donation from Cl-At to ammonia.

On the Interplay between Charge-Shift Bonding and Halogen Bonding

The nature of halogen-bond interactions has been analysed from the perspective of the astatine element, which is potentially the strongest halogen-bond donor. Relativistic quantum calculations on complexes formed between halide anions and a series of Y₃ C-X (Y=F to X, X=I, At) halogen-bond donors disclosed unexpected trends, e. g., At₃ C-At revealing a weaker donating ability than I₃ C-I despite a stronger polarizability. All the observed peculiarities have their origin in a specific component of C-Y bonds: the charge-shift bonding. Descriptors of the Quantum Chemical Topology show that the halogen-bond strength can be quantitatively anticipated from the magnitude of charge-shift bonding operating in Y₃ C-X. The charge-shift mechanism weakens the ability of the halogen atom X to engage in halogen bonds. This outcome provides rationales for outlier halogen-bond complexes, which are at variance with the consensus that the halogen-bond strength scales with the polarizability of the halogen atom.

A new type of halogen bond involving multivalent astatine: an ab initio study

Theoretical studies on the dimers formed by CO with the halides of multivalent astatine as a Lewis-acid center are carried out to examine the typical characteristics of supervalent halogen bonds. Calculations at the MP2/aug-cc-pVTZ level reveal that the multiple nucleophilic sites of multivalent halide monomers can promote the formation of various types of halogen bonds, among which the most stable ones are At-halogen bond complexes with multivalent astatine as a Lewis acid center, followed by the π-halogen bond dimers, and the weakest ones are the X-halogen bonds. Compared with multivalent Cl-, Br-, and I-centers, At, as the heaviest halogen, exhibits the highest halogen-bond donating ability. We found that the electrostatic term and the dispersion term play an important role in the overall attractive interaction energy, and the smallest attraction term for all complexes is the polarization term (ΔEpol). Moreover, the tri and pentavalent halides analyzed here possess very "flexible" tautomerism in which the transformation occurs during the formation of the dimers. AIM theory and NBO analysis are also employed here.