Time‐Dependent Density Functional Theory Molecular Dynamics Simulations of Liquid Water Radiolysis

The early stages of the Coulomb explosion of a doubly ionized water molecule immersed in liquid water are investigated with time‐dependent density functional theory molecular dynamics (TD–DFT MD) simulations. Our aim is to verify that the double ionization of one target water molecule leads to the formation of atomic oxygen as a direct consequence of the Coulomb explosion of the molecule. To that end, we used TD–DFT MD simulations in which effective molecular orbitals are propagated in time. These molecular orbitals are constructed as a unitary transformation of maximally localized Wannier orbitals, and the ionization process was obtained by removing two electrons from the molecular orbitals with symmetry 1B₁, 3A₁, 1B₂ and 2A₁ in turn. We show that the doubly charged H₂O²⁺ molecule explodes into its three atomic fragments in less than 4 fs, which leads to the formation of one isolated oxygen atom whatever the ionized molecular orbital. This process is followed by the ultrafast transfer of an electron to the ionized molecule in the first femtosecond. A faster dissociation pattern can be observed when the electrons are removed from the molecular orbitals of the innermost shell. A Bader analysis of the charges carried by the molecules during the dissociation trajectories is also reported.     

A Time‐Dependent Density Functional Theory Investigation on the Origin of Red Chemiluminescence

Is the resonance‐based anionic keto form of oxyluciferin the chemical origin of multicolor bioluminescence? Can it modulate green into red luminescence? There is as yet no definitive answer from experiment or theory. The resonance‐based anionic keto forms of oxyluciferin have been proposed as a cause of multicolor bioluminescence in the firefly. We model the possible structures by adding sodium or ammonium cations and investigating the ground‐ and excited‐state geometries as well as the electronic absorption and emission spectra. A role for the resonance structures is obvious in the gas phase. The absorption and emission spectra of the two structures are quite different—one in the blue and another in the red. The differences in the spectra of the models are small in aqueous solution, with all the absorption and emission spectra in the yellow–green region. The resonance‐based anionic keto form of oxyluciferin may be one origin of the red‐shifted luminescence but is not the exclusive explanation for the variation from green (≈530 nm) to red (≈635 nm). We study the geometries, absorption, and emission spectra of the possible protonated compounds of keto(?1) in the excited states. A new emitter keto(?1)′‐H is considered.     

Theoretical Investigation on the Origin of Yellow‐Green Firefly Bioluminescence by Time‐Dependent Density Functional Theory

The question whether the emitter of yellow‐green firefly bioluminescence is the enol or keto‐constrained form of oxyluciferin (OxyLH₂) still has no definitive answer from experiment or theory. In this study, Arg220, His247, adenosine monophosphate (AMP), Water324, Phe249, Gly343, and Ser349, which make the dominant contributions to color tuning of the fluorescence, are selected to simulate the luciferase (Luc) environment and thus elucidate the origin of firefly bioluminescence. Their respective and compositive effects on OxyLH₂ are considered and the electronic absorption and emission spectra are investigated with B3LYP, B3PW91, and PBE1KCIS methods. Comparing the respective effects in the gas and aqueous phases revealed that the emission transition is prohibited in the gas phase but allowed in the aqueous phase. For the compositive effects, the optimized geometry shows that OxyLH₂ exists in the keto(?1) form when Arg220, His247, AMP, Water324, Phe249, Gly343, and Ser349 are all included in the model. Furthermore, the emission maximum wavelength of keto(?1)+Arg+His+AMP+H₂O+Phe+Gly+Ser is close to the experimental value (560 nm). We conclude that the keto(?1) form of OxyLH₂ is a possible emitter which can produce yellow‐green bioluminescence because of the compositive effects of Arg220, His247, AMP, Water324, Phe249, Gly343, and Ser349 in the luciferase environment. Moreover, AMP may be involved in enolization of the keto(?1) form of OxyLH₂. Water324 is indispensable with respect to the environmental factors around luciferin (LH₂).     

Excited States of Xanthene Analogues: Photofragmentation and Calculations by CC2 and Time‐Dependent Density Functional Theory

Action spectroscopy has emerged as an analytical tool to probe excited states in the gas phase. Although comparison of gas‐phase absorption properties with quantum‐chemical calculations is, in principle, straightforward, popular methods often fail to describe many molecules of interest—such as xanthene analogues. We, therefore, face their nano‐ and picosecond laser‐induced photofragmentation with excited‐state computations by using the CC2 method and time‐dependent density functional theory (TDDFT). Whereas the extracted absorption maxima agree with CC2 predictions, the TDDFT excitation energies are blueshifted. Lowering the amount of Hartree–Fock exchange in the DFT functional can reduce this shift but at the cost of changing the nature of the excited state. Additional bandwidth observed in the photofragmentation spectra is rationalized in terms of multiphoton processes. Observed fragmentation from higher‐lying excited states conforms to intense excited‐to‐excited state transitions calculated with CC2. The CC2 method is thus suitable for the comparison with photofragmentation in xanthene analogues.     

Conjugation‐Promoted Reaction of Open‐Cage Fullerene: A Density Functional Theory Study

Density functional theory calculations are performed to study the addition mechanism of e‐rich moieties such as triethyl phosphite to a carbonyl group on the rim of a fullerene orifice. Three possible reaction channels have been investigated. The obtained results show that the reaction of a carbonyl group on a fullerene orifice with triethyl phosphite most likely proceeds along the classical Abramov reaction; however, the classical product is not stable and is converted into the experimental product. An attack on a fullerene carbonyl carbon will trigger a rearrangement of the phosphate group to the carbonyl oxygen as the conversion transition state is stabilized by fullerene conjugation. This work provides a new insight on the reactivity of open‐cage fullerenes, which may prove helpful in designing new switchable fullerene systems.     

Density Functional Characterization of the Electronic Structure and Visible‐Light Absorption of Cr‐Doped Anatase TiO2

To evaluate the electronic and optical properties of Cr‐doped anatase TiO₂, three possible Cr‐doped TiO₂ models, including Cr at a Ti site (model I), Cr at a Ti site with an oxygen vacancy compensation (model II), and an interstitial Cr site (model III), are studied by means of first principles density functional theory calculations. In model I, the splitting behavior of the Cr 3d states and the insulating properties are successfully depicted by the GGA+U method, from which it is proposed that Cr at a Ti site should exist as Cr⁴⁺ instead of the generally believed Cr³⁺. As a result, the electron transitions between these impurity states, the conduction band (CB), and the valence band (VB), as well as the d–d transitions between occupied and unoccupied Cr 3d states, provide a reasonable explanation for the experimentally observed major and minor absorption bands. In models II and III, the impurity states and associated optical transition processes—as well as the corresponding electron configurations—are examined.     

The Structure of (SCN)x: A Study Using Molecular and Solid‐State Density Functional Theory Calculations

A riddle solved! Despite its simple formula, the structure of the (SCN)_x polymer has remained elusive since its first synthesis in 1929. From energetics as well as NMR chemical shifts, based on DFT calculations, we have strong evidence that it is indeed a tangle of linear chains, made up from N‐linked S₂C₂N five‐membered rings.

     

Theoretical Studies on Muti‐Hydroxyimides as Highly Efficient Catalysts for Aerobic Oxidation

N,N‐dihydroxypyromellitimide (NDHPI) and N,N′,N′′‐trihydroxyisocyanuric acid (THICA) have been recently demonstrated to act as better carbon‐radical‐producing catalysts than the popular N‐hydroxyphthalimide (NHPI). To gain a mature understanding of these particular catalysts, herein their geometrical, electronic, and thermochemical properties, as well as their catalytic activities, have been systemically investigated by a theoretical analysis. It appears that THICA, unlike NDHPI and NHPI, is unsuitable for solvent‐free catalysis or catalysis in aprotic solvents due to its favorable coexistent planar conformer. Besides, the more remarkable catalytic efficiencies of NDHPI and THICA compared to NHPI can be ascribed to the lower barriers and the endothermicity in the H‐abstraction processes by their radicals, especially by their multi‐radicals which show stronger electron‐withdrawing effects. Furthermore, the generation of THICA radicals would be much feasible at high temperature without co‐catalysts. This study provides a new perspective towards the rational design of reactive hydroxyimide organocatalysts for industrial applications.     

Manganese,Iron, Cobalt,and Nickel Oxo‐, Peroxo‐, and Superoxoclusters: A Density Functional Theory Study

The 3d‐transition‐metal dioxo‐, peroxo‐, and superoxoclusters with the general composition MO₂, M(O₂), and MOO (M=Mn, Fe, Co, and Ni) were studied by DFT by the B1LYP functional. The dioxides in their ground states represent the global minima for the M+O₂ system. Both ground‐state dioxides and the lowest‐energy peroxides are in their (d‐only) highest spin states. The ⁶A₁ state of Co(O₂) exceeds the d‐only spin‐multiplicity value (quartet), being nearly isoenergetic with the ⁴A₁ state of Co(O₂). The energy gain on transforming the peroxides to the corresponding dioxides decreases in the order Mn(O₂)>Fe(O₂)>Co(O₂)>Ni(O₂) and varies in the range 0.27–1.8 eV. The dissociation energy to M+O₂ for all studied peroxides is less than 1 eV being the lowest (0.47 eV) for Mn(O₂). The Mn and Fe peroxides need less than 0.3 eV to rupture one of the MO bonds to form the corresponding superoxide. Mn and Fe superoxides are less stable than the corresponding peroxides; the superoxide of Co is more stable than its peroxide, while Ni superoxide is unstable—its energy is above the limit of dissociation to Ni+O₂. According to the electrostatic potential maps, the oxygen atoms in the peroxides are more nucleophilic than those in the dioxides and superoxides, in which the terminal oxygen atom is more nucleophilic than the M‐bonded oxygen atom. This result differs from the expectations based on charge‐distribution analysis.     

New Insights into the Band‐Gap Narrowing of (N,P)‐Codoped TiO2 from Hybrid Density Functional Theory Calculations

The electronic properties of anatase‐TiO₂ codoped by N and P at different concentrations have been investigated via generalized Kohn–Sham theory with the Heyd–Scuseria–Ernzerhof (HSE06) hybrid functional for exchange‐correlation in the context of density functional theory. At high doping concentrations, we find that the high photocatalytic activity of (N, P)‐codoped anatase TiO₂ vis‐à‐vis the N‐monodoped case can be rationalized by a double‐hole‐mediated coupling mechanism [Yin et al., Phys. Rev. Lett. 2011, 106, 066801] via the formation of an effective N? P bond. On the other hand, Ti³⁺ and Ti⁴⁺ ions’ spin double‐exchange results in more substantial gap narrowing for larger separations between N and P atoms. At low doping concentrations, double‐hole‐coupling is dominant, regardless of the N? P distance.     

Local Control Theory in Trajectory Surface Hopping Dynamics Applied to the Excited‐State Proton Transfer of 4‐Hydroxyacridine

The application of local control theory combined with nonadiabatic ab initio molecular dynamics to study the photoinduced intramolecular proton transfer reaction in 4‐hydroxyacridine was investigated. All calculations were performed within the framework of linear‐response time‐dependent density functional theory. The computed pulses revealed important information about the underlying excited‐state nuclear dynamics highlighting the involvement of collective vibrational modes that would normally be neglected in a study performed on model systems constrained to a subset of the full configuration space. This study emphasizes the strengths of local control theory for the design of pulses that can trigger chemical reactions associated with the population of a given molecular excited state. In addition, analysis of the generated pulses can help to shed new light on the photophysics and photochemistry of complex molecular systems.     

A Density Functional Theory Study on the Effect of Zero‐Point Energy Corrections on the Methanation Profile on Fe(100)

The thermodynamics and kinetics of the surface hydrogenation of adsorbed atomic carbon to methane, following the reaction sequence C+4 H?CH+3 H?CH₂+2 H?CH₃+H?CH₄, are studied on Fe(100) by means of density functional theory. An assessment is made on whether the adsorption energies and overall energy profile are affected when zero‐point energy (ZPE) corrections are included. The C, CH and CH₂ species are most stable at the fourfold hollow site, while CH₃ prefers the twofold bridge site. Atomic hydrogen is adsorbed at both the twofold bridge and fourfold hollow sites. Methane is physisorbed on the surface and shows neither orientation nor site preference. It is easily desorbed to the gas phase once formed. The incorporation of ZPE corrections has a very slight, if any, effect on the adsorption energies and does not alter the trends with regards to the most stable adsorption sites. The successive addition of hydrogen to atomic carbon is endothermic up to the addition of the third hydrogen atom resulting in the methyl species, but exothermic in the final hydrogenation step, which leads to methane. The overall methanation reaction is endothermic when starting from atomic carbon and hydrogen on the surface. Zero‐point energy corrections are rarely provided in the literature. Since they are derived from C? H bonds with characteristic vibrations on the order of 2500–3000 cm^?1, the equivalent ZPE of 1/2 hν is on the order of 0.2–0.3 eV and its effect on adsorption energy can in principle be significant. Particularly in reactions between CH_x and H, the ZPE correction is expected to be significant, as additional C? H bonds are formed. In this instance, the methanation reaction energy of +0.77 eV increased to +1.45 eV with the inclusion of ZPE corrections, that is, less favourable. Therefore, it is crucial to include ZPE corrections when reporting reactions involving hydrogen‐containing species.     

Density Functional Theoretical Analysis of the Molecular Structural Effects on Raman Spectra of β‐Carotene and Lycopene

The molecular structural and Raman spectroscopic characteristics of β‐carotene and lycopene are investigated by density functional calculations. The effects of molecular structure and solvent environment on the Raman spectra are analyzed by comparing the calculated and measured results. It is found that the B3LYP/6‐31G(d) method can predict the reasonable result for β‐carotene, but the ν1 Raman activities of lycopene overflow at all the used theoretical methods because of the longer conjugation length. The calculated results indicate that the rotation of β‐rings in β‐carotene impedes the delocalization of π‐electrons, shortens the effective conjugation length, and results in higher frequency and lower activity of the ν1 mode in β‐carotene than lycopene. The measured ν1 bands of β‐carotene and lycopene shift respectively to higher and lower frequencies in solution compared with that in crystals since the crystal packing forces can lead to different conformational variations in the carotenoids molecules. The polarized continuum model theoretical analysis suggests that solvent has slight (significant) effects on the Raman frequencies (intensities) of both carotenoids.     

1，2，5-噻重氮和1，4-二正戊氧基苯环化的自由四氮杂卟啉及其金属镁配合物的结构和性质：密度泛函理论计算

用DFT方法计算分析了1，2，5-噻重氮和1，4-二正戊氧基苯环化的自由四氮杂卟啉及其金属镁配合物的分子和电子结构，理论计算的键参数和单晶结构测定结果一致。进一步对1，2，5-噻重氮和1，4-二正戊氧基苯环化的自由四氮杂卟啉的红外光谱进行了正则坐标分析和光谱模拟，以及用TD-DFT方法对1，2，5-噻重氮和1，4-二正戊氧基苯环化的四氮杂卟啉金属镁配合物的电子吸收光谱进行了分析和谱峰归属，比较了四氮杂卟啉环上取代基的电子性质对四氮杂卟啉衍生物光谱性质的影响。     

HNO与O自由基反应的密度泛函理论研究

用密度泛函理论(DFT)和从头算方法，对HNO与O自由基反应进行了研究。在(U)B3LYP/6-311G^**和(U)B3LYP/aug-cc-pVTZ水平下优化了反应通道上各驻点(反应物、中间体、过渡态及产物)的几何构型。在(U)QCISD/aug-cc-pVTZ水平下计算了各物种的单点能，并对总能量进行了零点能校正。研究结果表明，HNO与O自由基反应过程中存在O → N、O → O和O → H进攻的竞争机制，且存在着多条反应通道。采用过渡态理论计算了600～2 000 K温度范围内3条慢反应通道的速率常数。求得lnk和1/T之间的线性关系。3种通道的阿累尼乌斯指前因子分别为1.469 × 10¹⁰、1.22 × 10¹⁰(1.06 × 10¹⁰)和2.26 × 10¹³。