Puzzling Lack of Temperature Dependence of the ... - ACS Publications

Jan 23, 2017 - open-shell actinide systems have been carried out successfully with ab initio ... electrostatic environment in the embedded system lead...
1 downloads 0 Views 1MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Letter pubs.acs.org/JPCL

Puzzling Lack of Temperature Dependence of the PuO2 Magnetic Susceptibility Explained According to Ab Initio Wave Function Calculations Frédéric Gendron and Jochen Autschbach* Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States S Supporting Information *

ABSTRACT: The electronic structure and the magnetic properties of solid PuO2 are investigated by wave function theory calculations, using a relativistic complete active space (CAS) approach including spin−orbit coupling. The experimental magnetic susceptibility is well reproduced by calculations for an embedded PuO812− cluster model. The calculations indicate that the surprising lack of temperature dependence of the magnetic susceptibility χ of solid PuO2 can be rationalized based on the properties of a single Pu4+ ion in the cubic ligand field of the surrounding oxygen ions. Below ∼300 K, the only populated state is the nonmagnetic ground state, leading to standard temperature-independent paramagnetism (TIP). Above 300 K, there is an almost perfect cancellation of temperature-dependent contributions to χ that depends delicately on the mixing of ion levels in the electronic states, their relative energies, and the magnetic coupling between them.

T

rationalize the magnetic properties of PuO2. However, an energy gap of 123 meV is inconsistent with TIP at high temperatures. Rather, the magnetic susceptibility should become temperature-dependent at around ∼300 K. We also note that with the magnetic transition moment between A1 and T1 obtained from CF theory, the low-temperature TIP susceptibility estimated with a 123 meV energy gap would be more than double the observed value. Clearly, some of the interpretations of the experimental data are at odds with each other. An explanation of the unexpected magnetic susceptibility of PuO2 by theory has proven challenging. Neither CF9−12 nor density functional theory (DFT) calculations were fully satisfactory, and DFT in different flavors tends to produce magnetic GSs.13−18 DFT19,20 and DFT-based dynamical meanfield theory (DMFT)21−23 have predicted Pu−O covalency, which has been used to explain satellite peaks in the core-level photoemission spectra.23 In other words, donation bonding would take place between the ligands and the metal. A DFT− impurity model combination with semiempirical parameters was able to produce an only weakly T-dependent susceptibility24 due to canceling T-dependent terms, although the susceptibility was sensitive to the model parameters. As an (often superior) alternative to DFT, many studies of open-shell actinide systems have been carried out successfully with ab initio multiconfigurational wave function methods.25 Among these is the Complete Active Space Self-Consistent Field (CASSCF)26 approach and its varieties, which allows one (i) to construct all possible configurations within a chosen

he PuO2 solid has been well characterized experimentally. It crystallizes in the cubic CaF2 structure with eightcoordinate octahedral plutonium sites and four-coordinate oxygen sites. The Pu4+ ion has four valence electrons, with the the ability to give rise to a variety of magnetic and nonmagnetic electronic states. Solid PuO2 affords a nonmagnetic electronic ground state (GS), which was characterized both by magnetic susceptibility measurements1,2 and by 17O NMR spectroscopy.3,4 239Pu NMR measurements also suggested that the magnetic moments of the Pu centers vanish.5 A very interesting observation for PuO2 is a temperature-independent paramagnetism (TIP) over a very large range of temperatures, namely, an essentially constant susceptibility χ of 536 × 10−3 cm3 mol−1 between 1 and 1000 K.2 The surprising TIP up to such high temperatures has since attracted much attention, and its origin remains a subject of debate. TIP at lower temperatures can be explained with the help of crystal field (CF) theory. The ground level of the Pu4+ ion is 5 I4. The 9-fold degenerate 5I4 level is split by the cubic CF into one singlet A1, one doublet E, and two triplets T1 and T2.6 The GS of PuO2 is expected to be the nonmagnetic singlet A1. If there is no mixing between the ground and excited ion levels, the TIP susceptibility can only arise from a magnetic coupling between A1 and T12 because the magnetic moment matrix elements involving A1 and the other states of the 5I4 level are zero. The experimental magnetic susceptibility combined with magnetic transition moments from CF theory lead to the estimate that the T1 excited state (ES) should be 284 meV above the GS.2 (In terms of wavenumbers, 1 meV corresponds to 8.07 cm−1.) Neutron inelastic scattering (NIS) experiments later measured this energy gap at 123 meV.7,8 The NIS experiment was thought to confirm the localized nature of the 5f states and, therefore, the validity of using a CF model to © 2017 American Chemical Society

Received: December 17, 2016 Accepted: January 18, 2017 Published: January 23, 2017 673

DOI: 10.1021/acs.jpclett.6b02968 J. Phys. Chem. Lett. 2017, 8, 673−678

Letter

The Journal of Physical Chemistry Letters

spatially three-fold degenerate spin-quintets (5T1g + 25T2g), one spatially two-fold degenerate spin-quintet (5Eg), and two spatially nondegenerate spin-quintets (5A1g + 5A2g). The SO coupling strongly mixes these states among each other as well as with states derived from excited SR triplet terms of Pu4+. The SO GS of PuO2 is the nondegenerate A1. The GS can also be characterized by considering the SO coupling for the Pu4+ ion first. The SR 5I term of the spherical ion is split by the SO interaction into 5IJ levels, with J going from 4 to 8. For a less than half filled f shell, the GS is the nine-fold degenerate J = 4 level. The 5I4 level is then split by the CF interaction into one singlet A1, one doublet E, and two triplets T1 and T2. The combined CF and SO interaction also mixes 5I4 with contributions from excited J = 4 levels (3F4, 3G4, and 3H4). Relative calculated electronic energies for the first four electronic states of the PuO812− cluster, without and with embedding, are collected in Table 1 for the different active

active orbital space and, optionally, a partitioning of the active space with limitations placed on the excitation levels and (ii) to include spin−orbit (SO) coupling via “state interaction” or from the onset. The CAS approach provides a reliable treatment of multiconfigurational character present in openshell actinide systems,25,27−36 and it can be used with allelectron relativistic Hamiltonians. Furthermore, this computational strategy can be applied to solids as long as a finite cluster model is appropriate. For instance, Ungur et al. were able to describe the electronic structures of several fluorite structures of UO2 by using a UO812− molecular cluster surrounded by four shells of point charges to describe the crystal environment.37 Herein, we offer an interpretation of the magnetism of PuO2 that relies only on the electronic structure of the Pu4+ ion surrounded by the nearest-neighbor oxygens. Relativistic CAS wave function results suggest a pseudo-TIP susceptibility of PuO2 at higher temperatures, that is, a near-perfect cancellation of temperature-dependent contributions to χ. The calculations utilize a PuO812− cluster surrounded by an embedding model for the crystal environment, as shown in the graphical abstract. CAS(n,m) and restricted active space (RAS) calculations with n electrons in m orbitals were carried out at the SCF level and with multistate second-order perturbation theory (PT2) for the treatment of dynamic electron correlation, as detailed in the Supporting Information (SI), using the Molcas program.38 Among the active spaces used in this work are (4,7) = four electrons in the seven 5f orbitals, (4,12) = (4,7) augmented with five 6d orbitals, and (12,16) = (4,12) augmented with four ligand-centered occupied orbitals. The ab initio results are conveniently analyzed by referring to the levels of the Pu4+ ion from which the various 5f4 states of PuO2 are derived. For illustration, an energy level diagram is shown in Figure 1. According to Hund’s rules, the spin-free (scalar relativistic, SR) ground term of the Pu4+ ion is 5I, that is, L = 6 and S = 2. The cubic CF in the solid splits 5I into three

Table 1. Relative Calculated Energies (meV) of the LowEnergy Electronic States of an Isolated and CrystalEmbedded PuO812− Cluster of PuO2. CAS/RAS-SCF-SO Results isolated PuO812−

embedded PuO812−

CAS(4,7)

CAS(4,7)

CAS(4,12)

RAS(12,16)

0 107 231 222

0 128 247 323

0 129 251 322

0 134 262 336

A1 T1 T2 E

spaces. The state ordering is A1 (GS), T1, T2, and E. The crystal electrostatic environment in the embedded system leads to a destabilization of the ESs, most notably for E. There are no experimental data for the ordering of the ESs. However, our embedded results are in good agreement with the state ordering calculated using LDA+DMFT, where T1, T2, and E were separated from the GS by 125, 226, and 319 meV, respectively.23 Assignments of the low-energy electronic states of PuO2 at the SCF-SO level are provided in Table 2, based on the Table 2. Assignment of the Low-Energy Electronic States of PuO2 with 5I4 Parentagea state

ΔE (meV)

A1 T1 T2 E

0 134 262 336

weight of 2S+1 LJ 85 86 86 83

5

I4 5 I4 5 I4 5 I4

weight of SR states 49 5T2g(1), 33 5Eg, 3 5T2g(2) 44 5T2g (1), 26 5Eg, 4 5T1g, 8 5T2g(2), 2 5A2g 41 5T2g(1), 27 5T1g, 8 5A1g, 9 5T2g(2) 7 5T2g(1), 24 5Eg, 15 5T1g, 6 5A1g, 23 5 T2g(2), 8 5A2g

a

RAS(12,16)SCF-SO calculations for the embedded cluster model. Weights are in %. The notations used for the contributions refer to the labels used in Figure 1.

embedded cluster model. As expected, on the basis of the qualitative discussion of level splitting and mixing, the 5I4 level contributions are dominant, with about 83−86% weight in the wave functions, with the remaining contributions from excited J = 4 levels. Considering the SO coupling a posteriori, the nondegenerate GS results from an admixture between the lowest SR states 5T2g (1) (49%) and 5Eg (33%) and numerous other small contributions.

Figure 1. Relative energies, in terms of wavenumbers (cm−1), of the low-energy SR and SO states, with the CF treated before (left) and after (right) the SO coupling. CAS(4,7)SCF calculations. E = 0 corresponds to the SR GS. The energies of SR+CF and SR+SO+CF are drawn to scale according to the ab initio calculations for the embedded PuO812− cluster. 674

DOI: 10.1021/acs.jpclett.6b02968 J. Phys. Chem. Lett. 2017, 8, 673−678

Letter

The Journal of Physical Chemistry Letters

excludes the 6d shell but includes the four ligand-centered orbitals produces nearly identical state energies and magnetic data (SI). A relatively weak antibonding interaction between the 5f Pu and 2p O orbitals is visible for the Pu 5fa2u NO, with an occupation of 0.158. There is a ligand-centered a2u NO (occupation = 1.998) featuring the corresponding bonding interaction, but the contribution from the metal orbital is seen to be small in Figure 2. The fact that SO coupling shifts electron density at the metal center from the formally occupied 5ft2u and 5ft1u also into the 5fa2u orbital, which exhibits some Pu−O antibonding character, means that the covalency involving the 5f shell is quite weak. The calculated 5f and 6d shell occupations and the weak metal−oxygen hybridization in the SO GS suggest that the GS electronic structure of PuO2 can be effectively rationalized by using a single-ion picture. Recent LDA+DMFT calculations found strong 5f−2p hybridizations, leading to 5f shell occupation numbers much exceeding four.23,24 In this context, it is important to mention that KS DFT-based calculations probably overestimate the orbital mixing due to the KS delocalization error.40,41 Our SCF calculations do not include dynamic electron correlation, and therefore, they may inherit an opposite “localization error” from Hartree−Fock theory and underestimate any covalent effects. Consequently, the effect of dynamic correlation was investigated at the PT2 level for CAS(4,7) and CAS(4,12). The resulting relative calculated electronic energies are collected in Table S2 in the SI, and the corresponding NOs are shown in Figures S3 and S5. The introduction of the dynamic correlation leads to a slight energetic destabilization of the ESs when compared to the corresponding SCF-SO results. However, there is also a noticeable symmetry breaking of the PT2 wave functions, which leads to an undesired energetic splitting of the degenerate states. For instance, the three components of the first excited triplet T1 are split by more than 13 meV with CAS(4,12)PT2-SO. Furthermore, the comparison between the NOs obtained from the SCF-SO and PT2-SO calculations reveals that the dynamic correlation only leads to a minor increase of the 5f−2p mixing. For example, with the CAS(4,12) active space, some minor additional 2p oxygen character shows up in the 5ft1u and 6deg orbitals. Overall, the introduction of the correlation effects via PT2 does not increase the metal−ligand orbital mixing to a degree that would warrant a reclassification of the system as covalent. Because of the spurious symmetry breaking with PT2, only the magnetic properties calculated at the SCF level are discussed in the following section. We note in this context that CASSCF is usually considered to be suitable for magnetic susceptibility calculations.37,42,43 The calculated magnetic moment matrix elements μij for the states of Tables 1 and 2 are given in Table 3. Due to the cubic

Of particular interest for this study is the energy gap between A1 and T1. This gap is found to be only weakly dependent on the size of the active space. When the active space includes 5f, 6d, and four formally doubly occupied ligand-centered orbitals of symmetry a2u and t2g, that is, RAS(12,16), the excited triplet T1 is calculated to be 134 meV above the GS. This result is in good agreement with the INS experiment (ΔE = 123 meV)8 and at the upper end of the range determined in previous computational studies (97−134 meV).10−12,23,24,39 Natural orbitals (NOs) and occupation numbers for the SO GS of PuO2 calculated with the RAS(12,16) space are shown in Figure 2 (see Figures S2 and S4 of the SI for the CAS(4,7) and

Figure 2. Selected NOs and occupations for the SO GS of PuO2. RAS(12,16)SCF-SO calculations for the embedded cluster model. Isosurfaces at ±0.03 au.

CAS(4,12) active spaces). The occupations of the 5f orbitals do not depend strongly on the size of the active space (see SI). The combined 5f shell occupation of the SO GS obtained from the relativistic ab initio calculations is nf = 3.99, that is, very close to the idealized Pu4+ 5f occupation. The four unpaired electrons are predominantly distributed over the 5ft2u and 5ft1u orbitals, with combined occupations of 3 × 0.726 = 2.178 and 3 × 0.552 = 1.656, respectively. This occupation pattern arises from mixing of the SR states 5T2g(1) and 5Eg (Table 2). The SR 5 T2g(1) state has a dominant configuration (5ft2u)3(5ft1u)1 (72%), while 5Eg has a dominant configuration (5ft2u)2(5ft1u)2 (80%). When SO coupling mixes these configurations, the combined 5ft2u occupation ends up between 2 and 3 (closer to 2 in our case), while the 5ft1u occupation is between 1 and 2. At least some of the 5fa2u occupation is also due to this particular state mixing by the SO coupling; the 5Eg state contains a configuration (5ft1u)2(5ft2u)1(5fa2u)1 with 17% weight, and there are contributions to the 5fa2u occupation from SR spin-triplets. The 6d shell of the Pu center remains basically unoccupied; the occupation numbers of the corresponding metal-centered NOs do not exceed 0.002, and the metal character of the doubly occupied ligand-centered orbitals of the same symmetry is very small. Indeed, calculations with a (12,11) active space that

Table 3. Calculated Magnetic Moment Matrix Elements μij between PuO2 SO States of 5I4 Parentagea μij

A1

T1

T2

E

A1 T1 T2 E

0.000 0.623 0.000 0.000

0.623 0.928 0.270 0.694

0.000 0.270 1.124 0.662

0.000 0.694 0.662 0.000

a RAS(12,16)SCF-SO results are in units of μB. For degenerate states, the listed value corresponds to the root-mean-square. The use of these tabulated μij values in eq 1 leads to curves (a)−(c) of Figure 3.

675

DOI: 10.1021/acs.jpclett.6b02968 J. Phys. Chem. Lett. 2017, 8, 673−678

Letter

The Journal of Physical Chemistry Letters site symmetry of Pu, the magnetic moment vector components are equivalent, and therefore, only the calculated data for one magnetic axis are discussed in the following section. The GS of PuO2 does not exhibit a magnetic moment as it is not degenerate. The excited non-Kramers doublet E is also calculated to be nonmagnetic (see the SI for further details on the E state). The triplet states T1 and T2 exhibit nonzero magnetic moments with calculated root-mean-square ⟨μ⟩ values of 0.928 and 1.124 μB, respectively. We note in passing that SO coupling heavily influences the magnetic moment matrix elements of PuO2, as expected. With the SO GS being nonmagnetic, off-diagonal magnetic moment matrix elements between the GS and the low-energy ESs are predominantly responsible for the TIP of PuO2 for temperatures up to about 300 K. From Table 3, it is seen that for the states with 5I4 parentage only the magnetic coupling between A1 and T1 is nonzero, as already mentioned. We also point out that large magnetic moment matrix elements are calculated between T1 and E and between T2 and E, along with a small magnetic coupling between T1 and T2. The magnetic moment matrix elements μij of Table 3, along with the state energies listed in Table 1, can be used to calculate the magnetic susceptibility χ of PuO2 from the van Vleck equation χ=

Figure 3. Magnetic susceptibility χ (×10−3 cm3 mol−1) as a function of T (K) calculated from eq 1 (a) using only the states A1 and T1 in the equation, (b) using only the states A1, T1, and T2, (c) using all states with 5I4 parentage, and (d) using all calculated electronic states (Molcas result for the susceptibility). The experimental (Expt.) data points were extracted from graphical material of ref 2. The experimental susceptibility of χ = 0.536 × 10−3 cm−3 mol−1 is represented by the dashed line.

would be consistent with pure 5I4 CF states if the A1/T1 energy gap were 284 meV. However, the ab initio calculations agree well both with the refined experimental A1/T1 energy gap of 123 meV and with the experimental TIP susceptibility. The reason for the excellent agreement in our calculation is that the ab initio calculated magnetic transition moment of 0.623 μB (Table 3) entering the susceptibility is much smaller than the CF estimate of 0.861 determined from the 5I4 CF states (further details are provided in the SI). The reduction by a factor of 0.72 is consistent with the reduced weights of the 5I4 configurations in the wave functions, around 85% (Table 2), because 0.852 = 0.72. The susceptibility based on all of the calculated states is also influenced by magnetic coupling of the GS with higher ESs (vide infra). However, these couplings do not involve the J = 4 levels that mix into the GS, and therefore, the weights of the 5I4 states in A1 and T1 best rationalize the reduction of the magnetic transition moment. Above 300 K, χ(a) becomes temperature-dependent because of the thermal population of the T1 state and decreases to 0.362 × 10−3 cm3 mol−1 at T = 1000 K. Such a T dependence of the susceptibility was not observed experimentally, but it is consistent with a previous CF analysis performed by Coliareti-Tosti et al.10 When contributions from the second triplet T2 are included in eq 1, we obtain curve (b) of Figure 3. Because the magnetic moment matrix elements between A1 and T2 vanish (see Table 3), the calculated TIP in χ(b) at low temperature is the same as in χ(a). Magnetic coupling between the T1 and the T2 states causes a less pronounced decrease of χ(b) at higher T when compared to χ(a). The influence of the introduction of the excited doublet E in eq 1 is shown by curve (c) in Figure 3. In this case, the Van Vleck equation includes all of the states of Pu4+ 5I4 parentage. Because the magnetic coupling between A1 and E vanishes, the TIP segment of χ(c) is not affected by the E state and remains in excellent agreement with the experiment. Just as importantly, however, χ(c) remains practically constant above 300 K and reaches 0.569 × 10−3 cm3 mol−1 at T = 1000 K. The reasons for this behavior are the large magnetic couplings between E and the triplets T1 and T2 (Table 3). The T-dependent Curie contribution from the thermal population of T1 above 300 K is

⎡ 1 μ0 μB 2 ∑ e−βEλ⎢β ∑ |⟨ψλa|L̂u + geSû |ψλa ′⟩|2 ⎢⎣ a , a ′ Q0 λ +2

∑∑ λ ′≠ λ a , a ′

|⟨ψλa|L̂u + geSû |ψλ ′ a ′⟩|2 ⎤ ⎥ ⎥⎦ Eλ ′ − Eλ

(1)

Here, χ is the susceptibility for an applied magnetic field in direction u coinciding with a principal magnetic axis. For PuO2, the susceptibility is isotropic. The summations λ, λ′ go over the set of available electronic states. In principle, this should be the complete set of states of the system, but in practice, the sum is restricted to the states that are actually calculated or a subset thereof. Further, Q0 = ∑λ,a e−βEλ is the partition function and β = 1/(kT). The indices a,a′ count the components of degenerate states. The factors μ0 and μB are the vacuum permeability and Bohr magneton, respectively. The first “Curie” term inside of the brackets on the right-hand side of eq 1 is proportional to β and therefore intrinsically temperature-dependent. This term results from contributions of populated degenerate magnetic states. The second “linear response” (LR) term is not intrinsically temperature-dependent and involves magnetic coupling between different electronic states that describe how state λ responds to the presence of the magnetic field. For a system with a nonmagnetic GS, like PuO2, at sufficiently low temperatures where the population of the ESs is negligible, only the coupling between the GS and ESs in the second term contribute to the susceptibility and TIP is observed. At higher temperatures, Curie terms from populated magnetic ESs and LR contributions from magnetic or nonmagnetic ESs may contribute to the susceptibility. Figure 3 shows several calculated χ(T) that were obtained from eq 1 and the data in Tables 1 and 3. Curve (a) of Figure 3 was calculated by using only the GS A1 and the first ES T1. The magnetic coupling between A1 and T1 gives a TIP χ(a) = 0.562 × 10−3 cm3 mol−1, between 0 K and room temperature. This calculated value of χ(a) is in excellent agreement with the experiment. We remind the reader that the experimental χ 676

DOI: 10.1021/acs.jpclett.6b02968 J. Phys. Chem. Lett. 2017, 8, 673−678

Letter

The Journal of Physical Chemistry Letters then counterbalanced in χ(c) by a contribution arising from the magnetic coupling between T1, T2, and E. The net result is a pseudo-TIP susceptibility up to very high temperature. The ab initio calculated χ(T) from all calculated states with the RAS(12,16) active space is also shown in Figure 3 (curve (d)). This calculated χ is in aceptable agreement with the experimental data. At low temperature, the calculated χ overestimates the experimental value by 0.17 × 10−3 cm3 mol−1. Such a deviation from experiment is not untypical for ab initio susceptibility calculations for actinide systems.36,42,44 The difference between curves (c) and (d) reveals that the calculated TIP of PuO2 is also influenced by magnetic coupling between the GS and ESs derived from other Pu ion levels at higher energies. For instance, nonzero transition magnetic moments were calculated between the SO GS and two triplet states arising from the excited 5I5 ion level. Further, above ∼600 K, a minor temperature dependence is caused by magnetic coupling between the populated T1 state and states derived from excited ion levels. The fact that the TIP susceptibility is overestimated, compared to the experiment, when these states are included suggests (i) that the magnetic coupling between these states and the low-energy electronic states is overestimated, (ii) that the energies of these states are predicted to be too low, or (iii) a combination thereof that may be caused by the various approximations in the calculations. In this case, curve (c) of Figure 3 obtained by including only the 5I4 parentage states in the susceptibility calculation would be the most realistic. Curve (c) is also the one that agrees best with the experiment over the 1000 K temperature range. To summarize, the calculated energy gap between the GS and the first ES and the magnetic susceptibility are in good agreement with experimental data. The puzzling experimental finding that the susceptibility is essentially constant up to 1000 K has a single-ion interpretation based on the ab initio calculations as follows: Up to ∼300 K, magnetic coupling between the nonmagnetic GS A1 and the magnetic first ES T1 is responsible for a standard TIP. Compared to a simple CF model, the A1/T1 magnetic transition moment is significantly reduced due to mixing of states derived from higher-energy ion levels into the low-energy states of 5I4 parentage. Above 300 K, a pseudo-TIP is caused by an almost perfect cancellation of temperature-dependent contributions to the susceptibility that arise from thermal population of the low-energy ESs, as previous results from semiempirical DFT-impurity model calculations suggested.24 While interactions between metal centers in the solid, for example, via electron correlations, cannot be ruled out as contributing to the observed magnetic behavior, the calculations indicate relatively weak covalency between Pu and the oxygen ligands that would help facilitate such interactions.



ORCID

Frédéric Gendron: 0000-0002-1896-3978 Jochen Autschbach: 0000-0001-9392-877X Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors acknowledge support from the U.S. Department of Energy, Office of Basic Energy Sciences, Heavy Element Chemistry program, under Grant DE-SC0001136 (formerly DE-FG02-09ER16066). We thank the Center for Computational Research (CCR) at the University at Buffalo for providing computational resources and Prof. Paul Bagus for stimulating discussions on the topic of the PuO2 electronic structure.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.6b02968. Additional information and data regarding the embedding calculations and plots of orbitals calculated at the SCF-SO and PT2-SO levels (PDF)



REFERENCES

(1) Candela, G. A.; Hutchison, C. A.; Burton, W. B. Magnetic susceptibility of Uranium (IV) and Plutonium (IV) Ion in Cubic Fields. J. Chem. Phys. 1959, 30, 246. (2) Raphael, G.; Lallement, R. Susceptibilité Magnétique de PuO2. Solid State Commun. 1968, 6, 383−385. (3) Tokunaga, Y.; Sakai, H.; Fujimoto, T.; Kambe, S.; Walstedt, R. E.; Ikushima, K.; Yasuoka, H.; Aoki, D.; Homma, Y.; Haga, Y.; Matsuda, T. D.; Ikeda, S.; Yamamoto, E.; Nakamura, A.; Shiokawa, Y.; Nakajima, K.; Arai, Y.; Onuki, Y. NMR Studies of Actinide Dioxides. J. Alloys Compd. 2007, 444-445, 241−245. (4) Martel, L.; Magnani, N.; Vigier, J.-F.; Boshoven, J.; Selfslag, C.; Farnan, I.; Griveau, J.-C.; Somers, J.; Fanghänel, T. High-Resolution Solid-State Oxygen-17 NMR of Actinide-Bearing Compounds: An Insight into the 5f Chemistry. Inorg. Chem. 2014, 53, 6928−6933. (5) Yasuoka, H.; Koutroulakis, G.; Chudo, H.; Richmond, S.; Veirs, D. K.; Smith, A. I.; Bauer, E. D.; Thompson, J. D.; Jarvinen, G. D.; Clark, D. L. Observation of 239Pu Nuclear Magnetic Resonance. Science 2012, 336, 901−904. (6) Lea, K. R.; Leask, M. J. M.; Wolf, W. P. The Raising of Angular Momentum Degeneracy of f-electron Terms by Cubic Crystal Fields. J. Phys. Chem. Solids 1962, 23, 1381−1405. (7) Kern, S.; Loong, C.-K.; Goodman, G. L.; Cort, B.; Lander, G. H. Crystal-Field Spectroscopy of PuO2: Further Complications in Actinides Dioxides. J. Phys.: Condens. Matter 1990, 2, 1933−1940. (8) Kern, S.; Robinson, R. A.; Nakotte, H.; Lander, G. H.; Cort, B.; Watson, P.; Vigil, F. A. Crystal-field transition in PuO2. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 104−106. (9) Santini, P.; Carretta, S.; Amoretti, G.; et al. Multipolar Interactions in f-electron Systems: The Paradigm of Actinide Dioxides. Rev. Mod. Phys. 2009, 81, 807. (10) Colarieti-Tosti, M.; Eriksson, O.; Nordstrom, L.; Wills, J.; Brooks, M. S. S. Crystal-field levels and magnetic susceptibility in PuO2. Phys. Rev. B: Condens. Matter Mater. Phys. 2002, 65, 195102. (11) Magnani, N.; Santini, P.; Amoretti, G.; Caciuffo, R. Perturbative Approach to J mixing in f-Electron Systems: Application to Actinide Dioxides. Phys. Rev. B: Condens. Matter Mater. Phys. 2005, 71, 054405. (12) Magnani, N.; Amoretti, G.; Carretta, S.; Santini, P.; Caciuffo, R. Unified Crystal-Field Picture for Actinide Dioxides. J. Phys. Chem. Solids 2007, 68, 2020−2023. (13) Prodan, I. D.; Scuseria, G. E.; Sordo, J. A.; Kudin, K. N.; Martin, R. L. Lattice Defects and Magnetic Ordering in Plutonium Oxides: A Hybrid Density-Functional-Theory Study of Strongly Correlated Materials. J. Chem. Phys. 2005, 123, 014703. (14) Wen, X.-D.; Martin, R. L.; Henderson, T. M.; Scuseria, G. E. Density Functional Theory Studies of the Electronic Structure of Solid State Actinide Oxides. Chem. Rev. 2013, 113, 1063−1096. (15) Suzuki, M.-T.; Magnani, N.; Oppeneer, P. M. Microscopic Theory of the Insulating electronic Ground States of the Actinide Dioxides AnO2 (An = U, Np, Pu, Am, and Cm). Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, 195146.

AUTHOR INFORMATION

Corresponding Author

*E-mail: jochena@buffalo.edu. 677

DOI: 10.1021/acs.jpclett.6b02968 J. Phys. Chem. Lett. 2017, 8, 673−678

Letter

The Journal of Physical Chemistry Letters

(36) Mounce, A. M.; Yasuoka, H.; Koutroulakis, G.; Lee, J. A.; Cho, H.; Gendron, F.; Zurek, E.; Scott, B. L.; Trujillo, J. A.; Slemmons, A. K.; Cross, J. N.; Thompson, J. D.; Kozimor, S. A.; Bauer, E. D.; Autschbach, J.; Clark, D. L. Nuclear Magnetic Resonance Measurements and Electronic Structure of Pu(IV) in [(Me)4N]2PuCl6. Inorg. Chem. 2016, 55, 8371−8380. (37) Chibotaru, L. F.; Ungur, L. Ab initio calculation of anisotropic magnetic properties of complexes. I. Unique definition of pseudospin Hamiltonians and their derivation. J. Chem. Phys. 2012, 137, 064112− 22. (38) Aquilante, F.; Autschbach, J.; Carlson, R. K.; Chibotaru, L. F.; Delcey, M. G.; De Vico, L.; et al. Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table. J. Comput. Chem. 2016, 37, 506−541. (39) Zhou, F.; Ozoliņ s.̌ Self-consistent Density Functional Calculations of the Crystal Field Levels in Lanthanide and Actinide Dioxides. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, 075124. (40) Cohen, A. J.; Mori-Sánchez, P.; Yang, W. Insights into current limitations of density functional theory. Science 2008, 321, 792−794. (41) Autschbach, J.; Srebro, M. Delocalization error and ‘functional tuning’ in Kohn-Sham calculations of molecular properties. Acc. Chem. Res. 2014, 47, 2592−2602. (42) Gendron, F.; Le Guennic, B.; Autschbach, J. Magnetic properties and electronic structures of Ar3UIV−L complexes with Ar = C5(CH3)4H−or C5H5− and L = CH3, NO, and Cl. Inorg. Chem. 2014, 53, 13174−13187. (43) Jung, J.; Huang, G.; Calvez, G.; Daiguebonne, C.; Guillou, O.; Cador, O.; Caneschi, A.; Roisnel, T.; Le Guennic, B.; Bernot, K.; Yi, X. Analysis of the electrostatics in DyIII single-molecule magnets: the case study of Dy(Murex)3. Dalton Trans. 2015, 44, 18270−18275. (44) Autillo, M.; Guerin, L.; Bolvin, H.; Moisy, P.; Berthon, C. Magnetic Susceptibility of Actinide(III) Cations: an Experimental and Theoretical Study. Phys. Chem. Chem. Phys. 2016, 18, 6515−6525.

(16) Nakamura, H.; Machida, M.; Kato, M. Effects of Spin-Orbit Coupling and Strong Correlation on the Paramagnetic Insulating State in Plutonium Dioxide. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 82, 155131. (17) Sun, B.; Zhang, P.; Zhao, X.-G. First-Principles Local Density Approximation + U and Generalized Gradient Approxiamtion + U Study of Plutonium Dioxides. J. Chem. Phys. 2008, 128, 084705. (18) Jomard, G.; Amadon, B.; Bottin, F.; Torrent, M. Structural, Thermodynamic, and Electronic Properties of Plutonium Oxides from First Principles. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 78, 075125. (19) Prodan, I. D.; Scuseria, G. E.; Martin, R. L. Assessment of Metageneralized Gradient Approximation and Screened Coulomb Hybrid Density Functionals on Bulk Actinide Oxides. Phys. Rev. B: Condens. Matter Mater. Phys. 2006, 73, 045104. (20) Prodan, I. D.; Scuseria, G. E.; Martin, R. L. Covalency in the Actinide Dioxides: Systematic Study of the Electronic Properties using Screened Hybrid Density Functional Theory. Phys. Rev. B: Condens. Matter Mater. Phys. 2007, 76, 033101. (21) Yin, Q.; Kutepov, A.; Haule, K.; Kotliar, G.; et al. Electronic Correlation and Transport Poperties of Nuclear Fuel Materials. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, 195111. (22) Yin, Q.; Savrasov, S. Y. Orgin of Low Thermal Conductivity in Nuclear Fuels. Phys. Rev. Lett. 2008, 100, 225504. (23) Kolorenč, J.; Shick, A. B.; Lichtenstein, A. I. Electronic Structure and Core-Level Spectra of Light Actinide Dioxides in the Dynamical mean-field Theory. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, 085125. (24) Shick, A. B.; Kolorenč, J.; Havela, L.; Gouder, T.; Caciuffo, R. Nonmagnetic Ground State of PuO2. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 89, 041109. (25) Kovács, A.; Konings, R. J. M.; Gibson, J. K.; Infante, I.; Gagliardi, L. Quantum Chemical Calculations and Experimental Investigations of Molecular Actinide Oxides. Chem. Rev. 2015, 115, 1725−1759. (26) Roos, B. O.; Taylor, P. R.; Siegbahn, P. E. M. A Complete Active Space SCF method (CASSCF) using a density matrix formulated super-CI approach. Chem. Phys. 1980, 48, 157−173. (27) Notter, F.-P.; Dubillard, S.; Bolvin, H. A theoretical study of the ecited states of AmOn+ 2 , n = 1, 2, 3. J. Chem. Phys. 2008, 128, 164315. (28) La Macchia, G.; Infante, I.; Raab, J.; Gibson, J. K.; Gagliardi, L. A theoretical study of the ground state and lowest excited states of PuO0/+/+2 and PuO20/+/+2. Phys. Chem. Chem. Phys. 2008, 10, 7278− 7283. (29) Danilo, C.; Vallet, V.; Flament, J.-P.; Wahlgren, U. Effects of the first hydration sphere and the bulk solvent on the spectra of the f 2 isoelectronic actinide compounds: U4+, NpO+2 , and PuO2+ 2 . Phys. Chem. Chem. Phys. 2010, 12, 1116−1130. (30) Kovács, A.; Infante, I. Theoretical Study of the Electronic Spectra of Neutral and Cationic NpO abd NpO2. J. Chem. Phys. 2015, 143, 074305. (31) Ruipérez, F.; Danilo, C.; Réal, F.; Flament, J.-P.; Vallet, V.; Wahlgren, U. An ab Initio Theoretical Study of the Electronic Structures of UO+2 and [UO2(CO3)3]5−. J. Phys. Chem. A 2009, 113, 1420−1428. (32) Gendron, F.; Páez-Hernández, D.; Notter, F.-P.; Pritchard, B.; Bolvin, H.; Autschbach, J. Magnetic properties and electronic structure of neptunylVI complexes: Wavefunctions, orbitals, and crystal-field models. Chem. - Eur. J. 2014, 20, 7994−8011. (33) Notter, F.-P.; Bolvin, H. Optical and magnetic properties of the series: A theoretical study. J. Chem. Phys. 2009, 130, 5f1 AnXq− 6 184310−11. (34) Chibotaru, L. F.; Ceulemans, A.; Bolvin, H. Unique Definition of the Zeeman-Splitting g Tensor of a Kramers Doublet. Phys. Rev. Lett. 2008, 101, 033003−4. (35) Hernandez, D. P.; Bolvin, H. Magnetic properties of a fourfold degenerate state: Np4+ ion diluted in Cs2 ZrCl6. J. Electron Spectrosc. Relat. Phenom. 2014, 194, 74−80. 678

DOI: 10.1021/acs.jpclett.6b02968 J. Phys. Chem. Lett. 2017, 8, 673−678