Spin-flip excitations induced by dehydrogenation in a magnetic single-molecule junction

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.

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.

Local decomposition of hybridization functions: Chemical insight into correlated molecular adsorbates

Hybridization functions are an established tool for investigating the coupling between a correlated subsystem (often a single transition metal atom) and its uncorrelated environment (the substrate and any ligands present). The hybridization function can provide valuable insight into why and how strong correlation features such as the Kondo effect can be chemically controlled in certain molecular adsorbates. To deepen this insight, we introduce a local decomposition of the hybridization function, based on a truncated cluster approach, enabling us to study individual effects on this function coming from specific parts of the systems (e.g., the surface, ligands, or parts of larger ligands). It is shown that a truncated-cluster approach can reproduce the Co 3d and Mn 3d hybridization functions from periodic boundary conditions in Co(CO)₄/Cu(001) and MnPc/Ag(001) qualitatively well. By locally decomposing the hybridization functions, it is demonstrated at which energies the transition metal atoms are mainly hybridized with the substrate or with the ligand. For the Kondo-active 3d_x²-y² orbital in Co(CO)₄/Cu(001), the hybridization function at the Fermi energy is substrate-dominated, so we can assign its enhancement compared with ligand-free Co to an indirect effect of ligand-substrate interactions. In MnPc/Ag(001), the same is true for the Kondo-active orbital, but for two other orbitals, there are both direct and indirect effects of the ligand, together resulting in such strong screening that their potential Kondo activity is suppressed. A local decomposition of hybridization functions could also be useful in other areas, such as analyzing the electrode self-energies in molecular junctions.

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.