PubMed Health⌕ Search

Biomedical subjects

Markus Reiher

Publications and source records attributed to Markus Reiher.

At least 19 recordsLinked to original sources

Trinuclear copper(II) complexes derived from Schiff-base ligands based on a 6-amino-6-deoxyglucopyranoside: structural and magnetic characterization.

The trinuclear copper(II) complexes ([CuL1)(mu-ac)Cu(mu-ac)CuL1) (1) and ([CuL2)(mu-ac)Cu(mu-ac)CuL2) (2) of the tridentate aminosaccharide-derived Schiff-base ligands H2L1 [6-N-(salicylidene)amino-6-deoxy-1,2,3-tri-O-methyl-alpha-D-glucopyranoside] and H2L2 [6-N-(3,5-di-tert-butylsalicylidene)amino-6-deoxy-1,2,3-tri-O-methyl-alpha-D-glucopyranoside] were synthesized and structurally characterized. The trinuclear complex units can be described as two terminal copper-ligand moieties bridged by a central copper acetate moiety, with the Cu centers arranged in a triangular fashion. IR and UV/vis spectroscopic studies strongly indicate that the trinuclear structure is maintained in a methanolic solution. The temperature dependence of the magnetic susceptibility of both complexes shows a moderate antiferromagnetic coupling and can be well interpreted by applying a symmetric Cua-Cub-Cua' model with linear spin topology. The fit of the magnetic data affords coupling constants J of -34 and -24 cm(-1) for 1 and 2, respectively [H = -J(SaSb + SbSa')]. For mu-alkoxo-mu-acetato-bridged copper(II) complexes with a large dihedral angle between the adjacent coordination planes, as found in 1 and 2, such an antiferromagnetic coupling is unusual. However, density functional theory calculations of 2 using BP86, B3LYP*, and B3LYP density functionals confirmed a symmetric doublet ground state.

Journal Article↗

Can Raman optical activity separate axial from local chirality? A theoretical study of helical deca-alanine.

Motivated by experimental work on the distinction of protein secondary structure motifs by Raman optical activity (ROA) spectroscopy, we demonstrate using density functional theory that axial chirality in structures with different local chirality can be filtered out by ROA spectroscopy. To this purpose, two diastereomers of right-handed helical deca-alanine, the (all-S) and the (R,S,R,S,R,S,R,S,R,S) form, are compared. Furthermore, we suggest to interpret calculated ROA spectra of large molecules in terms of vibrational bands rather than individual peaks. This is due to the non-homogeneous effect of the harmonic approximation as well as of the chosen electronic structure method onto the vibrational frequencies, which in a dense region of many vibrations will strongly determine the shape of the spectrum. In addition, the calculated ROA spectrum of (all-S)-deca-alanine is compared to the experimental spectrum of poly-(L)-alanine in solution.

Alanine↗

Intrinsic dinitrogen activation at bare metal atoms.

There is ongoing interest in metal complexes which bind dinitrogen and facilitate either its reduction or oxidation under mild conditions. In nature, the enzyme nitrogenase catalyzes this process, and dinitrogen fixation occurs under mild and ambient conditions at a metal-sulfur cluster in the center of the MoFe protein, but the mechanism of this process remains largely unknown. In the last few years, new important discoveries have been made in this field. In this review are discussed recent findings on the interaction of N(2) with metal atoms and metal-atom dimers from all groups of the periodic table as provided by gas-phase as well as matrix-isolation experiments. Intrinsic dinitrogen activation at such bare metal atoms is then related to corresponding processes at complexes, clusters, and surfaces.

Journal Article↗

Exact decoupling of the Dirac Hamiltonian. III. Molecular properties.

Recent advances in the theory of the infinite-order Douglas-Kroll-Hess (DKH) transformation of the Dirac Hamiltonian require a fresh and unified view on the calculation of atomic and molecular properties. It is carefully investigated how the four-component Dirac Hamiltonian in the presence of arbitrary electric and magnetic potentials is decoupled to two-component form. In order to cover the whole range of electromagnetic properties on the same footing, a consistent description within the DKH theory is presented. Subtle distinctions are needed between errors arising from any finite-order DKH scheme and effects due to oversimplified and thus approximate decoupling strategies for the Dirac operator, which will, though being numerically negligible in most cases, still be visible in the infinite-order limit of the two-component treatment. Special focus is given to the issue, whether the unitary DKH transformations to be applied to the Dirac Hamiltonian should depend on the property under investigation or not. It is explicitly shown that up to third order in the external potential the transformed property operator is independent of the chosen parametrization of the unitary transformations of the generalized DKH scheme. Since the standard DKH protocol covers the transformation of one-electron integrals only, the presentation is developed for one-electron properties for the sake of brevity. Nevertheless, all findings for the calculation of one-electron properties within a two-component framework presented here also hold for two-electron properties as well.

Journal Article↗

Exact decoupling of the Dirac Hamiltonian. IV. Automated evaluation of molecular properties within the Douglas-Kroll-Hess theory up to arbitrary order.

In Part III [J. Chem. Phys. 124, 064102 (2005)] of this series of papers on exact decoupling of the Dirac Hamiltonian within transformation theory, we developed the most general account on how to treat magnetic and electric properties in a unitary transformation theory on the same footing. In this paper we present an implementation of a general algorithm for the calculation of magnetic as well as electric properties within the framework of Douglas-Kroll-Hess theory. The formal and practical principles of this algorithm are described. We present the first high-order Douglas-Kroll-Hess results for property operators. As for model properties we propose to use the well-defined radial moments, i.e., expectation values of r(k), which can be understood as terms of the Taylor-series expansion of any property operator. Such moments facilitate a rigorous comparison of methods free of uncertainties which may arise in a direct comparison with experiment. This is important in view of the fact that various approaches to two-component molecular properties may yield numerically very small terms whose approximate or inaccurate treatment would not be visible in a direct comparison to experimental data or to another approximate computational reference. Results are presented for various degrees of decoupling of the model properties within the Douglas-Kroll-Hess scheme.

Journal Article↗

Construction of environment states in quantum-chemical density-matrix renormalization group calculations.

The application of the quantum-chemical density-matrix renormalization group (DMRG) algorithm is cumbersome for complex electronic structures with many active orbitals. The high computational cost is mainly due to the poor convergence of standard DMRG calculations. A factor which affects the convergence behavior of the calculations is the choice of the start-up procedure. In this start-up step matrix representations of operators have to be calculated in a guessed many-electron basis of the DMRG environment block. Different possibilities for the construction of these basis states exist, and we first compare four procedures to approximate the environment states using Slater determinants explicitly. These start-up procedures are applied to DMRG calculations on a sophisticated test system: the chromium dimer. It is found that the converged energies and the rate of convergence depend significantly on the choice of the start-up procedure. However, since already the most simple start-up procedure, which uses only the Hartree-Fock determinant, is comparatively good, Slater determinants, in general, appear not to be a good choice as approximate environment basis states for convergence acceleration. Based on extensive test calculations it is demonstrated that the computational cost can be significantly reduced if the number of total states m is successively increased. This is done in such a way that the environment states are built up stepwise from system states of previous truncated DMRG sweeps for slowly increasing m values.

Journal Article↗

Spin states in polynuclear clusters: the [Fe2O2] core of the methane monooxygenase active site.

The ability to provide a correct description of different spin states of mono- and polynuclear transition metal complexes is essential for a detailed investigation of reactions that are catalyzed by such complexes. We study the energetics of different total and local spin states of a dinuclear oxygen-bridged iron(IV) model for the intermediate Q of the hydroxylase component of methane monooxygenase by means of spin-unrestricted Kohn-Sham density functional theory. Because it is known that the spin state total energies depend systematically on the density functional, and that this dependence is intimately connected to the exact exchange admixture of present-day hybdrid functionals, we compare total energies, local and total spin values, and Heisenberg coupling constants calculated with the established functionals BP86 and B3LYP as well as with a modified B3LYP version with an exact exchange admixture ranging from 0 to 24%. It is found that exact exchange enhances local spin polarization. As the exact exchange admixture increases, the high-spin state is energetically favored, although the Broken-Symmetry state always is the ground state. Instead of the strict linear variation of the energy splittings observed for mononuclear complexes, a slightly nonlinear dependence is found. The Heisenberg coupling constants J(Fe1Fe2) --evaluated according to three different proposals from the literature -- are found to vary from -129 to -494 cm(-1) accordingly. The experimental finding that intermediate Q has an antiferromagnetic ground state is thus confirmed.

Binding Sites↗

Theoretical study of catalytic dinitrogen reduction under mild conditions.

Density functional theory results on key steps of Schrock's catalytic cycle are presented. The quantum chemical modeling of the dinitrogen-reducing reaction steps is conducted without simplifying the bulky HIPT [HIPT = 3,5-(2,4,6-iPr(3)C(6)H(2))(2)C(6)H(3)] substituents at the triamidoamine ligand.

Journal Article↗

Synthesis of chiral self-assembling rhombs and their characterization in solution, in the gas phase, and at the liquid-solid interface.

Chiral, enantiopure metallo-supramolecular rhombs self-assemble in solution through coordination of bis-pyridyl-substituted ligands with (en)M(NO3)2 (en = ethylenediamine, M = Pd(II), Pt(II)). Characterization by NMR and CD spectroscopy in solution and by ESI-FT-ICR mass spectrometry in the gas phase suggests that an equilibrium exists in water/methanol of a major 2:2 complex and a minor 3:3 complex of ligands and metal corners. In the gas phase, doubly charged 2:2 complexes fragment into two identical singly charged halves followed by metal-mediated C-H and C-C bond activation reactions within the ethylenediamine ligands. Electrochemical scanning tunneling microscopy (EC-STM) provides in situ imaging of the complexes even with submolecular resolution. Flat-lying rhombs are deposited under potential control from an aqueous electrolyte on a Cu(100) electrode surface precovered by a tetragonal pattern of chloride anions from the supporting electrolyte. Chirality induces the formation of only one domain orientation. Density functional calculations help to interpret the STM images.

Journal Article↗

Nitrogen fixation under mild ambient conditions: part I--the initial dissociation/association step at molybdenum triamidoamine complexes.

In several recent studies Schrock and collaborators demonstrated for the first time how molecular dinitrogen can be catalytically transformed under mild and ambient conditions to ammonia by a molybdenum triamidoamine complex. In this work, we investigate the geometrical and electronic structures involved in this process of dinitrogen activation with quantum chemical methods. Density functional theory (DFT) has been employed to calculate the coordination energies of ammonia and dinitrogen relevant for the dissociation/association step in which ammonia is substituted by dinitrogen. In the DFT calculations the triamidoamine chelate ligand has been modeled by a systematic hierarchy of increasingly complex substituents at the amide nitrogen atoms. The most complex ligand considered is an experimentally known ligand with an HMT = 3,5-(2,4,6-Me3C6H2)2C6H3 substituent. Several assumptions by Schrock and collaborators on key reaction steps are confirmed by our calculations. Additional information is provided on many species not yet observed experimentally. Particular attention is paid to the role of the charge of the complexes. The investigation demonstrates that dinitrogen coordination is enhanced for the negatively charged metal fragment, that is, coordination is more favorable for the anionic metal fragment than for the neutral species. Coordination of N2 is least favorable for the cationic metal fragment. Furthermore, ammonia abstraction from the cationic complex is energetically unfavorable, while NH3 abstraction is less difficult from the neutral and easily feasible from the anionic low-spin complex.

Journal Article↗

Relativistic DMRG calculations on the curve crossing of cesium hydride.

Over the past few years, it has been shown in various studies on small molecules with only a few electrons that the density-matrix renormalization group (DMRG) method converges to results close to the full configuration-interaction limit for the total electronic energy. In order to test the capabilities of the method for molecules with complex electronic structures, we performed a study on the potential-energy curves of the ground state and the first excited state of 1sigma+ symmetry of the cesium hydride molecule. For cesium relativistic effects cannot be neglected, therefore we have used the generalized arbitrary-order Douglas-Kroll-Hess protocol up to tenth order, which allows for a complete decoupling of the Dirac Hamiltonian. Scalar-relativistic effects are thus fully incorporated in the calculations. The potential curves of the cesium hydride molecule feature an avoided crossing between the ground state and the first excited state, which is shown to be very well described by the DMRG method. Compared to multireference configuration-interaction results, the potential curves hardly differ in shape, for both the ground state and the excited state, but the total energies from the DMRG calculations are in general consistently lower. However, the DMRG energies are as accurate as corresponding coupled cluster energies at the equilibrium distance, but convergence to the full configuration-interaction limit is not achieved.

Journal Article↗

Theoretical study on the spin-state energy splittings and local spin in cationic [Re]-Cn-[Re] complexes.

Total spin-state energy splittings are calculated for mono- and dications of the formula {[Re]-Cn-[Re]}z+ where [Re] = eta5-(C5Me5)Re(NO)(PPh3). Cn is an even-numbered carbon chain with n ranging from 4 to 20, and z is 1 or 2. These complexes are experimentally known, and their potential role as molecular electronic devices initiated this work. We have considered the different total spin states monocation/doublet, monocation/quartet, dication/singlet, and dication/triplet. Data obtained for two density functionals BP86 and B3LYP were compared to verify the internal consistency of the results. In both ionization states, the low-spin state is the ground state, but the spin-state splittings decrease as the chain gets longer. For the dications, the splitting reaches a nearly constant value of about 10 kJ/mol with BP86 and about 4 kJ/mol with B3LYP when there are at least 14 carbon atoms in the chain, whereas for the monocations, no constant value appears to be reached asymptotically, not even if 20 carbon atoms are in the chain. For monocations, the splittings range from 138 kJ/mol (n = 4) to 68 kJ/mol (n = 20) with BP86 and from 134 kJ/mol (n = 4) to 73 kJ/mol (n = 20) with B3LYP and are thus considerably higher than those of the dications. The spin-state splittings are qualitatively mirrored by the energy splitting between the highest-occupied molecular orbital with beta spin (HOMObeta) and the lowest-unoccupied molecular orbital with alpha spin (LUMOalpha) as obtained in the low-spin state. Furthermore, the HOMOalpha-LUMOalpha gaps decrease as the carbon chain lengthens. In addition, the local distribution of the ŝz expectation value is analyzed for the monocation/doublet, the monocation/quartet, and the dication/triplet state using a modified Löwdin partitioning scheme. In the monocation/doublet and the dication/triplet state, the electron spin is distributed mainly on the metal centers and slightly delocalized onto the carbon chain. In the monocation/quartet state for chain lengths of more than 8 carbon atoms, the electron spin is mainly localized on selected atoms of the chain and not on the metal centers. In all cases, the spin delocalization onto the chain increases as the chain gets longer.

Journal Article↗

Basis set and density functional dependence of vibrational Raman optical activity calculations.

We present an extensive investigation of the dependence of the scattering intensity difference of right and left circularly polarized light observed in vibrational Raman optical activity (VROA) on the choice of basis set and exchange-correlation functional. These dependencies are investigated for five molecules for which accurate experimental data are available: (S)-methyloxirane, (R)-epichlorhydrin, (S)-glycidol, (M)-spiro[2,2]pentane-1,4-diene, and (M)-sigma-[4]-helicene. Calculations are presented using the SVWN exchange-correlation functional (LDA), the BLYP exchange-correlation functional, and the B3LYP hybrid functional, using six different basis sets: the cc-pVDZ, cc-pVTZ, aug-cc-pVDZ, aug-cc-pVTZ, Sadlej's polarized basis set, and a minimal VROA basis set recently proposed by Zuber and Hug. It is demonstrated that results from pure gradient-corrected and hybrid functionals are comparable and that the aug-cc-pVDZ and aug-cc-pVTZ basis sets yield similar results. Furthermore, the combination of the small basis set by Zuber and Hug with an accurate force field represents the best compromise between computational accuracy and computational efficiency.

Journal Article↗

Understanding the template preorganization step of an artificial arginine receptor.

A biomimetic complex which mimics the arginine-phosphonate diester interaction of the arginine fork is investigated with respect to structure and energetics of stable configurations. Within this work, we provide knowledge on local minima of the isolated system obtained from first-principles calculations. Non-negligible solvation effects are studied in a microsolvation approach. The interactions which govern the structural patterns of molecular recognition in this tweezer-guest complex can be significantly modulated by the action of hydrogen bond accepting and donating solvent molecules, such as dimethyl sulfoxide or water, which were present in experimental investigations on this system. Different tweezer-guest structures are evaluated with respect to their temperature-dependent thermodynamical properties as products of the first association reaction step of the bisphosphonate tweezer template and the guanidinium moiety.

Arginine↗

Comparison of density functionals for differences between the high- (5T2g) and low- (1A1g) spin states of iron(II) compounds. IV. Results for the ferrous complexes [Fe(L)('NHS4')].

Previous work testing density functionals for use in calculating high-spin-low-spin energy differences, deltaE(HL), for iron(II) spin-crossover transitions has tended to conclude that only properly reparametrized hybrid functionals can predict deltaE(HL) since it seems to depend critically on a correct description of the electron pairing energy governed by the exchange term. Exceptions to this rule are the previous three papers (I, II, and III in the present series of papers) where it was found that modern generalized gradient approximations (GGAs) and meta-GGAs could do as well as hybrid functionals, if not better, for this type of problem. In the present paper, we extend these previous studies to five more molecules which are too large to treat with high-quality ab initio calculations, namely, the series [Fe(L)('NHS(4)')], where NHS(4)=2.2'-bis(2-mercaptophenylthio)diethylamine dianion, and L=NH(3), N(2)H(4), PMe(3), CO, and NO(+). Since we know of no reliable experimental estimate of deltaE(HL), we content ourselves with a comparison against the experimentally determined ground-state spin symmetry including, in so far as possible, finite-temperature effects. Together with the results of Papers I, II, and III, this paper provides a test of a large number of functionals against the high-spin/low-spin properties of a diverse set of Fe(II) compounds, making it possible to draw some particularly interesting conclusions. Trends among different classes of functionals are discussed and it is pointed out that there is at least one functional, namely, the OLYP generalized gradient approximation, which is able to give a reasonably good description of the delicate spin energetics of Fe(II) coordination compounds without resorting to hybrid functionals which require the relatively more expensive calculation of a Hartree-Fock-type exchange term.

Journal Article↗