Democritos National Simulation Center

Trieste, Italy

Democritos National Simulation Center

Trieste, Italy
SEARCH FILTERS
Time filter
Source Type

Pinilla C.,Universidad del Norte, Colombia | Pinilla C.,University College London | Acuna-Rojas M.,Universidad del Norte, Colombia | Seriani N.,Abdus Salam International Center For Theoretical Physics | And 2 more authors.
Computational Materials Science | Year: 2017

A classical force field for MgSiO3 polymorphs is presented and tested for the perovskite, post-perovskite and enstatite phases. In this force field each ion has a fixed partial charge and a fixed polarizability and the potential energy is minimised with respect to the ions’ dipole moments at each step of molecular dynamics by iterating them to self-consistency. The potential parameters are obtained by fitting to forces, stresses and energies calculated by density functional theory. The phase transition from perovskite to post-perovskite is predicted to take place at a pressure of 60 GPa at zero temperature. The potential reproduces well structural and thermodynamic properties as well as defects formation in a wide range of pressures, and therefore will be useful for the investigation of MgSiO3 based materials in a geophysical context. © 2016 Elsevier B.V.


Raji A.T.,Abdus Salam International Center For Theoretical Physics | Raji A.T.,University of Cape Town | Mazzarello R.,RWTH Aachen | Scandolo S.,Abdus Salam International Center For Theoretical Physics | And 4 more authors.
Solid State Communications | Year: 2011

The electronic structures of hexagonal closed-packed (h.c.p) titanium containing a vacancy and krypton impurity atoms at various insertion sites are calculated by first-principles methods in the framework of the density-functional theory (DFT). The density of states (DOS) for titanium containing a vacancy defect shows resonance-like features. Also, the bulk electron density decreases from ∼0.15 3 to ∼0.05 3 at the vacancy centre. Electronic structure calculations have been performed to investigate what underlies the krypton site preference in titanium. The DOS of the nearest-neighbour (NN) titanium atoms to the octahedral krypton appears to be less distorted (relative to pure titanium) when compared to the NN titanium atoms to the tetrahedral krypton. The electronic density deformation maps show that polarization of the titanium atoms is stronger when the krypton atom is located at the tetrahedral site. Since krypton is a closed-shell atom, thus precluding any bonding with the titanium atoms, we may conclude that the polarization of the electrons in the vicinity of the inserted krypton atoms and the distortion of the DOS of the NN titanium atoms to the krypton serve to indicate which defect site is preferred when a krypton atom is inserted into titanium. Based on these considerations, we conclude that the substitutional site is the most favourable one, and the octahedral is the preferred interstitial site, in agreement with recent DFT calculations of the energetics of krypton impurity sites© 2011 Elsevier B.V. All rights reserved.


Raji A.T.,University of Cape Town | Mazzarello R.,RWTH Aachen | Scandolo S.,Abdus Salam International Center For Theoretical Physics | Scandolo S.,Democritos National Simulation Center | And 3 more authors.
Physical Review B - Condensed Matter and Materials Physics | Year: 2011

We report on the migration of monovacancies, divacancies, and substitutional krypton atoms in hcp titanium using the nudged elastic band method, in the framework of ab initio density functional theory. The divacancy migration energy barrier is found to be lower than that of a monovacancy. The migration of substitutional krypton in titanium is governed by the vacancy mechanism, if there is an excess of nonequilibrium vacancies after the implantation. However, if thermal vacancies dominate, krypton atom migration is expected to proceed via the dissociation mechanism. We also present ab initio calculations on the formation of clusters of multiple substitutional impurity krypton atoms and the interactions between krypton impurities and vacancies in titanium. We have found that, analogous to multiple interstitial krypton clusters, clusters of substitutional krypton atoms are energetically unstable, but are stabilized by vacancies. However, the efficiency of stabilization by vacancies depends strongly on the spatial distribution of the vacancies within the clusters. This study indicates the possibility of inert gases to nucleate in voids created by vacancies in ion implantation processes. © 2011 American Physical Society.


Pinilla C.,Abdus Salam International Center For Theoretical Physics | Irani A.H.,Abdus Salam International Center For Theoretical Physics | Seriani N.,Abdus Salam International Center For Theoretical Physics | Scandolo S.,Abdus Salam International Center For Theoretical Physics | Scandolo S.,Democritos National Simulation Center
Journal of Chemical Physics | Year: 2012

A novel all-atom, dissociative, and polarizable force field for water is presented. The force field is parameterized based on forces, stresses, and energies obtained form ab initio calculations of liquid water at ambient conditions. The accuracy of the force field is tested by calculating structural and dynamical properties of liquid water and the energetics of small water clusters. The transferability of the force field to dissociated states is studied by considering the solvation of a proton and the ionization of water at extreme conditions of pressure and temperature. In the case of the solvated proton, the force field properly describes the presence of both Eigen and Zundel configurations. In the case of the pressure-induced ice VIII/ice X transition and the temperature-induced transition to a superionic phase, the force field is found to describe accurately the proton symmetrization and the melting of the proton sublattice, respectively. © 2012 American Institute of Physics.


Dupuy N.,CNRS Institute of Mineralogy, Materials Physics and Cosmochemistry | Bouaouli S.,University Pierre and Marie Curie | Mauri F.,CNRS Institute of Mineralogy, Materials Physics and Cosmochemistry | Sorella S.,International School for Advanced Studies | And 2 more authors.
Journal of Chemical Physics | Year: 2015

We study the ionization energy, electron affinity, and the π → π- (1La) excitation energy of the anthracene molecule, by means of variational quantum Monte Carlo (QMC) methods based on a Jastrow correlated antisymmetrized geminal power (JAGP) wave function, developed on molecular orbitals (MOs). The MO-based JAGP ansatz allows one to rigorously treat electron transitions, such as the HOMO → LUMO one, which underlies the 1La excited state. We present a QMC optimization scheme able to preserve the rank of the antisymmetrized geminal power matrix, thanks to a constrained minimization with projectors built upon symmetry selected MOs. We show that this approach leads to stable energy minimization and geometry relaxation of both ground and excited states, performed consistently within the correlated QMC framework. Geometry optimization of excited states is needed to make a reliable and direct comparison with experimental adiabatic excitation energies. This is particularly important in π-conjugated and polycyclic aromatic hydrocarbons, where there is a strong interplay between low-lying energy excitations and structural modifications, playing a functional role in many photochemical processes. Anthracene is an ideal benchmark to test these effects. Its geometry relaxation energies upon electron excitation are of up to 0.3 eV in the neutral 1La excited state, while they are of the order of 0.1 eV in electron addition and removal processes. Significant modifications of the ground state bond length alternation are revealed in the QMC excited state geometry optimizations. Our QMC study yields benchmark results for both geometries and energies, with values below chemical accuracy if compared to experiments, once zero point energy effects are taken into account. © 2015 AIP Publishing LLC.


Costanzo F.,University of Padua | Costanzo F.,Democritos National Simulation Center | Silvestrelli P.L.,University of Padua | Silvestrelli P.L.,Democritos National Simulation Center | And 2 more authors.
Archives of Metallurgy and Materials | Year: 2012

Hydrogen is frequently touted as the "fuel of the future" because of its huge potential as clean energy source, although the large-scale adoption of this technology has yet to be realized. One of the remaining barriers to the utilization of hydrogen energy is an efficient and inexpensive means of hydrogen storage. In this work we investigate the nature of this process by first principle calculation. In particular, we study the way in which the H2 molecule can interact with graphene sheet through physisorption and chemisorption mechanism. The first mechanism involves the condensation of the hydrogen molecule on the graphene as a result of weak van der Waals forces, while the chemisorption mechanism involves the preliminary dissociation of the H2 molecule and the subsequent reaction of hydrogen atoms with the unsatured C-C bonds to form C-H bonds. To study carefully the possible physisorbed configurations on the graphene sheet, we take in to account van der Waals (vdW) interactions in DFT using the new method (DFT/vdW-WF) recently developed in our group and based on the concept of maximally localized Wannier functions. There are three possible way in which the H2 molecule can adapt to the structure of graphene: the hollow, the bridge and the top site called H, B and T configurations, respectively. We find the hollow site to be most stable physisorbed state with a binding energy of -50 meV. This value, in agreement with experimental results, is also compared with other vdW-correction methods as described in the following paper. Diffusion of the physisorbed configurations on the graphene sheet and activated reaction pathways in which the molecule starts from a physisorbed configuration to end up in a chemisorbed configurations have also been studied.


Zen A.,London Center for Nanotechnology | Zen A.,University College London | Sorella S.,International School for Advanced Studies | Sorella S.,Democritos National Simulation Center | And 6 more authors.
Physical Review B - Condensed Matter and Materials Physics | Year: 2016

Diffusion Monte Carlo (DMC) simulations for fermions are becoming the standard for providing high-quality reference data in systems that are too large to be investigated via quantum chemical approaches. DMC with the fixed-node approximation relies on modifications of the Green's function to avoid singularities near the nodal surface of the trial wave function. Here we show that these modifications affect the DMC energies in a way that is not size consistent, resulting in large time-step errors. Building on the modifications of Umrigar et al. and DePasquale et al. we propose a simple Green's function modification that restores size consistency to large values of the time step, which substantially reduces time-step errors. This algorithm also yields remarkable speedups of up to two orders of magnitude in the calculation of molecule-molecule binding energies and crystal cohesive energies, thus extending the horizons of what is possible with DMC. © 2016 American Physical Society.


Rossi M.,University of Milan | Vitali E.,University of Milan | Vitali E.,Democritos National Simulation Center | Galli D.E.,University of Milan | Reatto L.,University of Milan
Journal of Physics Condensed Matter | Year: 2010

Defects are believed to play a fundamental role in the supersolid state of 4He. We have studied solid 4He in two dimensions (2D) as a function of the number of vacancies nv, up to 30, inserted in the initial configuration at ρ = 0.0765 Å-2, close to the melting density, with the exact zero-temperature shadow path integral ground state method. The crystalline order is found to be stable also in the presence of many vacancies and we observe two completely different regimes. For small nv, up to about 6, vacancies form a bound state and cause a decrease of the crystalline order. At larger nv, the formation energy of an extra vacancy at fixed density decreases by one order of magnitude to about 0.6K. It is no longer possible to recognize vacancies in the equilibrated state because they mainly transform into quantum dislocations and crystalline order is found almost independently of how many vacancies have been inserted in the initial configuration. The one-body density matrix in this latter regime shows a non-decaying large distance tail: dislocations, that in 2D are point defects, turn out to be mobile, their number is fluctuating, and they are able to induce exchanges of particles across the system mainly triggered by the dislocation cores. These results indicate that the notion of the incommensurate versus the commensurate state loses meaning for solid 4He in 2D, because the number of lattice sites becomes ill defined when the system is not commensurate. Crystalline order is found to be stable also in 3D in the presence of up to 100 vacancies. © 2010 IOP Publishing Ltd.


Devaux N.,CNRS Institute of Mineralogy, Materials Physics and Cosmochemistry | Casula M.,CNRS Institute of Mineralogy, Materials Physics and Cosmochemistry | Decremps F.,CNRS Institute of Mineralogy, Materials Physics and Cosmochemistry | Sorella S.,International School for Advanced Studies | Sorella S.,Democritos National Simulation Center
Physical Review B - Condensed Matter and Materials Physics | Year: 2015

The cerium α-γ phase transition is characterized by means of a many-body Jastrow-correlated wave function, which minimizes the variational energy of the first-principles scalar-relativistic Hamiltonian, and includes correlation effects in a nonperturbative way. Our variational ansatz accurately reproduces the structural properties of the two phases, and proves that even at temperature T=0K the system undergoes a first-order transition, with ab initio parameters which are seamlessly connected to the ones measured by experiment at finite T. We show that the transition is related to a complex rearrangement of the electronic structure, with a key role played by the p-f hybridization. The underlying mechanism unveiled by this work can hold in many Ce-bearing compounds, and more generally in other f-electron systems. © 2015 American Physical Society.


Dagrada M.,CNRS Institute of Mineralogy, Materials Physics and Cosmochemistry | Casula M.,CNRS Institute of Mineralogy, Materials Physics and Cosmochemistry | Saitta A.M.,CNRS Institute of Mineralogy, Materials Physics and Cosmochemistry | Sorella S.,International School for Advanced Studies | And 2 more authors.
Journal of Chemical Theory and Computation | Year: 2014

We report an extensive theoretical study of the protonated water dimer H5O2 + (Zundel ion) by means of the highly correlated variational Monte Carlo and lattice regularized Monte Carlo approaches. This system represents the simplest model for proton transfer (PT), and a correct description of its properties is essential in order to understand the PT mechanism in more complex aqueous systems. Our Jastrow correlated AGP wave function ensures an accurate treatment of electron correlation. By exploiting the advantage of contracting the primitive basis set over atomic hybrid orbitals, we are able to limit dramatically the number of variational parameters with a systematic control on the numerical precision, a crucial ingredient in order to simulate larger systems. For both energetics and geometrical properties, our QMC results are found to be in excellent agreement with state-of-the-art coupled cluster CCSD(T) techniques. A comparison with density functional theory in the PBE approximation points to the crucial role of electron correlation for a correct description of the PT in the dimer. We prove that the QMC framework used in this work is able to resolve the tiny energy differences (∼0.3 kcal/mol) and structural variations involved in proton transfer reactions. Our approach combines these features and a favorable N 4 scaling with the number of particles which paves the way to the simulation of more realistic PT models. A test calculation on a larger protonated water cluster is carried out. The QMC approach used here represents a promising candidate to provide the first high-level ab initio description of PT in water. © 2014 American Chemical Society.

Loading Democritos National Simulation Center collaborators
Loading Democritos National Simulation Center collaborators