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.
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.
The adsorption of 20 amino acids (AAs) on the (101) surface of anatase titanium dioxide (TiO2) has been investigated under the scheme of density functional theory. Through the analysis of adsorption geometries, amino group and side chains of AAs have been identified as the major side to adsorb on TiO2, while the carboxyl group prefers to stay outside to avoid the repulsion between negatively charged oxygen from TiO2 and AAs. On the surface, two-coordinated oxygen is the major site to stabilize AAs through O–H interactions. The above conclusion does not change when it is in the aqueous solution based on the calculations with AAs surrounded by explicit water molecules. The above knowledge is helpful in predicting how AAs and even peptides adsorb on inorganic materials.
Graphical abstract The adsorption of 20 amino acids (AAs) on the (101) surface of anatase titanium dioxide (TiO2) has been investigated under the scheme of density functional theory.
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.
The C2 fragmentation energies of the most stable isolated-pentagon-rule (IPR) isomers of the C80 and C82 fullerenes were evaluated with second-order Møller-Plesset (MP2) theory, density-functional theory (DFT) and the semiempirical self-consistent charge density-functional tight-binding (SCC-DFTB) method. Zero-point energy, ionization energy and empirical C2 corrections were included in the calculation of fragmentation energies for comparison with experimental C2 fragmentation energies of the fullerene cations. In the case of the most probable Stone-Wales pathway of C2 fragmentation of C80, the calculated \(D_{0} {\left( {{\text{C}}_{{{\text{80}}}} ^{ + } } \right)}\) agree well with experimental data, whereas in the case of C82 fragmentation, the calculated \(D_{0} {\left( {{\text{C}}_{{{\text{82}}}} ^{ + } } \right)}\) exceed by up to 1.2 eV the experimental ones, which suggests that other IPR isomers may be present in sufficient amounts in experimental samples. Computer-intensive MP2 calculations and DFT calculations with larger basis sets do not yield much improved C2 fragmentation energies, compared to those reported earlier with B3LYP/3-21G. On the other hand, semiempirical approaches such as SCC-DFTB, which are orders of magnitude less intensive, yield satisfactory fragmentation energies for higher fullerenes and may become a method of choice for routine calculations of fullerenes and carbon nanotubes.
Figure C2 fragmentation energies of C80 and C82 fullerenes have been calculated with B3LYP/6-31G* model chemistry, with semiempirical self-consistent-charge density-functional tight-binding (SCC-DFTB) method and with the more rigorous MP2 method. The influence of basis set extension and level of theory on the resulting fragmentation energies is discussed
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.
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.
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.
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.
The 1A1 ground and the first 1B2 excited states of the methylenecyclopropene (triafulvene) are described by localized wave functions, based on 20 structures valence bond structures. The results are compared to CASSCF(4,4) calculations for both the energetics and the dipole moment. Additional calculations with partial electronic delocalization are presented, and it is shown that the dipole moment modification does not correspond to a situation where the antiaromatic situation prevails (with 4n electrons in the cycle). Part of the analysis uses a “trust factor” that helps to decide if a wave function is appropriate to describe a given state. The trust factor compares the VB wave function to the CASSCF’s with their overlap. Finally, the valence bond density is used to produce density maps that illustrate the electron transfer upon excitation.
Graphical Abstract A projector-based method compares CASSCF wave functions to local wave functions, including Lewis structures as shown in the picture. A “trust factor” (τ) is obtained. Both the ground state and the first excited state of the methylenecyclopropene are discussed
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
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.
Cinnamaldehyde (CAL) belongs to the group of aromatic α,β-unsaturated aldehydes; the selective hydrogenation of CAL plays an important role in the fine chemical and pharmaceutical industries. Using Ptn clusters as catalytic models, we studied the selective hydrogenation reaction mechanism for CAL catalyzed by Ptn (n?=?6, 10, 14, 18) clusters by means of B3LYP in density functional theory at the 6-31+?G(d) level (the LanL2DZ extra basis set was used for the Pt atom). The rationality of the transition state was proved by vibration frequency analysis and intrinsic reaction coordinate computation. Moreover, atoms in molecules theory and nature bond orbital theory were applied to discuss the interaction among orbitals and the bonding characteristics. The results indicate that three kinds of products, namely 3-phenylpropyl aldehyde, 3-phenyl allyl alcohol and cinnamyl alcohol, are produced in the selective hydrogenation reaction catalyzed by Ptn clusters; each pathway possesses two reaction channels. Ptn clusters are more likely to catalyze the activation and hydrogenation of the C?=?O bond in CAL molecules, eventually producing cinnamic alcohol, which proves that Ptn clusters have a strong reaction selectivity to catalyze CAL. The reaction selectivity of the catalyzer cluster is closely related to the size of the Ptn cluster, with Pt14 clusters having the greatest reaction selectivity.
Graphical Abstract The reaction mechanism for the selective hydrogenation reaction ofcinnamaldehyde catalyzed by Ptn clusters was studied by densityfunctional theory. The reactionselectivity of cluster catalyzer was concluded to be closely related to the size of Ptn clusters, with Pt14 clusters having the greatest reaction selectivity
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.
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.
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.
This paper inquires the C60 capabilities to contain radio-iodide (131I2) molecules. The encapsulation conditions are investigated applying first principles method to simulate with geometric optimizations and molecular dynamics at 310 K and atmospheric pressure. We find that the n131I2@C60 system, where n?=?1, 2, 3…, is stable if the content does not exceed three molecules of radio-iodide. The application of density functional theory allows us to determine that, the nanocapsules content limit is related with the amount of charge that is transferred from the iodine 131I2 molecules to the carbon atoms in the fullerene surface. The Mulliken population analysis reveals that the excess of charge increases the repulsive forces between atoms and the bond length average in the C60 structure. The weakened bonds easily break and will critically damage the encapsulation properties. Additionally, we test the interaction nanocapsules with different amounts of radioactive iodine diatomic molecules content with calcium atoms, and find that only the fullerene containing one radioactive iodine diatomic molecule was able to interact with up to nine atoms of calcium without disrupting or cracking. Other fullerenes with two and three radio iodine diatomic molecules cannot resist the interaction with a single calcium atom without cracking or being broken.
Carbohydrates have drawn considerable interest from researchers recently due to their affinity for CO2. However, most of the research in this field has focused on peracetylated derivatives. Compared with acetylated carbohydrates, which have already been studied in depth, methyl d-glucopyranoside derivatives are more stable and could have additional applications. Thus, in the present work, ab initio calculations were performed to elucidate the characteristics of the interactions of methylglucoside derivatives with CO2, and to investigate how the binding energy (ΔE) is affected by isomerization or the introduction of various acyl groups. Four methyl d-glucopyranosides (each with two anomers) bearing acetyl, propionyl, butyryl, and isobutyryl moieties, respectively, were designed as substrates, and the 1:1 complexes of a CO2 molecule with each of these sugar substrates were modeled. The results indicate that ΔE is mainly influenced by interaction distance and the number of negatively charged donors or interacting pairs in the complex; the structure of the acyl group present in the substrate is a secondary influence. Except in the case of methyl 2-O-acetyl-d-glucopyranose, the ΔE values of the α- and β-anomers of each methylglucoside were found to be almost the same. Therefore, we would expect the CO2 affinities of the four derivatives studied here to be as strong as or even stronger than that of peracetylated d-glucopyranose.
Graphical Abstract The binding energy between methyl d-glucopyranoside derivatives with various substituted acyl groups and CO2 are evaluated by ab initio calculations. The strong interaction between these methyl dglucopyranoside derivatives and CO2 showed the potential of their application for CO2 capture
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.
