Source author record

Emilio Artacho

Emilio Artacho 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

36works
10topics
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

36 published item(s)

preprint2025arXiv

First-Principles Nanocapacitor Simulations of the Optical Dielectric Constant in Water Ice

We introduce a combined density functional theory (DFT) and non-equilibrium Green's function (NEGF) framework to compute the capacitance of nanocapacitors and directly extract the dielectric response of a sub-nanometer dielectric under bias. We identify that at the nanoscale conventional capacitance evaluations based on stored charge per unit voltage suffer from an ill-posed partitioning of electrode and dielectric charge. This partitioning directly impacts the geometric definition of capacitance through the capacitor width, which in turn makes the evaluation of dielectric response uncertain. This ambiguous separation further induces spurious interfacial polarizability when analyzed via maximally localized Wannier functions. Focusing on crystalline ice, we develop a robust charge-separation protocol that yields unique capacitance-derived polarizability and dielectric constants, unequivocally demonstrating that confinement neither alters ice's intrinsic electronic response nor its insensitivity to proton order. Our results lay the groundwork for rigorous interpretation of capacitor measurements in low-dimensional dielectric materials.

preprint2022arXiv

Bragg's Additivity Rule and Core and Bond model studied by real-time TDDFT electronic stopping simulations: the case of water vapor

The electronic stopping power ($S_e$) of water vapor (H$_2$O), hydrogen (H$_2$) and oxygen (O$_2$) gases for protons in a broad range of energies, centered in the Bragg peak, was calculated using real-time time-dependent density functional theory (rt-TDDFT) simulations with Gaussian basis sets. This was done for a kinetic energy of incident protons ($E_k$) ranging from 1.56 keV/amu to 1.6 MeV/amu. $S_e$ was calculated as the average over geometrically pre-sampled short ion trajectories. The average $S_e(E_k)$ values were found to rapidly converge with 25-30 pre-sampled, 2 nm-long ion trajectories. The rt-TDDFT $S_e(E_k)$ curves were compared to experimental and SRIM data, and used to validate the Bragg's Additivity Rule (BAR). Discrepancies were analyzed in terms of basis set effects and omitted nuclear stopping at low energies. At variance with SRIM, we found that BAR is applicable to our rt-TDDFT simulations of $\mathrm{2H_2+O_2}$ $\mathrm{\rightarrow 2H_2O}$ without scaling for $E_k>40$ keV/amu. The hydrogen and oxygen Core and Bond (CAB) contributions to electronic stopping were calculated and found to be slightly smaller than SRIM values as a result of a red-shift in our rt-TDDFT $S_e(E_k)$ curves and a re-distribution of weights due to some bond contributions being neglected in SRIM.

preprint2022arXiv

Inelastic scattering of electrons in water from first principles: cross sections and inelastic mean free path for use in Monte Carlo track-structure simulations of biological damage

Modelling the inelastic scattering of electrons in water is fundamental, given their crucial role in biological damage. In Monte Carlo track-structure codes used to assess biological damage, the energy loss function, from which cross sections are extracted, is derived from different semi-empirical optical models. Only recently, first ab-initio results for the energy loss function and cross-sections in water became available. For benchmarking purpose, in this work, we present ab-initio linear-response time-dependent density functional theory calculations of the energy loss function of liquid water. We calculated the inelastic scattering cross sections, inelastic mean free paths, and electronic stopping powers and compared our results with recent calculations and experimental data showing a good agreement. In addition, we provide an in-depth analysis of the contributions of different molecular orbitals, species, and orbital angular momenta to the total energy loss function. Moreover, we present single-differential cross sections computed for each molecular orbital channel, which should prove useful for Monte-Carlo track-structure simulations.

preprint2022arXiv

Semi-empirical and Linear-Scaling DFT Methods to Characterize duplex DNA and G-quadruplexes in Presence of Interacting Small Molecules

The computational study of DNA and its interaction with ligands is a highly relevant area of research, with significant consequences for developing new therapeutic strategies. However, the computational description of such large and complex systems requires considering interactions of different types. All these considerations imply a real challenge for computational chemistry. Using quantum methods for the entire system requires significant computational resources, with improvements in parallelization and optimization of theoretical strategies. Computational methods, such as LS-DFT and DLPNO-CCSD(T), may allow performing ab initio QM calculations, including explicitly the electronic structure for large biological systems, at a reasonable computing time. In this work, we study the interaction of small molecules and cations with DNA (duplex-DNA and G-quadruplexes), comparing different computational methods: a linear-scaling DFT (LS-DFT) at LMKLL/DZDP level of theory, semi-empirical methods (PM6-DH2 and PM7), mixed QM/MM, and DLPNO-CCSD(T). Our goal is to demonstrate the adequacy of LS-DFT to treat the different types of interactions present in DNA-dependent systems. We show that LMKLL/DZDP using SIESTA can yield very accurate geometries and energetics in all the different systems considered in this work: duplex DNA (dDNA), phenanthroline intercalating dDNA, G-quadruplexes, and Metal-G-tetrads considering alkaline metals of different sizes. As far as we know, this is the first time that full G-quadruplex geometry optimizations have been carried out using a DFT method thanks to its linear-scaling capabilities. Moreover, we show that LS-DFT provides high-quality structures, and some semi-empirical Hamiltonian can also yield suitable geometries. However, DLPNO-CCSD(T) and LS-DFT are the only methods that accurately describe interaction energies for all the systems considered in our study.

preprint2021arXiv

Local orbital formulation of the Floquet theory of projectile electronic stopping

A recently proposed theoretical framework for the description of electronic quantum friction for constant-velocity nuclear projectiles traversing periodic crystals is here implemented using a local basis representation. The theory requires a change of reference frame to that of the projectile, and a basis set transformation for the target basis functions to a gliding basis is presented, which is time-periodic but does not displace in space with respect to the projectile, allowing a local-basis Floquet impurity-scattering formalism to be used. It is illustrated for a one-dimensional single-band tight-binding model, as the simplest paradigmatic example, displaying the qualitative behaviour of the formalism. The time-dependent non-orthogonality of the gliding basis requires care in the proper (simplest) definition of a local projectile perturbation. The Fermi level is tilted with a slope given by the projectile velocity, which complicates integration over occupied states. It is solved by a recurrent application of the Lippmann-Schwinger equation, in analogy with previous non-equilibrium treatment of electron ballistic transport. Aiming towards a first-principles mean-field-like implementation, the final result is the time-periodic particle density in the region around the projectile, describing the stroboscopically stationary perturbation cloud around the projectile, out of which other quantities can be obtained, such as the electronic stopping power.

preprint2021arXiv

Manifold curvature and Ehrenfest forces with a moving basis

Known force terms arising in the Ehrenfest dynamics of quantum electrons and classical nuclei, due to a moving basis set for the former, can be understood in terms of the curvature of the manifold hosting the quantum states of the electronic subsystem. Namely, the velocity-dependent terms appearing in the Ehrenfest forces on the nuclei acquire a geometrical meaning in terms of the intrinsic curvature of the manifold, while Pulay terms relate to its extrinsic curvature.

preprint2020arXiv

Efficient ab initio calculation of electronic stopping in disordered systems via geometry pre-sampling: application to liquid water

Knowledge of the electronic stopping curve for swift ions, $S_e(v)$, particularly around the Bragg peak, is important for understanding radiation damage. Experimentally, however, the determination of such feature for light ions is very challenging, especially in disordered systems such as liquid water and biological tissue. Recent developments in real-time time-dependent density functional theory (rt-TDDFT) have enabled the calculation of $S_e(v)$ along nm-sized trajectories. However, it is still a challenge to obtain a meaningful statistically averaged $S_e(v)$ that can be compared to observations. In this work, taking advantage of the correlation between the local electronic structure probed by the projectile and the distance from the projectile to the atoms in the target, we devise a trajectory pre-sampling scheme to select, geometrically, a small set of short trajectories to accelerate the convergence of the averaged $S_e(v)$ computed via rt-TDDFT. For protons in liquid water, we first calculate the reference probability distribution function (PDF) for the distance from the proton to the closest oxygen atom, $ϕ_R(r_{p{\rightarrow}O})$, for a trajectory of a length similar to those sampled experimentally. Then, short trajectories are sequentially selected so that the accumulated PDF reproduces $ϕ_R(r_{p{\rightarrow}O})$ to increasingly high accuracy. Using these pre-sampled trajectories, we demonstrate that the averaged $S_e(v_p)$ converges in the whole velocity range with less than eight trajectories, while other averaging methods using randomly and uniformly distributed trajectories require approximately ten times the computational effort. This allows us to compare the $S_e(v_p)$ curve to experimental data, and assess widely used empirical tables based on Bragg's rule.

preprint2020arXiv

Entropic bonding of the type 1 pilus from experiment and simulation

The type 1 pilus is a bacterial filament consisting of a long coiled proteic chain of subunits joined together by non-covalent bonding between complementing $β$-strands. Its strength and structural stability are critical for its anchoring function in uropathogenic Escherichia coli bacteria. The pulling and unravelling of the FimG subunit of the pilus was recently studied by atomic force microscopy (AFM) experiments and steered molecular dynamics (SMD) simulations [A. Alonso-Caballero et al., Nature Commun. 9, 2758 (2018)]. In this work we perform a quantitative comparison between experiment and simulation, showing a good agreement in the underlying work values for the unfolding. The simulation results are then used to estimate the free energy difference for the detachment of FimG from the complementing strand of the neighbouring subunit in the chain, FimF. Finally, we show that the large free energy difference for the unravelling and detachment of the subunits which leads to the high stability of the chain is entirely entropic in nature.

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.

preprint2020arXiv

The CECAM Electronic Structure Library and the modular software development paradigm

First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (\esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include "heavy-duty" ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers.

preprint2019arXiv

Dynamical continuum simulation of condensed matter from first-principles

Macroscale continuum mechanics simulations rely on material properties stemming from the microscale, which are normally described using phenomenological equations of state (EOS). A method is proposed for the automatic generation of first-principles unconstrained EOSs using a Gaussian process on a set of ab initio molecular dynamics simulations, thereby closing the continuum equations. We illustrate it on a hyperelasticity simulation of bulk silicon using density-functional theory (DFT), following the dynamics of shock waves after a cylindrical region is instantaneously heated.

preprint2019arXiv

Electrostatics and domains in ferroelectric superlattices

The electrostatics arising in ferroelectric/dielectric two-dimensional heterostructures and superlatitices is revisited here within a simplest Kittel model, in order to define a clear paradigmatic reference for domain formation. The screening of the depolarizing field in isolated ferroelectric or polar thin films via the formation of 180$^{\circ}$ domains is well understood, whereby the width of the domains $w$ grows as the square-root of the film thickness $d$, following Kittel's law, for thick enough films ($w\ll d$). This behavior is qualitatively unaltered when the film is deposited on a dielectric substrate, sandwiched between dielectrics, and even in a superlattice setting, with just a suitable renormalisation of Kittel's length. As $d$ decreases, $w(d)$ deviates from Kittel's law, reaching a minimum and then diverging onto the mono-domain limit for thin enough films, always assuming a given spontaneous polarization $P$ of the ferrolectric, only modified by linear response to the depolarizing field. In most cases of experimental relevance $P$ would vanish before reaching that thin-film regime. This is not the case for superlattices. Unlike single films, for which the increase of the dielectric constant of the surrounding medium pushes the deviation from the Kittel's regime to lower values of $d$, there is a critical value of the relative thickness of ferroelectric/dielectric films in superlattices beyond which that behavior is reversed, and which defines the separation between strong and weak ferroelectric coupling in superlattices.

preprint2019arXiv

Floquet theory for the electronic stopping of projectiles in solids

A general theoretical framework for the study of electronic stopping of particle projectiles in crystalline solids is proposed. It neither relies on perturbative or linear response approximations, nor on an ideal metal host. Instead, it exploits the discrete translational symmetries in a space-time diagonal determined by a projectile with constant velocity moving along a trajectory with crystalline periodicity. This allows for the characterisation of (stroboscopically) stationary solutions, by means of Floquet theory for time-periodic systems. Previous perturbative and non-linear jellium models are recovered from this general theory. An analysis of the threshold velocity effect in insulators is presented based on Floquet quasi-energy conservation.

preprint2016arXiv

Continuous melting through a hexatic phase in confined bilayer water

Liquid water is not only of obvious importance but also extremely intriguing, displaying many anomalies that still challenge our understanding of such an a priori simple system. The same is true when looking at nanoconfined water: The liquid between constituents in a cell is confined to such dimensions, and there is already evidence that such water can behave very differently from its bulk counterpart. A striking finding has been reported from computer simulations for two-dimensionally confined water: The liquid displays continuous or discontinuous melting depending on its density. In order to understand this behavior, we have analyzed the melting exhibited by a bilayer of nanoconfined water by means of molecular dynamics simulations. At high density we observe the continuous melting to be related to the phase change of the oxygens only, with the hydrogens remaining liquid-like throughout. Moreover, we find an intermediate hexatic phase for the oxygens between the liquid and a triangular solid ice phase, following the Kosterlitz-Thouless-Halperin-Nelson-Young theory for two-dimensional melting. The liquid itself tends to maintain the local structure of the triangular ice, with its two layers being strongly correlated, yet with very slow exchange of matter. The decoupling in the behavior of the oxygens and hydrogens gives rise to a regime in which the complexity of water seems to disappear, resulting in what resembles a simple monoatomic liquid. This intrinsic tendency of our simulated water may be useful for understanding novel behaviors in other confined and interfacial water systems.

preprint2016arXiv

Enhanced configurational entropy in high-density nanoconfined bilayer ice

A novel kind of crystal order in high-density nanoconfined bilayer ice is proposed from molecular dynamics and density-functional theory simulations. A first-order transition is observed between a low-temperature proton-ordered solid and a high-temperature proton-disordered solid. The latter is shown to possess crystalline order for the oxygen positions, arranged on a close-packed triangular lattice with AA stacking. Uniquely amongst the ice phases, the triangular bilayer is characterized by two levels of disorder (for the bonding network and for the protons) which results in a configurational entropy twice that of bulk ice.

preprint2016arXiv

Structural and configurational properties of nanoconfined monolayer ice from first principles

Understanding the structural tendencies of nanoconfined water is of great interest for nanoscience and biology, where nano/micro-sized objects may be separated by very few layers of water. Here we investigate the properties of ice confined to a quasi-2D monolayer by a featureless, chemically neutral potential, in order to characterize its intrinsic behaviour. We use density-functional theory simulations with a non-local van der Waals density functional. An ab initio random structure search reveals all the energetically competitive monolayer configurations to belong to only two of the previously-identified families, characterized by a square or honeycomb hydrogen-bonding network, respectively. We discuss the modified ice rules needed for each network, and propose a simple point dipole 2D lattice model that successfully explains the energetics of the square configurations. All identified stable phases for both networks are found to be non-polar (but with a topologically non-trivial texture for the square) and, hence, non-ferroelectric, in contrast to previous predictions from a five-site empirical force-field model. Our results are in good agreement with very recently reported experimental observations.

preprint2015arXiv

Ab initio calculation of the shock Hugoniot of bulk silicon

We describe a simple annealing procedure to obtain the Hugoniot locus (states accessible by a shock wave) for a given material in a computationally efficient manner. We apply this method to determine the Hugoniot locus in bulk silicon from ab initio molecular dynamics with forces from density-functional theory, up to 70 GPa. The fact that shock waves can split into multiple waves due to phase transitions or yielding is taken into account here by specifying the strength of any preceding waves explicitly based on their yield strain. Points corresponding to uniaxial elastic compression along three crystal axes and a number of post-shock phases are given, including a plastically-yielded state, approximated by an isotropic stress configuration following an elastic wave of predetermined strength. The results compare well to existing experimental data for shocked silicon.

preprint2015arXiv

Electronic stopping power in a narrow band gap semiconductor from first principles

The direction and impact parameter dependence of electronic stopping power, along with its velocity threshold behavior, is investigated in a prototypical small band gap semiconductor. We calculate the electronic stopping power of H in Ge, a semiconductor with relatively low packing density, using time-evolving time-dependent density-functional theory. The calculations are carried out in channeling conditions with different impact parameters and in different crystal directions, for projectile velocities ranging from 0.05 to 0.6 atomic units. The satisfactory comparison with available experiments supports the results and conclusions beyond experimental reach. The calculated electronic stopping power is found to be different in different crystal directions; however, strong impact parameter dependence is observed only in one of these directions. The distinct velocity threshold observed in experiments is well reproduced, and its non-trivial relation with the band gap follows a perturbation theory argument surprisingly well. This simple model is also successful in explaining why different density functionals give the same threshold even with substantially different band gaps.

preprint2015arXiv

Model of two-dimensional electron gas formation at ferroelectric interfaces

The formation of a two-dimensional electron gas at oxide interfaces as a consequence of polar discontinuities has generated an enormous amount of activity due to the variety of interesting effects it gives rise to. Here we study under what circumstances similar processes can also take place underneath ferroelectric thin films. We use a simple Landau model to demonstrate that in the absence of extrinsic screening mechanisms a monodomain phase can be stabilized in ferroelectric films by means of an electronic reconstruction. Unlike in the LaAlO$_3$/SrTiO$_3$ heterostructure, the emergence with thickness of the free charge at the interface is discontinuous. This prediction is confirmed by performing first principles simulations of free standing slabs of PbTiO$_3$. The model is also used to predict the response of the system to an applied electric field, demonstrating that the two-dimensional electron gas can be switched on and off discontinuously and in a non-volatile fashion. Furthermore, the reversal of the polarization can be used to switch between a two-dimensional electron gas and a two-dimensional hole gas, which should, in principle, have very different transport properties. We discuss the possible formation of polarization domains and how such configuration competes with the spontaneous accumulation of free charge at the interfaces.

preprint2015arXiv

Two-dimensional electron gas at the PbTiO$_3$/SrTiO$_3$ interface: an ab initio study

In the polar catastrophe scenario, polar discontinuity accounts for the driving force of the formation of a two-dimensional electron gas (2DEG) at the interface between polar and non-polar insulators. In this paper, we substitute the usual, non-ferroelectric, polar material with a ferroelectric thin film and use the ferroelectric polarization as the source for polar discontinuity. We use ab initio simulations to systematically investigate the stability, formation and properties of the two-dimensional free-carrier gases formed in PbTiO$_3$/SrTiO$_3$ heterostructures under realistic mechanical and electrical boundary conditions. Above a critical thickness, the ferroelectric layers can be stabilized in the out-of-plane monodomain configuration due to the electrostatic screening provided by the free-carriers. Our simulations also predict that the system can be switched between three stable configurations (polarization up, down or zero), allowing the non-volatile manipulation of the free charge density and sign at the interface. Furthermore, the link between ferroelectric polarization and free charge density demonstrated by our analysis constitutes compelling support for the polar catastrophe model that is used to rationalize the formation of 2DEG at oxide interfaces.

preprint2014arXiv

Effects of stoichiometric doping in superconducting Bi-O-S compounds

Newly discovered Bi-O-S compounds remain an enigma in attempts to understand their electronic properties. A recent study of Bi$_{4}$O$_{4}$S$_{3}$ has shown it to be a mixture of two phases, Bi$_{2}$OS$_{2}$ and Bi$_{3}$O$_{2}$S$_{3}$, the latter being superconducting [W. A. Phelan et al., J. Am. Chem. Soc. 135, 5372 (2013)]. Using density functional theory, we explore the electronic structure of both the phases and the effect of the introduction of stacking faults. Our results demonstrate that the S$_{2}$ layers dope the bismuth-sulphur bands. The bands at the Fermi level are of clear two-dimensional character. One band manifold is confined to the two adjacent, square-lattice bismuth-sulphur planes, a second manifold is confined to the square lattice of sulphur dimers. We show that the introduction of defects in the stacking does not influence the electronic structure. Finally, we also show that spin-orbit coupling does not have any significant effect on the states close to the Fermi level at the energy scale considered.

preprint2014arXiv

Origin of two-dimensional electron gases at oxide interfaces: insights from theory

The response of oxide thin films to polar discontinuities at interfaces and surfaces has generated an enormous activity due to the variety of interesting effects it gives rise to. A case in point is the discovery of the electron gas at the interface between LaAlO3 and SrTiO3, which has since been shown to be quasi-two-dimensional, switchable, magnetic and/or superconducting. Despite these findings, the origin of the two-dimensional electron gas is highly debated and several possible mechanisms remain. Here we review the main proposed mechanisms and attempt to model expected effects in a quantitative way with the ambition of better constraining what effects can/cannot explain the observed phenomenology. We do it in the framework of a phenomenological model for understanding electronic and/or redox screening of the chemical charge in oxide heterostructures. We also discuss the effect of intermixing, both conserving and non-conserving the total stoichiometry.

preprint2013arXiv

Knock-on damage in bilayer graphene: indications for a catalytic pathway

We study by high-resolution transmission electron microscopy the structural response of bilayer graphene to electron irradiation with energies below the knock-on damage threshold of graphene. We observe that one type of divacancy, which we refer to as the butterfly defect, is formed for radiation energies and doses for which no vacancies are formed in clean monolayer graphene. By using first principles calculations based on density-functional theory, we analyze two possible causes related with the presence of a second layer that could explain the observed phenomenon: an increase of the defect stability or a catalytic effect during its creation. For the former, the obtained formation energies of the defect in monolayer and bilayer systems show that the change in stability is negligible. For the latter, ab initio molecular dynamics simulations indicate that the threshold energy for direct expulsion does not decrease in bilayer graphene as compared with monolayer graphene, and we demonstrate the possibility of creating divacancies through catalyzed intermediate states below this threshold energy. The estimated cross section agrees with what is observed experimentally. Therefore, we show the possibility of a catalytic pathway for creating vacancies under electron radiation below the expulsion threshold energy.

preprint2013arXiv

Optimal finite-range atomic basis sets for liquid water and ice

Finite-range numerical atomic orbitals are the basis functions of choice for several first principles methods, due to their flexibility and scalability. Generating and testing such basis sets, however, remains a significant challenge for the end user. We discuss these issues and present a new scheme for generating improved polarization orbitals of finite range. We then develop a series of high-accuracy basis sets for the water molecule, and report on their performance in describing the monomer and dimer, two phases of ice, and liquid water at ambient and high density. The tests are performed by comparison with plane-wave calculations, and show the atomic orbital basis sets to exhibit an excellent level of transferability and consistency. The highest-order bases (quadruple-zeta) are shown to give accuracies comparable to a plane-wave kinetic energy cutoff of at least ~1000 eV for quantities such as energy differences and ionic forces, as well as achieving significantly greater accuracies for total energies and absolute pressures.

preprint2013arXiv

Room temperature compressibility and diffusivity of liquid water from first principles

The isothermal compressibility of water is essential to understand its anomalous properties. We compute it by ab initio molecular dynamics simulations of 200 molecules at five densities, using two different van der Waals density functionals. While both functionals predict compressibilities within ~30% of experiment, only one of them accurately reproduces, within the uncertainty of the simulation, the density dependence of the self-diffusion coefficient in the anomalous region. The discrepancies between the two functionals are explained in terms of the low- and high-density structures of the liquid.

preprint2013arXiv

Species Fractionation in Atomic Chains from Mechanically Stretched Alloys

Bettini et al. [Nature Nanotech 1, 182 (2006)] reported the first experimental realization of linear atomic chains (LACs) composed of different atoms (Au and Ag). Different contents of Au and Ag were observed in the chains from what found in the bulk alloys, which rises the question of what is the wire composition if in equilibrium with a bulk alloy. In this work we address the thermodynamic driving force for species fractionation in LACs under tension, and we present density-functional theory results for Ag-Au chain alloys. A pronounced stabilization of wires with an alternating Ag-Au sequence is observed, which could be behind the experimentally observed Au enrichment in LACs from alloys of high Ag content.

preprint2012arXiv

Electrochemical ferroelectric switching: The origin of polarization reversal in ultrathin films

Against expectations, robust switchable ferroelectricity has been recently observed in ultrathin (1 nm) ferroelectric films exposed to air [V. Garcia $et$ $al.$, Nature {\bf 460}, 81 (2009)]. Based on first-principles calculations, we show that the system does not polarize unless charged defects or adsorbates form at the surface. We propose electrochemical processes as the most likely origin of this charge. The ferroelectric polarization of the film adapts to the bound charge generated on its surface by redox processes when poling the film. This, in turn, alters the band alignment at the bottom electrode interface, explaining the observed tunneling electroresistance. Our conclusions are supported by energetics calculated for varied electrochemical scenarios.

preprint2012arXiv

Electronic stopping power in gold: The role of d electrons and the H/He anomaly

The electronic stopping power of H and He moving through gold is obtained to high accuracy using time-evolving density-functional theory, thereby bringing usual first-principles accuracies into this kind of strongly coupled, continuum non-adiabatic processes in condensed matter. The two key unexplained features of what observed experimentally have been reproduced and understood: (i) The non-linear behaviour of stopping power versus velocity is a gradual crossover as excitations tail into the d-electron spectrum; and (ii) the low-velocity H/He anomaly (the relative stopping powers are contrary to established theory) is explained by the substantial involvement of the d electrons in the screening of the projectile even at the lowest velocities where the energy loss is generated by s-like electron-hole pair formation only.

preprint2012arXiv

One-dimensional half-metallic interfaces of two-dimensional honeycomb insulators

We study zigzag interfaces between insulating compounds that are isostructural to graphene, specifically II-VI, III-V and IV-IV two-dimensional (2D) honeycomb insulators. We show that these one-dimensional interfaces are polar, with a net density of excess charge that can be simply determined by using the ideal (integer) formal valence charges, regardless of the predominant covalent character of the bonding in these materials. We justify this finding on fundamental physical grounds, by analyzing the topology of the formal polarization lattice in the parent bulk materials. First principles calculations elucidate an electronic compensation mechanism not dissimilar to oxide interfaces, which is triggered by a Zener-like charge transfer between interfaces of opposite polarity. In particular, we predict the emergence of one dimensional electron and hole gases (1DEG), which in some cases are ferromagnetic half-metallic.

preprint2012arXiv

Van der Waals interaction in magnetic bilayer graphene nanoribbons

We study the interaction energy between two graphene nanoribbons by first principles calculations, including van der Waals interactions and spin polarization. For ultranarrow zigzag nanoribbons, the direct stacking is even more stable than Bernal, competing in energy for wider ribbons. This behavior is due to the magnetic interaction between edge states. We relate the reduction of the magnetization in zigzag nanoribbons with increasing ribbon width to the structural changes produced by the magnetic interaction, and show that when deposited on a substrate, zigzag bilayer ribbons remain magnetic for larger widths.

preprint2011arXiv

Proposal of a one-dimensional electron gas in the steps at the LaAlO$_3$-SrTiO$_3$ interface

The two-dimensional electron gas (2DEG) at the interface between LaAlO$_3$ (LAO) and SrTiO$_3$ (STO) has become one of the most fascinating and highly-debated oxide systems of recent times. Here we propose that a one-dimensional electron gas (1DEG) can be engineered at the step edges of the LAO/STO interface. These predictions are supported by first principles calculations and electrostatic modeling which elucidate the origin of the 1DEG as an electronic reconstruction to compensate a net surface charge in the step edge. The results suggest a novel route to increasing the functional density in these electronic interfaces.

preprint2011arXiv

Surface defects and conduction in polar oxide heterostructures

The polar interface between LaAlO$_{3}$ and SrTiO$_{3}$ has shown promise as a field effect transistor, with reduced (nanoscale) feature sizes and potentially added functionality over conventional semiconductor systems. However, the mobility of the interfacial two-dimensional electron gas (2DEG) is lower than desirable. Therefore to progress, the highly debated origin of the 2DEG must be understood. Here we present a case for surface redox reactions as the origin of the 2DEG, in particular surface O vacancies, using a model supported by first principles calculations that describes the redox formation. In agreement with recent spectroscopic and transport measurements, we predict a stabilization of such redox processes (and hence Ti 3$d$ occupation) with film thickness beyond a critical value, which can be smaller than the critical thickness for 2D electronic conduction, since the surface defects generate trapping potentials that will affect the interface electron mobility. Several other recent experimental results, such as lack of core level broadening and shifts, find natural explanation. Pristine systems will likely require changed growth conditions or modified materials with a higher vacancy free energy.

preprint2010arXiv

Density, structure and dynamics of water: the effect of Van der Waals interactions

It is known that ab initio molecular dynamics (AIMD) simulations of liquid water, based on the generalized gradient approximation (GGA) to density functional theory (DFT), yield structural and diffusive properties in reasonable agreement with experiment only if artificially high temperatures are used in the simulations. The equilibrium density, at normal conditions, of DFT water has been recently shown by Schmidt et al. [J. Phys. chem. B, 113, 11959 (2009)] to be underestimated by different GGA functionals for exchange and correlation, and corrected by the addition of interatomic pair potentials to describe van derWaals (vdW) interactions. In this contribution we present a DFTAIMD study of liquid water using several GGA functionals as well as the van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401(2004)]. As expected, we find that the density of water is grossly underestimated by GGA functionals. When a vdW-DF is used, the density improves drastically and the experimental diffusivity is reproduced without the need of thermal corrections. We analyze the origin of the density differences between all the functionals. We show that the vdW-DF increases the population of non-H-bonded interstitial sites, at distances between the first and second coordination shells. However, it excessively weakens the H-bond network, collapsing the second coordination shell. This structural problem is partially associated to the choice of GGA exchange in the vdW-DF. We show that a different choice for the exchange functional is enough to achieve an overall improvement both in structure and diffusivity.

preprint2009arXiv

Topology of the polarization field in ferroelectric nanowires from first principles

The behaviour of the cross-sectional polarization field is explored for thin nanowires of barium titanate from first-principles calculations. Topological defects of different winding numbers have been obtained, beyond the known textures in ferroelectric nanostructures. They result from the inward accommodation of the polarization patterns imposed at the surface of the wire by surface and edge effects. Close to a topological defect the polarization field orients out of the basal plane in some cases, maintaining a close to constant magnitude, whereas it virtually vanishes in other cases.

preprint2006arXiv

Resistive and rectifying effects of pulling gold atoms at thiol-gold nano-contacts

We investigate, by means of first-principles calculations, structural and transport properties of junctions made of symmetric dithiolated molecules placed between Au electrodes. As the electrodes are pulled apart, we find that it becomes energetically favorable that Au atoms migrate to positions between the electrode surface and thiol terminations, with junction structures alternating between symmetric and asymmetric. As a result, the calculated $\emph{IV}$ curves alternate between rectifying and non-rectifying behaviors as the electrodes are pulled apart, which is consistent with recent experimental results.

preprint1999arXiv

Linear-scaling ab-initio calculations for large and complex systems

A brief review of the SIESTA project is presented in the context of linear-scaling density-functional methods for electronic-structure calculations and molecular-dynamics simulations of systems with a large number of atoms. Applications of the method to different systems are reviewed, including carbon nanotubes, gold nanostructures, adsorbates on silicon surfaces, and nucleic acids. Also, progress in atomic-orbital bases adapted to linear-scaling methodology is presented.