Heats of formation of xenon fluorides and the fluxionality of XeF(6) from high level electronic structure calculations

Noble Gas Anions: An Overview of Strategies and Bonding Motifs

This review article aims to provide an overview of the strategies employed to prepare noble gas anions under different environments and experimental conditions, and of the bonding motifs typically occurring in these species. Observed systems include anions fixed into synthesized salts, detected in the gas phase or in high-pressure devices. The major role of the theoretical calculations is also highlighted, not only in support of the experiments, but also as effective in predicting still unreported species. The chemistry of noble gas anions overall appears as a varied and rich paint, offering fascinating opportunities for both experimentalists and theoreticians.

Crystal Structures of Xenon(VI) Salts: XeF5Ni(AsF6)3, XeF5AF6 (A = Nb, Ta, Ru, Rh, Ir, Pt, Au), and XeF5A2F11 (A = Nb, Ta)

Experiments on the preparation of the new mixed cations XeF₅M(AF₆)₃ (M = Cu, Ni; A = Cr, Nb, Ta, Ru, Rh, Re, Os, Ir, Pt, Au, As), XeF₅M(SbF₆)₃ (M = Sn, Pb), and XeF₅M(BF₄)_x(SbF₆)_3-x (x = 1, 2, 3; M = Co, Mn, Ni, Zn) salts were successful only in the preparation of XeF₅Ni(AsF₆)₃. In other cases, mixtures of different products, mostly XeF₅AF₆ and XeF₅A₂F₁₁ salts, were obtained. The crystal structures of XeF₅Ni(AsF₆)₃, XeF₅TaF₆, XeF₅RhF₆, XeF₅IrF₆, XeF₅Nb₂F₁₁, XeF₅Ta₂F₁₁, and [Ni(XeF₂)₂](IrF₆)₂ were determined for the first time on single crystals at 150 K by X-ray diffraction. The crystal structures of XeF₅NbF₆, XeF₅PtF₆, XeF₅RuF₆, XeF₅AuF₆, and (Xe₂F₁₁)₂(NiF₆) were redetermined by the same method at 150 K. The crystal structure of XeF₅RhF₆ represents a new structural type in the family of XeF₅AF₆ salts, which crystallize in four different structural types. The XeF₅A₂F₁₁ salts (M = Nb, Ta) are not isotypic and both represent a new structure type. They consist of [XeF₅]⁺ cations and dimeric [A₂F₁₁]^- anions. The crystal structure of [Ni(XeF₂)₂](IrF₆)₂ is a first example of a coordination compound in which XeF₂ is coordinated to the Ni²⁺ cation.

Prediction of stable radon fluoride molecules and geometry optimization using first-principles calculations

Noble gases possess extremely low reactivity because their valence shells are closed. However, previous studies have suggested that these gases can form molecules when they combine with other elements with high electron affinity, such as fluorine. Radon is a naturally occurring radioactive noble gas, and the formation of radon-fluorine molecules is of significant interest owing to its potential application in future technologies that address environmental radioactivity. Nevertheless, because all isotopes of radon are radioactive and the longest radon half-life is only 3.82 days, experiments on radon chemistry have been limited. Here, we study the formation of radon molecules using first-principles calculations; additionally, possible compositions of radon fluorides are predicted using a crystal structure prediction approach. Similar to xenon fluorides, di-, tetra-, and hexafluorides are found to be stabilized. Coupled-cluster calculations reveal that RnF₆ stabilizes with O_h point symmetry, unlike XeF₆ with C_3v symmetry. Moreover, we provide the vibrational spectra of our predicted radon fluorides as a reference. The molecular stability of radon di-, tetra-, and hexafluoride obtained through calculations may lead to advances in radon chemistry.

A Computational Study on Closed-Shell Molecular Hexafluorides MF6 (M=S, Se, Te, Po, Xe, Rn, Cr, Mo, W, U) - Molecular Structure, Anharmonic Frequency Calculations, and Prediction of the NdF6 Molecule

Quantum chemical methods were used to study the molecular structure and anharmonic IR spectra of the experimentally known closed-shell molecular hexafluorides MF₆ (M=S, Se, Te, Xe, Mo, W, U). First, the molecular structures and harmonic frequencies were investigated using Density Functional Theory (DFT) with all-electron basis sets and explicitly considering the influence of spin-orbit coupling. Second, anharmonic frequencies and IR intensities were calculated with the CCSD(T) coupled cluster method and compared, where available, with IR spectra recorded by us. These comparisons showed satisfactory results. The anharmonic IR spectra provide means for identifying experimentally too little studied or unknown MF₆ molecules with M=Cr, Po, Rn. To the best of our knowledge, we predict the NdF₆ molecule for the first time and show it to be a true local minimum on the potential energy surface. We used intrinsic bond orbital (IBO) analyses to characterize the bonding situation in comparison with the UF₆ molecule.

Predicting the structure and NMR coupling constant 1J(129Xe-19F) of XeF6 using quantum mechanics methods

The XeF₆ molecule exists as a monomer in the gas phase and as the (XeF₆)₄ tetramer in solution. Herein we used distinct quantum mechanics methods to study the conformational equilibrium for the XeF₆ monomer, which is represented mainly by O_h and C_3v symmetric geometries, and for the (XeF₆)₄ structure found in condensate phases. The NMR ¹J(¹²⁹Xe-¹⁹F) coupling constant is predicted using our own NMR-DKH basis set, designed for NMR properties. The C_3v conformer of XeF₆ was stable only with HF, CCSD, and hybrid DFT functionals with at least 28% exact HF exchange. Increasing the % of HF exchange improves the description of the geometry and the O_h→C_3v equilibrium. The BMK, BHandHLYP and LC-ωPBE functionals produce results in excellent agreement with experiments and high-level calculations for the XeF₆ molecule. When it comes to the ¹J(¹²⁹Xe-¹⁹F) coupling constant, the (XeF₆)₄ structure must be considered. For that compound, BHandHLYP leads to the best structure, and BMK leads to the best coupling constant; therefore, the generalized protocol BMK/NMR-DKH//BHandHLYP/def2-SVP is recommended to study the XeF₆ molecule in the gas phase and solution.

Noble Gas Bonding Interactions Involving Xenon Oxides and Fluorides

Noble gas (or aerogen) bond (NgB) can be outlined as the attractive interaction between an electron-rich atom or group of atoms and any element of Group-18 acting as an electron acceptor. The IUPAC already recommended systematic nomenclature for the interactions of groups 17 and 16 (halogen and chalcogen bonds, respectively). Investigations dealing with noncovalent interactions involving main group elements (acting as Lewis acids) have rapidly grown in recent years. They are becoming acting players in essential fields such as crystal engineering, supramolecular chemistry, and catalysis. For obvious reasons, the works devoted to the study of noncovalent Ng-bonding interactions are significantly less abundant than halogen, chalcogen, pnictogen, and tetrel bonding. Nevertheless, in this short review, relevant theoretical and experimental investigations on noncovalent interactions involving Xenon are emphasized. Several theoretical works have described the physical nature of NgB and their interplay with other noncovalent interactions, which are discussed herein. Moreover, exploring the Cambridge Structural Database (CSD) and Inorganic Crystal Structure Database (ICSD), it is demonstrated that NgB interactions are crucial in governing the X-ray packing of xenon derivatives. Concretely, special attention is given to xenon fluorides and xenon oxides, since they exhibit a strong tendency to establish NgBs.

Noble-Gas Chemistry More than Half a Century after the First Report of the Noble-Gas Compound

Recent development in the synthesis and characterization of noble-gas compounds is reviewed, i.e., noble-gas chemistry reported in the last five years with emphasis on the publications issued after 2017. XeF2 is commercially available and has a wider practical application both in the laboratory use and in the industry. As a ligand it can coordinate to metal centers resulting in [M(XeF2)x]n+ salts. With strong Lewis acids, XeF2 acts as a fluoride ion donor forming [XeF]+ or [Xe2F3]+ salts. Latest examples are [Xe2F3][RuF6]·XeF2, [Xe2F3][RuF6] and [Xe2F3][IrF6]. Adducts NgF2·CrOF4 and NgF2·2CrOF4 (Ng = Xe, Kr) were synthesized and structurally characterized at low temperatures. The geometry of XeF6 was studied in solid argon and neon matrices. Xenon hexafluoride is a well-known fluoride ion donor forming various [XeF5]+ and [Xe2F11]+ salts. A large number of crystal structures of previously known or new [XeF5]+ and [Xe2F11]+ salts were reported, i.e., [Xe2F11][SbF6], [XeF5][SbF6], [XeF5][Sb2F11], [XeF5][BF4], [XeF5][TiF5], [XeF5]5[Ti10F45], [XeF5][Ti3F13], [XeF5]2[MnF6], [XeF5][MnF5], [XeF5]4[Mn8F36], [Xe2F11]2[SnF6], [Xe2F11]2[PbF6], [XeF5]4[Sn5F24], [XeF5][Xe2F11][CrVOF5]·2CrVIOF4, [XeF5]2[CrIVF6]·2CrVIOF4, [Xe2F11]2[CrIVF6], [XeF5]2[CrV2O2F8], [XeF5]2[CrV2O2F8]·2HF, [XeF5]2[CrV2O2F8]·2XeOF4, A[XeF5][SbF6]2 (A = Rb, Cs), Cs[XeF5][BixSb1-xF6]2 (x = ~0.37-0.39), NO2XeF5(SbF6)2, XeF5M(SbF6)3 (M = Ni, Mg, Zn, Co, Cu, Mn and Pd) and (XeF5)3[Hg(HF)]2(SbF6)7. Despite its extreme sensitivity, many new XeO3 adducts were synthesized, i.e., the 15-crown adduct of XeO3, adducts of XeO3 with triphenylphosphine oxide, dimethylsulfoxide and pyridine-N-oxide, and adducts between XeO3 and N-bases (pyridine and 4-dimethylaminopyridine). [Hg(KrF2)8][AsF6]2·2HF is a new example of a compound in which KrF2 serves as a ligand. Numerous new charged species of noble gases were reported (ArCH2+, ArOH+, [ArB3O4]+, [ArB3O5]+, [ArB4O6]+, [ArB5O7]+, [B12(CN)11Ne]-). Molecular ion HeH+ was finally detected in interstellar space. The discoveries of Na2He and ArNi at high pressure were reported. Bonding motifs in noble-gas compounds are briefly commented on in the last paragraph of this review.

Covalent and Non-covalent Noble Gas Bonding Interactions in XeF n Derivatives (n = 2-6): A Combined Theoretical and ICSD Analysis

A noble gas bond (also known in the literature as aerogen bond) can be defined as the attractive interaction between any element of group-18 acting as a Lewis acid and any electron rich atom of group of atoms, thus following the IUPAC recommendation available for similar π,σ-hole interactions involving elements of groups 17 (halogens) and 16 (chalcogens). A significant difference between noble gas bonding (NgB) and halogen (HaB) or chalcogen (ChB) bonding is that whilst the former is scarcely found in the literature, HaB and ChB are very common and their applications in important fields like catalysis, biochemistry or crystal engineering have exponentially grown in the last decade. This article combines theory and experiment to highlight the importance of non-covalent NgBs in the solid state of several xenon fluorides [XeF_n]^m+ were the central oxidation state of Xe varies from +2 to +6 and the number of fluorine atoms varies from n = 2 to 6. The compounds with an odd number of fluorine atoms (n = 3 and 5) are cationic (m = 1). The Inorganic Crystal Structural Database (ICSD) strongly evidences the relevance of NgBs in the solid state structures of xenon derivatives. The ability of Xe compounds to participate in π,σ-hole interactions has been studied using different types of electron donors (Lewis bases and anions) using DFT calculations (PBE1PBE-D3/def2-TZVP) and the molecular electrostatic potential (MEP) surfaces.

Theoretical Study of the Reactions of H Atoms with CH3I and CH2I2

High level ab initio methods have been used to provide reliable kinetic data for the H + CH₃I and H + CH₂I₂ gas-phase reactions. The (H, I)-abstraction and I-substitution reaction pathways were identified. The structures were determined on the potential energy surface at the MP2/aug-cc-pVTZ level of theory. The energetics was then refined using the coupled cluster theory. For the iodinated species, the spin-orbit coupling was calculated using the MRCI approach. The core valence and the scalar relativistic corrections were considered. Thermal rate constants were reported using the canonical transition-state theory (TST) and compared to computed values with the canonical variational transition-state theory (CVT) using the zero curvature tunneling (ZCT) and the small curvature tunneling (SCT) corrections over a wide temperature range (250-2500 K) to show the importance of quantum tunneling effects at low temperatures. They are given by the following expressions for the overall reactions using the CVT/SCT method: k_H+CH3I( T) = 1.07 × 10^-17 × T^2.13 exp(2.68 (kJ mol^-1)/ RT) and k_H+CH2I2( T) = 5.73 × 10^-21 × T^2.97 exp(3.15 (kJ mol^-1)/ RT). The I-abstraction is predicted to be the major pathway for both H + CH₃I and H + CH₂I₂ reactions. The obtained kinetic parameters for the H + CH₃I reaction are in excellent agreement with their experimental counterparts over the temperature range 300-750 K. On the basis of our calculated reaction enthalpies, a new evaluation of the standard enthalpy of formation at 298 K of CH₂I and CHI₂ has been provided. Obtained values are Δ_f H°_298K (CH₂I) = 219.5 kJ mol^-1 and Δ_f H°_298K(CHI₂) = 296.3 kJ mol^-1.

Reactivity of Hydrogen Peroxide with Br and I Atoms

The reaction mechanisms of Br and I atoms with H₂O₂ have been investigated using DFT and high-level ab initio calculations. The H-abstraction and OH-abstraction channels were highlighted. The geometries of the stationary points were optimized at the B3LYP/aug-cc-pVTZ level of theory, and the energetics were recalculated with the coupled cluster theory. Spin-orbit coupling for each halogenated species was also explicitly computed by employing the MRCI level of theory. Thermochemistry for HOBr and HOI has been revised and updated standard enthalpies of formation at 298 K for HOBr and HOI are the following: Δ_fH°_298K(HOBr) = (-66.2 ± 4.6) kJ mol^-1 and Δ_fH°_298K(HOI) = (-66.8 ± 4.7) kJ mol^-1. The rate constants have been estimated using transition state theory (TST), canonical variational transition state theory (CVT), and CVT with small curvature tunneling (CVT/SCT) over a wide temperature range (250-2500 K). For the direct abstraction mechanism, the overall rate constant at 300 K was predicted to be 2.58 × 10^-16 and 7.42 × 10^-25 cm³ molecule^-1s^-1 for the Br + H₂O₂ and I + H₂O₂ reactions, respectively. The modified Arrhenius parameters have been estimated for the overall reactions: k_Br+H2O2(T) = 4.80 × 10^-26 T^4.31 exp(-5.51 (kJ mol^-1)/RT) and k_I+H2O2(T) = 3.41 × 10^-23 T^3.29 exp(-56.32 (kJ mol^-1)/RT).

Matrix-Isolation and Quantum-Chemical Analysis of the C3v Conformer of XeF6, XeOF4, and Their Acetonitrile Adducts

A joint experimental-computational study of the molecular structure and vibrational spectra of the XeF₆ molecule is reported. The vibrational frequencies, intensities, and in particular the isotopic frequency shifts of the vibrational spectra for ¹²⁹XeF₆ and ¹³⁶XeF₆ isotopologues recorded in the neon matrix agree very well with those obtained from relativistic coupled-cluster calculations for XeF₆ in the C_3v structure, thereby strongly supporting the observation of the C_3v conformer of the XeF₆ molecule in the neon matrix. A C_3v transition state connecting the C_3v and O_h local minima is located computationally. The calculated barrier of 220 cm^-1 between the C_3v minima and the transition state corroborates the experimental observation of the C_3v conformer and the absence of the O_h conformer in solid noble gas matrices. For comparison matrix-isolation spectra have also been recorded and analyzed for the ¹²⁹XeOF₄ and the ¹³⁶XeOF₄ isotopologues. The matrix-isolation complexation shifts obtained for the XeF₆·NCCH₃ relative to those of free matrix isolated XeF₆ and CH₃CN are in good agreement with those reported for crystalline XeF₆·NCCH₃.

A theoretical study of stabilities, reactivities and bonding properties of XKrOH (X = F, Cl, Br and I) as potential new krypton compounds using coupled cluster, MP2 and DFT calculations

DFT andAb initiocalculations were employed to disclose the conceivable existence of new noble gas molecules, XKrOH.

Infrared spectra of NgBeS (Ng = Ne, Ar, Kr, Xe) and BeS2 in noble-gas matrices

Laser-ablated beryllium atom has been codeposited at 4 K with hydrogen sulfide in excess noble gas matrices. Four noble-gas compounds NgBeS (Ng = Ne, Ar, Kr, Xe) and the BeS(2) molecule are identified on the basis of the S-34 isotopic substitution, DFT and CCSD(T) theoretical predictions, and a comparison of noble-gas substitution. The agreement between the experimental and calculated vibrational frequencies supports the identification of these molecules. The dissociation energies are calculated at 1.6, 12.6, 10.7, and 13.4 kcal/mol for NeBeS, ArBeS, KrBeS, and XeBeS, respectively, at the CCSD(T) level. The BeS Lewis acid molecule favors strong chemical binding between the Be and Ng atoms.

Freezing in resonance structures for better packing: XeF2 becomes (XeF+)(F-) at large compression

Recent high-pressure experiments conducted on xenon difluoride (XeF(2)) suggested that this compound undergoes several phase transitions up to 100 GPa, becoming metallic above 70 GPa. In this theoretical study, in contrast to experiment, we find that the ambient pressure molecular structure of xenon difluoride, of I4/mmm symmetry, remains the most stable one up to 105 GPa. In our computations, the structures suggested from experiment have either much higher enthalpies than the I4/mmm structure or converge to that structure upon geometry optimization. We discuss these discrepancies between experiment and calculation and point to an alternative interpretation of the measured cell vectors of XeF(2) at high pressure. At pressures exceeding those studied experimentally, above 105 GPa, the I4/mmm structure transforms to one of Pnma symmetry. The Pnma phase contains bent FXeF molecules, with unequal Xe-F distances, and begins to bring other fluorines into the coordination sphere of the Xe. Further compression of this structure up to 200 GPa essentially results in self-dissociation of XeF(2) into an ionic solid (i.e., [XeF](+)F(-)), similar to what is observed for nitrous oxide (N(2)O) at high pressure.

Two- and three-dimensional extended solids and metallization of compressed XeF2

The application of pressure, internal or external, transforms molecular solids into extended solids with more itinerant electrons to soften repulsive interatomic interactions in a tight space. Examples include insulator-to-metal transitions in O(2), Xe and I(2), as well as molecular-to-non-molecular transitions in CO(2) and N(2). Here, we present new discoveries of novel two- and three-dimensional extended non-molecular phases of solid XeF(2) and their metallization. At approximately 50 GPa, the transparent linear insulating XeF(2) transforms into a reddish two-dimensional graphite-like hexagonal layered structure of semiconducting XeF(4). Above 70 GPa, it further transforms into a black three-dimensional fluorite-like structure of the first observed metallic XeF(8) polyhedron. These simultaneously occurring molecular-to-non-molecular and insulator-to-metal transitions of XeF(2) arise from the pressure-induced delocalization of non-bonded lone-pair electrons to sp(3)d(2) hybridization in two-dimensional XeF(4) and to p(3)d(5) in three-dimensional XeF(8) through the chemical bonding of all eight valence electrons in Xe and, thereby, fulfilling the octet rule at high pressures.

Heats of formation of XeF(3)(+), XeF(3)(-), XeF(5)(+), XeF(7)(+), XeF(7)(-), and XeF(8) from high level electronic structure calculations

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for XeF(3)(+), XeF(3)(-), XeF(5)(+), XeF(7)(+), XeF(7)(-), and XeF(8) from coupled cluster theory (CCSD(T)) calculations with effective core potential correlation-consistent basis sets for Xe and including correlation of the nearest core electrons. Additional corrections are included to achieve near chemical accuracy of +/-1 kcal/mol. Vibrational zero point energies were computed at the MP2 level of theory. Unlike the other neutral xenon fluorides, XeF(8) is predicted to be thermodynamically unstable with respect to loss of F(2) with the reaction calculated to be exothermic by 22.3 kcal/mol at 0 K. XeF(7)(+) is also predicted to be thermodynamically unstable with respect to the loss of F(2) by 24.1 kcal/mol at 0 K. For XeF(3)(+), XeF(5)(+), XeF(3)(-), XeF(5)(-), and XeF(7)(-), the reactions for loss of F(2) are endothermic by 14.8, 37.8, 38.2, 59.6, and 31.9 kcal/mol at 0 K, respectively. The F(+) affinities of Xe, XeF(2), XeF(4), and XeF(6) are predicted to be 165.1, 155.3, 172.7, and 132.5 kcal/mol, and the corresponding F(-) affinities are 6.3, 19.9, 59.1, and 75.0 kcal/mol at 0 K, respectively.

Roaming is the dominant mechanism for molecular products in acetaldehyde photodissociation

Reaction pathways that bypass the conventional saddle-point transition state (TS) are of considerable interest and importance. An example of such a pathway, termed "roaming," has been described in the photodissociation of H(2)CO. In a combined experimental and theoretical study, we show that roaming pathways are important in the 308-nm photodissociation of CH(3)CHO to CH(4) + CO. The CH(4) product is found to have extreme vibrational excitation, with the vibrational distribution peaked at approximately 95% of the total available energy. Quasiclassical trajectory calculations on full-dimensional potential energy surfaces reproduce these results and are used to infer that the major route to CH(4) + CO products is via a roaming pathway where a CH(3) fragment abstracts an H from HCO. The conventional saddle-point TS pathway to CH(4) + CO formation plays only a minor role. H-atom roaming is also observed, but this is also a minor pathway. The dominance of the CH(3) roaming mechanism is attributed to the fact that the CH(3) + HCO radical asymptote and the TS saddle-point barrier to CH(4) + CO are nearly isoenergetic. Roaming dynamics are therefore not restricted to small molecules such as H(2)CO, nor are they limited to H atoms being the roaming fragment. The observed dominance of the roaming mechanism over the conventional TS mechanism presents a significant challenge to current reaction rate theory.

Xe129 chemical shift by the perturbational relativistic method: Xenon fluorides

(129)Xe nuclear shielding tensor is calculated at the leading-order, one-electron Breit-Pauli perturbation theory (BPPT) level for the xenon fluorides XeF(+), XeF(2), XeF(3) (+), and XeF(4) that cover the large nuclear magnetic resonance chemical shift range of this nucleus. BPPT is found to improve the shift range and relative shifts as compared to the nonrelativistic (NR) theory. While the full BPPT expansion consists of 16 relativistic terms, 5 of them are responsible for the entire chemical shift and shielding anisotropy. The remaining terms are practically isotropic, corelike contributions that are significant for the absolute shielding constant but cancel for the relative chemical shifts. The five principal terms are due to the spin-orbit-modified wave function allowing the Fermi contact and spin-dipole hyperfine interactions to be coupled to the orbital Zeeman interaction, as well as three distinct scalar relativistic modifications of the NR paramagnetic shielding: wave function change due to mass-velocity and Darwin interactions and the relativistic modification of the orbital hyperfine interaction. A very good agreement with the experimental shifts is obtained for XeF(2) and the particularly challenging XeF(+) species when both the NR and the five main relativistic terms are calculated at electron-correlated ab initio levels of theory. The performance of density-functional theory (DFT) with different pure and hybrid exchange-correlation functionals (with increasing exact exchange admixture) is tested against the ab initio data for each individual contribution. It is shown that DFT has difficulties in the description of paramagnetic shielding, already and especially in the NR contribution, which causes a large discrepancy of DFT results with experiment for xenon fluorides. In contrast, the DFT errors for the relativistic terms cancel out to the extent that a fairly good approximation of the total relativistic shift and anisotropy contributions may be obtained. A combination of high-level ab initio NR calculation with hybrid DFT estimates of the five main BPPT terms is proposed for reasonable estimates of xenon chemical shift in molecules. For the difficult cases such as the present XeF(+) and XeF(3) (+) cations, correlated ab initio calculations are unavoidable throughout. None of the other currently available relativistic methods, either at the fully relativistic or a variationally stable quasirelativistic levels of theory, surpasses the quality of the present approach for Xe shifts in these systems.

Energy-consistent relativistic pseudopotentials and correlation consistent basis sets for the 4d elements Y–Pd

Scalar-relativistic pseudopotentials and corresponding spin-orbit potentials of the energy-consistent variety have been adjusted for the simulation of the [Ar]3d(10) cores of the 4d transition metal elements Y-Pd. These potentials have been determined in a one-step procedure using numerical two-component calculations so as to reproduce atomic valence spectra from four-component all-electron calculations. The latter have been performed at the multi-configuration Dirac-Hartree-Fock level, using the Dirac-Coulomb Hamiltonian and perturbatively including the Breit interaction. The derived pseudopotentials reproduce the all-electron reference data with an average accuracy of 0.03 eV for configurational averages over nonrelativistic orbital configurations and 0.1 eV for individual relativistic states. Basis sets following a correlation consistent prescription have also been developed to accompany the new pseudopotentials. These range in size from cc-pVDZ-PP to cc-pV5Z-PP and also include sets for 4s4p correlation (cc-pwCVDZ-PP through cc-pwCV5Z-PP), as well as those with extra diffuse functions (aug-cc-pVDZ-PP, etc.). In order to accurately assess the impact of the pseudopotential approximation, all-electron basis sets of triple-zeta quality have also been developed using the Douglas-Kroll-Hess Hamiltonian (cc-pVTZ-DK, cc-pwCVTZ-DK, and aug-cc-pVTZ-DK). Benchmark calculations of atomic ionization potentials and 4d(m-2)5s(2)-->4d(m-1)5s(1) electronic excitation energies are reported at the coupled cluster level of theory with extrapolations to the complete basis set limit.