Development of minimized mixing molecular orbital method for designing organic ferromagnets

Predicting the high spin stability of organic radicals correctly for designing organic ferromagnets remains a significant challenge. We have developed a method with an index (L_min) for predicting the high spin stability of conjugated organic radicals at the restricted open‐shell Hartree–Fock level. Unitary transformations were performed for localizing the coefficients of nonbonding molecular orbitals, and subsequently the localized coefficients were used to calculate L_min that indicates the high spin stability of conjugated organic radicals. This method can be combined with the elongation method to treat huge high spin open‐shell systems. Thus, this method is useful for designing organic ferromagnets. © 2015 Wiley Periodicals, Inc.     

Effects of electron–electron interaction on the band structure in polydiacetylenes

The ground‐state band structure of polydiacetylenes is theoretically studied with the extensional Su–Schriffer–Heeger model supplemented by electron–electron interactions. The results show the following. First, the interval of valence bands (conduction bands) increases because of the electron–electron interactions. Second, the effect of the on‐site Coulomb energy (U) is different from that of the nearest neighbor Coulomb repulsion (V); the competition between U and V shows that U makes the bandwidth narrower and the gap broader, whereas V makes the bandwidth broader and the gap narrower. There is a critical value of U/V. Third, the whole band width (E_w) decreases when the U/V ratio is less than 1.0 and increases when the U/V ratio is greater than 1.0 at V = 2.0 eV. Thus, the ground‐state band structure is sensitive to the U/V ratio. The results also show that electron–electron interactions can play an important role in the band structure of polydiacetylenes. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 1656–1661, 2000     

Noncollinear spin density functional theory for spin‐frustrated and spin‐degenerate systems

We have implemented ab initio linear combinations of Gaussian‐type orbital calculations with generalized localized spin density approximation (GLSDA) for a dimer of equilateral H₃ as a model of the noncollinear magnetic clusters. It has been found that the GLSDA solution with the three‐dimensional noncollinear spin structure is, contrary to prior band calculations by other groups, the ground state near the O_h conformation. Further computational results are compared to that of ab initio generalized Hartree–Fock. The difference between them and the influence of the correlation correction were discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

An instability condition for the Hartree–Fock solution of the infinite one‐dimensional system with two‐crossing bands. I. Singlet‐instability check of metallic carbon nanotube

An instability condition is derived for the Hartree–Fock solution so that it can be applied to the system in which the highest occupied and the lowest unoccupied bands cross at the in‐between point in the Brillouin zone. The instability check developed here is further applied to a metallic single‐walled carbon nanotube having the two‐crossing bands toward prediction of its instability. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 574–582, 2000     

Quantum chemical investigation of lanthanum doping effects in tetragonal and rhombohedral PZT materials

Structural and electronic properties of lead–zirconate–titanate (PZT) materials doped with a lanthanum (La) impurity are studied using a quantum‐chemical approach based on the Hartree–Fock theory. Performed geometry optimization in the defective crystals shows that the atomic movements are predominantly outward with respect to the impurity position. It is found that the La impurity enhances a covalent character in the chemical bonding between the Ti and O atoms, as well as the Zr and O atoms situated in the neighborhood of the defect despite the fact that the La‐O interaction remains purely ionic. The occurrence of local energy levels within the band gap of the material is analyzed in light of the available experimental data on La concentration influence upon dielectric and piezoelectric properties in these crystals. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Seminumerical calculation of the Hartree-Fock exchange matrix: application to two-component procedures and efficient evaluation of local hybrid density functionals

A two-component extension of the seminumerical procedure for the calculation of the Hartree-Fock (HF) exchange matrix recently presented by Neese et al. (Chem Phys 2009, 356, 98) was implemented into the program system TURBOMOLE. It is demonstrated that this allows for efficient self-consistent treatment of spin-orbit coupling at HF and hybrid density functional theory level. One-component HF calculations were performed to study the accuracy of integration grids and the exploitation of the molecular point group symmetry. The efficiency was tested, and for one-component calculations compared to the implementation realized by Neese. It was further demonstrated that local hybrid density functionals can be evaluated with this technique. The "prototype" of this class of functionals, Lh-BLYP, was applied to an organic molecule with more than 150 atoms.     

A quantum‐chemical study of phosphor impurity in BaTiO3 crystal

Structural and electronic effects produced by a P impurity in the cubic and tetragonal BaTiO₃ lattices are investigated using a simple quantum chemical computer code based on the Hartree–Fock methodology. The obtained atomic displacements due to the defect present in otherwise pure crystal are mainly toward the impurity atom, thus reducing the interatomic distances within the defective region. It is also found that the phosphor produces some redistribution of electron density from the defect‐neighboring atoms toward the chemical bonds thus diminishing the charges on atoms. We also observe a local energy level in the band‐gap of material being composed mainly of P 3s atomic orbital. The level finds itself close to the top of the upper valence band, in no case contributing into the n‐type conductivity in BaTiO₃. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Excited states of OsO4: A comprehensive time‐dependent relativistic density functional theory study

A large number of scalar as well as spinor excited states of OsO₄, in the experimentally accessible energy range of 3–11 eV, have been captured by time‐dependent relativistic density functional linear response theory based on an exact two‐component Hamiltonian resulting from the symmetrized elimination of the small component. The results are grossly in good agreement with those by the singles and doubles coupled‐cluster linear response theory in conjunction with relativistic effective core potentials. The simulated‐excitation spectrum is also in line with the available experiment. Furthermore, combined with detailed analysis of the excited states, the nature of the observed optical transitions is clearly elucidated. It is found that a few scalar states of ³T₁ and ³T₂ symmetries are split significantly by the spin‐orbit coupling. The possible source for the substantial spin‐orbit splittings of ligand molecular orbitals is carefully examined, leading to a new interpretation on the primary valence photoelectron ionization spectrum of OsO₄. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Generalized spin orbital GW theory for spin‐frustrated and spin‐degenerate systems

We proposed the ab initio linear combination of Gaussian type orbital (LCGTO) generalized spin orbital GW (GSO–GW) method and calculated triangular hydrogen molecules as models of the noncolinear magnetic clusters. A remarkable improvement of ionized potentials (IPs) by the GW procedure for GHFS solutions is observed in comparison with calculational results by full CI. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 369–374, 2001     

Relativistic Gaussian basis sets for molecular calculations: Fully optimized single‐family exponent basis sets for H?Hg

Relativistic single‐family exponent Gaussian basis sets for molecular calculations are presented for the 80 atoms ₁H through ₈₀Hg. The exponent parameters shared by Gaussian basis functions of all symmetry species are fully optimized. Two nucleus models of uniformly charged sphere and Gaussian charge distribution are considered and two kinds of basis sets are generated accordingly. The total energy errors are less than 2 mhartree in any atoms. Some of the present basis sets include small variational collapse (or prolapse), but test calculations show that they could be reliably applied to molecular calculations. © 2005 Wiley Periodicals, Inc. J Comput Chem 27: 48–52, 2006     

6‐31G* basis set for third‐row atoms

Medium basis sets based upon contractions of Gaussian primitives are developed for the third‐row elements Ga through Kr. The basis functions generalize the 6‐31G and 6‐31G* sets commonly used for atoms up to Ar. A reexamination of the 6‐31G* basis set for K and Ca developed earlier leads to the inclusion of 3d orbitals into the valence space for these atoms. Now the 6‐31G basis for the whole third‐row K through Kr has six primitive Gaussians for 1s, 2s, 2p, 3s, and 3p orbitals, and a split‐valence pair of three and one primitives for valence orbitals, which are 4s, 4p, and 3d. The nature of the polarization functions for third‐row atoms is reexamined as well. The polarization functions for K, Ca, and Ga through Kr are single set of Cartesian d‐type primitives. The polarization functions for transition metals are defined to be a single 7f set of uncontracted primitives. Comparison with experimental data shows good agreement with bond lengths and angles for representative vapor‐phase metal complexes. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 976–984, 2001     

An improved generator coordinate Hartree–Fock method applied to the choice of contracted Gaussian basis sets for first‐row diatomic molecules

Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first‐row atoms, generated with an improved generator coordinate Hartree–Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B₂, C₂, BeO, CN⁻, LiF, N₂, CO, BF, NO⁺, O₂, and F₂. At the Hartree–Fock (HP), second‐order Møller–Plesset (MP2), fourth‐order Møller–Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between our total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 15–23, 2000     

Calculation of quasiparticle energy of molecular systems by the GW method

The quasiparticle energy of the H₂ molecule is calculated by using the GW method, in which the self‐energy operator fully depends on the frequency. The initial Green function G₀ is constructed from the wave function obtained by the Hartree–Fock approximation (HFA) and local density approximation (LDA) in the framework of the density functional theory (DFT). From the results obtained we have shown that the wave function from the DFT–LDA is more effective than that from the HFA for G₀. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 348–353, 2001     

Multiplicity,instability, and SCF convergence problems in Hartree–Fock solutions

We present a study of the instability and convergence of Hartree–Fock (HF) ab initio solutions for the diatomic systems H₂, LiH, CH, C₂, and N₂. In our study, we consider real molecular orbitals (MOs) and analyze the classes of single‐determinant functions associated to Hartree–Fock–Roothaan (HFR) and Hartree–Fock–Pople–Nesbet (HFPN) equations. To determine the multiple HF solutions, we used either an SCF iterative procedure with aufbau and non‐aufbau ordering rules or the algebraic method (AM). Stability conditions were determined using TICS and ASDW stability matrices, derived from the maximum and minimum method of functions (MMF). We examined the relationship between pure SCF convergence criterion with the aufbau ordering rule, and the classification of the HF solution as an extremum point in its respective class of functions. Our results show that (i) in a pure converged SCF calculation, with the aufbau ordering rule, the solutions are not necessarily classified as a minimum of the HF functional with respect to the TICS or ASDW classes of solutions, and (ii) for all studied systems, we obtained local minimum points associated only with the aufbau rule and the solutions of lower energies. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 600–610, 2000     

Theoretical study of one‐carbon unit transfer between methyl‐AICA and N1‐methyl‐N1‐acryloyl‐formamide

The mechanism of one‐carbon unit transfer between 1‐methyl‐5‐amino‐4‐carboxamide imidazole (M‐AICA) and N¹‐methyl‐N¹‐acryloyl‐formamide (the model molecule of 10‐f‐H4F) is investigated by the Hartree–Fock and DFT methods, respectively, at the 6‐31G* basis level. There are two different channels for the proton transfer, resulting in two reaction pathways with different properties. The results indicate that both channels can complete the reaction, but path a is slightly favored due to its lower active energy barrier. Furthermore, the influence of 4‐carboxamindde in M‐AICA is also discussed. This group can stabilize the reactant and intermediates, and reduce the active energy barrier through the intermolecular hydrogen bond. The intermolecular hydrogen bond results in an enlarged conjugation system and makes the transition states more stable. Our results are in agreement with experiments. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

H‐doped PbTiO3: Structure and electronic properties

The geometry and electronic properties of the interstitial H atom in the tetragonal PbTiO₃ crystal have been studied using an advanced quantum chemical computer code developed for the modeling of crystals. The inserted H atom was found to bind to one of the O atoms and to form the hydroxyl, O? H group, with the inter‐atomic distance equal to 0.93 Å and 1.00 Å for the hydroxyls containing O atom in the dimerized and nondimerized Ti? O? Ti chains, respectively. Atomic displacements in the vicinity of O? H complex are calculated and analyzed in relation to the H‐produced changes upon the atomic charges in defective region. The role of H impurity on the ferroelectric polarization in the tetragonal PbTiO₃ is discussed in terms of the results obtained in our research and those presented in the other studies on this subject. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007