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.
In the present work, a series of simulation tools were used to determine structure-activity relationships for the endomorphins (EMs) and derive μ-pharmacophore models for these peptides. Potential lowest energy conformations were determined in vacuo by systematically varying the torsional angles of the Tyr1-Pro2 (ω1) and Pro2-Trp3/Phe3 (ω2) as tuning parameters in AM1 calculations. These initial models were then exposed to aqueous conditions via molecular dynamics simulations. In aqueous solution, the simulations suggest that endomorphin conformers strongly favor the trans/trans pair of the ω1/ω2 amide bonds. From two-dimensional probability distributions of the ring-to-ring distances with respect to the pharmacophoric angles for EMs, a selectivity range of μ1 is ca. 8.3 ~ 10.5 Å for endomorphin-2 and selectivity range of μ2 is ca. 10.5 ~ 13.0 Å for endomorphin-1 were determined. Four-component μ-pharmacophore models are proposed for EMs and are compared to the previously published δ- and κ-pharmacophore models.
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.
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.
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.
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.
There are four types of aromaticity criteria: energetic, electronic, magnetic and geometric. The delocalization, density and degeneracy-based index of aromaticity, D3BIA, is an electronic aromaticity index from QTAIM that is not reference dependent and can be used for aromatic, homoaromatic, sigma aromatic and other aromatic systems with varying ring size containing hetereoatoms or not. We used B3LYP, MP2 and MP3 methods to search for linear relations between well-known aromaticity indices and D3BIA for a series of acenes. We found that the D3BIA versus FLU correlation exceeded 91 % and reasonably good correlations exist between D3BIA and HOMA and between D3BIA and PDI. Previous works have shown that D3BIA can be used for homoaromatic systems and tetrahedrane derivatives (sigma aromaticity), but no previous work has validated D3BIA for benzenoid systems. This is the first time we have shown that D3BIA can be used successfully for benzenoid systems, for example, acenes. This work supports and validates the use of D3BIA in classical aromatic systems.
The dual role of the ionic liquid 1-butyl-3-methyl-imidazolium trifluoroacetic acid ([C4mim]TFA) as an extractant for thiophene (TH) and a catalyst for the oxidation of TH was explored at the molecular level by performing density functional theory (DFT) calculations. The calculated interaction energies demonstrated why [C4mim]TFA is a better extractant for thiophene sulfone (THO2) than for TH. Two pathways were proposed for the oxidation of TH to THO2 with [C4mim]TFA acting as a catalyst. In the dominant pathway, a peracid is formed which then oxidizes TH to the sulfoxide and sulfones. The presence of [C4mim]TFA was found to greatly reduce the barrier to the oxidative desulfurization (ODS) of TH using H2O2 as an oxidant.
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.
In the development of quantum computing and communications, improvements in materials capable of single photon emission are of great importance. Advances in single photon emission have been achieved experimentally by introducing nitrogen-vacancy (N-V) centers on diamond nanostructures. However, theoretical modeling of the anisotropic effects on the electronic properties of these materials is almost nonexistent. In this study, the electronic band structure and density of states of diamond nanowires with N-V defects were analyzed through first principles approach using the density functional theory and the supercell scheme. The nanowires were modeled on two growth directions [001] and [111]. All surface dangling bonds were passivated with hydrogen (H) atoms. The results show that the N-V introduces multiple trap states within the energy band gap of the diamond nanowire. The energy difference between these states is influenced by the growth direction of the nanowires, which could contribute to the emission of photons with different wavelengths. The presence of these trap states could reduce the recombination rate between the conduction and the valence band, thus favoring the single photon emission.
In the domain of production of flavor delivery systems based on carbohydrate matrices, Hot Melt Extrusion (HME) is one of the processes used which implies a relatively slow water evaporation step compared to Spray Drying or other common processes. The HME process of generating liquid syrup is stable only at industrial scales. The drying kinetic associated to the syrup formation is followed by collecting the vapored water. This experiment is exploited to estimate the diffusion coefficient of liquid water in a continuously concentrating carbohydrate melt ranging from 30.4 to 10.2%w/w water content and at temperatures ranging from 100 to 116 °C. A combination of Fick, Arnold and Vogel-Tammann-Fulcher models was used to derive the viscosity of the solution, the gradient thickness at the surface of the stirred solution and the diffusion coefficient of water in the liquid under boiling conditions. It is here demonstrated that under continuous boiling conditions, the water concentration gradient depth from the surface of the solution has a dimension that is comparable to the size of the hydrated carbohydrate molecules and the diffusion coefficient of water in boiling carbohydrate solutions ranges from 5.10−9 cm.s−2 for low temperature and high water content to 10−11 cm.s−2 for high temperature and low water content.
Unknown force-field parameters for metal organic beryllium complexes used in emitting and electron transporting layers of OLED structures are determined. These parameters can be used for the predictive atomistic simulations of the structure and properties of amorphous organic layers containing beryllium complexes. The parameters are found for the AMBER force field using a relaxed scan procedure and quantum-mechanical DFT calculations of potential energy curves for specific internal (angular) coordinates in a series of three Be complexes (Bebq2; Be(4-mpp)2; Bepp2). The obtained parameters are verified in calculations of some molecular and crystal structures available from either quantum-mechanical DFT calculations or experimental data.
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.
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 characterization of the seleno-sulfide-bromo systems and the isomerization process on the [H, S, Se, Br] potential energy surface were investigated using state-of-the-art theoretical methods. The CCSD(T) and the MP2 levels of theory were employed along with the series of correlation consistent basis sets extrapolated to the complete basis set (CBS) limit in the optimization of the geometrical parameters and computation of electronic energies. The relative stability, in kcal mol?1, at the CCSD(T)/CBS follows the trend: HSSeBr (0) > HSeSBr (9.51) > SSeHBr (24.02) > SeSHBr (25.42). This order was observed in the previous study of the [H, S, Se, Cl] species. The structural parameters and vibrational frequencies of the [H, S, Se, Br] species are reported. This research work should be helpful to experimentalists in order to gain insights into these novel heteroatom molecules.
Graphical abstract Relative energy profile (in kcal mol?1) using the CCSD(T)/CBS and MP2/CBS (in parentheses) method of the stationary states on the [H, S, Se, Br] PES.
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
Ionic liquids (ILs) constitute a fast growing class of compounds finding multiple applications in science and technology. Morpholinium-based ILs (MBILs) and their mixtures with polar molecular co-solvents are interesting as sustainable electrolyte systems for electrochemistry. We investigate local structures of protic and apropic morpholinium cations in acetonitrile (ACN) using semi-empirical molecular dynamics (MD) simulations. An impact of an anion (acetate) on the cation solvation regularities is discussed. Unlike oxygen, nitrogen of the morpholine ring is a strong electrophilic binding center. This site is responsible for the interactions of the cation with the solvent and with the anion. In protic MBILs, the role of nitrogen is delegated to the proton, which is linked to nitrogen. The acetate anion weakens solvation of the cation due to occupation of space near nitrogen or proton. The analysis reveals a favorable solvation of MBILs in ACN, which is a prerequisite for a new high-performance electrolyte system. The reported structural data were validated through point-to-point comparison with the MP2 post-Hartree-Fock theory and density functional theory.
A mechanistic investigation using Becke3LYP density functional theory (DFT) was carried out on the palladium-catalyzed amidition of bromobenzene and tBu-isocyanide. The whole catalytic cycle consists of five steps: oxidative addition, migratory insertion, anion exchange, reductive elimination, and hydrogen migration. The rate-determining step is oxidative addition, with a small Gibbs free energy of 14.6 kcal mol?1. In the migratory insertion step, tBu-isocyanide provides an important source of carboxy and amino groups to establish the amide group. For anion exchange, path 1a is suggested as the most favorable pathway with the help of the base, and water provides a source of oxygen which is perfectly in line with experimental observations. Finally, in the hydrogen migration step, we illustrate that the six-membered ring path is energetically favored due to the assisting influence of water. In addition, our calculations indicate that using dimethyl sulfoxide as a solvent does not change the rate-determining step.
