Researcher profile

Emilio Artacho

Emilio Artacho contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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