Next Article in Journal
Modelling the Effect of Zero-Field Splitting on the 1H, 13C and 29Si Chemical Shifts of Lanthanide and Actinide Compounds
Next Article in Special Issue
Five 2,6-Di(pyrazol-1-yl)pyridine-4-carboxylate Esters, and the Spin States of their Iron(II) Complexes
Previous Article in Journal
Improved 19F{1H} Saturation Transfer Difference Experiments for Sensitive Detection to Fluorinated Compound Bound to Proteins
Previous Article in Special Issue
Abrupt Spin Crossover Behavior in a Linear N1,N2-Triazole Bridged Trinuclear Fe(II) Complex
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interfacial Spin Manipulation of Nickel-Quinonoid Complex Adsorbed on Co(001) Substrate

by
Indukuru Ramesh Reddy
1,
Peter M. Oppeneer
2 and
Kartick Tarafder
1,*
1
Department of Physics, National Institute of Technology Karnataka, Srinivasnagar, Surathkal, Mangalore 575025, India
2
Department of Physics and Astronomy, Uppsala University, Box 516, S-75120 Uppsala, Sweden
*
Author to whom correspondence should be addressed.
Magnetochemistry 2019, 5(1), 2; https://doi.org/10.3390/magnetochemistry5010002
Submission received: 29 November 2018 / Revised: 20 December 2018 / Accepted: 21 December 2018 / Published: 24 December 2018
(This article belongs to the Special Issue Spin-Crossover Beyond the Immediate Tribute)

Abstract

:
We studied the structural, electronic, and magnetic properties of a recently synthesized Ni(II)-quinonoid complex upon adsorption on a magnetic Co(001) substrate. Our density functional theory + U (DFT+U) calculations predict that the molecule undergoes a spin-state switching from low-spin S = 0 in the gas phase to high-spin S 1 when adsorbed on the Co(001) surface. A strong covalent interaction of the quinonoid rings and surface atoms leads to an increase of the Ni–O(N) bond lengths in the chemisorbed molecule that support the spin-state switching. Our DFT+U calculations show that the molecule is ferromagnetically coupled to the substrate. The Co surface–Ni center exchange mechanism was carefully investigated. We identified an indirect exchange interaction via the quinonoid ligands that stabilizes the molecule’s spin moment in ferromagnetic alignment with the Co surface magnetization.

Graphical Abstract

1. Introduction

Molecular spintronics has gained enormous attention in the past two decades because it holds promise for advanced multifunctional spin devices [1,2]. Especially spin-bearing molecules with one or more transition metal (TM) ions having a partially filled d shell with d 4 to d 7 configurations have drawn special attention as for many such complexes the spin properties can be manipulated externally. For spin-crossover complexes, it has been shown that it is possible to switch reversibly between different spin states via an external perturbation such as pressure, temperature, light irradiation, electric or magnetic fields [3,4]. Many magnetic molecules have been synthesized and intensively investigated both theoretically and experimentally in the past years [5,6,7,8]. The progress made recently in the application of spin-crossover complexes in functional devices has been reviewed lately [9,10]. The understanding of the underlying electronic and magnetic interactions on the single molecule level is, however, still sparse. Efficient and deterministic spin detection, manipulation and transportation in organic materials are yet to be achieved, properties that form an integral part of any advanced spin device.
Growing molecules on magnetic as well as inert surfaces with atomic precision is one of the approaches to develop such devices. Immobilized magnetic metal–organic molecules with tailorable magnetic properties can be deposited on solid substrates to construct building blocks for molecular spintronic devices [11,12]. In this process, the molecules exhibit magnetic stability at higher temperatures due to the substrate induced magnetic exchange coupling and/or substrate induced magnetic anisotropy [13,14,15]. Moreover, new magnetochemical properties arise due to the molecule–substrate interaction that can also offer an advantage over traditional spintronics [16,17,18,19,20]. It has been observed that magnetic molecules reacting on a metallic surface may lose their magnetic properties by coupling with the electrodes, even if they remain intact [21]. Hence, a proper understanding of molecule–substrate interactions is essential and needs to be achieved well before its technological application in data storage, sensor, quantum computing, or other single molecular spintronic device applications [22].
An isolated molecule that is at a relatively large distance from the surface/electrode or in the gas phase has discrete molecular energy levels. The bonding interactions may occur at the interface when molecules are brought near to the surface and adsorbed on it. The interaction can vary from weak (in physisorption) to strong bonding interactions (chemisorption) depending on the physical and chemical nature of the molecule and surface. The electronic structure of adsorbed molecules may change considerably. In physisorption, the molecule–substrate separation is relatively large, interaction is relatively weak, and the discrete molecular electronic levels are weakly broadened. On the other hand, electronic levels of chemisorbed adsorbates are significantly changed. Atoms of the adsorbed molecule in proximity with surface atoms form direct chemical bonds with surface atoms, which results in a much stronger substrate–molecule interaction. Particularly, in the case of magnetic substrate, new spin-resolved interface electronic states may appear near the Fermi energy due to the π d hybridization [23,24] that may completely change the electronic structure of the molecule and also affect the molecular geometry. To study and understand the chemical and magnetic exchange interactions between the adsorbed metal–organic molecules and magnetic substrates is therefore a fundamental focus. Such study may reveal an effective spin manipulation mechanism of the adsorbed molecule. The chemical environment between the adsorbed molecule and substrate determines the type and strength of magnetic exchange interactions which can be controlled externally by changing chemical environment [25,26]. Efficient chemical modification can lead to a significant change in π d hybridization across the interface and provides a potential path for molecular magnetism that is stable at room temperature ferromagnetism [13,23]. Hence, the interface magnetic properties that arises due to interface chemistry of molecule and surface, demand a careful study with the help of advanced computational and sophisticated experimental tools (XMCD and SP-STM) [27] to provide an insight into the spin-interface chemistry [28,29].
In general, non-planar molecules would bind weakly to a metal surface due to steric hindrance, whereas planar molecules could bind more strongly or self-assemble on the metal surfaces. Therefore, planar molecules are suitable candidates to design molecular devices by depositing them on matching metal surfaces. A number of theoretical as well as experimental observations have been reported on the magnetic interactions of transition-metal porphyrins (TM-P) [13,23], transition-metal phthalocyanine (TM-Pc) [12,30,31], and cyclohexane-based organic molecules [32] with metal substrates. It has been observed that planar TM-P and TM-Pc molecules are very promising as they show a stable magnetic behavior due to exchange coupling to the magnetic substrate. Moreover, their magnetic properties can be tailored by engineering the spin states of the metal center using on-surface magnetochemistry, i.e., chemical perturbation by coordination in the free ligand position or through interlayer modified exchange interaction [25,29,33,34,35,36]. Spin-state switching is also observed for a planar metal–organic molecule due to stretching the metal–ligand distance [37], similar to the bond-length induced switching of spin-crossover molecules [3]. While this focuses primarily on bond-length changes, low-spin to high-spin state switching in addition brings about thermodynamic changes such as entropy increase, which could play a role in the adsorption process.
Recently, Kar et al. designed a methanol-triggered vapochromism and spin-switching in a Ni(II)-quinonoid (NiQ) complex [Ni(HL) 2 , H 2 L = 4-methylamino-6-methyliminio-3- oxocyclohexa-1, 4-dien-1-olate], which is a planar molecule. They showed that the addition of methanol (MeOH) to the molecule leads to a magnetochemical reaction, in which the central Ni atom of the molecule changes it co-ordination from square-planar to octahedral by axial ligation of two MeOH groups at the Ni site, accompanied by a temperature-robust spin transition from S = 0 to S = 1 [38]. This observation suggests that the Ni(II)-quinonoid could be a potential candidate material for application in spintronics devices. However, for applications, a detailed understanding of the molecular adsorption mechanism on a metal surface is essential and, in addition, a major question is if magnetic interface interactions are sufficient to trigger the spin-switching.
The main objective of the present investigation was therefore to study the geometric and electronic structure of the multifunctional NiQ molecule when adsorbed on a Co(001) substrate, the emergent magnetic interactions at the spin-interface, and, specifically, whether an on-surface spin-state transition can be achieved. Thereto, we carried out a detailed theoretical investigation using density functional + U theory to unravel the nature of adsorption, electronic structure, and the magnetic state of the adsorbed NiQ molecule. Our investigation predicts that the NiQ molecule undergoes an on-surface spin-state transition to a high spin state when adsorbed on a Co(001) surface.

2. Computational Methodology

We employed density functional theory (DFT) calculations to study the magnetic interaction of a NiQ molecule with a ferromagnetic cobalt surface. The calculations were performed using the plane wave, pseudo-potential method, as implemented in Vienna Ab-initio Simulation Package (VASP) [39] with projector augmented plane wave (PAW) potential [40]. For the exchange-correlation function, the generalized gradient approximation (GGA) was used with the Perdew–Burke–Ernzerhof parameterization [41]. We further used the DFT+U technique to capture the strong electron–electron correlation effect that exists in the partially filled 3 d shell of Ni and the missing correlation effect beyond the plain GGA. This technique was proven to be very promising to achieve the precise spin state of the molecule [24,42,43,44]. The on-site Coulomb and exchange parameters U and J used in this study were chosen as 5.0 eV and 1.0 eV for Ni atom. These values of the U and J parameters have previously been found to provide an accurate description of the low and high spin states of bi-stable metalorganic molecules [44,45,46,47]. A 3 × 3 × 1 Monkhorst–Pack k-point grid was used for reciprocal space sampling. A sufficiently large plane-wave kinetic energy cutoff 500 eV was considered to obtain the desire accuracy. The convergence criterion was set to be 10 5 eV for the self-consistent electronic minimization. We modeled the Ni(II)-quinonoid molecule on top of three atomic Co layers within a large super-cell of the Co(001) magnetic substrate (we used 8 × 5 lateral supercell) and applied boundary conditions to maintain the periodicity of the surface. We optimized the geometry of the molecule on the Co surface by relaxing the atoms’ positions using the conjugate gradient method until all residual inter-atomic forces were minimized up to 0.01 eV/Å. Two bottom layers of Co atoms of the substrate were kept fixed during the optimization process to minimize the computational effort.

3. Results and Discussion

3.1. Optimized Geometry

We modeled the square-planar NiQ molecule structure according to the reported data [38]. To start with, we optimized the geometry of the isolated NiQ molecule in the gas phase. The optimized structure of NiQ is shown in Figure 1. The molecule remains planar with average Ni–O(N) bond lengths of the molecule being 1.868 Å, which is close to the experimental value (1.866 Å).
Next, we placed the optimized NiQ molecule on the Co(001) substrate in two different configurations, namely the HOLLOW and TOP adsorption sites. In the HOLLOW configuration, the molecule adsorbed on the Co(001) surface in such a way that the central Ni atom of the molecule is placed on top of a second layer Co atom (i.e., hollow site) of the Co(001) substrate (see Figure 2). The optimized distance between the surface and Ni atom is 1.778 Å. In this configuration, the Ni–O and Ni–N bond lengths have been stretched in such a way that the O and N atoms of the molecule come closer to surface Co atoms and form direct chemical bonds with the surface. The average Ni–O(N) bond length is increased by 0.204 Å compared to its gas phase structure.
In TOP configuration (see Figure 3), the molecule is placed on the Co (001) surface such that the Ni atom sits exactly on top of a surface Co atom. The optimized molecule substrate separation in this configuration is 2.4864 Å. On optimization, the ligands connected to the Ni atom of the molecule are impetus to stretch and orient themselves in such a way that the central O atoms of the molecule, which were initially at the hollow site, move toward the top site. Therefore, the bond length between Ni and O is increased approximately by 0.26 Å compared to the gas phase, whereas the Ni–N bond length does not change significantly. This is mainly because the nearest C atoms, which were already bonded with surface Co atoms, prevent the further stretching of Ni–N bonds. The average bond length is increased by 0.128 Å. Our calculations show that the HOLLOW configuration is energetically more favorable, whereas the molecule is structurally distorted more in the TOP configuration. The magnetic moment of the Ni atom, Ni–O(N) bond length in gas phase, and in different chemisorbed configurations are shown in Table 1.
Our calculations predict a zero moment ( S = 0 ) for the gas-phase NiQ molecule, consistent with experiments [38]. We note, with regard to the obtained S = 0 state, that we verified that the zero-moment property is robust against variations of the Coulomb U value between 2 and 8 eV. The molecular spin moment changes when the molecule is chemisorbed on the Co(001) surface, to 1.561 and 1.337 μ B for HOLLOW and TOP site, respectively. Thus, an on-surface spin-state switching to a high-spin state can be induced by adsorption of NiQ on the Co substrate. Compared to the reported magnetochemical reaction [38] in which two MeOH groups bond to Ni in the axial ligand positions, in the on-surface reaction, the molecule bonds one-sided to the top-layer surface atoms. The spin moment on the molecule is less than what would be expected for the high-spin S = 1 state. This is however not unexpected due to the extended hybridization of the quinonoid atoms to the surface atoms and the structural distortion.
The physisorption of this class of molecules is also possible on the Co surface [42]. Therefore, we looked for such physisorption of NiQ on Co surface. The distance between the NiQ molecule and Co(001) surface was successively changed in both HOLLOW and TOP configurations and the total energy of the systems in each situation calculated. The computed energy versus distance curve for the HOLLOW site absorption is shown in Figure 4. The single energy minima at 1.8 Å clearly indicates that only chemisorption of the molecule on Co(001) is possible. A weak physisorption of the molecule was observed at a distance 3.5 Å in TOP configuration, however the calculated energy barrier was very small and the magnetic interaction was found to be negligible in this configuration. The details of geometric structure and magnetic behavior in physisorption in TOP geometry are given in Appendix A.

3.2. Electronic Structure and Magnetic Properties

The distorted molecular geometry in chemisorption leads to a significant change in the electronic energy levels and their spin occupation, which further affects the magnetic properties of the molecule. In the square planar Ni(II)-quinonoid molecule with d 8 configuration, the four ligands (two oxygen and two nitrogen atoms) are bonded to the central Ni atom in the basal x y plane and provide a crystal-field effect on the Ni 3 d -orbitals. The Coulomb repulsion along with strong σ -type hybridization leads to a strong ligand field (LF) that splits the d-orbitals of Ni in degenerate t 2 g and e g levels. In the absence of an axial co-ordination, one of the e g orbitals, i.e., d z 2 is relatively lower in energy compare to the other, d x 2 y 2 . A stronger LF effect is imposed on the orbitals when the Ni–O and Ni–N bond lengths are shorter due to the strong hybridization and Coulomb repulsion. As a result, all d-electrons occupy the lower, LF-splitting dominated energy levels and yield a low-spin state ( S = 0 ) in the gas phase, where all d-orbitals are completely occupied, except d x 2 y 2 which is completely empty (see Figure 5a).
To investigate the charge transfer between molecule and substrate, we have calculated the Bader charges for the chemisorbed NiQ molecule in HOLLOW and TOP configurations. The molecule gains 2.267 and 1.959 electron charges from the Co surface in the HOLLOW and TOP configuration, respectively. The Bader charge on the Ni atom increases however only by 0.169 and 0.090, respectively.
In the case of molecular chemisorption, electrons in the outer-plane π -orbitals of the quinonoid strongly interact with the surface state electrons and form chemical bonds between the molecule and surface atoms. As a result, the Ni–O(N) bond lengths are stretched, which leads to a significant reduction of the LF effect and subsequently reduces the d x 2 y 2 orbital energy. In this situation, the intra-atomic exchange interaction becomes important and the electrons occupy the Ni- 3 d orbitals according to Hund’s rule. Hence, the d x 2 y 2 orbital level, which was completely unoccupied before, is now partially occupied (see Figure 5b,c). This explains the change in the spin state of a chemisorbed NiQ molecule governed by the shift in the energy level of the d x 2 y 2 orbital, as shown in Figure 5.
The magnetic exchange interaction energies E e x in different chemisorption configurations of the molecule are calculated using E e x = E A F M E F M , where E A F M and E F M are the total energies of anti-parallel and parallel spin alignment between the central Ni atom of the molecule and substrate Co atoms, respectively. Our calculation shows that the molecule is ferromagnetically coupled with the substrate with coupling energy of 68 and 88 meV for HOLLOW and TOP configurations, respectively.
To understand the spin coupling mechanism, we carefully investigated the atom projected density of states (PDOS) of the systems. In the HOLLOW configuration, the Ni-ligating O and N atoms are located on top of the surface Co atoms and form chemical bonds with the substrate Co atoms underneath (Figure 7). We observed that there are overlapping peaks of Co, N, and O projected densities of states marked with vertical arrows in the Figure 6a, indicating strong hybridization between out-of-plane orbitals ( d z 2 , d x z , d y z , and p z ) of Co, O, and N atoms.The total charge density plot, shown in Figure 7 (left) as light green hypersurface, also confirms that the orbital overlapping is only between surface Co and molecular ligand’s orbitals as there is no direct charge density overlap between the central Ni atom of the molecule and the Co surface atoms. This clearly indicates that the molecule interacts with the substrate only through its ligands, which in turn strongly favors magnetic coupling through the indirect super exchange interaction.
We computed a total 1.561 μ B spin moment on the molecule in HOLLOW position, whereas 1.178 μ B spin moment is located on the central Ni atoms (Table 1). A small spin-polarization is also observed on the ligand atoms, which further indicates that the exchange coupling path between the magnetic center of the molecule and the spins of Co surface atoms goes through the molecular ligands. To obtain more insight in the nature and path of the magnetic exchange interaction between magnetic centers, we computed the partial magnetization density of the system, shown in Figure 8, at the energy region where the hybridization between atomic orbitals of Co surface Co and molecular orbitals from N, O, and Ni are dominant, which is predicted from the atom projected DOS (see the overlapping peaks in Figure 6a). It is evident from the magnetization density plot that the magnetic exchange interaction between Ni and surface Co atoms is mainly mediated through ligands of the molecule. The negative magnetic moments on O and N atoms, shown through the light blue isosurface in this figure, indicate a 90 ferromagnetic indirect exchange mechanism with the Kanamori–Goodenough superexchange rule. A very weak direct exchange interaction could also be possible when one considers some of the noticeable PDOS peaks of Ni and Co atoms that have overlap, as indicated with vertical bars in the atom projected DOS shown in Figure 6a. However, we did not find any direct charge or spin-density overlap between Ni and Co atoms in the charge density and magnetization plots.
In the case of chemisorption of the molecule in the TOP configuration, the direct O and N atoms bonding to the Ni atom are in the hollow position with respect to the surface atoms and do not directly bond to the surface Co atoms; however, there are other atoms of the quinonoid ligand that chemically bond to the surface. The total charge density plot in Figure 7 (right) clearly indicates the non-overlapping nature of the orbital wave functions between the central Ni atom and the surface. However, quinonoid ligand atoms of the molecule strongly interact with the surface Co atoms. There are overlapping peaks present in the atom-projected PDOS plot (Figure 6b) marked with vertical arrows, which indicates the strong hybridization between out-of-plane orbitals ( d z 2 ) of surface Co with O and N p orbitals. The magnetization density calculated for the TOP configuration is shown in Figure 9. A small negative magnetization density is present on the ligand atoms, as shown by the light blue isosurface that indicates an indirect exchange interaction between the Ni spin center and the magnetic Co surface atoms, mediated through the ligand.

4. Conclusions

To conclude, we have demonstrated that the adsorption of the NiQ molecule on a Co(001) surface leads to a spin crossover from S = 0 in the gas phase to a nearly high-spin S = 1 state when chemisorbed on the Co surface. Chemisorption at a HOLLOW site is found to be energetically most favorable. The covalent interactions between the molecule’s ligands and surface lead to a stretching of the Ni–O and Ni–N bond lengths that lead to the spin-crossover of the non-magnetic NiQ into a spin bearing molecule. An indirect superexchange path through the molecular ligand stabilizes the Ni spin of the molecule parallel to the substrate magnetization direction. In view of the on-going efforts to realize bi-stable molecules with reversible spin state that are atomically precise anchored to substrates, a future research direction is to explore the possibility of achieving spin-switching on non-magnetic metallic substrates. Our investigation not only sheds light on the interactions of a NiQ molecule with a magnetic Co(001) substrate, but also suggests a possible pathway to explore planar metal-quinonoid molecules as spin-bearing interface molecule in future spintronic devices.

Author Contributions

K.T. conceived of the presented idea and designed the study. I.R.R. performed the computations. K.T. and P.M.O. supervised this work. All authors discussed the results and contributed to the final manuscript.

Funding

We gratefully acknowledge financial support from the Indo-Swedish Research Collaboration, funded through the Swedish Research Council (VR) and the Department of Science and Technology, India (DST-SERB, Project No. SB/FTP/PS-032/2014), and the Knut and Alice Wallenberg Foundation (Contract No. 2015.0060).

Acknowledgments

We acknowledge computer time received from the Swedish National Infrastructure for Computing (SNIC).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DFTDensity Functional Theory
DFT + U Density Functional Theory + U
XMCDX-ray magnetic circular dichroism
SP-STMSpin-polarized scanning tunneling microscopy
TMTransition metal
TM-Ptransition-metal porphyrin
TM-PcTransition-metal phthalocyanine
NiQNi(II)-quinonoid
MeOHMethanol
VASPVienna Ab-initio Simulation Package
PAWProjector augmented plane wave
GGAGeneralized gradient approximation
LFLigand field
DOSDensity of states
PDOSPartial density of states

Appendix A. Physisorption of NiQ on Co(001)

We studied theoretically the possible physisorption of the NiQ molecule on Co(001) in the TOP configuration. The total energy of the molecule in TOP configuration versus the molecule–substrate distance is shown in Figure A1. The shallow energy minima at 3.5 Å clearly indicates that physisorption of the molecule is possible. We note however that we have not included Van der Waals interactions (see, e.g., [48,49]) in the present calculation, but these would be required for an accurate description of the physisorption. The shallow minimum obtained at 3.5 Å is hence mainly driven by higher-order electrostatic interactions [42]. The optimized geometry of physisorbed NiQ on the Co surface is shown in Figure A1 (right). In the physisorption state, the NiQ molecule preserves mostly its gas-phase properties. The Ni–O and Ni–N bond lengths, the average Ni–O(N) bond length (1.8764 Å), and the magnetic moment on the Ni atom are very close to those of the gas-phase molecule. The computed charge and magnetization densities are shown in Figure A2. The charge density plot shows that the there is no orbital overlap between the molecule and the substrate, and the magnetization density plot exemplifies the non-magnetic nature of the physisorbed molecule on substrate.
Figure A1. (Left) DFT+U calculated total energy ( Δ E) of NiQ on Co as function of the Ni atom-substrate distance. The curve clearly shows a shallow energy minima at 3.5 Å which corresponds to physisorption of the NiQ molecule on Co(001) substrate in TOP configuration. (Right) top view of the physisorbed NiQ molecule on the surface of Co(001), in TOP configuration. Ni, O, N, C, H and Co atoms are represented with orange, red, blue, black, magenta and green colored balls. Second layer (hollow site) Co atoms are represented by grey colored balls.
Figure A1. (Left) DFT+U calculated total energy ( Δ E) of NiQ on Co as function of the Ni atom-substrate distance. The curve clearly shows a shallow energy minima at 3.5 Å which corresponds to physisorption of the NiQ molecule on Co(001) substrate in TOP configuration. (Right) top view of the physisorbed NiQ molecule on the surface of Co(001), in TOP configuration. Ni, O, N, C, H and Co atoms are represented with orange, red, blue, black, magenta and green colored balls. Second layer (hollow site) Co atoms are represented by grey colored balls.
Magnetochemistry 05 00002 g0a1
Figure A2. Calculated charge density (left) and magnetization density (right) of physisorbed NiQ molecule on Co(001) in the TOP configuration. The computed magnetization density shows that there is no magnetization density present on the molecule. The bright yellow hypersurface depicts the positive magnetization density.
Figure A2. Calculated charge density (left) and magnetization density (right) of physisorbed NiQ molecule on Co(001) in the TOP configuration. The computed magnetization density shows that there is no magnetization density present on the molecule. The bright yellow hypersurface depicts the positive magnetization density.
Magnetochemistry 05 00002 g0a2

References

  1. Bogani, L.; Wernsdorfer, W. Molecular spintronics using single-molecule magnets. Nat. Mater. 2008, 7, 179. [Google Scholar] [CrossRef] [PubMed]
  2. Moodera, J.S.; Koopmans, B.; Oppeneer, P.M. On the path toward organic spintronics. MRS Bull. 2014, 39, 578–581. [Google Scholar] [CrossRef] [Green Version]
  3. Gütlich, P.; Goodwin, H.A. Spin Crossover in Transition Metal Compounds I; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  4. Hao, H.; Zheng, X.; Song, L.; Wang, R.; Zeng, Z. Electrostatic spin crossover in a molecular junction of a single-molecule magnet Fe2. Phys. Rev. Lett. 2012, 108, 017202. [Google Scholar] [CrossRef] [PubMed]
  5. Coronado, E.; Palacio, F.; Veciana, J. Molecule-Based Magnetic Materials. Angew. Chem. Int. Ed. 2003, 42, 2570–2572. [Google Scholar] [CrossRef] [PubMed]
  6. Gatteschi, D.; Sessoli, R.; Villain, J. Molecular Nanomagnets; Oxford University Press: Oxford, UK, 2006. [Google Scholar]
  7. Benelli, C.; Gatteschi, D. Introduction to Molecular Magnetism—From Transition Metals to Lanthanides; Wiley–VCH: Weinheim, Germany, 2015. [Google Scholar]
  8. Banerjee, H.; Chakraborty, S.; Saha-Dasgupta, T. Design and Control of Cooperativity in Spin-Crossover in Metal-Organic Complexes: A Theoretical Overview. Inorganics 2017, 5, 47. [Google Scholar] [CrossRef]
  9. Kumar, K.S.; Ruben, M. Emerging trends in spin crossover (SCO) based functional materials and devices. Coord. Chem. Rev. 2017, 346, 176–205. [Google Scholar] [CrossRef]
  10. Molnár, G.; Rat, S.; Salmon, L.; Nicolazzi, W.; Bousseksou, A. Spin crossover nanomaterials: From fundamental concepts to devices. Adv. Mater. 2018, 30, 1703862. [Google Scholar] [CrossRef] [PubMed]
  11. Auwärter, W.; Écija, D.; Klappenberger, F.; Barth, J.V. Porphyrins at interfaces. Nat. Chem. 2015, 7, 105–120. [Google Scholar] [CrossRef] [PubMed]
  12. Gottfried, J.M. Surface chemistry of porphyrins and phthalocyanines. Surf. Sci. Rep. 2015, 70, 259–379. [Google Scholar] [CrossRef]
  13. Scheybal, A.; Ramsvik, T.; Bertschinger, R.; Putero, M.; Nolting, F.; Jung, T.A. Induced magnetic ordering in a molecular monolayer. Chem. Phys. Lett. 2005, 411, 214–220. [Google Scholar] [CrossRef]
  14. Gambardella, P.; Stepanow, S.; Dmitriev, A.; Honolka, J.; de Groot, F.M.F.; Lingenfelder, M.; Gupta, S.S.; Sarma, D.D.; Bencok, P.; Stanescu, S.; et al. Supramolecular control of the magnetic anisotropy in two-dimensional high-spin Fe arrays at a metal interface. Nat. Mater. 2009, 8, 189–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Miller, J.S.; Gatteschi, D. Molecule-based magnets. Chem. Soc. Rev. 2011, 40, 3065–3066. [Google Scholar] [CrossRef]
  16. Raman, K.V.; Kamerbeek, A.M.; Mukherjee, A.; Atodiresei, N.; Sen, T.K.; Lazić, P.; Caciuc, V.; Michel, R.; Stalke, D.; Mandal, S.K.; et al. Interface-engineered templates for molecular spin memory devices. Nature 2013, 493, 509. [Google Scholar] [CrossRef] [PubMed]
  17. Djeghloul, F.; Ibrahim, F.; Cantoni, M.; Bowen, M.; Joly, L.; Boukari, S.; Ohresser, P.; Bertran, F.; Le Févre, P.; Thakur, P. Direct observation of a highly spin-polarized organic spinterface at room temperature. Sci. Rep. 2013, 3, 1272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Steil, S.; Großmann, N.; Laux, M.; Ruffing, A.; Steil, D.; Wiesenmayer, M.; Mathias, S.; Monti, O.L.A.; Cinchetti, M.; Aeschlimann, M. Spin-dependent trapping of electrons at spinterfaces. Nat. Phys. 2013, 9, 242–247. [Google Scholar] [CrossRef]
  19. Bairagi, K.; Bellec, A.; Repain, V.; Chacon, C.; Girard, Y.; Garreau, Y.; Lagoute, J.; Rousset, S.; Breitwieser, R.; Hu, Y.C. Tuning the magnetic anisotropy at a molecule-metal interface. Phys. Rev. Lett. 2015, 114, 247203. [Google Scholar] [CrossRef] [PubMed]
  20. Barraud, C.; Bouzehouane, K.; Deranlot, C.; Fusil, S.; Jabbar, H.; Arabski, J.; Rakshit, R.; Kim, D.J.; Kieber, C.; Boukari, S. Unidirectional spin-dependent molecule-ferromagnet hybridized states anisotropy in cobalt phthalocyanine based magnetic tunnel junctions. Phys. Rev. Lett. 2015, 114, 206603. [Google Scholar] [CrossRef] [PubMed]
  21. Sanvito, S. Molecular spintronics. Chem. Soc. Rev. 2011, 40, 3336–3355. [Google Scholar] [CrossRef]
  22. Rocha, A.R.; Garcia-Suarez, V.M.; Bailey, S.W.; Lambert, C.J.; Ferrer, J.; Sanvito, S. Towards molecular spintronics. Nat. Mater. 2005, 4, 335. [Google Scholar] [CrossRef]
  23. Wende, H.; Bernien, M.; Luo, J.; Sorg, C.; Ponpandian, N.; Kurde, J.; Miguel, J.; Piantek, M.; Xu, X.; Eckhold, P.; et al. Substrate-induced magnetic ordering and switching of iron porphyrin molecules. Nat. Mater. 2007, 6, 516. [Google Scholar] [CrossRef]
  24. Lach, S.; Altenhof, A.; Tarafder, K.; Schmitt, F.; Ali, M.E.; Vogel, M.; Sauther, J.; Oppeneer, P.M.; Ziegler, C. Metal-organic hybrid interface states of a ferromagnet/organic semiconductor hybrid junction as basis for engineering spin injection in organic spintronics. Adv. Funct. Mater. 2012, 22, 989–997. [Google Scholar] [CrossRef]
  25. Wäckerlin, C.; Tarafder, K.; Siewert, D.; Girovsky, J.; Hählen, T.; Iacovita, C.; Kleibert, A.; Nolting, F.; Jung, T.A.; Oppeneer, P.M.; et al. On-surface coordination chemistry of planar molecular spin systems: Novel magnetochemical effects induced by axial ligands. Chem. Sci. 2012, 3, 3154–3160. [Google Scholar] [CrossRef]
  26. Wäckerlin, C.; Tarafder, K.; Girovsky, J.; Nowakowski, J.; Hählen, T.; Shchyrba, A.; Siewert, D.; Kleibert, A.; Nolting, F.; Oppeneer, P.M.; et al. Ammonia Coordination Introducing a Magnetic Moment in an On-Surface Low-Spin Porphyrin. Angew. Chem. Int. Ed. 2013, 125, 4666–4669. [Google Scholar] [CrossRef] [Green Version]
  27. Cornia, A.; Mannini, M.; Sainctavit, P.; Sessoli, R. Chemical strategies and characterization tools for the organization of single molecule magnets on surfaces. Chem. Soc. Rev. 2011, 40, 3076–3091. [Google Scholar] [CrossRef] [PubMed]
  28. Iacovita, C.; Rastei, M.V.; Heinrich, B.W.; Brumme, T.; Kortus, J.; Limot, L.; Bucher, J.P. Visualizing the spin of individual cobalt-phthalocyanine molecules. Phys. Rev. Lett. 2008, 101, 116602. [Google Scholar] [CrossRef]
  29. Hermanns, C.F.; Tarafder, K.; Bernien, M.; Krüger, A.; Chang, Y.M.; Oppeneer, P.M.; Kuch, W. Magnetic coupling of porphyrin molecules through graphene. Adv. Mater. 2013, 25, 3473–3477. [Google Scholar] [CrossRef] [PubMed]
  30. Lodi Rizzini, A.; Krull, C.; Balashov, T.; Mugarza, A.; Nistor, C.; Yakhou, F.; Sessi, V.; Klyatskaya, S.; Ruben, M.; Stepanow, S.; et al. Exchange Biasing Single Molecule Magnets: Coupling of TbPc2 to Antiferromagnetic Layers. Nano Lett. 2012, 12, 5703–5707. [Google Scholar] [CrossRef] [PubMed]
  31. Girovsky, J.; Nowakowski, J.; Ali, M.E.; Baljozovic, M.; Rossmann, H.R.; Nijs, T.; Aeby, E.A.; Nowakowska, S.; Siewert, D.; Srivastava, G.; et al. Long-range ferrimagnetic order in a two-dimensional supramolecular Kondo lattice. Nat. Commun. 2017, 8, 15388. [Google Scholar] [CrossRef] [Green Version]
  32. Gruber, M.; Berndt, R. Manipulation of Cyclohexene-Based Organic Molecules on Various Metallic Substrates. J. Phys. Chem. C 2016, 120, 18642–18650. [Google Scholar] [CrossRef]
  33. Flechtner, K.; Kretschmann, A.; Steinrück, H.P.; Gottfried, J.M. NO-Induced Reversible Switching of the Electronic Interaction between a Porphyrin-Coordinated Cobalt Ion and a Silver Surface. J. Am. Chem. Soc. 2007, 129, 12110–12111. [Google Scholar] [CrossRef]
  34. Bernien, M.; Miguel, J.; Weis, C.; Ali, M.E.; Kurde, J.; Krumme, B.; Panchmatia, P.M.; Sanyal, B.; Piantek, M.; Srivastava, P.; et al. Tailoring the nature of magnetic coupling of Fe-porphyrin molecules to ferromagnetic substrates. Phys. Rev. Lett. 2009, 102, 047202. [Google Scholar] [CrossRef] [PubMed]
  35. Ballav, N.; Wäckerlin, C.; Siewert, D.; Oppeneer, P.M.; Jung, T.A. Emergence of On-Surface Magnetochemistry. J. Phys. Chem. Lett. 2013, 4, 2303–2311. [Google Scholar] [CrossRef] [Green Version]
  36. Herper, H.C.; Bernien, M.; Bhandary, S.; Hermanns, C.F.; Krüger, A.; Miguel, J.; Weis, C.; Schmitz-Antoniak, C.; Krumme, B.; Bovenschen, D.; et al. Iron porphyrin molecules on Cu (001): Influence of adlayers and ligands on the magnetic properties. Phys. Rev. B 2013, 87, 174425. [Google Scholar] [CrossRef]
  37. Bhandary, S.; Brena, B.; Panchmatia, P.M.; Brumboiu, I.; Bernien, M.; Weis, C.; Krumme, B.; Etz, C.; Kuch, W.; Wende, H.; et al. Manipulation of spin state of iron porphyrin by chemisorption on magnetic substrates. Phys. Rev. B 2013, 88, 024401. [Google Scholar] [CrossRef] [Green Version]
  38. Kar, P.; Yoshida, M.; Shigeta, Y.; Usui, A.; Kobayashi, A.; Minamidate, T.; Matsunaga, N.; Kato, M. Methanol-Triggered Vapochromism Coupled with Solid-State Spin Switching in a Nickel (II)-Quinonoid Complex. Angew. Chem. Int. Ed. 2017, 129, 2385–2389. [Google Scholar] [CrossRef]
  39. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169. [Google Scholar] [CrossRef]
  40. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953. [Google Scholar] [CrossRef]
  41. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. [Google Scholar] [CrossRef]
  42. Ali, M.E.; Sanyal, B.; Oppeneer, P.M. Tuning the magnetic interaction between manganese porphyrins and ferromagnetic Co substrate through dedicated control of the adsorption. J. Phys. Chem. C 2009, 113, 14381–14383. [Google Scholar] [CrossRef]
  43. Tarafder, K.; Kanungo, S.; Oppeneer, P.M.; Saha-Dasgupta, T. Pressure and temperature control of spin-switchable metal-organic coordination polymers from ab initio calculations. Phys. Rev. Lett. 2012, 109, 077203. [Google Scholar] [CrossRef]
  44. Ali, M.E.; Sanyal, B.; Oppeneer, P.M. Electronic Structure, Spin-States, and Spin-Crossover Reaction of Heme-Related Fe-Porphyrins: A Theoretical Perspective. J. Phys. Chem. B 2012, 116, 5849–5859. [Google Scholar] [CrossRef] [PubMed]
  45. Lebègue, S.; Pillet, S.; Ángyán, J.G. Modeling spin-crossover compounds by periodic DFT+U approach. Phys. Rev. B 2008, 78, 024433. [Google Scholar] [CrossRef]
  46. Panchmatia, P.M.; Ali, M.E.; Sanyal, B.; Oppeneer, P.M. Halide Ligated Iron Porphines: A DFT+U and UB3LYP Study. J. Phys. Chem. A 2010, 114, 13381–13387. [Google Scholar] [CrossRef] [PubMed]
  47. Maldonado, P.; Kanungo, S.; Saha-Dasgupta, T.; Oppeneer, P.M. Two-step spin-switchable tetranuclear Fe(II) molecular solid: Ab initio theory and predictions. Phys. Rev. B 2013, 88, 020408. [Google Scholar] [CrossRef]
  48. Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D.C.; Lundqvist, B.I. Van der Waals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92, 246401. [Google Scholar] [CrossRef] [PubMed]
  49. Tkatchenko, A.; Scheffler, M. Accurate Molecular Van Der Waals Interactions from Ground-State Electron Density and Free-Atom Reference Data. Phys. Rev. Lett. 2009, 102, 073005. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The optimized structure of the nickel quinonoid molecule. Ni, O, N, C, and H atoms are represented through the orange, red, blue, black, and magenta colored balls, respectively.
Figure 1. The optimized structure of the nickel quinonoid molecule. Ni, O, N, C, and H atoms are represented through the orange, red, blue, black, and magenta colored balls, respectively.
Magnetochemistry 05 00002 g001
Figure 2. Top and side view of the chemisorbed NiQ molecule on the surface of Co(001), in the HOLLOW configuration. The NiQ atoms are represented in the same way as in Figure 1, and the Co atoms are represented by green and gray colored spheres. To make the Co positions in the second layer (hollow site) clearer, these Co atoms are represented with gray color and their bonds with other layer Co atoms are not shown.
Figure 2. Top and side view of the chemisorbed NiQ molecule on the surface of Co(001), in the HOLLOW configuration. The NiQ atoms are represented in the same way as in Figure 1, and the Co atoms are represented by green and gray colored spheres. To make the Co positions in the second layer (hollow site) clearer, these Co atoms are represented with gray color and their bonds with other layer Co atoms are not shown.
Magnetochemistry 05 00002 g002
Figure 3. Top and side view of the chemisorbed NiQ molecule on the surface of Co(001), in the TOP configuration. The color code of the atoms is done in the same way as in Figure 2.
Figure 3. Top and side view of the chemisorbed NiQ molecule on the surface of Co(001), in the TOP configuration. The color code of the atoms is done in the same way as in Figure 2.
Magnetochemistry 05 00002 g003
Figure 4. Calculated relative total energy ( Δ E) versus the molecule–surface distance. The total-energy curve clearly displays its energy minima at 1.8 Å, which corresponds to chemisorption of the NiQ molecule on the Co(001) substrate in HOLLOW configuration.
Figure 4. Calculated relative total energy ( Δ E) versus the molecule–surface distance. The total-energy curve clearly displays its energy minima at 1.8 Å, which corresponds to chemisorption of the NiQ molecule on the Co(001) substrate in HOLLOW configuration.
Magnetochemistry 05 00002 g004
Figure 5. d-orbital projected DOS of Ni atom in: (a) gas phase NiQ; (b) chemisorbed in HOLLOW configuration; and (c) chemisorbed in TOP configuration. Spin-up and spin-down DOS are shown by positive and negative values, respectively.
Figure 5. d-orbital projected DOS of Ni atom in: (a) gas phase NiQ; (b) chemisorbed in HOLLOW configuration; and (c) chemisorbed in TOP configuration. Spin-up and spin-down DOS are shown by positive and negative values, respectively.
Magnetochemistry 05 00002 g005
Figure 6. Atom-projected partial DOS of chemisorbed NiQ molecule on Co(001): (a) NiQ adsorbed in HOLLOW configuration; and (b) NiQ adsorbed in TOP configuration. Spin-up and spin-down DOS are shown by positive and negative values, respectively.
Figure 6. Atom-projected partial DOS of chemisorbed NiQ molecule on Co(001): (a) NiQ adsorbed in HOLLOW configuration; and (b) NiQ adsorbed in TOP configuration. Spin-up and spin-down DOS are shown by positive and negative values, respectively.
Magnetochemistry 05 00002 g006
Figure 7. Computed charge density (depicted by the light green hypersurface) of NiQ: (left) in HOLLOW configuration; and (right) in TOP configuration. The charge densities reveal that hybridization between the surface and molecule occurs only through the quinonoid ligands.
Figure 7. Computed charge density (depicted by the light green hypersurface) of NiQ: (left) in HOLLOW configuration; and (right) in TOP configuration. The charge densities reveal that hybridization between the surface and molecule occurs only through the quinonoid ligands.
Magnetochemistry 05 00002 g007
Figure 8. Partial magnetization density computed for an energy window indicated with two arrow marks in the spin down channel of the PDOS plot for the HOLLOW configuration (Figure 6a). The magnetization density clearly shows the ferromagnetic coupling between molecule and surface, and also that an indirect superexchange mechanism is more dominant than the direct exchange interaction. The bright yellow and light blue hypersurfaces depict the positive and negative magnetization densities, respectively.
Figure 8. Partial magnetization density computed for an energy window indicated with two arrow marks in the spin down channel of the PDOS plot for the HOLLOW configuration (Figure 6a). The magnetization density clearly shows the ferromagnetic coupling between molecule and surface, and also that an indirect superexchange mechanism is more dominant than the direct exchange interaction. The bright yellow and light blue hypersurfaces depict the positive and negative magnetization densities, respectively.
Magnetochemistry 05 00002 g008
Figure 9. Computed total magnetization density of NiQ molecule chemisorbed in the TOP configuration. The magnetization density on ligand N and C atoms points to a ferromagnetic coupling between the molecule and substrate through a direct exchange path. The bright yellow and light blue hypersurfaces depict the positive and negative magnetization densities, respectively.
Figure 9. Computed total magnetization density of NiQ molecule chemisorbed in the TOP configuration. The magnetization density on ligand N and C atoms points to a ferromagnetic coupling between the molecule and substrate through a direct exchange path. The bright yellow and light blue hypersurfaces depict the positive and negative magnetization densities, respectively.
Magnetochemistry 05 00002 g009
Table 1. The Ni–O and Ni–N bond lengths (in Å), magnetic moment on Ni, and total molecular moment (in μ B ) of NiQ in gas phase, and chemisorbed in HOLLOW and TOP configurations.
Table 1. The Ni–O and Ni–N bond lengths (in Å), magnetic moment on Ni, and total molecular moment (in μ B ) of NiQ in gas phase, and chemisorbed in HOLLOW and TOP configurations.
Gas phaseHOLLOWTOP
Ni–O11.85622.11512.1156
Ni–O21.85582.11582.1166
Ni–N11.87822.02901.8778
Ni–N21.87992.02951.8778
Average1.86832.07241.9967
Magnetic Ni moment0.0001.1781.133
Magnetic moment0.0001.5611.337

Share and Cite

MDPI and ACS Style

Reddy, I.R.; Oppeneer, P.M.; Tarafder, K. Interfacial Spin Manipulation of Nickel-Quinonoid Complex Adsorbed on Co(001) Substrate. Magnetochemistry 2019, 5, 2. https://doi.org/10.3390/magnetochemistry5010002

AMA Style

Reddy IR, Oppeneer PM, Tarafder K. Interfacial Spin Manipulation of Nickel-Quinonoid Complex Adsorbed on Co(001) Substrate. Magnetochemistry. 2019; 5(1):2. https://doi.org/10.3390/magnetochemistry5010002

Chicago/Turabian Style

Reddy, Indukuru Ramesh, Peter M. Oppeneer, and Kartick Tarafder. 2019. "Interfacial Spin Manipulation of Nickel-Quinonoid Complex Adsorbed on Co(001) Substrate" Magnetochemistry 5, no. 1: 2. https://doi.org/10.3390/magnetochemistry5010002

APA Style

Reddy, I. R., Oppeneer, P. M., & Tarafder, K. (2019). Interfacial Spin Manipulation of Nickel-Quinonoid Complex Adsorbed on Co(001) Substrate. Magnetochemistry, 5(1), 2. https://doi.org/10.3390/magnetochemistry5010002

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop