Stability, energetics, and magnetic states of cobalt adatoms on graphene

Same effect of biquadratic exchange interaction and Heisenberg linear interaction in a spin spiral

Monolayer (ML) NiCl₂ exhibits a strong biquadratic exchange interaction between the first neighboring magnetic atoms (B₁), as demonstrated by the spin spiral model in J. Ni et al., Phys. Rev. Lett., 2021, 127, 247204. This interaction is crucial for stabilizing the ferromagnetic collinear order within the ML NiCl₂. However, they neither point out the role of B₁ nor discuss the dispersion relation from spin orbit coupling (SOC) in the spin spiral. As we have done in this work, these parameters might theoretically potentially be derived directly by fitting the calculated spin spiral dispersion relation. Here, we draw attention to the fact that B₁ is equivalent to half of J₃ in Heisenberg linear interactions and that the positive B₁ partially counteracts the negative J₃'s impact on the spin spiral to make the ML NiCl₂ ferromagnetic. The comparatively small J₃ + 1/2B₁ from the spin spiral led us to believe that J₃ could be substituted by B₁, yet it still exists and plays a crucial function in magnetic semiconductors or insulators. The dispersion relation, which we also obtain from SOC, displays weak antiferromagnetic behavior in the spin spiral.

Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry

The accurate prediction of reaction mechanisms in heterogeneous (surface) catalysis is one of the central challenges in computational chemistry. Quantum Monte Carlo methods─Diffusion Monte Carlo (DMC) in particular─are being recognized as higher-accuracy, albeit more computationally expensive, alternatives to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the broader adoption of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that it is possible to significantly reduce this computational cost by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like adsorption (binding) energies and mapping potential energy surfaces. We analyze the cancellation and convergence of many-body finite-size errors in two well-known adsorbate/slab systems, H₂O/LiH(001) and CO/Pt(111). Based on this analysis, we identify strategies for obtaining binding energies in the thermodynamic limit that optimally utilize error cancellation to balance accuracy and computational efficiency. Using one such strategy, we then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for a wide range of density functionals. Our accurate and inexpensive DMC calculations are found to unambiguously recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to use this DMC method to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.

Visualization of Transition Metal Decoration on h-BN Surface

Functional h-BN (hexagonal boron nitride) has been prepared via the incorporation of transition metal (TM) impurities like nanoparticles and single atoms. Herein, scanning transmission electron microscopy (STEM) combined with density functional theory (DFT) was employed to study Ta-, Co-, Ni-, and Ir-decorated h-BN monolayers to provide an overview of their preferential site occupancies and morphological evolutions on h-BN. Ta, Ni, Ir, and Co single atoms are all positioned on the nitrogen of h-BN; however DFT predicts the occupancy site can vary with their spin state. In terms of microstructural evolution, Co, Ni, and Ir atoms form 3D nanoclusters while Ta atoms are well dispersed and thus the single Ta atom can be decorated on h-BN. This study highlights on TM/h-BN interaction dynamics and presents an avenue for designing nanostructures for electrocatalytic application.

Current-Induced One-Dimensional Diffusion of Co Adatoms on Graphene Nanoribbons

One-dimensional diffusion of Co adatoms on graphene nanoribbons has been induced and investigated by means of scanning tunnelling microscopy (STM). To this end, the nanoribbons and the Co adatoms have been imaged before and after injecting current pulses into the nanoribbons, with the STM tip in direct contact with the ribbon. We observe current-induced motion of the Co atoms along the nanoribbons, which is approximately described by a distribution expected for a thermally activated one-dimensional random walk. This indicates that the nanoribbons reach temperatures far beyond 100 K, which is well above the temperature of the underlying Au substrate. This model system can be developed further for the study of electromigration at the single-atom level.

Regional Embedding Enables High-Level Quantum Chemistry for Surface Science

Compared to common density functionals, ab initio wave function methods can provide greater reliability and accuracy, which could prove useful when modeling adsorbates or defects of otherwise periodic systems. However, the breaking of translational symmetry necessitates large supercells that are often prohibitive for correlated wave function methods. As an alternative, this paper introduces the regional embedding approach, which enables correlated wave function treatments of only a target fragment of interest through small, fragment-localized orbital spaces constructed using a simple overlap criterion. Applications to the adsorption of water on lithium hydride, hexagonal boron nitride, and graphene substrates show that regional embedding combined with focal-point corrections can provide converged CCSD(T) (coupled-cluster) adsorption energies with very small fragment sizes.

The effect of Hubbard-like interaction on molecular magnetism of TM-coronene complex (TM = Fe and Co)

The adsorption of two transition metal adatoms, Fe and Co added to a coronene molecule were studied by means of a full potential local orbital method in the framework of relativistic density functional theory. A sequence of fixed spin moment calculations based on the Hubbard-like interaction (U) were carried out on the basis of the generalized gradient approximation (GGA) of the exchange-correlation functional. We found a transition from low-spin to high-spin state upon increasing the transition metal-coronene distance. Furthermore, the magnetism of both Fe-coronene and Co-coronene complexes revealed sensitivity to the magnitude of the U values. In our GGA+U (U  =  2 eV) calculations, a low-spin ground state was found with in-plane magnetic anisotropy for both Fe-coronene and Co-coronene, whereas the GGA+U (U  =  4 eV) calculations resulted in high-spin state with out-of-plane and in-plane magnetic anisotropy for Fe-coronene and Co-coronene, respectively.

Strong current polarization and perfect negative differential resistance in few-FeN4-embedded zigzag graphene nanoribbons

The FeN₄-embedded graphene nanoribbons possess intrinsic current polarization and negative differential resistance behaviours, the performance of which can be enhanced easily.

3d Transition Metal Adsorption Induced the valley-polarized Anomalous Hall Effect in Germanene

Based on DFT + U and Berry curvature calculations, we study the electronic structures and topological properties of 3d transition metal (TM) atom (from Ti to Co) adsorbed germanene (TM-germanene). We find that valley-polarized anomalous Hall effect (VAHE) can be realized in germanene by adsorbing Cr, Mn, or Co atoms on its surface. A finite valley Hall voltage can be easily detected in their nanoribbon, which is important for valleytronics devices. Moreover, different valley-polarized current and even reversible valley Hall voltage can be archived by shifting the Fermi energy of the systems. Such versatile features of the systems show potential in next generation electronics devices.

High-coverage stable structures of 3d transition metal intercalated bilayer graphene

Alkali-metal intercalated graphite and graphene have been intensively studied for decades, where alkali-metal atoms are found to form ordered structures at the hollow sites of hexagonal carbon rings. Using first-principles calculations, we have predicted various stable structures of high-coverage 3d transition metal (TM) intercalated bilayer graphene (BLG) stabilized by the strain. Specifically, with reference to the bulk metal, Sc and Ti can form stable TM-intercalated BLG without strain, while the stabilization of Fe, Co, and Ni intercalated BLG requires the biaxial strain of over 7%. Under the biaxial strain ranging from 0% to 10%, there are four ordered sandwich structures for Sc with the coverage of 0.25, 0.571, 0.684, and 0.75, in which the Sc atoms are all distributed homogenously instead of locating at the hollow sites. According to the phase diagram, a homogenous configuration of C8Ti3C8 with the coverage of 0.75 and another inhomogeneous structure with the coverage of 0.692 were found. The electronic and magnetic properties as a function of strain were also analyzed to indicate that the strain was important for the stabilities of the high-coverage TM-intercalated BLG.

Growth mode and structures of magnetic Mn clusters on graphene

Lowest-energy geometries of Mn clusters on graphene. Blue and pink balls represent Mn atoms with negative and positive magnetic moments.