A post-calculation correction is established for PM7 band gaps of transition-metal oxides. The correction is based on the charge on the metal cation of interest, as obtained from MOPAC PM7 calculations. Application of the correction reduces the average error in the PM7 band gap from ~3 eV to ~1 eV. The residual error after correction is shown to be uncorrelated to the Hartree–Fock method upon which PM7 is based.
Graphical Abstract Comparison between calculated band gaps and experimental band gaps for binary oxides. The orange crosses are for corrected PM7 band gaps. Blue squares are uncorrected values. The orange crosses fall closer to the diagonal dashed line, showing an overall improvement of the accuracy of calculated values
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.
In this study, the doped defects in nitromethane crystals were investigated using first-principles calculations for the first time. We introduce dopant atoms in the interstitial sites of the nitromethane lattice, aiming to study the effects of element-doping on the structural properties, electronic properties, and sensitivity characteristics. The obtained results show that doped defects obviously affect the neighboring nitromethane molecules. The modification of electronic properties shows that the band gaps are significantly influenced by doped defects. Partial density of states and population analysis further reveal the mechanism for sensitivity control of nitromethane. It is shown that the new electronic states were introduced in the forbidden bands and the doped defects resulted in charge redistributions in the systems.
Aminopeptidase N (APN) is a zinc-dependent ectopeptidase involved in cell proliferation, secretion, invasion, and angiogenesis, and is widely recognized as an important cancer target. However, the mechanisms whereby ligands leave the active site of APN remain unknown. Investigating ligand dissociation processes is quite difficult, both in classical simulation methods and in experimental approaches. In this study, random acceleration molecular dynamics (RAMD) simulation was used to investigate the potential dissociation pathways of ligand from APN. The results revealed three pathways (channels A, B and C) for ligand release. Channel A, which matches the hypothetical channel region, was the most preferred region for bestatin to dissociate from the enzyme, and is probably the major channel for the inner bound ligand. In addition, two alternative channels (channels B and C) were shown to be possible pathways for ligand egression. Meanwhile, we identified key residues controlling the dynamic features of APN channels. Identification of the dissociation routes will provide further mechanistic insights into APN, which will benefit the development of more promising APN inhibitors.
The effect of alkali metal oxides MnO (M?=?Li, Na, K; n?=?2, 3, 4) on the geometric, electronic, and linear and nonlinear optical properties of the Mg12O12 nanocage was investigated by density-functional-based methods. According to the computational results, these alkali metal oxides are adsorbed on the Mg12O12 nanocage because this adsorption reduces its energy gap. The static first hyperpolarizability (β0) of the nanocage is dramatically increased in the presence of the alkali metal oxides, with the greatest increase seen in the presence of the superalkalis (i.e., M3O; M?=?Li, Na, and K). The highest first hyperpolarizability (β0?≈?600,000 a.u.) was calculated for K3O@Mg12O12, which was considerably more than that for Mg12O12. The thermodynamic properties and relative stabilities of these inorganic compounds are discussed.
High-level ab initio calculations on the complexes between noble gas atoms (He, Ne, Ar, Kr, and Xe) and dihalogen molecules (F2, Cl2, Br2, and I2) reveal trends, both in interaction energies and the energy difference between the linear and T-shaped structures, that can be explained well in terms of dispersion interactions enhanced by polar flattening of the halogens. The partial discrepancies with experimental findings are discussed.
Herein we report a study of the switchable [3]rotaxane reported by Huang et al. (Appl Phys Lett 85(22):5391–5393, 1) that can be mounted to a surface to form a nanomechanical, linear, molecular motor. We demonstrate the application of semiempirical electronic structure theory to predict the average and instantaneous force generated by redox-induced ring shuttling. Detailed analysis of the geometric and electronic structure of the system reveals technical considerations essential to success of the approach. The force is found to be in the 100–200 pN range, consistent with published experimental estimates.
A series of three star-shaped compounds containing both donor (carbazole) and acceptor (2,4,6-triphenyl-1,3,5-triazine) moieties linked through various linking bridges was studied theoretically at the linear response TD-DFT level of theory to describe their absorption and fluorescence spectra. The concept of a localized charge-transfer excited state has been applied successfully to explain the observed strong solvatochromic effect in the emission spectra of the studied molecules, which can be utilized for the fabrication of color tunable solution-processable OLEDs. The concept is in particularly applicable to donor–acceptor species with a C3 symmetry point group where the static dipole moment changes dramatically upon electronic excitation. An important peculiarity of the studied molecules is that they are characterized by non-zero values of the HOMO and LUMO orbitals in the same common part of molecular space that provides a large electric dipole transition moment for both light absorption and emission.
In this article, we explore the capacity of formed Schiff base complexes to trap metal atoms or ions, using their aromatic ends. The intrinsic geometry of each complex defines the process of substitution. Two cases were studied; one involving a trans Schiff base complex and the other considering how a salen ligand, with nickel systems traps chromium. We also assessed the nature of the new bonds and the frontier molecular orbitals.
The chelating properties of diethanoldithiocarbamate (DEDC) and π-electron flow from the nitrogen atom to the sulfur atom via a plane-delocalized π-orbital system (quasi ring) was studied using a density functional theory method. The molecular structure of DEDC and its complexes with Zn(II), Cd(II), and Hg(II) were also considered. First, the geometries of this ligand and DEDC-Zn(II), DEDC-Cd(II), and DEDC-Hg(II) were optimized, and the formation energies of these complexes were then calculated based on the electronic energy, or sum of electronic energies, with the zero point energy of each species. Formation energies indicated the DEDC-Zn(II) complex as the most stable complex, and DEDC-Cd(II) as the least stable. Structural data showed that the N1–C2 π-bond was localized in the complexes rather than the ligand, and a delocalized π-bond over S7–C2–S8 was also present. The stability of DEDC-Zn(II), DEDC-Cd(II), and DEDC-Hg(II) complexes increased in the presence of the non-specific effects of the solvent (PCM model), and their relative stability did not change. There was π-electron flow or resonance along N1–C2–S7 and along S7–C2–S8 in the ligand. The π-electron flow or resonance along N1–C2–S7 was abolished when the metal interacted with sulfur atoms. Energy belonging to van der Waals interactions and non-covalent delocalization effects between the metal and sulfur atoms of the ligand was calculated for each complex. The results of nucleus-independent chemical shift (NICS) indicated a decreasing trend as Zn(II)?<?Cd(II)?<?Hg(II) for the aromaticity of the quasi-rings. Finally, by ignoring van der Waals interactions and non-covalent delocalization effects between the metal and sulfur atoms of the ligand, the relative stability of the complexes was changed as follows:
The factors that explain the competition between intramolecular NO linkage photoisomerization and NO photorelease in five ruthenium nitrosyl complexes were investigated. By applying DFT-based methods, it was possible to characterize the ground states and lowest triplet potential energy surfaces of these species, and to establish that both photoisomerization and photorelease processes can occur in the lowest triplet state of each species. This work highlights the crucial role of the sideways-bonded isomer, a metastable state also known as the MS2 isomer, in the photochemical loss of NO, while the results obtained also indicate that the population of the triplet state of this isomer is compulsory for both processes and show how photoisomerization and photorelease interfere.
Spin-polarized first-principles total-energy calculations have been performed to investigate the possible chain reaction of acetylene molecules mediated by hydrogen abstraction on hydrogenated hexagonal boron nitride monolayers. Calculations have been done within the periodic density functional theory (DFT), employing the PBE exchange correlation potential, with van der Waals corrections (vdW-DF). Reactions at two different sites have been considered: hydrogen vacancies on top of boron and on top of nitrogen atoms. As previously calculated, at the intermediate state of the reaction, when the acetylene molecule is attached to the surface, the adsorption energy is of the order of ?0.82 eV and ?0.20 eV (measured with respect to the energy of the non interacting molecule-substrate system) for adsorption on top of boron and nitrogen atoms, respectively. After the hydrogen abstraction takes place, the system gains additional energy, resulting in adsorption energies of ?1.52 eV and ?1.30 eV, respectively. These results suggest that the chain reaction is energetically favorable. The calculated minimum energy path (MEP) for hydrogen abstraction shows very small energy barriers of the order of 5 meV and 22 meV for the reaction on top of boron and nitrogen atoms, respectively. Finally, the density of states (DOS) evolution study helps to understand the chain reaction mechanism.
The conversion of 2-phenylbenzimidazole using o-phenylenediamine and benzaldehyde can be improved significantly under β-cyclodextrin (β-CD). The density functional theory (DFT) method was applied to study the whole process. According to energy parameters (binding energy, deformation energy) and structural deformation, entry models and the reaction process can be pinpointed, with o-phenylenediamine embedding β-CD from a wide rim, and then benzaldehyde passing into the inclusion from the narrow rim. Subsequently, natural bonding orbital (NBO), Mulliken charge, frontier orbital, FuKui function and nuclear magnetic resonance (NMR) methods were employed to reveal the mechanism of electron transfer. The results illustrate that β-CD plays a catalytic role in synthesis reaction mechanism on the secondary side, improving the reactivity and selectivity of the process.
Graphical Abstract Density functional theory study of the effects of β-cyclodextrin in synthesis of 2-phenylbenzimidazole via benzaldehyde and o-phenylenediamine
The solubility advantage (SA) of meloxicam cocrystalized with mono- and dicarboxylic acids was expressed in terms of equilibrium constants involving active pharmaceutical ingredient and coformer in aqueous solutions. It is argued that SA can be quantified by concentration of pairs formed in water. The pH and concentration of dissolved components is included explicitly in the model. The alternative behavior of mono- and dicarboxylic acids was emphasized and addressed to different structural motifs. The structural and energetic properties of meloxicam and its complexes with carboxylic acids were characterized, including tautmerism and dissociation in aqueous media. In particular, performed in silico modeling confirmed experimental observation that meloxicam dissolved in water or modest acidic solutions is expected to be a mixture of anionic form in equilibrium with at least five neutral isomers. Tautomer-related diversity of pairs formation and the possibility of salt formation is also discussed.
The low bending rigidity of graphene facilitates the formation of folds into the structure. This curvature change affects the reactivity and electron transport of the sheet. One novel extension of this is the intercalation of small molecules into these folds. We construct a model incorporating a single-walled carbon nanotube into a sheet of folded graphene. Variational calculus techniques are employed to determine the minimum energy structure and the resulting curves are shown to agree well with molecular dynamics study.
Graphical Abstract Using calculus of variations, the elastic bending energy and van der Waals energy are minimised giving rise to Euler-Lagrange equation for which analytical solutions are derived to determine the optimal curved sturctures of graphene wrapped around carbon nanotubes . Overall agreement between the analytical solutions (with different values of bending rigidities) and results from molecular dynamics simulations (grey) is shown here for (6,6), (8,8) and (10,10) armchair nanotubes, respectively.
Density functional theory (DFT) was utilized to elucidate the reaction mechanisms of and the key factors that influence the Ni(0)-catalyzed cross-dimerization and -trimerization of trimethylsilylacetylene (R1) and diphenylacetylene (R2). Calculated results revealed that the electron-donating ability of the ligand plays a crucial role in determining the regionselectivity of this tandem reaction. The use of strongly electron-donating ligands favors the formation of cross-dimer intermediates, whereas cross-trimer products can easily be synthesized using weakly electron-donating ligands. A simple method of estimating the electron-donating abilities of different ligands based on the Mulliken charge distribution of the ligand–ligand pair was employed. The present theoretical results allow us to elucidate the reaction mechanisms for and to identify the factors that exert the greatest influence on the ligand-controlled cross-dimerization and -trimerization of trimethylsilylacetylene and diphenylacetylene. Guidelines for the design of novel ligand systems with Ni(0) catalysts are also proposed.
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.
