Catalytic fields illustrate topology of the optimal charge distribution of a molecular environment reducing the activation energy for any process involving barrier crossing, like chemical reaction, bond rotation etc. Until now, this technique has been successfully applied to predict catalytic effects resulting from intermolecular interactions with individual water molecules constituting the first hydration shell, aminoacid mutations in enzymes or Si→Al substitutions in zeolites. In this contribution, hydrogen to fluorine (H→F) substitution effects for two model reactions have been examined indicating qualitative applicability of the catalytic field concept in the case of systems involving intramolecular interactions.
The functionalization of graphene with transition metals is of great interest due to its wide range of applications, such as hydrogen storage, spintronics, information storage, etc. Due to its magnetic property adsorption of Mn atom on graphene has a high consequence on the electronic properties of graphene. The increase in size of the graphene sheet with hydrogen termination has a high impact on the transformation of electronic properties of the graphene sheet. Hence in this work, we investigate the size as well as change in structural and electronic properties of pristine/defective graphene sheets on adsorption of Mn atom using density functional theory methods. From the results obtained a higher adsorption energy value of 3.04 eV is found for Mn adatom on the defected graphene sheet than the pristine, 1.85 eV. It is subject to the coverage effect which decreases on increasing number of carbon atoms. Moreover, a decrease in energy gap is observed in pristine and defected graphene sheets with a high number of carbon atoms. The density of states illustrates the significant effect for hydrogen termination in the conduction band of the Mn adsorbed graphene sheet with low carbon atoms.
MP2/aug-cc-pVTZ calculations were performed for complexes linked by hydrogen bonds. Three types of proton donating species were taken into account: H2O, CCl3H, and H3O+. These calculations are supported by the natural bond orbital (NBO) method and the quantum theory of atoms in molecules (QTAIM) approach. Numerous correlations between parameters of H-bonded systems were found. The most important are those which show the response of the system on the H-bond formation; for example, the increase of polarization of the A-H bond correlates with the strength of the hydrogen bond. Similar relationships were found for the σ-hole bonds while the π-hole bonds do not follow the trends known for the hydrogen bonds.
Quantum chemical computations were used for prediction of the structure and color of alizarin complex with alkali metal hydroxides in methanolic solutions. The color prediction relying on the single Gaussian-like band once again proved the usefulness of the PBE0 density functional due to the observed smallest color difference between computed and experimentally derived values. It was found that the alkali metal hydroxide molecules can bind to the two oxygen atoms of both hydroxyl groups of alizarin or to one of these atoms and the oxygen atom from the keto group in a complex with three methanol molecules. This means that two electronic transitions need to be taken into account when considering the spectra of the studied complexes. The resulting bond lengths and angles are correlated with the properties of the alkali metal atoms. The molar mass, the atomic radius, and the Pauling electronegativity of studied metals are quite accurate predictors of the geometric properties of hydroxide complexes with alizarin in methanol solution.
Graphical abstract The spectra of the neutral and monoanionic form of alizarin together with color changes resulting from addition of different metal hydroxides and represented in CIE color space
Magnetic shielding constants for an isolated fullerene C60, cucurbituril CB[9], and the host-guest complex C60@CB[9] were calculated as a function of separation of the monomers. Our results in the gas phase and water indicate a significant variation of the magnetic properties for all atoms of the monomers in the complex and after liberation of fullerene C60 from the interior of the CB[9] cavity. The interaction between the two monomers results in a charge transfer that collaborates with a redistribution of electron density to deshield the monomers.
To explore the adsorption mechanism of NO, NH3, N2 on a carbon surface, and the effect of basic and acidic functional groups, density functional theory was employed to investigate the interactions between these molecules and carbon surfaces. Molecular electrostatic potential, Mulliken population analyses, reduced density gradient, and Mayer bond order analyses were used to clarify the adsorption mechanism. The results indicate that van der Waals interactions are responsible for N2 physisorption, and N2 is the least likely to adsorb on a carbon surface. Modification of carbon materials to decorate basic or acidic functional groups could enhance the NH3 physisorption because of hydrogen bonding or electrostatic interactions, however, NO physisorption on a carbon surface is poor. Zig-zag sites are more reactive than armchair sites when these gas molecules absorb on the edge sites of carbon surface.
The local and condensed Fukui functions as well as the principle of hard and soft acids and bases were used to study the addition of free radicals to the exocyclic and endocyclic double bonds of seven monocyclic monoterpenes of formula C10H16. The results obtained showed that, in general, the most reactive double bond was the one with the most substituents on the double-bonded carbon atoms, and that the reaction of a double bond with an electrophile is a soft–soft interaction. The effects of substituents on the double-bonded carbon atoms and the stabilization of the monoterpenes were interpreted by invoking hyperconjugated structures, which led us to propose a simple rule: the larger the value of the Fukui function for the double bond, the greater the hyperconjugative stabilization and the susceptibility of the double bond to electrophilic attack. In general, our results are in good accordance with relevant experimental and theoretical results published in the literature.
Bond critical points (BCPs) in the quantum theory of atoms in molecules (QTAIM) are shown to be a consequence of the molecular topology, symmetry, and the Poincaré-Hopf relationship, which defines the numbers of critical points of different types in a scalar field. BCPs can be induced by a polarizing field or by addition of a single non-bonded atom to a molecule. BCPs and their associated bond paths are therefore suggested not to be a suitable means of identifying chemical bonds, or even attractive intermolecular interactions.
We propose a new pathway for the adsorption of benzyl alcohol on the surface of TiO2 and the formation of interfacial surface complex (ISC). The reaction free energies and reaction kinetics were thoroughly investigated by density functional calculations. The TiO2 surfaces were modeled by clusters consisting of 4 Ti atoms and 18 O atoms passivated by H, OH group and H2O molecules. Compared with solid-state calculations utilizing the periodicity of the materials, such cluster modeling allows inclusion of the high-order correlation effects that seem to be essential for the adsorption of organic molecules onto solid surfaces. The effects of both acidity and solvation are included in our calculations, which demonstrate that the new pathway is competitive with a previous pathway. The electronic structure calculations based on the relaxed ISC structures reveal that the chemisorption of benzyl alcohol on the TiO2 surface greatly alters the nature of the frontier molecular orbitals. The resulted reduced energy gap in ISC matches the energy of visible light, showing how the adsorption of benzyl alcohol sensitizes the TiO2 surface.
Graphical Abstract The chemisorption of benzyl alcohol on TiO2 surface greatly alters the nature of the frontier molecular orbitals and the formed interfacial surface complex can be sensitized by visible light
Molecular simulations were carried out to study the sodium dodecyl sulfate (SDS) surfactant with the interleukin 8 (IL8) protein as a model to investigate the influence of amphiphilic molecules on proteins. Simulations for an SDS micelle with an IL8 protein show that both aggregates, which were initially separated, eventually approach each other to form a single complex. The results showed that the protein was attached to the SDS micelle by the charged positive amino acids whereas less contacts were observed for the negatively charged amino acids. Structural protein properties, such as amino acid contacts and pair correlation functions were conducted between the micelle and the protein groups and they showed greater interactions between the surfactant headgroups and the positively charged residues in the protein. Moreover, hydrogen bonds were also calculated between both structures and a greater number of bonds among the SDS headgroups and the charged positive amino acids in the protein was found.
Graphical Abstract Attachment of the interleukin 8 protein with a sodium dodecyl sulfate (SDS) micelle is given by the charged positive amino acids as indicated by the interaction of those amino acids with the SDS headgroups.
Coarse-grained force field (CGFF) methods were applied to study the self-assembly of sodium dodecyl sulfate with fragrance additives. The CGFF parameters were parameterized and validated using experimental and all-atom simulation data. Direct molecular dynamics simulations were carried out to characterize the initial aggregation, partitioning of fragrances, and chemical potentials of the surfactant and fragrance molecules in aggregates of different sizes. The equilibrium critical micelle concentrations (CMCs) and micelle size distributions, which could not be obtained by direct simulation, were predicted using the calculated chemical potentials in combination with a thermodynamic model. The predicted partitioning of fragrances, CMCs, micelle sizes, and micelle structures agree well with previously reported experimental data.
Phthalocyanines represent a crucial class of organic compounds with high technological appeal. By doping the center of these systems with metals, one obtains the so-called metal-phthalocyanines, whose property of being an effective electron donor allows for potentially interesting uses in organic electronics. In this sense, investigating optical and electronic structure changes in the phthalocyanine profiles in the presence of different metals is of fundamental importance for evaluating the appropriateness of the resulting system as far as these uses are concerned. In the present work, we carry out this kind of effort for phthalocyanines doped with different metals, namely, copper, nickel, and magnesium. Density functional theory was applied to obtain the absorption spectra, and electronic and structural properties of the complexes. Our results suggest that depending on the dopant, a different level of change is achieved. Moreover, electrostatic potential energy mapping shows how the charge distribution can be affected by solar radiation. Our contribution is crucial in describing the best possible candidates for use in different organic photovoltaic applications.
Graphical Abstract Representation of meta-phthalocyanine systems. All calculations of this work are based on varying metal position along z axis, considering the z-axis has its zero point matching with the center of phthalocyanine cavityconsidering.
The catalytic coupling reaction mechanism of the transformation from 4-nitrobenzenethiol (4-NBT) to 4,4′-dimercaptoazobenzene (4,4′-DMAB) on a silver cluster was studied by density functional theory. Reactants, intermediates, transition states and products were optimized with the B3LYP method using the 6-311 + G(d,p) basis set (Ag using the pseudo potential basis set of LanL2DZ). Transition states and intermediates were confirmed by the corresponding vibration analysis and intrinsic reaction coordinates (IRC). Consistent with literature reports, the key point of the transformation from 4-NBT absorbed on the surface of Ag5 clusters to 4,4′-DMAB is the elimination of two O atoms on the amino group. Meanwhile, the catalytic coupling reaction of 4-nitrobenzenethiol on a silver cluster is easy to carry out under irradiation. The possibility of “inter system channeling” (ISC) between different potential energy surfaces in the coupling reaction of 4-NBT is further discussed. The irradiation has an auxiliary catalytic effect on the coupling reaction. Our research results can explain the observed experimental phenomena.
Graphical abstract Catalytic coupling reaction mechanism of the transformation from 4-nitrothiophenol (4-NBT) to 4,4′-dimercaptoazobenzene (4,4′-DMAB) on silver clusters studied by density functional theory
In this article, we explore, both theoretically and experimentally, the general reactivity of alkyl hydrogeno-phenylphosphinates with alcohols. We show that alcohol molecules act exclusively as nucleophilic species, and add to alkyl hydrogeno-phenylphosphinates, leading to pentacoordinated intermediates. These intermediates are shown to subsequently competitively undergo alcohol eliminations and/or Berry pseudorotations. This offers several possible routes for racemizations and/or alcohol exchange reactions. Transition standard Gibbs free energies predicted from DFT calculations for the overall alcohol exchange mechanism are shown to be compatible with those experimentally measured in case ethanol reacts with ethyl hydrogeno-phenylphosphinate (134.5～136.0 kJ mol?1 at 78 °C).
The family of human epidermal growth factor receptors (HER) is involved in tumor cell growth. Homodimerization and heterodimerization of the HER family are important for activation of these receptors. The structures of homodimer conformation are well characterized, while the structures of heterodimer conformations, especially between HER1 and HER2, are not completely understood. In this study, two models of possible asymmetric HER1/HER2 kinase domains were built. Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) methods were applied to examine the possibility of these two-heterodimer interactions. From our results, it could be concluded that the HER2 kinase domain prefers to serve as the receiver rather than the activator. Key binding residues of this dimer complex at N lobe of HER2 is ALA683 and at C lobe of HER1 are GLU914, GLU917, and ASP930. This study will be useful in allowing us to predict and be able to control activity of this enzyme in disease in the future.
Graphical abstract A model of the asymmetric dimer of HER1–HER2 heterodimer with key intereaction residues in (a) HER1A and (b) HER2R by molecular dynamic simulation
Vitamin C is one of the most abundant exogenous antioxidants in the cell, and it is of the utmost importance to elucidate its mechanism of action against radicals. In this study, the reactivity of vitamin C toward OH and \( {HO}_2/{O}_2^{-} \) radicals in aqueous medium was analyzed by ab initio molecular dynamics using CPMD code. The simulations led to results similar to those of static studies or experiments for the pair of \( {HO}_2/{O}_2^{-} \) radicals but bring new insights for the reactivity with hydroxyl radical: the reaction takes place before the formation of an adduct and consists of two steps: first an electron is transferred to hydroxyl radical and then the ascorbyl radical loses a proton.
One of the central assumptions when a particle moves through a window in microporous materials is that interaction of the diffusing particle with the silicon (Si) and aluminum (Al) atoms of the framework can be neglected, as the presence of bulkier oxygen in the host structure is thought to hinder close proximity of the diffusing particle to Si and Al. We examine this assumption, exploring the diffusion path and cross-checking the bottleneck associated with the diffusion process. Our study reveals that short-range interactions between the diffusing species and Si/Al of the host have a significant effect on the diffusion process. Guest–host interaction energy increases significantly if interaction between Si and Al atoms with the diffusing species is considered. The self-diffusion coefficient (D) decreases significantly in the linear regime, whereas in the anomalous regime, surprisingly, D increases. The increase in D is due to a decrease in the activation energy in the anomalous regime, whereas in the linear regime, activation energy increases, thus D decreases.
Graphical abstracta Interaction energies (Ea) for different LJ potential for guest–guest interactions (σgg) along the diffusion path; b correspondingdiffusivity values
A new compound based on the D-π-A concept, where D = dimethylamino-phenyl and A = naphthoic acid, separated by an imine motif, was designed, synthesized and characterized. The spectral, energetics, and structural characteristics of the compound were studied thoroughly theoretically by density functional theory (DFT) in the gas and aqueous phases and experimentally (steady-state absorption) in aqueous media with various degrees of polarity and hydrogen bonding ability. This compound shows high sensitivity to the polarity, basicity and proton affinity of the environment. Based on DFT, TD-DFT and NBO analysis, the compound exists in the ground-state with both intermolecular and intramolecular hydrogen bond conformations in association with the –COOH, with latter isomer calculated to be more stable. Furthermore, structural changes via intermolecular solute–solvent interactions, dictate electronic modifications and spectral changes.
