A perfectly planar Al13+ cluster (CI) and a quasi-planar Al13+ cluster (CII) have been found for the first time. Both clusters have a triangular core surrounded by a set of ten Al atoms in the form of a ring. These cationic clusters have substantial aromatic character. The planar CI cluster has local antiaromatic patches within global aromatic sea. It is doubly aromatic having both σ and π aromatic character. The quasi-planar CII cluster is also aromatic but it has more σ-delocalization.
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.
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.
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
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).
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.
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.
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
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.
The methylxanthines have attracted interest due to the changes on their biological activities and physicochemical properties in terms of the number and position of the methyl groups present in the xanthine moiety. We report a theoretical study of the influence of the methyl substituent in the basicity and reactivity of xanthine and its methylated derivatives. Our results provide that when the xanthine increases the number of methyl substituents, the gas phase basicity increases (reactivity to proton increases), and the global hardness decreases. The result is in agreement with the maximum hardness principle (MHP) that states, “at equilibrium, chemical systems are as hard as possible” (Pearson, R.G., J. Chem. Educ., 1987, 64, 561–567, and Parr R.G., Chattaraj P.K., J. Am. Chem. Soc. 1991, 113, 1854–1855).
Efficient design of ionic compounds requires a systematic understanding of cation–anion interactions. Weakening of electrostatic attraction is essential to increase the liquid range of the ionic compound and decrease its melting point. Here, we report simulations of the closest-approach cation–anion distances in a variety of ion pairs containing the tetrakis(pentafluorophenyl)borate (TFPB–) anion. Small alkali cations (Li+, Na+) penetrate the TFPB– core, whereas K+ and larger organic cations do not. In the latter case, the shortest possible distance from the cations to the boron atom of TFPB– ranges from 0.50 nm to 0.63 nm. TFPB– was shown to be substantially rigid, providing a steric hindrance to thermodynamically efficient cation–anion coordination. Our results prove that TFPB– is more efficient for electrostatic charge confinement than the tetraoctylammonium cation, whereas the perfluorophenyl group is more efficient than linear alkyl chains. These simulations will motivate development of TFPB–-based ionic liquids with low phase transition points.
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.
This work investigated interactions between calcium cations (Ca2+) and three common types of oxygen-based functional groups of concrete superplasticizers using density functional theory (DFT) calculations and all-atom molecular dynamics (MD) simulations. The three common types of oxygen-based functional groups were modeled as three hypothetical, low-molecular-weight organic molecules, each containing a methyl-terminated oxyethylene dimer and an adsorbing head of two oxygen-based functional groups, and are referred to as carboxylate, sulfonate, and phosphate groups, respectively, following the usual terminology in the field of concrete admixtures. Our DFT results show that the binding strength of the three groups with calcium cations follows (from high to low) phosphate>carboxylate>sulfonate, and both the electrophilic attack and the chemical reactivity of the three groups contribute significantly to the binding strength. The MD simulation results indicate that the adsorption of the three small molecules on the calcite (1 0 4) surface in aqueous solution shares a similar pattern in the sense that just two oxygen atoms of two adjacent anchor groups adsorb on the calcium atoms on the top layer of the crystal. The adsorption strength among the three types of functional groups follows the same order as the binding strength obtained from DFT calculations; both results corroborate a similar rule-of-thumb established by experiments. Furthermore, interactions of the three types of groups with water molecules suggest that strong hydrogen-bonding interactions exist in those systems.
We employ computational methods to investigate the possibility of using electron-donating or electron-withdrawing substituents to reduce the free-energy barriers of the thermal isomerizations that limit the rotational frequencies achievable by synthetic overcrowded alkene-based molecular motors. Choosing as reference systems one of the fastest motors known to date and two variants thereof, we consider six new motors obtained by introducing electron-donating methoxy and dimethylamino or electron-withdrawing nitro and cyano substituents in conjugation with the central olefinic bond connecting the two (stator and rotator) motor halves. Performing density functional theory calculations, we then show that electron-donating (but not electron-withdrawing) groups at the stator are able to reduce the already small barriers of the reference motors by up to 18 kJ mol?1. This result outlines a possible strategy for improving the rotational frequencies of motors of this kind. Furthermore, exploring the origin of the catalytic effect, it is found that electron-donating groups exert a favorable steric influence on the thermal isomerizations, which is not manifested by electron-withdrawing groups. This finding suggests a new mechanism for controlling the critical steric interactions of these motors.
Graphical Abstract The introduction of electron-donating groups in one of the fastest rotary molecular motors known to date is found to reduce the free-energy barriers of the thermal steps that limit the rotational frequencies by up to 18 kJ mol?1.
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.
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.
Theoretical calculations for the first tri-iron-based extended metal atom chain (EMAC) molecule are reported. The studied triple-high-spin (S?=?6) complex exhibits ferromagnetic ordering (according to Ising and spin-projection approximations), which renders it unique among all previously prepared and theoretically calculated EMAC compounds. This ordering originates from the prevailing ferromagnetic nearest-neighbor interactions, while the magnetic superexchange between terminal Fe2+ sites is weaker and antiferromagnetic. Calculations indicate that this linear chain system based on a tri-iron core shows potential for the development of spin-frustrated behavior, which could be achieved through rational modification of the equatorial and axial ligands.
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.
