The design of organic catalysis for epoxidation by hydrogen peroxide

The potential of various organic species to catalyze epoxidation of ethene by hydrogen peroxide is explored with B3LYP/6-31G* DFT calculations. Electronic Supplementary Materials Supplementary material is available for this article at http://dx.doi.org/10.1007/s00894-005-0044-4.Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Theoretical 49Ti NMR chemical shifts

⁴⁹Ti chemical shifts for a total of 20 titanium complexes are reported, and several levels of theory are evaluated in order to identify a reliable approach for the calculation of titanium NMR data. The popular B3LYP/6–31G(d)//B3LYP/6–31G(d) proves to give very good agreement with experimental data over a range from 1,400 to −1,300 ppm. The MP2/6–31G(d)//MP2/6–31G(d) level computes even smaller average deviations but fails for TiI₄. This behavior together with its huge demand for computational resources requires careful handling of this theoretical level. In addition, NMR data for five titanium fulvene (or related) complexes are given. Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Umpolung catalysts: comparative assessments on reactivities

Umpolung catalysis is studied by a sequence of model reactions (CPCM in THF, B3LYP/6-31G*) with different aldehydes and catalysts. We involved addition of the catalyst to the aldehyde and 1,2-H-migration to form a carbanionic d¹-species, which is the crucial intermediate according to the Lapworth- and Breslow-mechanisms. Cyanide, N-methylthiazol-2-ylidene, and a glycol-based phosphite perform as umpolung catalysts, formaldehyde, acetaldehyde, benzaldehyde, and acrolein are substrates in this study. In these aldehyde substrates, alkyl-substitution disfavors but π-conjugation favors formation of the carbanionic d¹-intermediate. The nucleophilic carbene, N-methylthiazol-2-ylidene, is the strongest umpolung catalyst, while the phosphite is about as active as cyanide. Figure Transitions structure for the umpolung of formaldehyde with a glycol phosphinite catalyst Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Does the donor–acceptor concept work for designing synthetic metals?

Homopolymers of quinoxaline (QX), benzothiadiazole (BT), benzobisthiadiazole (BBT), thienopyrazine (TP), thienothiadiazole (TT), and thienopyrazinothiadiazole (TTP) and copolymers of these acceptors with thiophene (TH) and pyrrole (PY) were investigated with density functional theory. Theoretical band-gap predictions reproduce experimental data well. For all but six copolymers, band-gap reductions with respect to either homopolymer are obtained. Four of the acceptors, BBT, TP, TT, and TTP, give rise to copolymers with band gaps that are smaller than that of polyacetylene. BBT and TTP copolymers with PY in 1:2 stoichiometry are predicted to be synthetic metals. Band-gap reductions result from upshifts of HOMO energies and much smaller upshifts of LUMO values. The smallest band gaps are predicted with TTP, since changes in LUMO energies upon copolymerization are particularly small. The consequence of the small interactions between LUMO levels of donor and acceptor are vanishingly small conduction bandwidths. Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

An experimental and theoretical study of the enantioselective deprotonation of cyclohexene oxide with isopinocampheyl-based chiral lithium amides

The mechanism of the enantioselective deprotonation of cyclohexene oxide with isopinocampheyl-based chiral lithium amide was studied by quantum chemical calculations. The transition states of eight molecules were fully optimized at the ab initio HF/3-21G and density functional B3LYP/3-21G levels with Gaussian 98. The activation energies were calculated at the B3LYP/6-31+G(3df,2p)//B3LYP/3-21G level. We found the theoretical evaluation to be consistent with the experimental data. At the best case, an enantiomeric excess of up to 95% for (R)-2-scyclohexen-1-ol was achieved with (−)-N, N-diisopinocampheyl lithium amide. Enantioselective deprotonation of cyclohexene oxide Electronic Supplementary Material Supplementary material is available for this article at Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday.     

Boron nitride cages from B12N12 to B36N36: square–hexagon alternants vs boron nitride tubes

The structures and stabilities of square–hexagon alternant boron nitrides (B _x N _x , x=12–36) vs their tube isomers containing octagons, decagons and dodecagons have been computed at the B3LYP density functional level of theory with the correlation-consistent cc-pVDZ basis set of Dunning. It is found that octagonal B₂₀N₂₀ and B₂₄N₂₄ tube structures are more stable than their square–hexagon alternants by 18.6 and 2.4 kcal mol⁻¹, respectively, while the square–hexagon alternants of other cages are more stable. Trends in stability as a function of cluster size are discussed.Figure The octagonal B₂₀N₂₀ and B₂₄N₂₄ tube structures are more stable than their square-hexagon alternant cagesDedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Numerical simulation and graphical analysis of in vitro benign tumor growth: application of single-particle state bosonic matter equation with length scaling

We describe the application of a non-linear single-particle state bosonic condensate equation to simulate multicellular tumor growth by treating it as a coupling of two classical wave equations with real components. With one component representing the amplitude of the cells in their volume growth phase and the other representing the amplitude of the cells in their proliferation or mitosis phase, the two components of the coupled equation feed each other during the time evolution and are coupled together through diffusion and other linear and non-linear terms. The features of quiescent and necrotic cells, which result from poor nutrient diffusion into a tumor, have been found to correspond quite well to experimental data when they are modeled as depending on higher cell density. Classical hallmarks of benign tumor growth, such as the initial rapid growth, followed by a dramatic collapse in the proliferating cell count and a strong re-growth thereafter appear quite encouragingly in the theoretical results. A tool for graphical analysis of the tumor simulation results has been developed to provide morphological information about tumors at various growth stages. The model and the graphical analysis can be extended further to create an effective tool to predict/monitor tumor growth. 1 Screen shot from the graphical analysis tool showing simulation results after ten days: clustering of cells of the tumor (up); cell density profile (down) Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Synthesis, DFT calculations, linear and nonlinear optical properties of binuclear phthalocyanine gallium chloride

The axially substituted binuclear GaCl/GaCl phthalocyanine 1 with an unsymmetrical pattern of substitution has been prepared and its nonlinear optical (NLO) properties determined. The resulting binuclear complex retains approximately the same transition energies of monomeric (RO)₈PcGaCl as far as the linear optical spectrum is concerned, although 1 has a double concentration of central atoms per molecule and an enlarged conjugated ligand. The lack of significant spectral shifts in passing from mononuclear to binuclear complexes has been rationalized theoretically by means of density functional theory calculations. The purpose of the present study is to determine whether binuclearity affects the optical limiting behavior of 1 with respect to monomeric (RO)₈PcGaCl in the NLO regime determined by nanosecond laser pulses. Figure Bis axially substituted binuclear phthalocyanine: synthesis, DFT calculations and NLO properties Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday.     

Periodic trends and easy estimation of relative stabilities in 11-vertex nido-p-block-heteroboranes and -borates

Density functional theory computations were carried out for 11-vertex nido-p-block-hetero(carba)boranes and -borates containing silicon, germanium, tin, arsenic, antimony, sulfur, selenium and tellurium heteroatoms. A set of quantitative values called “estimated energy penalties” was derived by comparing the energies of two reference structures that differ with respect to one structural feature only. These energy penalties behave additively, i.e., they allow us to reproduce the DFT-computed relative stabilities of 11-vertex nido-heteroboranes in general with good accuracy and to predict the thermodynamic stabilities of unknown structures easily. Energy penalties for neighboring heteroatoms (HetHet and HetHet′) decrease down the group and increase along the period (indirectly proportional to covalent radii). Energy penalties for a five- rather than four-coordinate heteroatom, [Het_5k(1) and Het_5k(2)], generally, increase down group 14 but decrease down group 16, while there are mixed trends for group 15 heteroatoms. The sum of HetHet′ energy penalties results in different but easily predictable open-face heteroatom positions in the thermodynamically most stable mixed heterocarbaboranes and -borates with more than two heteroatoms. Figure Correlation of HetHet′ and HetC increments with covalent radii of group 15 heteroatoms Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Bond torsion affects the product distribution in the photoreaction of retinal model chromophores

Ab initio molecular dynamics (MD) calculations have been performed to study the photoisomerization of a 3-double-bond retinal model chromophore, the all-trans-4, 6-dimethylpenta-3, 5-dieniminium cation, and the possible influence of non-planar distortions on the product distribution. In total, 171 trajectories have been generated for four different conformations of the structure, a planar one and three in which the C4–C5 and the C5=C6 bonds were increasingly twisted out of plane. Starting geometries randomly distributed about the equilibrium geometry were generated by zero-point energy sampling; trajectories were calculated using CASSCF-BOMD methodology and were followed until the photoproduct and its configuration could be assigned. For the latter, two different approaches were applied, one involving the CASSCF configuration vectors, the other an analysis of the MD at the first possible hopping event. Isomerization was found to occur almost exclusively about the central C3=C4 double bond in the case of the planar model compound. Twisting the conjugated π-system shifts the isomerization site from the central double bond to the terminal C5=C6 double bond. With both the C4–C5 and the C5=C6 bonds twisted by 20°, about 35% of the trajectories lead to the configurationally inverted 5-cis product. The results are discussed with reference to the highly selective and efficient photo-induced isomerization of the retinal chromophore in rhodopsin. Figure Product distribution in the MD simulations of models 1–4. The percentage of trajectories that lead to either C3=C4 or C5=C6 rotation is given beside the bar graphs. The green and the red portions of the bars represent the productive and the unproductive events, respectively, with respect to that particular rotation estimated from the torsion angles at the first close approach of the energy surfaces Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday.     

Theoretical study of AlH(2+)n (n=1-7) dications

The structures and stability of 1–7 dications were calculated at the ab initio MP2/aug-cc-pVTZ level of theory. The dications AlH²⁺ 1 and 2 were characterized to be unstable thermodynamically. However, these and the stable dications, 3–7 have considerable kinetic barriers for deprotonation. Each of the structures 3–7 contains one or more two-electron three-center (2e–3c) bonds. Aluminum atoms of these dications carry most of the positive charges, as indicated by NBO charge calculations.Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Probing the nature of hydrogen bonds in DNA base pairs

Energy decomposition analyses based on the block-localized wave-function (BLW-ED) method are conducted to explore the nature of the hydrogen bonds in DNA base pairs in terms of deformation, Heitler–London, polarization, electron-transfer and dispersion-energy terms, where the Heitler–London energy term is composed of electrostatic and Pauli-exchange interactions. A modest electron-transfer effect is found in the Watson–Crick adenine–thymine (AT), guanine–cytosine (GC) and Hoogsteen adenine-thymine (H-AT) pairs, confirming the weak covalence in the hydrogen bonds. The electrostatic attraction and polarization effects account for most of the binding energies, particularly in the GC pair. Both theoretical and experimental data show that the GC pair has a binding energy (−25.4 kcal mol⁻¹ at the MP2/6-31G** level) twice that of the AT (−12.4 kcal mol⁻¹) and H-AT (−12.8 kcal mol⁻¹) pairs, compared with three conventional N-H···O(N) hydrogen bonds in the GC pair and two in the AT or H-AT pair. Although the remarkably strong binding between the guanine and cytosine bases benefits from the opposite orientations of the dipole moments in these two bases assisted by the π-electron delocalization from the amine groups to the carbonyl groups, model calculations demonstrate that π-resonance has very limited influence on the covalence of the hydrogen bonds. Thus, the often adopted terminology “resonance-assisted hydrogen bonding (RHAB)” may be replaced with “resonance-assisted binding” which highlights the electrostatic rather than electron-transfer nature of the enhanced stabilization, as hydrogen bonds are usually regarded as weak covalent bonds. Figure Electron density difference (EDD) maps for the GC pair: a shows the polarization effect (isodensity 1.2×10⁻³ a.u.); b shows the charge transfer effect (isodensity 2×10⁻⁴ a.u.) Dedicated to Professor Paul von Ragué Schleyer on the occasion of his 75th birthday     

Chemical carcinogens in non-enzymatic cytosine deamination: 3-isocyanatoacrylonitrile

Uracil has long been known as the main product of nitrosative cytosine deamination in aqueous solution. Recent mechanistic studies of cytosinediazonium ion suggest that the cation formed by its dediazoniation can ring-open to N-protonated (Z,s-cis)-3-isocyanatoacrylonitrile 7. Stereochemical preferences are discussed of the 3-isocyanatoacrylonitriles (Z,s-cis)-10, (E,s-cis)-11, (Z,s-trans)-12, and (E,s-trans)-13. The electronic structures of 7 and 10–13 have been analyzed and a rationale is provided for the thermodynamic preference for (Z,s-cis)-10. It is shown that s-cis/s-trans-interconversion occurs via C−N rotation–inversion paths with barriers below 3 kcal mol⁻¹. The proton affinities of 3-isocyanatoacrylonitrile 10 and water are nearly identical and, thus, 3-isocyanatoacrylonitriles can and should be formed in aqueous media from 7 along with 3-aminoacrylonitriles 9. The results highlight the relevance of the chemistry of 3-isocyanatoacrylonitriles for the understanding of the chemical toxicology of nitrosation of the nucleobase cytosine. Electronic Supplementary Material Supplementary Material is available for this article at Dedicated to professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Modeling study of the influences of the aromatic transitions and the local environment on the far-UV rotational strengths in TEM-1 β-lactamase

Rotational strengths in the far-UV of TEM-1 β-lactamase have been investigated with two theoretical models based on the matrix method. The first model excludes, and a second includes, effects of the local electrostatic interactions on the chromophore energies. Special attention is given to the contributions of the aromatic side-chain chromophores, and the mechanisms of generation of rotational strengths are analyzed. The sensitivity of the computational models with respect to the structural changes of the protein are discussed. Figure Structure of TEM-1 β-lactamase. Both domains—α and αβ, the secondary structural elements and the aromatic and disulfide chromophores are shown     

Theoretical study of the hydrogen-bonded complexes serotonin–water/hydrogen peroxide

Eight H-bonded complexes between serotonin (5-hydroxy-tryptamine) and water/hydrogen peroxide were studied at the B3LYP and HF levels of theory, using the 6-31+G(d) basis set. A thermodynamic analysis was performed in order to find the most stable complex. The calculated bonding parameters showed that the most stable H-bonded complex is formed between serotonin and hydrogen peroxide by means of the intermolecular H-bond –H₂N...H–OOH. Fig. a Theoretical study of the hydrogen-bonded supersystems serotonin-water/hydrogen peroxide     

A hydrophobic similarity analysis of solvation effects on nucleic acid bases

We investigate the changes in the solvation properties of the natural nucleic acid bases due to the formation of the canonical Watson–Crick hydrogen-bonded complexes. To this end, the changes in the free energy of solvation of the bases induced upon hydrogen-bonded dimerization are analyzed by means of the hydrophobic similarity index, which relies on the atomic contributions to the free energy of solvation determined by the partitioning method implemented in the framework of the MST continuum model. Such an index is also used to examine the hydrophobic similarity between the canonical nucleic acid bases and a series of highly apolar analogues, which have been designed as potential candidates to expand the genetic alphabet. The ability of these analogues to be incorporated into modified DNA duplexes can be related to the large reduction in the hydrophilicity of the natural bases upon formation of the canonical hydrogen-bonded dimers. The results illustrate the suitability of the hydrophobic similarity index to rationalize the role played by solvation in molecular recognition. Proceedings of “Modeling Interactions in Biomolecules II”, Prague, September 5th–9th, 2005.     

Structures of inclusion complexes of halogenbenzoic acids and α-cyclodextrin based on AM1 calculations

Semi-empirical quantum mechanics calculations using AM1 (Austin Method 1) were carried out for various host-guest combinations of α-cyclodextrin and mono-halogen benzoic acids. The energetically favorable inclusion structures were identified. The AM1 results show that α-cyclodextrin complexes with mono-halogen benzoic acid acids (where the halogen is chlorine, bromide, iodine) as guest compounds are more stable in the “head first” position than in the “tail-first” position for meta and para isomers while ortho mono-halogen benzoic acids complexes with α-cyclodextrin are more stable in “tail-first” position. The calculated structures were found to be in good agreement with those obtained from crystalographic databases.      

A QCAR-approach to materials modeling

Little is known about the relationship between the function and structure of materials. Materials (solids with a function) are complex entities and a better knowledge of the parameters that contribute to function is desirable. Here, we present modeling approaches that correlate chemical composition with function of heterogeneous catalysts. The complete composition space of the mixed oxides of Ni–Cr–Mn and of Ni–Co–Mo–Mn (10% spacing) have been measured for the oxidation of propene to acroleine. The data have been collected, visualized and modeled. Different mathematical approaches such as Support Vector Machines, multilevel B-splines approximation and Kriging have been applied to model this relationship. High-throughput screening data of ternary and quaternary composition spreads are approximated to locate catalysts of high activity within the search space. For quaternary systems, slice plots offer a good tool for visualization of the results. Using these approximation techniques, the composition of the most active catalysts can be predicted. The study documents that distinct relationships between chemical composition and catalytic function exist and can be described by mathematical models. Visualization of a ternary catalyst system and its approximation using slice plots     

A molecular-dynamics simulation study of diffusion of a single model carbonic chain on a graphite (001) surface

Molecular-dynamics simulations have been used to study the diffusion of a short single model carbonic chain on the graphite (001) surface. The calculated diffusion coefficient (D) first increases, then decreases with increasing chain length (N). This abnormal behavior is similar to polymer lateral diffusion at the solid–liquid interface. Furthermore, we have studied the relation between the mean-square gyration radius and N. Figure Log–log plot of the self-diffusion coefficient D versus the chain length N. The error bars are the standard deviation measured in three repeated simulations     

SOMMER: self-organising maps for education and research

SOMMER is a publicly available, Java-based toolbox for training and visualizing two- and three-dimensional unsupervised self-organizing maps (SOMs). Various map topologies are implemented for planar rectangular, toroidal, cubic-surface and spherical projections. The software allows for visualization of the training process, which has been shown to be particularly valuable for teaching purposes. Spread of a spherical self-organizing map (SOM) in a three-dimensional data space