Precise Control of Local Spin States in an Adsorbed Magnetic Molecule with an STM Tip: Theoretical Insights from First-Principles-Based Simulation

Differences in Kondo Splitting of Surface Quantum Systems Induced by Two Distinct Magnetic Tips: A Joint Method of DFT and HEOM

The interactions between a magnetic tip and local spin impurities initiate unconventional Kondo phenomena, such as asymmetric suppression or even splitting of the Kondo peak. However, a lack of realistic theoretical models and comprehensive explanations for this phenomenon persists due to the complexity of the interactions. This research employs a joint method of density functional theory (DFT) and hierarchical equation of motion (HEOM) to simulate and contrast the modulation of the spin state and Kondo behavior in the Fe/Cu(100) system with two distinct magnetic tips. A cobalt tip, possessing a larger magnetic moment, incites greater atomic displacement of the iron atom, more notable alterations in electronic structure, and enhanced charge transfer with the environment compared with the control process utilizing a nickel tip. Furthermore, the Kondo resonance undergoes asymmetric splitting as a result of the ferromagnetic correlation between the iron atom and the magnetic tip. The Co tip's higher spin polarization results in a wider spacing between the splitting peaks. This investigation underscores the precision of the DFT + HEOM approach in predicting complex quantum phenomena and explaining the underlying physical principles. This provides valuable theoretical support for developing more sophisticated quantum regulation experiments.

Unraveling the Nature of Spin Coupling in a Metal-Free Diradical: Theoretical Distinction of Ferromagnetic and Antiferromagnetic Interactions

Metal-free diradicals based on polycyclic aromatic hydrocarbons are promising candidates for organic spintronics due to their stable magnetism and tunable spin coupling. However, distinguishing and elucidating the origins of ferromagnetic and antiferromagnetic interactions in these systems remain challenging. Here, we investigate the 2-OS diradical molecule sandwiched between gold electrodes using a combined density functional theory and hierarchical equations of motion approach. We find that the dihedral angle between the radical moieties controls the nature and strength of the intramolecular spin coupling, transitioning smoothly from antiferromagnetic to ferromagnetic as the angle increases. Distinct features in the inelastic electron tunneling spectra are identified that can discern the two coupling regimes, including spin excitation steps whose energies directly reveal the exchange coupling constant. Mechanical stretching of the junction is predicted to modulate the spectral line shapes by adjusting the hybridization of the molecular radicals with the electrodes. Our work elucidates the electronic origin of tunable spin interactions in 2-OS and provides spectroscopic fingerprints for characterizing magnetism in metal-free diradicals.

Unconventional Surface Doping Effect on the Spin State of an Adsorbed Magnetic Molecule

Magnetic molecules adsorbed on two-dimensional (2D) substrates have attracted broad attention because of their potential applications in quantum device applications. Experimental observations have demonstrated substantial alteration in the spin excitation energy of iron phthalocyanine (FePc) molecules when adsorbed on nitrogen-doped graphene substrates. However, the underlying mechanism responsible for this notable change remains unclear. To shed light on this, we employ an embedding method and ab initio quantum chemistry calculations to investigate the effects of surface doping on molecular properties. Our study unveils an unconventional chemical bonding at the interface between the FePc molecule and the N-doped graphene. This bonding interaction, stronger than non-covalent interactions, significantly modifies the magnetic anisotropy energy of the adsorbed molecule, consistent with experimental observations. These findings provide valuable insights into the electronic and magnetic properties of molecules on 2D substrates, offering a promising pathway for precise manipulation of molecular spin states.

Voltage-induced modulation of the magnetic exchange in binuclear Fe(III) complex deposited on Au(111) surface

We present a complete computational study devoted to the deposition of a magnetic binuclear complex on a metallic surface, aimed to obtain insight into the interaction of magnetically coupled complexes with their supporting substrates, as well as their response to external electrical stimuli applied through a surface-molecule-STM molecular junction-like architecture. Our results not only show that the deposition is favorable in two of the four studied orientations, but also, that the magnetic coupling is only slightly perturbed once the complex is adsorbed. We observe that the effects of the applied bias voltage on the magnetic coupling strongly depend on the molecule orientation with respect to the surface and the voltage polarity. Further analysis shows that this behavior is attributable to the stabilization/destabilization of the d-type singly occupied orbitals of the iron centers, reinforced by the strong local electric fields and induced charge densities only present in certain orientations of the deposited molecule and applied voltage polarity.

Tweezer-like magnetic tip control of the local spin state in the FeOEP/Pb(111) adsorption system: a preliminary exploration based on first-principles calculations

The magnetic interactions between the spin-polarized scanning tunnelling microscopy (SP-STM) tip and the localized spin impurities lead to various forms of the Kondo effect. Although these intriguing phenomena enrich Kondo physics, detailed theoretical simulations and explanations are still lacking due to the rather complex formation mechanisms. Here, by combining density functional theory (DFT), complete active space self-consistent field (CASSCF) theory, and hierarchical equations of motion (HEOM) methods, we perform first-principles-based simulation to elaborate the regulation process of the magnetic Co-tip on the spin state and transport behaviour of FeOEP/Pb(111) system. Compared with the non-magnetic tip, the stronger interaction between the magnetic tip and FeOEP molecule results in a more drastic deformation of the molecular structure with more electron transfer from the local environment to Fe-3d orbitals. The magnetic anisotropy of FeOEP changes very drastically from positive values in the tunnelling region to negative values in the contact region. The ferromagnetic electron correlation between the magnetic tip and the molecule induces an asymmetric Kondo line-shape near the Fermi level. This work highlights that the DFT + CASSCF + HEOM approach can not only predict complex quantum phenomena and explain underlying physical mechanisms, but also facilitate the design of more fascinating quantum control experiments.

Hierarchical equations of motion approach for accurate characterization of spin excitations in quantum impurity systems

Recent technological advancement in scanning tunneling microscopes has enabled the measurement of spin-field and spin-spin interactions in single atomic or molecular junctions with an unprecedentedly high resolution. Theoretically, although the fermionic hierarchical equations of motion (HEOM) method has been widely applied to investigate the strongly correlated Kondo states in these junctions, the existence of low-energy spin excitations presents new challenges to numerical simulations. These include the quest for a more accurate and efficient decomposition for the non-Markovian memory of low-temperature environments and a more careful handling of errors caused by the truncation of the hierarchy. In this work, we propose several new algorithms, which significantly enhance the performance of the HEOM method, as exemplified by the calculations on systems involving various types of low-energy spin excitations. Being able to characterize both the Kondo effect and spin excitation accurately, the HEOM method offers a sophisticated and versatile theoretical tool, which is valuable for the understanding and even prediction of the fascinating quantum phenomena explored in cutting-edge experiments.

Unveiling the Decisive Factor for the Sharp Transition in the Scanning Tunneling Spectroscopy of a Single Nickelocene Molecule

Scanning tunneling microscopy (STM) has been utilized to realize the precise measurement and control of local spin states. Experiments have demonstrated that when a nickelocene (Nc) molecule is attached to the apex of an STM tip, the dI/dV spectra exhibit a sharp or a smooth transition when the tip is displaced toward the substrate. However, what leads to the two distinct types of transitions remains unclear, and more intriguingly, the physical origin of the abrupt change in the line shape of dI/dV spectra remains unclear. To clarify these intriguing issues, we perform first-principles-based simulations on the STM tip control process for the Cu tip/Nc/Cu(100) junction. In particular, we find that the suddenly enhanced hybridization between the d orbitals on the Ni ion and the metallic bands in the substrate leads to Kondo correlation overwhelming spin excitation, which is the main cause of the sharp transition in the dI/dV spectra observed experimentally.

Controlling the Spin States of FeTBrPP on Au(111)

Spin-flip excitations of iron porphyrin molecules on Au(111) are investigated with a low-temperature scanning tunneling microscope. The molecules adopt two distinct adsorption configurations on the surface that exhibit different magnetic anisotropy energies. Density functional theory calculations show that the different structures and excitation energies reflect unlike occupations of the Fe 3d levels. We demonstrate that the magnetic anisotropy energy can be controlled by changing the adsorption site, the orientation, or the tip-molecule distance.

Unravelling the robustness of magnetic anisotropy of a nickelocene molecule in different environments: a first-principles-based study

Recent scanning tunneling spectroscopy with single metallocene molecule-functionalized tips have proved to be a powerful tool to probe and control individual spins and spin-spin exchange interactions due to the robustness of the magnetic properties of the metallocene molecule in different surroundings. However, accurate prediction of such robustness at a first-principles-based level by the conventional density functional theory (DFT) has remained challenging. In this paper, we have performed a detailed investigation of the evolution of electronic and magnetic properties of a nickelocene molecule (NiCp₂) in different environments, i.e., free-standing, adsorbed on Cu(100) and as a functionalized tip apex. Using an embedding method, which combines DFT and the complete active space self-consistent field (CASSCF) method recently developed, we demonstrate that the nickelocene molecule almost preserves its spin and magnetic anisotropy upon adsorption on Cu(100), and also in the position of the tip apex. In particular, the cyclic π* orbital of the Cp rings could hybridize with the singly occupied d_π orbitals of the Ni center of the molecule, protecting these orbitals from external states. Hence the molecular spin maintains S = 1, the same as in the free-standing case, and its magnetic anisotropy is also robust with energies of 3.56, 3.34, and 3.51 meV in free-standing, adsorbed on Cu(100), and functionalized tip apex states, respectively, in good agreement with previous theoretical and experimental results. This work thus provides a first-principles-based understanding of the relevant experiments. Such agreement between theoretical simulations and experimental measurements highlights the potential usefulness of the method for investigating the local electronic and spin states of organometallic molecule-surface composite systems.

Origin of Asymmetric Splitting of Kondo Peak in Spin-Polarized Scanning Tunneling Spectroscopy: Insights from First-Principles-Based Simulations

The spin-polarized scanning tunneling microscope (SP-STM) has served as a versatile tool for probing and manipulating the spintronic properties of atomic and molecular devices with high precision. The interplay between the local spin state and its surrounding magnetic environment significantly affects the transport behavior of the device. Particularly, in the contact regime, the strong hybridization between the SP-STM tip and the magnetic atom or molecule could give rise to unconventional Kondo resonance signatures in the differential conductance (dI/dV) spectra. This poses challenges for the simulation of a realistic tip control process. By combining the density functional theory and the hierarchical equations of motion methods, we achieve first-principles-based simulation of the control of a Ni-tip/Co/Cu(100) junction in both the tunneling and contact regimes. The calculated dI/dV spectra reproduce faithfully the experimental data. A cotunneling mechanism is proposed to elucidate the physical origin of the observed unconventional Kondo signatures.

Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.

Precise Spin Manipulation of Single Molecule Positioning on Graphene by Coordination Chemistry

Precise spin manipulation of single molecules is crucial for future molecular spintronics. However, it has been a formidable challenge due to the complexities of the strong molecule-substrate coupling as well as the response of the molecule to external stimulus. Here we demonstrate by density functional theory calculations that precise spin manipulation can be achieved by extra CO and NO molecules coordination to transition metal phthalocyanine (TMPc) (TM = Co, Fe, Mn) molecules deposited on metal-supported graphene; the spins of TMPc molecules are switched from S to S - ¹/₂ (|S - 1|) after NO (CO) coordination. With the aid of a combination of molecular orbitals (MO) theory and recently developed principal interacting spin-orbital (PISO) analysis, the impacts of NO and CO coordinations on both adsorption configuration and spin polarization of TMPc are well elucidated. We reveal the different coordination geometries that CO always coordinates axially to the TM center with a linear geometry, while NO prefers a bent geometry, which can be attributed to the competition between the σ- and π-type interactions according to the PISO analysis. Particularly, the NO-MnPc complex adopts a bent geometry deviating from the prediction by the existing Enemark-Feltham formalism. In addition, MO analysis suggests that during the CO coordination, the simultaneous existence of σ-donation and π-back-donation promotes electrons flowing from the d_z² to partially occupied d_π (d_xz and d_xz) orbitals with subsequent reordering of the TM d-orbitals, resulting in the spin transition of S → |S - 1|. In comparison, given that NO is regarded as NO^- when it adopts a bent geometry coordinating to the TM center, the complete (CoPc) or partial (FePc and MnPc) quenching of the molecular spins caused by NO coordination is attributed to the electron transfer from TM to NO. These theoretical findings provide important insights into relevant experiments and offer an effective design strategy to realize underlying single-molecular spintronics devices integrated with two-dimensional materials.

Probing Magnetism in Artificial Metal-Organic Complexes Using Electronic Spin Relaxometry

Single spins are considered as a versatile candidate for miniaturizing information devices down to the nanoscale. To engineer the spin's properties, metal-organic frameworks provide a promising route which in turn requires thorough understanding of the metal-molecule interaction. Here, we investigate the magnetic robustness of a single iron (Fe) atom in artificially built Fe-tetracyanoethylene (TCNE) complexes by using low-temperature scanning tunneling microscopy (STM). We find that the magnetic anisotropy and spin relaxation dynamics of the Fe atom within the complexes remain unperturbed in comparison to well-isolated Fe atoms. Density functional theory (DFT) calculations support our experimental findings, suggesting that the 3d orbitals of the Fe atom remain largely undisturbed while the 4s and 4p orbitals are rearranged in the process of forming a complex. To precisely determine the location of the spin center within the complex, we utilize STM-based spin relaxometry, mapping out the spatial dependence of spin relaxation with subnanometer resolution. Our work suggests that the magnetic properties of atoms can remain unchanged while being embedded in a weakly bound molecular framework.

Hierarchical equations of motion method based on Fano spectrum decomposition for low temperature environments

The hierarchical equations of motion (HEOM) method has become one of the most popular methods for the studies of the open quantum system. However, its applicability to systems at ultra-low temperatures is largely restrained by the enormous computational cost, which is caused by the numerous exponential functions required to accurately characterize the non-Markovian memory of the reservoir environment. To overcome this problem, a Fano spectrum decomposition (FSD) scheme has been proposed recently [Cui et al., J. Chem. Phys. 151, 024110 (2019)], which expands the reservoir correlation functions using polynomial-exponential functions and hence greatly reduces the size of the memory basis set. In this work, we explicitly establish the FSD-based HEOM formalisms for both bosonic and fermionic environments. The accuracy and efficiency of the FSD-based HEOM are exemplified by the calculated low-temperature dissipative dynamics of a spin-boson model and the dynamic and static properties of a single-orbital Anderson impurity model in the Kondo regime. The encouraging numerical results highlight the practicality and usefulness of the FSD-based HEOM method for general open systems at ultra-low temperatures.

Electronic and magnetic properties of CoPc and FePc molecules on graphene: the substrate, defect, and hydrogen adsorption effects

Transition metal phthalocyanines (TMPcs) are particularly appealing for spintronic processing and data storage devices due to their structural simplicity and functional flexibility. To realize effective control of the spins in TMPc-based systems, it is necessary to quantify how the structural and chemical environment of the molecule affects its spin center. Herein we perform a detailed investigation of the electronic and spintronic properties of vertically stacked heterostructures formed by CoPc or FePc adsorbed on a monolayer of graphene under the influences of the gold substrate, vacancies in graphene, and extra atomic hydrogen coordination on the TMPc. By using density functional theory (DFT), we reveal that both the TMPc molecules prefer the carbon-top position on graphene, and the existence of the Au substrate enhances the stability of the adsorption, while this enhanced adsorption will not modify the molecular magnetism, keeping it the same value as in the free standing case. Moreover, with the aid of a combination of DFT and ab initio wavefunction-based calculations, our results indicate that the magnetic anisotropy of the FePc-graphene complex can be actively tuned by the Au substrate. Our calculations also show that defects in graphene including single and double vacancies can modify the magnetism of these heterostructures. In particular, the spin state of FePc can be tuned from S = 1 to S = 2 with such defect engineering. Further spin state tunability can be achieved from a hydrogenation process, with the coordination of one extra hydrogen on the Co-top site for CoPc and the pyridinic N site for FePc, respectively, tuning their spin states from S = 1/2 to S = 0 and from S = 1 to S = 2. These findings may prove to be instrumental for rational design of future molecular spintronics devices integrated with two-dimensional materials.

Manipulation of spin and magnetic anisotropy in bilayer magnetic molecular junctions

The Kondo effect and magnetic anisotropy in bilayer TMPc/TMPc/Pb(111) junctions can be actively tuned by changing the intermediate decoupling layer.