Laurent Joubert

On the physical role of exchange in the formation of an intramolecular bond path between two electronegative atoms

In this paper, we present a detailed energetic decomposition of intramolecular O···X interactions (X being O, S, or a halogen atom) based on the interacting quantum atoms approach of Pendás and co-workers. The nature of these interactions (repulsive or attractive, more or less electrostatic) is discussed in the framework of Baders atoms in molecules theory, a particular emphasis being put on delocalization (measured by delocalization indexes and in terms of the source function) and on the exchange contributions. Notably, the concept of exchange channels introduced by Pendás and collaborators provides means of rationalizing and predicting the presence of bond critical points, enhancing the physical meaning of bond paths.

Insights into the chemical meanings of the reaction electronic flux

Abstract The negative derivative of the chemical potential with respect to the reaction coordinate is called reaction electronic flux and has recently focused a wide interest to better understand chemical reactions at molecular level. After much consideration, it is now well accepted that positive REF values are associated with spontaneous processes, while negative REF ones translate unspontaneous phenomena. These characteristics of the REF are based on a thermodynamic analogy and have been shown right through computational results. In this paper, we develop two analytical expressions of the REF in both the canonical and the grand canonical ensembles. The connection between both equations is established. They are then analyzed, and some arguments are put forward to support the alleged characteristic of the REF and its ability to properly discriminate spontaneous from unspontaneous phenomena.

Theoretical Study of the Decomposition Reactions in Substituted Nitrobenzenes

The influence of substituent nature and position on the unimolecular decomposition of nitroaromatic compounds was investigated using the density functional theory at a PBE0/6-31+G(d,p) level. As the starting point, the two main reaction paths for the decomposition of nitrobenzene were analyzed: the direct carbon nitrogen dissociation (C6H5 + NO2) and a two step mechanism leading to the formation of phenoxyl and nitro radicals (C6H5O + NO). The dissociation energy of the former reaction was calculated to be 7.5 kcal/mol lower than the activation energy of the second reaction. Then the Gibbs free energies were computed for 15 nitrobenzene derivatives characterized by different substituents (nitro, methyl, amino, carboxylic acid, and hydroxyl) in the ortho, meta, and para positions. In meta position, no significant changes appeared in the reaction energy profiles whereas ortho and para substitutions led to significant deviations in energies on the decomposition mechanisms due to the resonance effect of the nitro group without changing the competition between these mechanisms. In the case of para and meta substitutions, the carbon-nitro bond dissociation energy has been directly related to the Hammett constant as an indicator of the electron donor-acceptor effect of substituents.

On the influence of density functional approximations on some local Bader's atoms-in-molecules properties.

In this article, we assess the ability of various density functionals to predict accurate values for some basic properties of the bond critical points of about 50 small molecules, including the recently proposed reduced gradient variation rates and involving typical ionic and covalent bonds, agostic interactions, and van der Waals complexes. The relation between the computed deviations and the geometric variations are discussed, as well as the topology variations. The possible correlation of these descriptors to atomization energies is considered, and the relevance of an accurate QTAIM analysis for correct descriptions of potential energy surfaces is addressed. Finally, we provide typical margins of error for the evaluation of these quantities and discuss their consequences for computational applications.

The prediction of energies and geometries of hydrogen bonded DNA base-pairs via a topological electrostatic potential

We introduce an anisotropic model of electrostatic interaction based on topological atoms obtained from the gradient vector field partitioning of the electron density. High-order electrostatic moments are computed within the compact spherical tensor formalism from ab initio wave functions of rigid geometry-optimized monomers. When combined with a simple hard-sphere or Lennard-Jones potential this topological electrostatic interaction potential correctly predicts the geometries of 27 DNA base-pairs. With the Lennard-Jones(LJ) potential internuclear distances of frontier atoms between base-pairs differ only by 0.08 A compared to supermolecular calculations at B3LYP/6-31G(d,p) level. The discrepancy for angles involved in hydrogen bond contacts amounts to only 3.5°. The topological model globally reproduces the correct ranking in intermolecular interaction energy (6.5 ± 6.5 kJ mol−1). Subsequently we compare the interaction energy profiles of the topological model with distributed multipole analysis, Merz–Kollman charges and the natural population analysis at the B3LYP/6-311+G(2d,p) level. The convergence of the topological multipole expansion is somewhat more favorable than that based on DMA but the two models have similar basis set dependence. This work encourages the development of a topological intermolecular potential, which already predicts a reliable geometry of a DNA tetrad in its present form.

Role of Cationization and Multimers Formation for Diastereomers Differentiation by Ion Mobility-Mass Spectrometry

AbstractStereochemistry plays an important role in biochemistry, particularly in therapeutic applications. Indeed, enantiomers have different biological activities, which can have important consequences. Many analytical techniques have been developed in order to allow the identification and the separation of stereoisomers. Here, we focused our work on the study of small diastereomers using the coupling of traveling wave ion mobility and mass spectrometry (TWIMS-MS) as a new alternative for stereochemistry study. In order to optimize the separation, the formation of adducts between diastereomers (M) and different alkali cations (X) was carried out. Thus, monomers [M + X]+ and multimers [2M + X]+ and [3M + X]+ ions have been studied from both experimental and theoretical viewpoints. Moreover, it has been shown that the study of the multimer [2Y + M + Li]+ ion, in which Y is an auxiliary diastereomeric ligand, allows the diastereomers separation. The combination of cationization, multimers ions formation, and IM-MS is a novel and powerful approach for the diastereomers identification. Thus, by this technique, diastereomers can be identified although they present very close conformations in gaseous phase. This work presents the first TWIMS-MS separation of diastereomers, which present very close collision cross section thanks to the formation of multimers and the use of an auxiliary diastereomeric ligand. Figureᅟ

Improved convergence of the 'atoms in molecules' multipole expansion of electrostatic interaction

We study the anisotropic electrostatic interaction with high rank multipoles (ℓ = 5) for a set of 10 van der Waals complexes and 27 DNA base pairs. Multipoles are generated by the distributed multipole analysis (DMA) and by the topological quantum theory of ‘atoms in molecules’ (AIM). The convergence of the multipolar expansion of the interaction between topological atoms is improved by distributing the moments over extra off-nuclear sites via a shifting procedure. A clear theoretical distinction is made between partitioning and distributing. A substantial improvement in the convergence of the AIM multipole expansion is observed for the smaller van der Waals molecules. An AIM representation with extra sites on the bond midpoints performs as well as DMA with just nuclear sites. However, for the larger DNA base pairs no improvement follows from the introduction of extra sites with AIM. For these larger systems AIM and DMA expansions perform equivalently, provided that now DMA allows for extra sites. This work further encourages the development of a topological intermolecular force field.

A proposal for an extended dual descriptor: a possible solution when Frontier Molecular Orbital Theory fails

In this paper, we introduce new local descriptors in the framework of Conceptual Density Functional Theory. They can be considered as an extension of the dual descriptor [Morell et al., J. Phys. Chem. A, 2005, 109, 205]. These indices are particularly suited for the discrimination between electrophilic and nucleophilic sites inside a molecule. They are computed using the densities of the electronic excited states, giving a picture of the polarization of the electron density induced by the approach of a reactant. Links with the linear-response function are discussed, and the first examples of applications are given, highlighting how these new descriptors can be used in practice for reactivity studies. It has been found that this extension of the dual descriptor can handle tricky cases, such as nitrobenzene or isoquinoline, for which Frontier Molecular Orbital Theory fails.

On the prediction of thermal stability of nitroaromatic compounds using quantum chemical calculations

This work presents a new approach to predict thermal stability of nitroaromatic compounds based on quantum chemical calculations and on quantitative structure-property relationship (QSPR) methods. The data set consists of 22 nitroaromatic compounds of known decomposition enthalpy (taken as a macroscopic property related to explosibility) obtained from differential scanning calorimetry. Geometric, electronic and energetic descriptors have been selected and computed using density functional theory (DFT) calculation to describe the 22 molecules. First approach consisted in looking at their linear correlations with the experimental decomposition enthalpy. Molecular weight, electrophilicity index, electron affinity and oxygen balance appeared as the most correlated descriptors (respectively R(2)=0.76, 0.75, 0.71 and 0.64). Then multilinear regression was computed with these descriptors. The obtained model is a six-parameter equation containing descriptors all issued from quantum chemical calculations. The prediction is satisfactory with a correlation coefficient R(2) of 0.91 and a predictivity coefficient R(cv)(2) of 0.84 using a cross validation method.

Development of a QSPR model for predicting thermal stabilities of nitroaromatic compounds taking into account their decomposition mechanisms

The molecular structures of 77 nitroaromatic compounds have been correlated to their thermal stabilities by combining the quantitative structure–property relationship (QSPR) method with density functional theory (DFT). More than 300 descriptors (constitutional, topological, geometrical and quantum chemical) have been calculated, and multilinear regressions have been performed to find accurate quantitative relationships with experimental heats of decomposition (−ΔH). In particular, this work demonstrates the importance of accounting for chemical mechanisms during the selection of an adequate experimental data set. A reliable QSPR model that presents a strong correlation with experimental data for both the training and the validation molecular sets (R2 = 0.90 and 0.84, respectively) was developed for non-ortho-substituted nitroaromatic compounds. Moreover, its applicability domain was determined, and the model’s predictivity reached 0.86 within this applicability domain. To our knowledge, this work has produced the first QSPR model, developed according to the OECD principles of regulatory acceptability, for predicting the thermal stabilities of energetic compounds.