Source author record

Javier Junquera

Javier Junquera 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

23works
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

23 published item(s)

preprint2023arXiv

Emergent chirality in a polar meron to skyrmion phase transition

Polar skyrmions are predicted to emerge from the interplay of elastic, electrostatic and gradient energies, in contrast to the key role of the anti-symmetric Dzyalozhinskii-Moriya interaction in magnetic skyrmions. With the discovery of topologically-stable polar skyrmions, it is of both fundamental and practical interest to understand the microscopic nature and the possibility of temperature- and strain-driven phase transitions in ensembles of such polar skyrmions. Here, we explore the reversible transition from a skyrmion state (topological charge of -1) to a two-dimensional, tetratic lattice of merons (with topological charge of -1/2) upon varying the temperature and elastic boundary conditions in [(PbTiO3)16/(SrTiO3)16]8 lifted-off membranes. This topological phase transition is accompanied by a change in chirality, from zero-net chirality (in meronic phase) to net-handedness (in skyrmionic phase). To map these changes microscopically required developing new imaging methods. We show how scanning convergent beam electron diffraction provides a robust measure of the local polarization simultaneously with the strain state at sub-nm resolution, while also directly mapping the chirality of each skyrmion. Using this, we demonstrate strain as a crucial order parameter to drive isotropic-to-anisotropic structural transitions of chiral polar skyrmions to non-chiral merons, validated with X-ray reciprocal space mapping and theoretical phase-field simulations. These results revealed by our new measurement methods provide the first illustration of systematic control of rich variety of topological dipole textures by altering the mechanical boundary conditions, which may offer a promising way to control their functionalities in ferroelectric nanodevices using the local and spatial distribution of chirality and order.

preprint2022arXiv

Berezinskii-Kosterlitz-Thouless phases in ultra-thin PbTiO$_3$/SrTiO$_3$ superlattices

We study the emergence of Berezinskii-Kosterlitz-Thouless (BKT) phases in (PbTiO$_3$)$_3$/(SrTiO$_3$)$_3$ superlattices by means of second-principles simulations. Beyond a threshold tensile epitaxial strain of $ε= 0.25 \%$ the local dipole moments within the superlattices are confined to the film-plane, and thus the polarization can be effectively considered as two-dimensional. The analysis of the decay of the dipole-dipole correlation with the distance, together with the study of the density of defects and its distribution as function of temperature, supports the existence of a BKT phase in a range of temperature mediating the ordered ferroelectric (stable at low $T$), and the disordered paraelectric phase that appears beyond a critical temperature $T_{\rm BKT}$. This BKT phase is characterized by quasi-long-range order (whose signature is a power-law decay of the correlations with the distance), and the emergence of tightly bounded vortex-antivortex pairs whose density is determined by a thermal activation process. The proposed PbTiO$_{3}$/SrTiO$_{3}$ superlattice model and the imposed mechanical boundary conditions are both experimentally feasible, opening the door for the first experimental observation of these new topological phases in ferroelectric materials.

preprint2022arXiv

Melting of Polarization Vortices Crystals and Chiral Phase Transitions in Oxide Superlattices

We study the equilibrium arrangements of polarization vortices in (PbTiO$_3$)$_n$/(SrTiO$_3$)$_n$ superlattices by means of second-principles simulations. We find that, at low temperatures, polarization vortices organize in a regular arrangement in which clockwise and counter-clockwise vortices alternate positions, leading to a crystal-like structure with well defined handedness. This chiral crystal melts at a critical temperature $T_\mathrm{M}$ into a chiral liquid, where long-range order is lost but handedness is preserved. At even higher temperatures, $T_\mathrm{C}$, a second phase transition occurs, at which the chiral liquid of polarization vortices loses its handedness. Both phase transitions can be readily identified by the adequate choices of order parameters.

preprint2021arXiv

Order-disorder transitions in a polar vortex lattice

Order-disorder transitions are widely explored in various vortex structures in condensed matter physics, i.e., in the type-II superconductors and Bose-Einstein condensates. In this study, we have investigated the ordering of the polar vortex phase in the (PZT)n/(STO)n superlattice systems through phase-field simulations. An antiorder state is discovered for short periodicity superlattice on an SSO substrate, owing to the huge interfacial coupling between PZT and STO as well as the giant in-plane polarization in STO layers due to the large tensile strain. Increasing the periodicity leads to the anti-order to disorder transition, resulting from the loss of interfacial coupling and disappearance of the polarization in STO layers. On the other hand, for short periodicity superlattices, order-disorder-antiorder transition can be engineered by mediating the substrate strain, due to the delicate competition between the depoling effect, interfacial coupling, and strain effect. We envision this study to spur further interest towards the understanding of order-disorder transition in ferroelectric topological structures.

preprint2020arXiv

Comparison of band-fitting and Wannier-based model construction for WSe$_2$

Transition metal dichalcogenide materials $MX_2 (M=Mo,W;X=S,Se)$ are being thoroughly studied due to their novel two-dimensional structure, that is associated with exceptional optical and transport properties. From a computational point of view, Density Functional Theory simulations perform very well in these systems and are an indispensable tool to predict and complement experimental results. However, due to the time and length scales where even the most efficient DFT implementations can reach today, this methodology suffers of stringent limitations to deal with finite temperature simulations or electron-lattice coupling when studying excitation states: the unit cells required to study, for instance, systems with thermal fluctuations or large polarons would require a large computational power. Multi-scale techniques, like the recently proposed Second Principles Density Functional Theory, can go beyond these limitations but require the construction of tight-binding models for the systems under investigation. In this work, we compare two such methods to construct the bands of WSe$_2$. In particular, we compare the result of (i) Wannier-based model construction with (ii) the band fitting method of Liu and co-workers where the top of the valence band and the bottom of the conduction band are modeled by three bands symmetrized to have mainly Tungsten $d_{z^2}$, $d_{xy}$ and $d_{x^2-y^2}$ character. Our results emphasize the differences between these two approaches and how band-fitting model construction leads to an overestimation of the localization of the real-space basis in a tight-binding representation.

preprint2020arXiv

First-principles study of two-dimensional electron and hole gases at the head-to-head and tail-to-tail 180$^\circ$ domain walls in PbTiO$_{3}$ ferroelectric thin films

We study from first-principles the structural and electronic properties of head-to-head (HH) and tail-to-tail (TT) 180$^\circ$ domain walls in isolated free-standing PbTiO$_{3}$ slabs. For sufficiently thick domains ($n$ = 16 unit cells of PbTiO$_{3}$), a transfer of charge from the free surfaces to the domain walls to form localized electron (in the HH) and hole (in the TT) gases in order to screen the bound polarization charges is observed. The electrostatic driving force behind this electronic reconstruction is clearly visible from the perfect match between the smoothed free charge densities and the bound charge distribution, computed from a finite difference of the polarization profile obtained after the relaxation of the lattice degrees of freedom. The domain wall widths, of around six unit cells, are larger than in the conventional 180$^\circ$ neutral configurations. Since no oxygen vacancies, defects or dopant atoms are introduced in our simulations, all the previous physical quantities are the intrinsic limits of the system. Our results support the existence of an extra source of charge at the domains walls to explain the enhancement of the conductivity observed in some domains walls of prototypical, insulating in bulk, perovskite oxides.

preprint2020arXiv

SIESTA: recent developments and applications

A review of the present status, recent enhancements, and applicability of the SIESTA program is presented. Since its debut in the mid-nineties, SIESTA's flexibility, efficiency and free distribution has given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of SIESTA combines finite-support pseudo-atomic orbitals as basis sets, norm-conserving pseudopotentials, and a real-space grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here we describe the more recent implementations on top of that core scheme, which include: full spin-orbit interaction, non-repeated and multiple-contact ballistic electron transport, DFT+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, density-functional perturbation theory, efficient Van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as Wannier90 and the second-principles modelling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering SIESTA runs, and various postprocessing utilities. SIESTA has also been engaged in the Electronic Structure Library effort from its inception, which has allowed the sharing of various low level libraries, as well as data standards and support for them, in particular the PSML definition and library for transferable pseudopotentials, and the interface to the ELSI library of solvers. Code sharing is made easier by the new open-source licensing model of the program. This review also presents examples of application of the capabilities of the code, as well as a view of on-going and future developments.

preprint2017arXiv

The PSML format and library for norm-conserving pseudopotential data curation and interoperability

Norm-conserving pseudopotentials are used by a significant number of electronic-structure packages, but the practical differences among codes in the handling of the associated data hinder their interoperability and make it difficult to compare their results. At the same time, existing formats lack provenance data, which makes it difficult to track and document computational workflows. To address these problems, we first propose a file format (PSML) that maps the basic concepts of the norm-conserving pseudopotential domain in a flexible form and supports the inclusion of provenance information and other important metadata. Second, we provide a software library (libPSML) that can be used by electronic structure codes to transparently extract the information in the file and adapt it to their own data structures, or to create converters for other formats. Support for the new file format has been already implemented in several pseudopotential generator programs (including ATOM and ONCVPSP), and the library has been linked with Siesta and Abinit, allowing them to work with the same pseudopotential operator (with the same local part and fully non-local projectors) thus easing the comparison of their results for the structural and electronic properties, as shown for several example systems. This methodology can be easily transferred to any other package that uses norm-conserving pseudopotentials, and offers a proof-of-concept for a general approach to interoperability.

preprint2016arXiv

Efficient thermoelectric materials using nonmagnetic double perovskites with $d^0$/$d^6$ band filling

Efficient thermoelectric materials should present large Seebeck coefficient, high electrical conductivity and low thermal conductivity. An enhanced Seebeck coefficient can be obtained from materials where the Fermi level can be aligned with a large and narrow peak of the density of states, particularly when a substantial band valley degeneracy occurs. A high electrical conductivity comes as a consequence of large conductive hopping paths between the atoms of the material. Both physical quantities can be decoupled and optimized independently if their origins can be ascribed to different sets of bands. Based on these assumptions, double perovskites A$_2$BB'O$_6$ with $d^0/d^6$ filling for the B and B' metal cations, respectively, have been considered. They provide a desirable band structure with degenerate B-$t_{2g}$ / B'-$e_g$ bands above the Fermi level together with a low thermal conductivity. We have carried out first-principles simulations for various of these nonmagnetic double perovskites and showed that all of them present a large Seebeck coefficient (consequence of the localized and empty $t_{2g}$ states of the B-cation), and large electrical conductivity due to the more spread unoccupied $e_{g}$ band of the B' cation. We have seen that if they can be optimally doped, they could show a figure of merit comparable or even higher than the best $n$-type thermoelectric oxides, such as SrTiO$_{3}$. Different mechanisms to tune the band structure and enhance the thermoelectric figure of merit are explored, including epitaxial strain, hydrostatic pressure, chemical pressure, and external doping. A fully relaxed structure has also been studied, showing that a realistic calculation is necessary to make accurate predictions, but also proving that the main trends shown throughout the paper remain unchanged.

preprint2016arXiv

Giant Optical Polarization Rotation Induced by Spin-Orbit Coupling in Polarons

We have uncovered a giant gyrotropic magneto-optical response for doped ferromagnetic manganite La2/3Ca1/3MnO3 around the near room-temperature paramagnetic-to-ferromagnetic transition. At odds with current wisdom, where this response is usually assumed to be fundamentally fixed by the electronic band structure, we point to the presence of small polarons as the driving force for this unexpected phenomenon. We explain the observed properties by the intricate interplay of mobility, Jahn-Teller effect and spin-orbit coupling of small polarons. As magnetic polarons are ubiquitously inherent to many strongly correlated systems, our results provide an original, general pathway towards the generation of gigantic gyrotropic responses that can be harnessed for nonreciprocal devices that exploit the polarization of light.

preprint2015arXiv

Second-principles method including electron and lattice degrees of freedom

We present a first-principles-based (second-principles) scheme that permits large-scale materials simulations including both atomic and electronic degrees of freedom on the same footing. The method is based on a predictive quantum-mechanical theory, e.g., Density Functional Theory, and its accuracy can be systematically improved at a very modest computational cost. Our approach is based on dividing the electron density of the system into a reference part - typically corresponding to the system's neutral, geometry-dependent ground state - and a deformation part - defined as the difference between the actual and reference densities. We then take advantage of the fact that the bulk part of the system's energy depends on the reference density alone; this part can be efficiently and accurately described by a force field, thus avoiding explicit consideration of the electrons. Then, the effects associated to the difference density can be treated perturbatively with good precision by working in a suitably chosen Wannier function basis. Further, the electronic model can be restricted to the bands of interest. All these features combined yield a very flexible and computationally very efficient scheme. Here we present the basic formulation of this approach, as well as a practical strategy to compute model parameters for realistic materials. We illustrate the accuracy and scope of the proposed method with two case studies, namely, the relative stability of various spin arrangements in NiO and the formation of a two-dimensional electron gas at the interface between band insulators LaAlO$_3$ and SrTiO$_3$. We conclude by discussing ways to overcome the limitations of the present approach (most notably, the assumption of a fixed bonding topology), as well as its many envisioned possibilities and future extensions.

preprint2014arXiv

Atomically precise interfaces from non-stoichiometric deposition

Complex oxide heterostructures display some of the most chemically abrupt, atomically precise interfaces, which is advantageous when constructing new interface phases with emergent properties by juxtaposing incompatible ground states. One might assume that atomically precise interfaces result from stoichiometric growth, but here we show that the most precise control is obtained for non-stoichiometric growth where differing surface energies can be compensated by surfactant-like effects. For the precise growth of Sr$_{n+1}$Ti$_n$O$_{3n+1}$ Ruddlesden-Popper (RP) phases, stoichiometric deposition leads to the loss of the first RP rock-salt double layer, but growing with a strontium-rich surface layer restores the bulk stoichiometry and ordering of the subsurface RP structure. Our results dramatically expand the materials that can be prepared in epitaxial heterostructures with precise interface control---from just the $n=\infty$ end members (perovskites) to the entire RP family---enabling the exploration of novel quantum phenomena at a richer variety of oxide interfaces.

preprint2014arXiv

First-principles study of metal-induced gap states in metal/oxide interfaces and their relation with the complex band structure

We develop a simple model to compute the energy-dependent decay factors of metal-induced gap states in metal/insulator interfaces considering the collective behaviour of all the bulk complex bands in the gap of the insulator. The agreement between the penetration length obtained from the model (considering only bulk properties) and full first-principles simulations of the interface (including explicitly the interfaces) is good. The influence of the electrodes and the polarization of the insulator is analyzed. The method simplifies the process of screening materials to be used in Schootky barriers or in the design of giant tunneling electroresistance and magnetoresistance devices.

preprint2012arXiv

Domain walls in a perovskite oxide with two primary structural order parameters: first-principles study of BiFeO$_3$

We present a first-principles study of ferroelectric domain walls (FE-DWs) in multiferroic BiFeO$_3$ (BFO), a material in which the FE order parameter coexists with anti-ferrodistortive (AFD) modes involving rotations of the O$_6$ octahedra. We find that the energetics of the DWs are dominated by the capability of the domains to match their O$_6$ octahedra rotation patterns at the plane of the wall, so that the distortion of the oxygen groups is minimized. Our results thus indicate that, in essence, it is the discontinuity in the AFD order parameter, and not the change in the electric polarization, what decides which crystallographic planes are most likely to host BFO's FE-DWs. Such a result clearly suggests that the O$_6$ rotational patterns play a primary role in the FE phase of this compound, in contrast with the usual (implicit) assumption that they are subordinated to the FE order parameter. Our calculations show that, for the most favorable cases in BFO, the DW energy amounts to a several tens of mJ/m$^2$, which is higher than what was computed for other ferroelectric perovskites with no O$_6$ rotations. Interestingly, we find that the structure of BFO at the most stable DWs resembles the atomic arrangements that are characteristic of low-lying (meta)stable phases of the material. Further, we argue that our results for the DWs of bulk BFO are related with the nanoscale-twinned structures that Prosandeev et al. [Adv. Funct. Mats. (2012), doi: 10.1002/adfm.201201467] have recently predicted to occur in this compound, and suggest that BFO can be viewed as a polytypic material. Our work thus contributes to shape a coherent picture of the structural variants that BFO can present and the way in which they are related.

preprint2012arXiv

Ferroelectricity in ultrathin film capacitors

Going down to the limit of ultrathin films holds promise for a new generation of devices such as ferroelectric tunnel junctions or resistive memories. However, these length scales also make the devices sensitive to parasitic effects related to miniaturization, and a better understanding of what happens as size is reduced is of practical importance for the future development of these devices. As the experimental advances in materials preparation and characterization have come together with great progress in theoretical modeling of ferroelectrics, both theorists and experimentalists can finally probe the same length and time scales. This allows realtime feedback between theory and experiment, with new discoveries now routinely made both in the laboratory and on the computer. Throughout this chapter, we will highlight the recent advances in density functional theory based modeling and the role it played in our understanding of ultrathin ferroelectrics. We will begin with a brief introduction to ferroelectricity and ferroelectric oxides, followed by an overview of the major theoretical developments. We will then discuss some of the subtleties of ferroelectricity in perovskite oxides, before turning our attention to the main subject of the chapter -- ferroelectricity in ultrathin films. We will discuss in detail the influence of the mechanical, electrical and chemical boundary conditions on the stability of the polar state in a parallel plate capacitor geometry, introducing the notion of depolarization fields that tend to destabilize ferroelectricity. We will look at other ways in which a thin ferroelectric can preserve its polar state, focusing on ferroelectric domains and domain walls. Finally, we will briefly discuss artificially layered ferroelectrics and the potential they hold as tailor-made materials for electronic applications.

preprint2012arXiv

First-principles modeling of the thermoelectric properties of SrTiO$_3$/SrRuO$_3$ superlattices

Using a combination of first-principles simulations, based on the density functional theory and Boltzmann's semiclassical theory, we have calculated the transport and thermoelectric properties of the half-metallic two dimensional electron gas confined in single SrRuO$_3$ layers of SrTiO$_3$/SrRuO$_3$ periodic superlattices. Close to the Fermi energy we find that the semiconducting majority spin channel displays a very large in-plane component of the Seebeck tensor at room temperature, $S$ = 1500 $μ$V/K, and the minority spin channel shows good in-plane conductivity $σ$ = 2.5 (m$Ω$cm)$^{-1}$. However, contrary to the expectation of Hicks and Dresselhaus model about enhanced global thermoelectric properties due to the confinement of the metallic electrons, we find that the total power factor and thermoelectric figure of merit for reduced doping is too small for practical applications. The reason for this failure can be traced back on the electronic structure of the interfacial gas, which departs from the free electron behaviour on which the model was based. The evolution of the electronic structure, electrical conductivity, Seebeck coefficient, and power factor as a function of the chemical potential is explained by a simplified tight-binding model. We find that the electron gas in our system is composed by a pair of one dimensional electron gases orthogonal to each other. This reflects the fact the physical dimensionality of the electronic system can be even smaller than that of the spacial confinement of the carriers.

preprint2012arXiv

First-principles simulations on the structural and energetic properties of domains in PbTiO$_3$/SrTiO$_3$ superlattices

We report first-principles calculations, within the density functional theory, on the structural and energetic properties of 180$^\circ$ stripe domains in (PbTiO$_{3}$)$_{n}$/(SrTiO$_{3}$)$_{n}$ superlattices. For the explored periodicities ($n$=3 and 6) we find that the polydomain structures compete in energy with the monodomain phases. Our results suggest the progressive transition, as a function of $n$, from a strong to a weak electrostatic coupling regime between the SrTiO$_{3}$ and PbTiO$_{3}$ layers. Structurally, they display continuous rotation of polarization connecting 180$^\circ$ domains. A large offset between [100] atomic rows across the domain wall and huge strain gradients are observed. The domain wall energy is very isotropic, depending very weakly on the stripe orientation.

preprint2012arXiv

Highly-confined spin-polarized two-dimensional electron gas in SrTiO$_{3}$/SrRuO$_{3}$ superlattices

We report first principles characterization of the structural and electronic properties of (SrTiO$_{3}$)$_{5}$/(SrRuO$_{3}$)$_{1}$ superlattices. We show that the system exhibits a spin-polarized two-dimensional electron gas, extremely confined to the 4$d$ orbitals of Ru in the SrRuO$_{3}$ layer. Every interface in the superlattice behaves as a minority-spin half-metal ferromagnet, with a magnetic moment of $μ$ = 2.0 $μ_{\rm B}$/SrRuO$_{3}$ unit. The shape of the electronic density of states, half metallicity and magnetism are explained in terms of a simplified tight-binding model, considering only the $t_{2g}$ orbitals plus (i) the bi-dimensionality of the system, and (ii) strong electron correlations.

preprint2011arXiv

Band alignment at metal/ferroelectric interfaces: insights and artifacts from first principles

Based on recent advances in first-principles theory, we develop a general model of the band offset at metal/ferroelectric interfaces. We show that, depending on the polarization of the film, a pathological regime might occur where the metallic carriers populate the energy bands of the insulator, making it metallic. As the most common approximations of density functional theory are affected by a systematic underestimation of the fundamental band gap of insulators, this scenario is likely to be an artifact of the simulation. We provide a number of rigorous criteria, together with extensive practical examples, to systematically identify this problematic situation in the calculated electronic and structural properties of ferroelectric systems. We discuss our findings in the context of earlier literature studies, where the issues described in this work have often been overlooked. We also discuss formal analogies to the physics of polarity compensation at LaAlO3/SrTiO3 interfaces, and suggest promising avenues for future research.

preprint2011arXiv

Interplay of couplings between antiferrodistortive, ferroelectric, and strain degrees of freedom in monodomain PbTiO$_{3}$/SrTiO$_{3}$ superlattices

We report first-principles calculations on the coupling between epitaxial strain, polarization, and oxygen octahedra rotations in monodomain (PbTiO$_{3}$)$_{n}$/(SrTiO$_{3}$)$_{n}$ superlattices. We show how the interplay between (i) the epitaxial strain and (ii) the electrostatic conditions, can be used to control the orientation of the main axis of the system. The electrostatic constrains at the interface facilitate the rotation of the polarization and, as a consequence, we predict large piezoelectric responses at epitaxial strains smaller than those that would be required considering only strain effects. In addition, ferroelectric (FE) and antiferrodistortive (AFD) modes are strongly coupled. Usual steric arguments cannot explain this coupling and a covalent model is proposed to account for it. The energy gain due to the FE-AFD coupling decreases with the periodicity of the superlattice, becoming negligible for $n \ge 3$.

preprint2011arXiv

Nonlinear optical properties of TeO$_2$ crystalline phases from first principles

We have computed second and third nonlinear optical susceptibilities of two crystalline bulk tellurium oxide polymorphs: $α$-TeO$_{2}$ (the most stable crystalline bulk phase) and $γ$-TeO$_{2}$ (the crystalline phase that ressembles the more to the glass phase. Third order nonlinear susceptibilities of the crystalline phases are two orders of magnitude larger than $α$-SiO$_{2}$ cristoballite, thus extending the experimental observations on glasses to the case of crystalline compounds. While the electronic lone pairs of Te contribute to those large values, a full explanation of the anisotropy of the third order susceptibility tensor requires a detailed analysis of the structure, in particular the presence of helical chains, that seems to be linked to cooperative non-local polarizabilty effects. Our results demonstrate that first-principles simulations are a powerful predictive tool to estimate nonlinear optical susceptibilitites of materials.

preprint2010arXiv

Band selection and disentanglement using maximally-localized Wannier functions: the cases of Co impurities in bulk copper and the Cu (111) surface

We have adapted the maximally-localized Wannier function approach of [I. Souza, N. Marzari and D. Vanderbilt, Phys. Rev. B 65, 035109 (2002)] to the density functional theory based Siesta method [J. M. Soler et al., J. Phys.: Cond. Mat. 14, 2745 (2002)] and applied it to the study of Co substitutional impurities in bulk copper as well as to the Cu (111) surface. In the Co impurity case, we have reduced the problem to the Co d-electrons and the Cu sp-band, permitting us to obtain an Anderson-like Hamiltonian from well defined density functional parameters in a fully orthonormal basis set. In order to test the quality of the Wannier approach to surfaces, we have studied the electronic structure of the Cu (111) surface by again transforming the density functional problem into the Wannier representation. An excellent description of the Shockley surface state is attained, permitting us to be confident in the application of this method to future studies of magnetic adsorbates in the presence of an extended surface state.

preprint2010arXiv

First-principles study of competing ferroelectric and antiferroelectric instabilities in BaTiO3/BaO superlattices

We report a first-principles study of (BaTiO${_3}$)$_m$/(BaO)$_n$ superlattices for a wide range of periodicities $m/n$. We show that such a system develops a polar zone-center instability for sufficiently large \textit{m}/\textit{n} ratio, which can be understood, at least qualitatively, from a simple electrostatic model and should lead to a ferroelectric ground-state. However, the analysis of the phonon dispersion curves also points out the appearance of stronger antiferroelectric instabilities at the zone boundaries around $m=4$, before the critical ratio for ferroelectricity is reached and which still dominate beyond it. The dominant character of the anti-ferroelectric instability is explained from the depolarizing field which hardens the ferroelectric mode. This analysis allows us to predict that, (BaTiO${_3}$)$_m$/(BaO)$_n$ superlattices should present an antiferroelectric ground state for $m$ larger than 4, which should smoothly evolve to a multidomain structure for increasing $m$ values and only become ferroelectric for large $m$.