Tertius L. Fonseca

Static polarizabilities of doubly charged polyacetylene oligomers: basis set and electron correlation effects

Abstract Ab initio calculations, carried out with different basis sets, for the static longitudinal linear polarizability, αL, and second order hyperpolarizability, γL, of small doubly charged polyacetylene (PA) chains, are presented. The polarizabilities were calculated using the Hartree–Fock (HF) method while the electron correlation effects were included through the second-order Moller–Plesset perturbation theory (MP2). Positively and negatively charged bipolarons were studied. The results obtained for positive and negative chains show that the ionization state effect decreases more rapidly, as the chain length is increased, for αL than for γL. For both types of charged chains, the incorporation of the electron correlation increases the αL and γL values, as compared to the HF values. A comparison between the results obtained using the standard 6-31G basis set and augmented versions of this set, obtained by the addition of diffuse and polarization functions, shows that 6-31G basis set does not provide a good description of the negative chains studied here and that the addition of extra diffuse functions on the basis set is needed in order to obtain reliable estimates for polarizabilities, specially for γL.

Hyperpolarizabilities of the methanol molecule: A CCSD calculation including vibrational corrections.

In this work we present the results for hyperpolarizabilities of the methanol molecule including vibrational corrections and electron correlation effects at the CCSD level. Comparisons to random phase approximation results previously reported show that the electron correlation is in general important for both electronic contribution and vibrational corrections. The role played by the anharmonicities on the calculations of the vibrational corrections has also been analyzed and the obtained results indicate that the anharmonic terms are important for the dc-Pockels and dc-Kerr effects. For the other nonlinear optical properties studied the double-harmonic approximation is found to be suitable. Comparison to available experimental result in gas phase for the dc-second harmonic generation second hyperpolarizability shows a very good agreement with the electronic contribution calculated here while our total value is 14% larger than the experimental value.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7+/-0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm.

Theoretical study of the static first hyperpolarizability of azo-enaminone compounds

In this work the static electric properties of azo-enaminones, with special emphasis to the vector component of the first hyperpolarizability βvec, are determined at the Hartree–Fock (HF) level with the electron correlation (EC) effects included through the second-order Moller–Plesset perturbation theory (MP2). The ab initio results, in accordance with previous semiempirical calculations, show that appropriate choices of substituents to be incorporated to the molecular structure can have a marked influence on the first hyperpolarizability. An initial study about the changes on the βvec values of these compounds, as a result of the incorporation of different donor groups, indicates that this property increases as function of the donor group strength tending to a saturated value. A comparison between our HF and MP2 results, for all compounds studied here, show that the βvec values are strongly affected by the effects of the electron correlation correction.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense pi-pi* transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Dynamic (hyper)polarizabilities of the ozone molecule: Coupled cluster calculations including vibrational corrections

In this work, we present results for dynamical (hyper)polarizabilities of the ozone molecule with inclusion of vibrational corrections. Electronic contributions for dynamic properties were computed analytically at the single and double coupled cluster level through response theories for the frequencies 0, 0.0239, 0.0428, and 0.0656 hartree. In the static limit, the electronic contributions were also computed at the single and double coupled cluster with perturbative correction of connected triple excitations level by means of the finite-field method. It was found that the inclusion of connected triple excitations is important, especially for a reliable description of the hyperpolarizabilities. Vibrational corrections were calculated by means of the perturbation theoretical method. The zero-point vibrational average correction was found to be relevant only for the linear polarizability, representing approximately 8% of the corresponding electronic contribution. Results also showed that the pure vibrational correction is relevant for the dc-Pockels effect, dc-second harmonic generation, intensity dependent refractive index, and dc-Kerr effect nonlinear optical processes. The double-harmonic approximation is in general suitable to compute this correction, the anharmonicity being small for the dc-Kerr effect and negligible for the other processes.

A theoretical investigation of electric properties of L-arginine phosphate monohydrate including environment polarization effects

The dipole moment (μ), linear polarizability (α), and first hyperpolarizability (β(tot)) of the asymmetric unit of L-arginine phosphate (LAP) monohydrate crystal are investigated using the supermolecule approach in combination with an iterative electrostatic polarization scheme. Environment polarization effects are attained by assuring the convergence of the dipole moment of LAP embedded in the polarization field of the surrounding molecules whose atomic sites are treated as point charges. The results obtained show that in the presence of the embedding charges, the value of μ is increased by 9% but the static values of α and β(tot) are decreased, respectively, by 3% and 13%, as compared with the isolated situation. The MP2/6-311+G(d) model predicts for the in-crystal dipole moment the converged value of 33 D, in good concordance with the available experimental result of 32 D. Our estimates for the converged results of α and β(tot) are, respectively, 22.51×10(-24) and 5.01×10(-30) esu. Dispersion effects are found to have a small impact on the nonlinear optical responses of LAP in the visible region. In addition, MP2/6-311G results obtained for β(tot) by using isolated and embedded LAP dimers show that crystal packing effects have a significant contribution of the electrostatic interactions. Our results suggest that the role of the crystal environment is to minimize the effects of the intermolecular interactions in the electric properties. That is, μ and β(tot) gain a more additive character in the presence of the field of the embedding charges. This is specially marked for β(tot).

Vibrational effects on the dynamic electric properties of hydrogen peroxide

In this work we present a method based on the perturbation theoretic approach of Bishop and co-workers [J. Chem. Phys. 95, 2646 (1991); 97, 5255 (1992); 108, 10013 (1998)] to calculate the effect of torsional motion on the polarizability and hyperpolarizabilities of hydrogen peroxide. The frequency dependence has been evaluated using the time-dependent Hartree-Fock method. The results obtained show that the zero-point vibrational averaging contributions are small compared to the corresponding electronic contributions. In the static limit the pure vibrational contributions are very large, specially for beta and gamma. These contributions are significant for the hyperpolarizabilities even in the visible region, except for the second harmonic generation and third harmonic generation processes.

Ab initio polarizabilities calculations of singly charged polyacetylene oligomers

Abstract We present results for the static longitudinal linear polarizability and second-order hyperpolarizability of small polyacetylene chains bearing positively and negatively charged solitons, obtained through the second-order Moller–Plesset perturbation theory (MP2) method. Hartree–Fock (HF) calculations for these properties was performed only for negatively charged chains. The standard 6-31G basis set was used in all calculations. Our ab initio calculations showed that, regarding singly charged structures, only the second hyperpolarizability is affected by the ionization state. For both, positive and negative structures, it is shown that the electron correlation effect enhances the linear polarizability, and even more markedly the second hyperpolarizabilities.

Polarization effects on the electric properties of urea and thiourea molecules in solid phase

We present theoretical results for the dipole moment, linear polarizability, and first hyperpolarizability of the urea and thiourea molecules in solid phase. The in-crystal electric properties were determined by applying a supermolecule approach in combination with an iterative electrostatic scheme, in which the surrounding molecules are represented by point charges. It is found for both urea and thiourea molecules that the influence of the polarization effects is mild for the linear polarizability, but it is marked for the dipole moment and first hyperpolarizability. The replacement of oxygen atoms by sulfur atoms increases, in general, the electric responses. Our second-order Møller-Plesset perturbation theory based iterative scheme predicts for the in-crystal dipole moment of urea and thiourea the values of 7.54 and 9.19 D which are, respectively, increased by 61% and 58%, in comparison with the corresponding isolated values. The result for urea is in agreement with the available experimental result of 6.56 D. In addition, we present an estimate of macroscopic quantities considering explicit unit cells of urea and thiourea crystals including environment polarization effects. These supermolecule calculations take into account partially the exchange and dispersion effects. The results illustrate the role played by the electrostatic interactions on the static second-order nonlinear susceptibility of the urea crystal.