Relativity–Induced Bonding Pattern Change in Coinage Metal Dimers

Feb 20, 2018 - properties compared with the lighter dimers, including higher intrinsic bond energy, force constant, and .... The NR, SR, and SO energy...
0 downloads 6 Views 2MB Size
Article Cite This: Inorg. Chem. XXXX, XXX, XXX−XXX

pubs.acs.org/IC

Relativity−Induced Bonding Pattern Change in Coinage Metal Dimers M2 (M = Cu, Ag, Au, Rg) Wan-Lu Li, Jun-Bo Lu, Zhen-Ling Wang, Han-Shi Hu,* and Jun Li* Department of Chemistry and Laboratory of Organic Optoelectronics & Molecular Engineering of the Ministry of Education, Tsinghua University, Beijing 100084, China S Supporting Information *

ABSTRACT: 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.



atomic numbers,25−29 gold and mercury being prominent examples of such.10,30−32 Superheavy elements have especially significant relativistic effects. These elements can be studied by means of the gas-phase chromatography technique33 due to their strong volatilities, whereby they are deposited on the detector surfaces. Pershina and co-workers have performed a series of theoretical calculations to predict the adsorption behavior of many superheavy elements, so as to compare the character of the heavier element with that of lighter homologue in a chemical group.34−38 It is thus interesting to understand how and how much the relativistic effects can affect the conventional chemical bonding pattern and the periodic law.31 Recently, we studied the bond multiplicity change caused by relativistic effect in the group 6 diatomic molecules M2 (M = Cr, Mo, W, Sg) and provided an interesting example that relativity breaks the periodic trend when the element becomes heavier.39 Similar to the group 6 elements with [(n − 1)d5ns1] nonrelativistic configuration, the group 11 elements Cu, Ag, Au, and Roentgenium (Rg) possess [(n − 1)d10ns1] nonrelativistic configurations. Superheavy elements like Rg are synthetic elements not found in nature due to the short half-life period,25,40,41 and their chemical properties are challenging to predict due to enormous relativistic effects. More developed theoretical methods of relativistic quantum chemistry are required for accurate calculations.42 The basic chemical

INTRODUCTION The first recognized periodic table (PT) was developed in 19th century, which has been widely used in illustrating the periodic law or recurring trends in properties of the lighter elements since then.1 The PT has been continuously explored for elements with larger atomic number (Z), with current experimental PT covering element Og (Z = 118) and theoretical ones going much further (e.g., up to Z ≤ 172 by Pyykkö).2−5 Elements falling into the same group (column) generally present similar chemical properties, and to a certain degree d- or f-block elements sometimes tend to have comparable properties as well among some elements in the same period (row). The periodic law of lighter elements is due to nonrelativistic quantum mechanics, as explained by Bohr in 1920s.6−8 For relativistic effect in chemistry, Dirac stated that it would be of little use in the consideration of molecular structures or ordinary chemical reactions in 1929,9 and chemists by and large agreed with his assertion until the 1970s.10,11 Nevertheless, the relativistic effects turn out to be noticeable in all the elements and are significant for the ground states of molecules roughly beginning from the elements in the sixth row (Cs−Rn), especially for lanthanides, actinides, and the superheavy elements of the period table.12−19 For example, the so-called 6s2 inert pair effect refers to a tendency for this pair of 6s2 electrons to resist oxidation due to the stabilization resulting from the dominating direct relativistic ef fects.20−24 Related to relativistic effects, many unexpected and novel properties are brought into the chemistry of the heavier elements with high © XXXX American Chemical Society

Received: February 20, 2018

A

DOI: 10.1021/acs.inorgchem.8b00438 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 1. Relativistic effects on the atomic valence-shell orbital energy levels and radial distributions (D(r) = r2R(r)2) of the electron density. (left) The NR, SR, and SO energy levels of Cu, Ag, Au, and Rg all in comparable (n − 1)d10ns1) configuration by using the ZORA-PBE/TZ2P method (more computational details are in Supporting Information); (right) the radial distribution of 7s and 6d orbitals of Rg.

Figure 2. MO energy levels of Cu2, Ag2, Au2, and Rg2 at SR-ZORA PBE/TZ2P level with the HOMOs shown with two dots and all corresponding energy levels connected by short dash lines. The 3D MO contours (isovalue = 0.04 au) are shown in the right panel, and the MO compositions are listed in Table S1. The σ-MOs are singly degenerate, while the π- and δ-MOs are doubly degenerate.

rows of transition metals.49 In fact, the order of AO energies, (n − 1)d < ns, holds for neutral and positively charged elements of all groups 3−18 in all periods 1 to 6, while the somewhat unusual order ns < (n − 1)d often taught in introductory textbooks only holds in groups 1 and 2 and for the d elements of period 7. For Cu, Ag, Au, and Rg, as dominating direct/ indirect relativistic ef fects50 raise the energies of (n − 1)d orbitals and lower ns orbitals, the s−d gap becomes smaller for Au with respect to Ag and Cu. This order is even further reversed in superheavy element Rg. Such phenomenon has been noted in recent papers on the heavy d-elements, such as the particular stability of Hs8+O4 and Hs8+S4.51−53 Illustrated in Figure 1 are the nonrelativistic (NR), scalar-relativistic (SR), and spin−orbit (SO) energy levels of Cu, Ag, Au, and Rg and the electron density radial distribution of the 7s and 6d orbitals of Rg. Clearly, the energy level of the Rg 7s valence-shell orbital is relativistically stabilized so strongly that it stays energetically lower than the destabilized Rg 6d orbitals and becomes significantly contracted radially. Thus, energetically destabilized and radially expanded Rg 6d5/2 orbitals are relativistically put in

properties of Rg with predicted relativistic 6d97s2 electron configuration39,43 are deemed to differ from those of its lighter congeners, gold, silver, and copper to a certain extent, despite lacking of simple examples to directly show the differences.44,45 The chemical properties of Rg have attracted more interest than that of the two former elements Mt and Ds because of the strong relativistic contraction of the Rg 7s subshell.46−48 To explore how the relativistic effect might alter the chemical bonding picture, we here studied the coinage dimer systems Cu2, Ag2, Au2, and Rg2 for comparison and found that the chemical bonding of Rg2 is drastically different from that of its lighter congeners. The bonding pattern change of this series of simple metal dimer thus provides a textbook example of relativistic effects in periodicity.



RESULTS AND DISCUSSION 1. Different Molecular Orbital Pictures Arising from Distinctive Valence-Shell Energy Levels. According to the common atomic orbital (AO) shell effect, orbital energies of ns orbitals lie above the (n − 1)d ones in third-, fourth-, and fifth B

DOI: 10.1021/acs.inorgchem.8b00438 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

respectively, and the other MOs are in pairs to show bonding and antibonding cancelation without significantly contributing to the chemical bonding for the whole system. But when it comes to Au2, 6s participation in the σ bonding (HOMO−2 for Au2) is reduced to 80% because of the direct and indirect relativistic effects that cause significant s-d hybridization in gold. When moving down to Rg2, both 2σg (HOMO−4) and 1σu (HOMO−3) MOs dominated by 7s character are occupied, so that the net bonding for Rg−Rg single bond changes to involving HOMO−6. This orbital has bonding interaction between sdz2−sdz2 hybridized orbitals (with 52% 6d and 46% 7s character), with the LUMO acting as its corresponding antibonding orbital. In comparison, the NR and SR MO profiles of Rg2 are plotted in Figure S1. The MO contours of Rg2 based on the NR method are basically similar to these of the Au2 dimer. However, all the σ-type orbitals will be greatly stabilized when the SR effect is taken into account, because the direct relativistic ef fect tremendously lowers the energy of the 7s orbital. As a consequence, the bonding character of Rg−Rg is primarily derived from 6d orbitals rather than 7s orbitals, which is also confirmed in the multiconfigurational CASSCF and DMRG calculations below. It is worth noting that, although SO coupling lowers 7p1/2 energy slightly, 7p participation in bonding seems to be less important in Rg2 dimer. To investigate the possible multiconfigurational features, we performed CASSCF and DMRG calculations. The Löwdin natural orbital occupation number (NOON) of each orbital was determined. With the CASSCF (22e, 12o) calculations, NOONs of the s-σ bonding/antibonding paired orbitals of Cu2, Ag2, and Au2 are 1.94/0.06, 1.98/0.02, and 1.96/0.04, respectively, indicating little occupation of the antibonding σ*(s−s) orbital. While the low-energy natural orbitals are nearly fully occupied, for the d-dominated σ-type bonds of Rg2 the NOON is 1.90/0.11, implying that the ground state has a two-configuration feature. Thus, we performed DMRG-CI calculations with active space extending to (22e, 28o) to see how the NOONs change according to the inclusion of more virtual shells. The natural orbitals and NOONs of Rg2 are shown in Figure 3, while those for Cu2, Ag2, and Au2 are shown in Figures S2−S4. Compared with the CASSCF (22e, 12o) results, the Ag2 and Au2 dimers are not susceptible to the size of active space, where the outer-shell orbitals are practically unoccupied. However, in the Cu2 dimer, the NOONs obtained from the DMRG-CI (22e, 28o) method in the occupied region change a little compared with those from the calculations with the small active space, likely because the 3d orbitals of Cu are so contracted that the inclusion of 4p and 4d orbitals can recover some of the dynamic correlation. Moreover, because of the small HOMO− LUMO gap for Rg2, the single-reference DFT approach does not seem to be reliable enough when high accuracy is sought. From Figure 3 and Figures S2−S4, the effective bond orders (EBOs) are calculated, as listed in Table 2. Especially noteworthy is the EBO of Rg2 calculated to be 1.01, which includes 0.94 from 6d-σ, 0.03 from 7s-σ, 0.02 from 6d-π, and 0.01 from 6d-δ character. Overall, the chemical bonding of Rg2 is primarily a single bond from d-type orbitals, different from other lighter dimers mainly from s-type orbitals. 2. Geometric and Chemical Bonding Analyses. Geometries optimized from CCSD(T) and PBE methods are compared with the experimental data (Table 3). Results from CCSD(T) calculations are within 0.03 Å deviation compared to experimental data, which are reasonably better than those from

the valence region for chemical bonding, while low-energy 7s orbital becomes less valence-active, although 7s is still outside 6d in the valence region as shown by the radial distribution function.29,46,47 Different chemical bonding characters for the Rg2 dimer are therefore expected for the distinctive s2d9 configuration with open-d-shell (rather than s1d10) and inherent electronic structures caused by relativistic effects. The bonding pictures of these M2 dimers are simple and well-known from the qualitative molecular orbital (MO) model. Figure 2 shows the SR MO energy levels of the dimers, with the highest occupied molecular orbitals (HOMOs) illustrated by two dots. Take Cu2 as an example: the 3d orbitals form σ, π, δ bonding MOs (1σg, 1πu, 1δg) and σ*, π*, δ* antibonding MOs (1δu, 1πg, 1σu), while the 4s orbitals form a bonding 2σg and antibonding 2σu MO pair. The energy-level “bands” exert larger expansion from Cu2 to Rg2 due to the stronger orbital interactions between the two coinage metals.54 In contrast with Cu2, Ag2 and Au2 dimers with 1σg mainly formed between (n − 1)dσ orbitals, the 1σg MO of Rg2 has rather low energy because of remarkable mixing with the inert 7s orbital (see Table S1 of the Supporting Information). The HOMO−LUMO (LUMO = lowest unoccupied molecular orbital) gap for Rg2 dimers is rather small, as the LUMO becomes very low in energy due to Rg 7s stabilization, giving rise to multiconfigurational feature of the ground state (vide infra) and complicatedness in assigning the ground state when taking SO coupling relativistic effect into account. We then performed time-dependent density functional theory (TDDFT) calculations corrected with SO coupling effect to estimate the effective vertical singlet excitations of the Au2 and Rg2 dimers to qualitatively compare different electronic structures in the HOMO region. As shown in Table 1, the Table 1. Vertical Singlet Excitation Energies and Oscillator Strengths ( f) of Au2 and Rg2 at the SO-ZORA TD-SAOP/ TZ2P Level peak positionb ( f)

dimers Au2 Rg2

calcd expt55,56 calcd

2.34 (0.022), 2.89 (0.016), 3.28 (0.079), 3.86 (0.002) 2.58 (s), 3.13 (m), 3.30 (m), 3.96 (w)a 2.11 (0.002), 2.56 (0.019), 3.60 (0.004), 4.46 (0.002), 4.65 (0.002), 4.70 (0.050), 5.07 (0.010), 5.20 (0.004), 5.29 (0.001), 5.49 (0.008), 5.53 (0.005), 5.58 (0.002), 5.75 (0.034), 6.09 (0.005), 6.20 (0.102), 6.98 (0.008)

a

The relative intensities from experimental values. s = strong, m = medium, w = weak. bThe energies are in electronvolts, and the oscillator strengths are in atomic units.

theoretical excitation energies and oscillator strengths of Au2 using this approach are comparable with the experiment. We also calculated the excitation energies to predict the electronic spectra of Rg2 for comparison and found that more visiblerange excitations exist due to small HOMO−LUMO gap and the strong SO coupling effect. Moreover, the SO energies of the first excited states of Au2 and Rg2 were calculated to be 1.84 and 0.28 eV at the TD-SAOP/TZ2P level, respectively, suggesting much smaller HOMO−LUMO gap for Rg2 than for the lighter congeners, which are consistent with the fore MO discussions. The 12 frontier valence MOs derived from ns and (n − 1)d orbitals are also shown in the right panel of Figure 2. Especially noteworthy is the bonding pattern change from the lighter dimers to Rg2. In Cu2 and Ag2 molecules, there is a single bond with almost pure (more than 96% s-orbital) s−s σg orbital interactions as shown in Cu2 (HOMO−4) and Ag2 (HOMO), C

DOI: 10.1021/acs.inorgchem.8b00438 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 3. Löwdin natural valence orbitals (NOs) of Rg2 from DMRG-CI (22e, 28o), with Löwdin NOONs (isovalue = 0.04 au). For the degenerate π and δ NOs, only one set is shown. The first and third rows are the bonding NOs, and the second and fourth are the corresponding antibonding NOs.

with the canonical MOs and DMRG natural orbitals, strong orbital hybridization occurs in the case of the Rg2 molecule, with 6d (68%) dominating chemical bonding albeit nonnegligible 7s (32%). These results also give evidence to the change of bonding pattern for the Rg2 dimer. We further performed energy decomposition analysis-natural orbitals for chemical valence (EDA-NOCV) analyses for all dimers, and the results are listed in Table 4 to illustrate the orbital energy contribution to chemical bonding. Consistently, the primary bonding of Cu2, Ag2, and Au2 is derived from s-type σ-bonding, accounting for more than 85% of the total orbital interaction energy, with negligible d-type σ-bonding contribution. Differently, d-character turns out to be the dominant factor in the chemical bonding energy of Rg2, responsible for up to 84% of the orbital interaction. Meanwhile, the s-type σbonding descends dramatically to 11%. Overall, it gives convincing evidence that the single bond between the two Rg is mainly contributed by 6d orbitals, with some 7s hybridization. Interestingly, a little d-type π composition is also involved in the bonding of the Rg2 dimer with a contribution of ∼4%. To elucidate the inherent bonding interaction, the adiabatic force constant (ke) was calculated and shown in Table 5. The

Table 2. Effective Bond Orders of M2 (M = Cu, Ag, Au, Rg) Dimer from DMRG-CI (22e, 28o) Method (n − 1)d-σ (n − 1)d-π (n − 1)d-δ ns-σ total

Cu

Ag

Au

Rg

0.01 0 0 0.99 1.01

0 0 0 0.99 0.99

0 0 0 0.98 0.98

0.94 0.02 0.01 0.03 1.01

PBE. The general trend, albeit nonintuitive, that the M−M distances alter as Cu−Cu < Ag−Ag > Au−Au > Rg−Rg is reproduced, and it is also consistent with the atomic radii and the observed M−C distances of MCN (M = Cu, Ag, Au, Rg) system.57−59 In addition the natural bond orbital (NBO) results based on DFT/PBE and CCSD density matrices listed in Table 3 show the orbital hybridization of atoms in the M−M bond down the periodic table for the coinage metal elements. Results given by these two methods are quite consistent with each other. All the dimers have single σ-bonds while with different AO compositions. Again, NBO analyses reveal that, for the Cu, Ag, and Au dimers, the single bonds are mainly composed of ns-type orbitals, with little (n − 1)d participation. Consistent

Table 3. Calculated Geometric Parameters Compared with Experimental Results and NBO Analyses of Cu2, Ag2, Au2, and Rg2 Diatomic Moleculesa bond lengthb (Å) CCSD(T) Cu2 Ag2 Au2 Rg2

2.221 2.545 2.500 2.480

DFT/PBEd 2.222 2.567 2.523 2.513

(2.252) (2.607) (2.526) (2.561)

NBO analysesc ∑Rp60

expt46

DFT/PBE

CCSD

2.24 2.56 2.48 2.42

2.219 2.534 2.470

Cu [4s3d0.03] Ag [5s4d0.01] Au [6s5d0.06] Rg [6d7s0.59]

Cu [4s3d0.01] Ag [5s4d0.01] Au [6s5d0.05] Rg [6d7s0.47]

a

All calculations are at SR level. bTheoretical bond lengths from both CCSD(T) and DFT/PBE methods compared with the summation of Pyykkö covalent radii (∑Rp) and experimental values (expt). cNBO results analyzed from both DFT/PBE and CCSD density matrix. dThe small frozencore approximation was used for the values without parentheses, that is, [1s2-2p6] for Cu, [1s2-3d10] for Ag, [1s2-4d10] for Au, and [1s2-4f14] for Rg. The numbers in parentheses are obtained by using large core approximation for [1s2-3p6] for Cu, [1s2-4p6] for Ag, [1s2-4f14] for Au, and [1s2-5f14] for Rg. D

DOI: 10.1021/acs.inorgchem.8b00438 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry Table 4. EDA-NOCV Analyses for M2 (M = Cu, Ag, Au, Rg) from SR-PBE/TZ2Pa ΔEσ−s orb ΔEσ−d orb ΔEπ−d orb ΔErest orb ΔEorb a

Cu2

Ag2

Au2

Rg2

−34.4 (89%) −2.0 (5%) 0 −2.2 −38.6

−29.1 (91%) −2.0 (6%) 0 −1.0 −32.1

−46.2 (85%) −3.8 (7%) 0 −4.6 −54.6

−11.5 (11%) −90.3 (84%) −4.3 (4%) −1.4 −107.6

All energies are in kilocalories per mole. The electronic configuration of each fragment is identically taken as (n − 1)d10ns1 for comparison.

Table 5. Scalar-Relativistic Adiabatic Force Constants, Bond Energies and Bond Order Indices of M−M (M = Cu, Ag, Au and Rg) IBDE (eV) Cu2 Ag2 Au2 Rg2 a

bond order a

ke (N/cm)

PBE

CCSD(T)

expt

G-J

1.35 1.07 1.75 2.78

2.12 1.67 2.29 3.86

1.97 1.62 2.15 4.39c

2.0771 1.6672 2.2973

1.02 1.01 1.03 1.04

b

Table 6. Comparison of SO and SR Effects on Nuclear Shieldings (in ppm) of the M Atom in M2 Molecules at the ZORA-PBE/TZ2P Levela

Wiberg

b

NR Δ(SR‑NR)b Δ(SO‑NR)b Δ(SO‑SR)b

1.06 1.02 1.04 1.13

a b

c

Gopinathan-Jug bond orders. Wiberg bond orders. Rg (6d107s1) configuration was applied for the IBDE calculation. The bond dissociation energy (BDE) for Rg2 is 1.291 eV at the CCSD(T)/ QZVP level with Rg (6d97s2) electron configuration.45 SO coupling is not included in these calculations.

Cu

Ag

Au

Rg

2184 −45 77 121

4489 32 472 440

9129 2901 3976 3685

14 488 5102 35 184 30 082

Here no frozen-core approximation was used to ensure the accuracy. NR SO NR SO SR Δ(SR‑NR) = σSR iso − σiso , Δ(SO‑NR) = σiso − σiso , Δ(SO‑SR) = σiso − σiso.

the shielding constants for many heavy-element molecules.74−77 Both SR and SO coupling effects are evaluated herein to investigate the relativity-induced changes of the NMR shielding. Obviously, the chemical shifts resulting from the SR and SO coupling effects of Rg2 are dramatically larger than those of the other lighter congener dimers. More s-type orbital contributions will give rise to significant chemical shift caused by SO in NMR,78,79 compared with the nonrelativistic result, which is consistent with the fact that s-type orbitals will be stabilized by the direct relativistic ef fect. To explore the relativistic effects on the bonding character and covalency, electron localization function (ELF) analyses are further done for the coinage metal dimers. The ELF patterns in the bonding region on the basis of the SR level are depicted in Figure 4, also compared with the NR and SO coupled levels. It displays an apparent change in the Rg2 molecule from the NR to the SR level, where more electron density is situated in the bonding region between the two atoms. On account of the relativistic effects, Au2 also has a modest increase in electron density but with almost no change for Ag2 and Cu2. Moreover, SO coupling somewhat changes the ELF distribution of the covalent bond formation in Rg2, as an electron density reduction showing in the central yellow part and an increase showing in blue part slightly deviate from the central part of Figure 4, which is not found in the other congeners. 3. Determining the Ground State of Rg2. As the LUMO is energetically rather close to the HOMO in Rg2 (Figure 2), it is less straightforward to determine the electronic ground state. The actual difference of bonding character of 1 Σ g + ([1σ g 2 1π u 4 2σ g 2 1σ u 2 1δ g 4 ]1π g 4 1δ u 4 2σ u 0 ) and 3 Π u ([···] 1πg31δu42σu1) states is trivial, as the latter involves an electron transfer from d-π* to d-σ* orbitals, both of which are antibonding. The Mayer bond order index is calculated to be 1.07 and 0.90 for 1Σg+ and 3Πu, respectively, with the small decrease likely due to stronger bonding of σ-type than π-type. At the SR level, the 1Σg+ and 3Πu states are nearly degenerate, while the 3Δg ([···]1πg41δu32σu1) state originating from transferring one electron from HOMO to LUMO has much higher excitation energy (1.279 eV) than 1Σg+ (see Figure 5 and Table S2 for details).

force constants at first decline from Cu2 (1.35 N/cm) to Ag2 (1.07 N/cm) due to the increased Pauli repulsion in Ag(d10)− Ag(d10) between the less-contracted 4d semicore orbitals; the σ(s−s) bond of Cu2 is stronger than that of Ag2 as a result of significantly smaller Cu 3d orbitals due to the quantum primogenic effect.61−63 The force constants then increase markedly to 1.75 N/cm for Au2 and 2.78 N/cm for Rg2 as a result of increased s-d hybridization caused by relativistic effects. The various bond orders also present the same trend as the ke values. Specifically, Intrinsic Bond Dissociation Energies (IBDEs)64,65 are calculated in these coinage metal dimer cases with the same atomic electron configurations. IBDE of Rg2 reaches up to 4.391 eV, which is much higher than for the other molecules, partly because d10s1 is not its atomic ground-state configuration, which lies 3.1 eV above the d9s2 ground state at the CCSD(T) level. The bonding strength indicators mentioned above were calculated based on scalar relativistic calculations. While they seem to suggest that there is stronger bonding in Rg2 in comparison with the other dimers, the SO coupling effect to bonding strength is yet to be determined. The difference of bond energies and force constants were also predicted in M2−Ng (M = Cu−Au and Ng = Kr−Rn) systems,66 which can be explained by strong relativity as well as the lanthanide contraction in the case of Au that tends to destabilize the d and f shells and stabilize the s and p shells.12,18,30,32,67−69 Similar trends were evident in the previous calculations on the diatomic compound RgH, which shows that the strength of the roentgenium−hydrogen bond is doubled due to the relativistic effects.46,70 The bonding picture change from σ(s−s) bonding in Cu2, Ag2, and Au2 to σ(d−d) bonding in Rg2 implies that the nuclear shielding of the external magnetic field should be rather different for these two groups of molecules. We therefore calculated the coinage atom shielding constants for M2 using the ADF software with zeroth order regular approximation (ZORA), and the results are listed in Table 6. The SR- and SOZORA methods have been verified to be effective in simulating E

DOI: 10.1021/acs.inorgchem.8b00438 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 4. ELFs for Cu2, Ag2, Au2, and Rg2 calculated using the NR-, SR-, and SO-ZORA methods. The corresponding scale to each item is shown in the right panel.

the Ω state of 2u after SO splitting from 3Πu is competitive with 0g+ (from 1Σg+) as the ground state. In our CASSCF/CASPT2/ SO calculations, the energy difference between 2u and 0g+ is only ∼0.04 eV, which casts doubt on the true ground state of Rg2 due to methodology inaccuracy. Moreover, upon including the SO splitting, 3Δg will split into 3g, 2g, 1g states in energetic ascending order, where the lowest SO state 3g is 0.695 eV higher than 0g+. Since this calculation uses an approximation method to evaluate the energies of the SO states, more accurate methods are needed to deal with the dynamic and static electron correlation as well as the SO coupling effects to determine the true ground state for Rg2 in the future.



CONCLUSIONS

We have studied a series of coinage metal dimers M2 (M = Cu, Ag, Au, Rg) through quantum-chemical calculations to explore the relativistic effects on chemical bonding. This kind of metal dimer represents the simplest yet cogent case to explain the influence of relativistic effects on the change of bonding pattern. The relativistic effects are so important for the group11 superheavy element Rg that the 6d orbitals (with decreasing 7s contribution) play the prominent role in forming chemical bonding in the Rg2 molecule, which differs from the other lighter dimer species with ns−ns bonding. Such strong relativity-enhanced 6d7s hybridization seems to strengthen the Rg−Rg chemical bonding, as shown by the higher adiabatic force constants, intrinsic bond energies, and bond orders calculated at the scalar relativistic level. The calculated nuclear shieldings become much larger when SR and SO coupling effects are taken into account. The electron distribution in the superheavy Rg2 molecule is influenced by both SR and SO coupling relativistic effects. The small energy gap between the 1 + Σg and 3Πu states may result in a triplet-dominanted 0g+ ground state after SO coupling. The coinage metal dimers provide a simple and clear textbook example, where relativistic effects change the chemical bonding pattern and the periodic trends of elements in the same group of the periodic table.

Figure 5. Energy correlation diagram for Rg2: (a) SR state-averageCASPT2 energy levels (left) and (b) SO states from CASSCF/ CASPT2/SO (right).

The first SR-excited state 3Πu will split into 1u, 2u twofold degenerate states and 0u+, 0u− nondegenerate states upon firstorder SO splitting, whereas 1Σg+ turns out to be a 0g+ state being energetically modified only with the second-order SO splitting. The energy of the 3Πu state is much higher than 1Σg+ for Cu2, Ag2, and Au2; thus, the small SO coupling effect is not able to split the state enough so that the ground state is still singlet dominant beyond doubt. But in case of Rg2, 1Σg+ is only 0.42 eV more stable than 3Πu from SR CASPT2/CASSCF calculations (Table S2). For more accurate results, geometry optimization based on CCSD(T)/ QZVP calculations for these two electronic states was also performed; the complete basis set (CBS) extrapolated energies were then obtained by X−3 extrapolation, which serves as reference energy [CCSD(T)/CBS] for CCSDT(Q) calculations with δT+(Q) parameters at the QZVP level. The CCSDT(Q) results show that the energy difference between the SR 1Σg+ and 3Πu states is only 0.31 eV; as a consequence, F

DOI: 10.1021/acs.inorgchem.8b00438 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry



(15) Kaldor, U.; Wilson, S. Theoretical chemistry and physics of heavy and superheavy elements; Springer Science & Business Media, 2003; Vol. 11. (16) Thayer, J. S., Relativistic Effects and the Chemistry of the Heavier Main Group Elements. In Relativistic Methods for Chemists; Barysz, M., Ishikawa, Y., Eds.; Springer Netherlands: Dordrecht, The Netherlands, 2010; pp 63−97. (17) Schwarz, W. H. E. An Introduction to Relativistic Quantum Chemistry. In Relativistic Methods for Chemists; Barysz, M., Ishikawa, Y., Eds.; Springer: Dordrecht, Netherlands, 2010; pp 1−62. (18) Pyykkö, P. Relativistic Effects in Chemistry: More Common Than You Thought. Annu. Rev. Phys. Chem. 2012, 63, 45−64. (19) Autschbach, J. Perspective: Relativistic effects. J. Chem. Phys. 2012, 136, 150902. (20) Drago, R. S. Thermodynamic Evaluation of the Inert Pair Effect. J. Phys. Chem. 1958, 62, 353−357. (21) Sanderson, R. The inert-pair effect on electronegativity. Inorg. Chem. 1986, 25, 1856−1858. (22) Pyykkö , P. Semiempirical relativistic molecular structure calculations. In Methods in computational chemistry; Springer, 1988; pp 137−226. (23) Schwarz, W. H. E.; Wezenbeek, E. M. v.; Baerends, E. J.; Snijders, J. G. The origin of relativistic effects of atomic orbitals. J. Phys. B: At., Mol. Opt. Phys. 1989, 22, 1515−1530. (24) Baerends, E. J.; Schwarz, W. H. E.; Schwerdtfeger, P.; Snijders, J. G. Relativistic atomic orbital contractions and expansions: magnitudes and explanations. J. Phys. B: At., Mol. Opt. Phys. 1990, 23, 3225−3240. (25) Schädel, M. Chemistry of superheavy elements. Angew. Chem., Int. Ed. 2006, 45, 368−401. (26) Iliaš, M.; Kellö, V.; Urban, M. Relativistic effects in atomic and molecular properties. Acta Phys. Slovaca 2010, 60, 259−391. (27) Hu, H. S.; Qiu, Y. H.; Xiong, X. G.; Schwarz, W. H. E.; Li, J. On the maximum bond multiplicity of carbon: unusual C≣U quadruple bonding in molecular CUO. Chem. Sci. 2012, 3, 2786−2796. (28) Xu, W. H.; Jin, X.; Chen, M. H.; Pyykkö, P.; Zhou, M. F.; Li, J. Rare-earth monocarbonyls MCO: comprehensive infrared observations and a transparent theoretical interpretation for M = Sc; Y; La− Lu. Chem. Sci. 2012, 3, 1548−1554. (29) Schwerdtfeger, P. Relativity and Chemical Bonding. In The Chemical Bond; Wiley-VCH Verlag GmbH & Co. KGaA: Hoboken, NJ, 2014; pp 383−404. (30) Schwerdtfeger, P.; Dolg, M.; Schwarz, W. H. E.; Bowmaker, G. A.; Boyd, P. D. W. Relativistic effects in gold chemistry. I. Diatomic gold compounds. J. Chem. Phys. 1989, 91, 1762−1774. (31) Schwerdtfeger, P. Relativistic effects in gold chemistry. 2. The stability of complex halides of gold(III). J. Am. Chem. Soc. 1989, 111, 7261−7262. (32) Pyykkö, P. Theoretical chemistry of gold. II. Inorg. Chim. Acta 2005, 358, 4113−4130. (33) Gäggeler, H. W.; Türler, A. Chemistry of Superheavy Elements; Kluwer: Dordrecht, Netherlands, 2003. (34) Pershina, V.; Borschevsky, A.; Eliav, E.; Kaldor, U. Adsorption of inert gases including element 118 on noble metal and inert surfaces from ab initio Dirac−Coulomb atomic calculations. J. Chem. Phys. 2008, 129, 144106. (35) Pershina, V.; Borschevsky, A.; Eliav, E.; Kaldor, U. Atomic Properties of Element 113 and Its Adsorption on Inert Surfaces from ab Initio Dirac-Coulomb Calculations. J. Phys. Chem. A 2008, 112, 13712−13716. (36) Pershina, V.; Borschevsky, A.; Eliav, E.; Kaldor, U. Prediction of the adsorption behavior of elements 112 and 114 on inert surfaces from ab initio Dirac-Coulomb atomic calculations. J. Chem. Phys. 2008, 128, 024707. (37) Pershina, V.; Anton, J.; Jacob, T. Theoretical predictions of adsorption behavior of elements 112 and 114 and their homologs Hg and Pb. J. Chem. Phys. 2009, 131, 084713. (38) Pershina, V.; Borschevsky, A.; Anton, J.; Jacob, T. Theoretical predictions of trends in spectroscopic properties of gold containing

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.inorgchem.8b00438. Computational method and its corresponding references; tables of the MO compositions in percentage of M2 molecules; relative energies of spin−orbit states of Rg2 dimer; NR and SR MO correlation diagrams of Rg2; NO contours of Cu2, Ag2, and Au2 from DMRG-CI (22e, 28o) calculations (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. (H.S.H.) *E-mail: [email protected]. (J.L.) ORCID

Jun Li: 0000-0002-8456-3980 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors are grateful to Prof. Dr. W. H. E. Schwarz and Prof. Dr. P. Schwerdtfeger for fruitful discussion. This work was supported by the National Natural Science Foundation of China (Grant Nos. 21590792, 21433005, and 91426302). The calculations were partially performed by using supercomputers at the Computer Network Information Center, Chinese Academy of Sciences and Tsinghua National Laboratory for Information Science and Technology.



REFERENCES

(1) Mendeleev, D. On the relationship of the properties of the elements to their atomic weights. Zeitschrift für Chemie 1869, 12, 405− 406. (2) Wang, S. G.; Schwarz, W. H. E. Icon of chemistry: the periodic system of chemical elements in the new century. Angew. Chem., Int. Ed. 2009, 48, 3404−3415. (3) Schwarz, W. H. E.; Rich, R. L. Theoretical Basis and Correct Explanation of the Periodic System: Review and Update. J. Chem. Educ. 2010, 87, 435−443. (4) Schwarz, W. H. E. The Full Story of the Electron Configurations of the Transition Elements. J. Chem. Educ. 2010, 87, 444−448. (5) Pyykkö, P. A suggested periodic table up to Z⩽172, based on Dirac−Fock calculations on atoms and ions. Phys. Chem. Chem. Phys. 2011, 13, 161−168. (6) Bohr, N. The Structure of the Atom. Nature 1923, 112, 29−44. (7) Bohr, N.; Coster, D. Röntgenspektren und periodisches System der Elemente. Eur. Phys. J. A 1923, 12, 342−374. (8) Schwarz, W. H. E. 100th Anniversary of Bohr’s Model of the Atom. Angew. Chem., Int. Ed. 2013, 52, 12228−12238. (9) Dirac, P. A. Quantum mechanics of many-electron systems. Proc. R. Soc. London, Ser. A 1929, 123, 714−733. (10) Pitzer, K. S. Relativistic effects on chemical properties. Acc. Chem. Res. 1979, 12, 271−276. (11) Pyykkö, P.; Desclaux, J. P. Relativity and the periodic system of elements. Acc. Chem. Res. 1979, 12, 276−281. (12) Pyykkö, P. Relativistic effects in structural chemistry. Chem. Rev. 1988, 88, 563−594. (13) Kaltsoyannis, N. Relativistic effects in inorganic and organometallic chemistry. J. Chem. Soc., Dalton Trans. 1997, 1−12. (14) Kaltsoyannis, N. Recent developments in computational actinide chemistry. Chem. Soc. Rev. 2003, 32, 9−16. G

DOI: 10.1021/acs.inorgchem.8b00438 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

(61) Pyykkö, P. Dirac-Fock One-Centre Calculations Part 8. The 1Σ States of ScH, YH, LaH, AcH, TmH, LuH and LrH. Phys. Scr. 1979, 20, 647. (62) Kaupp, M. The role of radial nodes of atomic orbitals for chemical bonding and the periodic table. J. Comput. Chem. 2007, 28, 320−325. (63) Tang, Y.; Zhao, S.; Long, B.; Liu, J.-C.; Li, J. On the Nature of Support Effects of Metal Dioxides MO2 (M = Ti, Zr, Hf, Ce, Th) in Single-Atom Gold Catalysts: Importance of Quantum Primogenic Effect. J. Phys. Chem. C 2016, 120, 17514−17526. (64) Siegel, S.; Siegel, B. Intrinsic bond energies. J. Chem. Educ. 1963, 40, 143−145. (65) Gibson, J. K.; Hu, H.-S.; Van Stipdonk, M. J.; Berden, G.; Oomens, J.; Li, J. Infrared Multiphoton Dissociation Spectroscopy of a Gas-Phase Complex of Uranyl and 3-Oxa-Glutaramide: An Extreme Red-Shift of the [O = U=O]2+ Asymmetric Stretch. J. Phys. Chem. A 2015, 119, 3366−3374. (66) Jamshidi, Z.; Far, M. F.; Maghari, A. Binding of Noble Metal Clusters with Rare Gas Atoms: Theoretical Investigation. J. Phys. Chem. A 2012, 116, 12510−12517. (67) Pyykkö, P. Relativistic Effects on Periodic Trends. In The Effects of Relativity in Atoms, Molecules, and the Solid State; Springer: Boston, MA, 1991; pp 1−13. (68) Pyykkö, P. Theoretical Chemistry of Gold. Angew. Chem., Int. Ed. 2004, 43, 4412−4456. (69) Pyykkö, P. Theoretical chemistry of gold. III. Chem. Soc. Rev. 2008, 37, 1967−1997. (70) Thierfelder, C.; Schwerdtfeger, P.; Koers, A.; Borschevsky, A.; Fricke, B. Scalar relativistic and spin-orbit effects in closed-shell superheavy-element monohydrides. Phys. Rev. A: At., Mol., Opt. Phys. 2009, 80, 022501. (71) Rohlfing, E. A.; Valentini, J. J. UV laser excited fluorescence spectroscopy of the jet-cooled copper dimer. J. Chem. Phys. 1986, 84, 6560−6566. (72) Beutel, V.; Krämer, H. G.; Bhale, G. L.; Kuhn, M.; Weyers, K.; Demtröder, W. High-resolution isotope selective laser spectroscopy of Ag2 molecules. J. Chem. Phys. 1993, 98, 2699−2708. (73) Morse, M. D. Clusters of transition-metal atoms. Chem. Rev. 1986, 86, 1049−1109. (74) Aquino, F.; Govind, N.; Autschbach, J. Scalar Relativistic Computations of Nuclear Magnetic Shielding and g-Shifts with the Zeroth-Order Regular Approximation and Range-Separated Hybrid Density Functionals. J. Chem. Theory Comput. 2011, 7, 3278−3292. (75) Wolff, S. K.; Ziegler, T.; van Lenthe, E.; Baerends, E. J. Density functional calculations of nuclear magnetic shieldings using the zerothorder regular approximation (ZORA) for relativistic effects: ZORA nuclear magnetic resonance. J. Chem. Phys. 1999, 110, 7689−7698. (76) Krykunov, M.; Ziegler, T.; van Lenthe, E. Implementation of a Hybrid DFT Method for Calculating NMR Shieldings Using SlaterType Orbitals with Spin−Orbital Coupling Included. Applications to 187 Os, 195Pt, and 13C in Heavy-Metal Complexes. J. Phys. Chem. A 2009, 113, 11495−11500. (77) Autschbach, J.; Ziegler, T. Nuclear spin−spin coupling constants from regular approximate relativistic density functional calculations. I. Formalism and scalar relativistic results for heavy metal compounds. J. Chem. Phys. 2000, 113, 936−947. (78) Autschbach, J.; Ziegler, T. Relativistic Computation of NMR Shieldings and Spin-Spin Coupling Constants. In Encyclopedia of Nuclear Magnetic Resonance; John Wiley & Sons, Ltd: Hoboken, NJ, 2007; Vol. 9, pp 306−323. (79) Autschbach, J. Density functional theory applied to calculating optical and spectroscopic properties of metal complexes: NMR and optical activity. Coord. Chem. Rev. 2007, 251, 1796−1821.

dimers of the 6p and 7p elements and their adsorption on gold. J. Chem. Phys. 2010, 133, 104304. (39) Wang, Y. L.; Hu, H. S.; Li, W. L.; Wei, F.; Li, J. Relativistic Effects Break Periodicity in Group 6 Diatomic Molecules. J. Am. Chem. Soc. 2016, 138, 1126−1129. (40) Hofmann, S.; Ninov, V.; Heßberger, F.; Armbruster, P.; Folger, H.; Münzenberg, G.; Schött, H.; Popeko, A.; Yeremin, A.; Andreyev, A.; et al. The new element 111. Z. Phys. A: Hadrons Nucl. 1995, 350, 281−282. (41) Karol, P.; Nakahara, H.; Petley, B.; Vogt, E. On the discovery of the elements 110−112 (IUPAC Technical Report). Pure Appl. Chem. 2001, 73, 959−967. (42) Türler, A.; Pershina, V. Advances in the Production and Chemistry of the Heaviest Elements. Chem. Rev. 2013, 113, 1237− 1312. (43) Türler, A. Gas phase chemistry of superheavy elements. J. Nucl. Radiochem. Sci. 2004, 5, R19−R25. (44) Fricke, B., Superheavy elements a prediction of their chemical and physical properties. In Recent impact of physics on inorganic chemistry; Springer, 1975; pp 89−144. (45) Kullie, O.; Zhang, H.; Kolb, D. Relativistic and non-relativistic local-density functional, benchmark results and investigation on the dimers Cu2, Ag2, Au2, Rg2. Chem. Phys. 2008, 351, 106−110. (46) Liu, W. J.; van Wüllen, C. Spectroscopic constants of gold and eka-gold (element 111) diatomic compounds: The importance of spin−orbit coupling. J. Chem. Phys. 1999, 110, 3730−3735. (47) Pershina, V. Relativistic Electronic Structure Theory, Part 2: Applications. In Theoretical and Computational Chemistry; Schwerdtfeger, P., Ed.; Elsevier B. V.: Netherlands, 2004; Vol. 14. (48) Hancock, R. D.; Bartolotti, L. J.; Kaltsoyannis, N. Density Functional Theory-Based Prediction of Some Aqueous-Phase Chemistry of Superheavy Element 111. Roentgenium(I) Is the ‘Softest’ Metal Ion. Inorg. Chem. 2006, 45, 10780−10785. (49) Kramida, A.; Ralchenko, Y.; Reader, J. NIST-ASD Team, NIST Atomic Spectra Database Ionization Energies; National Institute of Standards and Technology: Gaithersburg, MD, 2014. http://physics. nist.gov/PhysRefData/ASD/ionEnergy.html. (50) Rose, S.; Grant, I.; Pyper, N. The direct and indirect effects in the relativistic modification of atomic valence orbitals. J. Phys. B: At. Mol. Phys. 1978, 11, 1171−1176. (51) Huang, W.; Xu, W.-H.; Schwarz, W. H. E.; Li, J. On the Highest Oxidation States of Metal Elements in MO4 Molecules (M = Fe, Ru, Os, Hs, Sm, and Pu). Inorg. Chem. 2016, 55, 4616−4625. (52) Huang, W.; Jiang, N.; Schwarz, W. H. E.; Yang, P.; Li, J. Diversity of Chemical Bonding and Oxidation States in MS4 Molecules of Group 8 Elements. Chem. - Eur. J. 2017, 23, 10580−10589. (53) Hu, S. X.; Li, W. L.; Lu, J. B.; Bao, J. L.; Yu, H. S.; Truhlar, D. G.; Gibson, J. K.; Marçalo, J.; Zhou, M.; Riedel, S.; Schwarz, W. H. E.; Li, J. On the Upper Limits of Oxidation States in Chemistry. Angew. Chem., Int. Ed. 2018, 57, 3242−3245. (54) Hay, P. J.; Thibeault, J. C.; Hoffmann, R. Orbital interactions in metal dimer complexes. J. Am. Chem. Soc. 1975, 97, 4884−4899. (55) Lecoultre, S.; Rydlo, A.; Félix, C.; Harbich, W. Gold dimer in neon: an absorption and fluorescence study. Eur. Phys. J. D 2009, 52, 187−190. (56) Lecoultre, S.; Rydlo, A.; Félix, C.; Buttet, J.; Gilb, S.; Harbich, W. UV−visible absorption of small gold clusters in neon: Aun (n = 1− 5 and 7−9). J. Chem. Phys. 2011, 134, 074302. (57) Lee, D.-K.; Lim, I. S.; Lee, Y. S.; Jeung, G.-H. Relativistic effects on the ground state properties of group 1 and group 11 cyanides estimated from quantum chemical calculations. Int. J. Mass Spectrom. 2008, 271, 22−29. (58) Zaleski-Ejgierd, P.; Patzschke, M.; Pyykkö, P. Structure and bonding of the MCN molecules, M = Cu, Ag, Au, Rg. J. Chem. Phys. 2008, 128, 224303. (59) Liao, M. S.; Schwarz, W. H. E. Effective radii of the monovalent coin metals. Acta Crystallogr., Sect. B: Struct. Sci. 1994, 50, 9−12. (60) Pyykkö, P.; Atsumi, M. Molecular Single-Bond Covalent Radii for Elements 1−118. Chem. - Eur. J. 2009, 15, 186−197. H

DOI: 10.1021/acs.inorgchem.8b00438 Inorg. Chem. XXXX, XXX, XXX−XXX