Using density functional theory (DFT) and molecular dynamics (MD), we studied the interaction of a titanium atom with a half of a C60 fullerene (i.e., C30), formed from the corannulene structure with a pentagonal base. We considered atmospheric pressure and 300 K. We found that the most stable adsorption of the titanium atom on C30 occurs in the concave surface of the molecule. Afterward, we investigated the interaction of the system C30-titanium with carbon monoxide and carbon dioxide molecules, respectively. We found that each of these molecules is chemisorbed, with no dissociation. The value of the adsorption energy for the carbon monoxide molecule varies from ?0.897 to ?1.673 eV, and for the carbon dioxide molecule, it is between ?1.065 and ?1.274 eV. These values depend on the initial orientation of these molecules with respect to TiC30.
For the first time, the structures, stabilities and electronic properties of alkaline-earth metal doped B44 fullerenes were investigated by means of density functional theory calculations. Our results reveal that M@B44 (M = Ca, Sr, Ba) possess endohedral configurations as their lowest energy structures, whereas the exohedral form is favored when metal is Be or Mg. The large binding energies and sizable HOMO–LUMO gap energies of Ca@B44, Sr@B44 and Ba@B44 suggest the considerable possibility to achieve these novel endohedral borofullerenes experimentally. Born-Oppenheimer molecular dynamics (BO-MD) simulations at various temperatures further confirmed the extreme dynamic stabilities of these endohedral complexes. Their bonding patterns were also analyzed in detail. Finally, we simulated their infrared absorption spectra and 11B nuclear magnetic resonance spectra to help future structural characterization.
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.
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.
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
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.
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 mechanistic details of N-heterocyclic olefin-catalyzed formation of cyclic carbonate from CO2 and propargylic alcohols were investigated by DFT calculations. Six mechanisms, four for the formation of five-membered cyclic carbonate (M-A, M-B, M-B’ and M-C), and two for six-membered cyclic carbonate (M-D and M-E), were fully investigated. The energy profiles in dichloromethane showed that M-B is the predominant reaction with the lowest barrier of 31.99 kcal mol?1, while M-C and M-D may be kinetically competitive to M-B. The very high activation energy of 45.37 kcal mol-1, 57.07 kcal mol-1 and 59.61 kcal mol?1 for M-A, M-B’ and M-E, respectively, suggest that they are of lesser importance in the overall mechanism.
Graphical abstract Formations of five-membered ring product and six-membered ring product are kinetically competitive, but five-membered ring product is thermodynamically more preferable.
The adsorption processes of elemental lead on carbonaceous surfaces which adsorbed CO/CO2/NO flue gases were investigated to understand the effects of CO/CO2/NO on elemental lead adsorption on carbonaceous surfaces with density functional theory. All calculations including optimizations, energies, and frequencies were conducted at B3PW91 density functional theory level, utilizing SDD basis set for lead and 6-31G(d) Pople basis set for other atoms. The results indicate that CO, CO2, and NO can promote the adsorption of elemental lead on the carbonaceous surface, but probably compete for adsorption sites with elemental lead. The promotion effects on adsorption can be attributed to active sites on the carbonaceous surface rather than flue gas adsorption on the carbonaceous surface. In addition, the adsorption order of three kinds of flue gas on the carbonaceous surface is CO2?>?NO?>?CO?>?Pb on average. Furthermore, the enhancement order of three kinds of flue gas on the elemental lead adsorption on carbonaceous surfaces is CO-CS?>?CO2-CS?>?NO-CS?>?CS in general. In particular, atomic charge and adsorption energy have good linear relationship in the process of elemental lead adsorption.
A density functional theory (DFT) study of cct-As, ccc, and cct-CO isomers of the ruthenium dihydride complex RuH2(CO)2(AsMe2Ph)2 is reported (see Scheme for the labeling isomer 34 structures of RuH2(CO)2(AsMe2Ph)2). Complex geometries and relative energies of different isomers have been calculated with both B3LYP and M06-2X functionals. The results show that the B3LYP calculated Boltzmann populations of cct-As, ccc, and cct-CO isomers are 65.5, 34.2, and 0.3%, respectively. These are in better agreement with the experimental data than those calculated at the M06-2X level. However, the calculations of 1H NMR chemical shifts were found to be better described with M06-2X than with B3LYP or with HF level of theories. In addition, a transition state between the two most stable isomers was determined through DFT/(B3LYP or M06-2X) calculations.
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.
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.
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.
Performance of 18 DFT functionals (B1B95, B3LYP, B3PW91, B97D, BHandHLYP, BMK, CAM-B3LYP, HSEh1PBE, M06-L, mPW1PW91, O3LYP, OLYP, OPBE, PBE1PBE, tHCTHhyb, TPSSh, wB97xD, VSXC) in combinations with six basis sets (cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, aug-cc-pVTZ, IGLO-II, and IGLO-III) and three methods for calculating magnetic shieldings (GIAO, CSGT, IGAIM) was tested for predicting 1H and 13C chemical shifts for 25 organic compounds, for altogether 86 H and 88 C atoms. Proton shifts varied between 1.03 ppm to 12.00 ppm and carbon shifts between 7.87 ppm to 209.28 ppm. It was found that the best method for calculating 13C shifts is PBE1PBE/aug-cc-pVDZ with CSGT or IGAIM approaches (mae?=?1.66 ppm), for 1H the best results were obtained with HSEh1PBE, mPW1PW91, PBE1PBE, CAM-B3LYP, and B3PW91 functionals with cc-pVTZ basis set and with CSGT or IGAIM approaches (mae?=?0.28 ppm). We found that often larger basis sets do not give better results for chemical shifts. The best basis sets for calculating 1H and 13C chemical shifts were cc-pVTZ and aug-cc-pVDZ, respectively. CSGT and IGAIM NMR approaches can perform really well and are in most cases better than popular GIAO approach.
Graphical Abstract Mean absolute errors for 1H and 13C chemical shifts and computational times of neutral toluene molecule with aug-cc-pVDZ basis set and CSGT approach
Mechanisms for the activation of water, ammonia, and other small molecules by the PCcarbeneP nickel pincer complex were studied computationally with the aid of density functional theory. The calculation results indicate that the strongly donating, nucleophilic carbene center can engage in a variety of heterolytic splitting of E?H (E=H, C, N, O) bonds, some of which are reversible. The cleavage of E?H bonds across the Ni=C bond represents a new mode of bond activation by ligand cooperativity in nickel pincer complex. On the basis of the calculations, we also demonstrate that reversible H2 activation across the Ir=C bond via the PCcarbeneP iridium pincer complex was observed in the experiments, while other E?H (E=C, N, O) bonds were not activated. Our calculations are in good agreement with experimental observations and could provide new insights into ligand cooperativity in nickel pincer complexes.
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
The aldol reaction in the presence of L-proline acting as an organocatalyst is a well-known example of asymmetric synthesis. Many theoretical and experimental studies have been carried out to probe the mechanism of this reaction. In this work, two levels of density functional theory in the gas phase and DMSO were used to elucidate the best pathways for this reaction, with the enamine and enol considered intermediates and L-proline considered either a reactant or a facilitator. The calculations indicated that both intermediates are formed simultaneously in the reaction medium. Interestingly, the formation of the enamine intermediate predominates in DMSO at room temperature, whereas the enol becomes the predominant intermediate upon the addition of water.
The cooperativity effects of the H-bonding interactions in HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane)???HMX???FA (formamide), HMX???HMX???H2O and HMX???HMX???HMX complexes involving the chair and chair–chair HMX are investigated by using the ONIOM2 (CAM-B3LYP/6–31++G(d,p):PM3) and ONIOM2 (M06-2X/6–31++G(d,p):PM3) methods. The solvent effect of FA or H2O on the cooperativity effect in HMX???HMX???HMX are evaluated by the integral equation formalism polarized continuum model. The results show that the cooperativity and anti-cooperativity effects are not notable in all the systems. Although the effect of solvation on the binding energy of ternary system HMX???HMX???HMX is not large, that on the cooperativity of H-bonds is notable, which leads to the mutually strengthened H-bonding interaction in solution. This is perhaps the reason for the formation of different conformation of HMX in different solvent. Surface electrostatic potential and reduced density gradient are used to reveal the nature of the solvent effect on cooperativity effect in HMX???HMX???HMX.
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.
