PubMed Health⌕ Search

Biomedical subjects

Stefan Grimme

Publications and source records attributed to Stefan Grimme.

At least 19 recordsLinked to original sources
PubMed2006

Time dependent density functional theory calculations for electronic circular dichroism spectra and optical rotations of conformationally flexible chiral donor-acceptor dyad.

Twelve conformations of a chiral donor-acceptor (charge-transfer) dyad and six conformations of its dimer complex were structurally optimized by using the Kohn-Sham density functional theory (BLYP/TZV2P) incorporating a recently developed empirical correction scheme that uses C6/R6 potentials for van der Waals interactions (DFT-D). Subsequent time-dependent DFT calculations with BH-LYP and B3-LYP functionals (with triple-zeta basis set) were performed to obtain theoretical circular dichroism (CD) spectra. The experimental CD spectra obtained independently were properly reproduced by averaging the calculated spectra of individual conformers according to a Boltzmann population derived from single-point SCS-MP2 energies. The optical rotations of the monomer were also calculated by using the same functionals with an aug-cc-pVDZ basis set. Dielectric continuum solvation models (COSMO) applied to correct the relative energies from the isolated molecule calculations resulted in conformer distributions that piled the same or even poorer level of agreement with the experimental CD spectrum. Our results clearly show the advantage of the DFT-D method for the geometry optimization of large systems with donor-acceptor interactions and the TD-DFT/BH-LYP calculations for reproducing the experimental CD spectra. As compared with the calculated optical rotations, the wealthy information embedded in the experimental/calculated CD spectra is requisite for the configurational and/or conformational analyses of relatively large and flexible chiral organic molecules in solution.

Journal Article↗
PubMed2006

Semiempirical GGA-type density functional constructed with a long-range dispersion correction.

A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.

Journal Article↗
PubMed2006

Experimental and theoretical characterization of the valence isomerization of Bi-2H-azirin-2-yls to diazabenzenes.

3,4-diazidocyclobutenes 16 were prepared from the corresponding dihalides. Some of these diazides, such as parent compound 16 d and phenyl-substituted derivatives 16 c,f, underwent spontaneous stereoselective electrocyclic ring opening below room temperature, whereas the tetraalkyl derivatives of 16 had to be heated to force the same reaction. In most cases, the resulting 1,4-diazidobuta-1,3-dienes 8 were isolated to study their photochemical transformation into bi-2H-azirin-2-yls 9 via intermediate mono-azirines 17. Except for starting materials with a low number of substituents such as 9 d and 9 f, title compounds 9 underwent a thermal valence isomerization which led exclusively to pyridazines 18 at surprisingly low temperatures. Based on quantum-chemical calculations for the parent bi-2H-azirinyl-2-yl 9 d at the UB3LYP/6-31+G(d) and MR-MP2/TZV(2df,2p) levels, the valence isomerization process is best explained by simultaneous homolytic cleavage of both C--N single bonds of 9 to generate energetically favorable N,N' diradicals 26, which cyclize to 18. The theoretical studies indicate also that one stereoisomer of 9, namely, the rac compound, should undergo valence isomerization more easily than the other, which is in conformity with different rates of these rearrangement reactions found experimentally. For the tetramethyl-bi-2H-azirin-2-yls 9 g, which are better models for the experimentally studied compounds, simultaneous homolytic cleavage of both C--N single bonds is also predicted by the calculations, although the intermediate diradicals 26 g are significantly higher in energy than those of the parent system 9 d.

Journal Article↗
PubMed2006

Metallosupramolecular chemistry with bis(benzene-o-dithiolato) ligands.

The bis(benzene-o-dithiol) ligands H(4)-1, H(4)-2, and H(4)-3 react with [Ti(OC(2)H(5))(4)] to give dinuclear triple-stranded helicates [Ti(2)L(3)](4)(-) (L = 1(4)(-), 2(4)(-), 3(4)(-)). NMR spectroscopic investigations revealed that the complex anions possess C(3) symmetry in solution. A crystal structure analysis for (PNP)(4)[Ti(2)(2)(3)] ((PNP)(4)[14]) confirmed the C(3) symmetry for the complex anion in the solid state. The complex anion in Li(PNP)(3)[Ti(2)(1)(3)] (Li(PNP)(3)[13]) does not exhibit C(3) symmetry in the solid state due to the formation of polymeric chains of lithium bridged complex anions. Complexes [13](4)(-) and [14](4)(-) were obtained as racemic mixtures of the Delta,Delta and Lambda,Lambda isomers. In contrast to that, complex (PNP)(4)[Ti(2)(3)(3)] ((PNP)(4)[15]) with the enantiomerically pure chiral ligand 3(4)(-) shows a strong Cotton effect in the CD spectrum, indicating that the chirality of the ligands leads to the formation of chiral metal centers. The o-phenylene diamine bridged bis(benzene-o-dithiol) ligand H(4)-4 reacts with Ti(4+) to give the dinuclear double-stranded complex Li(2)[Ti(2)(4)(2)(mu-OCH(3))(2)] containing two bridging methoxy ligands between the metal centers. The crystal structure analysis and the (1)H NMR spectrum of (Ph(4)As)(2)[Ti(2)(4)(2)(mu-OCH(3))(2)] ((Ph(4)As)(2)[(16]) reveal C(2) symmetry for the anion [Ti(2)(4)(2)(mu-OCH(3))(2)](2)(-). For a comparative study the dicatechol ligand H(4)-5, containing the same o-phenylene diamine bridging group as the bis(benzene-o-dithiol) ligands H(4)-4, was prepared and reacted with [TiO(acac)(2)] to give the dinuclear complex anion [Ti(2)(5)(2)(mu-OCH(3))(2)](2)(-). The molecular structure of (PNP)(2)[Ti(2)(5)(2)(mu-OCH(3))(2)] ((PNP)(2)[17]) contains a complex anion which is similar to [16](2)(-), with the exception that strong N-H...O hydrogen bonds are formed in complex anion [17](2)(-), while N-H...S hydrogen bonds are absent in complex anion [16](2)(-).

Journal Article↗
PubMed2006

Towards chemical accuracy for the thermodynamics of large molecules: new hybrid density functionals including non-local correlation effects.

Two hybrid density functionals that include a second-order perturbation correction for non-local correlation effects are tested for the full G3/05 test set. Very large AO basis sets including core-polarization/correlation functions have been employed that yield for the first time results quite close to the basis set limit for this set. The B2-PLYP functional and the new mPW2-PLYP approach with a modified exchange part give by far the lowest MAD over the whole G3/05 set ever reported for a DFT method (2.5 and 2.1 kcal mol(-1), respectively). The big improvement compared to common density functionals is further demonstrated by the reduction of the maximum and minimum errors (outliers) and by much smaller errors for complicated molecular systems.

Chemical Phenomena↗
PubMed2006

A theoretical study of the chiroptical properties of molecules with isotopically engendered chirality.

There has been a considerable interest in the chiroptical properties of molecules whose chirality is exclusively due to an isotopic substitution and numerous examples for the electronic circular dichroism (CD) spectra of isotopically chiral systems have been reported in literature. Four different explanations have been proposed for the mechanism as to how the isotopic substitution induces a chiral perturbation of the otherwise achiral electronic wave function; however, up to now no conclusive answer has been given about the dominating effect responsible for the experimental observations. In this study we will present, for the first time, fully quantum-mechanical calculations of the CD spectra of three different molecular systems with isotopically engendered chirality. As examples, we consider the spectra of organic molecules with ketone and alpha-diketone carbonyl and diene chromophores. The effect of vibronic couplings for the reorientation of the electric and magnetic transition dipole moments is taken into account within the Herzberg-Teller approximation. The ground and excited state geometries and vibrational normal modes are obtained with (time-dependent) density functional theory [(TD)DFT], while the vibronic coupling effects are calculated at the TDDFT and density functional theory/multireference configuration interaction (DFT/MRCI) levels of theory. Generally, the band shapes of the experimental CD spectra are reproduced very well, and also the absolute CD intensities from the simulations are of the right order of magnitude. The sign and the intensity of the CD band are determined by a delicate balance of the contributions of a large number of individual vibronic transitions, and it is found that the vibrational normal modes with a large displacement are dominant. The separation of the calculated CD spectrum into the different contributions due to the overlap of the in-plane and out-of-plane components (regarding the symmetry plane of the unsubstituted molecule) of the electric and magnetic transition dipole moments yields information about the influence of the vibronic coupling effects for the reorientation of the corresponding transition dipole moments. In conclusion, the calculations clearly show that vibronic effects are responsible or at least dominant for the chiroptical properties of isotopically chiral organic molecules.

Journal Article↗
PubMed2006

Consistent theoretical description of 1,3-dipolar cycloaddition reactions.

The cycloaddition reactions of 18 1,3-dipolar molecules to ethylene and acetylene have been reinvestigated by quantum chemical methods that are based on a second-order perturbation treatment of electron correlation. It is found that SCS-MP2 and the new perturbative B2-PLYP density functional provide accurate reaction barriers and outperform MP2 as well as standard density functionals such as B3-LYP. The new second-order based methods have the additional advantage that they perform better with increasing quality of the one-particle space, as is desired for a good quantum chemical method. The errors for the reaction enthalpies are in general larger than for the barriers when compared to CBS-QB3 literature values, which is related to strong changes in the electronic structures, but the deviations are again smaller than with MP2 or B3-LYP and are also more systematic. The results of a detailed basis set study suggest that properly polarized triple-zeta AO basis sets represent a good compromise between accuracy and computational speed. The combination of very inaccurate density functionals with small (double-zeta) basis sets, which yields good results for the initial part of the reactions due to error compensation, is not recommended.

Letter↗
PubMed2006

Oligoethylene chains terminated by ferrocenyl end groups: synthesis, structural properties, and two-dimensional self-assembly on surfaces.

A series of unsaturated long-chain-bridged diferrocenes Fc-(CH2)n-CH=CH-(CH2)n-Fc (4 a-e) was synthesized by means of olefin metathesis. Subsequent catalytic hydrogenation furnished the saturated alpha,omega-bis-ferrrocenyl oligoethylene products Fc-(CH2)m-Fc (5). Members of both series formed highly ordered laminar structures at the highly oriented pyrolytic graphite (HOPG) solid/liquid interface or on the Ag(110) surface, which were characterized by STM. Details of the structural features of these ordered physisorbed surface assemblies of 4 and 5 were analyzed by comparison with DFT calculations on model systems and with the characteristic packing modes of these systems in the crystal.

Journal Article↗
PubMed2006

New insights into the geometry of resorc[4]arenes: solvent-mediated supramolecular conformational and chiroptical control.

The conformations of inherently chiral resorc[4]arenes were studied by circular dichroism (CD) spectroscopy. Whereas in aprotic solvents the crown conformation (C4) is preferred, protic solvents favor the boat conformation (C2). As a result of electronic coupling of the lowest L(b) state of the resorcinol unit in the resorc[4]arene, the CD spectra show a strong dependence on the conformation of the macrocycle. For the first time the solvent dependence of the CD spectra was qualitatively analyzed and simulated by using theoretical methods. We have thus demonstrated not only that the conformation of the calixarene is dramatically manipulated by the solvent but also that the joint use of chiroptical measurements and theoretical calculations is a powerful and versatile tool for elucidating structural variations in supramolecular chemistry.

Journal Article↗
PubMed2006

Semiempirical hybrid density functional with perturbative second-order correlation.

A new hybrid density functional for general chemistry applications is proposed. It is based on a mixing of standard generalized gradient approximations (GGAs) for exchange by Becke (B) and for correlation by Lee, Yang, and Parr (LYP) with Hartree-Fock (HF) exchange and a perturbative second-order correlation part (PT2) that is obtained from the Kohn-Sham (GGA) orbitals and eigenvalues. This virtual orbital-dependent functional contains only two global parameters that describe the mixture of HF and GGA exchange (a(x)) and of the PT2 and GGA correlation (c), respectively. The parameters are obtained in a least-squares-fit procedure to the G297 set of heat of formations. Opposed to conventional hybrid functionals, the optimum a(x) is found to be quite large (53% with c=27%) which at least in part explains the success for many problematic molecular systems compared to conventional approaches. The performance of the new functional termed B2-PLYP is assessed by the G297 standard benchmark set, a second test suite of atoms, molecules, and reactions that are considered as electronically very difficult (including transition-metal compounds, weakly bonded complexes, and reaction barriers) and comparisons with other hybrid functionals of GGA and meta-GGA types. According to many realistic tests, B2-PLYP can be regarded as the best general purpose density functional for molecules (e.g., a mean absolute deviation for the two test sets of only 1.8 and 3.2 kcal/mol compared to about 3 and 5 kcal/mol, respectively, for the best other density functionals). Very importantly, also the maximum and minimum errors (outliers) are strongly reduced (by about 10-20 kcal/mol). Furthermore, very good results are obtained for transition state barriers but unlike previous attempts at such a good description, this definitely comes not at the expense of equilibrium properties. Preliminary calculations of the equilibrium bond lengths and harmonic vibrational frequencies for diatomic molecules and transition-metal complexes also show very promising results. The uniformity with which B2-PLYP improves for a wide range of chemical systems emphasizes the need of (virtual) orbital-dependent terms that describe nonlocal electron correlation in accurate exchange-correlation functionals. From a practical point of view, the new functional seems to be very robust and it is thus suggested as an efficient quantum chemical method of general purpose.

Journal Article↗
PubMed2005

Weak intermolecular interactions calculated with diffusion Monte Carlo.

The performance of fixed node diffusion Monte Carlo (FNDMC) for weakly interacting molecules is investigated. The effect of Gaussian basis sets on the asymptotic description of the molecular orbitals which is crucial for a successful importance sampling is analyzed for the example of the hydrogen atom. We find that accurate reference binding energies of the water, the ammonia, and the T-shaped as well as the parallel-displaced benzene dimer are correctly reproduced by FNDMC. The binding energies for the benzene dimers are -3.00(0.38) and -3.58(0.38) kcal/mol, respectively. The description of the methane dimer which has the smallest binding energy and a quite large intermolecular distance requires a more flexible basis set of diffuse quadruple-zeta quality in order to prevent sampling errors.

Journal Article↗
PubMed2005

Van der Waals complexes of polar aromatic molecules: unexpected structures for dimers of azulene.

Full geometry optimizations at the dispersion corrected DFT-BLYP/TZV2P level of theory have been performed for dimers of azulene that may serve as a model system for the van der Waals complexes of polar pi systems. The structures and binding energies for 11 dimers are investigated in detail. The DFT-D interaction energies have been successfully checked against results from the accurate SCS-MP2/aug-cc-pVTZ approach. Out of the nine investigated stacked complexes, eight have binding energies larger than 7.4 kcal/mol (SCS-MP2) that exceed the value of 7.1 kcal/mol for the best naphthalene dimer. T-shaped arrangements (CH...pi) are significantly less stable. Two out of the three best structures have an antiparallel alignment of the monomer dipole moments in the complex, although the best ones with a parallel orientation are only about 0.5 kcal/mol less strongly bound which points to a minor importance of dipole-dipole interactions to binding. Quite surprisingly, the energetically lowest structure (DeltaE = -9.2 kcal/mol) corresponds to a situation where the two seven-membered rings are located almost on top of each other (7-7) and the long molecular axes are rotated against each other by 130 degrees. The 7-7 structural motif is found also in other energetically low-lying structures, and the expected 5-7 (two-side) arrangement is less strongly bound by about 2 kcal/mol. This can be explained by the electrostatic potential of azulene that only partially reflects the charge separation according to the common 4n + 2 pi electron rule. General rules for predicting stable van der Waals complexes of polar pi systems are discussed.

Azulenes↗
PubMed2005

Catalytic enantioselective reactions driven by photoinduced electron transfer.

Photoinduced electron transfer is an essential step in the conversion of solar energy into chemical energy in photosystems I and II (ref. 1), and is also frequently used by chemists to build complex molecules from simple precursors. During this process, light absorption generates molecules in excited electronic states that are susceptible to accepting or donating electrons. But although the excited states are straightforward to generate, their short lifetimes makes it challenging to control electron transfer and subsequent product formation-particularly if enantiopure products are desired. Control strategies developed so far use hydrogen bonding, to embed photochemical substrates in chiral environments and to render photochemical reactions enantioselective through the use of rigid chiral complexing agents. To go beyond such stoichiometric chiral information transmission, catalytic turnover is required. Here we present a catalytic photoinduced electron transfer reaction that proceeds with considerable turnover and high enantioselectivity. By using an electron accepting chiral organocatalyst that enforces a chiral environment on the substrate through hydrogen bonding, we obtain the product in significant enantiomeric excess (up to 70%) and in yields reaching 64%. This performance suggests that photochemical routes to chiral compounds may find use in general asymmetric synthesis.

Journal Article↗
PubMed2005

Van der Waals interactions in aromatic systems: structure and energetics of dimers and trimers of pyridine.

Full geometry optimizations at the dispersion-corrected DFT-BLYP level of theory were carried out for dimers and trimers of pyridine. The DFT-D interaction energies were checked against results from single-point SCS-MP2/aug-cc-pVTZ calculations. Three stacked structures and a planar H-bonded dimer were found to be very close in energy (interaction energies in the range from -3.4 to -4.0 kcal mol(-1)). Two T-shaped geometries are higher lying, by about 1 kcal mol(-1), which is explained by the more favorable electrostatic interactions in the stacked and H-bonded arrangements. The DFT-D approach has proved to be a reliable and efficient tool to explore the conformational space of aromatic van der Waals complexes and furthermore provides interaction energies with errors of less than 10-20 % of DeltaE. Comparisons with previous results obtained by using only partially optimized model geometries strongly indicate that unconstrained optimizations are mandatory in such weakly bonded low-symmetry systems.

Hydrophobic and Hydrophilic Interactions↗