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).
Recently, a series of xanthone analogues has been identified as α-glucosidase inhibitors. To provide deeper insight into the three-dimensional (3D) structural requirements for the activities of these molecules, CoMFA and CoMSIA approaches were employed on 54 xanthones to construct 3D-QSAR models. Their bioactive conformations were first investigated by docking studies and optimized by subsequent molecular dynamics (MD) simulations using the homology modeled structure of the target protein. Based on the docking/MD-determined conformers, 3D-QSAR studies generated several significant models in terms of 47 molecules as the training set. The best model (CoMSIA-SHA) yielded q2 of 0.713, r2 of 0.967 and F of 140.250. The robustness of the model was further externally confirmed by a test set of the remaining molecules (q2 = 0.793, r2 = 0.902, and k = 0.905). Contour maps provided much information for future design and optimization of new compounds with high inhibitory activities towards α-glucosidase.
Density functional theory (B3LYP, B3LYP-D2 and wB97XD functionals) was used in finite models of zigzag carbon nanotubes (CNT), (n,0)×k with n?=?6–9 and k?=?2–4, to systematically investigate the effects of size on their structural and electronic properties. We found that the ratio between the length (Lt) and the diameter (dt) of the pristine CNT has to be larger than 2, i.e., Lt/dt?>?2, in order to provide the observed experimental trends of C=C bond distances, as well as to maintain the atomic charges nearly constant and zero around the center of the tube. Therefore, the concepts of useful length and volume were developed and tested for the encapsulation process of HCN and C2H2 into CNTs. The energies involved in these processes, as well as the changes in molecular structure and electronic properties of the dopants and the CNTs are discussed and rationalized by the amount of charge transferred between dopant and CNT.
A grid-based variational method was proposed and applied to the ground state energies of atoms from the first to the third period of the periodic table. The nonuniform grid in the radial coordinate was defined by a q-exponential sequence. Some unusual properties between the optimum q-parameters and the electronic energies or atomic numbers are described. The behavior of the electronic energy, with respect to the q-parameter, yields near Hartree-Fock accuracy with a relatively small number of integration points. A simple relationship between the optimum q-parameters and the atomic numbers was found, which allowed the determination of the optimum q-parameters for atoms of the same period from two results. The remarkable results provide a simple alternative route to reach accurate results. The consistent results also suggest that this is not a random or accidental effect, but some optimum condition achieved by using a q-exponential mesh grid.
The diffusion coefficients of 14 n-alkanes (ranging from methane to n-tetradecane) in liquid and supercritical methanol at infinite dilution (at a pressure of 10.5 MPa and at temperatures of 299 K and 515 K) were deduced via molecular dynamics simulations. Values for the radial distribution function, coordination number, and number of hydrogen bonds were then calculated to explore the local structure of each fluid. The flexibility of the n-alkane (as characterized by the computed dihedral distribution, end-to-end distance, and radius of gyration) was found to be a major influence and hydrogen bonding to be a minor influence on the local structure. Hydrogen bonding reduces the flexibility of the n-alkane, whereas increasing the temperature enhances its flexibility, with temperature having a greater effect than hydrogen bonding on flexibility.
Graphical abstract The flexibility of the alkane is a major influence and the hydrogen bonding is a minor influence on the first solvation shell; the coordination numbers of long-chain n-alkanes in the first solvation shell are rather low
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
Density functional theory (DFT) was used to study the cobalt(I)-catalyzed enantioselective intramolecular hydroacylation of ketones and alkenes. All intermediates and transition states were fully optimized at the M06/6-31G(d,p) level (LANL2DZ(f) for Co). The results demonstrated that the ketone and alkene present different reactivities in the enantioselective hydroacylation. In ketone hydroacylation catalyzed by the cobalt(I)–(R,R)-Ph-BPE complex, reaction channel “a” to (R)-phthalide was more favorable than channel “b” to (S)-phthalide. Hydrogen migration was both the rate-determining and chirality-limiting step, and this step was endothermic. In alkene hydroacylation catalyzed by the cobalt(I)–(R,R)-BDPP complex, reaction channel “c” leading to the formation of (S)-indanone was the most favorable, both thermodynamically and kinetically. Reductive elimination was the rate-determining step, but the chirality-limiting step was hydrogen migration, which occurred easily. The results also indicated that the alkene hydroacylation leading to (S)-indanone formation was more energetically favorable than the ketone hydroacylation that gave (R)-phthalide, both thermodynamically and kinetically.
Graphical abstract A DFT study demonstrated that the ketone and alkene in the cobalt(I)-catalyzed enantioselective intramolecular hydroacylation showed different reactivities
Seven models that related the features of molecular surface electrostatic potentials (ESPs) above the bond midpoints and rings, statistical parameters of ESPs to the experimental impact sensitivities h50 of eight strained cyclic explosives with the C–NO2 bonds were theoretically predicted at the DFT-B3LYP/6-311++G** level. One of the models was used to investigate the changes of h50 for the nitrocyclohydrocarbon frameworks in the H-bonded complexes of HF with nitrocyclopropane, nitrocyclobutane, nitrocyclopentane, and nitrocyclohexane. The results show that the correlation coefficients of the obtained models are small. When adding the effect of ring strain, the value of correlation coefficient is increased. According to the calculated h50, the sensitivities in the frameworks are increased after hydrogen bonding. As a global feature of molecules, surface electrostatic potential is more available to judge the sensitivity change than the trigger bond dissociation energy or ring strain energy in H-bonded complex.
Graphical Abstract A theoretical prediction of the relationships between the impact sensitivity and electrostatic potential in strained cyclic explosive and application to H-bonded complex of nitrocyclohydrocarbon?
Based on dissipative particle dynamics (DPD) methods and experimental data, we used an empirical relationship between the DPD temperature and the real temperature to build a model that describes the viscosity of molten TNT fluids. The errors in the predicted viscosity based on this model were no more than 2.3 %. We also studied the steady-state shear rheological behavior of molten TNT fluids containing nanoparticles (“nanofluids”). The dependence of the nanofluid viscosity on the temperature was found to satisfy an Arrhenius-type equation, η?=?AeB/T, where B, the flow activation energy, depends on particle content, size, and shape. We modified the Einstein-type viscosity model to account for the effects of nanoparticle solvation in TNT nanofluids. The resulting model was able to correctly predict the viscosities of suspensions containing nano- to microsized particles, and did not require any changes to the physical background of Einstein’s viscosity theory.
At high temperature, silicon oxycarbide (SiCO) exhibits excellent mechanical properties and thermal stability. The incorporation of boron in SiCO results in improved performance in creep temperatures. In this work, large-scale molecular dynamics calculations were applied to obtain amorphous SiCO structures containing boron. Phase separation of C–C, B–C and Si–O was achieved for three compositions, and silicon-centered mixed-bond tetrahedrons were reproduced successfully. As the boron content increases, the boron atoms tend to form B–C and B–Si bonds in the voids, which stretches the free carbon network in some instances, causing a increase in C–C distance. Young’s modulus remains stable at high temperature for the high-carbon case, which indicates that the free carbon network plays a critical role in the structural and thermal stability of SiBCO.
Graphical Abstract Three major typical structures in the cooling down process for silicon boron oxycarbide (Si5BC2O8). Bonds: red Si–O, blue Si–C, black C–C, green B–C, purple Si–B
The catalytic pyrolysis pathways of carbonyl compounds in coal were systematically studied using density functional theory (DFT), with benzaldehyde (C6H5CHO) employed as a coal-based model compound and ZnO, γ-Al2O3, and CaO as catalysts. The results show that the products of both pyrolysis and catalytic pyrolysis are C6H6 and CO. However, the presence of any of the catalysts changes the reaction pathway and reduces the energy barrier, indicating that these catalysts promote C6H5CHO decomposition.
Graphical abstract The presence of catalysts changes the reaction pathway and the energy barrier decreases in the order Ea (no catalyst)> Ea (CaO)> Ea (γ-Al2O3)> Ea (ZnO), indicating that these catalysts promote C6H5CHO decomposition.
The present study reports the geometries, electronic structures, growth behavior, and stabilities of neutral and ionized copper-doped germanium clusters containing 1–20 Ge atoms within the framework of linear combination of atomic orbitals density functional theory (DFT) under the spin-polarized generalized gradient approximation. It was found that Cu-capped Gen (or Cu-substituted Gen+1) and Cu-encapsulated Gen clusters mostly occur in the ground state at a particular cluster size (n). In order to explain the relative stabilities of the ground-state clusters, parameters such as the average binding energy per atom (BE), the embedding energy (EE), and the fragmentation energy (FE) of the clusters were calculated, and the resulting values are discussed. To explain the chemical stabilities of the clusters, parameters such as the energy gap between the highest occupied and the lowest unoccupied molecular orbitals (the HOMO–LUMO gap), the ionization energy (IP), the electron affinity (EA), the chemical potential (μ), the chemical hardness (η), and the polarizability were calculated, and the resulting values are also discussed. Natural atomic orbital (NAO) and natural bond orbital (NBO) analyses were also used to determine the electron-counting rule that should be applied to the most stable Ge10Cu cluster. Finally, the relevance of the calculated results to the design of Ge-based superatoms is discussed.
Figure Contributions of the valance orbitals of the Ge and Cu atom(s) to the HOMO of the ground-state icosahedral Ge10Cu cluster obtained from NBO analysis. The numbers below the clusters represent the occupancies of the HOMO orbitals
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.
Molecular dynamics simulations were performed to investigate the separation of trihalomethanes (THMs) from water using boron nitride nanosheets (BNNSs). The studied systems included THM molecules and a functionalized BNNS membrane immersed in an aqueous solution. An external pressure was applied to the z axis of the systems. Two functionalized BNNSs with large fluorinated-hydrogenated pore (F-H-pores) and small hydrogen-hydroxyl pore (H-OH-pores) were used. The pores of the BNNS membrane were obtained by passivating each nitrogen and boron atoms at the pore edges with fluorine and hydrogen atoms in the large pore or with hydroxyl and hydrogen atoms in the small pore. The results show that the BNNS with a small functionalized pore was impermeable to THM molecules, in contrast to the BNNS with a large functionalized pore. Using these membranes, water contaminants can be removed at lower cost.
Graphical Abstract A snapshot of the simulation system. The BNNS membrane with the large functionalized pore is located in the middle of the box. The size of the box is 3 × 3 × 5 nm3. Green chlorine, cyan carbon, red oxygen, white hydrogen
Multilayer-shaped compression and slide models were employed to investigate the complex sensitive mechanisms of cocrystal explosives in response to external mechanical stimuli. Here, density functional theory (DFT) calculations implementing the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) with the Tkatchenko-Scheffler (TS) dispersion correction were applied to a series of cocrystal explosives: diacetone diperoxide (DADP)/1,3,5-trichloro-2,4,6-trinitrobenzene (TCTNB), DADP/1,3,5-tribromo-2,4,6-trinitrobenzene (TBTNB) and DADP/1,3,5-triiodo-2,4,6-trinitrobenzene (TITNB). The results show that the GGA-PBE-TS method is suitable for calculating these cocrystal systems. Compression and slide models illustrate well the sensitive mechanism of layer-shaped cocrystals of DADP/TCTNB and DADP/TITNB, in accordance with the results from electrostatic potentials and free space per molecule in cocrystal lattice analyses. DADP/TCTNB and DADP/TBTNB prefer sliding along a diagonal direction on the a?c face and generating strong intermolecular repulsions, compared to DADP/TITNB, which slides parallel to the b?c face. The impact sensitivity of DADP/TBTNB is predicted to be the same as that of DADP/TCTNB, and the impact sensitivity of DADP/TBTNB may be slightly more insensitive than that of DADP and much more sensitive than that of TBTNB.
Density functional theory and its time-dependent extension (DFT, TDDFT) were employed to establish the feasibility of using a series of 4,4-difluoro-4-bora-3a,4a-diaza-s-indacenes (BODIPYs) in photodynamic therapy. Their absorption electronic spectra, singlet–triplet energy gaps, and spin–orbit matrix elements were computed and are discussed here. The effects of bromine substitution on the photophysical properties of BODIPY were elucidated. The investigated compounds were found to possess different excited triplet states that lie below the energy of the bright excited singlet state (S1 or S2), depending on the positions occupied by the bromine atoms. The computed spin–orbit matrix elements for the radiationless intersystem crossing Sn?→ ?Tm and the relative singlet–triplet energy gaps allowed the prediction of plausible nonradiative decay pathways for the production of singlet excited molecular oxygen, the key cytotoxic agent in photodynamic therapy.
Graphical Abstract The photophysical properties affected by the presence of bromine atoms in different positions of a BODIPY core have been here elucidated. In particular it has been found that SOC values strongly depend on the position of heavy atoms into the BODIPY core, suggesting positions 1 and 7 as the best ones to enhance the ISC kinetics
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.
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.
