Two complementary molecular energy decomposition schemes: the Mayer and Ziegler-Rauk methods in comparison

In the present paper we discuss and compare two different energy decomposition schemes: Mayer's Hartree-Fock energy decomposition into diatomic and monoatomic contributions [Chem. Phys. Lett. 382, 265 (2003)], and the Ziegler-Rauk dissociation energy decomposition [Inorg. Chem. 18, 1558 (1979)]. The Ziegler-Rauk scheme is based on a separation of a molecule into fragments, while Mayer's scheme can be used in the cases where a fragmentation of the system in clearly separable parts is not possible. In the Mayer scheme, the density of a free atom is deformed to give the one-atom Mulliken density that subsequently interacts to give rise to the diatomic interaction energy. We give a detailed analysis of the diatomic energy contributions in the Mayer scheme and a close look onto the one-atom Mulliken densities. The Mulliken density ρA has a single large maximum around the nuclear position of the atom A, but exhibits slightly negative values in the vicinity of neighboring atoms. The main connecting point between both analysis schemes is the electrostatic energy. Both decomposition schemes utilize the same electrostatic energy expression, but differ in how fragment densities are defined. In the Mayer scheme, the electrostatic component originates from the interaction of the Mulliken densities, while in the Ziegler-Rauk scheme, the undisturbed fragment densities interact. The values of the electrostatic energy resulting from the two schemes differ significantly but typically have the same order of magnitude. Both methods are useful and complementary since Mayer's decomposition focuses on the energy of the finally formed molecule, whereas the Ziegler-Rauk scheme describes the bond formation starting from undeformed fragment densities

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

One- and two-center energy components in the atoms in molecules theory

The energy decomposition scheme proposed in a recent paper has been realized by performing numerical integrations. The sample calculations carried out for some simple molecules show excellent agreement with the chemical picture of molecules, indicating that such an energy decomposition analysis can be useful from the point of view of connecting quantum mechanics with the genuine chemical concepts

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

Modelling energy consumption of relay-enabled MAC protocols in ad hoc networks

To enhance the throughput of ad hoc networks, dual-hop relay-enabled transmission schemes have recently been proposed. Since in ad hoc networks throughput is normally related to their energy consumption, it is important to examine the impact of using relay-enabled transmissions on energy consumption. In this paper, we present an analytical energy consumption model for dual-hop relay-enabled medium access control (MAC) protocols. Based on the recently reported relay-enabled distributed coordination function (rDCF), we have shown the efficacy of the proposed analytical model. This is a generalized model and can be used to predict energy consumption in saturated relay-enabled ad hoc networks via energy decomposition. This is helpful in designing MAC protocols for cooperative communications and it is shown that using a relay results not only in a better throughput but also better energy efficiency.

Decomposition kinetics of iron (III)-diclofenac compound

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Thermogravimetric studies of the decomposition kinetics of four different Hevea rubber clones using Ozawa's approach

Thermal degradation of natural rubber extracted from four different rubber clones of the Hevea brasiliensis species was investigated by thermogravimetry using Ozawa's approach to assess the kinetic parameters of the decomposition process. The results are discussed in terms of the order of reaction, kinetic parameters such as activation energy as a function of conversion degree, and thermal stability for each of the clones.

Insights into the nature of the coupling interactions between uracil corrosion inhibitors and copper : a DFT and molecular dynamics study

In this study, the nature of the coupling interactions between copper and uracil as well as its several derivatives has been systematically investigated employing the atoms in molecules (AIM) theory and energy decomposition analyses. The whole interaction process has been investigated through the analyses of the radial distribution functions of the Cu⋯X (X = S and O) contact on the basis of the ab initio molecular dynamics. No direct relationship between the adsorption strengths and inhibition efficiencies of the inhibitors has been observed. Additionally, the possibility of the methyl-substituted dithiouracil species to act as copper corrosion inhibitors has been tested.

Stabilization of an All-Metal Antiaromatic Molecule (Al4Li4) Using BH and C as Caps

It has been reported by Pati et al. (J. Am. Chem. Soc. 2005, 127, 3496) that coordination with a transition metal can stabilize the “antiaromatic”, all-metal compound Al4Li4. Here, we report that it can also be stabilized by capping with a main group element like C and its isoelectronic species BH. Our calculations of binding energy, nuclear independent chemical shift, energy decomposition analysis, and molecular orbital analysis support the capping-induced stability, reduction of bond length alternation, and increase of aromaticity of these BH/C-capped Al4Li4 systems. The interaction between px and py orbitals of BH/C and the HOMO and LUMO of Al4Li4 is responsible for the stabilization. Our calculations suggest that capping can introduce fluxionality at room temperature.

Contrasting Electronic Requirements for C-H Binding and C-H Activation in d(6) Half-Sandwich Complexes of Rhenium and Tungsten

A computational study of the interaction half-sandwich metal fragments (metal=Re/W, electron count=d(6)), containing linear nitrosyl (NO+), carbon monoxide (CO), trifluorophosphine (PF3), N-heterocyclic carbene (NHC) ligands with alkanes are conducted using density functional theory employing the hybrid meta-GGA functional (M06). Electron deficiency on the metal increases with the ligand in the order NHC < CO < PF3 < NO+. Electron-withdrawing ligands like NO+ lead to more stable alkane complexes than NHC, a strong electron donor. Energy decomposition analysis shows that stabilization is due to orbital interaction involving charge transfer from the alkane to the metal. Reactivity and dynamics of the alkane fragment are facilitated by electron donors on the metal. These results match most of the experimental results known for CO and PF3 complexes. The study suggests activation of alkane in metal complexes to be facile with strong donor ligands like NHC. (C) 2015 Wiley Periodicals, Inc.

A molecular model of braid-like DNA structure

The three-dimensional molecular models of DNA triple helices and triple-stranded brain-like structure were built up by molecular architecture, and their structural features and energy decomposition were examined. The results showed: (i) The base triplet is the element forming braid-like and triple helix DNA; (ii) Under specified conditions, DNA could form the triplet-stranded braid-like structure; (iii) DNA stability of the braid-like structure is less than that of the triple helix structure. (C) 1995 Academic Press Limited.

Influence of surface structures, subsurface carbon and hydrogen, and surface alloying on the activity and selectivity of acetylene hydrogenation on Pd surfaces: A density functional theory study

The selective hydrogenation of acetylene to ethylene on several Pd surfaces (Pd(111), Pd(100), Pd(211), and Pd(211)-defect) and Pd surfaces with subsurface species (carbon and hydrogen) as well as a number of Pd-based alloys (Pd-M/Pd(111) and Pd-M/Pd(211) (M = Cu, Ag and Au)) are investigated using density functional theory calculations to understand both the acetylene hydrogenation activity and the selectivity of ethylene formation. All the hydrogenation barriers are calculated, and the reaction rates on these surfaces are obtained using a two-step model. Pd(211) is found to have the highest activity for acetylene hydrogenation while Pd(100) gives rise to the lowest activity. In addition, more open surfaces result in over-hydrogenation to form ethane, while the close-packed surface (Pd(111)) is the most selective. However, we also find that the presence of subsurface carbon and hydrogen significantly changes the reactivity and selectivity of acetylene toward hydrogenation on Pd surfaces. On forming surface alloys of Pd with Cu, Ag and Au, the selectivity for ethylene is also found to be changed. A new energy decomposition method is used to quantitatively analyze the factors in determining the changes in selectivity. These surface modifiers are found to block low coordination unselective sites, leading to a decreased ethane production. (C) 2013 The Authors. Published by Elsevier Inc. All rights reserved.

Acrolein hydrogenation on Pt(211) and Au(211) surfaces:a density functional theory study

Partial hydrogenation of acrolein, the simplest alpha, beta-unsaturated aldehyde, is not only a model system to understand the selectivity in heterogeneous catalysis, but also technologically an important reaction. In this work, the reaction on Pt(211) and Au(211) surfaces is thoroughly investigated using density functional theory calculations. The formation routes of three partial hydrogenation products, namely propenol, propanal and enol, on both metals are studied. It is found that the pathway to produce enol is kinetically favoured on Pt while on Au the route of forming propenol is preferred. Our calculations also show that the propanal formation follows an indirect pathway on Pt(211). An energy decomposition method to analyze the barrier is utilized to understand the selectivities at Pt(211) and Au(211), which reveals that the interaction energies between the reactants involved in the transition states play a key role in determining the selectivity difference.

Synthèse totale du méthyl ester de la zincophorine et étude DFT du transfert d’hydrogène diastéréosélectif impliquant des intermédiaires radicalaires vicinaux à un ester

Cet ouvrage traite d’une méthodologie pour l’induction de centres stéréogènes au sein des motifs propionates par la chimie des radicaux et de son application à la synthèse totale du méthyl ester de la zincophorine. Cet ionophore, aux propriétés biologiques intéressantes, présente plusieurs défis synthétiques dont une séquence de type polypropionate anti, anti, anti, anti difficilement accessible ainsi qu’un tétrahydropyrane trans trisubstitué. Récemment, l’intérêt renouvelé pour ces composés polyéthers, en tant qu’agents anticancéreux, accentue l’importance de stratégies versatiles permettant l’accès à ces structures ainsi qu’à leurs analogues. Depuis quelques années, notre groupe s’intéresse à la synthèse d’unités propionates acycliques par une séquence réactionnelle contrôlée uniquement par le substrat. La première étape découle d’une aldolisation de Mukaiyama entre un aldéhyde alpha-chiral et un énoxysilane tétrasubstitué portant une liaison carbone-halogène, et où l’issue stéréochimique dépend de la nature monodentate ou bidentate de l’acide de Lewis employé. La seconde réaction-clé implique la formation d’un radical tertiaire, vicinal à un ester, pouvant être réduit diastéréosélectivement en présence d’hydrure d’étain. La première section décrit l’élaboration de motifs tétrahydropyranes trisubstitués et l’induction des centres stéréogènes vicinaux par une réduction radicalaire. Nous avons révélé que l’issue diastéréosélective de la réaction de cyclisation par une iodoéthérification était dictée par le groupement gamma-méthyle des esters alpha,beta-insaturés de départ. Nous avons ensuite démontré que les produits de la réaction radicalaire anti et syn pouvaient être obtenus sélectivement à partir d’un intermédiaire commun, respectivement en prenant appui sur l’effet exocyclique ou endocyclique lors de la réduction. Par une stratégie complémentaire, nous avons révélé que le précurseur radicalaire pouvait également être obtenu par une réaction de cycloéthérification en présence d’un énoxysilane tétrasubstitué. Une étude systématique des substituants du cycle a révélé que certaines relations stéréochimiques conduisaient à une perte de sélectivité au détriment du produit anti. La seconde section concerne l’étude DFT au niveau BHandHLYP/TZVP des intermédiaires radicalaires impliqués lors du transfert d’hydrogène. Par une étude de décomposition de l’énergie d’activation, nous avons été en mesure de rationaliser l’issue diastéréosélective de la réaction sur la base des énergies de distorsion (∆Ed‡) et d’interaction (∆Eint‡) requises pour accéder à la paire d’états de transition pro-anti et pro-syn. De plus, nous avons démontré qu’une analyse NBO permettait de relativiser l’impact des interactions stéréoélectroniques. Par la suite, l’évaluation des intermédiaires radicalaires borinates et aluminates nous a permis de révéler que l’encombrement stérique de la chaîne propionate était la cause principale de la formation sélective des produits anti ou syn lors d’une réduction en présence d’un acide de Lewis. La dernière section décrit l’élaboration de la séquence polypropionate de la zincophorine, et de ses isomères, à partir du fragment tétrahydropyrane substitué. Au cours de notre étude, nous avons identifié que le nombre de sites de chélation potentiels devait être limité à trois lors de l’aldolisation en condition Cram-chélate. De plus, nous avons démontré que les différents motifs acétates sont accessibles sélectivement par l’utilisation d’un énoxysilane encombré. Par ailleurs, nous avons révélé qu’une même séquence réactionnelle pouvait être employée pour la synthèse du fragment C17–C25 de son analogue naturel CP-78,545, et avons complété la plus récente synthèse totale du méthyl ester de la zincophorine.