Room temperature metathesis of aryl isocyanates and aromatic aldehydes catalyzed by group(IV) metal alkoxides: An experimental and computational study

Aromatic aldehydes and aryl isocyanates do not react at room temperature. However, we have shown for the first time that in the presence of catalytic amounts of group(IV) n-butoxide, they undergo metathesis at room temperature to produce imines with the extrusion of carbon dioxide. The mechanism of action has been investigated by a study of stoichiometric reactions. The insertion of aryl isocyanates into the metal n-butoxide occurs very rapidly. Reaction of the insertion product with the aldehyde is responsible for the metathesis. Among the n-butoxides of group(IV) metals, Ti((OBu)-Bu-n)(4) (8aTi) was found to be more efficient than Zr((OBu)-Bu-n)(4) (8aZr) and Hf((OBu)-Bu-n)(4) (8aHf) in carrying out metathesis. The surprisingly large difference in the metathetic activity of these alkoxides has been probed computationally using model complexes Ti(OMe)(4) (8bTi), Zr(OMe)(4) (8bZr) and Hf(OMe)(4) (8bHf) at the B3LYP/LANL2DZ level of theory. These studies indicate that the insertion product formed by Zr and Hf are extremely stable compared to that formed by Ti. This makes subsequent reaction of Zr and Hf complexes unfavorable.

Design, Synthesis, and DNA Binding Properties of Photoisomerizable Azobenzene−Distamycin Conjugates: An Experimental and Computational Study

Here, we present the synthesis, photochemical, and DNA binding properties of three photoisomerizable azobenzene−distamycin conjugates in which two distamycin units were linked via electron-rich alkoxy or electron-withdrawing carboxamido moieties with the azobenzene core. Like parent distamycin A, these molecules also demonstrated AT-specific DNA binding. Duplex DNA binding abilities of these conjugates were found to depend upon the nature and length of the spacer, the location of protonatable residues, and the isomeric state of the conjugate. The changes in the duplex DNA binding efficiency of the individual conjugates in the dark and with their respective photoirradiated forms were examined by circular dichroism, thermal denaturation of DNA, and Hoechst displacement assay with poly[d(A-T).d(T-A)] DNA in 150 mM NaCl buffer. Computational structural analyses of the uncomplexed ligands using ab initio HF and MP2 theory and molecular docking studies involving the conjugates with duplex d[(GC(AT)10CG)]2 DNA were performed to rationalize the nature of binding of these conjugates.

Selectivity in Garratt-Braverman Cyclization: An Experimental and Computational Study

Bispropargyl sulfones equipped with aromatic rings of dissimilar nature were synthesized. Under basic conditions, these sulfones isomerized to the bisallenic sulfones, creating a competitive scenario between two alternate Garratt-Braverman (GB) cyclization pathways. The observed product distribution ruled out the involvement of any ionic intermediate and supported the diradical mechanism with greater involvement of the electron-rich aromatic ring via the more nucleophilic radical. DFT-based calculations supported the diradical mechanism along with the observed selectivity.

Diastereospecific Coupling of Imines by Low-Valent Titanium: An Experimental and Computational Study

The reaction of phenylsilane (PhSiH3) and titanium(IV) isopropoxide Ti(OiPr)(4)] generates low-valent titanium alkoxides that reduce and reductively couple imines. The C-C coupling reaction is diastereospecific, with exclusive formation of the (+/-)-isomer. The yield is dependent on the concentration of titanium(IV) isopropoxide used. Using imines with varying electronic demand revealed that the coupling is most efficient for unsubstituted imines. The involvement of a trimeric titanium biradical species is invoked to explain the observed concentration dependence and diastereospecificity in the reaction. Multilayer ONIOM (''our own n-layered integrated molecular orbital and molecular mechanics'') calculations were carried out on the plausible intermediates involved by using the Guassian 09 suite of programs. Formation of the trimeric biradical species leading to the formation of the (+/-)-isomer is more favored than formation of the trimeric biradical species, which leads to the meso-isomer. Similar calculations with dimeric intermediates do not predict the (+/-)-selectivity. Mass spectral analysis of the reaction mixture indicates the presence of a trimeric titanium alkoxide species.

Metathesis of carbon dioxide and phenyl isocyanate catalysed by group(IV) metal alkoxides: An experimental and computational study

The insertion reactions of zirconium(IV) n-butoxide and titanium(IV) n-butoxide with a heterocumulene like carbodiimide, carbon dioxide or phenyl isocyanate are compared. Both give an intermediate which carries out metathesis at elevated temperatures by inserting a second heterocumulene in a head-to-head fashion. The intermediate metallacycle extrudes a new heterocumulene, different from the two that have inserted leading to metathesis. As the reaction is reversible, catalytic metathesis is feasible. In stoichiometric reactions heterocumulene insertion, metathesis and metathesis cum insertion products are observed. However, catalytic amounts of the metal alkoxide primarily led to metathesis products. It is shown that zirconium alkoxides promote catalytic metathesis (isocyanates, carbon dioxide) more efficiently than the corresponding titanium alkoxide. The difference in the metathetic activity of these alkoxides has been explained by a computational study using model complexes Ti(OMe)(4) (1bTi) and Zr(OMe)(4) (1bZr). The computation was carried out at the B3LYP/LANL2DZ level of theory.

A computational study of heat transfer, soot production and transport in an acetylene-air laminar diffusion flame in the field of a line vortex

The present work involves a computational study of soot formation and transport in case of a laminar acetylene diffusion flame perturbed by a co nvecting line vortex. The topology of the soot contours (as in an earlier experimental work [4]) have been investigated. More soot was produced when vortex was introduced from the air si de in comparison to a fuel side vortex. Also the soot topography was more diffused in case of the air side vortex. The computational model was found to be in good agreement with the ex perimental work [4]. The computational simulation enabled a study of the various parameters affecting soot transport. Temperatures were found to be higher in case of air side vortex as compared to a fuel side vortex. In case of the fuel side vortex, abundance of fuel in the vort ex core resulted in stoichiometrically rich combustion in the vortex core, and more discrete so ot topography. Overall soot production too was low. In case of the air side vortex abundan ce of air in the core resulted in higher temperatures and more soot yield. Statistical techniques like probability density fun ction, correlation coefficient and conditional probability function were introduced to explain the transient dependence of soot yield and transport on various parameters like temperature, a cetylene concentration.

A computational study of flame dynamics on soot production and transport for an acetylene-air laminar diffusion flame in the field of a convecting line vortex

The present work involves a computational study of soot formation and transport in case of a laminar acetylene diffusion flame perturbed by a convecting line vortex. The topology of the soot contours (as in an earlier experimental work [4]) have been investigated. More soot was produced when vortex was introduced from the air side in comparison to a fuel side vortex. Also the soot topography was more diffused in case of the air side vortex. The computational model was found to be in good agreement with the experimental work [4]. The computational simulation enabled a study of the various parameters affecting soot transport. Temperatures were found to be higher in case of air side vortex as compared to a fuel side vortex. In case of the fuel side vortex, abundance of fuel in the vort ex core resulted in stoichiometrically rich combustion in the vortex core, and more discrete soot topography. Overall soot production too was low. In case of the air side vortex abundance of air in the core resulted in higher temperatures and more soot yield. Statistical techniques like probability density function, correlation coefficient and conditional probability function were introduced to explain the transient dependence of soot yield and transport on various parameters like temperature, a cetylene concentration.

Charge-transfer-driven inclusion of neutral TCNQ molecules in the galleries of a layered double hydroxide: an experimental and computational study

The layered double hydroxides (LDH) or anionic clays are an important class of ion-exchange materials. They consist of positively charged brucite-like inorganic sheets with charge-compensating exchangeable anions in the interlamellar space. Here we show how neutral TCNQ (7,7,8,8-tetracyanoquinodimethane) molecules can be included within the galleries of an LDH. To do so, we exploit the fact that TCNQ is a good electron acceptor that forms donor acceptor complexes with a variety of donors. The electron donor aniline was intercalated into a Mg-Al LDH as p-aminobenzoate (AB) ions by a conventional ion-exchange reaction. We show here that neutral TCNQ molecules may be driven into the galleries of the layered solid by charge-transfer complex formation with the intercalated p-aminobenzoate anions. We use diffraction and spectroscopic measurements in combination with molecular dynamics simulations and quantum chemical calculations to establish the nature of interactions and arrangement of the charge-transfer complex within the galleries of the layered double hydroxide. Electrostatic interactions between the TCNQ molecules and the anchored AB ions, subsequent to charge transfer, are the driving force for the inclusion of TCNQ molecules in the galleries of the LDH.

Computational Study of the Migration of Rhenium from One Enantioface of an Olefin to the Other Facilitated by (C-H)center dot center dot center dot Re Interactions

The migration of a metal atom in a metal olefin complex from one pi face of the olefin to the opposite pi face has been rarely documented. Gladysz and co-workers showed that such a movement is indeed possible in monosubstituted chiral Re olefin complexes, resulting in diastereomerization. Interestingly, this isomerization occurred without dissociation, and on the basis of kinetic isotope effects, the involvement of a trans C-H bond was indicated. Either oxidative addition or an agostic interaction of the vinylic C-H(D) bond with the metal could account for the experimentally observed kinetic isotope effect. In this study we compute the free energy of activation for the migration of Re from one enantioface of the olefin to the other through various pathways. On the basis of DFT calculations at the B3LYP level we show that a trans (C-H)center dot center dot center dot Re interaction and trans C-H oxidative addition provide a nondissociative path for the diastereomerization. The trans (C-H)center dot center dot center dot Re interaction path is computed to be more favorable by 2.3 kcal mol(-1) than the oxidative addition path. While direct experimental evidence was not able to discount the migration of the metal through the formation of a eta(2)-arene complex (conducted tour mechanism), computational results at the B3LYP level show that it is energetically more expensive. Surprisingly, a similar analysis carried out at the M06 level computes a lower energy path for the conducted tour mechanism and is not consistent with the experimental isotope effects observed. Metal-(C-H) interactions and oxidative additions of the metal into C-H bonds are closely separated in energy and might contribute to unusual fluxional processes such as this diastereomerization.

Computational Study of the Formation of Short Centrosymmetric N-center dot center dot center dot S Supramolecular Synthon and Related Weak Interactions in Crystalline 1,2,4-Triazoles

A comprehensive analysis of the crystal packing and the energetic features of a series of four biologically active molecules belonging to the family of substituted 4-(benzylideneamino)-3-(4-fluoro-3-phenoxyphenyl)-1H-1,2,4-triazole-5-(4 H)-thione derivatives have been performed based on the molecular conformation and the supramolecular packing. This involves the formation of a short centrosymmetric R-2(2)(8) NH...S supramolecular synthon in the solid state, including the presence of CH...S, CH...O, CH...N, CH...F, CH...Cl, CF...FC, CCl...ClC, and CH...pi intermolecular interactions along with pp stacking to evaluate the role of noncovalent interactions in the crystal. The presence of such synthons has a substantial contribution toward the interaction energy (-18 to -20 kcal/mol) as obtained from the PIXEL calculation, wherein the Coulombic and polarization contribution are more significant than the dispersion contribution. The geometrical characteristics of such synthons favor short distance, and the population of related molecules having these geometries is rare as has been obtained from the Cambridge Structural Database (CSD). Furthermore, their interaction energies have been compared with those present in our molecules in the solid state. The topological characteristics of the NH...S supramolecular synthon, in addition to related weak interactions, CH...N, CH...Cl, CF...FC, and CCl...ClC, have been estimated using the quantum theory of atoms in molecules (QTAIM). In addition, an analysis of the Hirshfeld surface and associated fingerprint plots of these four molecules also have provided a platform for the evaluation of the contribution of different atom...atom contacts, which contribute toward the packing of the molecules in solids.

A Computational Study of the Factors Influencing the PVC-Triggering Ability of a Cluster of Early Afterdepolarization-Capable Myocytes

Premature ventricular complexes (PVCs), which are abnormal impulse propagations in cardiac tissue, can develop because of various reasons including early afterdepolarizations (EADs). We show how a cluster of EAD-generating cells (EAD clump) can lead to PVCs in a model of cardiac tissue, and also investigate the factors that assist such clumps in triggering PVCs. In particular, we study, through computer simulations, the effects of the following factors on the PVC-triggering ability of an EAD clump: (1) the repolarization reserve (RR) of the EAD cells; (2) the size of the EAD clump; (3) the coupling strength between the EAD cells in the clump; and (4) the presence of fibroblasts in the EAD clump. We find that, although a low value of RR is necessary to generate EADs and hence PVCs, a very low value of RR leads to low-amplitude EAD oscillations that decay with time and do not lead to PVCs. We demonstrate that a certain threshold size of the EAD clump, or a reduction in the coupling strength between the EAD cells, in the clump, is required to trigger PVCs. We illustrate how randomly distributed inexcitable obstacles, which we use to model collagen deposits, affect PVC-triggering by an EAD clump. We show that the gap-junctional coupling of fibroblasts with myocytes can either assist or impede the PVC-triggering ability of an EAD clump, depending on the resting membrane potential of the fibroblasts and the coupling strength between the myocyte and fibroblasts. We also find that the triggering of PVCs by an EAD clump depends sensitively on factors like the pacing cycle length and the distribution pattern of the fibroblasts.

Enhancement in Bath Mixing and Plume area in a New Degassing Process - A Computational Fluid Dynamic Study

Reaction between the various species in slag and metal phase is usually mass transfer controlled. There have been continuous efforts to increase the reaction efficiency in slag-metal system, especially during decarburization of steel to produce the ultra low carbon steel (ULCS) in secondary steelmaking. It has been found that the surface reaction is a dominant factor in the final stage of decarburization. In the initial stage, the inner site reaction is major factor in the refining process. The mixing of bath affects the later reaction. However, the former reaction (surface reaction) is affected by the plume size area at the top of the metal surface. Therefore, a computational study has been made to understand the fluid dynamics of a new secondary steelmaking process called Revolutionary Degasser Activator (REDA) to study the bath mixing and plume area. REDA process has been considered as it is claimed that this process can reduce the carbon content in steel below 10ppm in a less time than the other existing processes such as RH and Tank degasser. This study shows that both bath mixing and plume area are increased in REDA process facilitating it to give the desired carbon content in less time. Qualitative comments are made on slag-metal reaction system based on this finding.

A numerical study of laminar axisymmetric plumes due to the combined buoyancy of heat and mass diffusion

This paper presents the results of a computational study of laminar axisymmetric plumes generated by the simultaneous diffusion of thermal energy and chemical species. Species concentrations are assumed small. The plume is treated as a boundary layer. Boussinesq approximations are incorporated and the governing conservation equations of mass, momentum, energy and species are suitably non-dimensionalised. These equations are solved using one time-step-forward explicit finite-difference method. Upwind differencing is employed for convective terms. The results thus obtained are explained in terms of the basic physical mechanisms that govern these flows. They show many interesting aspects of the complex interaction of the two buoyant mechanisms.

New insight into the architecture of oxy-anion pocket in unliganded conformation of GAT domains: A MD-simulation study

Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N-terminal ``glutaminase'' (GAT) and C-terminal ``synthetase'' domain. The enzyme is identified as a potential target for anticancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H-bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)-substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non-catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site.

Implicit and multigrid procedures for steady-state computations with upwind algorithms

A computational study for the convergence acceleration of Euler and Navier-Stokes computations with upwind schemes has been conducted in a unified framework. It involves the flux-vector splitting algorithms due to Steger-Warming and Van Leer, the flux-difference splitting algorithms due to Roe and Osher and the hybrid algorithms, AUSM (Advection Upstream Splitting Method) and HUS (Hybrid Upwind Splitting). Implicit time integration with line Gauss-Seidel relaxation and multigrid are among the procedures which have been systematically investigated on an individual as well as cumulative basis. The upwind schemes have been tested in various implicit-explicit operator combinations such that the optimal among them can be determined based on extensive computations for two-dimensional flows in subsonic, transonic, supersonic and hypersonic flow regimes. In this study, the performance of these implicit time-integration procedures has been systematically compared with those corresponding to a multigrid accelerated explicit Runge-Kutta method. It has been demonstrated that a multigrid method employed in conjunction with an implicit time-integration scheme yields distinctly superior convergence as compared to those associated with either of the acceleration procedures provided that effective smoothers, which have been identified in this investigation, are prescribed in the implicit operator.