Theoretical study of one- and two-photon absorption spectra of azoaromatic compounds

In this study, the one- and two-photon absorption spectra of seven azoaromatic compounds (five pseudostilbenes-type and two aminoazobenzenes) were theoretically investigated using the density functional theory combined with the response functions formalism. The equilibrium molecular structure of each compound was obtained at three different levels of theory: Hartree-Fock, density functional theory (DFT), and Moller-Plesset 2. The effect of solvent on the equilibrium structure and the electronic transitions of the compounds were investigated using the polarizable continuum model. For the one-photon absorption, the allowed pi ->pi(*) transition energy showed to be dependent on the molecular structures and the effect of solvent, while the n ->pi(*) and pi ->pi(*)(n) transition energies exhibited only a slight dependence. An inversion between the bands corresponding to the pi ->pi(*) and n ->pi(*) states due to the effect of solvent was observed for the pseudostilbene-type compounds. To characterize the allowed two-photon absorption transitions for azoaromatic compounds, the response functions formalism combined with DFT using the hybrid B3LYP and PBE0 functionals and the long-range corrected CAM-B3LYP functional was employed. The theoretical results support the previous findings based on the three-state model. The model takes into account the ground and two electronic excited states and has already been used to describe and interpret the two-photon absorption spectrum of azoaromatic compounds. The highest energy two-photon allowed transition for the pseudostilbene-type compounds shows to be more effectively affected (similar to 20%) by the torsion of the molecular structure than the lowest allowed transition (similar to 10%). In order to elucidate the effect of the solvent on the two-photon absorption spectra, the lowest allowed two-photon transition (dipolar transition) for each compound was analyzed using a two-state approximation and the polarizable continuum model. The results obtained reveal that the effect of solvent increases drastically the two-photon cross-section of the dipolar transition of the pseudostilbene-type compounds. In general, the features of both one- and two-photon absorption spectra of the azoaromatic compounds are well reproduced by the theoretical calculations.

Polarized–unpolarized ground state of small polycyclic aromatic hydrocarbons

Do polyacenes, circumacenes, periacenes, nanographenes, and graphene nanoribbons show a spin polarized ground state? In this work, we present monodeterminantal (Hartree–Fock (HF) and density functional theory (DFT) types), and multideterminantal calculations (Møller–Plesset and Coupled Cluster), for several families of unsaturated organic molecules (n-Acenes, n-Periacenes and n-Circumacenes). All HF calculations and many DFT show a spin-polarized (antiferromagnetic) ground state, in agreement with previous calculations. Nevertheless, the multideterminantal calculations, carried out with perturbative and variational wavefunctions, show that the more stable state is obtained starting from the unpolarized HF wavefunction. The trend of the stabilization of wavefunctions (polarized or unpolarized) with respect to exchange and correlation potentials, and to the number of benzene rings, has been analyzed. A study of the spin (〈Ŝ2〉) and the spin density on the carbon atoms has also been carried out.

Electric field gradients in (111)In-doped (Hf/Zr)(3)Al(2) and (Hf/Zr)(4)Al(3) mixed compounds: ab initio calculations, perturbed angular correlation measurements and site preference

The quadrupolar hyperfine interactions of in-diffused (111)In -> (111)Cd probes in polycrystalline isostructural Zr(4)Al(3) and Hf(4)Al(3) samples containing small admixtures of the phases (Zr/Hf)(3)Al(2) were investigated. A strong preference of (111)In solutes for the contaminant (Zr/Hf)(3)Al(2) minority phases was observed. Detailed calculations of the electric field gradient (EFG) at the Cd nucleus using the full-potential augmented plane wave + local orbital formalism allowed us to assign the observed EFG fractions to the various lattice sites in the (Zr/Hf)(3)Al(2) compounds and to understand the preferential site occupation of the minority phases by the (111)In atoms. The effects of the size of the supercell and relaxation around the oversized In and Cd probe atoms were investigated in detail.

The role of electron localization in the atomic structure of transition-metal 13-atom clusters: the example of Co(13), Rh(13), and Hf(13)

The crystalline structure of transition-metals (TM) has been widely known for several decades, however, our knowledge on the atomic structure of TM clusters is still far from satisfactory, which compromises an atomistic understanding of the reactivity of TM clusters. For example, almost all density functional theory (DFT) calculations for TM clusters have been based on local (local density approximation-LDA) and semilocal (generalized gradient approximation-GGA) exchange-correlation functionals, however, it is well known that plain DFT fails to correct the self-interaction error, which affects the properties of several systems. To improve our basic understanding of the atomic and electronic properties of TM clusters, we report a DFT study within two nonlocal functionals, namely, the hybrid HSE (Heyd, Scuseria, and Ernzerhof) and GGA + U functionals, of the structural and electronic properties of the Co(13), Rh(13), and Hf(13) clusters. For Co(13) and Rh(13), we found that improved exchange-correlation functionals decrease the stability of open structures such as the hexagonal bilayer (HBL) and double simple-cubic (DSC) compared with the compact icosahedron (ICO) structure, however, DFT-GGA, DFT-GGA + U, and DFT-HSE yield very similar results for Hf(13). Thus, our results suggest that the DSC structure obtained by several plain DFT calculations for Rh(13) can be improved by the use of improved functionals. Using the sd hybridization analysis, we found that a strong hybridization favors compact structures, and hence, a correct description of the sd hybridization is crucial for the relative energy stability. For example, the sd hybridization decreases for HBL and DSC and increases for ICO in the case of Co(13) and Rh(13), while for Hf(13), the sd hybridization decreases for all configurations, and hence, it does not affect the relative stability among open and compact configurations.

Basis set convergence in Hartree-Fock calculations of some diatomic molecules containing first and second-row atoms

Investiga-se a convergência de conjuntos de bases em direção ao limite numérico da energia Hartree-Fock (HF) total para as seqüências hierárquicas dos conjuntos de bases XZP e ccpVXZ. Para as duas hierarquias, melhoramentos significativos são obtidos com cada incremento em X. Para estimar o limite do conjunto de base completo, uma forma exponencial foi usada. Entre as várias aproximações consideradas aqui, uma extrapolação exponencial de três parâmetros aplicada aos resultados TZP, QZP e 5ZP deu os limites do conjunto de bases mais precisos. Em adição, energias HF dos orbitais moleculares ocupados mais altos de algumas moléculas diatômicas foram calculadas com o conjunto 5ZP e comparadas com as correspondentes obtidas com o conjunto cc-pV5Z e com um método numérico HF.

Counterpoise-corrected geometries and harmonic frequencies of N-body clusters: application to (HF)n (n=3,4)

A study was conducted on the methods of basis set superposition error (BSSE)-free geometry optimization and frequency calculations in clusters larger than a dimer. In particular, three different counterpoise schemes were critically examined. It was shown that the counterpoise-corrected supermolecule energy can be easily obtained in all the cases by using the many-body partitioning of energy

Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.

NMR investigation and theoretical calculations of the effect of solvent on the conformational analysis of 4',7-di-hydroxy-8-prenylflavan

The NMR conformational study of 4',7-di-hydroxy-8-prenylflavan 1 was carried out in acetone-d6, DMSO-d6 and CDCl3 which enabled the proposition of three conformations, namely 1a, 1b and 1c, differing in the position of the prenyl group. Geometry optimizations performed using AM1 method showed that 1a (deltaHf = -86.2 kcal/mol) is as stable as 1b (deltaHf = -85.1 kcal/mol) and 1c (deltaHf = -85.4 kcal/mol). When the solvent was included, the calculations showed that the solute-solvent interactions could be explained either in the light of the electronic intermolecular delocalization or the electrostatic character between solute and solvent. Theoretical calculations (HF/6-31G*, deltaFT/BLYP/6-31G*, and deltaFT/B3LYP/6-31G*) showed that the combination of these types of interactions present in each solute-solvent system, dependent on the chemical properties of the solvent, lead to different spatial arrangements of the prenyl group, which in turn determined the conformation of 1.

Estudo das ligações de hidrogênio para dímeros formados pelas moléculas de H2O, NH3, HF, HCl e HBr através de cálculos baseados em primeiros princípios

Hydrogen bond energies of fifteen dimers were calculated using the large basis set 6-311++G(3df,3pd), at Hartree-Fock (HF) level including Møller-Plesset (MP2) calculations. The procedure for obtaining such energies were based on the dimer's energy rise provoked by increasing in intermolecular distance of the system component units. Deviations from a strictly linear hydrogen bond were investigated and rotational barriers were also computed allowing the calculation of the second order attractive interactions. In order to provide a more objective definition of hydrogen bond, a lower energy limit was proposed in place of the merely empirical parameters employed in the classical definition

Covalent anion of Li(HF)3- system: an ab initio study

In this paper, we report the stability of the Li(HF)3- molecular anion calculated at the MP2/6-31++G** and CCSD(T)/6-31++G** level of theory. Five possible conformers of Li(HF)3- molecular anions have been determined employing ab initio MP2 method with 6-31++G** basis set. The most stable conformer of five Li(HF)3- anions is in a cyclic ring structure Li(HF)3-(1). From our calculations we show that the molecule is stable towards electron attachment, with an electron adiabatic electron affinity (AEA) of 199.5 meV (233.1 meV with zero point energy correction) and 471.3 meV at the MP2 and CCSD(T) levels, respectively. In addition we present vertical detachment energies of 230.2 meV and 795.8 meV at the MP2, CCSD(T), respectively. The importance of the latter has to do with the ability of experimental detection of this value.

Ionization potentials and radii of atoms and ions of element 104 (unnilquadium) and of hafnium (2+) derived from multiconfiguration Dirac-Fock calculations

Multiconfiguration relativistic Dirac-Fock (MCDF) values have been computed for the first four ionization potentials (IPs) of element 104 (unnilquadium) and of the other group 4 elements (Ti, Zr, and Hf). Factors were calculated that allowed correction of the systematic errors between the MCDF IPs and the experimental IPs. Single "experimental" IPs evaluated in eV (to ± 0.1 eV) for element 104 are: [104(0),6.5]; [104( 1 + ),14.8]; [104(2 + ),23.8]; [104(3 + ),31.9]. Multiple experimental IPs evaluated in eV for element 104 are: [(0-2+ ),21.2±0.2]; [(0-3+ ),45.1 ±0.2]; [(0-4+ ),76.8±0.3].Our MCDF results track 11 of the 12 experimental single IPs studied for group 4 atoms and ions. The exception is Hf( 2 + ). We submit our calculated IP of 22.4 ± 0.2 eV as much more accurate than the value of 23.3 eV derived from experiment.

Counterpoise-corrected geometries and harmonic frequencies of N-body clusters: application to (HF)n (n=3,4)

A study was conducted on the methods of basis set superposition error (BSSE)-free geometry optimization and frequency calculations in clusters larger than a dimer. In particular, three different counterpoise schemes were critically examined. It was shown that the counterpoise-corrected supermolecule energy can be easily obtained in all the cases by using the many-body partitioning of energy

Force constants and dipole moment derivatives of molecules from perturbed Hartree-Fock calculations II. Applications to limited basis set SCF-MO wavefunctions

The perturbed Hartree–Fock theory developed in the preceding paper is applied to LiH, BH, and HF, using limited basis‐set SCF–MO wavefunctions derived by previous workers. The calculated values for the force constant ke and the dipole‐moment derivative μ(1) are (experimental values in parentheses): LiH, ke  =  1.618(1.026)mdyn/Å,μ(1)  =  −18.77(−2.0±0.3)D/ÅBH,ke  =  5.199(3.032)mdyn/Å,μ(1)  =  −1.03(−)D/Å;HF,ke  =  12.90(9.651)mdyn/Å,μ(1)  =  −2.15(+1.50)D/Å. The values of the force on the proton were calculated exactly and according to the Hellmann–Feynman theorem in each case, and the discrepancies show that none of the wavefunctions used are close to the Hartree–Fock limit, so that the large errors in ke and μ(1) are not surprising. However no difficulties arose in the perturbed Hartree–Fock calculation, so that the application of the theory to more accurate wavefunctions appears quite feasible.

The molecular structure of hexamethyldigermane determined by gas-phase electron diffraction with theoretical calculations for (CH3)3M-M(CH3)3 Where M = C, Si, and Ge

Gas-phase electron diffraction (GED) data together with results from ab initio molecular orbital calculations (HF and MP2/6-311+G(d,p)) have been used to determine the structure of hexamethyldigermane ((CH3)3Ge-Ge(CH3)3). The equilibrium symmetry is D3d, but the molecule has a very low-frequency, largeamplitude, torsional mode (φCGeGeC) that lowers the thermal average symmetry. The effect of this largeamplitude mode on the interatomic distances was described by a dynamic model which consisted of a set of pseudoconformers spaced at even intervals. The amount of each pseudoconformer was obtained from the ab initio calculations (HF/6-311+G(d,p)). The results for the principal distances (ra) and angles (∠h1) obtained from the combined GED/ab initio (with estimated 1σ uncertainties) are r(Ge-Ge) ) 2.417(2) Å, r(Ge-C) ) 1.956(1) Å, r(C-H) ) 1.097(5) Å, ∠GeGeC ) 110.5(2)°, and ∠GeCH ) 108.8(6)°. Theoretical calculations were performed for the related molecules ((CH3)3Si-Si(CH3)3 and (CH3)3C-C(CH3)3).