Precise characterization of archetypal systems of aqueous hydrogen-bonding networks is essential for developing accurate potential functions and universal models of water. The structures of water clusters (H2O)n (n = 2-9) have been verified recently through size-specific infrared spectroscopy with a vacuum ultraviolet free electron laser (VUV-FEL) and quantum chemical studies. For (H2O)10, the pentagonal prism and butterfly motifs were proposed to be important building blocks and were observed in previous experiments. Here we report the size-specific infrared spectra of (H2O)10 via a joint experimental and theoretical study. Well-resolved spectra provide a unique signature for the coexistence of pentagonal prism and butterfly motifs. These (H2O)10 motifs develop from the dominant structures of (H2O)n (n = 8, 9) clusters. This work provides an intriguing prelude to the diverse structure of liquid water and opens avenues for size-dependent measurement of larger systems to understand the stepwise formation mechanism of hydrogen-bonding networks.
We present here a set of scalar-relativistic norm-conserving 4f-in-core pseudopotentials, together with complementary valence-shell Gaussian basis sets, for the lanthanide (Ln) series (Ce-Lu). The Goedecker, Teter, and Hutter (GTH) formalism is adopted with the generalized gradient approximation (GGA) for the exchange-correlation Perdew-Burke-Ernzerhof (PBE) functional. The 4f-in-core pseudopotentials are built through attributing 4f-subconfiguration 4fn (n = 1-14) for Ln (Ln = Ce-Lu) into the atomic core region, making it possible to circumvent the difficulty of the description of the open 4fn valence shell. A wide variety of computational benchmarks and tests have been carried out on lanthanide systems including Ln3+-containing molecular complexes, aqueous solutions, and bulk solids to validate the accuracy, reliability, and efficiency of the optimized 4f-in-core GTH pseudopotentials and basis sets. The 4f-in-core GTH pseudopotentials successfully replicate the main features of lanthanide structural chemistry and reaction energetics, particularly for nonredox reactions. The chemical bonding features and solvation shells, hydrolysis energetics, acidity constants, and solid-state properties of selected lanthanide systems are also discussed in detail by utilizing these new 4f-in-core GTH pseudopotentials. This work bridges the idea of keeping highly localized 4f electrons in the atomic core and efficient pseudopotential formalism of GTH, thus providing a highly efficient approach for studying lanthanide chemistry in multi-scale modeling of constituent-wise and structurally complicated systems, including electronic structures of the condensed phase and first-principles molecular dynamics simulations.
Actinoid tetroxide molecules AnO4 (An = Ac - Cm) are investigated with the ab initio density matrix renormalization group (DMRG) approach. Natural orbital shapes are used to read out the oxidation state (OS) of the f-elements, and the atomic orbital energies and radii are used to explain the trends. The highest OSs reveal a "volcano"-type variation: For An = Ac - Np, the OSs are equal to the number of available valence electrons, that is, AcIII , ThIV , PaV , UVI , and NpVII . Starting with plutonium as the turning point, the highest OSs in the most stable AnO4 isomers then decrease as PuV , AmV , and CmIII , indicating that the 5f-electrons are hard to be fully oxidized off from Pu onward. The variations are related to the actinoid contraction and to the 5f-covalency characteristics. Combined with previous work on OSs, we review their general trends throughout the periodic table, providing fundamental understanding of OS-relevant phenomena.
Versatile graphene-like two-dimensional materials with s-, p- and d-block elements have aroused significant interest because of their extensive applications while there is a lack of such materials with f-block elements. Herein we report a unique one composed of the f-block element moiety of uranyl (UO2 2+) through a global-minimum structure search. Its geometry is found to be similar to that of graphene with a honeycomb-like hexagonal unit composed of six uranyl ligands, where each uranyl is bridged by two superoxido groups and a pair of hydroxyl ligands. All the uranium and bridging oxygen atoms form an extended planar 2D structure, which shows thermodynamic, kinetic and thermal stabilities due to σ/π bonding as well as electrostatic interactions between ligands. Each superoxido ligand has one unpaired (2pπ*)1 electron and is antiferromagnetically coupled through uranyl bridges with 2pπ*-5f δ -2pπ* superexchange interactions, forming a rare type of one-dimensional Heisenberg chain with p-orbital antiferromagnetism, which might become valuable for application in antiferromagnetic spintronics.
Lanthanide/actinide separation is a worldwide challenge for atomic energy and nuclear waste treatment. Separation of americium (Am), a critical actinide element in the nuclear fuel cycle, from lanthanides (Ln) is highly desirable for minimizing the long-term radiotoxicity of nuclear waste, yet it is extremely challenging given the chemical similarity between trivalent Am(III) and Ln(III). Selective oxidation of Am(III) to a higher oxidation state (OS) could facilitate this separation, but so far, it is far from satisfactory for practical application as a result of the unstable nature of Am in a high OS. Herein, we find a novel strategy to generate stable pentavalent Am (Am(V)) through coordination of Am(III) with a diglycolamide ligand and oxidation with Bi(V) species in the presence of an organic solvent. This strategy leads to efficient stabilization of Am(V) and an extraordinarily high separation factor (>104) of Am from Ln through one single contact in solvent extraction, thereby opening a new avenue to study the high-OS chemistry of Am and fulfill the crucial task of Ln/Am separation in the nuclear fuel cycle. The synergistic coordination and oxidation process is found to occur in the organic solvent, and the mechanism has been well elucidated by quantum-theoretical modeling.
The 24 trioxide halide molecules MO3X of the manganese group (M = Mn-Bh; X = F-Ts), which are iso-valence-electronic with the famous MnO4- ion, have been quantum-chemically investigated by quasi-relativistic density-functional and ab initio correlated approaches. Geometric and electronic structures, valence and oxidation numbers, vibrational and electronic spectral properties, energetic stabilities of the monomers in the gas phase, and the decay mode of MnO3F have been investigated. The light Mn-3d species are most strongly electron-correlated, indicating that the concept of a closed-shell Lewis-type single-configurational structure [Mn+7(d0) O-2(p6)3 F-(p6)] reaches its limits. The concept of real-valued spin orbitals φ(r)·α and φ(r)·ß breaks down for the heavy Bh-6d, At-6p and Ts-7p elements because of the dominating spin-orbit coupling. The vigorous decomposition of MnO3F at ambient conditions starts by the autocatalyzed release of n O2 and the formation of MnmO3m-2nFm clusters, triggered by the electron-depleted "oxylic" character of the oxide ligands in MnO3X.
We have developed a new set of norm-conserving pseudopotentials and companion Gaussian basis sets for the actinide (An) series (Ac-Lr) using the Goedecker, Teter, and Hutter (GTH) formalism with the Perdew, Burke, and Ernzerhof (PBE) exchange-correlation functional of generalized gradient approximation. To test the accuracy and reliability of the newly parameterized An-GTH pseudopotentials and basis sets, a variety of benchmarks on actinide-containing molecules were carried out and compared to all-electron and available experimental results. The new pseudopotentials include both medium- ([Xe]4f14) and large-core ([Xe]4f145d10) options that successfully reproduce the structures and energetics, particularly redox processes. The medium-core size set, in particular, reproduces all-electron calculations over multiple oxidation states from 0 to VII, whereas the large-core set is suitable only for the early series elements and low oxidation states. The underlying reason for these transferability issues is discussed in detail. This work fills a critical void in the literature for studying the chemistry of 5f-block elements in the condensed phase.
To resolve the fleeting structures of lanthanide Ln3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln3+ ions (Nd3+, Dy3+, Er3+, Lu3+) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. We performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln-O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L3- and L1-edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln-O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln3+ aqua ions.
A complete set of pseudopotentials and accompanying basis sets for all lanthanide elements are presented based on the relativistic, norm-conserving, separable, dual-space Gaussian-type pseudopotential protocol of Goedecker, Teter, and Hutter (GTH) within the generalized gradient approximation (GGA) and the exchange-correlation functional of Perdew, Burke, and Ernzerhof (PBE). The corresponding basis sets have been molecularly optimized (MOLOPT) using a contracted form with a single set of Gaussian exponents for the s, p, and d states. The f states are uncontracted explicitly with Gaussian exponents. Moreover, the Hubbard U values for each lanthanide element, for DFT+U calculations, are also tabulated, allowing for the proper treatment of the strong on-site Coulomb interactions of localized 4f electrons. The accuracy and reliability of our GTH pseudopotentials and companion basis sets optimized for lanthanides is illustrated by a series of test calculations on lanthanide-centered molecules, and solid-state systems, with the most common oxidation states. We anticipate that these pseudopotentials and basis sets will enable larger-scale density functional theory calculations and ab initio molecular dynamics simulations of lanthanide molecules in either gas or condensed phases, as well as of solid state lanthanide-containing materials, allowing further exploration of the chemical and physical properties of lanthanide systems.
Single transition-metal site catalysts with s-, p-, or d-block atom anchor for nitrogen fixation have been extensively studied, and yet the studies of the f-block atom anchor are rarely reported. Thus, we investigate the feasibility of using a newly synthesized U-Co complex featuring a single CoI site coordinated by tetrakis(phophinoamide) and an UIV anchor for N2-to-NH3 conversion by theoretical modeling. We characterize the evolution of oxidation states of U and Co along the reaction pathways from ab initio density matrix renormalization group (DMRG) calculations, and we find that the variation of the Co â U dative bond is correlated with the changes of oxidation states. Both uranium and cobalt can serve as electron reservoirs to facilitate breaking the N-N bond. Our study demonstrates the viability of metal â metal dative bonds, particularly the df-d one, for the reduction of N2 to NH3, and thus, this opens up a new avenue to the rational design of efficient catalyst for nitrogen fixation.
Heteronuclear transition-metal-main-group element carbonyl anion complexes of AFe(CO)3- (A = Ge, Sn, and Pb) are prepared using a laser vaporization supersonic ion source in the gas phase, which were studied by mass-selected infrared (IR) photodissociation spectroscopy. The geometric and electronic structures of the experimentally observed species are identified by a comparison of the measured and calculated IR spectra. These anion complexes have a 2A1 doublet electronic ground state and feature an A[triple bond, length as m-dash]Fe triply bonded C3v structure with all of the carbonyl ligands bonded at the iron center. Bonding analyses of AFe(CO)3- (A = C, Si, Ge, Sn, Pb, and Fl) indicate that the complexes are triply bonded between the valence np atomic orbitals of bare group-14 atoms and the hybridized 3d and 4p atomic orbitals of iron.
Lanthanide elements typically exhibit a +III oxidation state (OS) in chemical compounds with a few in +IV or even +V OS. Although lanthanides with +II OS have been observed recently in organometallic compounds, +I OS is extremely rare. Using a joint photoelectron spectroscopy and quantum theoretical study, we have found two low OS lanthanides in doped boron clusters, PrB3- and PrB4-. These two clusters are shown to have planar structures, in which the Pr atom is bonded to the aromatic boron clusters via two Pr-B σ bonds. Chemical bonding and electronic structure analyses reveal that the Pr atom is in a very low OS in the two boride clusters: +II in PrB3- and +I in PrB4-. The current finding suggests that there should exist a whole class of boride complexes featuring rather low-valent lanthanides and expands the frontier of lanthanide chemistry.
Recycle of thorium is an essential process in the thorium-uranium closed fuel cycle of molten salt reactor (MSR). Pyrochemical treatment of spent nuclear fuel using chloride molten salts as medium has been considered as a promising method. In this article, we performed molecular dynamics simulations on the ThCl4 LiCl molten salts using a polarizable force field parameterized by us from first-principles calculations. The microscopic structures and macroscopic properties at different compositions were investigated using the developed force field to understand the structure/property relationship in the mixture. The differences between ThCl4 LiCl and ThF4 LiF MSs are compared to understand the behaviors of Th4+ in the fluoride-chloride mixed media. In the molten fluorides, the coordination number of Th4+ is larger, and the resulting more shared anions lead to lower ThF dissociation barrier and shorter lifetime of the Th4+ first solvation shell. Our results also indicate the Pauling's structural rules for crystals can be used to rationalize the local structures in molten salts. © 2018 Wiley Periodicals, Inc.
We report a supramolecular strategy to promote radical-mediated Fenton oxidation by the rational design of a folded host-guest complex based on cucurbit[8]uril (CB[8]). In the supramolecular complex between CB[8] and a derivative of 1,4-diketopyrrolo[3,4-c]pyrrole (DPP), the carbonyl groups of CB[8] and the DPP moiety are brought together through the formation of a folded conformation. In this way, the electrostatic effect of the carbonyl groups of CB[8] is fully applied to highly improve the reactivity of the DPP radical cation, which is the key intermediate of Fenton oxidation. As a result, the Fenton oxidation is extraordinarily accelerated by over 100 times. It is anticipated that this strategy could be applied to other radical reactions and enrich the field of supramolecular radical chemistry in radical polymerization, photocatalysis, and organic radical battery and holds potential in supramolecular catalysis and biocatalysis.
Model systems of the FeMo cofactor of nitrogenase have been explored extensively in catalysis to gain insights into their ability for nitrogen fixation that is of vital importance to the human society. Here we investigate the trigonal pyramidal borane-ligand Fe complex by first-principles calculations, and find that the variation of oxidation state of Fe along the reaction path correlates with that of the reverse-dative Fe â B bonding. The redox-flexibility of the reverse-dative Fe â B bonding helps to provide an electron reservoir that buffers and stabilizes the evolution of Fe oxidation state, which is essential for forming the key intermediates of N2 activation. Our work provides insights for understanding and optimizing homogeneous and surface single-atom catalysts with reverse-dative donating ligands for efficient dinitrogen fixation. The extension of this kind of molecular catalytic active center to heterogeneous catalysts with surface single-clusters is also discussed.
The periodic table provides a fundamental protocol for qualitatively classifying and predicting chemical properties based on periodicity. While the periodic law of chemical elements had already been rationalized within the framework of the nonrelativistic description of chemistry with quantum mechanics, this law was later known to be affected significantly by relativity. We here report a systematic theoretical study on the chemical bonding pattern change in the coinage metal dimers (Cu2, Ag2, Au2, Rg2) due to the relativistic effect on the superheavy elements. Unlike the lighter congeners basically demonstrating ns- ns bonding character and a 0g+ ground state, Rg2 shows unique 6d-6d bonding induced by strong relativity. Because of relativistic spin-orbit (SO) coupling effect in Rg2, two nearly degenerate SO states, 0g+ and 2u, exist as candidate of the ground state. This relativity-induced change of bonding mechanism gives rise to various unique alteration of chemical properties compared with the lighter dimers, including higher intrinsic bond energy, force constant, and nuclear shielding. Our work thus provides a rather simple but clear-cut example, where the chemical bonding picture is significantly changed by relativistic effect, demonstrating the modified periodic law in heavy-element chemistry.
Quantum chemical calculations using ab initio methods at the CCSD(T) level and density functional theory have been carried out for the title molecules. The electronic structures of the molecules were analyzed with a variety of charge and energy decomposition methods. The equilibrium geometries of the M2O2 rhombic clusters exhibit very short distances between the transannular metal atoms M = Be, Mg, Ca. The calculated distances are close to standard values between double and triple bonds, but there are no chemical M-M bonds. The metal atoms M carry large positive partial charges, which are even bigger than in diatomic MO. The valence electrons of M are essentially shifted toward oxygen in M2O2, which makes it possible that there is practically no electronic charge in the region between the metal atoms. The bond dissociation energies for fragmentation of M2O2 into two metal oxides MO are very large. The metal-oxide bonds in the rhombic clusters are shorter and stronger than in diatomic MO. A detailed analysis of the electronic structure suggests that there is no significant direct M-M interaction in the M2O2 rhombic clusters, albeit weak three-center M-O-M bonding.
The concept of oxidation state (OS) is based on the concept of Lewis electron pairs, in which the bonding electrons are assigned to the more electronegative element. This approach is useful for keeping track of the electrons, predicting chemical trends, and guiding syntheses. Experimental and quantum-chemical results reveal a limit near +8 for the highest OS in stable neutral chemical substances under ambient conditions. OS=+9 was observed for the isolated [IrO4 ]+ cation in vacuum. The prediction of OS=+10 for isolated [PtO4 ]2+ cations is confirmed computationally for low temperatures only, but hasn't yet been experimentally verified. For high OS species, oxidation of the ligands, for example, of O-2 with formation of . O-1 and O-O bonds, and partial reduction of the metal center may be favorable, possibly leading to non-Lewis type structures.
The experimentally known highest oxidation state of iron has been determined to be Fe(vi) so far. Here we report a combined matrix-isolation infrared spectroscopic and theoretical study of two interconvertible iron oxide anions: a dioxoiron peroxide complex [(η2-O2)FeO2]- with a C2v-structure and a tetroxide FeO4- with a D2d tetrahedral structure, which are formed by co-condensation of laser-ablated iron atoms and electrons with O2/Ar mixtures at 4 K. Quantum chemistry theoretical studies indicate that the Jahn-Teller distorted tetroxide FeO4- anion is a d1 species with hereto the highest iron formal oxidation state Fe(vii).
