BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrea Zen

Andrea Zen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phphysics.comp-phcond-mat.str-elcond-mat.mtrl-sciquant-phcond-mat.softphysics.atom-ph

Catalog footprint

What is connected

14works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

14 published item(s)

preprint2022arXiv

General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides

The O vacancy (Ov) formation energy, $E_\textrm{Ov}$, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g. the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute $E_\textrm{Ov}$ accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods such as local coupled cluster with single, double and perturbative triple excitations [CCSD(T)]. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rocksalt MgO, producing the first accurate and well-converged determinations of $E_\textrm{Ov}$ with this method. These reference values are used to benchmark exchange-correlation functionals in DFT and we find that all studied functionals underestimate $E_\textrm{Ov}$, with the average error decreasing along the rungs of Jacob's ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Long range ionic and short range hydration effects govern strongly anisotropic clay nanoparticle interactions

The aggregation of clay particles in aqueous solution is a ubiquitous everyday process of broad environmental and technological importance. However, it is poorly understood at the all-important atomistic level since it depends on a complex and dynamic interplay of solvent-mediated electrostatic, hydrogen-bonding, and dispersion interactions. With this in mind we have performed an extensive set of classical molecular dynamics simulations (included enhanced sampling simulations) on the interactions between model kaolinite nanoparticles in pure and salty water. Our simulations reveal highly anisotropic behaviour in which the interaction between the nanoparticles varies from attractive to repulsive depending on the relative orientation of the nanoparticles. Detailed analysis reveals that at large separation (>1.5 nm) this interaction is dominated by electrostatic effects whereas at smaller separations the nature of the water hydration structure becomes critical. This study highlights an incredible richness in how clay nanoparticles interact, which should be accounted for in e.g. coarse grained models of clay nanoparticle aggregation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

TurboRVB: a many-body toolkit for {\it ab initio} electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for {\it ab initio} Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC), and Diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions. The electronic wave function (WF) is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal product with spin-singlet pairing, or as a Pfaffian, including both singlet and triplet correlations. This wave function can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) ansatz, first proposed by L. Pauling and P. W. Anderson in quantum chemistry and condensed matter physics, respectively. The RVB ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization is made possible thanks to state-of-the-art stochastic algorithms for energy minimization. The code has been efficiently parallelized by using a hybrid MPI-OpenMP protocol, that is also an ideal environment for exploiting the computational power of modern GPU accelerators.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Boosting the accuracy and speed of quantum Monte Carlo: size-consistency and time-step

Diffusion Monte Carlo (DMC) simulations for fermions are becoming the standard to provide 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 function to avoid singularities near the nodal surface of the trial wavefunction. 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 {\em et al.} and of DePasquale {\em et al.} we propose a simple Green function modification that restores size-consistency to large values of time-step; substantially reducing the time-step errors. The new 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insight into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that in the last few decades have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state of the art computational methods, by reviewing simulations of e.g. ice nucleation or crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insight into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that in doing so the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that by improving (i.) existing interatomic potentials; and (ii.) currently available enhanced sampling methods, the community can move towards accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evidence for Stable Square Ice from Quantum Monte Carlo

Recent experiments on ice formed by water under nanoconfinement provide evidence for a two-dimensional (2D) `square ice' phase. However, the interpretation of the experiments has been questioned and the stability of square ice has become a matter of debate. Partially this is because the simulation approaches employed so far (force fields and density functional theory) struggle to accurately describe the very small energy differences between the relevant phases. Here we report a study of 2D ice using an accurate wave-function based electronic structure approach, namely Diffusion Monte Carlo (DMC). We find that at relatively high pressure square ice is indeed the lowest enthalpy phase examined, supporting the initial experimental claim. Moreover, at lower pressures a `pentagonal ice' phase (not yet observed experimentally) has the lowest enthalpy, and at ambient pressure the `pentagonal ice' phase is degenerate with a `hexagonal ice' phase. Our DMC results also allow us to evaluate the accuracy of various density functional theory exchange correlation functionals and force field models, and in doing so we extend the understanding of how such methodologies perform to challenging 2D structures presenting dangling hydrogen bonds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ice Formation on Kaolinite: Insights from Molecular Dynamics Simulations

The formation of ice affects many aspects of our everyday life as well as technologies such as cryotherapy and cryopreservation. Foreign substances almost always aid water freezing through heterogeneous ice nucleation, but the molecular details of this process remain largely unknown. In fact, insight into the microscopic mechanism of ice formation on different substrates is difficult to obtain even via state-of-the-art experimental techniques. At the same time, atomistic simulations of heterogeneous ice nucleation frequently face extraordinary challenges due to the complexity of the water-substrate interaction and the long timescales that characterize nucleation events. Here, we have investigated several aspects of molecular dynamics simulations of heterogeneous ice nucleation considering as a prototypical ice nucleating material the clay mineral kaolinite, which is of relevance in atmospheric science. We show via seeded molecular dynamics simulations that ice nucleation on the hydroxylated (001) face of kaolinite proceeds exclusively via the formation of the hexagonal ice polytype. The critical nucleus size is two times smaller than that obtained for homogeneous nucleation at the same supercooling. Previous findings suggested that the flexibility of the kaolinite surface can alter the time scale for ice nucleation within molecular dynamics simulations. However, we here demonstrate that equally flexible (or non flexible) kaolinite surfaces can lead to very different outcomes in terms of ice formation, according to whether or not the surface relaxation of the clay is taken into account. We show that very small structural changes upon relaxation dramatically alter the ability of kaolinite to provide a template for the formation of a hexagonal overlayer of water molecules at the water-kaolinite interface, and that this relaxation therefore determines the nucleation ability of this mineral.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Toward Accurate Adsorption Energetics on Clay Surfaces

Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ab-initio molecular dynamics simulation of liquid water by Quantum Monte Carlo

Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous Density Functional Theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab-initio simulations of complex chemical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Properties of Reactive Oxygen Species by Quantum Monte Carlo

The electronic properties of the oxygen molecule, in its singlet and triplet states, and of many small oxygen-containing radicals and anions have important roles in different fields of Chemistry, Biology and Atmospheric Science. Nevertheless, the electronic structure of such species is a challenge for ab-initio computational approaches because of the difficulties to correctly describe the statical and dynamical correlation effects in presence of one or more unpaired electrons. Only the highest-level quantum chemical approaches can yield reliable characterizations of their molecular properties, such as binding energies, equilibrium structures, molecular vibrations, charge distribution and polarizabilities. In this work we use the variational Monte Carlo (VMC) and the lattice regularized Monte Carlo (LRDMC) methods to investigate the equilibrium geometries and molecular properties of oxygen and oxygen reactive species. Quantum Monte Carlo methods are used in combination with the Jastrow Antisymmetrized Geminal Power (JAGP) wave function ansatz, which has been recently shown to effectively describe the statical and dynamical correlation of different molecular systems. In particular we have studied the oxygen molecule, the superoxide anion, the nitric oxide radical and anion, the hydroxyl and hydroperoxyl radicals and their corresponding anions, and the hydrotrioxyl radical. Overall, the methodology was able to correctly describe the geometrical and electronic properties of these systems, through compact but fully-optimised basis sets and with a computational cost which scales as $N^3-N^4$, where $N$ is the number of electrons. This work is therefore opening the way to the accurate study of the energetics and of the reactivity of large and complex oxygen species by first principles.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Static and dynamical correlation in diradical molecules by Quantum Monte Carlo using the Jastrow Antisymmetrized Geminal Power ansatz

Diradical molecules are essential species involved in many organic and inorganic chemical reactions. The computational study of their electronic structure is often challenging, because a reliable description of the correlation, and in particular of the static one, requires multi-reference techniques. The Jastrow correlated Antisymmetrized Geminal Power (JAGP) is a compact and efficient wave function ansatz, based on the valence-bond representation, which can be used within Quantum Monte Carlo (QMC) approaches. The AGP part can be rewritten in terms of molecular orbitals, obtaining a multi-determinant expansion with zero-seniority number. In the present work we demonstrate the capability of the JAGP ansatz to correctly describe the electronic structure of two diradical prototypes: the orthogonally twisted ethylene, C2H4, and the methylene, CH2, representing respectively a homosymmetric and heterosymmetric system. On the other hand, we show that the simple ansatz of a Jastrow correlated Single Determinant (JSD) wave function is unable to provide an accurate description of the electronic structure in these diradical molecules, both at variational level and, more remarkably, in the fixed-nodes projection schemes showing that a poor description of the static correlation yields an inaccurate nodal surface. The suitability of JAGP to correctly describe diradicals with a computational cost comparable with that of a JSD calculation, in combination with a favorable scalability of QMC algorithms with the system size, opens new perspectives in the ab initio study of large diradical systems, like the transition states in cycloaddition reactions and the thermal isomerization of biological chromophores.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Molecular properties by Quantum Monte Carlo: an investigation on the role of the wave function ansatz and the basis set in the water molecule

Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely: the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudo potential, and the basis set for QMC calculations. We also introduce a new strategy for the definition of the atomic orbitals involved in the Jastrow - Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Optimized structure and vibrational properties by error affected potential energy surfaces

The precise theoretical determination of the geometrical parameters of molecules at the minima of their potential energy surface and of the corresponding vibrational properties are of fundamental importance for the interpretation of vibrational spectroscopy experiments. Quantum Monte Carlo techniques are correlated electronic structure methods promising for large molecules, which are intrinsically affected by stochastic errors on both energy and force calculations, making the mentioned calculations more challenging with respect to other more traditional quantum chemistry tools. To circumvent this drawback in the present work we formulate the general problem of evaluating the molecular equilibrium structures, the harmonic frequencies and the anharmonic coefficients of an error affected potential energy surface. The proposed approach, based on a multidimensional fitting procedure, is illustrated together with a critical evaluation of systematic and statistical errors. We observe that the use of forces instead of energies in the fitting procedure reduces the the statistical uncertainty of the vibrational parameters by one order of magnitude. Preliminary results based on Variational Monte Carlo calculations on the water molecule demonstrate the possibility to evaluate geometrical parameters, harmonic and anharmonic coefficients at this level of theory with an affordable computational cost and a small stochastic uncertainty (<0.07% for geometries and <0.7% for vibrational properties).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite temperature electronic simulations beyond the Born-Oppenheimer approximation

We introduce a general technique to compute finite temperature electronic properties by a novel covariant formulation of the electronic partition function. By using a rigorous variational upper bound to the free energy we are led to the evaluation of a partition function that can be computed stochastically by sampling electronic wave functions and atomic positions (assumed classical). In order to achieve this target we show that it is extremely important to consider the non trivial geometry of the space defined by the wave function ansatz. The method can be extended to any technique capable to provide an energy value over a given wave function ansatz depending on several variational parameters and atomic positions. In particular we can take into account electronic correlation, by using the standard variational quantum Monte Carlo method, that has been so far limited to zero temperature ground state properties. We show that our approximation reduces correctly to the standard Born-Oppenheimer (BO) one at zero temperature and to the correct high temperature limit. At large enough temperatures this method allows to improve the BO, providing lower values of the electronic free energy, because within this method it is possible to take into account the electron entropy. We test this new method on the simple hydrogen molecule, where at low temperature we recover the correct BO low temperature limit. Moreover, we show that the dissociation of the molecule is possible at a temperature much smaller than the BO prediction. Several extension of the proposed technique are also discussed, as for instance the calculation of critical (magnetic, superconducting) temperatures, or transition rates in chemical reactions.

0 citations0 reviewsNo code linkedOpen access