Source author record

J. J. Palacios

J. J. Palacios 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

26works
4topics
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

26 published item(s)

preprint2020arXiv

Consistency between ARPES and STM measurements on SmB$_6$

Strongly correlated topological surface states are promising platforms for next-generation quantum applications, but they remain elusive in real materials. The correlated Kondo insulator SmB$_6$ is one of the most promising candidates, with theoretically predicted heavy Dirac surface states supported by transport and scanning tunneling microscopy (STM) experiments. However, a puzzling discrepancy appears between STM and angle-resolved photoemission (ARPES) experiments on SmB$_6$. Although ARPES detects spin-textured surface states, their velocity is an order of magnitude higher than expected, while the Dirac point -- the hallmark of any topological system -- can only be inferred deep within the bulk valence band. A significant challenge is that SmB$_6$ lacks a natural cleavage plane, resulting in ordered surface domains limited to 10s of nanometers. Here we use STM to show that surface band bending can shift energy features by 10s of meV between domains. Starting from our STM spectra, we simulate the full spectral function as an average over multiple domains with different surface potentials. Our simulation shows excellent agreement with ARPES data, and thus resolves the apparent discrepancy between large-area measurements that average over multiple band-shifted domains and atomically-resolved measurements within a single domain.

preprint2019arXiv

Directional bonding explains high conductance values of atomic contacts in bcc metals

Atomic-sized junctions of iron, created by controlled rupture, present unusually high values of conductance compared to other metals. This result is counter-intuitive since, at the nanoscale, body-centered cubic metals are expected to exhibit lower coordination than face-centered cubic metals. In this work, classical molecular dynamics simulations of contact rupture, using an interatomic potential that accounts for directional bonding, yield highly-coordinated stable structures before rupture, unlike an isotropic bonding potential, which results in the expected stable single-atom contacts. Density functional theory electronic transport calculations show that conductance values of these highly coordinated and highly stable structures, can explain the experimentally measured values for conductance of body-centered cubic atomic contacts, thus revealing the important role of directional bonding in these metals.

preprint2016arXiv

Effective mass theory for the anisotropic exciton in 2D crystals: Application to phosphorene

We present a theoretical study of the exciton binding energy for anisotropic two-dimensional crystals. We obtain analytical expressions from variational wave functions in different limits of the screening length to exciton size ratio and compare them with numerical solutions, both variational and exact. As an example, we apply these results to phosphorene, a monolayer of black phosphorous. Aided by density functional theory calculations for the evaluation of the two-dimensional polarizability, our analytical solution for the exciton binding energy gives a result which is very close to the numerical ones and, in turn, is comparable to the experimental value, as recently reported.

preprint2016arXiv

Local electroexfoliation of graphene with a STM tip

Graphite surfaces can be manipulated by several methods to create graphene structures of different shapes and sizes. Scanning tunneling microscopy (STM) can be used to create these structures either through mechanical contact between the tip and the surface or through electro-exfoliation. In the latter, the mechanisms involved in the process of exfoliation with an applied voltage are not fully understood. Here we show how a graphite surface can be locally exfoliated in a systematic manner by applying an electrostatic force with a STM tip at the edge of a terrace, forming triangular flakes several nanometers in length. We demonstrate, through experiments and simulations, how these flakes are created by a two-step process: first a voltage ramp must be applied at the edge of the terrace, and then the tip must be scanned perpendicularly to the edge. Ab-initio electrostatic calculations reveal that the presence of charges on the graphite surface weakens the interaction between layers allowing for exfoliation at voltages in the same range as those used experimentally. Molecular dynamics simulations show that a force applied locally on the edge of a step produces triangular flakes such as those observed under STM. Our results provide new insights towards surface modification that can be extended to other layered materials.

preprint2015arXiv

A theoretical study of the dynamics of atomic hydrogen on graphene bilayers

We present a theoretical study of the dynamics of H atoms adsorbed on graphene bilayers with Bernal stacking. First, through extensive density functional theory calculations, including van der Waals interactions, we obtain the activation barriers involved in the desorption and migration processes of a single H atom. These barriers, along with attempt rates and the energetics of H pairs, are used as input parameters in kinetic Monte Carlo simulations to study the time evolution of an initial random distribution of adsorbed H atoms. The simulations reveal that, at room temperature, H atoms occupy only one sublattice before they completely desorb or form clusters. This sublattice selectivity in the distribution of H atoms may last for sufficiently long periods of time upon lowering the temperature down to 0 C. The final fate of the H atoms, namely, desorption or cluster formation, depends on the actual relative values of the activation barriers which can be tuned by doping. In some cases a sublattice selectivity can be obtained for periods of time experimentally relevant even at room temperature. This result shows the possibility for observation and applications of the ferromagnetic state associated with such distribution.

preprint2015arXiv

Modeling contact formation between atomic-sized gold tips via molecular dynamics

The formation and rupture of atomic-sized contacts is modelled by means of molecular dynamics simulations. Such nano-contacts are realized in scanning tunnelling microscope and mechanically controlled break junction experiments. These instruments routinely measure the conductance across the nano-sized electrodes as they are brought into contact and separated, permitting conductance traces to be recorded that are plots of conductance versus the distance between the electrodes. One interesting feature of the conductance traces is that for some metals and geometric configurations a jump in the value of the conductance is observed right before contact between the electrodes, a phenomenon known as jump-to-contact. This paper considers, from a computational point of view, the dynamics of contact between two gold nano-electrodes. Repeated indentation of the two surfaces on each other is performed in two crystallographic orientations of face-centred cubic gold, namely (001) and (111). Ultimately, the intention is to identify the structures at the atomic level at the moment of first contact between the surfaces, since the value of the conductance is related to the minimum cross-section in the contact region. Conductance values obtained in this way are determined using first principles electronic transport calculations, with atomic configurations taken from the molecular dynamics simulations serving as input structures.

preprint2015arXiv

Resonant transport and electrostatic effects in single-molecule electrical junctions

In this contribution we demonstrate structural control over a transport resonance in HS(CH$_{2}$)$_{n}$[1,4 - C$_{6}$H$_{4}$](CH$_{2}$)$_{n}$SH (n = 1, 3, 4, 6) metal - molecule - metal junctions, fabricated and tested using the scanning tunnelling microscopy-based $I(z)$ method. The Breit-Wigner resonance originates from one of the arene $π$-bonding orbitals, which sharpens and moves closer to the contact Fermi energy as $n$ increases. Varying the number of methylene groups thus leads to a very shallow decay of the conductance with the length of the molecule. We demonstrate that the electrical behaviour observed here can be straightforwardly rationalized by analyzing the effects caused by the electrostatic balance created at the metal-molecule interface. Such resonances offer future prospects in molecular electronics in terms of controlling charge transport over longer distances, and also in single molecule conductance switching if the resonances can be externally gated.

preprint2015arXiv

Susbtrate-Induced Stabilization and Reconstruction of Zigzag Edges in Graphene Nanoislands on Ni(111)

We combine experimental observations by scanning tunneling microscopy (STM) and density functional theory (DFT) to reveal the most stable edge structures of graphene on Ni(111) as well as the role of stacking-driven activation and suppression of edge reconstruction. Depending on the position of the outermost carbon atoms relative to hollow and on-top Ni sites, zigzag edges have very different energies. Triangular graphene nanoislands are exclusively bound by the more stable zigzag hollow edges. In hexagonal nanoislands, which are constrained by geometry to alternate zigzag hollow and zigzag top edges along their perimeter, only the hollow edge is stable whereas the top edges spontaneously reconstruct into the (57) pentagon-heptagon structure. Atomically-resolved STM images are consistent with either top-fcc or top-hcp epitaxial stacking of graphene and Ni sites, with the former being favored by DFT. Finally, we find that there is a one-to-one relationship between the edge type, graphene stacking, and orientation of the graphene islands.

preprint2014arXiv

Anomalous exchange interaction between intrinsic spins in conducting graphene systems

We address the nature and possible observable consequences of singular one-electron states that appear when strong defects are introduced in the metallic family of graphene, namely, metallic carbon nanotubes and nanotori. In its simplest form, after creating two defects on the same sublattice, a state may emerge at the Fermi energy presenting very unusual properties: It is unique, normalizable, and features a wave function equally distributed around both defects. As a result, the exchange coupling between the magnetic moments generated by the two defects is anomalous. The intrinsic spins couple ferromagnetically, as expected, but do not present an antiferromagnetic excited state at any distance. We propose the use of metallic carbon nanotubes as a novel electronic device based on this anomalous coupling between spins which can be useful for the robust transmission of magnetic information at large distances.

preprint2014arXiv

Hydrogenation-induced ferromagnetism on graphite surfaces

We calculate the electronic structure and magnetic properties of hydrogenated graphite surfaces using van der Waals density functional theory (DFT) and model Hamiltonians. We find, as previously reported, that the interaction between hydrogen atoms on graphene favors adsorption on different sublattices along with an antiferromagnetic coupling of the induced magnetic moments. On the contrary, when hydrogenation takes place on the surface of graphene multilayers or graphite (Bernal stacking), the interaction between hydrogen atoms competes with the different adsorption energies of the two sublattices. This competition may result in all hydrogen atoms adsorbed on the same sublattice and, thereby, in a ferromagnetic state for low concentrations. Based on the exchange couplings obtained from the DFT calculations, we have also evaluated the Curie temperature by mapping this system onto an Ising-like model with randomly located spins. Remarkably, the long-range nature of the magnetic coupling in these systems makes the Curie temperature size dependent and larger than room temperature for typical concentrations and sizes.

preprint2014arXiv

Isolation and characterization of few-layer black phosphorus

Isolation and characterization of mechanically exfoliated black phosphorus flakes with a thickness down to two single-layers is presented. A modification of the mechanical exfoliation method, which provides higher yield of atomically thin flakes than conventional mechanical exfoliation, has been developed. We present general guidelines to determine the number of layers using optical microscopy, Raman spectroscopy and transmission electron microscopy in a fast and reliable way. Moreover, we demonstrate that the exfoliated flakes are highly crystalline and that they are stable even in free-standing form through Raman spectroscopy and transmission electron microscopy measurements. A strong thickness dependence of the band structure is found by density functional theory calculations. The exciton binding energy, within an effective mass approximation, is also calculated for different number of layers. Our computational results for the optical gap are consistent with preliminary photoluminescence results on thin flakes. Finally, we study the environmental stability of black phosphorus flakes finding that the flakes are very hydrophilic and that long term exposure to air moisture etches black phosphorus away. Nonetheless, we demonstrate that the aging of the flakes is slow enough to allow fabrication of field-effect transistors with strong ambipolar behavior. Density functional theory calculations also give us insight into the water-induced changes of the structural and electronic properties of black phosphorus.

preprint2014arXiv

Theory of projections with non-orthogonal basis sets: Partitioning techniques and effective Hamiltonians

Here we present a detailed account of the fundamental problems one encounters in projection theory when non-orthogonal basis sets are used for representation of the operators. In particular, we re-examine the use of projection operators in connection with the calculation of projected (or reduced) Green's functions and associated physical quantities such as the local density of states (LDOS), local charge, and conductance. The unavoidable ambiguity in the evaluation of the LDOS and charge is made explicit with the help of simple examples of metallic nanocontacts while the conductance, within certain obvious limits, remains invariant against the type of projection. We also examine the procedure to obtain effective Hamiltonians from reduced Green's functions. For completeness we include a comparison with results obtained with block-orthogonal basis sets where both direct and dual spaces are used.

preprint2013arXiv

Kondo effect and spin quenching in high-spin molecules on metal substrates

Using a state-of-the art combination of density functional theory and impurity solver techniques we present a complete and parameter-free picture of the Kondo effect in the high-spin ($S=3/2$) coordination complex known as Manganese Phthalocyanine adsorbed on the Pb(111) surface. We calculate the correlated electronic structure and corresponding tunnel spectrum and find an asymmetric Kondo resonance, as recently observed in experiments. Contrary to previous claims, the Kondo resonance stems from only one of three possible Kondo channels with origin in the Mn 3d-orbitals, its peculiar asymmetric shape arising from the modulation of the hybridization due to strong coupling to the organic ligand. The spectral signature of the second Kondo channel is strongly suppressed as the screening occurs via the formation of a many-body singlet with the organic part of the molecule. Finally, a spin-1/2 in the 3d-shell remains completely unscreened due to the lack of hybridization of the corresponding orbital with the substrate, hence leading to a spin-3/2 underscreened Kondo effect.

preprint2013arXiv

Topologically protected quantum transport in locally exfoliated bismuth at room temperature

We report electrical conductance measurements of Bi nanocontacts created by repeated tip-surface indentation using a scanning tunneling microscope at temperatures of 4 K and 300 K. As a function of the elongation of the nanocontact we measure robust, tens of nanometers long plateaus of conductance G0 = 2e^2/h at room temperature. This observation can be accounted for by the mechanical exfoliation of a Bi(111) bilayer, a predicted QSH insulator, in the retracing process following a tip-surface contact. The formation of the bilayer is further supported by the additional observation of conductance steps below G0 before break-up at both temperatures. Our finding provides the first experimental evidence of the possibility of mechanical exfoliation of Bi bilayers, of the existence of the QSH phase in a two-dimensional crystal, and, most importantly, of the observation of the QSH phase at room temperature.

preprint2013arXiv

Understanding the structure of the first atomic contact in Gold

We have studied experimentally the phenomena of jump-to-contact (JC) and jump-out-of-contact (JOC) in gold electrodes. JC can be observed at the first contact when the two metals approach each other while JOC occurs in the last contact before breaking. When the indentation depth between the electrodes is limited to a certain value of conductance, a highly reproducible behaviour in the evolution of the conductance can be obtained for hundreds of cycles of formation and rupture. Molecular dynamics simulations of this process show how the two metallic electrodes are shaped into tips of a well-defined crystallographic structure formed through a mechanical annealing mechanism. We report a detailed analysis of the atomic configurations obtained before contact and rupture of these stable structures and obtained their conductance using first-principlesquantum transport calculations. These results help us understand the values of conductance obtained experimentally in the JC and JOC phenomena and improve our understanding of atomic-sized contacts and the evolution of their structural characteristics.

preprint2012arXiv

A critical analysis of vacancy-induced magnetism in mono and bilayer graphene

The observation of intrinsic magnetic order in graphene and graphene-based materials relies on the formation of magnetic moments and a sufficiently strong mutual interaction. Vacancies are arguably considered the primary source of magnetic moments. Here we present an in-depth density functional theory study of the spin-resolved electronic structure of (monoatomic) vacancies in graphene and bilayer graphene. We use two different methodologies: supercell calculations with the SIESTA code and cluster-embedded calculations with the ALACANT package. Our results are conclusive: The vacancy-induced extended $π$ magnetic moments, which present long-range interactions and are capable of magnetic ordering, vanish at any experimentally relevant vacancy concentration. This holds for $σ$-bond passivated and un-passivated reconstructed vacancies, although, for the un-passivated ones, the disappearance of the $π$ magnetic moments is accompanied by a very large magnetic susceptibility. Only for the unlikely case of a full $σ$-bond passivation, preventing the reconstruction of the vacancy, a full value of 1$μ_B$ for the $π$ extended magnetic moment is recovered for both mono and bilayer cases. Our results put on hold claims of vacancy-induced ferromagnetic or antiferromagnetic order in graphene-based systems, while still leaving the door open to $σ$-type paramagnetism.

preprint2012arXiv

Manipulation of the magnetic ground state of individual manganese phtalocyanines by coordination to CO molecules

We show that the magnetic state of individual manganese phthalocyanine (MnPc) molecules on a Bi(110) surface is modified when the Mn2+ center coordinates to CO molecules adsorbed on top. Using scanning tunneling spectroscopy we identified this change in magnetic properties from the broadening of a Kondo-related zero-bias anomaly when the CO-MnPc complex is formed. The original magnetic state can be recovered by selective desorption of individual CO molecules. First principle calculations show that the CO molecule reduces the spin of the adsorbed MnPc from S = 1 to S = 1=2 and strongly modifies the respective screening channels, driving a transition from an underscreened Kondo state to a state of mix-valence.

preprint2012arXiv

Mechanical annealing of metallic electrodes at the atomic scale

The process of creating an atomically defined and robust metallic tip is described and quantified using measurements of contact conductance between gold electrodes and numerical simulations. Our experiments show how the same conductance behavior can be obtained for hundreds of cycles of formation and rupture of the nanocontact by limiting the indentation depth between the two electrodes up to a conductance value of approximately $5G_{0}$ in the case of gold. This phenomenon is rationalized using molecular dynamics simulations together with density functional theory transport calculations which show how, after repeated indentations (mechanical annealing), the two metallic electrodes are shaped into tips of reproducible structure. These results provide a crucial insight into fundamental aspects relevant to nano-tribology or scanning probe microscopies.

preprint2012arXiv

Spin-filtered edge states in graphene

Spin orbit coupling changes graphene, in principle, into a two-dimensional topological insulator, also known as quantum spin Hall insulator. One of the expected consequences is the existence of spin-filtered edge states that carry dissipationless spin currents and undergo no back-scattering in the presence of non-magnetic disorder, leading to quantization of conductance. Whereas, due to the small size of spin orbit coupling in graphene, the experimental observation of these remarkable predictions is unlikely, the theoretical understanding of these spin-filtered states is shedding light on the electronic properties of edge states in other two-dimensional quantum spin Hall insulators. Here we review the effect of a variety of perturbations, like curvature, disorder, edge reconstruction, edge crystallographic orientation, and Coulomb interactions on the electronic properties of these spin filtered states.

preprint2010arXiv

Critical comparison of electrode models in density functional theory based quantum transport calculations

We study the performance of two different electrode models in quantum transport calculations based on density functional theory: Parametrized Bethe lattices and quasi-one dimensional wires or nanowires. A detailed account of implementation details in both cases is given. From the systematic study of nanocontacts made of representative metallic elements, we can conclude that parametrized electrode models represent an excellent compromise between computational cost and electronic structure definition as long as the aim is to compare with experiments where the precise atomic structure of the electrodes is not relevant or defined with precision. The results obtained using parametrized Bethe lattices are essentially similar to the ones obtained with quasi one dimensional electrodes for large enough sections of these, adding a natural smearing to the transmission curves that mimics the true nature of polycrystalline electrodes. The latter are more demanding from the computational point of view, but present the advantage of expanding the range of applicability of transport calculations to situations where the electrodes have a well-defined atomic structure, as is case for carbon nanotubes, graphene nanoribbons or semiconducting nanowires. All the analysis is done with the help of codes developed by the authors which can be found in the quantum transport toolbox Alacant and are publicly available.

preprint2010arXiv

Hydrogenated Graphene Nanoribbons for Spintronics

We show how hydrogenation of graphene nanoribbons at small concentrations can open new venues towards carbon-based spintronics applications regardless of any especific edge termination or passivation of the nanoribbons. Density functional theory calculations show that an adsorbed H atom induces a spin density on the surrounding $π$ orbitals whose symmetry and degree of localization depends on the distance to the edges of the nanoribbon. As expected for graphene-based systems, these induced magnetic moments interact ferromagnetically or antiferromagnetically depending on the relative adsorption graphene sublattice, but the magnitude of the interactions are found to strongly vary with the position of the H atoms relative to the edges. We also calculate, with the help of the Hubbard model, the transport properties of hydrogenated armchair semiconducting graphene nanoribbons in the diluted regime and show how the exchange coupling between H atoms can be exploited in the design of novel magnetoresistive devices.

preprint2010arXiv

Origin of the quasi-universality of the graphene minimal conductivity

It is a fact that the minimal conductivity $σ_0$ of most graphene samples is larger than the well-established universal value for ideal graphene $4e^2/πh$; in particular, larger by a factor $\gtrsimπ$. Despite intense theoretical activity, this fundamental issue has eluded an explanation so far. Here we present fully atomistic quantum mechanical estimates of the graphene minimal conductivity where electron-electron interactions are considered in the framework of density functional theory. We show the first conclusive evidence of the dominant role on the minimal conductivity of charged impurities over ripples, which have no visible effect. Furthermore, in combination with the logarithmic scaling law for diffusive metallic graphene, we ellucidate the origin of the ubiquitously observed minimal conductivity in the range $8e^2/h > σ_0 \gtrsim 4e^2/h$.

preprint2010arXiv

Spin orbit in curved graphene ribbons

We study the electronic properties of electrons in flat and curved zigzag graphene ribbons using a tight-binding model within the Slater Koster approximation. We find that curvature dramatically enhances the action of spin orbit effects in graphene ribbons and has a strong effect on the spin orientation of the edge states: whereas spins are normal to the surface in the case of flat ribbons, this is no longer the case in the case of curved ribbons. We find that for the edge states, the spin density lies always in the plane perpendicular to the ribbon axis, and deviate strongly from the normal to the ribbon, even for very small curvature and the small spin orbit coupling of carbon. We find that curvature results also in an effective second neighbor hopping that modifies the electronic properties of zigzag graphene ribbons. We discuss the implications of our finding in the spin Hall phase of curved graphene Ribbons.

preprint2009arXiv

Spin-transfer torque on a single magnetic adatom

We theoretically show how the spin orientation of a single magnetic adatom can be controlled by spin polarized electrons in a scanning tunneling microscope configuration. The underlying physical mechanism is spin assisted inelastic tunneling. By changing the direction of the applied current, the orientation of the magnetic adatom can be completely reversed on a time scale that ranges from a few nanoseconds to microseconds, depending on bias and temperature. The changes in the adatom magnetization direction are, in turn, reflected in the tunneling conductance.

preprint2007arXiv

Performance limits of graphene-ribbon-based field effect transistors

The performance of field effect transistors based on an single graphene ribbon with a constriction and a single back gate are studied with the help of atomistic models. It is shown how this scheme, unlike that of traditional carbon-nanotube-based transistors, reduces the importance of the specifics of the chemical bonding to the metallic electrodes in favor of the carbon-based part of device. The ultimate performance limits are here studied for various constriction and metal-ribbon contact models. In particular we show that, even for poorly contacting metals, properly taylored constrictions can give promising values for both the on-conductance and the subthreshold swing.

preprint1997arXiv

Bulk charge distributions on integer and fractional quantum Hall plateaus

We discuss the charge distributions across the bulk of a two-dimensional electron gas system which is on an integer or fractional quantum Hall plateau. Our analysis is based on a relation, derived from the long wavelength limit of the bulk density-density response function, between the induced charge and leading derivatives of a slowly varying Hall potential. We use the Wiener-Hopf method to solve the coupled response and Poisson equations. Unlike the integer case, treated in previous work, the induced charge in the fractional case can alternate in sign.