Asymptotic behavior and interpretation of virtual states: The effects of confinement and of basis sets

A real-space high order finite difference method is used to analyze the effect of spherical domain size on the Hartree-Fock (and density functional theory) virtual eigenstates. We show the domain size dependence of both positive and negative virtual eigenvalues of the Hartree-Fock equations for small molecules. We demonstrate that positive states behave like a particle in spherical well and show how they approach zero. For the negative eigenstates, we show that large domains are needed to get the correct eigenvalues. We compare our results to those of Gaussian basis sets and draw some conclusions for real-space, basis-sets, and plane-waves calculations. (C) 2016 AIP Publishing LLC.

Intramolecular hydrogen bond: Can it be part of the basis set of valence internal coordinates in normal mode analysis?

It has been shown earlier1] that the relaxed force constants (RFCs) could be used as a measure of bond strength only when the bonds form a part of the complete valence internal coordinates (VIC) basis. However, if the bond is not a part of the complete VIC basis, its RFC is not necessarily a measure of bond strength. Sometimes, it is possible to have a complete VIC basis that does not contain the intramolecular hydrogen bond (IMHB) as part of the basis. This means the RFC of IMHB is not necessarily a measure of bond strength. However, we know that IMHB is a weak bond and hence its RFC has to be a measure of bond strength. We resolve this problem of IMHB not being part of the complete basis by postulating `equivalent' basis sets where IMHB is part of the basis at least in one of the equivalent sets of VIC. As long as a given IMHB appears in one of the equivalent complete VIC basis sets, its RFC could be used as a measure of bond strength parameter.

A comparative study of the energetics, structures, and mechanisms of the HCN H HNC and LiCN <-> LiNC isomerizations

The potential energy surfaces of the HCN<->HNC and LiCN<->LiNC isomerization processes were determined by ab initio theory using fully optimized triple-zeta double polarization types of basis sets. Both the MP2 corrections and the QCISD level of calculations were performed to correct for the electron correlation. Results show that electron correlation has a considerable influence on the energetics and structures. Analysis of the intramolecular bond rearrangement processes reveals that, in both cases, H (or Li+) migrates in an almost elliptic path in the plane of the molecule. In HCN<->HNC, the migrating hydrogen interacts with the in-plane pi,pi* orbitals of CN, leading to a decrease in the C-N bond order. In LiCN<->LiNC, Li+ does not interact with the corresponding pi,pi* orbitals of CN.

Ab initio molecular orbital calculations on ion-molecule and ion pair-molecule complexes of the water-lithium cyanide system

Ab initio molecular orbital (MO) calculations with the 3-21G and 6-31G basis sets were performed on a series of ion-molecule and ion pair-molecule complexes for the H2O + LiCN system. Stabilisation energies (with counter-poise corrections), geometrical parameters, internal force constants and harmonic vibrational frequencies were evaluated for 16 structures of interest. Although the interaction energies are smaller, the geometries and relative stabilities of the monohydrated contact ion pair are reminiscent of those computed for the complexes of the individual ions. Thus, interaction of the oxygen lone pair with lithium leads to a highly stabilised C2v structure, while the coordination of water to the cyanide ion involves a slightly non-linear hydrogen bond. Symmetrical bifurcated structures are computed to be saddle points on the potential energy surface, and to have an imaginary frequency for the rocking mode of the water molecule. On optimisation the geometries of the solvent shared ion pair structures (e.g. Li+cdots, three dots, centered OH2cdots, three dots, centered CN−) revealed a proton transfer from the water molecule leading to hydrogen bonded forms such as Li-O-Hcdots, three dots, centered HCN. The variation in the force constants and harmonic frequencies in the various structures considered are discussed in terms of ion-molecular and ion pair-molecule interactions.

AB-initio study of alpha-silyl and beta-silyl substituent effects on vinyl radicals

The geometries of alpha- and beta-silyl substituted vinyl radicals and of alpha,beta-disilylvinyl radical have been optimised with the STO-3G and the STO-3G* basis sets. The relative stabilities of various conformers have been determined at the UMP2/6-31G* level. The stabilisation of vinyl radicals through alpha-silyl substitution is larger than that due to corresponding alkyl groups. The presence of an alpha-silyl group also leads to a tendency towards linearisation of the vinyl radical centre and a corresponding reduction in the inversion barrier. In marked contrast, the beta-silyl effect is negligible. The geometric, conformational and energetic consequences are insignificant. Overall, the silyl substituent effect at vinyl radicals is very different from that computed earlier for the vinyl cations, but qualitatively similar to that found in carbanions.

SCF and electron correlation studies on structure, force constants and vibrational spectra of borane monoammoniate complex

Effects of basis set and electron correlation on the equilibrium geometry, force constants and vibrational spectra of BH3NH3 have been studied. A series of basis sets ranging from double zeta to triple zeta including polarization and diffuse functions have been utilized. All the SCF based calculations overestimate the dative B-N bond distance and considerable improvement occurs when the treatment for electron correlation is introduced. Detailed vibrational analysis for BH3NH3 has been carried out. The mean absolute percentage deviation of the ab initio predicted vibration frequencies of (BH3NH3)-B-11 from the experiment is about 10% for the SCF based calculations and the MP2 method shows better agreement, the overall deviation being 5-6%. The ground state effective force constants of BH3NH3 were obtained using RECOVES procedure. The RECOVES sets of force constants are found to be highly satisfactory for the prediction of the vibrational frequencies of different isotopomers of BH3NH3. The mean absolute percentage deviation of the calculated frequencies of different isotopomers from the experiment is much less than 1%. The RECOVES-MP2/augDZP set of force constants was found to be the best set among the different sets for this molecule. Theoretical infrared intensities are in fair agreement with the observed spectral features.

[n]peristylanes and [n]oxa[n]peristylanes (n = 3-6): A theoretical study

Theoretical studies at the HF and Becke3LYP levels using 6-31G* basis sets were carried out on a series of [n]peristylanes and [n]oxa[n]peristylanes (n = 3-6) to understand their structure and energetics. The structures of the [3]- and [4]peristylanes (1, 2) and their era-derivatives (5, 6) were calculated to have the anticipated high symmetry, C-nv. In contrast, a C-s structure (9) at HF/6-31G* and another (25) at the Becke3LYP/6-31G* level were calculated for the [5]oxa[5]peristylane. The energy difference between them is extremely small even though there are major differences in the structures indicating every soft potential energy surface: On the other hand, the potential energy surface of [6]oxa[6]peristylane is not as soft. Similar structures were also calculated for the top rings. Calculations on the seco-compounds 11-14 and 15-19 (Table 4) indicate that there is no unusual strain involved in the formation of 27 from 19. The Li+ interaction energies of the [n]oxa[n]peristylanes are 61.7 (n = 3), 72.8 (n = 4), 84.2 (n = 5) and 91.7 (n = 6) kcal mol(-1) at the Becke3LYP/6-3IG* level. Dramatic differences between the C-C bond lengths obtained from the solid state X-ray diffraction studies and those from the calculations for the [n]oxa[n]peristylanes were also observed.

Ab initio studies of solvent effect on conformational equilibria and vibrational spectra of dipropionamide

Systematic ab initio molecular orbital studies of the conformational equilibria and vibrational spectra of dipropionamide using the basis sets 6-31g(d) and 6-31++G(d,p) have been carried out. The vibrational spectra of dipropionamide have been satisfactorily interpreted taking into account the agreement between the calculated frequencies, infrared and Raman band intensities and the shifts in the spectra of deuterated molecules with those observed. The previous assignments of most of the vibrational bands are well confirmed, a few bands need reassignment, however. The solvent effects were investigated by self-consistent reaction field theory using dipole and self-consistent isodensity polarized continuum model methods. The introduction of a dielectric medium has only a marginal effect on the conformational equilibria and vibrational spectra. However, the calculated changes in geometry and vibrational spectra on going from the gas phase to the solution phase are in accord with the increasing weight of the dipolar resonance structure in polar solvents. (C) 2002 Elsevier Science B.V. All rights reserved.

Ab initio, DFT and transition state theory calculations on 1,2-HF, HCI and CIF elimination reactions from CH2F-CH2Cl

This paper reports ab intio, DFT and transition state theory (TST) calculations on HF, HCI and CIF elimination reactions from CH2Cl-CH2F molecule. Both the ground state and the transition state for HX elimination reactions have been optimized at HF, MP2 and DFT calculations with 6-31G*, 6-31G** and 6-311++G** basis sets. In addition, CCSD(T) single point calculations were carried out with MP2/6-311++G** optimized geometry for more accurate determination of the energies of the minima and transition state, compared to the other methods employed here. Classical barriers are converted to Arrhenius activation energy by TST calculations for comparisons with experimental results. The pre-exponential factors, A, calculated at all levels of theory are significantly larger than the experimental values. For activation energy, E-a DFT gives good results for HF elimination, within 4-8 W mol(-1) from experimental values. None of the methods employed, including CCSD(T), give comparable results for HCI elimination reactions. However, rate constants calculated by CCSD(T) method are in very good agreement with experiment for HCI elimination and they are in reasonable agreement for HF elimination reactions. Due to the strong correlation between A and E., the rate constants could be fit to a lower A and E-a (as given by experimental fitting, corresponding to a tight TS) or to larger A and E-a (as given by high level ab initio calculations, corresponding to a loose TS). The barrier for CIF elimination is determined to be 607 U mol(-1) at HF level and it is unlikely to be important for CH2FCH2Cl. Results for other CH2X-CH2Y (X,Y = F/Cl) are included for comparison.

Thermal decomposition of haloethanols: single pulse shock tube and ab initio studies

This paper reports single pulse shock tube and ab initio studies on thermal decomposition of 2-fluoro and 2-chloroethanol at T=1000–1200 K. Both molecules have HX (X = F/Cl) and H2O molecular elimination channels. The CH3CHO formed by HX elimination is chemically active and undergoes secondary decomposition resulting in the formation of CH4, C2H6, and C2H4. A detailed kinetic simulation indicates that the formation of C2H4 could not be quantitatively explained as arising exclusively from secondary CH3CHO decomposition. Contributions from primary radical processes need to be considered to explain C2H4 quantitatively. Ab initio calculations on HX and H2O elimination reactions from the haloethanols at HF, MP2, and DFT levels with various basis sets up to 6/311++G**are reported. It is pointed out that due to strong correlations between A and Eα, comparison of these two parameters between experimental and theoretical results could be misleading.

Study on the structure and vibrational spectra of efavirenz conformers using DFT: Comparison to experimental data

Efavirenz, (S)-6-chloro-4-(cyclopropylethynyl)-1,4-dihydro-4-(trifluoromethyl)-2H-3 ,1-benzoxazin-2-one, is an anti HIV agent belonging to the class of the non-nucleoside inhibitors of the HIV-1 virus reverse transcriptase. A systematic quantum chemical study of the possible conformations, their relative stabilities and vibrational spectra of efavirenz has been reported. Structural and spectral characteristics of efavirenz have been studied by vibrational spectroscopy and quantum chemical methods. Density functional theory (DFT) calculations for potential energy curve, optimized geometries and vibrational spectra have been carried out using 6-311++G(d,p) basis sets and B3LYP functionals. Based on these results, we have discussed the correlation between the vibrational modes and the crystalline structure of the most stable form of efavirenz. A complete analysis of the experimental infrared and Raman spectra has been reported on the basis of wavenumber of the vibrational bands and potential energy distribution. The infrared and the Raman spectra of the molecule based on OFT calculations show reasonable agreement with the experimental results. The calculated HOMO and LUMO energies shows that charge transfer occur within the molecule. (C) 2011 Elsevier B.V. All rights reserved.

Cu-II-azide polynuclear complexes of Cu-4 building clusters with Schiff-base co-ligands: synthesis, structures, magnetic behavior and DFT studies

Three new copper-azido complexes Cu-4(N-3)(8)(L-1)(2)](n) (1), Cu-4(N-3)(6)(L-2)(2)(H2O)(2)] (2), and Cu-4(N-3)(6)(L-3)(2)](n) (3) L-1 is the imine resulting from the condensation of pyridine-2-carboxaldehyde with N-methylethylenediamine, HL2 and HL3 are the condensation products of 2-hydroxy-3-methoxybenzaldehyde with N,N-diethylethylenediamine and N-ethylethylenediamine respectively] have been synthesized by using 0.5 molar equivalents of the Schiff base ligands with Cu(NO3)(2)center dot 3H(2)O and an excess of NaN3. Single crystal X-ray structures show that the basic unit of these complexes contains very similar Cu-4(II) building blocks. While 1 and 3 have overall 1D structures, 2 forms discrete tetranuclear clusters due to blocking of two coordination sites on the tetranuclear cluster by water molecules. Magnetic susceptibility measurements over a wide range of temperatures exhibit the presence of both antiferromagnetic and ferromagnetic exchanges within the tetranuclear unit structures. Density functional theory calculations (using B3LYP functional and two different basis sets) have been performed on the complexes 1-3 to provide a qualitative theoretical interpretation of their overall magnetic behavior.

Cu-II-Azide Polynuclear Complexes of Three Different Building Clusters with the Same Schiff-Base Ligand: Synthesis, Structures, Magnetic Behavior, and Density Functional Theory Studies

Three copper-azido complexes Cu-4(N-3)(8)(L-1)(2)(MeOH)(2)](n) (1), Cu-4(N-3)(8)(L-1)(2)] (2), and Cu-5(N-3)(10)(L-1)(2)](n) (3) L-1 is the imine resulting from the condensation of pyridine-2-carboxaldehyde with 2-(2-pyridyl)ethylamine] have been synthesized using lower molar equivalents of the Schiff base ligand with Cu(NO3)(2)center dot 3H(2)O and an excess of NaN3. Single crystal X-ray structures show that the basic unit of the complexes 1 and 2 contains Cu-4(II) building blocks; however, they have distinct basic and overall structures due to a small change in the bridging mode of the peripheral pair of copper atoms in the linear tetranudear structures. Interestingly, these changes are the result of changing the solvent system (MeOH/H2O to EtOH/H2O) used for the synthesis, without changing the proportions of the components (metal to ligand ratio 2:1). Using even lower proportions of the ligand, another unique complex was isolated with Cu-5(II) building units, forming a two-dimensional complex (3). Magnetic susceptibility measurements over a wide range of temperature exhibit the presence of both antiferromagnetic (very weak) and ferromagnetic exchanges within the tetranuclear unit structures. Density functional theory calculations (using B3LYP functional, and two different basis sets) have been performed on the complexes 1 and 2 to provide a qualitative theoretical interpretation of their overall magnetic behavior.

Conformational Stability and Intramolecular Hydrogen Bonding in 1,2-Ethanediol and 1,4-Butanediol

The gas-phase infrared spectra of 1,2-ED and 1,4-BD have been, recorded at three different temperatures using a multipass gas cell of 6 m optical path length. DFT calculation has also been carried out using 6-311++G** and aug-cc-pVDZ basis sets to look for the existence of intramolecular hydrogen bonding, in them from the red shift and infrared absorption intensity enhancement of the bonded O-H band compared to that of the free O-H hand. Equilibrium population analysis With 10 conformers of 1,2-ED and 1,4-BD at experimental temperatures were-carried out for the reconstruction of the Observed vibrational spectra at that temperature,using standard statistical relationships. The most abundant conformer at experimental temperatures, was identified. In 1,2-ED a red shift of 45 cm(-1) in the intramolecularly interacting O-H stretching vibrational band position and no significant intensity enhancement compared to that of-the free O-H have been observed. On the contrary, in one of the hydrogen bonded conformers of 1,4-BD, a 124 cm(-1) red shift in the O-H stretching frequency and 8.5 times-intensity enhancement for the ``bonded'' O-H compared to that of the ``free'' O-H is seen. On the basis of this comparative study, we have concluded that strong intramolecular hydrogen bonding exists in 1,4-BD. But there appears, to be weak intramolecular hydrogen bonding in 1,2-ED at temperatures of 303, 313, and 323 Km the gas phase We have found that most stable hydrogen-bonded conformers of 1,4-BD are less populated than some of the non-hydrogen-bonded conformers. Even for the 1,4-BD, the relative population of the g'GG'Gt conformer, which has a strong intramolecular,hydrogen bond, is less than what is predicted. Perhaps the intramolecular hydrogen bond plays a less Significant role in the relative stability,of the various Conformers than what has been predicted from calculations and prevails in the literature.

Prospects of zero Schottky barrier height in a graphene-inserted MoS2-metal interface

A low Schottky barrier height (SBH) at source/drain contact is essential for achieving high drive current in atomic layer MoS(2-)channel-based field effect transistors. Approaches such as choosing metals with appropriate work functions and chemical doping are employed previously to improve the carrier injection from the contact electrodes to the channel and to mitigate the SBH between the MoS2 and metal. Recent experiments demonstrate significant SBH reduction when graphene layer is inserted between metal slab (Ti and Ni) and MoS2. However, the physical or chemical origin of this phenomenon is not yet clearly understood. In this work, density functional theory simulations are performed, employing pseudopotentials with very high basis sets to get insights of the charge transfer between metal and monolayer MoS2 through the inserted graphene layer. Our atomistic simulations on 16 different interfaces involving five different metals (Ti, Ag, Ru, Au, and Pt) reveal that (i) such a decrease in SBH is not consistent among various metals, rather an increase in SBH is observed in case of Au and Pt; (ii) unlike MoS2-metal interface, the projected dispersion of MoS2 remains preserved in any MoS2-graphene- metal system with shift in the bands on the energy axis. (iii) A proper choice of metal (e.g., Ru) may exhibit ohmic nature in a graphene-inserted MoS2-metal contact. These understandings would provide a direction in developing high-performance transistors involving heteroatomic layers as contact electrodes. (c) 2016 AIP Publishing LLC.