Source author record

Pablo García-Fernández

Pablo García-Fernández appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.str-el physics.comp-ph

Catalog footprint

What is connected

12works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

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.

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.

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

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

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$.

Pablo García-Fernández

What is connected

Connect this record

See the researcher in context

Building this map preview

12 published item(s)

Emergent chirality in a polar meron to skyrmion phase transition

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

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

Order-disorder transitions in a polar vortex lattice

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

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

SIESTA: recent developments and applications

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

Second-principles method including electron and lattice degrees of freedom

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

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

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