Source author record

Stefano Sanvito

Stefano Sanvito appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

62works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

62 published item(s)

preprint2026arXiv

Agentic Design of Compositional Descriptors via Autoresearch for Materials Science Applications

Autoresearch offers a flexible paradigm for automating scientific tasks, in which an AI agent proposes, implements, evaluates, and refines candidate solutions against a quantitative objective. Here, we use composition-based materials-property prediction to test whether such agents can perform a task beyond model selection and hyperparameter optimization: the design of input descriptors. We introduce Automat, an autoresearch framework where a coding agent based on a large language model generates composition-only descriptors for chemical compounds and evaluates them using a random forest workflow. The agent is restricted to information derivable from chemical formulas and iteratively proposes, implements, and tests chemically motivated descriptor strategies. We apply Automat, with OpenAI Codex using GPT-5.5 as the coding agent, to the prediction of experimental band gaps in inorganic materials and Curie temperatures in ferromagnetic compounds. In both tasks, Automat improves over fractional-composition, Magpie, and combined fractional-composition/Magpie baselines, while producing descriptor families that are chemically interpretable. These results provide a demonstration that autoresearch agents can generate competitive, task-specific materials descriptors without manual feature engineering during the run. They also reveal current limitations, including descriptor redundancy, sensitivity to greedy feature expansion, and the need for explicit complexity control, descriptor pruning, and more sophisticated search strategies.

preprint2022arXiv

A spectral-neighbour representation for vector fields: machine-learning potentials including spin

We introduce a translational and rotational invariant local representation for vector fields, which can be employed in the construction of machine-learning energy models of solids and molecules. This allows us to describe, on the same footing, the energy fluctuations due to the atomic motion, the longitudinal and transverse excitations of the vector field, and their mutual interplay. The formalism can then be applied to physical systems where the total energy is determined by a vector density, as in the case of magnetism. Our representation is constructed over the power spectrum of the combined angular momentum describing the local atomic positions and the vector field, and can be used in conjunction with different machine-learning schemes and data taken from accurate ab initio electronic structure theories. We demonstrate the descriptive power of our representation for a range of classical spin Hamiltonian and machine-learning algorithms. In particular, we construct energy models based on both linear Ridge regression, as in conventional spectral neighbour analysis potentials, and gaussian approximation. These are both built to represent a Heisenberg-type Hamiltonian including a longitudinal energy term and spin-lattice coupling.

preprint2022arXiv

Data-Driven Time Propagation of Quantum Systems with Neural Networks

We investigate the potential of supervised machine learning to propagate a quantum system in time. While Markovian dynamics can be learned easily, given a sufficient amount of data, non-Markovian systems are non-trivial and their description requires the memory knowledge of past states. Here we analyse the feature of such memory by taking a simple 1D Heisenberg model as many-body Hamiltonian, and construct a non-Markovian description by representing the system over the single-particle reduced density matrix. The number of past states required for this representation to reproduce the time-dependent dynamics is found to grow exponentially with the number of spins and with the density of the system spectrum. Most importantly, we demonstrate that neural networks can work as time propagators at any time in the future and that they can be concatenated in time forming an autoregression. Such neural-network autoregression can be used to generate long-time and arbitrary dense time trajectories. Finally, we investigate the time resolution needed to represent the system memory. We find two regimes: for fine memory samplings the memory needed remains constant, while longer memories are required for coarse samplings, although the total number of time steps remains constant. The boundary between these two regimes is set by the period corresponding to the highest frequency in the system spectrum, demonstrating that neural network can overcome the limitation set by the Shannon-Nyquist sampling theorem.

preprint2022arXiv

Fe- and Co-based magnetic tunnel junctions with AlN and ZnO spacers

AlN and ZnO, two wide band-gap semiconductors extensively used in the display industry, crystallise in the wurtzite structure, which can favour the formation of epitaxial interfaces to close-packed common ferromagnets. Here we explore these semiconductors as material for insulating barriers in magnetic tunnel junctions. In particular, the {\it ab initio} quantum transport code {\it Smeagol} is used to model the $X$[111]/$Y$[0001]/$X$[111] ($X=$ Co and Fe, $Y=$ AlN and ZnO) family of junctions. Both semiconductors display a valance-band top with $p$-orbital character, while the conduction band bottom exhibits $s$-type symmetry. The smallest complex-band decay coefficient in the forbidden energy-gap along the [0001] direction is associated with the $Δ_1$ symmetry, and connects across the band gap at the $Γ$ point in 2D Brillouin zones. This feature enables spin filtering and may result in a large tunnelling magnetoresistance. In general, we find that Co-based junctions present limited spin filtering and little magnetoresistance at low bias, since both spin sub-bands cross the Fermi level with $Δ_1$ symmetry. This contrasts the situation of Fe, where only the minority $Δ_1$ band is available. However, even in the case of Fe the magnitude of the magnetoresistance at low bias remains relatively small, mostly due to conduction away from the $Γ$ point and through complex bands with symmetry different than $Δ_1$. The only exception is for the Fe/AlN/Fe junction, where we predict a magnetoresitance of around 1,000\% at low bias.

preprint2022arXiv

Inversion of the chemical environment representations

Machine-learning generative methods for material design are constructed by representing a given chemical structure, either a solid or a molecule, over appropriate atomic features, generally called structural descriptors. These must be fully descriptive of the system, must facilitate the training process and must be invertible, so that one can extract the atomic configurations corresponding to the output of the model. In general, this last requirement is not automatically satisfied by the most efficient structural descriptors, namely the representation is not directly invertible. Such drawback severely limits our freedom of choice in selecting the most appropriate descriptors for the problem, and thus our flexibility to construct generative models. In this work, we present a general optimization method capable of inverting any local many-body descriptor of the chemical environment, back to a cartesian representation. The algorithm is then implemented together with the bispectrum representation of the local structure and demonstrated for a number of molecules. The scheme presented here, thus, represents a general approach to the inversion of structural descriptors, enabling the construction of efficient structural generative models.

preprint2022arXiv

Spin transfer torque in Mn$_3$Ga-based ferrimagnetic tunnel junctions from first principles

We report on first-principles calculations of spin-transfer torque (STT) in epitaxial magnetic tunnel junctions (MTJs) based on ferrimagnetic tetragonal Mn$_3$Ga electrodes, both as analyzer in an Fe/MgO stack, and also in an analogous stack with a second Mn$_3$Ga electrode (instead of Fe) as polarizer. Solving the ballistic transport problem (NEGF + DFT) for the nonequilibrium spin density in a scattering region extended to over 7.6 nm into the Mn$_3$Ga electrode, we find long-range spatial oscillations of the STT decaying on a length scale of a few tens of angstroms, both in the linear response regime and for finite bias. The oscillatory behavior of the STT in Mn$_3$Ga is robust against variations in the stack geometry and the applied bias voltage, which may affect the phase and the amplitude of the spatial oscillation, but the wave number is only responsive to variations in the longitudinal lattice constant of Mn$_3$Ga (for fixed in-plane geometry) without being commensurate with the lattice. Our interpretation of the long-range STT oscillations is based on the bulk electronic structure of Mn$_3$Ga, taking also into account the spin-filtering properties of the MgO barrier. Comparison to a fully Mn$_3$Ga-based stack shows similar STT oscillations, but a significant enhancement of both the TMR effect at the Fermi level and the STT at the interface, due to resonant tunneling for the mirror-symmetric junction with thinner barrier (three monoatomic layers). From the calculated energy dependence of the spin-polarized transmissions at 0 V, we anticipate asymmetric or symmetric TMR as a function of the applied bias voltage for the Fe-based and the all-Mn$_3$Ga stacks, respectively, which also both exhibit a sign change below 1 V. In the latter (symmetric) case we expect a TMR peak at zero, which is larger for the thinner barriers because of a spin-polarized resonant tunneling contribution.

preprint2022arXiv

Transition between large and small electron polaron at neutral ferroelectric domain walls in BiFeO$_3$

Ferroelectric domain walls are planes within an insulating material that can accumulate and conduct charge carriers, hence the interaction of the domain walls with the charge carriers can be important for photovoltaic and other electronic applications. By means of first principles calculations we predict a transition from a large two-dimensional electron polaron to a small polaron at the domain walls at a critical electron density, with polaron signatures in optical absorption and photoluminescence. We find that large and small polarons at the domain walls create different absorption peaks within the band gap that are not present in the case of pristine domain walls. These are an extended Drude peak in the case of large electron or hole polarons and a narrow mid-gap peak in the case of the small electron polaron.

preprint2021arXiv

On the conservation of the angular momentum in ultrafast spin dynamics

The total angular momentum of a close system is a conserved quantity, which should remain constant in time for any excitation experiment once the pumping signal has extinguished. Such conservation, however, is never satisfied in practice in any real-time first principles description of the demagnetization process. Furthermore, there is a growing experimental evidence that the same takes place in experiments. The missing angular momentum is usually associated to lattice vibrations, which are not measured experimentally and are never considered in real-time simulations. Here we critically analyse the issue and conclude that current state-of-the-art simulations violate angular momentum conservation already at the electronic level of description. This shortcoming originates from an oversimplified description of the spin-orbit coupling, which includes atomic contributions but neglects completely that of itinerant electrons. We corroborate our findings with time-dependent simulations using model tight-binding Hamiltonians, and show that indeed such conservation can be re-introduced by an appropriate choice of spin-orbit coupling. The consequences of our findings on recent experiments are also discussed.

preprint2020arXiv

Edge superconductivity in Multilayer WTe2 Josephson junction

WTe2, as a type-II Weyl semimetal, has 2D Fermi arcs on the (001) surface in the bulk and 1D helical edge states in its monolayer. These features have recently attracted wide attention in condensed matter physics. However, in the intermediate regime between the bulk and monolayer, the edge states have not been resolved owing to its closed band gap which makes the bulk states dominant. Here, we report the signatures of the edge superconductivity by superconducting quantum interference measurements in multilayer WTe2 Josephson junctions and we directly map the localized supercurrent. In thick WTe2 (~60 nm), the supercurrent is uniformly distributed by bulk states with symmetric Josephson effect ($\left|I_c^+(B)\right|=\left|I_c^-(B)\right|$). In thin WTe2 (10 nm), however, the supercurrent becomes confined to the edge and its width reaches up to 1.4 um and exhibits non-symmetric behavior $\left|I_c^+(B)\right|\neq \left|I_c^-(B)\right|$. The ability to tune the edge domination by changing thickness and the edge superconductivity establishes WTe2 as a promising topological system with exotic quantum phases and a rich physics.

preprint2020arXiv

Multiple Spin-Phonon Relaxation Pathways in a Kramer Single-Ion Magnet

We present a first-principles investigation of spin-phonon relaxation in a molecular crystal of Co(II) single-ion magnets. Our study combines electronic structure calculations with machine-learning force fields and unravels the nature of both the Orbach and the Raman relaxation channels in terms of atomistic processes. We find that although both mechanisms are mediated by the excited spin states, the low temperature spin dynamics is dominated by phonons in the THz energy range, which partially suppress the benefit of having a large magnetic anisotropy. This study also determines the importance of intra-molecular motions for both the relaxation mechanisms and paves the way to the rational design of a new generation of single-ion magnets with tailored spin-phonon coupling.

preprint2020arXiv

Role of longitudinal fluctuations in L$1_0$ FePt

L$1_0$ FePt is a technologically important material for a range of novel data storage applications. In the ordered FePt structure the normally non-magnetic Pt ion acquires a magnetic moment, which depends on the local field originating from the neighboring Fe atoms. In this work a model of FePt is constructed, where the induced Pt moment is simulated by using combined longitudinal and rotational spin dynamics. The model is parameterized to include a linear variation of the moment with the exchange field, so that at the Pt site the magnetic moment depends on the Fe ordering. The Curie temperature of FePt is calculated and agrees well with similar models that incorporate the Pt dynamics through an effective Fe-only Hamiltonian. By computing the dynamic correlation function the anisotropy field and the Gilbert damping are extracted over a range of temperatures. The anisotropy exhibits a power-law dependence with temperature with exponent $n\approx2.1$. This agrees well with what observed experimentally and it is obtained without including a two-ion anisotropy term as in other approaches. Our work shows that incorporating longitudinal fluctuations into spin dynamics calculations is crucial for understanding the properties of materials with induced moments.

preprint2020arXiv

The limit of spin lifetime in solid-state electronic spins

The development of spin qubits for quantum technologies requires their protection from the main source of finite-temperature decoherence: atomic vibrations. Here we eliminate one of the main barriers to the progress in this field by providing a complete first-principles picture of spin relaxation that includes up to two-phonon processes. Our method is based on machine learning and electronic structure theory and makes the prediction of spin lifetime in realistic systems feasible. We study a prototypical vanadium-based molecular qubit and reveal that the spin lifetime at high temperature is limited by Raman processes due to a small number of THz intra-molecular vibrations. These findings effectively change the conventional understanding of spin relaxation in this class of materials and open new avenues for the rational design of long-living spin systems.

preprint2019arXiv

Electronic spin-spin decoherence contribution in molecular qubits by quantum unitary spin dynamics

The realisation of quantum computers based on molecular electronic spins requires the design of qubits with very long coherence times, T2. Dephasing can proceed over several different microscopic pathways, active at the same time and in different regimes. This makes the rationalisation of the dephasing process not straightforward. Here we present a computational methodology able to address spin decoherence processes for a general ensemble of spins. The method consists in the propagation of the unitary quantum spin dynamics on a reduced Hilbert space. Then we study the dependence of spin dephasing over the magnetic dilution for a crystal of Vanadyl-based molecular qubits. Our results show the importance of long-range electronic spin-spin interactions and their effect on the shape of the spin-echo signal.

preprint2018arXiv

XDFT: an efficient first-principles method for neutral excitations in molecules

State-of-the-art methods for calculating neutral excitation energies are typically demanding and limited to single electron-hole pairs and their composite plasmons. Here we introduce excitonic density-functional theory (XDFT) a computationally light, generally applicable, first-principles technique for calculating neutral excitations based on generalized constrained DFT. In order to simulate an M-particle excited state of an N-electron system, XDFT automatically optimizes a constraining potential to confine N-M electrons within the ground-state Kohn-Sham valence subspace. We demonstrate the efficacy of XDFT by calculating the lowest single-particle singlet and triplet excitation energies of the well-known Thiel molecular test set, with results which are in excellent agreement with time-dependent DFT. Furthermore, going beyond the capability of adiabatic time-dependent DFT, we show that XDFT can successfully capture double excitations. Overall our method makes optical gaps, excition bindings and oscillator strengths readily accessible at a computational cost comparable to that of standard DFT. As such, XDFT appears as an ideal candidate to work within high-throughput discovery frameworks and within linear-scaling methods for large systems.

preprint2016arXiv

Observation of quasi-two-dimensional Dirac fermions in ZrTe5

Since the discovery of graphene, layered materials have attracted extensive interests owing to their unique electronic and optical characteristics. Among them, Dirac semimetal, one of the most appealing categories, has been a long-sought objective in layered systems beyond graphene. Recently, layered pentatelluride ZrTe5 was found to host signatures of Dirac semimetal. However, the low Fermi level in ZrTe5 strongly hinders a comprehensive understanding of the whole picture of electronic states through photoemission measurements, especially in the conduction band. Here, we report the observation of Dirac fermions in ZrTe5 through magneto-optics and magneto-transport. By applying magnetic field, we observe a square-root-B dependence of inter-Landau-level resonance and Shubnikov-de Haas (SdH) oscillations with non-trivial Berry phase, both of which are hallmarks of Dirac fermions. The angular-dependent SdH oscillations show a clear quasi-two-dimensional feature with highly anisotropic Fermi surface and band topology, in stark contrast to the 3D Dirac semimetal such as Cd3As2. This is further confirmed by the angle-dependent Berry phase measurements and the observation of bulk quantum Hall plateaus. The unique band dispersion is theoretically understood: the system is at the critical point between a 3D Dirac semimetal and a topological insulator phase. With the confined interlayer dispersion and reducible dimensionality, our work establishes ZrTe5 as an ideal platform for exploring exotic physical phenomena of Dirac fermions.

preprint2015arXiv

An in- and ex-situ TEM study into the oxidisation of titanium (IV) sulphide

The degradation of liquid dispersed titanium (IV) sulphide (TiS$_2$) is studied from the point of view of oxidisation. TiS$_2$ is a layered two-dimensional nanomaterial, with potential for energy storage applications. However, flakes in dispersion were observed to degrade. We examined two oxidisation routes: deionised water at room temperature and oxygen gas within the temperature range 150-350°C. Water was seen to slowly oxidise flakes inwards from the edge, forming an amorphous oxide, a result comparable to the state of flakes which degraded within the original dispersion. Oxygen gas was seen to rapidly oxidise entire flakes to a polycrystalline oxide, when heated 275 °C and above. Degradation was studied with scanning transmission electron microscopy (STEM), energy dispersed x-ray spectroscopy (EDX) and electron energy-loss spectroscopy (EELS). Density Functional Theory calculations have shown oxidisation by both water and gas molecules to be thermodynamically favourable.

preprint2015arXiv

Charge carrier mobility in hybrid halide perovskites

The charge transport properties of hybrid halide perovskites are investigated with a combination of density functional theory including van der Waals interaction and the Boltzmann theory for diffusive transport in the relaxation time approximation. We find the mobility of electrons to be in the range 5-10 cm$^2$V$^{-1}$s$^{-1}$ and that for holes within 1-5 cm$^2$V$^{-1}$s$^{-1}$, where the variations depend on the crystal structure investigated and the level of doping. Such results, in good agreement with recent experiments, set the relaxation time to about 1 ps, which is the time-scale for the molecular rotation at room temperature. For the room temperature tetragonal phase we explore two possible orientations of the organic cations and find that the mobility has a significant asymmetry depending on the direction of the current with respect to the molecular axis. This is due mostly to the way the PbI$_3$ octahedral symmetry is broken. Interestingly we find that substituting I with Cl has minor effects on the mobilities. Our analysis suggests that the carrier mobility is probably not a key factor in determining the high solar-harvesting efficiency of this class of materials.

preprint2015arXiv

Charge transfer energies of benzene physisorbed on a graphene sheet from constrained density functional theory

Constrained density functional theory (CDFT) is used to evaluate the energy level alignment of a benzene molecule as it approaches a graphene sheet. Within CDFT the problem is conveniently mapped onto evaluating total energy differences between different charge-separated states, and it does not consist in determining a quasi-particle spectrum. We demonstrate that the simple local density approximation provides a good description of the level aligmnent along the entire binding curve, with excellent agreement to experiments at an infinite separation and to GW calculations close to the bonding distance. The method also allows us to explore the effects due to the presence of graphene structural defects and of multiple molecules. In general all our results can be reproduced by a classical image charge model taking into account the finite dielectric constant of graphene.

preprint2015arXiv

Designing a fully-compensated half-metallic ferrimagnet

Recent experimental work on Mn2RuxGa demonstrates its potential as a compensated ferrimagnetic half-metal (CFHM).Here we present a set of high-throughput ab initio density functional theory calculations and detailed experimental characterisation, that enable us to correctly describe the nominal Mn2RuxGa thin films, in particular with regard to site-disorder and defects. We then construct models that accurately capture all the key features of the Mn-Ru-Ga system, including magnetic compensation and the spin gap at the Fermi level. We find that electronic doping is neccessary, which is achieved with a Mn/Ga ratio smaller than two. Our study shows how composition and substrate-induced biaxial strain can be combined to design the first room-temperature CFHM.

preprint2015arXiv

Gate-tunable quantum oscillations in ambipolar Cd3As2 thin films

Electrostatic doping in materials can lead to various exciting electronic properties, such as metal-insulator transition and superconductivity, by altering the Fermi level position or introducing exotic phases. Cd3As2, a three-dimensional (3D) analog of graphene with extraordinary carrier mobility, was predicted to be a 3D Dirac semimetal, a feature confirmed by recent experiments. However, most research so far has been focused on metallic bulk materials that are known to possess ultra-high mobility and giant magnetoresistance but limited carrier transport tunability. Here, we report on the first observation of a gate-induced transition from band conduction to hopping conduction in single-crystalline Cd3As2 thin films via electrostatic doping by solid electrolyte gating. The extreme charge doping enables the unexpected observation of p-type conductivity in a 50 nm-thick Cd3As2 thin film grown by molecular beam epitaxy. More importantly, the gate-tunable Shubnikov-de Haas (SdH) oscillations and the temperature-dependent resistance reveal a unique band structure and bandgap opening when the dimensionality of Cd3As2 is reduced. This is also confirmed by our first-principles calculations. The present results offer new insights towards nanoelectronic and optoelectronic applications of Dirac semimetals in general, and provide new routes in the search for the intriguing quantum spin Hall effect in low-dimension Dirac semimetals, an effect that is theoretically predicted but not yet experimentally realized.

preprint2015arXiv

Liquid exfoliation of solvent-stabilised black phosphorus: applications beyond electronics

Few layer black phosphorus is a new two-dimensional material which is of great interest for applications, mainly in electronics. However, its lack of stability severely limits our ability to synthesise and process this material. Here we demonstrate that high-quality, few-layer black phosphorus nanosheets can be produced in large quantities by liquid phase exfoliation in the solvent N-cyclohexyl-2-pyrrolidone (CHP). We can control nanosheet dimensions and have developed metrics to estimate both nanosheet size and thickness spectroscopically. When exfoliated in CHP, the nanosheets are remarkably stable unless water is intentionally introduced. Computational studies show the degradation to occur by reaction with water molecules only at the nanosheet edge, leading to the removal of phosphorus atoms and the formation of phosphine and phosphorous acid. We demonstrate that liquid exfoliated black phosphorus nanosheets are potentially useful in a range of applications from optical switches to gas sensors to fillers for composite reinforcement.

preprint2015arXiv

Raman Characterization of Platinum Diselenide Thin Films

Platinum diselenide (PtSe2) is a newly discovered 2D material which is of great interest for applications in electronics and catalysis. PtSe2 films were synthesized by thermally-assisted selenization of predeposited platinum films and scanning transmission electron microscopy revealed the crystal structure of these films to be 1T. Raman scattering of these films was studied as a function of film thickness, laser wavelength and laser polarization. Eg and A1g Raman active modes were identified using polarization measurements in the Raman setup. These modes were found to display a clear position and intensity dependence with film thickness, for multiple excitation wavelengths, and their peak positions agree with simulated phonon dispersion curves for PtSe2. These results highlight the practicality of using Raman spectroscopy as a prime characterization technique for newly-synthesized 2D materials.

preprint2015arXiv

Revealing the role of organic cations in hybrid halide perovskites CH3NH3PbI3

The hybrid halide perovskite CH$_{3}$NH$_{3}$PbI$_{3}$ has enabled solar cells to reach an efficiency of about 18\%, demonstrating a pace for improvements with no precedents in the solar energy arena. Despite such explosive progress, the microscopic origin behind the success of such material is still debated, with the role played by the organic cations in the light-harvesting process remaining unclear. Here van-der-Waals-corrected density functional theory calculations reveal that the orientation of the organic molecules plays a fundamental role in determining the material electronic properties. For instance, if CH$_{3}$NH$_{3}$ orients along a (011)-like direction, the PbI$_{6}$ octahedral cage will distort and the band gap will become indirect. Our results suggest that molecular rotations, with the consequent dynamical change of the band structure, might be at the origin of the slow carrier recombination and the superior conversion efficiency of CH$_{3}$NH$_{3}$PbI$_{3}$.

preprint2015arXiv

Single atom anisotropic magnetoresistance on a topological insulator surface

We demonstrate single atom anisotropic magnetoresistance on the surface of a topological insulator, arising from the interplay between the helical spin-momentum-locked surface electronic structure and the hybridization of the magnetic adatom states. Our first-principles quantum transport calculations based on density functional theory for Mn on Bi$_2$Se$_3$ elucidate the underlying mechanism. We complement our findings with a two dimensional model valid for both single adatoms and magnetic clusters, which leads to a proposed device setup for experimental realization. Our results provide an explanation for the conflicting scattering experiments on magnetic adatoms on topological insulator surfaces, and reveal the real space spin texture around the magnetic impurity.

preprint2015arXiv

Spin-valve Effect in NiFe/MoS2/NiFe Junctions

Two-dimensional (2D) layered transition metal dichalcogenides (TMDs) have been recently proposed as appealing candidate materials for spintronic applications owing to their distinctive atomic crystal structure and exotic physical properties arising from the large bonding anisotropy. Here we introduce the first MoS2-based spin-valves that employ monolayer MoS2 as the nonmagnetic spacer. In contrast with what expected from the semiconducting band-structure of MoS2, the vertically sandwiched-MoS2 layers exhibit metallic behavior. This originates from their strong hybridization with the Ni and Fe atoms of the Permalloy (Py) electrode. The spin-valve effect is observed up to 240 K, with the highest magnetoresistance (MR) up to 0.73% at low temperatures. The experimental work is accompanied by the first principle electron transport calculations, which reveal an MR of ~ 9% for an ideal Py/MoS2/Py junction. Our results clearly identify TMDs as a promising spacer compound in magnetic tunnel junctions and may open a new avenue for the TMDs-based spintronic applications.

preprint2015arXiv

The role of spin-orbit interaction in the ultrafast demagnetization of small iron clusters

The ultra-fast demagnetization of small iron clusters initiated by an intense optical excitation is studied with time-dependent spin density functional theory (TDSDFT). In particular we investigate the effect of the spin-orbit interaction on the onset of the demagnetization process. It is found that the initial rate of coherent spin loss is proportional to the square of the atomic spin-orbit coupling constant, $λ$. A simplified quantum spin model comprising spin-orbit interaction and a local time-dependent magnetic field is found to be the minimal model able to reproduce our {\it ab initio} results. The model predicts the $λ^2$ dependence of the onset rate of demagnetization when it is solved either numerically or analytically in the linear response limit. Our findings are supported by additional TDSDFT simulations of clusters made of Co and Ni.

preprint2014arXiv

Ab-initio transport across Bismuth Selenide surface barriers

We investigate the effect of potential barriers in the form of step edges on the scattering properties of Bi$_2$Se$_3$(111) topological surface states by means of large-scale ab-initio transport simulations. Our results demonstrate the suppression of perfect backscattering, while all other scattering processes, which do not entail a complete spin and momentum reversal, are allowed. Furthermore, we find that the spin of the surface state develops an out of plane component as it traverses the barrier. Our calculations reveal the existence of quasi-bound states in the vicinity of the surface barriers, which appear in the form of an enhanced density of states in the energy window corresponding to the topological state. For double barriers we demonstrate the formation of quantum well states. To complement our first-principles results we construct a two-dimensional low-energy effective model and show that band bending plays a significant role in the scattering process. Our findings are discussed in the context of a number of recent experimental works.

preprint2014arXiv

Efficient spin injection and giant magnetoresistance in Fe/MoS$_2$/Fe junctions

We demonstrate giant magnetoresistance in Fe/MoS$_2$/Fe junctions by means of \textit{ab-initio} transport calculations. We show that junctions incorporating either a mono- or a bi-layer of MoS$_2$ are metallic and that Fe acts as an efficient spin injector into MoS$_2$ with an efficiency of about 45\%. This is the result of the strong coupling between the Fe and S atoms at the interface. For junctions of greater thickness a maximum magnetoresistance of $\sim$300\% is obtained, which remains robust with the applied bias as long as transport is in the tunneling limit. A general recipe for improving the magnetoresistance in spin valves incorporating layered transition metal dichalcogenides is proposed.

preprint2014arXiv

Multiprobe quantum spin Hall bars

We analyze electron transport in multiprobe quantum spin Hall (QSH) bars using the Büttiker formalism and draw parallels with their quantum Hall (QH) counterparts. We find that in a QSH bar the measured resistance changes upon introducing side voltage probes, in contrast to the QH case. We also study four- and six-terminal geometries and derive the expressions for the resistances. For these our analysis is generalized from the single-channel to the multi-channel case and to the inclusion of backscattering originating from a constriction placed within the bar.

preprint2014arXiv

Phenomenology of Andreev reflection from first-principles transport theory

We study Andreev reflection in normal metal-superconductor junctions by using an extended Blonder-Tinkham-Klapwijk model combined with transport calculations based on density functional theory. Starting from a parameter-free description of the underlying electronic structure, we perform a detailed investigation of normal metal-superconductor junctions, as the separation between the superconductor and the normal metal is varied. The results are interpreted by means of transverse momentum resolved calculations, which allow us to examine the contributions arising from different regions of the Brillouin zone. Furthermore we investigate the effect of a voltage bias on the normal metal-superconductor conductance spectra. Finally, we consider Andreev reflection in carbon nanotubes sandwiched between normal and superconducting electrodes.

preprint2014arXiv

Proximity induced topological state in graphene

The appearance of topologically protected states at the surface of an ordinary insulator is a rare occurrence and to date only a handful of materials are known for having this property. An intriguing question concerns the possibility of forming topologically protected interfaces between different materials. Here we propose that a topological phase can be transferred to graphene by proximity with the three-dimensional topological insulator Bi$_2$Se$_3$. By using density functional and transport theory we prove that, at the verge of the chemical bond formation, a hybrid state forms at the graphene/Bi$_2$Se$_3$ interface. The state has Dirac-cone-like dispersion at the $Γ$ point and a well-defined helical spin-texture, indicating its topologically protected nature. This demonstrates that proximity can transfer the topological phase from Bi$_2$Se$_3$ to graphene.

preprint2014arXiv

Stretching of BDT-gold molecular junctions: thiol or thiolate termination?

It is often assumed that the hydrogen atoms in the thiol groups of a benzene-1,4-dithiol dissociate when Au-benzene-1,4-dithiol-Au junctions are formed. We demonstrate, by stability and transport properties calculations, that this assumption can not be made. We show that the dissociative adsorption of methanethiol and benzene-1,4-dithiol molecules on a flat Au(111) surface is energetically unfavorable and that the activation barrier for this reaction is as high as 1 eV. For the molecule in the junction, our results show, for all electrode geometries studied, that the thiol junctions are energetically more stable than their thiolate counterparts. Due to the fact that density functional theory (DFT) within the local density approximation (LDA) underestimates the energy difference between the lowest unoccupied molecular orbital and the highest occupied molecular orbital by several electron-volts, and that it does not capture the renormalization of the energy levels due to the image charge effect, the conductance of the Au-benzene-1,4-dithiol-Au junctions is overestimated. After taking into account corrections due to image charge effects by means of constrained-DFT calculations and electrostatic classical models, we apply a scissor operator to correct the DFT energy levels positions, and calculate the transport properties of the thiol and thiolate molecular junctions as a function of the electrodes separation.

preprint2014arXiv

Topological tuning in three-dimensional Dirac semimetals

We study with first-principles methods the interplay between bulk and surface Dirac fermions in three dimensional Dirac semimetals. By combining density functional theory with the coherent potential approximation, we reveal a topological phase transition in Na$_3$Bi$_{1-x}$Sb$_{x}$ and Cd$_3$[As$_{1-x}$P$_x$]$_2$ alloys, where the material goes from a Dirac semimetal to a trivial insulator upon changing Sb or P concentrations. Tuning the composition allows us to engineer the position of the bulk Dirac points in reciprocal space. Interestingly, the phase transition coincides with the reversal of the band ordering between the conduction and valence bands.

preprint2013arXiv

Ab-initio study on the possible doping strategies for MoS$_2$ monolayers

Density functional theory is used to systematically study the electronic and magnetic properties of doped MoS$_2$ monolayers, where the dopants are incorporated both via S/Mo substitution or as adsorbates. Among the possible substitutional dopants at the Mo site, Nb is identified as suitable p-type dopant, while Re is the donor with the lowest activation energy. When dopants are simply adsorbed on a monolayer we find that alkali metals shift the Fermi energy into the MoS$_2$ conduction band, making the system n-type. Finally, the adsorption of charged molecules is considered, mimicking an ionic liquid environment. We find that molecules adsorption can lead to both n- and p-type conductivity, depending on the charge polarity of the adsorbed species.

preprint2013arXiv

Atomistic Simulations of Highly Conductive Molecular Transport Junctions Under Realistic Conditions

We report state-of-the-art atomistic simulations combined with high-fidelity conductance calculations to probe structure-conductance relationships in Au-benzenedithiolate(BDT)-Au junctions under elongation. Our results demonstrate that large increases in conductance are associated with the formation of monatomic chains (MACs) of Au atoms directly connected to BDT. An analysis of the electronic structure of the simulated junctions reveals that enhancement in the s-like states in Au MACs causes the increases in conductance. Other structures also result in increased conductance but are too short-lived to be detected in experiment, while MACs remain stable for long simulation times. Examinations of thermally evolved junctions with and without MACs show negligible overlap between conductance histograms, indicating that the increase in conductance is related to this unique structural change and not thermal fluctuation. These results, which provide an excellent explanation for a recently observed anomalous experimental result [Bruot et al., Nature Nanotech., 2012, 7, 35-40], should aid in the development of mechanically responsive molecular electronic devices.

preprint2013arXiv

Detection of the electrostatic spin crossover effect in magnetic molecules

Scanning tunneling microscopy (STM) can be used to detect inelastic spin transitions in magnetic nano-structures comprising only a handful of atoms. Here we demonstrate that STM can uniquely identify the electrostatic spin crossover effect, whereby the exchange interaction between two magnetic centers in a magnetic molecule changes sign as a function of an external electric field. The fingerprint of such effect is a large drop in the differential conductance as the bias increases. Crucially in the case of a magnetic dimer the spin crossover transition inverts the order between the ground state and the first excited state, but does not change their symmetry. This means that at both sides of the conductance drop associated to the spin crossover transition there are two inelastic transition between the same states. The corresponding conductance steps split identically in a magnetic field and provide a unique way to identify the electrostatic spin crossover.

preprint2013arXiv

Dimensionality driven charge density wave instability in TiS$_2$

Density functional theory and density functional perturbation theory are used to investigate the electronic and vibrational properties of TiS$_2$. Within the local density approximation the material is a semi-metal both in the bulk and in the monolayer form. Most interestingly we observe a Kohn anomaly in the bulk phonon dispersion, which turns into a charge density wave instability when TiS$_2$ is thinned to less than four monolayers. Such charge density wave phase can be tuned by compressive strain, which appears to be the control parameter of the instability.

preprint2013arXiv

Dynamical Exchange Interaction From Time-Dependent Spin Density Functional Theory

We report on {\it ab initio} time-dependent spin dynamics simulations for a two-center magnetic molecular complex based on time-dependent non-collinear spin density functional theory. In particular, we discuss how the dynamical behavior of the {\it ab initio} spin-density in the time-domain can be mapped onto a model Hamiltonian based on the classical Heisenberg spin-spin interaction $J\vcr{S}_1\cdot \vcr{S}_2$. By analyzing individual localized-spin trajectories, extracted from the spin-density evolution, we demonstrate a novel method for evaluating the effective Heisenberg exchange coupling constant, $J$, from first principles simulations. We find that $J$, extracted in such a new dynamical way, agrees quantitatively to that calculated by the standard density functional theory broken-symmetry scheme.

preprint2013arXiv

Electric Field Effects on Armchair MoS2 Nanoribbons

{\it Ab initio} density functional theory calculations are performed to investigate the electronic structure of MoS$_2$ armchair nanoribbons in the presence of an external static electric field. Such nanoribbons, which are nonmagnetic and semiconducting, exhibit a set of weakly interacting edge states whose energy position determines the band-gap of the system. We show that, by applying an external transverse electric field, $E_\mathrm{ext}$, the nanoribbons band-gap can be significantly reduced, leading to a metal-insulator transition beyond a certain critical value. Moreover, the presence of a sufficiently high density of states at the Fermi level in the vicinity of the metal-insulator transition leads to the onset of Stoner ferromagnetism that can be modulated, and even extinguished, by $E_\mathrm{ext}$. In the case of bi-layer nanoribbons we further show that the band-gap can be changed from indirect to direct by applying a transverse field, an effect which might be of significance for opto-electronics applications.

preprint2013arXiv

Gate controlled spin pumping at a quantum spin Hall edge

We propose a four-terminal device designed to manipulate by all electrical means the spin of a magnetic adatom positioned at the edge of a quantum spin Hall insulator. We show that an electrical gate, able to tune the interface resistance between a quantum spin Hall insulator and the source and drain electrodes, can switch the device between two regimes: one where the system exhibits spin pumping and the other where the adatom remains in its ground state. This demonstrates an all-electrical route to control single spins by exploiting helical edge states of topological materials.

preprint2013arXiv

Origin of the n-type and p-type conductivity of MoS2 monolayers on a SiO2 substrate

Ab-initio density functional theory calculations are performed to study the electronic properties of a MoS2 monolayer deposited over a SiO2 substrate in the presence of interface impurities and defects. When MoS2 is placed on a defect-free substrate the oxide plays an insignificant role, since the conduction band top and the valence band minimum of MoS2 are located approximately in the middle of the SiO2 band-gap. However, if Na impurities and O dangling bonds are introduced at the SiO2 surface, these lead to localized states, which modulate the conductivity of the MoS2 monolayer from n- to p-type. Our results show that the conductive properties of MoS2 deposited on SiO2 are mainly determined by the detailed structure of the MoS2 /SiO2 interface, and suggest that doping the substrate can represent a viable strategy for engineering MoS2 -based devices.

preprint2013arXiv

Spin-filtering efficiency of ferrimagnetic spinels CoFe2O4 and NiFe2O4

We assess the potential of the ferrimagnetic spinel ferrites CoFe2O4 and NiFe2O4 to act as spin filtering barriers in magnetic tunnel junctions. Our study is based on the electronic structure calculated by means of first-principles density functional theory within different approximations for the exchange correlation energy. We show that, in agreement with previous calculations, the densities of states suggest a lower tunneling barrier for minority spin electrons, and thus a negative spin-filter effect. However, a more detailed analysis based on the complex band-structure reveals that both signs for the spin-filtering efficiency are possible, depending on the band alignment between the electrode and the barrier materials and depending on the specific wave-function symmetry of the relevant bands within the electrode.

preprint2013arXiv

Spin-pumping and inelastic electron tunneling spectroscopy in topological insulators

We demonstrate that a quantum spin Hall current, spontaneously generated at the edge of a two-dimensional topological insulator, acts as a source of spin-pumping for a magnetic impurity with uniaxial anisotropy. One can then manipulate the impurity spin direction by means of an electrical current without using either magnetic electrodes or an external magnetic field. Furthermore we show that the unique properties of the quantum spin Hall topological state have profound effects on the inelastic electron tunneling spectrum of the impurity. For low current densities inelastic spin-flip events do not contribute to the conductance. As a consequence the conductance steps, normally appearing at voltages corresponding to the spin excitations, are completely suppressed. In contrast an intense current leads to spin pumping and generates a transverse component of the impurity spin. This breaks the topological phase yielding to the conductance steps.

preprint2013arXiv

Structural Origins of Conductance Fluctuations in Gold-Thiolate Molecular Transport Junctions

We report detailed atomistic simulations combined with high-fidelity conductance calculations to probe the structural origins of conductance fluctuations in thermally evolving Au-benzene-1,4-dithiolate-Au junctions. We compare the behavior of structurally ideal junctions (electrodes with flat surfaces) to structurally realistic, experimentally representative junctions resulting from break junction simulations. The enhanced mobility of metal atoms in structurally realistic junctions results in significant changes to the magnitude and origin of the conductance fluctuations. Fluctuations are larger by a factor of 2-3 in realistic junctions compared to ideal junctions. Moreover, in junctions with highly deformed electrodes, the conductance fluctuations arise primarily from changes in the Au geometry, in contrast to results for junctions with non-deformed electrodes, where the conductance fluctuations are dominated by changes in the molecule geometry. These results provide important guidance to experimentalists developing strategies to control molecular conductance for device applications, and also to theoreticians invoking simplified structural models of junctions to predict their behavior.

preprint2012arXiv

Andreev reflection in two-dimensional topological insulators with either conserved or broken time-reversal symmetry

We investigate Andreev reflection in two-dimensional heterojunctions formed by a superconductor in contact with a topological insulator ribbon either possessing or breaking time-reversal symmetry. Both classes of topological insulators exhibit perfect Andreev reflection, which is robust against disorder. This is assigned to topologically protected edge states. In the time-reversal symmetric case we show that doping one of the ribbon edges with magnetic impurities suppresses one Andreev channel, while no such suppression is seen in the broken symmetry situation. Based on this observation we suggest a tabletop transport experiment able to distinguish between the two types of topological insulators, which does not involve the direct measurement of the material band structure.

preprint2012arXiv

Bias asymmetry in the conductance profile of magnetic ions on surfaces probed by scanning tunneling microscopy

The conductance profiles of magnetic transition metal atoms, such as Fe, Co and Mn, deposited on surfaces and probed by a scanning tunneling microscope (STM), provide detailed information on the magnetic excitations of such nano-magnets. In general the profiles are symmetric with respect to the applied bias. However a set of recent experiments has shown evidence for inherent asymmetries when either a normal or a spin-polarized STM tip is used. In order to explain such asymmetries here we expand our previously developed perturbative approach to electron-spin scattering to the spin- polarized case and to the inclusion of out of equilibrium spin populations. In the case of a magnetic STM tip we demonstrate that the asymmetries are driven by the non-equilibrium occupation of the various atomic spin-levels, an effect that reminds closely that electron spin-transfer. In contrast when the tip is not spin-polarized such non-equilibrium population cannot be build up. In this circumstance we propose that the asymmetry simply originates from the transition metal ion density of state, which is included here as a non-vanishing real component to the spin-scattering self-energy.

preprint2012arXiv

Coexistance of giant tunneling electroresistance and magnetoresistance in an all-oxide magnetic tunnel junction

We demonstrate with first-principles electron transport calculations that large tunneling magnetoresistance (TMR) and tunneling electroresistance (TER) effects can coexist in an all-oxide device. The TMR originates from the symmetry-driven spin filtering provided by the insulating BaTiO3 barrier to the electrons injected from SrRuO3. In contrast the TER is possible only when a thin SrTiO3 layer is intercalated at one of the SrRuO3/BaTiO3 interfaces. As the complex band-structure of SrTiO3 has the same symmetry than that of BaTiO3, the inclusion of such an intercalated layer does not negatively alter the TMR and in fact increases it. Crucially, the magnitude of the TER also scales with the thickness of the SrTiO3 layer. The SrTiO3 thickness becomes then a single control parameter for both the TMR and the TER effect. This protocol offers a practical way to the fabrication of four-state memory cells.

preprint2012arXiv

Electronic Transport Through EuO Spin Filter Tunnel Junctions

Epitaxial spin filter tunnel junctions based on the ferromagnetic semiconductor europium monoxide, EuO, are investigated by means of density functional theory. In particular, we focus on the spin transport properties of Cu(100)/EuO(100)/Cu(100) junctions. The dependence of the transmission coefficient and the current-voltage curves on the interface spacing and on the EuO thickness is explained in terms of the EuO density of states and the complex band structure. Furthermore we also discuss the relation between the spin transport properties and the Cu-EuO interface geometry. The level alignment of the junction is sensitively affected by the interface spacing, since this determines the charge transfer between EuO and the Cu electrodes. Our calculations indicate that EuO epitaxially grown on Cu can act as a perfect spin filter, with a spin polarization of the current close to 100%, and with both the Eu-5d conduction band and the Eu-4f valence band states contributing to the coherent transport. For epitaxial EuO on Cu a symmetry filtering is observed, with the Δ_1 states dominating the transmission. This leads to a transport gap larger than the fundamental EuO band gap. Importantly the high spin polarization of the current is preserved up to large bias voltages.

preprint2012arXiv

Magnetism and Antiferroelectricity in MgB$_6$

We report on a density functional theory study demonstrating the coexistence of weak ferromagnetism and antiferroelectricity in boron-deficient MgB6. A boron vacancy produces an almost one dimensional extended molecular orbital, which is responsible for the magnetic moment formation. Then, long-range magnetic order can emerge from the overlap of such orbitals above percolation threshold. Although there is a finite density of states at the Fermi level, the localized nature of the charge density causes an inefficient electron screening. We find that the Mg ions can displace from the center of their cubic cage, thus generating electrical dipoles. In the ground state these order in an antiferroelectric configuration. If proved experimentally, this will be the first material without d or f electrons displaying the coexistence of magnetic and electric order.

preprint2012arXiv

The ground state of a spin-crossover molecule calculated by diffusion Monte Carlo

Spin crossover molecules have recently emerged as a family of compounds potentially useful for implementing molecular spintronics devices. The calculations of the electronic properties of such molecules is a formidable theoretical challenge as one has to describe the spin ground state of a transition metal as the legand field changes. The problem is dominated by the interplay between strong electron correlation at the transition metal site and charge delocalization over the ligands, and thus it fits into a class of problems where density functional theory may be inadequate. Furthermore, the crossover activity is extremely sensitive to environmental conditions, which are difficult to fully characterize. Here we discuss the phase transition of a prototypical spin crossover molecule as obtained with diffusion Monte Carlo simulations. We demonstrate that the ground state changes depending on whether the molecule is in the gas or in the solid phase. As our calculation provides a solid benchmark for the theory we then assess the performances of density functional theory. We find that the low spin state is always over-stabilized, not only by the (semi-)local functionals, but even by the most commonly used hybrids (such as B3LYP and PBE0). We then propose that reliable results can be obtained by using hybrid functionals containing about 50% of exact-exchange.

preprint2012arXiv

Topological surface states scattering in Antimony

In this work we study the topologically protected states of the Sb(111) surface by using ab-initio transport theory. In the presence of a strong surface perturbation we obtain standing-wave states resulting from the superposition of spin-polarized surface states. By Fourier analysis, we identify the underlying two dimensional scattering processes and the spin texture. We find evidence of resonant transmission across surface barriers at quantum well states energies and evaluate their lifetimes. Our results are in excellent agreement with experimental findings. We also show that despite the presence of a step edge along a different high symmetry direction, not yet probed experimentally, the surface states exhibit unperturbed transmission around the Fermi energy for states with near to normal incidence.

preprint2011arXiv

A pertubative approach to the Kondo effect in magnetic atoms on nonmagnetic substrates

Recent experimental advances in scanning tunneling microscopy make the measurement of the conductance spectra of isolated and magnetically coupled atoms on nonmagnetic substrates possible. Notably these spectra are characterized by a competition between the Kondo effect and spin-flip inelastic electron tunneling. In particular they include Kondo resonances and a logarithmic enhancement of the conductance at voltages corresponding to magnetic excitations, two features that cannot be captured by second order perturbation theory in the electron-spin coupling. We have now derived a third order analytic expression for the electron-spin self-energy, which can be readily used in combination with the non-equilibrium Green's function scheme for electron transport at finite bias. We demonstrate that our method is capable of quantitative description the competition between Kondo resonances and spin-flip inelastic electron tunneling at a computational cost significantly lower than that of other approaches. The examples of Co and Fe on CuN are discussed in detail.

preprint2011arXiv

Current-induced energy barrier suppression for electromigration from first principles

We present an efficient method for evaluating current-induced forces in nanoscale junctions, which naturally integrates into the non-equilibrium Green's function formalism implemented within density functional theory. This allows us to perform dynamical atomic relaxation in the presence of an electric current while also evaluating the current-voltage characteristics. The central idea consists in expressing the system energy density matrix in terms of Green's functions. In order to validate our implementation we perform a series of benchmark calculations, both at zero and finite bias. Firstly we evaluate the current-induced forces acting over an Al nanowire and compare them with previously published results for fixed geometries. Then we perform structural relaxation of the same wires under bias and determine the critical voltage at which they break. We find that, while a perfectly straight wire does not break at any of the voltages considered, a zigzag wire is more fragile and snaps at 1.4 V, with the Al atoms moving against the electron flow. Finally we demonstrate the capability of our scheme to tackle the electromigration problem by studying the current-induced motion of a single Si atom covalently attached to the sidewall of a (4,4) armchair single-walled carbon nanotube. Our calculations indicate that if Si is attached along the current path, then current-induced forces can induce migration. In contrast, if the bonding site is away from the current path, then the adatom will remain stable regardless of the voltage. An analysis based on decomposing the total force into a wind and an electrostatic component, as well as on a detailed evaluation of the bond currents, shows that this remarkable electromigration phenomenon is due solely to the position-dependent wind force.

preprint2011arXiv

Effects of edge chemistry doping on graphene nanoribbon mobility

Doping of semiconductor is necessary for various device applications. Exploiting chemistry at its reactive edges was shown to be an effective way to dope an atomically thin graphene nanoribbon (GNR) for realizing new devices in recent experiments. The carrier mobility limited by edge doping is studied as a function of the GNR width, doping density, and carrier density by using ab initio density functional and parameterized tight binding simulations combined with the non-equilibrium Green's function formalism for quantum transport. The results indicate that for GNRs wider than about 4nm, the mobility scales approximately linearly with the GNR width, inversely proportional to the edge doping concentration and decreases for an increasing carrier density. For narrower GNRs, dependence of the mobility on the GNR width and carrier density can be qualitatively different.

preprint2011arXiv

Electric field control of valence tautomeric interconversion in Cobalt dioxolene

We demonstrate that the critical temperature for valence tautomeric interconversion in Cobalt dioxolene complexes can be significantly changed when a static electric field is applied to the molecule. This is achieved by effectively manipulating the redox potential of the metallic acceptor forming the molecule. Importantly our accurate density functional theory calculations demonstrate that already a field of 0.1 V/nm, achievable in Stark spectroscopy experiments, can produce a change in the critical temperature for the interconversion of 20 K. Our results indicate a new way for switching on and off the magnetism in a magnetic molecule. This offers the unique chance of controlling magnetism at the atomic scale by electrical means.

preprint2011arXiv

Electrical control of spin dynamics in finite one-dimensional systems

We investigate the possibility of the electrical control of spin transfer in monoatomic chains incorporating spin-impurities. Our theoretical framework is the mixed quantum-classical (Ehrenfest) description of the spin dynamics, in the spirit of the s-d-model, where the itinerant electrons are described by a tight-binding model while localized spins are treated classically. Our main focus is on the dynamical exchange interaction between two well-separated spins. This can be quantified by the transfer of excitations in the form of transverse spin oscillations. We systematically study the effect of an electrostatic gate bias V_g on the interconnecting channel and we map out the long-range dynamical spin transfer as a function of V_g. We identify regions of V_g giving rise to significant amplification of the spin transmission at low frequencies and relate this to the electronic structure of the channel.

preprint2011arXiv

Spin-flip inelastic electron tunneling spectroscopy in atomic chains

We present a theoretical study of the spin transport properties of mono-atomic magnetic chains with a focus on the spectroscopical features of the I-V curve associated to spin-flip processes. Our calculations are based on the s-d model for magnetism with the electron transport treated at the level of the non-equilibrium Green's function formalism. Inelastic spin-flip scattering processes are introduced perturbatively via the first Born approximation and an expression for the associated self-energy is derived. The computational method is then applied to describe the I-V characteristics and its derivatives of one dimensional chains of Mn atoms and the results are then compared to available experimental data. We find a qualitative and quantitative agreement between the calculated and the experimental conductance spectra. Significantly we are able to describe the relative intensities of the spin excitation features in the I-V curve, by means of a careful analysis of the spin transition selection rules associated to the atomic chains.

preprint2010arXiv

Large bias-dependent magnetoresistance in all-oxide magnetic tunnel junctions with a ferroelectric barrier

All-oxide magnetic tunnel junctions (MTJs) incorporating functional materials as insulating barriers have the potential of becoming the founding technology for novel multi-functional devices. We investigate, by first-principles density functional theory, the bias-dependent transport properties of an all-oxide SrRuO3/BaTiO3/SrRuO3 MTJ. This incorporates a BaTiO3 barrier which can be found either in a non-ferroic or in a ferroelectric state. In such an MTJ not only can the tunneling magnetoresistance reach enormous values, but also, for certain voltages, its sign can be changed by altering the barrier electric state. These findings pave the way for a new generation of electrically-controlled magnetic sensors.

preprint2010arXiv

Magnetic interaction of Co ions near the {10\bar{1}0} ZnO surface

Co-doped ZnO is the prototypical dilute magnetic oxide showing many of the characteristics of ferromagnetism. The microscopic origin of the long range order however remains elusive, since the conventional mechanisms for the magnetic interaction, such as super-exchange and double exchange, fail either at the fundamental or at a quantitative level. Intriguingly, there is a growing evidence that defects both in point-like or extended form play a fundamental role in driving the magnetic order. Here we explore one of such possibilities by performing {\it ab initio} density functional theory calculations for the magnetic interaction of Co ions at or near a ZnO \{10$\bar{1}$0\} surface. We find that extended surface states can hybridize with the $e$-levels of Co and efficiently mediate the magnetic order, although such a mechanism is effective only for ions placed in the first few atomic planes near the surface. We also find that the magnetic anisotropy changes at the surface from an hard-axis easy-plane to an easy axis, with an associated increase of its magnitude. We then conclude that clusters with high densities of surfacial Co ions may display blocking temperatures much higher than in the bulk.

preprint2009arXiv

Electron doping and magnetic moment formation in N- and C-doped MgO

The formation of the magnetic moment in C- and N-doped MgO is the result of a delicate interplay between Hund's coupling, hybridization and Jahn-Teller distortion. The balance depends on a number of environmental variables including electron doping. We investigate such a dependence by self-interaction corrected density functional theory and we find that the moment formation is robust with respect to electron doping. In contrast, the local symmetry around the dopant is more fragile and two different geometries can be stabilized. Crucially the magnetic moment is always extremely localized, making any carrier mediated picture of magnetism in d^0 magnets unlikely.

preprint2009arXiv

Electron transport across electrically switchable magnetic molecules

We investigate the electron transport properties of a model magnetic molecule formed by two magnetic centers whose exchange coupling can be altered with a longitudinal electric field. In general we find a negative differential conductance at low temperatures originating from the different scattering amplitudes of the singlet and triplet states. More interestingly, when the molecule is strongly coupled to the leads and the potential drop at the magnetic centers is only weakly dependent on the magnetic configuration, we find that there is a critical voltage V_C at which the current becomes independent of the temperature. This corresponds to a peak in the low temperature current noise. In such limit we demonstrate that the quadratic current fluctuations are proportional to the product between the conductance fluctuations and the temperature.

preprint2009arXiv

The magnetism of wurtzite CoO nanoclusters

The possibility that the apparent room temperature ferromagnetism, often measured in Co-doped ZnO, is due to uncompensated spins at the surface of wurtzite CoO nanoclusters is investigated by means of a combination of density functional theory and Monte Carlo simulations. We find that the critical temperature extracted from the specific heat systematically drops as the cluster size is reduced, regardless of the particular cluster shape. Furthermore the presence of defects, in the form of missing magnetic sites, further reduces $T_\mathrm{C}$. This suggests that even a spinodal decomposed phase is unlikely to sustain room temperature ferromagnetism in ZnO:Co.