Showing 1–22 of 22 results for author: Andrade, X

Dynamics of ballistic photocurrents driven by Coulomb scattering

Authors: Liang Z. Tan, Xavier Andrade, Sangeeta Rajpurohit, Alfredo A. Correa, Tadashi Ogitsu

Abstract: First principles real-time time dependent density functional theory (rt-TDDFT) calculations reveal the existence of ballistic photocurrents generated by Coulomb scattering, a form of photocurrent that has not previously been considered as a mechanism for the bulk photovoltaic effect. With monolayer GeS as an example, it is predicted that ballistic currents can exceed shift currents under experimen… ▽ More First principles real-time time dependent density functional theory (rt-TDDFT) calculations reveal the existence of ballistic photocurrents generated by Coulomb scattering, a form of photocurrent that has not previously been considered as a mechanism for the bulk photovoltaic effect. With monolayer GeS as an example, it is predicted that ballistic currents can exceed shift currents under experimentally accessible conditions. Moreover, these simulations reveal recombination pathways that work to modify shift current magnitudes beyond perturbation theory predictions. △ Less

Submitted 12 August, 2024; originally announced August 2024.

Multi-Objective Optimization of Consumer Group Autoscaling in Message Broker Systems

Authors: Diogo Landau, Nishant Saurabh, Xavier Andrade, Jorge G Barbosa

Abstract: Message brokers often mediate communication between data producers and consumers by adding variable-sized messages to ordered distributed queues. Our goal is to determine the number of consumers and consumer-partition assignments needed to ensure that the rate of data consumption keeps up with the rate of data production. We model the problem as a variable item size bin packing problem. As the rat… ▽ More Message brokers often mediate communication between data producers and consumers by adding variable-sized messages to ordered distributed queues. Our goal is to determine the number of consumers and consumer-partition assignments needed to ensure that the rate of data consumption keeps up with the rate of data production. We model the problem as a variable item size bin packing problem. As the rate of production varies, new consumer-partition assignments are computed, which may require rebalancing a partition from one consumer to another. While rebalancing a queue, the data being produced into the queue is not read leading to additional latency costs. As such, we focus on the multi-objective optimization cost of minimizing both the number of consumers and queue migrations. We present a variety of algorithms and compare them to established bin packing heuristics for this application. Comparing our proposed consumer group assignment strategy with Kafka's, a commonly employed strategy, our strategy presents a 90th percentile latency of 4.52s compared to Kafka's 217s with both using the same amount of consumers. Kafka's assignment strategy only improved the consumer group's performance with regards to latency with configurations that used at least 60% more resources than our approach. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 17 pages. arXiv admin note: substantial text overlap with arXiv:2206.11170

Kafka Consumer Group Autoscaler

Authors: Diogo Landau, Xavier Andrade, Jorge G. Barbosa

Abstract: Message brokers enable asynchronous communication between data producers and consumers in distributed environments by assigning messages to ordered queues. Message broker systems often provide with mechanisms to parallelize tasks between consumers to increase the rate at which data is consumed. The consumption rate must exceed the production rate or queues would grow indefinitely. Still, consumers… ▽ More Message brokers enable asynchronous communication between data producers and consumers in distributed environments by assigning messages to ordered queues. Message broker systems often provide with mechanisms to parallelize tasks between consumers to increase the rate at which data is consumed. The consumption rate must exceed the production rate or queues would grow indefinitely. Still, consumers are costly and their number should be minimized. We model the problem of determining the required number of consumers, and the partition-consumer assignments, as a variable item size bin packing variant. Data cannot be read when a queue is being migrated to another consumer. Hence, we propose the R-score metric to account for these rebalancing costs. Then, we introduce an assortment of R-score based algorithms, and compare their performance to established heuristics for the Bin Packing Problem for this application. We instantiate our method within an existing system, demonstrating its effectiveness. Our approach guarantees adequate consumption rates something the previous system was unable to at lower operational costs. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 10 pages, 10 figures, MSc dissertation work

ACM Class: D.4.1

Electron dynamics in extended systems within real-time time-dependent density functional theory

Authors: Alina Kononov, Cheng-Wei Lee, Tatiane Pereira dos Santos, Brian Robinson, Yifan Yao, Yi Yao, Xavier Andrade, Andrew David Baczewski, Emil Constantinescu, Alfredo A. Correa, Yosuke Kanai, Normand Modine, Andre Schleife

Abstract: Due to a beneficial balance of computational cost and accuracy, real-time time-dependent density functional theory has emerged as a promising first-principles framework to describe electron real-time dynamics. Here we discuss recent implementations around this approach, in particular in the context of complex, extended systems. Results include an analysis of the computational cost associated with… ▽ More Due to a beneficial balance of computational cost and accuracy, real-time time-dependent density functional theory has emerged as a promising first-principles framework to describe electron real-time dynamics. Here we discuss recent implementations around this approach, in particular in the context of complex, extended systems. Results include an analysis of the computational cost associated with numerical propagation and when using absorbing boundary conditions. We extensively explore the shortcomings for describing electron-electron scattering in real time and compare to many-body perturbation theory. Modern improvements of the description of exchange and correlation are reviewed. In this work, we specifically focus on the Qb@ll code, which we have mainly used for these types of simulations over the last years, and we conclude by pointing to further progress needed going forward. △ Less

Submitted 9 May, 2022; originally announced May 2022.

INQ, a modern GPU-accelerated computational framework for (time-dependent) density functional theory

Authors: Xavier Andrade, Chaitanya Das Pemmaraju, Alexey Kartsev, Jun Xiao, Aaron Lindenberg, Sangeeta Rajpurohit, Liang Z. Tan, Tadashi Ogitsu, Alfredo A. Correa

Abstract: We present INQ, a new implementation of density functional theory (DFT) and time-dependent DFT (TDDFT) written from scratch to work on graphical processing units (GPUs). Besides GPU support, INQ makes use of modern code design features and takes advantage of newly available hardware. By designing the code around algorithms, rather than against specific implementations and numerical libraries, we a… ▽ More We present INQ, a new implementation of density functional theory (DFT) and time-dependent DFT (TDDFT) written from scratch to work on graphical processing units (GPUs). Besides GPU support, INQ makes use of modern code design features and takes advantage of newly available hardware. By designing the code around algorithms, rather than against specific implementations and numerical libraries, we aim to provide a concise and modular code. The result is a fairly complete DFT/TDDFT implementation in roughly 12,000 lines of open-source C++ code representing a modular platform for community-driven application development on emerging high-performance computing architectures for the simulation of materials. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: 25 pages, 10 figures

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Authors: Nicolas Tancogne-Dejean, Micael J. T. Oliveira, Xavier Andrade, Heiko Appel, Carlos H. Borca, Guillaume Le Breton, Florian Buchholz, Alberto Castro, Stefano Corni, Alfredo A. Correa, Umberto De Giovannini, Alain Delgado, Florian G. Eich, Johannes Flick, Gabriel Gil, Adrián Gomez, Nicole Helbig, Hannes Hübener, René Jestädt, Joaquim Jornet-Somoza, Ask H. Larsen, Irina V. Lebedeva, Martin Lüders, Miguel A. L. Marques, Sebastian T. Ohlmann , et al. (9 additional authors not shown)

Abstract: Over the last years extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high-degree of precision. An appealing and challenging route towards engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a… ▽ More Over the last years extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high-degree of precision. An appealing and challenging route towards engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, providing an unique framework allowing to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory framework. The present article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, like the new theoretical framework of quantum electrodynamics density-functional formalism (QEDFT) for the description of novel light-matter hybrid states. Those advances, and other being released soon as part of the Octopus package, will enable the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (QED-materials). △ Less

Submitted 17 December, 2019; originally announced December 2019.

Journal ref: J. Chem. Phys. 152, 124119 (2020)

Effect of chemical disorder on the electronic stopping of solid solution alloys

Authors: Edwin E. Quashie, Rafi Ullah, Xavier Andrade, Alfredo A. Correa

Abstract: The electronic stopping power of nickel-based equiatomic solid solutions alloys NiCr, NiFe and NiCo for protons and alpha projectiles is investigated in detail using real-time time-dependent density functional theory over a wide range of velocities. Recently developed numerical electronic structure methods are used to probe fundamental aspects of electron-ion coupling non-perturbatively and in a f… ▽ More The electronic stopping power of nickel-based equiatomic solid solutions alloys NiCr, NiFe and NiCo for protons and alpha projectiles is investigated in detail using real-time time-dependent density functional theory over a wide range of velocities. Recently developed numerical electronic structure methods are used to probe fundamental aspects of electron-ion coupling non-perturbatively and in a fully atomistic context, capturing the effect of the atomic scale disorder. The effects of particular electronic band structures and density of states reflect in the low velocity limit behavior. We compare our results for the alloys with those of a pure nickel target to understand how alloying affects the electronic stopping. We discover that NiCo and NiFe have similar stopping behavior as Ni while NiCr has an asymptotic stopping power that is more than a factor of two larger than its counterparts for velocities below 0.1 a.u.. We show that the low-velocity limit of electronic stopping power can be manipulated by controlling the broadening of the d-band through the chemical disorder. In this regime, the Bragg's additive rule for the stopping of composite materials also fails for NiCr. △ Less

Submitted 16 July, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: 9 pages, 7 figures

Journal ref: Acta Mater., 196, 575 (2020)

Non-linear conductivity of metals from real-time quantum simulations

Authors: Xavier Andrade, Sebastien Hamel, Alfredo A. Correa

Abstract: We simulate bulk materials under strong currents by following in real-time the dynamics of the electrons under an electric field. By changing the intensity of the electric field, our method can model, for the first time, non-linear effects in the conductivity from first principles. To illustrate our approach, we show calculations that predict that liquid aluminum exhibits negative-differential con… ▽ More We simulate bulk materials under strong currents by following in real-time the dynamics of the electrons under an electric field. By changing the intensity of the electric field, our method can model, for the first time, non-linear effects in the conductivity from first principles. To illustrate our approach, we show calculations that predict that liquid aluminum exhibits negative-differential conductivity for current densities of the order of $10^{12}-10^{13}~\mathrm{A/cm^2}$. We find that the change in the non-linear conductivity emerges from a competition between the accumulation of charge around the nuclei that increases the scattering of the conduction electrons, and a decreasing scattering cross-section at high currents. △ Less

Submitted 1 February, 2017; originally announced February 2017.

Comments: 8 pages, 8 figures

Self-interaction effects on charge-transfer collisions

Authors: Edwin E. Quashie, Bidhan C. Saha, Xavier Andrade, Alfredo A. Correa

Abstract: In this article, we investigate the role of the self-interaction error in the simulation of collisions using time-dependent density functional theory (TDDFT) and Ehrenfest dynamics. We compare many different approximations of the exchange and correlation potential, using as a test system the collision of $\mathrm{H^+ + CH_4}$ at $30~\mathrm{eV}$. We find that semi-local approximations, like PBE, a… ▽ More In this article, we investigate the role of the self-interaction error in the simulation of collisions using time-dependent density functional theory (TDDFT) and Ehrenfest dynamics. We compare many different approximations of the exchange and correlation potential, using as a test system the collision of $\mathrm{H^+ + CH_4}$ at $30~\mathrm{eV}$. We find that semi-local approximations, like PBE, and even hybrid functionals, like B3LYP, produce qualitatively incorrect predictions for the scattering of the proton. This discrepancy appears because the self-interaction error allows the electrons to jump too easily to the proton, leading to radically different forces with respect to the non-self-interacting case. From our results, we conclude that using a functional that is self-interaction free is essential to properly describe charge-transfer collisions between ions and molecules in TDDFT. △ Less

Submitted 5 April, 2017; v1 submitted 8 September, 2016; originally announced September 2016.

Comments: 9 pages, 8 figures

Journal ref: Phys. Rev. A 95, 042517 (2017)

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Authors: Xavier Andrade, David A. Strubbe, Umberto De Giovannini, Ask Hjorth Larsen, Micael J. T. Oliveira, Joseba Alberdi-Rodriguez, Alejandro Varas, Iris Theophilou, Nicole Helbig, Matthieu Verstraete, Lorenzo Stella, Fernando Nogueira, Alán Aspuru-Guzik, Alberto Castro, Miguel A. L. Marques, Ángel Rubio

Abstract: Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of… ▽ More Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems. △ Less

Submitted 22 January, 2015; originally announced January 2015.

Comments: 28 pages, 16 figures

Journal ref: Phys. Chem. Chem. Phys. 17, 31371-31396 (2015)

A sparse-sampling approach for the fast computation of matrices: application to molecular vibrations

Authors: Jacob N. Sanders, Xavier Andrade, Alán Aspuru-Guzik

Abstract: This article presents a new method to compute matrices from numerical simulations based on the ideas of sparse sampling and compressed sensing. The method is useful for problems where the determination of the entries of a matrix constitutes the computational bottleneck. We apply this new method to an important problem in computational chemistry: the determination of molecular vibrations from elect… ▽ More This article presents a new method to compute matrices from numerical simulations based on the ideas of sparse sampling and compressed sensing. The method is useful for problems where the determination of the entries of a matrix constitutes the computational bottleneck. We apply this new method to an important problem in computational chemistry: the determination of molecular vibrations from electronic structure calculations, where our results show that the overall scaling of the procedure can be improved in some cases. Moreover, our method provides a general framework for bootstrapping cheap low-accuracy calculations in order to reduce the required number of expensive high-accuracy calculations, resulting in a significant 3x speed-up in actual calculations. △ Less

Submitted 17 October, 2014; originally announced October 2014.

Comments: 9 pages, 7 figures

More accurate and efficient bath spectral densities from super-resolution

Authors: Thomas Markovich, Samuel M. Blau, John Parkhill, Christoph Kreisbeck, Jacob N. Sanders, Xavier Andrade, Alán Aspuru-Guzik

Abstract: Quantum transport and other phenomena are typically modeled by coupling the system of interest to an environment, or bath, held at thermal equilibrium. Realistic bath models are at least as challenging to construct as models for the quantum systems themselves, since they must incorporate many degrees of freedom that interact with the system on a wide range of timescales. Owing to computational lim… ▽ More Quantum transport and other phenomena are typically modeled by coupling the system of interest to an environment, or bath, held at thermal equilibrium. Realistic bath models are at least as challenging to construct as models for the quantum systems themselves, since they must incorporate many degrees of freedom that interact with the system on a wide range of timescales. Owing to computational limitations, the environment is often modeled with simple functional forms, with a few parameters fit to experiment to yield semi-quantitative results. Growing computational resources have enabled the construction of more realistic bath models from molecular dynamics (MD) simulations. In this paper, we develop a numerical technique to construct these atomistic bath models with better accuracy and decreased cost. We apply a novel signal processing technique, known as super-resolution, combined with a dictionary of physically-motivated bath modes to derive spectral densities from MD simulations. Our approach reduces the required simulation time and provides a more accurate spectral density than can be obtained via standard Fourier transform methods. Moreover, the spectral density is provided as a convenient closed-form expression which yields an analytic time-dependent bath kernel. Exciton dynamics of the Fenna-Matthews-Olsen light-harvesting complex are simulated with a second order time-convolutionless master equation, and spectral densities constructed via super-resolution are shown to reproduce the dynamics using only a quarter of the amount of MD data. △ Less

Submitted 16 July, 2013; originally announced July 2013.

Comments: 8 pages, 3 figures

Real-space density functional theory on graphical processing units: computational approach and comparison to Gaussian basis set methods

Authors: Xavier Andrade, Alán Aspuru-Guzik

Abstract: We discuss the application of graphical processing units (GPUs) to accelerate real-space density functional theory (DFT) calculations. To make our implementation efficient, we have developed a scheme to expose the data parallelism available in the DFT approach; this is applied to the different procedures required for a real-space DFT calculation. We present results for current-generation GPUs from… ▽ More We discuss the application of graphical processing units (GPUs) to accelerate real-space density functional theory (DFT) calculations. To make our implementation efficient, we have developed a scheme to expose the data parallelism available in the DFT approach; this is applied to the different procedures required for a real-space DFT calculation. We present results for current-generation GPUs from AMD and Nvidia, which show that our scheme, implemented in the free code Octopus, can reach a sustained performance of up to 90 GFlops for a single GPU, representing a significant speed-up when compared to the CPU version of the code. Moreover, for some systems our implementation can outperform a GPU Gaussian basis set code, showing that the real-space approach is a competitive alternative for DFT simulations on GPUs. △ Less

Submitted 30 August, 2013; v1 submitted 12 June, 2013; originally announced June 2013.

Comments: 16 pages, 19 figures

Journal ref: J. Chem. Theory Comput. 2013

A survey of the parallel performance and the accuracy of Poisson solvers for electronic structure calculations

Authors: Pablo García-Risueño, Joseba Alberdi-Rodriguez, Micael J. T. Oliveira, Xavier Andrade, Michael Pippig, Javier Muguerza, Agustin Arruabarrena, Angel Rubio

Abstract: We present an analysis of different methods to calculate the classical electrostatic Hartree potential created by charge distributions. Our goal is to provide the reader with an estimation on the performance ---in terms of both numerical complexity and accuracy--- of popular Poisson solvers, and to give an intuitive idea on the way these solvers operate. Highly parallelisable routines have been im… ▽ More We present an analysis of different methods to calculate the classical electrostatic Hartree potential created by charge distributions. Our goal is to provide the reader with an estimation on the performance ---in terms of both numerical complexity and accuracy--- of popular Poisson solvers, and to give an intuitive idea on the way these solvers operate. Highly parallelisable routines have been implemented in the first-principle simulation code Octopus to be used in our tests, so that reliable conclusions about the capability of methods to tackle large systems in cluster computing can be obtained from our work. △ Less

Submitted 14 November, 2012; v1 submitted 9 November, 2012; originally announced November 2012.

Comments: 48 pages, 8 figures. A supporting info of 14 pages and 4 images

Compressed sensing for multidimensional electronic spectroscopy experiments

Authors: J. N. Sanders, S. Mostame, S. K. Saikin, X. Andrade, J. R. Widom, A. H. Marcus, A. Aspuru-Guzik

Abstract: Compressed sensing is a processing method that significantly reduces the number of measurements needed to accurately resolve signals in many fields of science and engineering. We develop a two-dimensional (2D) variant of compressed sensing for multidimensional electronic spectroscopy and apply it to experimental data. For the model system of atomic rubidium vapor, we find that compressed sensing p… ▽ More Compressed sensing is a processing method that significantly reduces the number of measurements needed to accurately resolve signals in many fields of science and engineering. We develop a two-dimensional (2D) variant of compressed sensing for multidimensional electronic spectroscopy and apply it to experimental data. For the model system of atomic rubidium vapor, we find that compressed sensing provides significantly better resolution of 2D spectra than a conventional discrete Fourier transform from the same experimental data. We believe that by combining powerful resolution with ease of use, compressed sensing can be a powerful tool for the analysis and interpretation of ultrafast spectroscopy data. △ Less

Submitted 16 July, 2012; originally announced July 2012.

Comments: 5 pages, 4 figures

Application of compressed sensing to the simulation of atomic systems

Authors: X. Andrade, J. N. Sanders, A. Aspuru-Guzik

Abstract: Compressed sensing is a method that allows a significant reduction in the number of samples required for accurate measurements in many applications in experimental sciences and engineering. In this work, we show that compressed sensing can also be used to speed up numerical simulations. We apply compressed sensing to extract information from the real-time simulation of atomic and molecular systems… ▽ More Compressed sensing is a method that allows a significant reduction in the number of samples required for accurate measurements in many applications in experimental sciences and engineering. In this work, we show that compressed sensing can also be used to speed up numerical simulations. We apply compressed sensing to extract information from the real-time simulation of atomic and molecular systems, including electronic and nuclear dynamics. We find that for the calculation of vibrational and optical spectra the total propagation time, and hence the computational cost, can be reduced by approximately a factor of five. △ Less

Submitted 29 May, 2012; originally announced May 2012.

Comments: 7 pages, 5 figures

Prediction of the derivative discontinuity in density functional theory from an electrostatic description of the exchange and correlation potential

Authors: Xavier Andrade, Alán Aspuru-Guzik

Abstract: We propose a new approach to approximate the exchange and correlation (XC) functional in density functional theory. The XC potential is considered as an electrostatic potential, generated by a fictitious XC density, which is in turn a functional of the electronic density. We apply the approach to develop a correction scheme that fixes the asymptotic behavior of any approximated XC potential for fi… ▽ More We propose a new approach to approximate the exchange and correlation (XC) functional in density functional theory. The XC potential is considered as an electrostatic potential, generated by a fictitious XC density, which is in turn a functional of the electronic density. We apply the approach to develop a correction scheme that fixes the asymptotic behavior of any approximated XC potential for finite systems. Additionally, the correction procedure gives the value of the derivative discontinuity; therefore it can directly predict the fundamental gap as a ground-state property. △ Less

Submitted 25 July, 2011; v1 submitted 21 July, 2011; originally announced July 2011.

Comments: 5 pages, 4 figures

Basis set effects on the hyperpolarizability of CHCl_3: Gaussian-type orbitals, numerical basis sets and real-space grids

Authors: Fernando Vila, David Strubbe, Yoshinari Takimoto, Xavier Andrade, Angel Rubio, S. G. Louie, J. J. Rehr

Abstract: Calculations of the hyperpolarizability are typically much more difficult to converge with basis set size than the linear polarizability. In order to understand these convergence issues and hence obtain accurate ab initio values, we compare calculations of the static hyperpolarizability of the gas-phase chloroform molecule (CHCl_3) using three different kinds of basis sets: Gaussian-type orbitals,… ▽ More Calculations of the hyperpolarizability are typically much more difficult to converge with basis set size than the linear polarizability. In order to understand these convergence issues and hence obtain accurate ab initio values, we compare calculations of the static hyperpolarizability of the gas-phase chloroform molecule (CHCl_3) using three different kinds of basis sets: Gaussian-type orbitals, numerical basis sets, and real-space grids. Although all of these methods can yield similar results, surprisingly large, diffuse basis sets are needed to achieve convergence to comparable values. These results are interpreted in terms of local polarizability and hyperpolarizability densities. We find that the hyperpolarizability is very sensitive to the molecular structure, and we also assess the significance of vibrational contributions and frequency dispersion. △ Less

Submitted 30 March, 2010; originally announced March 2010.

Journal ref: J. Chem. Phys. 133, 034111 (2010)

A modified Ehrenfest formalism for efficient large-scale ab initio molecular dynamics

Authors: Xavier Andrade, Alberto Castro, David Zueco, J. L. Alonso, Pablo Echenique, Fernando Falceto, Angel Rubio

Abstract: We present in detail the recently derived ab-initio molecular dynamics (AIMD) formalism [Phys. Rev. Lett. 101 096403 (2008)], which due to its numerical properties, is ideal for simulating the dynamics of systems containing thousands of atoms. A major drawback of traditional AIMD methods is the necessity to enforce the orthogonalization of the wave-functions, which can become the bottleneck for… ▽ More We present in detail the recently derived ab-initio molecular dynamics (AIMD) formalism [Phys. Rev. Lett. 101 096403 (2008)], which due to its numerical properties, is ideal for simulating the dynamics of systems containing thousands of atoms. A major drawback of traditional AIMD methods is the necessity to enforce the orthogonalization of the wave-functions, which can become the bottleneck for very large systems. Alternatively, one can handle the electron-ion dynamics within the Ehrenfest scheme where no explicit orthogonalization is necessary, however the time step is too small for practical applications. Here we preserve the desirable properties of Ehrenfest in a new scheme that allows for a considerable increase of the time step while keeping the system close to the Born-Oppenheimer surface. We show that the automatically enforced orthogonalization is of fundamental importance for large systems because not only it improves the scaling of the approach with the system size but it also allows for an additional very efficient parallelization level. In this work we provide the formal details of the new method, describe its implementation and present some applications to some test systems. Comparisons with the widely used Car-Parrinello molecular dynamics method are made, showing that the new approach is advantageous above a certain number of atoms in the system. The method is not tied to a particular wave-function representation, making it suitable for inclusion in any AIMD software package. △ Less

Submitted 15 April, 2009; v1 submitted 15 December, 2008; originally announced December 2008.

Comments: 28 pages, 5 figures, published in a special issue of J. Chem. Theory Comp. in honour of John Perdew

Journal ref: Journal of Chemical Theory and Computation 5 (2009) 728-742

Cluster-surface and cluster-cluster interactions: Ab initio calculations and modeling of van der Waals forces

Authors: Silvana Botti, Alberto Castro, Xavier Andrade, Angel Rubio, Miguel A. L. Marques

Abstract: We present fully ab-initio calculations of van der Waals coefficients for two different situations: i) the interaction between hydrogenated silicon clusters; and ii) the interactions between these nanostructures and a non metallic surface (a silicon or a silicon carbide surface). The methods used are very efficient, and allow the calculation of systems containing hundreds of atoms. The results o… ▽ More We present fully ab-initio calculations of van der Waals coefficients for two different situations: i) the interaction between hydrogenated silicon clusters; and ii) the interactions between these nanostructures and a non metallic surface (a silicon or a silicon carbide surface). The methods used are very efficient, and allow the calculation of systems containing hundreds of atoms. The results obtained are further analyzed and understood with the help of simple models. These models can be of interest for molecular dynamics simulations of silicon nanostructures on surfaces, where they can give a very fast yet sufficiently accurate determination of the van der Waals interaction at large separations. △ Less

Submitted 18 June, 2008; originally announced June 2008.

Comments: Phys. Rev. B

Efficient formalism for large scale ab initio molecular dynamics based on time-dependent density functional theory

Authors: J. L. Alonso, Xavier Andrade, Pablo Echenique, Fernando Falceto, Diego Prada-Gracia, Angel Rubio

Abstract: A new "on the fly" method to perform Born-Oppenheimer ab initio molecular dynamics (AIMD) is presented. Inspired by Ehrenfest dynamics in time-dependent density functional theory, the electronic orbitals are evolved by a Schroedinger-like equation, where the orbital time derivative is multiplied by a parameter. This parameter controls the time scale of the fictitious electronic motion and speeds… ▽ More A new "on the fly" method to perform Born-Oppenheimer ab initio molecular dynamics (AIMD) is presented. Inspired by Ehrenfest dynamics in time-dependent density functional theory, the electronic orbitals are evolved by a Schroedinger-like equation, where the orbital time derivative is multiplied by a parameter. This parameter controls the time scale of the fictitious electronic motion and speeds up the calculations with respect to standard Ehrenfest dynamics. In contrast to other methods, wave function orthogonality needs not be imposed as it is automatically preserved, which is of paramount relevance for large scale AIMD simulations. △ Less

Submitted 29 August, 2008; v1 submitted 17 October, 2007; originally announced October 2007.

Comments: 5 pages, 3 color figures, revtex4 package

Journal ref: Physical Review Letters 101 (2008) 096403

A time-dependent density functional theory scheme for efficient calculations of dynamic (hyper)polarizabilities

Authors: Xavier Andrade, Silvana Botti, Miguel A. L. Marques, Angel Rubio

Abstract: We present an efficient perturbative method to obtain both static and dynamic polarizabilities and hyperpolarizabilities of complex electronic systems. This approach is based on the solution of a frequency dependent Sternheimer equation, within the formalism of time-dependent density functional theory, and allows the calculation of the response both in resonance and out of resonance. Furthermore… ▽ More We present an efficient perturbative method to obtain both static and dynamic polarizabilities and hyperpolarizabilities of complex electronic systems. This approach is based on the solution of a frequency dependent Sternheimer equation, within the formalism of time-dependent density functional theory, and allows the calculation of the response both in resonance and out of resonance. Furthermore, the excellent scaling with the number of atoms opens the way to the investigation of response properties of very large molecular systems. To demonstrate the capabilities of this method, we implemented it in a real-space (basis-set free) code, and applied it to benchmark molecules, namely CO, H2O, and paranitroaniline (PNA). Our results are in agreement with experimental and previous theoretical studies, and fully validate our approach. △ Less

Submitted 4 April, 2007; v1 submitted 25 January, 2007; originally announced January 2007.

Comments: 9 pages, 4 figures