Source author record

T. Thonhauser

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

26 published item(s)

preprint2016arXiv

Chemistry in confined spaces: Reactivity of the Zn-MOF-74 channels

Using infrared spectroscopy combined with ab initio methods we study reactions of H$_2$O and CO inside the confined spaces of Zn-MOF-74 channels. Our results show that, once the water dissociation reaction H$_2$O$\;\rightarrow\;$OH+H takes place at the metal centers, the addition of 40 Torr of CO at 200 $^{\circ}$C starts the production of formic acid via OH+H+CO$\;\rightarrow\;$HCO$_2$H. Our detailed analysis shows that the overall reaction H$_2$O+CO$\;\rightarrow\;$HCO$_2$H takes place in the confinement of MOF-74 without an external catalyst, unlike the same reaction on flat surfaces. This discovery has several important consequences: It opens the door to a new set of catalytic reactions inside the channels of the MOF-74 system, it suggests that a recovery of the MOF's adsorption capacity is possible after it has been exposed to water (which in turn stabilizes its crystal structure), and it produces the important industrial feedstock formic acid.

preprint2016arXiv

Cluster Assisted Water Dissociation Mechanism in MOF-74 and Controlling it Using Helium

We show that the water dissociation reaction H$_2$O$\rightarrow$OH+H in the confined environment of MOF-74 channels can be precisely controlled by the addition of the noble gas He. Elucidating the entire reaction process with ab initio methods and infrared (IR) spectroscopy, we prove that the interaction between water molecules is critical to the formation of water clusters, which reduce the dissociation barrier by up to 37% and thus influence the reaction significantly. Our time-resolved IR measurements confirm that the formation of these clusters can be suppressed by introducing He gas, providing unprecedented control over water dissociation rates. Since the water dissociation reaction is the cause of the structural instability of MOF-74 in the presence of water, our finding of the reaction mechanism lays the groundwork for designing water stable versions of MOF-74 as well as understanding water related phenomena in MOFs in general.

preprint2016arXiv

Suppressing diborane production during the hydrogen release of metal borohydrides: The example of alloyed Al(BH$_4$)$_3$

Aluminum borohydride (Al(BH$_4$)$_3$) is an example of a promising hydrogen storage material with exceptional hydrogen densities by weight and volume and a low hydrogen desorption temperature. But, unfortunately, its production of diborane (B$_2$H$_6$) gases upon heating to release the hydrogen restricts its practical use. To elucidate this issue, we investigate the properties of a number of metal borohydrides with the same problem and find that the electronegativity of the metal cation is not the best descriptor of diborane production. We show that, instead, the closely related formation enthalpy is a better descriptor and we find that diborane production is an exponential function thereof. We conclude that diborane production is sufficiently suppressed for formation enthalpies of $-$80 kJ/mol BH$_4$ or lower, providing specific design guidelines to tune existing metal borohydrides or synthesize new ones. We then use first-principles methods to study the effects of Sc alloying in Al(BH$_4$)$_3$. Our results for the thermodynamic properties of the Al$_{1-x}$Sc$_x$(BH$_4$)$_3$ alloy clearly show the stabilizing effect of Sc alloying and thus the suppression of diborane production. We conclude that stabilizing Al(BH$_4$)$_3$ and similar borohydrides via alloying or other means is a promising route to suppress diborane production and thus develop viable hydrogen storage materials.

preprint2016arXiv

Understanding and Controlling Water Stability of MOF-74

Metal organic framework (MOF) materials in general, and MOF-74 in particular, have promising properties for many technologically important processes. However, their instability under humid conditions severely restricts practical use. We show that this instability and the accompanying reduction of the CO$_2$ uptake capacity of MOF-74 under humid conditions originate in the water dissociation reaction H$_2$O$\rightarrow$OH+H at the metal centers. After this dissociation, the OH groups coordinate to the metal centers, explaining the reduction in the MOF's CO$_2$ uptake capacity. This reduction thus strongly depends on the catalytic activity of MOF-74 towards the water dissociation reaction. We further show that-while the water molecules themselves only have a negligible effect on the crystal structure of MOF-74-the OH and H products of the dissociation reaction significantly weaken the MOF framework and lead to the observed crystal structure breakdown. With this knowledge, we propose a way to suppress this particular reaction by modifying the MOF-74 structure to increase the water dissociation energy barrier and thus control the stability of the system under humid conditions.

preprint2015arXiv

Spin signature of nonlocal-correlation binding in metal organic frameworks

We develop a proper nonempirical spin-density formalism for the van der Waals density functional (vdW-DF) method. We show that this generalization, termed svdW-DF, is firmly rooted in the single-particle nature of exchange and we test it on a range of spin systems. We investigate in detail the role of spin in the nonlocal-correlation driven adsorption of H$_2$ and CO$_2$ in the linear magnets Mn-MOF74, Fe-MOF74, Co-MOF74, and Ni-MOF74. In all cases, we find that spin plays a significant role during the adsorption process despite the general weakness of the molecular-magnetic responses. The case of CO$_2$ adsorption in Ni-MOF74 is particularly interesting, as the inclusion of spin effects results in an increased attraction, opposite to what the diamagnetic nature of CO$_2$ would suggest. We explain this counter-intuitive result, tracking the behavior to a coincidental hybridization of the O $p$ states with the Ni $d$ states in the down-spin channel. More generally, by providing insight on nonlocal correlation in concert with spin effects, our nonempirical svdW-DF method opens the door for a deeper understanding of weak nonlocal magnetic interactions.

preprint2014arXiv

Positional disorder in ammonia borane at ambient conditions

We solve a long-standing experimental discrepancy of NH$_3$BH$_3$, which---as a molecule---has a threefold rotational axis, but in its crystallized form at room temperature shows a fourfold symmetry about the same axis, creating a geometric incompatibility. To explain this peculiar experimental result, we study the dynamics of this system with ab initio Car-Parrinello molecular dynamics and nudged-elastic-band simulations. We find that rotations, rather than spatial static disorder, at angular velocities of 2 rev/ps---a time scale too small to be resolved by standard experimental techniques---are responsible for the fourfold symmetry.

preprint2014arXiv

Structural, elastic, thermal, and electronic response of small-molecule-loaded metal organic framework materials

We combine infrared spectroscopy, nano-indentation measurements, and \emph{ab initio} simulations to study the evolution of structural, elastic, thermal, and electronic responses of the metal organic framework MOF-74-Zn when loaded with H$_2$, CO$_2$, CH$_4$, and H$_2$O. We find that the molecular adsorption in this MOF triggers remarkable responses in all of these properties of the host material, with specific signatures for each of the guest molecules. With this comprehensive study we are able to clarify and correlate the underlying mechanisms regulating these responses with changes of the physical and chemical environment. Our findings suggest that metal organic framework materials in general, and MOF-74-Zn in particular, can be very promising materials for novel transducers and sensor applications, including highly selective small-molecule detection in gas mixtures.

preprint2014arXiv

Tuning the hydrogen desorption of Mg(BH$_4$)$_2$ through Zn alloying

We study the effect of Zn alloying on the hydrogen desorption properties of Mg(BH$_4$)$_2$ using $\it{ab\ initio}$ simulations. In particular, we investigate formation/reaction enthalpies/entropies for a number of compounds and reactions at a wide range of temperatures and Zn concentrations in Mg$_{1-x}$Zn$_x$(BH$_4$)$_2$. Our results show that the thermodynamic stability of the resulting material can be significantly lowered through Zn alloying. We find that e.g. the solid solution Mg$_{2/3}$Zn$_{1/3}$(BH$_4$)$_2$ has a reaction enthalpy for the complete hydrogen desorption of only 25.3 kJ/mol H$_2$$-$a lowering of 15 kJ/mol H$_2$ compared to the pure phase and a corresponding lowering in critical temperature of 123 K. In addition, we find that the enthalpy of mixing is rather small and show that the decrease in reaction enthalpy with Zn concentration is approximately linear.

preprint2014arXiv

van der Waals density functionals built upon the electron-gas tradition: Facing the challenge of competing interactions

The theoretical description of sparse matter attracts much interest, in particular for those ground-state properties that can be described by density functional theory (DFT). One proposed approach, the van der Waals density functional (vdW-DF) method, rests on strong physical foundations and offers simple yet accurate and robust functionals. A very recent functional within this method called vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412] stands out in its attempt to use an exchange energy derived from the same plasmon-based theory from which the nonlocal correlation energy was derived. Encouraged by its good performance for solids, layered materials, and aromatic molecules, we apply it to several systems that are characterized by competing interactions. These include the ferroelectric response in PbTiO$_3$, the adsorption of small molecules within metal-organic frameworks (MOFs), the graphite/diamond phase transition, and the adsorption of an aromatic-molecule on the Ag(111) surface. Our results indicate that vdW-DF-cx is overall well suited to tackle these challenging systems. In addition to being a competitive density functional for sparse matter, the vdW-DF-cx construction presents a more robust general purpose functional that could be applied to a range of materials problems with a variety of competing interactions.

preprint2014arXiv

van der Waals forces in density functional theory: The vdW-DF method

A density functional theory (DFT) that accounts for van der Waals (vdW) interactions in condensed matter, materials physics, chemistry, and biology is reviewed. The insights that led to the construction of the Rutgers-Chalmers van der Waals Density Functional (vdW-DF) are presented with the aim of giving a historical perspective, while also emphasising more recent efforts which have sought to improve its accuracy. In addition to technical details, we discuss a range of recent applications that illustrate the necessity of including dispersion interactions in DFT. This review highlights the value of the vdW-DF method as a general-purpose method, not only for dispersion bound systems, but also in densely packed systems where these types of interactions are traditionally thought to be negligible.

preprint2013arXiv

High-throughput screening of small-molecule adsorption in MOF

Using high-throughput screening coupled with state-of-the-art van der Waals density functional theory, we investigate the adsorption properties of four important molecules, H_2, CO_2, CH_4, and H_2O in MOF-74-M with M = Be, Mg, Al, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Sr, Zr, Nb, Ru, Rh, Pd, La, W, Os, Ir, and Pt. We show that high-throughput techniques can aid in speeding up the development and refinement of effective materials for hydrogen storage, carbon capture, and gas separation. The exploration of the configurational adsorption space allows us to extract crucial information concerning, for example, the competition of water with CO_2 for the adsorption "pockets." We find that only a few noble metals---Rh, Pd, Os, Ir, and Pt---favor the adsorption of CO_2 and hence are potential candidates for effective carbon-capture materials. Our findings further reveal significant differences in the binding characteristics of H_2, CO_2, CH_4, and H_2O within the MOF structure, indicating that molecular blends can be successfully separated by these nano-porous materials.

preprint2013arXiv

NMR study of small molecule adsorption in MOF-74-Mg

We calculate the carbon nuclear magnetic resonance (NMR) shielding for CO_2 and the hydrogen shieldings for both H_2 and H_2O inside the metal organic framework MOF-74-Mg. Our ab initio calculations are at the density functional theory level using the van der Waals including density functional vdW-DF. The shieldings are obtained while placing the small molecules throughout the structure, including the calculated adsorption site for various loading scenarios. We then explore relationships between loading, rotational and positional characteristics, and the NMR shieldings for each adsorbate. Our NMR calculations show a change in the shielding depending on adsorbate, position, and loading in a range that is experimentally observable. We further provide a simple model for the energy and the NMR shieldings throughout the cavity of the MOF. By providing this mapping of shielding to position and loading for these adsorbates, we argue that NMR probes could be used to provide additional information about the position at which these small molecules bind within the MOF, as well as the loading of the adsorbed molecule.

preprint2013arXiv

When metal organic frameworks turn into linear magnets

We investigate the existence of linear magnetism in the metal organic framework materials MOF-74-Fe, MOF-74-Co, and MOF-74-Ni, using first-principles density functional theory. MOF-74 displays regular quasi-linear chains of open-shell transition metal atoms, which are well separated. Our results show that within these chains-for all three materials-ferromagnetic coupling of significant strength occurs. In addition, the coupling in-between chains is at least one order of magnitude smaller, making these materials almost perfect 1D magnets at low temperature. The inter-chain coupling is found to be anti-ferromagnetic, in agreement with experiments. While some quasi-1D materials exist that exhibit linear magnetism-mostly complex oxides, polymers, and a few other rare material-they are typically very difficult to synthesize. The significance of our finding is that MOF-74 is very easy to synthesize and it is likely the simplest realization of the 1D Ising model in nature. MOF-74 could thus be used in future experiments to study 1D magnetism at low temperature.

preprint2012arXiv

A theoretical study of the hydrogen-storage potential of (H_2)_4CH_4 in metal organic framework materials and carbon nanotubes

The hydrogen-methane compound (H_2)_4CH_4---or for short H4M---is one of the most promising hydrogen-storage materials. This van der Waals compound is extremely rich in molecular hydrogen: 33.3 mass%, not including the hydrogen bound in CH_4; including it, we reach even 50.2 mass%. Unfortunately, H4M is not stable under ambient pressure and temperature, requiring either low temperature or high pressure. In this paper, we investigate the properties and structure of the molecular and crystalline forms of H4M, using ab initio methods based on van der Waals DFT (vdW-DF). We further investigate the possibility of creating the pressures required to stabilize H4M through external agents such as metal organic framework (MOF) materials and carbon nanotubes, with very encouraging results. In particular, we find that certain MOFs can create considerable pressure for H4M in their cavities, but not enough to stabilize it at room temperature, and moderate cooling is still necessary. On the other hand, we find that all investigated carbon nanotubes can create the high pressures required for H4M to be stable at room temperature, with direct implications for new and exciting hydrogen-storage applications.

preprint2012arXiv

Analyzing the frequency shift of physiadsorbed CO2 in metal organic framework materials

Combining first-principles density functional theory simulations with IR and Raman experiments, we determine the frequency shift of vibrational modes of CO2 when physiadsorbed in the iso-structural metal organic framework materials Mg-MOF74 and Zn-MOF74. Surprisingly, we find that the resulting change in shift is rather different for these two systems and we elucidate possible reasons. We explicitly consider three factors responsible for the frequency shift through physiabsorption, namely (i) the change in the molecule length, (ii) the asymmetric distortion of the CO$_2$ molecule, and (iii) the direct influence of the metal center. The influence of each factor is evaluated separately through different geometry considerations, providing a fundamental understanding of the frequency shifts observed experimentally.

preprint2012arXiv

Diffusion of Small Molecules in Metal Organic Framework Materials

Ab initio simulations are combined with in situ infrared spectroscopy to unveil the molecular transport of H$_2$, CO$_2$, and H$_2$O in the metal organic framework MOF-74-Mg. Our study uncovers---at the atomistic level---the major factors governing the transport mechanism of these small molecules. In particular, we identify four key diffusion mechanisms and calculate the corresponding diffusion barriers, which are nicely confirmed by time-resolved infrared experiments. We also answer a long-standing question about the existence of secondary adsorption sites for the guest molecules, and we show how those sites affect the macroscopic diffusion properties. Our findings are important to gain a fundamental understanding of the diffusion processes in these nano-porous materials, with direct implications for the usability of MOFs in gas sequestration and storage applications.

preprint2012arXiv

Molecular Biology at the Quantum Level: Can Modern Density Functional Theory Forge the Path?

Recent years have seen vast improvements in the ability of rigorous quantum-mechanical methods to treat systems of interest to molecular biology. In this review article, we survey common computational methods used to study such large, weakly bound systems, starting from classical simulations and reaching to quantum chemistry and density functional theory. We sketch their underlying frameworks and investigate their strengths and weaknesses when applied to potentially large biomolecules. In particular, density functional theory---a framework that can treat thousands of atoms on firm theoretical ground---can now accurately describe systems dominated by weak van der Waals interactions. This newfound ability has rekindled interest in using this tried-and-true approach to investigate biological systems of real importance. In this review, we focus on some new methods within density functional theory that allow for accurate inclusion of the weak interactions that dominate binding in biological macromolecules. Recent work utilizing these methods to study biologically-relevant systems will be highlighted, and a vision for the future of density functional theory within molecular biology will be discussed.

preprint2012arXiv

Structure and energetics of a ferroelectric organic crystal of phenazine and chloranilic acid

We report first-principles calculations for a ferroelectric organic crystal of phenazine and chloranilic acid molecules. Weak intermolecular interactions are properly treated by using a second version of van der Waals density functional known as vdW-DF2 [K. Lee et al., Phys. Rev. B 82, 081101 (2010)]. Lattice constants, total energies, spontaneous electric polarizations, phonon modes and frequencies, and the energy barrier of proton transfer are calculated and compared with PBE and experiments whenever possible. We show that the donation of one proton from a chloranilic acid molecule to a neighboring phenazine molecule is energetically favorable. This proton transfer is the key structural change that breaks the centrosymmetry and leads to the ferroelectric structure. However, there is no unstable phonon associated with the proton transfer, and an energy barrier of 8 meV is found between the paraelectric and ferroelectric states.

preprint2012arXiv

Tuning the Gate Opening Pressure of Metal Organic Frameworks (MOFs) for the Selective Separation of Hydro-carbons

Separation of hydrocarbons is one of the most energy demanding processes. The need to develop materials for the selective adsorption of hydrocarbons, under reasonable conditions, is therefore of paramount importance. This work unveils unexpected hydrocarbon selectivity in a flexible Metal Organic Framework (MOF), based on differences in their gate opening pressure. We show selectivity dependence on both chain length and specific framework-gas interaction. Combining Raman spectroscopy and theoretical van der Waals Density Functional (vdW-DF) calculations, the separation mechanisms governing this unexpected gate opening behavior are revealed.

preprint2011arXiv

Ab initio energetics and kinetics study of H_2 and CH_4 in the SI Clathrate Hydrate

We present ab initio results at the density functional theory level for the energetics and kinetics of H_2 and CH_4 in the SI clathrate hydrate. Our results complement a recent article by some of the authors [G. Roman-Perez et al., Phys. Rev. Lett. 105, 145901 (2010)] in that we show additional results of the energy landscape of H_2 and CH_4 in the various cages of the host material, as well as further results for energy barriers for all possible diffusion paths of H_2 and CH_4 through the water framework. We also report structural data of the low-pressure phase SI and the higher-pressure phases SII and SH.

preprint2011arXiv

Theory of Orbital Magnetization in Solids

In this review article, we survey the relatively new theory of orbital magnetization in solids-often referred to as the "modern theory of orbital magnetization"-and its applications. Surprisingly, while the calculation of the orbital magnetization in finite systems such as atoms and molecules is straight forward, in extended systems or solids it has long eluded calculations owing to the fact that the position operator is ill-defined in such a context. Approaches that overcome this problem were first developed in 2005 and in the first part of this review we present the main ideas reaching from a Wannier function approach to semi-classical and finite-temperature formalisms. In the second part, we describe practical aspects of calculating the orbital magnetization, such as taking k-space derivatives, a formalism for pseudopotentials, a single k-point derivation, a Wannier interpolation scheme, and DFT specific aspects. We then show results of recent calculations on Fe, Co, and Ni. In the last part of this review, we focus on direct applications of the orbital magnetization. In particular, we will review how properties such as the nuclear magnetic resonance shielding tensor and the electron paramagnetic resonance g-tensor can elegantly be calculated in terms of a derivative of the orbital magnetization.

preprint2011arXiv

Van der Waals interactions in the ground state of Mg(BH4)2 from density functional theory

In order to resolve an outstanding discrepancy between experiment and theory regarding the ground-state structure of Mg(BH4)2, we examine the importance of long-range dispersive interactions on the compound's thermodynamic stability. Careful treatment of the correlation effects within a recently developed nonlocal van der Waals density functional (vdW-DF) leads to a good agreement with experiment, favoring the α-Mg(BH4)2 phase (P6122) and a closely related Mn(BH4)2-prototype phase (P3112) over a large set of polymorphs at low temperatures. Our study demonstrates the need to go beyond (semi)local density functional approximations for a reliable description of crystalline high-valent metal borohydrides.

preprint2011arXiv

vdW-DF Study of energetic, structural, and vibrational properties of small water clusters and ice Ih

We present results for a density functional theory study of small water clusters and hexagonal ice Ih, using the van der Waals density functional (vdW-DF). In particular, we examine energetic, structural, and vibrational properties of these systems. Our results exhibit excellent agreement with both experiment and quantum-chemistry calculations and show a systematic and consistent improvement over standard exchange-correlation functionals---making vdW-DF a promising candidate for resolving longstanding difficulties with density functional theory in describing water. In addition, by comparing our vdW-DF results with quantum-chemistry calculations and other standard exchange-correlation functionals, we shed light on the question of why standard functionals often fail to describe these systems accurately.

preprint2011arXiv

Wannier-based calculation of the orbital magnetization in crystals

We present a first-principles scheme that allows the orbital magnetization of a magnetic crystal to be evaluated accurately and efficiently even in the presence of complex Fermi surfaces. Starting from an initial electronic-structure calculation with a coarse ab initio k-point mesh, maximally localized Wannier functions are constructed and used to interpolate the necessary k-space quantities on a fine mesh, in parallel to a previously-developed formalism for the anomalous Hall conductivity [X.Wang, J. Yates, I. Souza, and D. Vanderbilt, Phys. Rev. B 74, 195118 (2006)]. We formulate our new approach in a manifestly gauge-invariant manner, expressing the orbital magnetization in terms of traces over matrices in Wannier space. Since only a few (e.g., of the order of 20) Wannier functions are typically needed to describe the occupied and partially occupied bands, these Wannier matrices are small, which makes the interpolation itself very efficient. The method has been used to calculate the orbital magnetization of bcc Fe, hcp Co, and fcc Ni. Unlike an approximate calculation based on integrating orbital currents inside atomic spheres, our results nicely reproduce the experimentally measured ordering of the orbital magnetization in these three materials.

preprint2010arXiv

An ab-initio converse NMR approach for pseudopotentials

We extend the recently developed converse NMR approach [T. Thonhauser, D. Ceresoli, A. Mostofi, N. Marzari, R. Resta, and D. Vanderbilt, J. Chem. Phys. \textbf{131}, 101101 (2009)] such that it can be used in conjunction with norm-conserving, non-local pseudopotentials. This extension permits the efficient ab-initio calculation of NMR chemical shifts for elements other than hydrogen within the convenience of a plane-wave pseudopotential approach. We have tested our approach on several finite and periodic systems, finding very good agreement with established methods and experimental results.

preprint2007arXiv

Van der Waals density functional: Self-consistent potential and the nature of the van der Waals bond

We derive the exchange-correlation potential corresponding to the nonlocal van der Waals density functional [M. Dion, H. Rydberg, E. Schroder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)]. We use this potential for a self-consistent calculation of the ground state properties of a number of van der Waals complexes as well as crystalline silicon. For the latter, where little or no van der Waals interaction is expected, we find that the results are mostly determined by semilocal exchange and correlation as in standard generalized gradient approximations (GGA), with the fully nonlocal term giving little effect. On the other hand, our results for the van der Waals complexes show that the self-consistency has little effect at equilibrium separations. This finding validates previous calculations with the same functional that treated the fully nonlocal term as a post GGA perturbation. A comparison of our results with wave-function calculations demonstrates the usefulness of our approach. The exchange-correlation potential also allows us to calculate Hellmann-Feynman forces, hence providing the means for efficient geometry relaxations as well as unleashing the potential use of other standard techniques that depend on the self-consistent charge distribution. The nature of the van der Waals bond is discussed in terms of the self-consistent bonding charge.