Tag: Hirshfeld-I

Predicting Partial Atomic Charges in Siliceous Zeolites

Authors: Jarod J. Wolffis, Danny E. P. Vanpoucke, Amit Sharma, Keith V. Lawler, and Paul M. Forster
Journal: Microporous Mesoporous Mater. 277, 184-196 (2019)
doi: 10.1016/j.micromeso.2018.10.028
IF(2017): 3.649
export: bibtex
pdf: <MicroporousMesoporousMater>


Partial charges in zeolites for force fields.
Graphical Abstract: Partial charges in zeolites for force fields.


Partial atomic charge, which determines the magnitude of the Coulombic non-bonding interaction, represents a critical parameter in molecular mechanics simulations. Partial charges may also be used as a measure of physical properties of the system, i.e. covalency, acidic/catalytic sites, etc. A range of methods, both empirical and ab initio, exist for calculating partial charges in a given solid, and several of them are compared here for siliceous (pure silica) zeolites. The relationships between structure and the predicted partial charge are examined. The predicted partial charges from different methods are also compared with related experimental observations, showing that a few of the methods offer some guidance towards identifying the T-sites most likely to undergo substitution or for proton localization in acidic framework forms. Finally, we show that assigning unique calculated charges to crystallographically unique framework atoms makes an appreciable difference in simulating predicting N2 and O2 adsorption with common dispersion-repulsion parameterizations.

Linker Functionalization in MIL-47(V)-R Metal–Organic Frameworks: Understanding the Electronic Structure

Authors: Danny E. P. Vanpoucke
Journal: J. Phys. Chem. C 121(14), 8014-8022 (2017)
doi: 10.1021/acs.jpcc.7b01491
IF(2017): 4.484
export: bibtex
pdf: <J.Phys.Chem.C>
Graphical Abstract: Evolution of the electronic band structure of MIL-47(V) upon OH-functionalization of the BDC linker.
Graphical Abstract: Evolution of the electronic band structure of MIL-47(V) upon OH-functionalization of the BDC linker. The π-orbital of the BDC linker splits upon functionalisation, and the split-off π-band moves up into the band gap, effectively reducing the latter.


Metal–organic frameworks (MOFs) have gained much interest due to their intrinsic tunable nature. In this work, we study how linker functionalization modifies the electronic structure of the host MOF, more specifically, the MIL-47(V)-R (R = −F, −Cl, −Br, −OH, −CH3, −CF3, and −OCH3). It is shown that the presence of a functional group leads to a splitting of the π orbital on the linker. Moreover, the upward shift of the split-off π-band correlates well with the electron-withdrawing/donating nature of the functional groups. For halide functional groups the presence of lone-pair back-donation is corroborated by calculated Hirshfeld-I charges. In the case of the ferromagnetic configuration of the host MIL-47(V+IV) material a half-metal to insulator transition is noted for the −Br, −OCH3, and −OH functional groups, while for the antiferromagnetic configuration only the hydroxy group results in an effective reduction of the band gap.

First-Principles Study of Antisite Defect Configurations in ZnGa2O4:Cr Persistent Phosphors

Authors: Arthur De Vos, Kurt Lejaeghere, Danny E. P. Vanpoucke, Jonas J. Joos, Philippe F. Smet, and Karen Hemelsoet
Journal: Inorg. Chem. 55(5), 2402-2412 (2016)
doi: 10.1021/acs.inorgchem.5b02805
IF(2016): 4.857
export: bibtex
pdf: <Inorg.Chem>
Graphical Abstract: (left) Ball-and-stick model of zinc gallate (right) density of states of Cr doped zinc gallate.
Graphical Abstract: First-principles simulations on zinc gallate solid phosphors (ZGO) containing a chromium dopant and antisite defects (left) rationalize the attractive interactions between the various elements. A large number of antisite pair configurations is investigated and compared with isolated antisite defects. Defect energies point out the stability of the antisite defects in ZGO. Local structural distortions are reported, and charge transfer mechanisms are analyzed based on theoretical density of states (right) and Hirshfeld-I charges.


Zinc gallate doped with chromium is a recently developed near-infrared emitting persistent phosphor, which is now extensively studied for in vivo bioimaging and security applications. The precise mechanism of this persistent luminescence relies on defects, in particular, on antisite defects and antisite pairs. A theoretical model combining the solid host, the dopant, and/or antisite defects is constructed to elucidate the mutual interactions in these complex materials. Energies of formation as well as dopant, and defect energies are calculated through density-functional theory simulations of large periodic supercells. The calculations support the chromium substitution on the slightly distorted octahedrally coordinated gallium site, and additional energy levels are introduced in the band gap of the host. Antisite pairs are found to be energetically favored over isolated antisites due to significant charge compensation as shown by calculated Hirshfeld-I charges. Significant structural distortions are found around all antisite defects. The local Cr surrounding is mainly distorted due to a ZnGa antisite. The stability analysis reveals that the distance between both antisites dominates the overall stability picture of the material containing the Cr dopant and an antisite pair. The findings are further rationalized using calculated densities of states and Hirshfeld-I charges.

Convergence of Atomic Charges with the Size of the Enzymatic Environment

Authors: Danny E. P. Vanpoucke, Julianna Oláh, Frank De Proft, Veronique Van Speybroeck, and Goedele Roos
Journal: J. Chem. Inf. Model. 55(3), 564-571 (2015)
doi: 10.1021/ci5006417
IF(2015): 3.657
export: bibtex
pdf: <J.Chem.Inf.Model.> 
Graphical Abstract: The influence of the cluster size and water presence on the atomic charge of active and inactive sites in Biomolecules.
Graphical Abstract: Graphical Abstract: The influence of the cluster size and water presence on the atomic charge of active and inactive sites in Bio-molecules.


Atomic charges are a key concept to give more insight into the electronic structure and chemical reactivity. The Hirshfeld-I partitioning scheme applied to the model protein human 2-cysteine peroxiredoxin thioredoxin peroxidase B is used to investigate how large a protein fragment needs to be in order to achieve convergence of the atomic charge of both, neutral and negatively charged residues. Convergence in atomic charges is rapidly reached for neutral residues, but not for negatively charged ones. This study pinpoints difficulties on the road towards accurate modeling of negatively charged residues of large bio-molecular systems in a multiscale approach.

Quasi-1D physics in metal-organic frameworks: MIL-47(V) from first principles

Authors: Danny E. P. Vanpoucke, Jan W. Jaeken, Stijn De Baerdemacker, Kurt Lejaeghere
and Veronique Van Speybroeck
Journal: Beilstein J. Nanotechnol. 5, 1738-1748 (2014)
doi: 10.3762/bjnano.5.184
IF(2014): 2.670
export: bibtex
pdf: <Beilstein> (open access)
Graphical Abstract: (left) Spin density of anti-ferromagnetic MIL-47(V) with ferromagnetic chains. (right) Electronic band structure and density of states.
Graphical Abstract: The MIL-47(V) MOF has one unpaired electron per V site. As a result, different spin configurations are possible, several of which lead to an anti-ferromagnetic state. The spin density of an antiferromagnetic state, containing only ferromagnetic chains is shown on the left. On the right, the electronic band structure of the same system is presented.


The geometric and electronic structure of the MIL-47(V) metal-organic framework (MOF) is investigated by using ab initio density functional theory (DFT) calculations. Special focus is placed on the relation between the spin configuration and the properties of the MOF. The ground state is found to be antiferromagnetic, with an equilibrium volume of 1554.70 Å3. The transition pressure of the pressure-induced large-pore-to-narrow-pore phase transition is calculated to be 82 MPa and 124 MPa for systems with ferromagnetic and antiferromagnetic chains, respectively. For a mixed system, the transition pressure is found to be a weighted average of the ferromagnetic and antiferromagnetic transition pressures. Mapping DFT energies onto a simple-spin Hamiltonian shows both the intra- and inter-chain coupling to be antiferromagnetic, with the latter coupling constant being two orders of magnitude smaller than the former, suggesting the MIL-47(V) to present quasi-1D behavior. The electronic structure of the different spin configurations is investigated and it shows that the band gap position varies strongly with the spin configuration. The valence and conduction bands show a clear V d-character. In addition, these bands are flat in directions orthogonal to VO6 chains, while showing dispersion along the the direction of the VO6 chains, similar as for other quasi-1D materials.

Tetravalent Doping of CeO2: The impact of valence electron character on group IV dopant influence

Authors: Danny E. P. Vanpoucke, Stefaan Cottenier, Veronique Van Speybroeck, Isabel Van Driessche, and Patrick Bultinck
Journal: J. Am. Ceram. Soc. 97(1), 258-266 (2014)
doi: 10.1111/jace.12650
IF(2014): 2.610
export: bibtex
pdf: <J.Am.Ceram.Soc.> <arXiv>


Fluorite CeO2 doped with group IV elements is studied within the density functional theory (DFT) and DFT + U framework. Concentration-dependent formation energies are calculated for Ce1−xZxO2 (Z = C, Si, Ge, Sn, Pb, Ti, Zr, Hf) with 0 ≤ x ≤ 0.25 and a roughly decreasing trend with ionic radius is observed. The influence of the valence and near valence electronic configuration is discussed, indicating the importance of filled d and f shells near the Fermi level for all properties investigated. A clearly different behavior of group IVa and IVb dopants is observed: the former are more suitable for surface modifications and the latter are more suitable for bulk modifications. For the entire set of group IV dopants, there exists an inverse relation between the change, due to doping, of the bulk modulus, and the thermal expansion coefficients. Hirshfeld-I atomic charges show that charge-transfer effects due to doping are limited to the nearest-neighbor oxygen atoms.

Cover Image of Journal of Computational Chemistry : Extending Hirshfeld-I

Authors: Danny E. P. Vanpoucke
Journal: J. Comput. Chem. 34(5), i-ii (2013)
doi: 10.1002/jcc.23239
IF(2013): 3.601
export: bibtex
pdf: <J.Comput.Chem.>


The image shows an isosurface of Hirshfeld-I “atoms in molecules” for Ti-doped CeO2, taken at an electron density of 0.03e/Å3, as presented by Danny E. P. Vanpoucke, Patrick Bultinck, and Isabel Van Driessche on page 405. The cubic Ce0.75Ti0.25O2 unit cell is shown along the 111 direction. The different atoms are still clearly distinguishable at this iso-surface level, and show the Ti atom in the corners to be much smaller than the Ce atoms on the sides. In this issue, this implementation of the Hirshfeld- I method for solids is published back to back with a Comment from Thomas A. Manz and the authors’ Reply.

Cover of Journal of Computational Chemistry: Volume 34, Issue 5, feb. 15, 2013

Reply to ‘Comment on “Extending Hirshfeld-I to bulk and periodic materials” ‘

Authors: Danny E. P. Vanpoucke, Isabel Van Driessche, and Patrick Bultinck
Journal: J. Comput. Chem. 34(5), 422-427 (2013)
doi: 10.1002/jcc.23193
IF(2013): 3.601
export: bibtex
pdf: <J.Comput.Chem.> <arXiv>
Graphical Abstract: Hirshfeld-I atoms-in-molecules atoms in Ti doped CeO2. Graphical Abstract:The issues raised in the preceding comment are addressed. It is shown why Hirshfeld-I is, from a theoretical point of view, a good method for defining AIM and obtaining charges. Charges for a set of ionic systems are calculated using our presented method and shown to be chemically feasable. Comparison of pseudo-density to all-electron based results shows the pseudo-densities to be sufficient to obtain all-electron quality results. Timing results for systems containing hundreds of atoms.


The issues raised in the comment by Manz are addressed through the presentation of calculated atomic charges for NaF, NaCl, MgO, SrTiO3, and La2Ce2O7, using our previously presented method for calculating Hirshfeld-I charges in solids (Vanpoucke et al., J. Comput. Chem. doi: 10.1002/jcc.23088). It is shown that the use of pseudovalence charges is sufficient to retrieve the full all-electron Hirshfeld-I charges to good accuracy. Furthermore, we present timing results of different systems, containing up to over 200 atoms, underlining the relatively low cost for large systems. A number of theoretical issues are formulated, pointing out mainly that care must be taken when deriving new atoms in molecules methods based on “expectations” for atomic charges.

Extending Hirshfeld-I to bulk and periodic materials

Authors: Danny E. P. Vanpoucke, Patrick Bultinck, and Isabel Van Driessche,
Journal: J. Comput. Chem. 34(5), 405-417 (2013)
doi: 10.1002/jcc.23088
IF(2013): 3.601
export: bibtex
pdf: <J.Comput.Chem.> <arXiv>
Graphical Abstract: Hirshfeld-I atoms-in-molecules carbon atoms in a graphene sheet. Graphical Abstract: The Hirshfeld-I method is extended to solids, allowing for the partitioning of a solid density into constituent atoms. The use of precalculated density grids makes the implementation code independent, and the use of pseudo-potential based electron density distributions is shown to give qualitatively the same results as all electron densities. Results for some simple solids/periodic systems like cerium oxide and graphene are presented.


In this work, a method is described to extend the iterative Hirshfeld-I method, generally used for molecules, to periodic systems. The implementation makes use of precalculated pseudopotential-based electron density distributions, and it is shown that high-quality results are obtained for both molecules and solids, such as ceria, diamond, and graphite. The use of grids containing (precalculated) electron densities makes the implementation independent of the solid state or quantum chemical code used for studying the system. The extension described here allows for easy calculation of atomic charges and charge transfer in periodic and bulk systems. The conceptual issue of obtaining reference densities for anions is discussed, and the delocalization problem for anionic reference densities originating from the use of a plane wave basis set is identified and handled.

New Functionalized Metal-Organic Frameworks MIL-47-X (X = -Cl, -Br, -CH3, -CF3, -OH, -OCH3): Synthesis, Characterization and CO2 Adsorption Properties

Authors: Shyam Biswas, Danny E. P. Vanpoucke, Toon Verstraelen, Matthias Vandichel, Sarah Couck, Karen Leus, Ying-Ya Liu, Michel Waroquier, Veronique Van Speybroeck, Joeri F. M. Denayer, and Pascal Van Der Voort
Journal: J. Phys. Chem. C 117(44), 22784-22796 (2013)
doi: 10.1021/jp406835n
IF(2013): 4.835
export: bibtex
pdf: <J.Phys.Chem.C>


Six new functionalized vanadium hydroxo terephthalates [VIII(OH)(BDC-X)]·n(guests) (MIL-47(VIII)-X-AS) (BDC = 1,4-benzenedicarboxylate; X = −Cl, −Br, −CH3, −CF3, −OH, −OCH3; AS = as-synthesized) along with the parent MIL-47 were synthesized under rapid microwave-assisted hydrothermal conditions (170 °C, 30 min, 150 W). The unreacted H2BDC-X and/or occluded solvent molecules can be removed by thermal activation under vacuum, leading to the empty-pore forms of the title compounds (MIL-47(VIV)-X). Except pristine MIL-47 (+III oxidation state), the vanadium atoms in all the evacuated functionalized solids stayed in the +IV oxidation state. The phase purity of the compounds was ascertained by X-ray powder diffraction (XRPD), diffuse reflectance infrared Fourier transform (DRIFT) spectroscopy, Raman, thermogravimetric (TG), and elemental analysis. The structural similarity of the filled and empty-pore forms of the functionalized compounds with the respective forms of parent MIL-47 was verified by cell parameter determination from XRPD data. TGA and temperature-dependent XRPD (TDXRPD) experiments in an air atmosphere indicate high thermal stability in the 330–385 °C range. All the thermally activated compounds exhibit significant microporosity (SBET in the 305–897 m2 g–1 range), as verified by the N2 and CO2 sorption analysis. Among the six functionalized compounds, MIL-47(VIV)-OCH3 shows the highest CO2 uptake, demonstrating the determining role of functional groups on the CO2 sorption behavior. For this compound and pristine MIL-47(VIV), Widom particle insertion simulations were performed based on ab initio calculated crystal structures. The theoretical Henry coefficients show a good agreement with the experimental values, and calculated isosurfaces for the local excess chemical potential indicate the enhanced CO2 affinity is due to two effects: (i) the interaction between the methoxy group and CO2 and (ii) the collapse of the MIL-47(VIV)-OCH3 framework.