skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Origin of Magnetic Anisotropy in Nickelocene Molecular Magnet and Resilience of Its Magnetic Behavior
The robustness of nickelocene’s (NiCp2, Cp = cyclopentadienyl) magnetic anisotropy and addressability of its spin states make this molecular magnet attractive as a spin sensor. However, microscopic understanding of its magnetic anisotropy is still lacking, especially when NiCp2 is deposited on a surface to make quantum sensing devices. Quantum chemical calculations of such molecule/solid-state systems are limited to density functional theory (DFT) or DFT+U (Hubbard correction to DFT). We investigate the magnetic behavior of NiCp2 using the spin-flip variant of the equation-of-motion coupled-cluster (EOM-SF-CC) method and use the EOM-SF-CC results to benchmark SF-TD-DFT. Our first-principle calculations agree well with experimentally derived magnetic anisotropy and susceptibility values. The calculations show that magnetic anisotropy in NiCp2 originates from a large spin–orbit coupling (SOC) between the triplet ground state and the third singlet state, whereas the coupling with lower singlet excited states is negligible. We also considered a set of six ring-substituted NiCp2 derivatives and a model system of the NiCp2/MgO(001) adsorption complex, for which we used SF-TD-DFT method. To gain insight into the electronic structure of these systems, we analyze spinless transition density matrices and their natural transition orbitals (NTOs). The NTO analysis of SOCs explains how spin states and magnetic properties are retained upon modification of the NiCp2 coordination environment and upon its adsorption on a surface. Such resilience of the NiCp2 magnetic behavior supports using NiCp2 as a spin-probe molecule by functionalization of the tip of a scanning tunneling microscope.  more » « less
Award ID(s):
2221453
PAR ID:
10534618
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Chemical Society
Date Published:
Journal Name:
Journal of physical chemistry C
ISSN:
1932-7447
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We present a new implementation for computing spin–orbit couplings (SOCs) within a time-dependent density-functional theory (TD-DFT) framework in the standard spin-conserving formulation as well in the spin–flip variant (SF-TD-DFT). This approach employs the Breit–Pauli Hamiltonian and Wigner–Eckart’s theorem applied to the reduced one-particle transition density matrices, together with the spin–orbit mean-field treatment of the two-electron contributions. We use a state-interaction procedure and compute the SOC matrix elements using zero-order non-relativistic states. Benchmark calculations using several closed-shell organic molecules, diradicals, and a single-molecule magnet illustrate the efficiency of the SOC protocol. The results for organic molecules (described by standard TD-DFT) show that SOCs are insensitive to the choice of the functional or basis sets, as long as the states of the same characters are compared. In contrast, the SF-TD-DFT results for small diradicals (CH 2 , [Formula: see text], SiH 2 , and [Formula: see text]) show strong functional dependence. The spin-reversal energy barrier in a Fe(III) single-molecule magnet computed using non-collinear SF-TD-DFT (PBE0, ωPBEh/cc-pVDZ) agrees well with the experimental estimate. 
    more » « less
  2. Heteroaromatic species are commonly found in complex gaseous mixtures, from tobacco smoke to petroleum and asphaltene combustion products. At high temperatures, C–H bond rupture produces various dehydro radical isomers. We have used the spin–flip formulation of equation-of-motion coupled cluster theory with single and double substitutions (EOM-SF-CCSD) to characterize the energies and wave functions of the lowest lying singlet and triplet states of the diradical (2,3), (2,4), (2,5), and (3,4) di-dehydro isomers of pyrrole, furan, and thiophene. In all cases, these diradicals are minima on the broken-symmetry ωB97X-D/cc-pVDZ potential energy surface. In most cases, the diradical geometries distort to enhance through-space or through-bond coupling in the singlet states and to avoid Coulombic or exchange repulsion in the triplet states. EOM-SF-CCSD results indicate that all diradical isomers are two-configurational, closed shell singlet states. The only exceptions to this are for (2,3) and (2,4) thiophene and (2,3) pyrrole, which each contain more than two configurations. In all cases, the leading term in the multiconfigurational diradical wave function doubly occupies the symmetric radical σ orbital, indicative of either through-space or 1,3 through-bond coupling. We utilized the nucleus-independent chemical shift (NICS) approach to qualitatively assess aromaticity and find that this property varies and may be related to the energetic splittings in these diradical isomers. 
    more » « less
  3. 31 P NMR spectroscopy and the study of nuclear spin singlet relaxation phenomena are of interest in particular due to the importance of phosphorus-containing compounds in physiology. We report the generation and measurement of relaxation of 31 P singlet order in a chemically equivalent but magnetically inequivalent case. Nuclear magnetic resonance singlet state lifetimes of 31 P pairs have heretofore not been reported. Couplings between 1 H and 31 P nuclei lead to magnetic inequivalence and serve as a mechanism of singlet state population conversion within this molecule. We show that in this molecule singlet relaxation occurs at a rate significantly faster than spin–lattice relaxation, and that anticorrelated chemical shift anisotropy can account for this observation. Calculations of this mechanism, with the help of molecular dynamics simulations and ab initio calculations, provide excellent agreement with the experimental findings. This study could provide guidance for the study of 31 P singlets within other compounds, including biomolecules. 
    more » « less
  4. Abstract Large separation of magnetic levels and slow relaxation in metal complexes are desirable properties of single‐molecule magnets (SMMs). Spin‐phonon coupling (interactions of magnetic levels with phonons) is ubiquitous, leading to magnetic relaxation and loss of memory in SMMs and quantum coherence in qubits. Direct observation of magnetic transitions and spin‐phonon coupling in molecules is challenging. We have found that far‐IR magnetic spectra (FIRMS) of Co(PPh3)2X2(Co‐X; X=Cl, Br, I) reveal rarely observed spin‐phonon coupling as avoided crossings between magnetic andu‐symmetry phonon transitions. Inelastic neutron scattering (INS) gives phonon spectra. Calculations using VASP and phonopy programs gave phonon symmetries and movies. Magnetic transitions among zero‐field split (ZFS) levels of theS=3/2 electronic ground state were probed by INS, high‐frequency and ‐field EPR (HFEPR), FIRMS, and frequency‐domain FT terahertz EPR (FD‐FT THz‐EPR), giving magnetic excitation spectra and determining ZFS parameters (D, E) andgvalues. Ligand‐field theory (LFT) was used to analyze earlier electronic absorption spectra and give calculated ZFS parameters matching those from the experiments. DFT calculations also gave spin densities inCo‐X, showing that the larger Co(II) spin density in a molecule, the larger its ZFS magnitude. The current work reveals dynamics of magnetic and phonon excitations in SMMs. Studies of such couplings in the future would help to understand how spin‐phonon coupling may lead to magnetic relaxation and develop guidance to control such coupling. 
    more » « less
  5. The examination and optimized preparation of nuclear spin singlet order has enabled the development of new types of applications that rely on potentially long-term polarization storage. Lifetimes several orders of magnitude longer than T 1 have been observed. The efficient creation of such states relies on special pulse sequences. The extreme cases of very large and very small magnetic equivalence received primary attention, while relatively little effort has been directed towards studying singlet relaxation in the intermediate regime. The intermediate case is of interest as it is relevant for many spin systems, and would also apply to heteronuclear systems in very low magnetic fields. Experimental evidence for singlet–triplet leakage in the intermediate regime is sparse. Here we describe a pulse sequence for efficiently creating singlets in the intermediate regime in a broad-band fashion. Singlet lifetimes are studied with a specially synthesized molecule over a wide range of magnetic fields using a home-built sample-lift apparatus. The experimental results are supplemented with spin simulations using parameters obtained from ab initio calculations. This work indicates that the chemical shift anisotropy (CSA) mechanism is relatively weak compared to singlet–triplet leakage for the proton system observed over a large magnetic field range. These experiments provide a mechanism for expanding the scope of singlet NMR to a larger class of molecules, and provide new insights into singlet lifetime limiting factors. 
    more » « less