Numerically Exact Study of Polarizabilities and Hyperpolarizabilities of Correlated Conjugated Organic Models

Polarizabilities and Hyperpolarizabilities of conjugated organic chains are calculated using correlated model Hamiltonians. While correlations reduce the Polarizabilities and extend the range of linear response, the Hyperpolarizabilities essentially are unaffected by the same. This explains the apparently large Hyperpolarizabilities of conjugated electronic systems.

Exact static polarizabilities of correlated finite model systems

A finite-field method for calculating exact polarizabilities of correlated conjugated model systems within the valence bond (VB) framework is presented. The correlations reduce the polarizabilities from their noninteracting values and extend the range of linearity to higher external fields. The large nonlinear polarizabilities observed in strongly correlated conjugated organic molecules cannot be directly attributed to electron correlations. The method described can be employed to calculate static polarizabilities for any desired state of a correlated system.

Electric and magnetic polarizabilities of hexagonal Ln(2)CuTiO(6)(Ln = Y, Dy, Ho, Er, and Yb)

We investigated the rare-earth transition-metal oxide series, Ln(2)CuTiO(6) (Ln = Y, Dy, Ho, Er, and Yb), crystallizing in the hexagonal structure with noncentrosymmetric P6(3)cm space group for possible occurrences of multiferroic properties. Our results show that while these compounds, except Ln = Y, exhibit a low-temperature antiferromagnetic transition due to the ordering of the rare-earth moments, the expected ferroelectric transition is frustrated by the large size difference between Cu and Ti at the B site. Interestingly, this leads these compounds to attain a rare and unique combination of desirable paraelectric properties with high dielectric constants, low losses, and weak temperature and frequency dependencies. First-principles calculations establish these exceptional properties result from a combination of two effects. A significant difference in the MO5 polyhedral sizes for M = Cu and M = Ti suppress the expected cooperative tilt pattern of these polyhedra, required for the ferroelectric transition, leading to relatively large values of the dielectric constant for every compound investigated in this series. Additionally, it is shown that the majority contribution to the dielectric constant arises from intermediate-frequency polar vibrational modes, making it relatively stable against any temperature variation. Changes in the temperature stability of the dielectric constant among different members of this series are shown to arise from changes in relative contributions from soft polar modes.

1. Ultrasonic studies of binary liquid structure in the critical region. Theory and experiment for the 2,6-lutidine/water system. 2. Hartree-Fock calculations of electric polarizabilities of some simple atoms and molecules, and their practicality. 3. Calculation of vibrational transition probabilities in collinear atom-diatom and diatom-diatom collisions with Lennard-Jones interaction

Part 1. Many interesting visual and mechanical phenomena occur in the critical region of fluids, both for the gas-liquid and liquid-liquid transitions. The precise thermodynamic and transport behavior here has some broad consequences for the molecular theory of liquids. Previous studies in this laboratory on a liquid-liquid critical mixture via ultrasonics supported a basically classical analysis of fluid behavior by M. Fixman (e. g., the free energy is assumed analytic in intensive variables in the thermodynamics)--at least when the fluid is not too close to critical. A breakdown in classical concepts is evidenced close to critical, in some well-defined ways. We have studied herein a liquid-liquid critical system of complementary nature (possessing a lower critical mixing or consolute temperature) to all previous mixtures, to look for new qualitative critical behavior. We did not find such new behavior in the ultrasonic absorption ascribable to the critical fluctuations, but we did find extra absorption due to chemical processes (yet these are related to the mixing behavior generating the lower consolute point). We rederived, corrected, and extended Fixman's analysis to interpret our experimental results in these more complex circumstances. The entire account of theory and experiment is prefaced by an extensive introduction recounting the general status of liquid state theory. The introduction provides a context for our present work, and also points out problems deserving attention. Interest in these problems was stimulated by this work but also by work in Part 3.

Part 2. Among variational theories of electronic structure, the Hartree-Fock theory has proved particularly valuable for a practical understanding of such properties as chemical binding, electric multipole moments, and X-ray scattering intensity. It also provides the most tractable method of calculating first-order properties under external or internal one-electron perturbations, either developed explicitly in orders of perturbation theory or in the fully self-consistent method. The accuracy and consistency of first-order properties are poorer than those of zero-order properties, but this is most often due to the use of explicit approximations in solving the perturbed equations, or to inadequacy of the variational basis in size or composition. We have calculated the electric polarizabilities of H₂, He, Li, Be, LiH, and N₂ by Hartree-Fock theory, using exact perturbation theory or the fully self-consistent method, as dictated by convenience. By careful studies on total basis set composition, we obtained good approximations to limiting Hartree-Fock values of polarizabilities with bases of reasonable size. The values for all species, and for each direction in the molecular cases, are within 8% of experiment, or of best theoretical values in the absence of the former. Our results support the use of unadorned Hartree-Pock theory for static polarizabilities needed in interpreting electron-molecule scattering data, collision-induced light scattering experiments, and other phenomena involving experimentally inaccessible polarizabilities.

Part 3. Numerical integration of the close-coupled scattering equations has been carried out to obtain vibrational transition probabilities for some models of the electronically adiabatic H₂-H₂ collision. All the models use a Lennard-Jones interaction potential between nearest atoms in the collision partners. We have analyzed the results for some insight into the vibrational excitation process in its dependence on the energy of collision, the nature of the vibrational binding potential, and other factors. We conclude also that replacement of earlier, simpler models of the interaction potential by the Lennard-Jones form adds very little realism for all the complication it introduces. A brief introduction precedes the presentation of our work and places it in the context of attempts to understand the collisional activation process in chemical reactions as well as some other chemical dynamics.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.

Investigation of molecular polarizabilities and derivatives in halomethanes

The research undertaken was to obtain absolute Raman intensities for the symmetric stretching vibrations of the methyl halides, CH3X with (X=F, CI, Br), by experiment and theory. The intensities were experimentally measured using the Ar+ ion gas laser as excitation source, a Spex 14018 double monochromator and a RCA C-31034 photomultiplier tube as detector. These intensities arise from changes in the derivative of the polarizability (8 a'), with respect to vibration along a normal coordinate (8qi). It was intended that these derivatives obtained with respect to normal coordinates would be converted to derivatives with respect to internal coordinates, for a quantitative comparison with theory. Theoretical numerical polarizability derivatives for the stretching vibrations are obtained using the following procedure. A vibration was simulated in the molecule by increasi.ng and decreasing the respective bond by the amount ±o.oosA for the C-H bonds and ±o.oIA for the C-X (X=F, CI, Br) bond. The derivative was obtained by taking the difference in the polarizability for the equilibrium geometry and the geometry when a particular bond is changed. This difference, when divided by the amount of change in each bond and the number of bonds present results in the derivative of the polarizability with respect to internal coordinate i.e., !1u/!1r. These derivatives were obtained by two methods: I} ab initio molecular orbital calculation and 2} theory of atoms in molecules (AIM) analysis. Due to errors in the experimental setup only a qualitative analysis of the results was undertaken relative to the theory. Theoretically it is predicted that the symmetric carbonhalogen stretch vibrations are more intense than the respective carbon-hydrogen stretch, but only for the methyl chloride and bromide. The carbon fluorine stretch is less intense than the carbon-hydrogen stretch, a fact which is attributed to the small size and high electronegativity of the fluorine atom. The experimental observations are seen to agree qualitatively with the theory results. It is hoped that when the experiment is repeated, a quantitative comparison can be made. The analysis by the theory of atoms in molecules, along with providing polarizabilities and polarizability derivatives, gives additional information outlined below. The theory provides a pictorial description of the main factors contributing to the molecular polarizability and polarizability derivative. These contributions are from the charge transfer and atomic dipole terms i.e., transfer of charge from one atom to another and the reorganization of atomic electronic charge distribution due to presence of an electric field. The linear relationship between polarizability and molecular volume was also observed.

Quantum Monte Carlo study of electrostatic polarizabilities of H and He atoms /

The infinitesimal differential quantum Monte Carlo (QMC) technique is used to estimate electrostatic polarizabilities of the H and He atoms up to the sixth order in the electric field perturbation. All 542 different QMC estimators of the nonzero atomic polarizabilities are derived and used in order to decrease the statistical error and to obtain the maximum efficiency of the simulations. We are confident that the estimates are "exact" (free of systematic error): the two atoms are nodeless systems, hence no fixed-node error is introduced. Furthermore, we develope and use techniques which eliminate systematic error inherent when extrapolating our results to zero time-step and large stack-size. The QMC results are consistent with published accurate values obtained using perturbation methods. The precision is found to be related to the number of perturbations, varying from 2 to 4 significant digits.

On the correlation between electric polarizabilities and the ionization potential of atoms

It is found that the electric dipole polarizabilities of neutral atoms correlate very strongly with their first ionization potential within the groups of elements with the same angular momenta of the outermost electrons. As the latter values are known very accurately, this allows a very good (<30%) prediction of various atomic polarizabilities.

Additional compact formulas for vibrational dynamic dipole polarizabilities and hyperpolarizabilities

Compact expressions, complete through second order in electrical and/or mechanical anharmonicity, are given for the dynamic dipole vibrational polarizability and dynamic first and second vibrational hyperpolarizabilities. Certain contributions not previously formulated are now included

Finite field treatment of vibrational polarizabilities and hyperpolarizabilities: on the role of the Eckart conditions, their implementation, and their use in characterizing key vibrations

In the finite field (FF) treatment of vibrational polarizabilities and hyperpolarizabilities, the field-free Eckart conditions must be enforced in order to prevent molecular reorientation during geometry optimization. These conditions are implemented for the first time. Our procedure facilities identification of field-induced internal coordinates that make the major contribution to the vibrational properties. Using only two of these coordinates, quantitative accuracy for nuclear relaxation polarizabilities and hyperpolarizabilities is achieved in π-conjugated systems. From these two coordinates a single most efficient natural conjugation coordinate (NCC) can be extracted. The limitations of this one coordinate approach are discussed. It is shown that the Eckart conditions can lead to an isotope effect that is comparable to the isotope effect on zero-point vibrational averaging, but with a different mass-dependence

Determination of vibrational polarizabilities and hyperpolarizabilities using field-induced coordinates

An analytical set of field-induced coordinates is defined and is used to show that the vibrational degrees of freedom required to completely describe nuclear relaxation polarizabilities and hyperpolarizabilities is reduced from 3N-6 to a relatively small number. As this number does not depend upon the size of the molecule, the process provides computational advantages. A method is provided to separate anharmonic contributions from harmonic contributions as well as effective mechanical from electrical anharmonicity. The procedures are illustrated by Hartree-Fock calculations, indicating that anharmonicity can be very important

High precision calculation of multipolar dynamic polarizabilities and two- and three-body dispersion coefficients of atomic hydrogen

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)