Broken-symmetry unrestricted hybrid density functional calculations on nickel dimer and nickel hydride

In the present work we investigate the adequacy of broken-symmetry unrestricted density functional theory for constructing the potential energy curve of nickel dimer and nickel hydride, as a model for larger bare and hydrogenated nickel cluster calculations. We use three hybrid functionals: the popular B3LYP, Becke's newest optimized functional Becke98, and the simple FSLYP functional (50% Hartree-Fock and 50% Slater exchange and LYP gradient-corrected correlation functional) with two basis sets: all-electron (AE) Wachters+f basis set and Stuttgart RSC effective core potential (ECP) and basis set. We find that, overall, the best agreement with experiment, comparable to that of the high-level CASPT2, is obtained with B3LYP/AE, closely followed by Becke98/AE and Becke98/ECP. FSLYP/AE and B3LYP/ECP give slightly worse agreement with experiment, and FSLYP/ECP is the only method among the ones we studied that gives an unacceptably large error, underestimating the dissociation energy of Ni(2) by 28%, and being in the largest disagreement with the experiment and the other theoretical predictions. We also find that for Ni(2), the spin projection for the broken-symmetry unrestricted singlet states changes the ordering of the states, but the splittings are less than 10 meV. All our calculations predict a deltadelta-hole ground state for Ni(2) and delta-hole ground state for NiH. Upon spin projection of the singlet state of Ni(2), almost all of our calculations: Becke98 and FSLYP both AE and ECP and B3LYP/AE predict (1)(d(A)(x(2)-y(2)d(B)(x(2)-y(2)) or (1)(d(A)(xy) (d)(B)(xy)) ground state, which is a mixture of (1)Sigma(g) (+) and (1)Gamma(g). B3LYP/ECP predicts a (3)(d(A)(x(2)-y(2))d(B)(xy) (mixture of (3)Sigma(g) (-) and (3)Gamma(u)) ground state virtually degenerate with the (1)(d(A)(x(2)-y(2)d(B)(x)(2)-y(2)/(1)(d(A)(xy)D(B)(xy) state. The doublet delta-hole ground state of NiH predicted by all our calculations is in agreement with the experimentally predicted (2)Delta ground state. For Ni(2), all our results are consistent with the experimentally predicted ground state of 0(g) (+) (a mixture of (1)Sigma(g) (+) and (3)Sigma(g) (-)) or 0(u) (-) (a mixture of (1)Sigma(u) (-) and (3)Sigma(u) (+)).     

Generalized approximation to the reaction path: the formic acid dimer case

A set of mass-weighted internal coordinates was derived and applied to the double proton transfer reaction in the formic acid dimer (FAD). The coordinate set was obtained starting from the Hirschfelder "mobile" by an optimization procedure consisting of a sequence of kinematic rotations. In FAD, the optimization procedure leads to three coordinates that do change significantly along the reaction path. These coordinates span the reaction space, whereas the remaining modes are treated in a harmonic approximation. The effect that the dimer dissociative motion has on the ground and excited vibrational states dynamics was explored. In the frequency region corresponding to the symmetric OH-stretch vibration four doublets have been identified with splittings of 2.76, 0.07, 0.60, and 4.03 cm(-1).     

Compliant fields for molecular interactions: Water dimer and formic acid dimer

The redundancy-free internal valence compliance constants of open-chain water dimer and formic acid cyclic dimer have been determined by the combined use of the CNDO /Force method and the compliance constant formalism. The final compliant fields of these dimers have been refined with the help of experimental frequency data.     

Geometries and stabilities of l-ascorbic acid dimer and its derivatives: a computational investigation by density functional methods

Geometries, relative stabilities, and hydrogen bonds of l-ascorbic acid (LAA) and d-Erythroascorbate (DEAA) dimers as well as their S- and Se-substituted isomers in gas phase and water solvent are studied using density functional method. Furthermore, the hydrogen bond lengths in LAA and DEAA dimers are generally increased along with the binding dissociation energy of the dimers being decreased as apex O atoms in the five-membered C5 rings of LAA and DEAA dimers are substituted by S and Se atoms in gas phase and water solvent. Interestingly, one LAA dimer and its S- or Se-substituted isomer with four hydrogen bonds in gas and water solvent are the three-centers structures. In addition, the chemical bonding and charge distributions of all the dimers are discussed. A good agreement with available experimental results is reached.     

An experimental and theoretical study of molecular structure and vibrational spectra of 3- and 4-pyridineboronic acid molecules by density functional theory calculations

The experimental and theoretical vibrational spectra of 3- and 4-pyridineboronic acids (abbreviated as p3 and p4) were studied. The Fourier transform Raman and Fourier transform infrared spectra of p3 and p4 molecules were recorded in the solid phase. The structural and spectroscopic analysis of the p3 and p4 acids were made by using density functional harmonic calculations. Both p3 and p4 only one form was most stable using B3LYP level with the 6-31G(d), 6-31G(d,p), 6-311G(d), 6-311G(d,p) and 6-311++G(d,p) basis sets. Selected experimental bands were assigned based on the scaled theoretical wavenumbers. Finally, geometric parameters, infrared and Raman bands and intensities were compared with experimental data of the molecules.     

A density functional theory study for the hydrogen-bonded nucleic acid base pair: cytosine dimer

Theoretical investigation for the geometric and energetic properties, rotational constants, harmonic vibrational frequencies, and binding energies of nucleic acid base pair, cytosine dimer, are carried out by using the density functional theory method. The dimer structures resulting from both the keto and the enol (cis/trans) tautomers are investigated in the present study. Various isomers are considered to find the stable structures of the cytosine dimer. The planar cytosine dimer, K-K3 with C2h symmetry, resulting from nonplanar keto tautomers, is found to be thermodynamically most stable out of the four different stable isomers and having the highest binding energy value, 19.51 kcal/mol (including basis set superposition error correction). The vibrational frequency analysis also suggests a red shift of 367.97 cm(-1) for the hydrogen-bonding K-K3 symmetric dimer with two hydrogen bond lengths, each of length 1.913 angstroms. Moreover, charge distribution (ChelpG charges), Laplacian electronic density distribution, and the dimerization equilibrium for the most stable dimer, K-K3, have also been investigated using the same method and the basis set.     

An examination of density functional theories on isomerization energy calculations of organic molecules

Long-range corrected (LC) density functional theories (DFTs) were applied to the isomerization energy calculations of organic molecules to make clear why conventional DFTs including B3LYP have given poor isomerization reaction energies. Combining with local response dispersion (LRD) method, we performed LC-DFT calculations for the benchmark set of isomerization reactions. Consequently, we found that LC-DFT?+?LRD methods give accurate reaction energies equivalent to up-to-date DFTs containing many semi-empirical parameters. This result indicates that long-range exchange and intramolecular dispersion correlation interactions, which have been neglected in conventional DFTs, play prominent roles in isomerization reactions. However, we also found that these interactions are not sufficient to give accurate isomerization energies especially for cyclization reactions. Considering that Gaussian-attenuated LC-DFTs (LCgau-DFTs) give better isomerization reaction energies than LC-DFTs, we suggested that the isomerization energies will be further improved by correcting the short-range part of exchange functionals in DFT with keeping the whole long-range exchange interactions.     

Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules

We present the implementation of a variational finite element solver in the HelFEM program for benchmark calculations on diatomic systems. A basis set of the form is used, where (μ, ν, φ) are transformed prolate spheroidal coordinates, B _n(μ) are finite element shape functions, and are spherical harmonics. The basis set allows for an arbitrary level of accuracy in calculations on diatomic molecules, which can be performed at present with either nonrelativistic Hartree-Fock (HF) or density functional (DF) theory. Hundreds of DFs at the local spin density approximation (LDA), generalized gradient approximation (GGA), and the meta-GGA level can be used through an interface with the Libxc library; meta-GGA and hybrid DFs are not available in other fully numerical diatomic program packages. Finite electric fields are also supported in HelFEM , enabling access to electric properties. We introduce a powerful tool for adaptively choosing the basis set by using the core Hamiltonian as a proxy for its completeness. HelFEM and the novel basis set procedure are demonstrated by reproducing the restricted open-shell HF limit energies of 68 diatomic molecules from the first to the fourth period with excellent agreement with literature values, despite requiring orders of magnitude fewer parameters for the wave function. Then, the electric properties of the BH and N₂ molecules under finite field are studied, again yielding excellent agreement with previous HF limit values for energies, dipole moments, and dipole polarizabilities, again with much more compact wave functions than what were needed for the literature references. Finally, HF, LDA, GGA, and meta-GGA calculations of the atomization energy of N₂ are performed, demonstrating the superb accuracy of the present approach.     

Symmetric double proton tunneling in formic acid dimer: a diabatic basis approach

A model of double proton tunneling in formic acid dimer is developed using a reaction surface Hamiltonian. The surface includes the symmetric OH stretch plus the in-plane stretch and bend interdimer vibrations. The surface Hamiltonian is coupled to a bath of five A1g and B3g normal modes obtained at the D2h transition state structure. Eigenstates are calculated using Davidson and block-Davidson iterative methods. Strong mode specific effects are found in the tunneling splittings for the reaction surface, where splittings are enhanced upon excitation of the interdimer bend motion. The results are interpreted within the framework of a diabatic representation of reaction surface modes. The splitting patterns observed for the reaction surface eigenstates are only slightly modified upon coupling to the bath states. Splitting patterns for the bath states are also determined. It is found that predicting these splittings is greatly complicated by subtle mixings with the inter-dimer bend states.     

Exploring odd-even effects of simple oligomer-like DRCNnT series: a study based on density functional theory/time-dependent density functional theory calculations

Odd-even effects of short-circuit current density and power conversion efficiency (PCE) are an interesting phenomenon in some organic solar cells. Although some explanations have been given, why they behave in such a way is still an open question. In the present work, we investigate a set of acceptor-donor-acceptor simple oligomer-like small molecules, named the DRCNnT (n = 5-9) series, to give an insight into this phenomenon because the solar cells based on them have high PCE (up to 10.08%) and show strong odd-even effects in experiments. By modeling the DRCNnT series and using density functional theory, we have studied the ground-state electronic structures of the DRCNnT (n = 5-9) series in condensed phase. The calculated results reproduce the experimental trends well. Furthermore, we find that the exciton-binding energies of the DRCNnT series may be one of the key parameters to explain this phenomenon because they also show odd-even effects. In addition, by studying the effects of alkyl branch and terminal group on odd-even effects of dipole moment, we find that eliminating one or two alkyl branches does not break the odd-even effects of dipole moments, but eliminating one or two terminal groups does. Finally, we conclude that removing one alkyl branch close to the terminal group of DRCN5T can decrease highest occupied molecular orbital (HOMO) energy (thus increasing open circuit voltage) and increase dipole moment (thus enhancing charge separation and short-circuit current). This could be a new and simple method to increase the PCE of DRCN5T-based solar cells.     

Hydrogen bonding versus coordination of adsorbate molecules on Ti-silicalites: a density functional theory study

The present study discusses the results of theoretical calculations obtained at the B3LYP/ 6-31G level on the structural, electronic, and energetic properties of Ti-silicalites. Particularly, the relevance of 5T cluster models, either H- or OH-terminated, in large-scale calculations has been critically considered. It was shown that an open surface structure with one OH group and a closed-bulk structure with no bonded OH group at the Ti site are responsible for the observed UV-vis properties of Ti-silicalite materials. Both water and methanol can preferably interact with Ti-silicalites through the H-bonding mechanism, while ammonia can form either H-bonded or coordination complexes. The calculations support the existence of highly dispersed Ti sites in a tetrahedral environment only in Ti-silicalites because an increase in the coordination number of the Ti site by next-neighbor lattice oxygens is the energetically less favorable process.     

Dissociation of sulfur oxoacids by two water molecules studied using ab initio and density functional theory calculations

Using ab initio [SCS‐MP2 and CCSD(T)] and density functional theory (M062X) calculations, we have studied the geometries and energies of sulfur oxoacids H₂S_mO₆ (m = 2–4) and their monohydrated and dihydrated clusters. When including the results from previously reported disulfuric acid (H₂S₂O₇) cases, the gas phase acidity is ordered as H₂S₂O₆ < H₂S₃O₆ < H₂S₂O₇ < H₂S₄O₆. The intramolecular H‐bonding, which may indicate the degree of structural flexibility in this molecular series, is an important factor for the order of the gas phase acidity. All these sulfur oxoacids show dissociated (or deprotonated) geometries with only two water molecules, although the energies of the dissociated conformers are ranked differently. All of the dissociated conformers form a unique H‐bonding network structure in which the protonated first water (H₃O⁺) is triply H‐bonded to each oxygen atom of two SO₃ moieties as well as the second water, which in turn is H‐bonded to a SO₃ moiety. H₂S₃O₆ has the best molecular flexibility for adopting such an H‐bonding network structure, and thereby all the low‐lying conformers of H₂S₃O₆(H₂O)₂ are dissociated. In contrast, the least flexible H₂S₂O₆ forms such a structure with a high strain, and dissociation of H₂S₂O₆(H₂O)₂ is found from the third lowest conformer. Although the gas phase acidity of H₂S₄O₆ is the highest in this series, the lowest dissociated conformer and the lowest undissociated conformer of H₂S₄O₆(H₂O)₂ are very close in energy. This is because forming the H‐bonding network structure is somewhat difficult due to the large distance between the two SO₃ moieties.     

Mesophase formation in a system of top-shaped hard molecules: density functional theory and Monte Carlo simulation

We present the phase diagram of a system of mesogenic top-shaped molecules based on the Parsons-Lee density functional theory and Monte Carlo simulation. The molecules are modeled as a hard spherocylinder with a hard sphere embedded in its center. The stability of five different phases is studied, namely, isotropic, nematic, smectic A, smectic C, and columnar phases. The positionally ordered phases are investigated only for the case of parallel alignment. It is found that the central spherical unit destabilizes the nematic with respect to the isotropic phase, while increasing the length of the cylinder has the opposite effect. Also, the central hard sphere has a strong destabilizing effect on the smectic A phase, due the inefficient packing of the molecules into layers. For large hard sphere units the smectic A phase is completely replaced by a smectic C structure. The columnar phase is first stabilized with increasing diameter of the central unit, but for very large hard sphere units it becomes less stable again. The density functional results are in good agreement with the simulations.     

Modeling proton transfer in imidazole-like dimers: a density functional theory study

A detailed theoretical study of proton transfer reaction in protonated imidazole, 1,2,3-triazole, and tetrazole dimers, the basic components of polymeric membrane used in proton exchange membranes fuel cells, has been carried out. In particular, several approaches based on density functional theory have been considered and their results compared with those provided by post-HF methods. From a computational point of view, these molecules appear to be a very challenging playground also for robust and recent functionals. Indeed none of the considered approaches provide results in close agreement with the reference post-HF data and a combined BMK//B3LYP model seems the only approach able to reproduce both the energetic and the structural features. From a chemical point of view, two new mechanisms of proton transfer in tetrazole dimers have been investigated and found to be more favorable than that previously hypothesized in literature. At the same time, the theoretical results show a direct connection between the obtained proton transfer barrier and the charge localized on the transferred hydrogen.     

Modeling a halogen dance reaction mechanism: A density functional theory study

Since the discovery of the halogen dance (HD) reaction more than 60 years ago, numerous insights into the mechanism have been unveiled. To date however, the reaction has not been investigated from a theoretical perspective. Density functional theory (DFT) was used to model the potential energy surface linking the starting reagents to the lithiated products for each step in the mechanism using a thiophene substrate. It was found that the lithium‐halogen exchange mechanism is critical to understand the HD mechanism in detail and yielded the knowledge that S_N2 transition states (TS) are favored over the four‐center type for the lithium‐bromine exchange steps. The overall driving force for the HD is thermodynamics, while the kinetic factors tightly control the reaction path through temperature. The S_N2 lithium‐bromide TS are barrierless, except the second, which is the limiting step. Finally, the model for the HD is discovered to be a pseudo‐clock type, due to a highly favorable bromide catalysis step and the reformation of 2‐bromothiophene. © 2016 Wiley Periodicals, Inc.     

Substituted 3-phenylpropenoates and related analogs: electron ionization mass spectral fragmentation and density functional theory calculations

Analysis of ethyl 3-(2-chlorophenyl)propenoate by electron ionization mass spectrometry showed the distinct loss of an ortho chlorine. To characterize the structural requisites for the observed mass fragmentation, a series of 30 halogen-substituted 3-phenylpropenoate-related structures were examined. All ester-containing alkene derivatives exhibited loss of the distinctive chlorine from the 2-position of the phenyl ring. Analogous derivatives with the halogen (chlorine or bromine) in the para position did not evidence selective halogen loss. Results demonstrated that substituted 3-phenylpropenoates and their analogs fragment via the formation of a previously reported benzopyrylium intermediate. To understand the correlation between the intramolecular radical substitution and the abundance and selectivity of the chlorine (or other halogen) displacement, density functional theory calculations were performed to determine the charge on the principal cation involved in the chlorine loss (in the ortho, meta, and para positions), the charge for the neutral radical (noncation), the excess alpha-electron density on the relevant atom and the energy to form the cation from the neutral atom (ionization energy). Results showed that the selectivity and extent of halogen displacement correlated highly to the electrophilicity of the radical cation as well as the neutral radical. These data further support the proposed fragmentation mechanism involving intramolecular radical elimination.     

A characterization study on 2,6-dimethyl-4-nitropyridine N-oxide by density functional theory calculations

This study deals with the identification of a title compound, 2,6-dimethyl-4-nitropyridine N-oxide by means of theoretical calculations. The optimized molecular structures, vibrational frequencies, corresponding vibrational assignments, thermodynamic properties and atomic charges of the title compound in the ground state were evaluated using density functional theory (DFT) with the standard B3LYP/6-311G(d,p) method and basis set combination for the first time. Theoretical vibrational spectra were interpreted with the aid of normal coordinate analysis based on scaled density functional force field. The results show that the optimized geometric parameters (bond lengths and bond angles) and vibrational frequencies were observed to be in good agreement with the available experimental results. Based on the results of comparison between experimental results and theoretical data, the chosen calculation level is powerful approach for understanding the molecular structures and vibrational spectra of the 2,6-dimethyl-4-nitropyridine N-oxide. Moreover, we not only simulated frontier molecular orbitals (FMO) and molecular electrostatic potential (MEP) but also determined the transition state and energy band gap. Based on the investigations, the title compound is found to be useful to bond metallically and interact intermolecularly. Infrared intensities and Raman activities were also reported.     

Structure of the ammonia dimer studied by density functional theory

Self-consistent Kohn–Sham density functional calculations have been carried out to study the structure of the ammonia dimer. The local-density approximation yields unusually large binding energy and short internitrogen distance compared with the experimental and more accurate theoretical data. The results from the Becke–Perdew gradient-corrected functionals are generally in good agreement with those at the SCF MP 2 level when the geometry is fully optimized with various large basis sets. With our best estimation, the staggered quasi-linear structure (C_s) is 0.6 kcal/mol lower in energy than the symmetric cyclic one (C_2h). The hydrogen-bonded N—H bond in the staggered quasi-linear structure is found to be 0.008 Å longer than the N—H bond in ammonia. In our calculations, we could not find the minima on the energy surface corresponding to the two asymmetric cyclic structures suggested by microwave spectra and coupled pair functional calculations. © 1994 John Wiley & Sons, Inc.     

Molecular structure, vibrational spectral studies of pyrazole and 3,5-dimethyl pyrazole based on density functional calculations

In this work, the experimental and theoretical vibrational spectra of pyrazole (PZ) and 3,5-dimethyl pyrazole (DMP) have been studied. FTIR and FT-Raman spectra of the title compounds in the solid phase are recorded in the region 4000-400 cm(-1) and 4000-50 cm(-1), respectively. The structural and spectroscopic data of the molecules in the ground state are calculated using density functional methods (B3LYP) with 6-311+G** basis set. The vibrational frequencies are calculated and scaled values are compared with experimental FTIR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete vibrational assignments are performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanical (SM) method. 13C and 1H NMR chemical shifts results are compared with the experimental values.