Ab Initio Design of Molecular Qubits with Electric Field Control.

Current scalable quantum computers require large footprints and complex interconnections due to the design of superconducting qubits. While this architecture is competitive, molecular qubits offer a promising alternative due to their atomic scale and tuneable properties through chemical design. The use of electric fields to precisely, selectively and coherently manipulate molecular spins with resonant pulses has the potential to solve the experimental limitations of current molecular spin manipulation techniques such as electron paramagnetic resonance (EPR) spectroscopy. EPR can only address a macroscopic ensemble of molecules, defeating the inherent benefits of molecule-based quantum information. Hence, numerous experiments have been performed using EPR in combination with electric fields to demonstrate coherent spin manipulation. In this work, we explore the underlying theory of spin-electric coupling in lanthanide molecules, and outline ab initio methods to design molecules with enhanced electric field responses. We show how structural distortions arising from electric fields generate coupling elements in the crystal field Hamiltonian within a Kramers doublet ground state and demonstrate the impact of molecular geometry on this phenomenon. We use perturbation theory to rationalize the magnetic and electric field orientation dependence of the spin-electric coupling. We use pseudo-symmetry point groups to decompose molecular distortions to understand the role that symmetry has on spin-electric coupling. Finally, we present an analytical electric field model of structural perturbations that provides large savings in computational expense and allows for the investigation of experimentally accessible electric field magnitudes which cannot be accessed using common ab initio methods.

1.3 Å Crystal Structure of E. coli Peptidyl-Prolyl Isomerase B with Uniform Substitution of Valine by (2S,3S)-4-Fluorovaline Reveals Structure Conservation and Multiple Staggered Rotamers of CH2F Groups.

(2S,3S)-4-Fluorovaline (FVal) is an analogue of valine, where a single CH3 group is substituted by a CH2F group. In the absence of valine, E. coli valyl-tRNA synthetase uses FVal as a substitute, enabling the production of proteins uniformly labeled with FVal. Here, we describe the production and analysis of E. coli peptidyl-prolyl isomerase B where all 16 valine residues have been replaced by FVal synthesized with a 13C-labeled CH2F group. Although the melting temperature is lower by about 11 °C relative to the wild-type protein, the three-dimensional protein structure is almost completely conserved, as shown by X-ray crystallography. The CH2F groups invariably populate staggered rotamers. Most CH2F groups populate two different rotamers. The increased space requirement of fluorine versus hydrogen does not prohibit rotamers that position fluorine next to a backbone carbonyl carbon. 19F NMR spectra show a signal dispersion over 25 ppm. The most high-field shifted 19F resonances correlate with large 3JHF coupling constants, confirming the impact of the γ-gauche effect on the signal dispersion. The present work is the second experimental verification of the effect and extends its validity to fluorovaline. The abundance of valine in proteins and structural conservation with FVal renders this valine analogue attractive for probing proteins by 19F NMR spectroscopy.

δ-Bonding modulates the electronic structure of formally divalent nd1 rare earth arene complexes.

Landmark advances in rare earth (RE) chemistry have shown that divalent complexes can be isolated with non-Aufbau 4f n {5d/6s}1 electron configurations, facilitating remarkable bonding motifs and magnetic properties. We report a series of divalent bis-tethered arene complexes, [RE(NHAriPr6 )2] (2RE; RE = Sc, Y, La, Sm, Eu, Tm, Yb; NHAriPr6 = {N(H)C6H3-2,6-(C6H2-2,4,6-iPr3)2}). Fluid solution EPR spectroscopy gives g iso < 2.002 for 2Sc, 2Y, and 2La, consistent with formal nd1 configurations, calculations reveal metal-arene δ-bonding via mixing of nd(x 2-y 2) valence electrons into arene π* orbitals. Experimental and calculated EPR and UV-Vis-NIR spectroscopic properties for 2Y show that minor structural changes markedly alter the metal d(x 2-y 2) contribution to the SOMO. This contrasts 4f n {5d/6s}1 complexes where the valence d-based electron resides in a non-bonding orbital. Complexes 2Sm, 2Eu, 2Tm, and 2Yb contain highly-localised 4f n+1 ions with no appreciable metal-arene bonding by density functional calculations. These results show that the physicochemical properties of divalent rare earth arene complexes with both formal nd1 and 4f n+1 configurations are nuanced, may be controlled through ligand modification, and require a multi-pronged experimental and theoretical approach to fully rationalise.

Reactivity of a triamidoamine terminal uranium(VI)-nitride with 3d-transition metal metallocenes.

Reactions between [(TrenTIPS)UVIN] (1, TrenTIPS = {N(CH2CH2NSiPri3)3}3-) and [MII(η5-C5R5)2] (M/R = Cr/H, Mn/H, Fe/H, Ni/H) were intractable, but M/R = Co/H or Co/Me afforded [(TrenTIPS)UVN-(η1:η4-C5H5)CoI(η5-C5H5)] (2) and [(TrenTIPS)UIV-NH2] (3), respectively. For M/R = V/H [(TrenTIPS)UIV-NVIV(η5-C5H5)2] (4), was isolated. Complexes 2-4 evidence one-/two-electron uranium reductions, nucleophilic nitrides, and partial N-atom transfer.

Tris-Silanide f-Block Complexes: Insights into Paramagnetic Influence on NMR Chemical Shifts.

The paramagnetism of f-block ions has been exploited in chiral shift reagents and magnetic resonance imaging, but these applications tend to focus on 1H NMR shifts as paramagnetic broadening makes less sensitive nuclei more difficult to study. Here we report a solution and solid-state (ss) 29Si NMR study of an isostructural series of locally D 3h -symmetric early f-block metal(III) tris-hypersilanide complexes, [M{Si(SiMe3)3}3(THF)2] (1-M; M = La, Ce, Pr, Nd, U); 1-M were also characterized by single crystal and powder X-ray diffraction, EPR, ATR-IR, and UV-vis-NIR spectroscopies, SQUID magnetometry, and elemental analysis. Only one SiMe3 signal was observed in the 29Si ssNMR spectra of 1-M, while two SiMe3 signals were seen in solution 29Si NMR spectra of 1-La and 1-Ce. This is attributed to dynamic averaging of the SiMe3 groups in 1-M in the solid state due to free rotation of the M-Si bonds and dissociation of THF from 1-M in solution to give the locally C 3v -symmetric complexes [M{Si(SiMe3)3}3(THF) n ] (n = 0 or 1), which show restricted rotation of M-Si bonds on the NMR time scale. Density functional theory and complete active space self-consistent field spin-orbit calculations were performed on 1-M and desolvated solution species to model paramagnetic NMR shifts. We find excellent agreement of experimental 29Si NMR data for diamagnetic 1-La, suggesting n = 1 in solution and reasonable agreement of calculated paramagnetic shifts of SiMe3 groups for 1-M (M = Pr and Nd); the NMR shifts for metal-bound 29Si nuclei could only be reproduced for diamagnetic 1-La, showing the current limitations of pNMR calculations for larger nuclei.

Intramolecular bridging strategies to suppress two-phonon Raman spin relaxation in dysprosocenium single-molecule magnets.

Dy(III) bis-cyclopentadienyl (Cp) sandwich compounds exhibit extremely strong single-ion magnetic anisotropy which imbues them with magnetic memory effects such as magnetic hysteresis, and has put them at the forefront of high-performance single-molecule magnets (SMMs). Owing to the great success of design principles focused on maximising the anisotropy barrier, ever higher Ueff values have been reported leading to significant slow down of single-phonon Orbach spin relaxation. However, anisotropy-based SMM design has largely ignored two-phonon Raman spin relaxation, which is still limiting the temperatures at which a memory effect can be observed. In this work, we study the suppression of Raman relaxation through covalent bridging of the Cp ligands by alkyl chains, testing the hypothesis that increasing the rigidity of the ligand framework results in a blue shift of low frequency vibrations in the first coordination sphere of the Dy(III) ion. This reshaping of the vibrational low-energy density of states (DOS) results in lower occupation of pseudo-acoustic phonons available to drive Raman relaxation at low temperatures. We simulate Orbach and Raman spin relaxation in a series of zero-, mono-, di- and tri-bridged [Dy(Cpttt)2]+ analogues fully ab initio, using a quantum mechanics (QM)/molecular mechanics (MM) condensed phase embedding protocol in a periodic solvent matrix as a generic and experimentally testable environment model that can include (pseudo-)acoustic phononic degrees of freedom. We show that this approach can simulate magnetic relaxation dynamics in the condensed phase for the existing non-bridged [Dy(Cpttt)2]+ compound with quantitative experimental accuracy. Subsequently, we find a significant slowing down of Raman relaxation can be achieved for the singly-bridged SMM, while the introduction of further bridges leads to faster relaxation. A key result being that we find the two-phonon Raman rates correlate with the purity of the first-excited Kramers doublet in terms of its mJ = ±13/2 content. Even though the bridging design principle is successful at progressively reshaping the low-energy DOS, the introduction of linker atoms in the equatorial plane successively degrades magnetic anisotropy, suggesting the importance of refined design of the linker chemistry. The accuracy of our results emphasises the value of a generic periodic solvent embedding model, such that it permits the modelling of molecular spin dynamics in the condensed phase without knowledge of a crystal structure. This allows the study of hypothetical molecules or aggregates under real-world conditions, which we expect to have utility beyond the field of molecular magnetism.

Coercive Fields Exceeding 30 T in the Mixed-Valence Single-Molecule Magnet (CpiPr5)2Ho2I3.

Mixed-valence dilanthanide complexes of the type (CpiPr5)2Ln2I3 (CpiPr5 = pentaisopropylcyclopentadienyl; Ln = Gd, Tb, Dy) featuring a direct Ln-Ln σ-bonding interaction have been shown to exhibit well-isolated high-spin ground states and, in the case of the Tb and Dy variants, a strong axial magnetic anisotropy that gives rise to a large magnetic coercivity. Here, we report the synthesis and characterization of two new mixed-valence dilanthanide compounds in this series, (CpiPr5)2Ln2I3 (1-Ln; Ln = Ho, Er). Both compounds feature a Ln-Ln bonding interaction, the first such interaction in any molecular compounds of Ho or Er. Like the Tb and Dy congeners, both complexes exhibit high-spin ground states arising from strong spin-spin coupling between the lanthanide 4f electrons and a single σ-type lanthanide-lanthanide bonding electron. Beyond these similarities, however, the magnetic properties of the two compounds diverge. In particular, 1-Er does not exhibit observable magnetic blocking or slow magnetic relaxation, while 1-Ho exhibits magnetic blocking below 28 K, which is the highest temperature among Ho-based single-molecule magnets, and a spin reversal barrier of 556(4) cm-1. Additionally, variable-field magnetization data collected for 1-Ho reveal a coercive field of greater than 32 T below 8 K, more than 6-fold higher than observed for the bulk magnets SmCo5 and Nd2Fe14B, and the highest coercive field reported to date for any single-molecule magnet or molecule-based magnetic material. Multiconfigurational calculations, supported by far-infrared magnetospectroscopy data, reveal that the stark differences in magnetic properties of 1-Ho and 1-Er arise from differences in the local magnetic anisotropy of the lanthanide centers.

Ligand Effects on the Spin Relaxation Dynamics and Coherent Manipulation of Organometallic La(II) Potential Qudits.

We present pulsed electron paramagnetic resonance (EPR) studies on three La(II) complexes, [K(2.2.2-cryptand)][La(Cp')3] (1), [K(2.2.2-cryptand)][La(Cp″)3] (2), and [K(2.2.2-cryptand)][La(Cptt)3] (3), which feature cyclopentadienyl derivatives as ligands [Cp' = C5H4SiMe3; Cp″ = C5H3(SiMe3)2; Cptt = C5H3(CMe3)2] and display a C3 symmetry. Long spin-lattice relaxation (T1) and phase memory (Tm) times are observed for all three compounds, but with significant variation in T1 among 1-3, with 3 being the slowest relaxing due to higher s-character of the SOMO. The dephasing times can be extended by more than an order of magnitude via dynamical decoupling experiments using a Carr-Purcell-Meiboom-Gill (CPMG) sequence, reaching 161 µs (5 K) for 3. Coherent spin manipulation is performed by the observation of Rabi quantum oscillations up to 80 K in this nuclear spin-rich environment (1H, 13C, and 29Si). The high nuclear spin of 139La (I = 7/2), and the ability to coherently manipulate all eight hyperfine transitions, makes these molecules promising candidates for application as qudits (multilevel quantum systems featuring d quantum states; d >2) for performing quantum operations within a single molecule. Application of HYSCORE techniques allows us to quantify the electron spin density at ligand nuclei and interrogate the role of functional groups to the electron spin relaxation properties.

Relativistic Quantum Chemical Investigation of Actinide Covalency Measured by Electron Paramagnetic Resonance Spectroscopy.

We investigate actinide covalency effects in two [AnCptt3] (An = Th, U) complexes recently studied with pulsed electron paramagnetic resonance spectroscopy, using the Hyperion package to obtain relativistic hyperfine coupling constants from relativistic multiconfigurational wave functions. 1H and 13C HYSCORE simulations using the computed parameters show excellent agreement with the experimental data, highlighting the accuracy of modern relativistic ab initio methods. The extent of covalency indicated from the calculations on [ThCptt3] is in agreement with the original report based on traditional spectral fitting methods, while the covalency in [UCptt3] is found to be previously overestimated. The latter is due to the paramagnetic spin-orbit effect that arises naturally in a relativistic theory of hyperfine coupling and yet was not accounted for in the original study, thus highlighting the necessity of relativistic approaches for the interpretation of magnetic resonance data pertaining to actinides.

Halobenzene Adducts of a Dysprosocenium Single-Molecule Magnet.

Dysprosium complexes with strong axial crystal fields are promising candidates for single-molecule magnets (SMMs), which could be used for high-density data storage. Isolated dysprosocenium cations, [Dy(CpR)2]+ (CpR = substituted cyclopentadienyl), have recently shown magnetic hysteresis (a memory effect) above the temperature of liquid nitrogen. Synthetic efforts have focused on reducing strong transverse ligand fields in these systems as they are known to enhance magnetic relaxation by spin-phonon mechanisms. Here we show that equatorial coordination of the halobenzenes PhX (X = F, Cl, Br) and o-C6H4F2 to the cation of a recently reported dysprosocenium complex [Dy(Cpttt)(Cp*)][Al{OC(CF3)3}4] (Cpttt = C5H2tBu3-1,2,4; Cp* = C5Me5) reduces magnetic hysteresis temperatures compared to that of the parent cation. We find that this is due to increased effectiveness of both one- (Orbach) and two-phonon (Raman) relaxation mechanisms, which correlate with the electronegativity and number of interactions with the halide despite κ1-coordination of a single halobenzene having a minimal effect on the metrical parameters of [Dy(Cpttt)(Cp*)(PhX-κ1-X)]+ cations vs the isolated [Dy(Cpttt)(Cp*)]+ cation. We observe unusual divergent behavior of relaxation rates at low temperatures in [Dy(Cpttt)(Cp*)(PhX)][Al{OC(CF3)3}4], which we attribute to a phonon bottleneck effect. We find that, despite the transverse fields introduced by the monohalobenzenes in these cations, the interactions are sufficiently weak that the effective barriers to magnetization reversal remain above 1000 cm-1, being only ca. 100 cm-1 lower than for the parent complex, [Dy(Cpttt)(Cp*)][Al{OC(CF3)3}4].

Synthesis and Magnetic Properties of Bis-Halobenzene Decamethyldysprosocenium Cations.

The decamethyldysprosocenium cation, [Dy(Cp*)2]+ (Cp* = {C5Me5}), was a target single-molecule magnet (SMM) prior to the isolation of larger dysprosocenium cations, which have recently shown magnetic memory effects up to 80 K. However, the relatively short Dy···Cp*centroid distances of [Dy(Cp*)2]+, together with the reduced resonance of its vibrational modes with electronic states compared to larger dysprosocenium cations, could lead to more favorable SMM behavior. Here, we report the synthesis and magnetic properties of a series of solvated adducts containing bis-halobenzene decamethyldysprosocenium cations, namely [Dy(Cp*)2(PhX-κ-X)2][Al{OC(CF3)3}4] (X = F or Cl) and [Dy(Cp*)2(C6H4F2-κ2-F,F)(C6H4F2-κ-F)][Al{OC(CF3)3}4]. These complexes were prepared by the sequential reaction of [Dy(Cp*)2(µ-BH4)]∞ with allylmagnesium chloride and [NEt3H][Al{OC(CF3)3}4], followed by recrystallization from parent halobenzenes. The complexes were characterized by powder and single crystal X-ray diffraction, NMR and ATR-IR spectroscopy, elemental analysis, and SQUID magnetometry; experimental data were rationalized by a combination of density functional theory and ab initio calculations. We find that bis-halobenzene adducts of the [Dy(Cp*)2]+ cation exhibit highly bent Cp*···Dy···Cp* angles; these cations are also susceptible to decomposition by C-X (X = F, Cl, Br) activation and displacement of halobenzenes by O-donor ligands. The effective energy barrier to reversal of magnetization measured for [Dy(Cp*)2(PhF-κ-F)2][Al{OC(CF3)3}4] (930(6) cm-1) sets a new record for SMMs containing {Dy(Cp*)2} fragments, though all SMM parameters are lower than would be predicted for an isolated [Dy(Cp*)2]+ cation, as expected due to transverse ligand fields introduced by halobenzenes and the large deviation of the Cp*···Dy···Cp* angle from linearity promoting magnetic relaxation.

Carbene Complexes of Plutonium: Structure, Bonding, and Divergent Reactivity to Lanthanide Analogs.

Organoplutonium chemistry was established in 1965, yet structurally authenticated plutonium-carbon bonds remain rare being limited to π-bonded carbocycle and σ-bonded isonitrile and hydrocarbyl derivatives. Thus, plutonium-carbenes, including alkylidenes and N-heterocyclic carbenes (NHCs), are unknown. Here, we report the preparation and characterization of the diphosphoniomethanide-plutonium complex [Pu(BIPMTMSH)(I)(µ-I)]2 (1Pu, BIPMTMSH = (Me3SiNPPh2)2CH) and the diphosphonioalkylidene-plutonium complexes [Pu(BIPMTMS)(I)(DME)] (2Pu, BIPMTMS = (Me3SiNPPh2)2C) and [Pu(BIPMTMS)(I)(IMe4)2] (3Pu, IMe4 = C(NMeCMe)2), thus disclosing non-actinyl transneptunium multiple bonds and transneptunium NHC complexes. These Pu-C double and dative bonds, along with cerium, praseodymium, samarium, uranium, and neptunium congeners, enable lanthanide-actinide and actinide-actinide comparisons between metals with similar ionic radii and isoelectronic 4f5 vs 5f5 electron-counts within conserved ligand fields over 12 complexes. Quantum chemical calculations reveal that the orbital-energy and spatial-overlap terms increase from uranium to neptunium; however, on moving to plutonium the orbital-energy matching improves but the spatial overlap decreases. The bonding picture that emerges is more complex than the traditional picture of the bonding of lanthanides being ionic and early actinides being more covalent but becoming more ionic left to right. Multiconfigurational calculations on 2M and 3M (M = Pu, Sm) account for the considerably more complex UV/vis/NIR spectra for 5f5 2Pu and 3Pu compared to 4f5 2Sm and 3Sm. Supporting the presence of Pu═C double bonds in 2Pu and 3Pu, 2Pu exhibits metallo-Wittig bond metathesis involving the highest atomic number element to date, reacting with benzaldehyde to produce the alkene PhC(H)═C(PPh2NSiMe3)2 (4) and "PuOI". In contrast, 2Ce and 2Pr do not react with benzaldehyde to produce 4.

Metal-carbon bonding in early lanthanide substituted cyclopentadienyl complexes probed by pulsed EPR spectroscopy.

We examine lanthanide (Ln)-ligand bonding in a family of early Ln3+ complexes [Ln(Cptt)3] (1-Ln, Ln = La, Ce, Nd, Sm; Cptt = C5H3tBu2-1,3) by pulsed electron paramagnetic resonance (EPR) methods, and provide the first characterization of 1-La and 1-Nd by single crystal XRD, multinuclear NMR, IR and UV/Vis/NIR spectroscopy. We measure electron spin T1 and Tm relaxation times of 12 and 0.2 µs (1-Nd), 89 and 1 µs (1-Ce) and 150 and 1.7 µs (1-Sm), respectively, at 5 K: the T1 relaxation of 1-Nd is more than 102 times faster than its valence isoelectronic uranium analogue. 13C and 1H hyperfine sublevel correlation (HYSCORE) spectroscopy reveals that the extent of covalency is negligible in these Ln compounds, with much smaller hyperfine interactions than observed for equivalent actinide (Th and U) complexes. This is corroborated by ab initio calculations, confirming the predominant electrostatic nature of the metal-ligand bonding in these complexes.

Isolation of a Bent Dysprosium Bis(amide) Single-Molecule Magnet.

The isolation of formally two-coordinate lanthanide (Ln) complexes is synthetically challenging, due to predominantly ionic Ln bonding regimes favoring high coordination numbers. In 2015, it was predicted that a near-linear dysprosium bis(amide) cation [Dy{N(SiiPr3)2}2]+ could provide a single-molecule magnet (SMM) with an energy barrier to magnetic reversal (Ueff) of up to 2600 K, a 3-fold increase of the record Ueff for a Dy SMM at the time; this work showed a potential route to SMMs that can provide high-density data storage at higher temperatures. However, synthetic routes to a Dy complex containing only two monodentate ligands have not previously been realized. Here, we report the synthesis of the target bent dysprosium bis(amide) complex, [Dy{N(SiiPr3)2}2][Al{OC(CF3)3}4] (1-Dy), together with the diamagnetic yttrium analogue. We find Ueff = 950 ± 30 K for 1-Dy, which is much lower than the predicted values for idealized linear two-coordinate Dy(III) cations. Ab initio calculations of the static electronic structure disagree with the experimentally determined height of the Ueff barrier, thus magnetic relaxation is faster than expected based on magnetic anisotropy alone. We propose that this is due to enhanced spin-phonon coupling arising from the flexibility of the Dy coordination sphere, in accord with ligand vibrations being of equal importance to magnetic anisotropy in the design of high-temperature SMMs.

Vibronic effects on the quantum tunnelling of magnetisation in Kramers single-molecule magnets.

Single-molecule magnets are among the most promising platforms for achieving molecular-scale data storage and processing. Their magnetisation dynamics are determined by the interplay between electronic and vibrational degrees of freedom, which can couple coherently, leading to complex vibronic dynamics. Building on an ab initio description of the electronic and vibrational Hamiltonians, we formulate a non-perturbative vibronic model of the low-energy magnetic degrees of freedom in monometallic single-molecule magnets. Describing their low-temperature magnetism in terms of magnetic polarons, we are able to quantify the vibronic contribution to the quantum tunnelling of the magnetisation, a process that is commonly assumed to be independent of spin-phonon coupling. We find that the formation of magnetic polarons lowers the tunnelling probability in both amorphous and crystalline systems by stabilising the low-lying spin states. This work, thus, shows that spin-phonon coupling subtly influences magnetic relaxation in single-molecule magnets even at extremely low temperatures where no vibrational excitations are present.

Accurate and Efficient Spin-Phonon Coupling and Spin Dynamics Calculations for Molecular Solids.

Molecular materials are poised to play a significant role in the development of future optoelectronic and quantum technologies. A crucial aspect of these areas is the role of spin-phonon coupling and how it facilitates energy transfer processes such as intersystem crossing, quantum decoherence, and magnetic relaxation. Thus, it is of significant interest to be able to accurately calculate the molecular spin-phonon coupling and spin dynamics in the condensed phase. Here, we demonstrate the maturity of ab initio methods for calculating spin-phonon coupling by performing a case study on a single-molecule magnet and showing quantitative agreement with the experiment, allowing us to explore the underlying origins of its spin dynamics. This feat is achieved by leveraging our recent developments in analytic spin-phonon coupling calculations in conjunction with a new method for including the infinite electrostatic potential in the calculations. Furthermore, we make the first ab initio determination of phonon lifetimes and line widths for a molecular magnet to prove that the commonplace Born-Markov assumption for the spin dynamics is valid, but such "exact" phonon line widths are not essential to obtain accurate magnetic relaxation rates. Calculations using this approach are facilitated by the open-source packages we have developed, enabling cost-effective and accurate spin-phonon coupling calculations on molecular solids.

Thermally Stable Terbium(II) and Dysprosium(II) Bis-amidinate Complexes.

The thermostable four-coordinate divalent lanthanide (Ln) bis-amidinate complexes [Ln(Piso)2] (Ln = Tb, Dy; Piso = {(NDipp)2CtBu}, Dipp = C6H3iPr2-2,6) were prepared by the reduction of parent five-coordinate Ln(III) precursors [Ln(Piso)2I] (Ln = Tb, Dy) with KC8; halide abstraction of [Ln(Piso)2I] with [H(SiEt3)2][B(C6F5)] gave the respective Ln(III) complexes [Ln(Piso)2][B(C6F5)]. All complexes were characterized by X-ray diffraction, ICP-MS, elemental analysis, SQUID magnetometry, UV-vis-NIR, ATR-IR, NMR, and EPR spectroscopy and ab initio CASSCF-SO calculations. These data consistently show that [Ln(Piso)2] formally exhibit Ln(II) centers with 4fn5dz21 (Ln = Tb, n = 8; Dy, n = 9) valence electron configurations. We show that simple assignments of the f-d coupling to either L-S or J-s schemes are an oversimplification, especially in the presence of significant crystal field splitting. The coordination geometry of [Ln(Piso)2] is intermediate between square planar and tetrahedral. Projecting from the quaternary carbon atoms of the CN2 ligand backbones shows near-linear C···Ln···C arrangements. This results in strong axial ligand fields to give effective energy barriers to magnetic reversal of 1920(91) K for the Tb(II) analogue and 1964(48) K for Dy(II), the highest values observed for mononuclear Ln(II) single-molecule magnets, eclipsing 1738 K for [Tb(C5iPr5)2]. We tentatively attribute the fast zero-field magnetic relaxation for these complexes at low temperatures to transverse fields, resulting in considerable mixing of mJ states.

Correction: Synthesis and characterization of heterometallic rings templated through alkylammonium or imidazolium cations.

Correction for 'Synthesis and characterization of heterometallic rings templated through alkylammonium or imidazolium cations' by Rajeh Alotaibi et al., Dalton Trans., 2023, 52, 7473-7481, https://doi.org/10.1039/D3DT00982C.

AtomAccess: A Predictive Tool for Molecular Design and Its Application to the Targeted Synthesis of Dysprosium Single-Molecule Magnets.

Isolated dysprosocenium cations, [Dy(CpR)2]+ (CpR = substituted cyclopentadienyl), have recently been shown to exhibit superior single-molecule magnet (SMM) properties over closely related complexes with equatorially bound ligands. However, gauging the crossover point at which the CpR substituents are large enough to prevent equatorial ligand binding, but small enough to approach the metal closely and generate strong crystal field splitting has required laborious synthetic optimization. We therefore created the computer program AtomAccess to predict the accessibility of a metal binding site and its ability to accommodate additional ligands. Here, we apply AtomAccess to identify the crossover point for equatorial coordination in [Dy(CpR)2]+ cations in silico and hence predict a cation that is at the cusp of stability without equatorial interactions, viz., [Dy(Cpttt)(Cp*)]+ (Cpttt = C5H2tBu3-1,2,4, Cp* = C5Me5). Upon synthesizing this cation, we found that it crystallizes as either a contact ion-pair, [Dy(Cpttt)(Cp*){Al[OC(CF3)3]4-κ-F}], or separated ion-pair polymorph, [Dy(Cpttt)(Cp*)][Al{OC(CF3)3}4]·C6H6. Upon characterizing these complexes, together with their precursors, yttrium and yttrium-doped analogues, we find that the contact ion-pair shows inferior SMM properties to the separated ion-pair, as expected, due to faster Raman and quantum tunneling of magnetization relaxation processes, while the Orbach region is relatively unaffected. The experimental verification of the predicted crossover point for equatorial coordination in this work tests the limitations of the use of AtomAccess as a predictive tool and also indicates that the application of this type of program shows considerable potential to boost efficiency in exploratory synthetic chemistry.

Structural Evolution of Paramagnetic Lanthanide Compounds in Solution Compared to Time- and Ensemble-Average Structures.

Anisotropy in the magnetic susceptibility strongly influences the paramagnetic shifts seen in nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI) experiments. A previous study on a series of C3-symmetric prototype MRI contrast agents showed that their magnetic anisotropy was highly sensitive to changes in molecular geometry and concluded that changes in the average angle between the lanthanide-oxygen (Ln-O) bonds and the molecular C3 axis due to solvent interactions had a significant impact on the magnetic anisotropy and, consequently, the paramagnetic shift. However, this study, like many others, was predicated on an idealized C3-symmetric structural model, which may not be representative of the dynamic structure in solution at the single-molecule level. Here, we address this by using ab initio molecular dynamics simulations to simulate how the molecular geometry, in particular the angles between the Ln-O bonds and the pseudo-C3 axis, evolves over time in the solution, mimicking typical experimental conditions. We observe large-amplitude oscillations in the O-Ln-C̃3 angles, and complete active space self-consistent field spin-orbit calculations show that this leads to similarly large oscillations in the pseudocontact (dipolar) paramagnetic NMR shifts. The time-averaged shifts show good agreement with experimental measurements, while the large fluctuations suggest that an idealized structure provides an incomplete description of the solution dynamics. Our observations have significant implications for modeling the electronic and nuclear relaxation times in this and other systems where the magnetic susceptibility is exquisitely sensitive to the molecular structure.