Anharmonicity contributions to the vibrational second hyperpolarizability of conjugated oligomers

Restricted Hartree-Fock 6-31G calculations of electrical and mechanical anharmonicity contributions to the longitudinal vibrational second hyperpolarizability have been carried out for eight homologous series of conjugated oligomers - polyacetylene, polyyne, polydiacetylene, polybutatriene, polycumulene, polysilane, polymethineimine, and polypyrrole. To draw conclusions about the limiting infinite polymer behavior, chains containing up to 12 heavy atoms along the conjugated backbone were considered. In general, the vibrational hyperpolarizabilities are substantial in comparison with their static electronic counterparts for the dc-Kerr and degenerate four-wave mixing processes (as well as for static fields) but not for electric field-induced second harmonic generation or third harmonic generation. Anharmonicity terms due to nuclear relaxation are important for the dc-Kerr effect (and for the static hyperpolarizability) in the σ-conjugated polymer, polysilane, as well as the nonplanar π systems polymethineimine and polypyrrole. Restricting polypyrrole to be planar, as it is in the crystal phase, causes these anharmonic terms to become negligible. When the same restriction is applied to polymethineimine the effect is reduced but remains quantitatively significant due to the first-order contribution. We conclude that anharmonicity associated with nuclear relaxation can be ignored, for semiquantitative purposes, in planar π-conjugated polymers. The role of zero-point vibrational averaging remains to be evaluated

The mapping of the local contributions of Fermi and Coulomb correlation into intracule and extracule density distributions

The contributions of the correlated and uncorrelated components of the electron-pair density to atomic and molecular intracule I(r) and extracule E(R) densities and its Laplacian functions ∇2I(r) and ∇2E(R) are analyzed at the Hartree-Fock (HF) and configuration interaction (CI) levels of theory. The topologies of the uncorrelated components of these functions can be rationalized in terms of the corresponding one-electron densities. In contrast, by analyzing the correlated components of I(r) and E(R), namely, IC(r) and EC(R), the effect of electron Fermi and Coulomb correlation can be assessed at the HF and CI levels of theory. Moreover, the contribution of Coulomb correlation can be isolated by means of difference maps between IC(r) and EC(R) distributions calculated at the two levels of theory. As application examples, the He, Ne, and Ar atomic series, the C2-2, N2, O2+2 molecular series, and the C2H4 molecule have been investigated. For these atoms and molecules, it is found that Fermi correlation accounts for the main characteristics of IC(r) and EC(R), with Coulomb correlation increasing slightly the locality of these functions at the CI level of theory. Furthermore, IC(r), EC(R), and the associated Laplacian functions, reveal the short-ranged nature and high isotropy of Fermi and Coulomb correlation in atoms and molecules

The relevance of the Laplacian of intracule and extracule density distributions for analyzing electron-electron interactions in molecules

A topological analysis of intracule and extracule densities and their Laplacians computed within the Hartree-Fock approximation is presented. The analysis of the density distributions reveals that among all possible electron-electron interactions in atoms and between atoms in molecules only very few are located rigorously as local maxima. In contrast, they are clearly identified as local minima in the topology of Laplacian maps. The conceptually different interpretation of intracule and extracule maps is also discussed in detail. An application example to the C2H2, C2H4, and C2H6 series of molecules is presented

Nuclear relaxation contribution to static and dynamic (infinite frequency approximation) nonlinear optical properties by means of electrical property expansions: Application to HF, CH4, CF4, and SF6

Electrical property derivative expressions are presented for the nuclear relaxation contribution to static and dynamic (infinite frequency approximation) nonlinear optical properties. For CF4 and SF6, as opposed to HF and CH4, a term that is quadratic in the vibrational anharmonicity (and not previously evaluated for any molecule) makes an important contribution to the static second vibrational hyperpolarizability of CF4 and SF6. A comparison between calculated and experimental values for the difference between the (anisotropic) Kerr effect and electric field induced second-harmonic generation shows that, at the Hartree-Fock level, the nuclear relaxation/infinite frequency approximation gives the correct trend (in the series CH4, CF4, SF6) but is of the order of 50% too small

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

One- and two-center physical space partitioning of the energy in the density functional theory

A conceptually new approach is introduced for the decomposition of the molecular energy calculated at the density functional theory level of theory into sum of one- and two-atomic energy components, and is realized in the "fuzzy atoms" framework. (Fuzzy atoms mean that the three-dimensional physical space is divided into atomic regions having no sharp boundaries but exhibiting a continuous transition from one to another.) The new scheme uses the new concept of "bond order density" to calculate the diatomic exchange energy components and gives them unexpectedly close to the values calculated by the exact (Hartree-Fock) exchange for the same Kohn-Sham orbitals

Índexos de localització i deslocalització derivats de la densitat bielectrònica: anàlisi i aplicacions en estructura molecular, reactivitat química i aromaticitat

La present tesi doctoral s'emmarca dins el camp de la química teòrica i computacional. El primer objectiu, que va ser el de partida, consistia a veure si els índexos bielectrònics derivats de la teoria AIM podien ser útils per estudiar la reorganització electrònica al llarg d'una reacció. De forma paral·lela es va plantejar un objectiu metodològic, el càlcul dels índexos mencionats a altres nivells de teoria, com podia ser DFT o CI. CI. Aquests dos objectius varen portar a molts d'altres, tot i que els principals són: Objectiu 1. Estudi de la reorganització de parells d'electrons al llarg d'una reacció. Objectiu 2. Efecte de la solvatació sobre l'estructura de parells d'electrons d'una molècula. Objectiu 3. Càlcul dels índexos de localització i deslocalització als nivells de teoria del funcional de la densitat i post-Hartree-Fock. Objectiu 4. Caracterització de l'enllaç per pont d'hidrogen per mitjà de l'estructura dels parells d'electrons. Objectiu 5. L'índex de deslocalització: nou criteri electrònic d'aromaticitat.

Spin polarization, orbital occupation and band gap opening in vanadium dioxide: the effect of screened Hartree–Fock exchange

The metal–insulator transition of VO2 so far has evaded an accurate description by density functional theory. The screened hybrid functional of Heyd, Scuseria and Ernzerhof leads to reasonable solutions for both the low-temperature monoclinic and high-temperature rutile phases only if spin polarization is excluded from the calculations. We explore whether a satisfactory agreement with experiment can be achieved by tuning the fraction of Hartree Fock exchange (a) in the density functional. It is found that two branches of locally stable solutions exist for the rutile phase for 12:5% 6 a 6 20%. One is metallic and has the correct stability as compared to the monoclinic phase, the other is insulating with lower energy than the metallic branch. We discuss these observations based on the V 3d orbital occupations and conclude that a ¼ 10% is the best possible choice for spin-polarized VO2 calculations.

Role of the Skyrme tensor force in heavy-ion fusion

We make use of the Skyrme effective nuclear interaction within the time-dependent Hartree-Fock framework to assess the effect of inclusion of the tensor terms of the Skyrme interaction on the fusion window of the 16O–16O reaction. We find that the lower fusion threshold, around the barrier, is quite insensitive to these details of the force, but the higher threshold, above which the nuclei pass through each other, changes by several MeV between different tensor parametrisations. The results suggest that eventually fusion properties may become part of the evaluation or fitting process for effective nuclear interactions.

New insights into the structure of exotic nuclei using the RISING active stopper

This conference paper outlines the operation and some of the preliminary physics results using the GSI RISING active stopper. Data are presented from an experiment using combined isomer and beta‐delayed gamma‐ray spectroscopy to study low‐lying spectral and decay properties of heavy‐neutron‐rich nuclei around A∼190 produced following the relativistic projectile fragmentation of 208Pb primary beam. The response of the RISING active stopper detector is demonstrated for both the implantation of heavy secondary fragments and in‐situ decay of beta‐particles. Beta‐delayed gamma‐ray spectroscopy following decays of the neutron‐rich nucleus 194Re is presented to demonstrate the experimental performance of the set‐up. The resulting information inferred from excited states in the W and Os daughter nuclei is compared with results from Skyrme Hartree‐Fock predictions of the evolution of nuclear shape.

Probing the (6,7)Li nucleon densities through a new break-up process approach

In this work we present a double folding optical model analysis of new near-barrier quasi-elastic experimental data for the (6,7)Li + (120)Sn systems. From the analysis, it was possible to confirm the ground-state nucleon densities assumed for the weakly bound (6,7)Li isotopes. The apparent discrepancies between the experimental densities and those based on Dirac-Hartree-Fock Bogoliubov (DHB) calculations were removed. A new approach that simulates the projectile break-up and a positive polarization from couplings of (6,7)Li bound states with the continuum was considered in the reaction mechanism. (C) 2010 Elsevier B.V. All rights reserved.

Extended analysis of the Kr V spectrum

The spectrum of four-times-ionized krypton (Kr V) has been observed in the 230-4900 angstrom wavelength range, resulting in 91 new classified lines. We were able to identify 21 new energy levels belonging to the 4s(2)4p5d, 4s(2)4p5s, 4s(2)4p6s, 4s(2)4p5p and 4s4p(2)4d configurations. Relativistic Hartree-Fock calculations were used to predict energy levels and transitions. (c) 2010 Elsevier Ltd. All rights reserved.

Finite temperature correction to the Thomas-Fermi approximation for a Bose-Einstein condensate: comparison between theory and experiment

We observe experimentally a deviation of the radius of a Bose-Einstein condensate from the standard Thomas-Fermi prediction, after free expansion, as a function of temperature. A modified Hartree-Fock model is used to explain the observations, mainly based on the influence of the thermal cloud on the condensate cloud.

Exploring new molecular species on the (1)[H, Se, F] singlet potential energy surface: Energetics, structures, IR spectra, and heats of formation

Structural, energetic, and vibrational properties of new molecular species, HSeF and HFSe, the associated transition state, and dissociation fragments are investigated using a state-of-the-art theoretical approach, CCSD(T)/CBS. HSeF is a normal covalently bonded molecule 38.98 kcal mol (1) more stable than the complex HF-Se, which shows an unusual structure with a central fluorine atom and a bond angle of 101.8 degrees.A barrier (Delta G(#)) of 49.01 kcal mol (1) separates the two species. Vibrational frequencies are also quite distinct. Heats of formation are evaluated for the diatomic fragments and HSeF. Final Delta(f)H values depend on the experimental accuracy of those of Se(g) and H(2)Se. (c) 2009 Elsevier B.V. All rights reserved.

Removal of metal ions from aqueous solution by chelating polymeric hydrogel

Polysaccharide natural seed coat from the tree Magonia pubescens, in the form of hydrogel was used to remove metals in aqueous solution. Swelling tests indicate that seed coat presents hydrogel behavior, with maximum water absorption of 292 g water/g. Adsorption experiments performed using Na(+), Mg(2+), K(+), Ca(2+), Cr(3+), Fe(3+) and Zn(2+) demonstrated that the polysaccharide structure has a high capacity to extract these ions from the aqueous solution. Scanning electron microscopy revealed significant morphological changes of the material before and after water contact. Differential scanning calorimetry measurements indicate a signal shift of the water evaporation temperature in the material with adsorbed zinc. X-ray photoelectron spectroscopy analysis combined with theoretical studies by the density functional theory and on Hartree-Fock (HF) level evidence that the metallic ions were adsorbed through coordination with hydroxyl groups of polysaccharide. In the case of Zn(2+) the lowest HF energy was observed for the tetracoordination mode, where Zn(2+) is coordinated by two hydroxyl groups and two water molecules.