Search

A Kohn--Sham-type computational scheme capable of treating systems with strong nondynamic correlation is examined. The experiment results of REKS calculation are compared with the available Kohn--Sham solutions for cases in which the exact density is known, as well as with results of conventional multireference ab initio methods.

Research is presented applying REKS-type (Filatov, M.; Shaik, S. Chem. Phys. Lett. 1999) density functional calculations to study the molecular and electronic structure of the triplet and the lowest singlet state of tetramethyleneethane diradical.

Research is done on measurements at the (14,10)CASSCF/6-31G** and the MR-(S)DCI/cc-pVTZ levels, which are employed by studying three possible modes of reaction between dioxygen and dihydrogen molecules at the ground triplet state and excited singlet state of O2. In conclusion, the measurements revealed that the lowest singlet state of O2 has very high barriers for reaction and therefore cannot mediate a facile oxidation of H2 in contrast to transition metal oxenide cation catalysts and mono-oxygenase enzymes and the fundamental difference is explained.

In the past, the design of light-driven single molecule rotary motors has been mainly guided by the modification of their ground-state conformational properties. Further progress in this field is thus likely to be achieved through a detailed understanding of light-induced dynamics of the system and the ways of modulating it by introducing chemical modifications. In the present theoretical work, the analysis of model organic chromophores and synthesized rotary motors is used for rationalizing the effect of electron-withdrawing heteroatoms (such as a cationic nitrogen) on the topography and branching plane of mechanistically relevant conical intersections. Such an analysis reveals how the character of rotary motion could be changed from a precessional motion to an axial rotational motion. These concepts are then used to design and build quantum chemical models of three distinct types of Schiff base rotary motors. One of these models, featuring the synthetically viable indanylidenepyrroline framework, has conical intersection structures consistent with an axial rotation not hindered by ground-state conformational barriers. It is expected that this type of motor should be capable of funneling the photon energy into specific rotary modes, thus achieving photoisomerization quantum efficiencies comparable to those seen in visual pigments. [ABSTRACT FROM AUTHOR]

The mechanisms of H2 oxidation by FeO+ were investigated. Density functional theory calculations using the B3LYP, BP86 and FT97 functionals with an extended basis set were used to examine the three mechanisms: addition-elimination, 'rebound' and oxene-insertion. The results showed that the addition-elimination mechanism is dominant at low temperature while the 'rebound' mechanism is dominant at high temperature. The intermediate temperatue was characterized by two-state reactivity where the two mechanism compete.

A new adaptive algorithm for penalty function optimization for minimum-energy three-states conical intersections (ME3CI) is suggested. The new algorithm differs from the original penalty function algorithm by (a) removing the redundancy in the target function, (b) using an adaptive increment for the penalty function weighting factor, and (c) using tighter convergence criteria for the energy gap. The latter was introduced to guarantee convergence to a true conical intersection rather than to a narrowly avoided crossing geometry. The new algorithm was tested in the optimization of the ME3CI geometries in butadiene and malonaldehyde, where all of the previously found true ME3CI geometries were recovered. The previously found butadiene's CI 3/2/1 turned out to be a narrowly avoided crossing. For butadiene, seven new ME3CI geometries have been located. Because of the removal of the redundancy and the use of the adaptive weighting factor, the convergence rate of the new algorithm is noticeably improved as compared to that of the previously proposed penalty function algorithm. The application to malonaldehyde and butadiene demonstrates that the three-state conical intersections may be more abundant and hence more involved in the photochemistry than previously thought. The recently developed mixed-reference spin flip (MRSF)-TDDFT method yields ME3CI geometries and relative energies quantitatively consistent with the previously reported calculations at a much reduced computational cost.

The computation of electron attachment energies (electron affinities) was implemented in connection with an ensemble density functional theory method, the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method. With the use of the extended Koopmans' theorem, the electron affinities and the respective Dyson orbitals are obtained directly for the neutral molecule, thus avoiding the necessity to compute the ionized system. Together with the EKT-SSR (extended Koopmans' theorem-SSR) method for ionization potentials, which was developed earlier, EKT-SSR for electron affinities completes the implementation of the EKT-SSR formalism, which can now be used for obtaining electron detachment as well as the electron attachment energies of molecules in the ground and excited electronic states. The extended EKT-SSR method was tested in the calculation of several closed-shell molecules. For the molecules in the ground states, the EKT-SSR energies of Dyson's orbitals are virtually identical to the energies of the unoccupied orbitals in the usual single-reference spin-restricted Kohn-Sham calculations. For the molecules in the excited states, EKT-SSR predicts an increase of the most positive electron affinity by approximately the amount of the vertical excitation energy. The electron affinities of a number of diradicals were calculated with EKT-SSR and compared with the available experimental data. With the use of a standard density functional (BH&HLYP), the EKT-SSR electron affinities deviate on average by ca. 0.2 eV from the experimental data. It is expected that the agreement with the experiment can be improved by designing density functionals parametrized for ionization energies.

The mixed-reference spin-flip (MRSF) time-dependent density functional theory (TDDFT) method eliminates the notorious spin contamination of SF-TDDFT, thus enabling identification of states of proper spin-symmetry for automatic geometry optimization and molecular dynamics simulations. Here, we analyze and optimize the MRSF-TDDFT in the calculations of the vertical excitation energies (VEEs) and the singlet-triplet (ST) gaps. The dependence of the obtained VEEs and ST gaps on the intrinsic parameters of the MRSF-TDDFT method is investigated, and prescriptions for the proper use of the method are formulated. For VEEs, MRSF-TDDFT displays similar or better accuracy than SF-TDDFT (ca. 0.5 eV), while considerably outperforming the LR-TDDFT for the ST gaps. As a result, a new functional of STG1X (dubbed here), especially for ST gaps is suggested on the basis of splitting between the components of the atomic multiplets.

A quantum theory of atoms in molecules (QTAIM) and stress tensor analysis was applied to analyze intramolecular interactions influencing the photoisomerization dynamics of a light-driven rotary molecular motor. For selected nonadiabatic molecular dynamics trajectories characterized by markedly different S 1 state lifetimes, the electron densities were obtained using the ensemble density functional theory method. The analysis revealed that torsional motion of the molecular motor blades from the Franck-Condon point to the S 1 energy minimum and the S 1 /S 0 conical intersection is controlled by two factors: greater numbers of intramolecular bonds before the hop-time and unusually strongly coupled bonds between the atoms of the rotor and the stator blades. This results in the effective stalling of the progress along the torsional path for an extended period of time. This finding suggests a possibility of chemical tuning of the speed of photoisomerization of molecular motors and related molecular switches by reshaping their molecular backbones to decrease or increase the degree of coupling and numbers of intramolecular bond critical points as revealed by the QTAIM/stress tensor analysis of the electron density. Additionally, the stress tensor scalar and vector analysis was found to provide new methods to follow the trajectories, and from this, new insight was gained into the behavior of the S 1 state in the vicinity of the conical intersection.

Various types of interactions between halogen (X) and π moiety (X-π interaction) including halogen bonding play important roles in forming the structures of biological, supramolecular, and nanomaterial systems containing halogens and aromatic rings. Furthermore, halogen molecules such as X 2 and CX 4 (X = Cl/Br) can be intercalated in graphite and bilayer graphene for doping and graphene functionalization/modification. Due to the X-π interactions, though recently highly studied, their structures are still hardly predictable. Here, using the coupled-cluster with single, double, and noniterative triple excitations (CCSD(T)), the Møller-Plesset second-order perturbation theory (MP2), and various flavors of density functional theory (DFT) methods, we study complexes of benzene (Bz) with halogen-containing molecules X 2 and CX 4 (X = Cl/Br) and analyze various components of the interaction energy using symmetry adapted perturbation theory (SAPT). As for the lowest energy conformers (S1), X 2 -Bz is found to have the T-shaped structure where the electropositive X atom-end of X 2 is pointing to the electronegative midpoint of CC bond of the Bz ring, and CX 4 -Bz has the stacked structure. In addition to this CX 4 -Bz (S1), other low energy conformers of X 2 -Bz (S2/S3) and CX 4 -Bz (S2) are stabilized primarily by the dispersion interaction, whereas the electrostatic interaction is substantial. Most of the density functionals show noticeable deviations from the CCSD(T) complete basis set (CBS) limit binding energies, especially in the case of strongly halogen-bonded conformers of X 2 -Bz (S1), whereas the deviations are relatively small for CX 4 -Bz where the dispersion is more important. The halogen bond shows highly anisotropic electron density around halogen atoms and the DFT results are very sensitive to basis set. The unsatisfactory performance of many density functionals could be mainly due to less accurate exchange. This is evidenced from the good performance by the dispersion corrected hybrid and double hybrid functionals. B2GP-PLYP-D3 and PBE0-TS(Tkatchenko-Scheffler)/D3 are well suited to describe the X-π interactions adequately, close to the CCSD(T)/CBS binding energies (within ∼1 kJ/mol). This understanding would be useful to study diverse X-π interaction driven structures such as halogen containing compounds intercalated between 2-dimensional layers.

Based on the normalized elimination of the small component relativistic formalism, a new approach to the calculation of hyperfine structure parameters of paramagnetic molecules is developed and implemented. The new method is tested in the calculation of the isotropic hyperfine structure constant for a series of open-shell molecules containing mercury. The results of calculations carried out in connection with ab initio methods of increasing complexity demonstrate the high accuracy of the formalism developed. In view of its computational simplicity, the new approach provides the basis for an efficient and accurate calculation of the HFS parameters of large molecules.