PubMed Health⌕ Search

Biomedical subjects

Jürg Hutter

Publications and source records attributed to Jürg Hutter.

13 recordsLinked to original sources
PubMed2006

Ionic liquids from Car-Parrinello simulations, part I: liquid AlCl3.

The properties of isolated AlCl3 clusters and the bulk system are investigated by means of static and dynamic electronic structure methods. We find important structural motifs with the edge connectivity dominant in a dimer and the corner connectivity dominant in a trimer. Furthermore, the trimer cluster exhibits an interesting ring structure with large cooperative effects relative to the dimer. Comparing the found structural motifs in isolated molecule calculations with the structure of the liquid allows us to determine the dominance of edge connectivity in the liquid. The size of the clusters present in the liquid indicates indeed that the dimer is the most abundant species, but there are also trimers, tetramers, and pentamers present. From the local dipole analysis both for the isolated clusters as well as for the liquid, further proof for the edge connectivity is given. However, all results point to the fact that there is also some small percentage of corner connectivity present that might be attributed to the most stable corner-connected cluster, namely the trimer. Importantly, we find that energetic considerations of isolated (static) clusters only do not represent the findings in liquid phase. Instead, a quantum cluster equilibrium approach or simulations are needed.

Journal Article↗
PubMed2006

Molecular dynamics simulation of liquid water: hybrid density functionals.

The structure, dynamical, and electronic properties of liquid water utilizing different hybrid density functionals were tested within the plane wave framework of first-principles molecular dynamics simulations. The computational approach, which employs modified functionals with short-ranged Hartree-Fock exchange, was first tested in calculations of the structural and bonding properties of the water dimer and cyclic water trimer. Liquid water simulations were performed at the state point of 350 K at the experimental density. Simulations included three different hybrid functionals, a meta-functional, four gradient-corrected functionals, and the local density and Hartree-Fock approximations. It is found that hybrid functionals are superior in reproducing the experimental structure and dynamical properties as measured by the radial distribution function and self-diffusion constant when compared to the pure density functionals. The local density and Hartree-Fock approximations show strongly over- and understructured liquids, respectively. Hydrogen bond analysis shows that the hybrid functionals give slightly smaller average numbers of hydrogen bonds than pure density functionals but similar hydrogen bond populations. The average molecular dipole moments in the liquid from the three hybrid functionals are lower than those of the corresponding pure density functionals.

Algorithms↗
PubMed2006

Simulating fluid-phase equilibria of water from first principles.

Efficient Monte Carlo algorithms and a mixed-basis set electronic structure program were used to compute from first principles the vapor-liquid coexistence curve of water. A water representation based on the Becke-Lee-Yang-Parr exchange and correlation functionals yields a saturated liquid density of 900 kg/m3 at 323 K and normal boiling and critical temperatures of 350 and 550 K, respectively. An analysis of the structural and electronic properties of the saturated liquid phase shows an increase of the asymmetry of the local hydrogen-bonded structure despite the persistence of a 4-fold coordination and decreases of the molecular dipole moment and of the spread of the lowest unoccupied molecular orbital with increasing temperature.

Journal Article↗
PubMed2005

Isobaric-isothermal monte carlo simulations from first principles: application to liquid water at ambient conditions.

A series of first-principles Monte Carlo simulations in the isobaric-isothermal ensemble were carried out for liquid water at ambient conditions (T=298 K and p=1 atm). The Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and norm-conserving Goedecker-Teter-Hutter (GTH) pseudopotentials were employed with the CP2 K simulation package to examine systems consisting of 64 water molecules. The fluctuations in the system volume encountered in simulations in the isobaric-isothermal ensemble require a reconsideration of the suitability of the typical charge-density cutoff and the regular grid-generation method previously used for the computation of the electrostatic energy in first-principles simulations in the microcanonical or canonical ensembles. In particular, it is noted that a much higher cutoff is needed and that the most computationally efficient method of creating grids can result in poor simulations. Analysis of the simulation trajectories using a very large charge-density cutoff at 1200 Ry and four different grid-generation methods point to a significantly underestimated liquid density of about 0.8 g cm-3 resulting in a somewhat understructured liquid (with a value of about 2.7 for the height of the first peak in the oxygen-oxygen radial distribution function) for BLYP-GTH water at ambient conditions. In addition, a simulation using a charge-density cutoff at 280 Ry yields a higher density of 0.9 g cm-3, showing the sensitivity of the simulation outcome to this parameter.

Journal Article↗
PubMed2005

A density-functional approach to polarizable models: a Kim-Gordon response density interaction potential for molecular simulations.

A combined linear-response-frozen electron-density model has been implemented in a molecular-dynamics scheme derived from an extended Lagrangian formalism. This approach is based on a partition of the electronic charge distribution into a frozen region described by Kim-Gordon theory [J. Chem. Phys. 56, 3122 (1972); J. Chem. Phys. 60, 1842 (1974)] and a response contribution determined by the instantaneous ionic configuration of the system. The method is free from empirical pair potentials and the parametrization protocol involves only calculations on properly chosen subsystems. We apply this method to a series of alkali halides in different physical phases and are able to reproduce experimental structural and thermodynamic properties with an accuracy comparable to Kohn-Sham density-functional calculations.

Journal Article↗
PubMed2005

Car-Parrinello molecular dynamics study of the initial dinitrogen reduction step in Sellmann-type nitrogenase model complexes.

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.

Hydrogen↗
PubMed2005

The influence of temperature and density functional models in ab initio molecular dynamics simulation of liquid water.

The performance of density functional theory methods for the modeling of condensed aqueous systems is hard to predict and validation by ab initio molecular simulation of liquid water is absolutely necessary. In order to assess the reliability of these tests, the effect of temperature on the structure and dynamics of liquid water has been characterized with 16 simulations of 20 ps in the temperature range of 280-380 K. We find a pronounced influence of temperature on the pair correlation functions and on the diffusion constant including nonergodic behavior on the time scale of the simulation in the lower temperature range (which includes ambient temperature). These observations were taken into account in a consistent comparison of a series of density functionals (BLYP, PBE, TPSS, OLYP, HCTH120, HCTH407). All simulations were carried out using an ab initio molecular dynamics approach in which wave functions are represented using Gaussians and the density is expanded in an auxiliary basis of plane waves. Whereas the first three functionals show similar behavior, it is found that the latter three functionals yield more diffusive dynamics and less structure.

Journal Article↗
PubMed2004

A photochemical activation scheme of inert dinitrogen by dinuclear Ru(II) and Fe(II) complexes.

A general photochemical activation process of inert dinitrogen coordinated to two metal centers is presented on the basis of high-level DFT and ab initio calculations. The central feature of this activation process is the occupation of an antibonding pi* orbital upon electronic excitation from the singlet ground state S0 to the first excited singlet state S1. Populating the antibonding LUMO weakens the triple bond of dinitrogen. After a vertical excitation, the excited complex may structurally relax in the S1 state and approaches its minimum structure in the S1 state. This excited-state minimum structure features the dinitrogen bound in a diazenoid form, which exhibits a double bond and two lone pairs localized at the two nitrogen atoms, ready to be protonated. Reduction and de-excitation then yield the corresponding diazene complex; its generation represents the essential step in a nitrogen fixation and reduction protocol. The consecutive process of excitation, protonation, and reduction may be rearranged in any experimentally appropriate order. The protons needed for the reaction from dinitrogen to diazene can be provided by the ligand sphere of the complexes, which contains sulfur atoms acting as proton acceptors. These protonated thiolate functionalities bring protons close to the dinitrogen moiety. Because protonation does not change the pi*-antibonding character of the LUMO, the universal and well-directed character of the photochemical activation process makes it possible to protonate the dinitrogen complex before it is irradiated. The pi*-antibonding LUMO plays the central role in the activation process, since the diazenoid structure was obtained by excitation from various occupied orbitals as well as by a direct two-electron reduction (without photochemical activation) of the complex; that is, the important bending of N2 towards a diazenoid conformation can be achieved by populating the pi*-antibonding LUMO.

Iron↗
PubMed2004

Solvent effects on electronic properties from Wannier functions in a dimethyl sulfoxide/water mixture.

We present an efficient implementation for the calculation of maximally localized Wannier functions (MLWFs) during parallel Car-Parrinello molecular dynamics simulations. The implementation is based on a block Jacobi method. The calculation of MLWFs results in only a moderate (10%-20%) increase in computer time. Consequently it is possible to calculate MLWFs routinely during Car-Parrinello simulations. The Wannier functions are then applied to derive molecular dipole moments of dimethyl sulfoxide (DMSO) in gas phase and aqueous solution. We observe a large increase of the local dipole moment from 3.97 to 7.39 D. This large solvent effect is caused by strong hydrogen bonding at the DMSO oxygen atom and methyl groups. Decomposing the dipole moment into local contributions from the S-O bond and the methyl groups is used to understand the electrostatic response of DMSO in aqueous solution. A scheme is given to derive charges on individual atoms from the MLWFs using the D-RESP methodology. The charges also display large solvent effects and give insight into the transferability of recent force field models for DMSO.

Journal Article↗
PubMed2004

Computational approaches to activity in rhodium-catalysed hydroformylation.

In this theoretical study on rhodium-catalysed hydroformylation we examine an unmodified hydridorhodium(I) carbonyl system a together with three variants modified by the model phosphane ligands PF3 (system b), PH3 (system c) and PMe3 (system d), which show increasing basicity on the Tolman chi parameter scale. The olefinic substrate for all systems is ethene. Based on the dissociative hydroformylation mechanism, static and dynamic quantum-mechanical approaches are made for preequilibria and the whole catalytic cycle. Agreement with experimental results was achieved with regard to the predominance of phosphane monocoordination in systems b-d, different sensitivity of unmodified and modified systems towards hydrogen pressure and the early location of the rate-determining step. Neither the catalytic cycle as a whole nor olefin insertion as an important selectivity-determining step gives a clear picture of activity differences among a-d. However, the crucial first catalytic step, association of ethene to the active species [HRhL3] (L=CO, PR3), may play the key role in the experimentally observed higher activity of a and systems with less basic phosphane ligands modelled by b.

Journal Article↗
PubMed2003

QM/MM Car-Parrinello molecular dynamics study of the solvent effects on the ground state and on the first excited singlet state of acetone in water.

We present a hybrid Car-Parrinello quantum mechanical/molecular mechanical (QM/MM) approach that is capable of treating the dynamics of molecular systems in electronically excited states in complex environments. The potential energy surface in the excited state is described either within the restricted open-shell Kohn-Sham (ROKS) formalism or within time-dependent density functional theory (TDDFT). As a test case, we apply this technique to the study of the solvent effects on the ground state and on the first excited singlet state of acetone in water. Our results demonstrate that for this system a purely classical description of the solvent is sufficient, since inclusion of the first solvent shell of 12 water molecules into the quantum system does not show a significant effect on this transition. The excited-state energies calculated with ROKS are red shifted by a constant value compared to the TDDFT results, while the relative variations of the excitation energy for different configurations are in very good agreement. The experimentally observed blue shift of the excitation energy in going from gas phase to condensed phase is well reproduced. Excited-state dynamics carried out with ROKS yield the relaxation of the solute and the rearrangement of the solvent structure on a picosecond timescale. The calculated Stokes shift is in reasonable agreement with experimental data.

Journal Article↗