Inelastic Neutron Scattering Measurement of the Ground State Tunneling Gap in Tb and Ho Analogues of a Dy Field-Induced Single-Molecule Magnet.

Recent advances in single-molecule magnet (SMM) research have placed great value on interpretation of inelastic neutron scattering (INS) data for rare earth (RE)-containing SMMs. Here, we present the synthesis of several rare earth complexes where combined magnetic and INS studies have been performed, supported by ab initio calculations. The reaction of rare earth nitrate salts with 2,2'-bipyridine (2,2'-bpy) and tetrahalocatecholate (X4Cat2-, X = Br, Cl) ligands in methanol (MeOH) afforded two new families of compounds [RE(2,2'-bpy)2(X4Cat)(X4CatH)(MeOH)] (X = Br and RE = Y, Eu, Gd, Tb, Dy, Ho, Yb for 1-RE; X = Cl and RE = Y, Tb, Dy, Ho, and Yb for 2-RE). Addition of triethylamine (Et3N) to the reaction mixture delivered Et3NH[RE(2,2'-bpy)2(Br4Cat)2] (3-RE, RE = Er and Yb). Interestingly, cerium behaves differently to the rest of the series, generating (2,2'-bpyH)2[Ce(Br4Cat)3(2,2'-bpy)] (4-Ce) with tetravalent Ce(IV) in contrast to the trivalent metal ions in 1-3. The static magnetic properties of 1-RE (RE = Gd, Tb, Dy and Ho) were investigated in conjunction with INS measurements on 1-Y, 1-Tb, and 1-Ho to probe their ground state properties and any crystal field excitations. To facilitate interpretation of the INS spectra and provide insight into the magnetic behavior, ab initio calculations were performed using the single-crystal X-ray diffraction structural data of 1-RE (RE = Tb, Dy and Ho). The ab initio calculations indicate ground doublets dominated by the maximal angular momentum projection states of Kramers type for 1-Dy and Ising type for 1-Tb and 1-Ho. Dynamic magnetic susceptibility measurements indicate that 1-Dy exhibits slow magnetic relaxation in the presence of a small applied magnetic field mainly through Raman pathways. Inelastic neutron scattering spectra exhibit distinct transitions corresponding to crystal field-induced tunneling gaps between the pseudo-doublet ground state components for 1-Tb and 1-Ho, which is one of the first direct experimental measurements with INS of such tunneling transitions in a molecular nanomagnet. The power of high-resolution INS is demonstrated with evidence of two distinct tunneling gaps measurable for the two crystallographically unique Tb coordination environments observed in the single crystal X-ray structure.

Effect of magnetic anisotropy on direct chiral discrimination in paramagnetic NMR spectroscopy.

We present and analyse a complete analytic expression for the generalized shielding polarizability third-rank tensor Φαßγ, whose isotropic average, the pseudoscalar [capital Phi, Greek, macron], is proportional to a chiral macroscopic electric polarization P measurable in a pulsed nuclear magnetic resonance (NMR) experiment probing molecules with an arbitrarily degenerate ground state, in solution. We have recently predicted that [capital Phi, Greek, macron] can be measurable at room temperature in chiral paramagnetic molecules, and identified strong magnetic anistropy (as featured e.g. in lanthanide complexes) as a crucial molecular property to achieve room temperature chiral discrimination using NMR spectroscopy. As previously proposed, the components of Φαßγ are obtained as analytical third derivatives of the electronic free energy. Here we present the explicit calculation of these derivatives, which provide working expressions for the explicit accurate ab initio calculation of Φαßγ. We apply our theory by performing ab initio multiconfigurational calculations of all contributions to Φαßγ, for a set of ten DyIII complexes, characterized by a strongly axial ground Kramers doublet, but also by thermally accessible excited Kramers doublets at room temperature. The results show that the thermally populated excited state contributions, while generally reducing the value of [capital Phi, Greek, macron] calculated on the assumption of a thermally isolated ground state, still confirm the room temperature detectability of this property for all ten studied complexes. Trends on the relative sign of dominant contributions are then discussed on the basis of a crystal field model electrostatic potential splitting a ground spin-orbit multiplet, which provides an insight into the properties of the generalized shielding polarizability tensor for open shell species.

CERES: An ab initio code dedicated to the calculation of the electronic structure and magnetic properties of lanthanide complexes.

We have developed and implemented a new ab initio code, Ceres (Computational Emulator of Rare Earth Systems), completely written in C++11, which is dedicated to the efficient calculation of the electronic structure and magnetic properties of the crystal field states arising from the splitting of the ground state spin-orbit multiplet in lanthanide complexes. The new code gains efficiency via an optimized implementation of a direct configurational averaged Hartree-Fock (CAHF) algorithm for the determination of 4f quasi-atomic active orbitals common to all multi-electron spin manifolds contributing to the ground spin-orbit multiplet of the lanthanide ion. The new CAHF implementation is based on quasi-Newton convergence acceleration techniques coupled to an efficient library for the direct evaluation of molecular integrals, and problem-specific density matrix guess strategies. After describing the main features of the new code, we compare its efficiency with the current state-of-the-art ab initio strategy to determine crystal field levels and properties, and show that our methodology, as implemented in Ceres, represents a more time-efficient computational strategy for the evaluation of the magnetic properties of lanthanide complexes, also allowing a full representation of non-perturbative spin-orbit coupling effects. © 2017 Wiley Periodicals, Inc.

Metal Cluster Electrides: A New Type of Molecular Electride with Delocalised Polyattractor Character.

Electrides are ionic substances containing isolated electrons. These confined electrons are topologically characterised by a quasi-atom, that is, a non-nuclear attractor (NNA) of the electron density. The electronic structure of the octahedral 4 A1g Li6+ and 5 A1g Be6 species shows that these species have a large number of NNAs. These NNAs have highly delocalised electron densities and, as a result, the chemical bonding pattern of these systems is reminiscent of that in solid metals, in which metal cations are surrounded by a "sea" of delocalised valence electrons. We propose the term metal cluster electrides to refer to this new class of compounds. In this study, we establish a computational protocol to identify, characterize, and design metal cluster electrides and we elucidate the intricate bonding patterns of this particular type of species.

Room Temperature Chiral Discrimination in Paramagnetic NMR Spectroscopy.

A recently proposed theory of chiral discrimination in NMR spectroscopy based on the detection of a molecular electric polarization P rotating in a plane perpendicular to the NMR magnetic field [A. D. Buckingham, J. Chem. Phys. 140, 011103 (2014)] is generalized here to paramagnetic systems. Our theory predicts new contributions to P, varying as the square of the inverse temperature. Ab initio calculations for ten Dy^{3+} complexes, at 293 K, show that, in strongly anisotropic paramagnetic molecules, P can be more than 1000 times larger than in diamagnetic molecules, making paramagnetic NMR chiral discrimination amenable to room temperature detection.

Configuration-averaged 4f orbitals in ab initio calculations of low-lying crystal field levels in lanthanide(iii) complexes.

A successful and commonly used ab initio method for the calculation of crystal field levels and magnetic anisotropy of lanthanide complexes consists of spin-adapted state-averaged CASSCF calculations followed by state interaction with spin-orbit coupling (SI-SO). Based on two observations valid for Ln(iii) complexes, namely: (i) CASSCF 4f orbitals are expected to change very little when optimized for different states belonging to the 4f electronic configuration, (ii) due to strong spin-orbit coupling the total spin is not a good quantum number, we show here via a straightforward analysis and direct calculation that the CASSCF/SI-SO method can be simplified to a single configuration-averaged HF calculation and one complete active space CI diagonalization, including spin-orbit coupling, on determinant basis. Besides its conceptual simplicity, this approach has the advantage that all spin states of the 4f(n) configuration are automatically included in the SO coupling, thereby overcoming one of the computational limitations of the existing CASSCF/SI-SO approach. As an example, we consider three isostructural complexes [Ln(acac)3(H2O)2], Ln = Dy(3+), Ho(3+), Er(3+), and find that the proposed simplified method yields crystal field levels and magnetic g-tensors that are in very good agreement with those obtained with CASSCF/SI-SO.

Ab initio-based determination of lanthanoid-radical exchange as visualised by inelastic neutron scattering.

Magnetic exchange coupling can modulate the slow magnetic relaxation in single-molecule magnets. Despite this, elucidation of exchange coupling remains a significant challenge for the lanthanoid(iii) ions, both experimentally and computationally. In this work, the crystal field splitting and 4f-π exchange coupling in the erbium-semiquinonate complex [ErTp2dbsq] (Er-dbsq; Tp- = hydro-tris(1-pyrazolyl)borate, dbsqH2 = 3,5-di-tert-butyl-1,2-semiquinone) have been determined by inelastic neutron scattering (INS), magnetometry, and CASSCF-SO ab initio calculations. A related complex with a diamagnetic ligand, [ErTp2trop] (Er-trop; tropH = tropolone), has been used as a model for the crystal field splitting in the absence of coupling. Magnetic and INS data indicate antiferromagnetic exchange for Er-dbsq with a coupling constant of Jex = -0.23 meV (-1.8 cm-1) (-2Jex formalism) and good agreement is found between theory and experiment, with the low energy magnetic and spectroscopic properties well modelled. Most notable is the ability of the ab initio modelling to reproduce the signature of interference between localised 4f states and delocalised π-radical states that is evident in the Q-dependence of the exchange excitation. This work highlights the power of combining INS with EPR and magnetometry for determination of ground state properties, as well as the enhanced capability of CASSCF-SO ab initio calculations and purposely developed ab initio-based theoretical models. We deliver an unprecedentedly detailed representation of the entangled character of 4f-π exchange states, which is obtained via an accurate image of the spin-orbital transition density between the 4f-π exchange coupled wavefunctions.

Direct observation of magnetoelastic coupling in a molecular spin qubit: new insights from crystal field neutron scattering data.

Single-molecule magnets are promising candidates for data storage and quantum computing applications. A major barrier to their use is rapid magnetic relaxation and quantum decoherence due to thermal vibrations. Here we report a reanalysis of inelastic neutron scattering (INS) data of the candidate qubit Na9[Ho(W5O18)2]·35D2O, wherein we demonstrate for the first time that magnetic relaxation times and mechanisms can be directly observed as crystal field (CF) peak broadening in INS spectra of a lanthanoid molecular system. The magnetoelastic coupling between the lower energy CF states and phonons (lattice vibrations) is determined by the simultaneous measurement of CF excitations and the phonon density of states, encoded within the same INS experiment. This directly results in the determination of relaxation coupling pathways that occur in this molecule. Such information is invaluable for the further advancement of SMMs and to date has only been obtained from techniques performed in external magnetic fields. Additionally, we determine a relaxation rate of quantum-tunnelling of magnetisation that is consistent with previously measured EPR spectroscopy data.