Barbara Kirchner

Find an error

Name:
Organization: Lehrstuhl für Theoretische Chemie der Universit?t Bonn , Germany
Department: Mulliken Centre for Theoretical Chemistry
Title: (PhD)

TOPICS

Co-reporter:Oldamur Hollóczki, Albrecht Berkessel, Julian Mars, Markus Mezger, Anton Wiebe, Siegfried R. Waldvogel, and Barbara Kirchner
ACS Catalysis March 3, 2017 Volume 7(Issue 3) pp:1846-1846
Publication Date(Web):January 31, 2017
DOI:10.1021/acscatal.6b03090
In the present contribution, we investigated catalytically active mixtures of 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP) and aqueous H2O2 by molecular dynamics simulations. It is clearly observable that the HFIP molecule strongly binds to the H2O2, which is necessary for the desired catalytic reaction to occur. Upon the addition of the substrate cyclooctene to the solution, this interaction is enhanced, which suggests that the catalytic activity is increased by the presence of the hydrocarbon. We could clearly observe the microheterogeneous structure of the mixture, which is the result of the separation of the hydroxyl groups, water, and H2O2 from the fluorinated alkyl moiety in the form of large domains, which span through large areas of the system. The hydrocarbon, however, does not fit into either one of these two microphases, and it forms separate aggregates in the macroscopically homogeneous liquid, creating thereby a triphilic mixture. The latter kinds of aggregates are mostly surrounded by the fluorous moieties, and therefore, the H2O2 has to move from the polar through the fluorous domain to be able to react with the cyclooctene. Accordingly, the present reaction should be described figuratively as a phase transfer or an interfacial reaction, rather than a homogeneous liquid-phase process.Keywords: 1,1,1,3,3,3-hexafluoro-2-propanol; cyclooctene; hydrocarbon; hydrogen peroxide;
Co-reporter:Alexander Korotkevich, Dzmitry S. Firaha, Agilio A.H. Padua, Barbara Kirchner
Fluid Phase Equilibria 2017 Volume 448(Volume 448) pp:
Publication Date(Web):25 September 2017
DOI:10.1016/j.fluid.2017.03.024
Deep eutectic solvents (DESs) are mixtures of ionic compounds and molecular hydrogen bond donors. Due to the many components and their different interacting subgroups, they give rise to a plethora of many different interactions which can be studied by ab initio molecular dynamics simulations, because within this method all the forces are calculated on the fly and no parametrization prior to the calculation is necessary. Since DESs can be applied in gas capture, for example for SO2 absorption, we performed ab initio molecular dynamics studies of both the pure choline chloride/glycerol DES and the same mixed with SO2. We identified the hydrogen bonding and other specific interactions between all components. With addition of SO2, we observed a decrease in the anion-OH group interplay, because the chloride anions form complexes with the SO2 molecules. Furthermore, the SO2 molecules are incorporated into the hydrophobic network and the interaction between the hydrogen bonds of all OH groups remain constant. This decrease of anion-OH interaction might be responsible for the more fluid state of the SO2-DESs mixture than the pure DES.
Co-reporter:Uli Fechner;Chris de Graaf;Andrew E. Torda
Journal of Cheminformatics 2016 Volume 8( Issue 1 Supplement) pp:
Publication Date(Web):2016 April
DOI:10.1186/s13321-016-0119-5
Co-reporter:Henry Weber and Barbara Kirchner
The Journal of Physical Chemistry B 2016 Volume 120(Issue 9) pp:2471-2483
Publication Date(Web):February 16, 2016
DOI:10.1021/acs.jpcb.6b00098
We carried out ab initio molecular dynamics simulations for the three cyano-based ionic liquids, 1-ethyl-3-methylimidazolium tetracyanoborate ([C2C1Im][B(CN)4]), 1-ethyl-3-methyl-imidazolium dicyanamide ([C2C1Im][N(CN)2]), and 1-ethyl-3-methylimidazolium thiocyanate ([C2C1Im][SCN]). We found that the [SCN]-based ionic liquid is much more prone to π–π stacking interactions as opposed to the other two ionic liquids, contrary to the fact that all liquids bear the same cation. Hydrogen bonding is strong in the dicyanamide- and the thiocyanate-based ionic liquids and it is almost absent in the tetracyanoborate liquid. The anion prefers to stay on-top of the imidazolium ring with the highest priority for the [N(CN)2]− anion followed by the [B(CN)4]− anion. We find that experimental viscosity trends cannot be correlated to the hydrogen bond dynamics which is fastest for [B(CN)4]− followed by [SCN]− and [N(CN)2]−. For the dynamics of the cation on-top of itself, we find the order of [B(CN)4]− followed by [N(CN)2]− and finally by [SCN]−. Interestingly, this trend correlates well with the viscosity, suggesting a relation between the cation–cation dynamics and the viscosity at least for these cyano-based ionic liquids. These findings, especially the apparent correlation between cation–cation dynamics and the viscosity, might be useful for the suggestion of better ionic liquids in electrolyte applications.
Co-reporter:Dzmitry S. Firaha ;Dr. Barbara Kirchner
ChemSusChem 2016 Volume 9( Issue 13) pp:1591-1599
Publication Date(Web):
DOI:10.1002/cssc.201600126

Abstract

One of the possible solutions to prevent global climate change is the reduction of CO2 emissions, which is highly desired for the sustainable development of our society. In this work, the chemical absorption of carbon dioxide in amino acid ionic liquids was studied through first-principles methods. The use of readily accessible and biodegradable amino acids as building blocks for ionic liquids makes them highly promising replacements for the widely applied hazardous aqueous solutions of amines. A detailed insight into the reaction mechanism of the CO2 absorption was obtained through state-of-the-art theoretical methods. This allowed us to determine the reason for the specific CO2 capacities found experimentally. Moreover, we have also conducted a theoretical design of ionic liquids to provide valuable insights into the precise tuning of the energetic and kinetic parameters of the CO2 absorption.

Co-reporter:Dzmitry S. Firaha ;Dr. Barbara Kirchner
ChemSusChem 2016 Volume 9( Issue 13) pp:
Publication Date(Web):
DOI:10.1002/cssc.201600837
Co-reporter:Martin Thomas
The Journal of Physical Chemistry Letters 2016 Volume 7(Issue 3) pp:509-513
Publication Date(Web):January 15, 2016
DOI:10.1021/acs.jpclett.5b02752
We present a new approach for calculating vibrational circular dichroism spectra by ab initio molecular dynamics. In the context of molecular dynamics, these spectra are given by the Fourier transform of the cross-correlation function of magnetic dipole moment and electric dipole moment. We obtain the magnetic dipole moment from the electric current density according to the classical definition. The electric current density is computed by solving a partial differential equation derived from the continuity equation and the condition that eddy currents should be absent. In combination with a radical Voronoi tessellation, this yields an individual magnetic dipole moment for each molecule in a bulk phase simulation. Using the chiral alcohol 2-butanol as an example, we show that experimental spectra are reproduced very well. Our approach requires knowing only the electron density in each simulation step, and it is not restricted to any particular electronic structure method.
Co-reporter:Barbara Kirchner and Frank Neese  
Physical Chemistry Chemical Physics 2015 vol. 17(Issue 22) pp:14268-14269
Publication Date(Web):31 Mar 2015
DOI:10.1039/C5CP90040A
A graphical abstract is available for this content
Co-reporter:Martin Thomas, Martin Brehm and Barbara Kirchner  
Physical Chemistry Chemical Physics 2015 vol. 17(Issue 5) pp:3207-3213
Publication Date(Web):15 Dec 2014
DOI:10.1039/C4CP05272B
We present the Voronoi tessellation of electron density data to obtain molecular dipole moments in bulk phase ab initio molecular dynamics simulations for the calculation of vibrational spectra. Opposed to the established scheme of maximally localized Wannier functions, this approach does not rely on computationally demanding localization procedures. Nevertheless, we show at the examples of methanol, benzene, and phenol that it provides infrared and Raman spectra of similar quality and is even superior in specific cases like the Raman spectra of benzene and phenol. We have also applied the Voronoi method to a mixture of the ionic liquid 1-ethyl-3-methylimidazolium acetate with water, and show that it is advantageous in systems with significant charge transfer.
Co-reporter:Dzmitry S. Firaha;Dr. Oldamur Hollóczki;Dr. Barbara Kirchner
Angewandte Chemie International Edition 2015 Volume 54( Issue 27) pp:7805-7809
Publication Date(Web):
DOI:10.1002/anie.201502296

Abstract

Ionic liquids (ILs), vary strongly in their interaction with CO2. We suggest simple theoretical approach to predict the CO2 absorption behavior of ILs. Strong interaction of the CO2 with the IL anions corresponds to chemical absorption whereas weak interaction indicates physical absorption. A predictive estimate with a clear distinction between physical and chemical absorption can be simply obtained according to geometries optimized in the presence of a solvation model instead of optimizing it only in gas phase as has been done to date. The resulting Gibbs free energies compare very well with experimental values and the energies were correlated with experimental capacities. Promising anions, for ionic liquids with reversible CO2 absorption properties can be defined by a reaction Gibbs free energy of absorption in the range of −30 to 16 kJ mol−1.

Co-reporter:Dzmitry S. Firaha;Dr. Oldamur Hollóczki;Dr. Barbara Kirchner
Angewandte Chemie 2015 Volume 127( Issue 27) pp:7916-7920
Publication Date(Web):
DOI:10.1002/ange.201502296

Abstract

Ionische Flüssigkeiten (ILs) variieren stark in ihrer Wechselwirkung mit CO2. Es wird ein einfaches theoretisches Verfahren vorgeschlagen, mit dem sich das Absorptionsverhalten von ILs vorhersagen lässt. Während starke Wechselwirkungen der IL-Anionen mit CO2 chemische Absorption bedeuten, weisen schwache auf physikalische Absorption hin. Eine Abschätzung mit klarer Trennung zwischen chemischer und physikalischer Absorption kann leicht erhalten werden, indem die Geometrie in Gegenwart eines Solvatationsmodells optimiert wird statt wie bisher nur in der Gasphase. Aus diesen solvatisierten Geometrien lassen sich Gibbs-Energien erhalten, die sehr gut mit experimentellen Werten übereinstimmen. Somit können die berechneten Gibbs-Energien mit vermessenen Kapazitäten korreliert werden. Vielversprechende Anionen für den Entwurf von ILs für reversible chemische CO2-Absorption können durch berechnete Gibbs-Reaktionsenergien definiert werden. Für diese Kandidaten lassen sich wiederum die Kapazitäten aus der eben erwähnten Korrelation abschätzen.

Co-reporter:Dr. Oldamur Hollóczki;Dr. Marina Macchiagodena;Henry Weber;Martin Thomas;Dr. Martin Brehm;Dr. Annegret Stark;Dr. Olga Russina;Dr. Alessro Triolo;Dr. Barbara Kirchner
ChemPhysChem 2015 Volume 16( Issue 15) pp:3325-3333
Publication Date(Web):
DOI:10.1002/cphc.201500473

Abstract

We present here the possibility of forming triphilic mixtures from alkyl- and fluoroalkylimidazolium ionic liquids, thus, macroscopically homogeneous mixtures for which instead of the often observed two domains—polar and nonpolar—three stable microphases are present: polar, lipophilic, and fluorous ones. The fluorinated side chains of the cations indeed self-associate and form domains that are segregated from those of the polar and alkyl domains. To enable miscibility, despite the generally preferred macroscopic separation between fluorous and alkyl moieties, the importance of strong hydrogen bonding is shown. As the long-range structure in the alkyl and fluoroalkyl domains is dependent on the composition of the liquid, we propose that the heterogeneous, triphilic structure can be easily tuned by the molar ratio of the components. We believe that further development may allow the design of switchable, smart liquids that change their properties in a predictable way according to their composition or even their environment.

Co-reporter:Dr. Martin Brehm;Henry Weber;Martin Thomas;Dr. Oldamur Hollóczki;Dr. Barbara Kirchner
ChemPhysChem 2015 Volume 16( Issue 15) pp:3271-3277
Publication Date(Web):
DOI:10.1002/cphc.201500471

Abstract

In the present computational work, we develop a new tool for our trajectory analysis program TRAVIS to analyze the well-known behavior of liquids to separate into microphases. The dissection of the liquid into domains of different subsets, for example, in the case of fluorinated ionic liquids with anionic and cationic head groups (forming together the polar domain), fluorous, and alkyl subsets is followed by radical Voronoi tessellation. This leads to useful average quantities of the subset neighbor count, that is, the domain count that gives the amount of particular domains in the liquid, the domain volume and surface, as well as the isoperimetric quotient, which provides a measure of the deviation of the domains from a spherical shape. Thus, the newly implemented method allows analysis of the domains in terms of their numbers and shapes on a qualitative and also quantitative basis. It is a simple, direct, and automated analysis that does not need evaluation of the structure beforehand in terms of, for example, first solvent shell minima. In the microheterogeneous ionic liquids that we used as examples, the polar subsets always form a single domain in all investigated liquids. Although the fluorous side chains are also more or less connected in one or, maximally, two domains, the alkyl phases are dispersed.

Co-reporter:Henry Weber
The Journal of Physical Chemistry C 2015 Volume 119(Issue 27) pp:15137-15149
Publication Date(Web):June 10, 2015
DOI:10.1021/acs.jpcc.5b02347
The adsorption behavior of the ionic liquids 1,3-dimethylimidazolium thiocyanate and 1,3-dimethylimidazolium tetracyanoborate on the anatase (101) surface is studied through theoretical methods. This includes a comparison of the PBE and PBE0 functionals as well as the D3 and D3M+ dispersion correction schemes for the calculation of energetic and structural properties. The PBE functional was found to underestimate interaction energies, and the D3 dispersion correction was observed to overbind. While the surface titanium atoms predominantly form contact with anion’s nitrogen atoms, the oxygen atoms of anatase interact in hydrogen-bond-like structures with the acidic hydrogen atoms of the imidazolium cation in addition to ubiquitous electrostatic and dispersive interactions. Ionic liquid ion pairs remain stable after adsorption and closely match the geometry of the isolated ions. A band shift of the density of states of TiO2, possibly resulting in higher values of the open-circuit voltage in complete solar cells, is attributed to the adsorption, which might be interesting for application of these ionic liquids in real dye-sensitized solar cell devices.
Co-reporter:Henry Weber
The Journal of Physical Chemistry C 2015 Volume 119(Issue 45) pp:25260-25267
Publication Date(Web):October 16, 2015
DOI:10.1021/acs.jpcc.5b08538
Structural properties of the ionic liquids 1-ethyl-3-methylimidazolium thiocyanate and tetracyanoborate at the anatase (101) surface are of crucial interest for energy harvesting and storage devices, but their investigation via molecular dynamics requires large simulation cells. Thus, two classical interaction potentials with and without polarization effects were parametrized in order to accurately model the interface between the liquid and solid. The parameters were fitted to match ab initio reference interaction energies. Application of the generated force field to model the interface is demonstrated. The adsorption profiles of the ionic liquid’s components reveal a very dense surface layer with an excess of cations, which has some possible implications for the use of the ionic liquids as electrolytes, for example in solar cells.
Co-reporter:Friedrich Malberg;Oldamur Hollóczki;Martin Thomas
Structural Chemistry 2015 Volume 26( Issue 5-6) pp:1343-1349
Publication Date(Web):2015 December
DOI:10.1007/s11224-015-0662-0
An increasing lack of single ion cation–anion associations (ion pairing) in ionic liquids suggests a structural motif that stands in contradiction to the single ion pair structure of their vapor phase, which was evidenced by different experimental and theoretical studies. Therefore, a structural rearrangement has to occur en route from the liquid to the vapor. In this study, we propose a detailed four-step evaporation mechanism for ionic liquids, providing a refined perspective on the theory of this process based on the connection between ion pairing and volatility. The process involves diffusion of ions from the bulk to the surface, where they float around until a well-defined ion pair is formed with a counterion, leading to the departure from the surface into the vacuum. To assess the validity of this scheme, we performed a series of classical and ab initio molecular dynamics simulations based on the most sophisticated methods and force fields available for ionic liquids.
Co-reporter:Oldamur Hollóczki, Friedrich Malberg, Tom Welton and Barbara Kirchner  
Physical Chemistry Chemical Physics 2014 vol. 16(Issue 32) pp:16880-16890
Publication Date(Web):04 Jul 2014
DOI:10.1039/C4CP01177E
In this paper we show by using static DFT calculations and classical molecular dynamics simulations that the charge transfer between ionic liquid ions plays a major role in the observed discrepancies between the overall mobility of the ions and the observed conductivities of the corresponding ionic liquids, while it also directly suppresses the association of oppositely charged ions, thus the ion pairing. Accordingly, in electrochemical applications of these materials it is important to consider this reduction of the total charges on the ions, which can greatly affect the performance of the given process or device in which the ionic liquid is used. By slightly shifting from the salt-like to a molecular liquid-like system via the decreased charges, the charge transfer also fluidizes the ionic liquid. We believe that this vital information on the molecular level structure of ionic liquids offers a better understanding of these materials, and allows us to improve the a priori design of ionic liquids for any given purpose.
Co-reporter:Eva Perlt, Marc Brüssel and Barbara Kirchner  
Physical Chemistry Chemical Physics 2014 vol. 16(Issue 15) pp:6997-7005
Publication Date(Web):14 Feb 2014
DOI:10.1039/C3CP54797C
We introduce an alternative ab initio molecular dynamics simulation as a unification of Hartree–Fock molecular dynamics and the floating orbital approach. The general scheme of the floating orbital molecular dynamics method is presented. Moreover, a simple but sophisticated guess for the orbital centers is provided to reduce the number of electronic structure optimization steps at each molecular dynamics step. The conservation of total energy and angular momentum is investigated in order to validate the floating orbital molecular dynamics approach with and without application of the initial guess. Finally, a water monomer and a water dimer are simulated, and the influence of the orbital floating on certain properties like the dipole moment is investigated.
Co-reporter:Dzmitry S. Firaha and Barbara Kirchner
Journal of Chemical & Engineering Data 2014 Volume 59(Issue 10) pp:3098-3104
Publication Date(Web):June 6, 2014
DOI:10.1021/je500166d
We present a first principles molecular dynamics study of carbon dioxide solvation by protic ionic liquids using ethylammonium nitrate as an example solvent. Microheterogeneity of the alkyl chains and the extended hydrogen bond network could be observed. Thus, the entire structure of the investigated protic ionic liquid mixed with CO2 closely resembles the one of the pure liquid. Our data indicates that CO2 most likely creates an energy loss due to entering the liquid via the too-small voids. But this is fully compensated by specific attractive interactions of CO2 with the cation and anions of ethylammonium nitrate. This result might serve as an explanation for the question of why the volume of the ionic liquid is not increasing through CO2 uptake. The CO2 cluster formation, which shows a structure similar to supercritical CO2, is guided by the dominance of the nonpolar groups in the CO2 solvation shell.
Co-reporter:Stefan Zahn, Martin Brehm, Marc Brüssel, Oldamur Hollóczki, Miriam Kohagen, Sebastian Lehmann, Friedrich Malberg, Alfonso Sanmartin Pensado, Matthias Schöppke, Henry Weber, Barbara Kirchner
Journal of Molecular Liquids 2014 192() pp: 71-76
Publication Date(Web):
DOI:10.1016/j.molliq.2013.08.015
Co-reporter:Martin Thomas;Martin Brehm;Dr. Oldamur Hollóczki;Dr. Barbara Kirchner
Chemistry - A European Journal 2014 Volume 20( Issue 6) pp:1622-1629
Publication Date(Web):
DOI:10.1002/chem.201303329

Abstract

The solvation of the carbene 1-ethyl-3-methylimidazole-2-ylidene in the ionic liquid 1-ethyl-3-methylimidazolium acetate was investigated by ab initio molecular dynamics simulations in order to reveal the interaction between these two highly important classes of materials: N-heterocyclic carbenes with superb catalytic activity and ionic liquids with advantageous properties as solvents and reaction media. In contrast to previously published data on analogous systems, no hydrogen bond is observed between the hypovalent carbon atom and the most acidic ring hydrogen atoms, as these interaction sites of the imidazolium ring are predominantly occupied by the acetate ions. Keeping the carbene away from the ring hydrogen atoms prevents stabilization of this reactive species, and hence any retarding effect on subsequent reactions, which explains the observed high reactivity of the carbene in acetate-based ionic liquids. Instead, the carbene exhibits a weaker interaction with the methyl group of the imidazolium cation by forming a hitherto unprecedented kind of C⋅⋅⋅HC hydrogen bond. This unexpected finding not only indicates a novel kind of hydrogen bond for carbenes, but also shows that such interaction sites of the imidazolium cation are not limited to the ring hydrogen atoms. Thus, the results give the solute–solvent interactions within ionic liquids a new perspective, and provide a further, albeit weak, site of interaction to tune in order to achieve the desired environment for any dissolved active ingredient.

Co-reporter:Friedrich Malberg, Martin Brehm, Oldamur Hollóczki, Alfonso S. Pensado and Barbara Kirchner  
Physical Chemistry Chemical Physics 2013 vol. 15(Issue 42) pp:18424-18436
Publication Date(Web):02 Sep 2013
DOI:10.1039/C3CP52966E
In this work we present a comprehensive temperature-dependence analysis of both the structural and the dynamic properties of a vaporized ionic liquid (1-ethyl-3-methylimidazolium ethylsulfate). This particular ionic liquid is known to be distillable from experimental studies and thus enables us to deepen the understanding of the evaporation mechanism of ionic liquids. We have used ab initio molecular dynamics of one ion pair at three different temperatures to accurately describe the interactions present in this model ionic liquid. By means of radial and spatial distribution functions a large impact on the coordination pattern at 400 K is shown which could explain the transfer of one ion pair from the bulk to the gas phase. Comparison of the free energy surfaces at 300 K and 600 K supports the idea of bulk phase-like and gas phase-like ion pairs. The different coordination patterns caused by the temperature, describing a loosening of the anion side chains, are also well reflected in the power spectra. The lifetime analysis of typical conformations for ionic liquids shows a characteristic behavior at 400 K (temperature close to the experimental evaporation temperature), indicating that conformational changes occur when the ionic liquid is evaporated.
Co-reporter:Oldamur Hollóczki, Dzmitry S. Firaha, Joachim Friedrich, Martin Brehm, Richard Cybik, Martin Wild, Annegret Stark, and Barbara Kirchner
The Journal of Physical Chemistry B 2013 Volume 117(Issue 19) pp:5898-5907
Publication Date(Web):April 8, 2013
DOI:10.1021/jp4004399
We present a theoretical study of carbene formation from the 1-ethyl-3-methylimidazolium acetate ionic liquid in the absence and presence of CO2 in gas and liquid phase. Although CO2 physisorption constitutes a precursory step of chemisorption (the CO2’s reaction with carbenes, which forms from cations via proton abstraction by anions), it also enables a very stable CO2–anion associate. However, this counteracts the chemical absorption by reducing the basicity of the anion and the electrophilicity of the CO2, which is reflected by charge transfer. Accordingly, the observable carbene formation in the gas phase is hindered in the presence of CO2. In the neat liquid, the carbene formation is also suppressed by the charge screening compared to the case of the gas phase; nevertheless, indications for carbene incidents appear. Interestingly, in the CO2-containing liquid we detect more carbene-like incidents than in the neat one, which is caused by the way CO2 is solvated. Despite the weakness of the CO2–cation interaction, the CO2–anion associate is distorted by cations, which can be seen in longer associate distances and reduced “binding” energies. While the single solvating anion is shifted away from CO2, many more solvating cations approach it compared to the case of the gas phase. This leads to the conclusion that while the ionic liquid effect stabilizes charged species, introducing neutral species such as CO2 provides an opposite trend, leading to an inverse ionic liquid effect with the facilitation of carbene formation and thus of chemical absorption.
Co-reporter:Barbara Kirchner Priv.-Doz. Dr.
ChemPhysChem 2007 Volume 8(Issue 1) pp:41-43
Publication Date(Web):23 NOV 2006
DOI:10.1002/cphc.200600476

Proton transfer: The ability of water molecules to dissociate into hydronium and hydroxide ions is one of the most important chemical properties of liquid water. Recent work on the electron binding energies of hydrated hydroxide and hydronium ions (see picture) is highlighted.

Co-reporter:Duc-Thanh Pham, Knud Gentz, Caroline Zörlein, Nguyen T. M. Hai, Sung-Lin Tsay, Barbara Kirchner, Simone Kossmann, Klaus Wandelt and Peter Broekmann  
New Journal of Chemistry 2006 vol. 30(Issue 10) pp:1439-1451
Publication Date(Web):15 Sep 2006
DOI:10.1039/B609421J
The surface redox-chemistry of adsorbed viologens is studied by means of cyclic voltammetry (CV) in combination with in situ scanning tunneling microscopy (STM). 1,1′-Dibenzyl-4,4′-bipyridinium molecules (DBV2+) adsorb on a chloride modified Cu(100) electrode surface under formation of a laterally well ordered 2D array of supramolecular cavitand ensembles. Each cavitand consists at least of 4 individual DBV2+ sub-units which are arranged in a certain circular manner making this supramolecular cavitand chiral. Both possible enantiomeric forms are found in two mirror domains at the surface. Reducing the di-cationic DBV2+(ads) species to the corresponding radical mono-cation DBV•+(ads) causes a phase transition from the pre-existing DBV2+(ads) cavitand phase to a stripe pattern following a nucleation and growth mechanism. DBV•+(ads) species are adsorbed with their main molecular axis parallel to the surface in a side-on adsorption geometry. Enhanced intermolecular π–π-interactions are identified as the main driving force for the formation of 1D oligomer and polymer chains as the characteristic structural motif of the DBV•+(ads) phase. These structural motifs are generally independent of the electronic and structural substrate properties. Chloride desorption through the viologen film is discussed as the reason for an order–disorder transition within the viologen film at even more negative potentials.
Co-reporter:Marcella Iannuzzi, Barbara Kirchner, Jürg Hutter
Chemical Physics Letters 2006 Volume 421(1–3) pp:16-20
Publication Date(Web):3 April 2006
DOI:10.1016/j.cplett.2005.08.155

Abstract

We introduce a density functional based embedding method for the study of molecular systems in condensed phase. Molecular subunits are treated using a standard Kohn–Sham method together with an embedding potential derived from orbital-free density functional theory, by using kinetic energy functionals. The method leads to a linear scaling electronic structure approach that maps naturally onto massively parallel computers. The application of the method for a molecular dynamics simulation of water at ambient conditions results in a liquid with unstructured second solvation shell.

Co-reporter:Barbara Kirchner Dr.;Markus Reiher Priv.-Doz. Dr.;Andreas Hille Dr.;Jürg Hutter Dr.;Bernd A. Hess Dr.
Chemistry - A European Journal 2005 Volume 11(Issue 2) pp:
Publication Date(Web):18 NOV 2004
DOI:10.1002/chem.200400709

We have studied reduction reactions for nitrogen fixation at Sellmann-type model complexes with Car–Parrinello simulation techniques. These dinuclear complexes are especially designed to emulate the so-called open-side FeMoco model. The main result of this work shows that in order to obtain the reduced species several side reactions have to be suppressed. These involve partial dissociation of the chelate ligands and hydrogen atom transfer to the metal center. Working at low temperature turns out to be one necessary pre-requisite in carrying out successful events. The successful events cannot be described by simple reaction coordinates. Complicated processes are involved during the initiation of the reaction. Our theoretical study emphasizes two experimental strategies which are likely to inhibit the side reactions. Clamping of the two metal fragments by a chelating phosphane ligand should prevent dissociation of the complex. Furthermore, introduction of tert-butyl substituents could improve the solubility and should thus allow usage of a wider range of (mild) acids, reductants, and reaction conditions.

Co-reporter:Barbara Kirchner
Physics Reports (March 2007) Volume 440(Issues 1–3) pp:1-111
Publication Date(Web):1 March 2007
DOI:10.1016/j.physrep.2006.11.005
It is the aim of this work to elucidate the usefulness and feasibility of the first-principles approach and to extend it to the regime of liquid molecular substances of complex structure. Physical and thermodynamic properties of complicated liquids are investigated by means of Car–Parrinello molecular dynamics (CPMD) and also with static quantum chemical methods. The connection between the dynamic and static approach is given by the quantum cluster equilibrium (QCE) theory. Since the QCE theory is not yet well established, a new implementation in the MD post-processing program PEACEMAKER is presented. It can be shown that it is by far more important to include cooperative effects rather than to concentrate the effort on the inclusion of weak dispersion forces not present in current density functionals. Traditionally, investigations of complicated liquids were also undertaken with the tools of simple liquids, because for some problems the size of the system does not allow for a more accurate description. Although linear-scaling techniques are simplifications from the point of view of quantum chemistry, they might be severe improvements when compared to traditional molecular dynamics simulations. For the interpretation of the liquid state the introduction of local properties is inevitable. New methods are presented for the calculation of local dipole moments and for the estimation of hydrogen bond energies in quantum mechanically nondecomposable systems. The latter also allows for the detection of hydrogen bonds in simulations through a wavefunction-based criterion instead of one which is solely grounded on the geometric structure of the atomic nuclei involved. The article then discusses prominent liquids which show properties that are not yet understood. Another part of the work analyzes the effect of solvent molecules on solutes and their reactions in the solvent. Finaly, neoteric solvents, such as ionic liquids are discussed.
Co-reporter:Friedrich Malberg, Martin Brehm, Oldamur Hollóczki, Alfonso S. Pensado and Barbara Kirchner
Physical Chemistry Chemical Physics 2013 - vol. 15(Issue 42) pp:NaN18436-18436
Publication Date(Web):2013/09/02
DOI:10.1039/C3CP52966E
In this work we present a comprehensive temperature-dependence analysis of both the structural and the dynamic properties of a vaporized ionic liquid (1-ethyl-3-methylimidazolium ethylsulfate). This particular ionic liquid is known to be distillable from experimental studies and thus enables us to deepen the understanding of the evaporation mechanism of ionic liquids. We have used ab initio molecular dynamics of one ion pair at three different temperatures to accurately describe the interactions present in this model ionic liquid. By means of radial and spatial distribution functions a large impact on the coordination pattern at 400 K is shown which could explain the transfer of one ion pair from the bulk to the gas phase. Comparison of the free energy surfaces at 300 K and 600 K supports the idea of bulk phase-like and gas phase-like ion pairs. The different coordination patterns caused by the temperature, describing a loosening of the anion side chains, are also well reflected in the power spectra. The lifetime analysis of typical conformations for ionic liquids shows a characteristic behavior at 400 K (temperature close to the experimental evaporation temperature), indicating that conformational changes occur when the ionic liquid is evaporated.
Co-reporter:Oldamur Hollóczki, Friedrich Malberg, Tom Welton and Barbara Kirchner
Physical Chemistry Chemical Physics 2014 - vol. 16(Issue 32) pp:NaN16890-16890
Publication Date(Web):2014/07/04
DOI:10.1039/C4CP01177E
In this paper we show by using static DFT calculations and classical molecular dynamics simulations that the charge transfer between ionic liquid ions plays a major role in the observed discrepancies between the overall mobility of the ions and the observed conductivities of the corresponding ionic liquids, while it also directly suppresses the association of oppositely charged ions, thus the ion pairing. Accordingly, in electrochemical applications of these materials it is important to consider this reduction of the total charges on the ions, which can greatly affect the performance of the given process or device in which the ionic liquid is used. By slightly shifting from the salt-like to a molecular liquid-like system via the decreased charges, the charge transfer also fluidizes the ionic liquid. We believe that this vital information on the molecular level structure of ionic liquids offers a better understanding of these materials, and allows us to improve the a priori design of ionic liquids for any given purpose.
Co-reporter:Eva Perlt, Marc Brüssel and Barbara Kirchner
Physical Chemistry Chemical Physics 2014 - vol. 16(Issue 15) pp:
Publication Date(Web):
DOI:10.1039/C3CP54797C
Co-reporter:Barbara Kirchner and Frank Neese
Physical Chemistry Chemical Physics 2015 - vol. 17(Issue 22) pp:NaN14269-14269
Publication Date(Web):2015/03/31
DOI:10.1039/C5CP90040A
A graphical abstract is available for this content
Co-reporter:Martin Thomas, Martin Brehm and Barbara Kirchner
Physical Chemistry Chemical Physics 2015 - vol. 17(Issue 5) pp:NaN3213-3213
Publication Date(Web):2014/12/15
DOI:10.1039/C4CP05272B
We present the Voronoi tessellation of electron density data to obtain molecular dipole moments in bulk phase ab initio molecular dynamics simulations for the calculation of vibrational spectra. Opposed to the established scheme of maximally localized Wannier functions, this approach does not rely on computationally demanding localization procedures. Nevertheless, we show at the examples of methanol, benzene, and phenol that it provides infrared and Raman spectra of similar quality and is even superior in specific cases like the Raman spectra of benzene and phenol. We have also applied the Voronoi method to a mixture of the ionic liquid 1-ethyl-3-methylimidazolium acetate with water, and show that it is advantageous in systems with significant charge transfer.
1H-IMIDAZOLIUM, 1-METHYL-3-(3,3,4,4,5,5,6,6,6-NONAFLUOROHEXYL)-, IODIDE
1H-Imidazolium, 1,1'-methylenebis[3-ethyl-, dibromide
2,4,6,8-TETRAPHOSPHATRICYCLO[3.3.1.02,8]NONA-3,6-DIENE
Phosphonium, trihexyltetradecyl-
2H-Imidazol-2-ylidene, 1-ethyl-1,3-dihydro-3-methyl-
1H-Imidazolium, 1,3-dimethyl-, thiocyanate
1-METHYL-3-(3,3,4,4,5,5,6,6,7,7,8,8,8-TRIDECAFLUOROOCTYL)IMIDAZOL-1-IUM;HEXAFLUOROPHOSPHATE
1-ethyl-3-methyl-1H-imidazol-3-ium bis[(pentafluoroethyl)sulfonyl]azanide