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

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

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

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

Variational calculation of static and dynamic vibrational nonlinear optical properties

The vibrational configuration interaction method used to obtain static vibrational (hyper)polarizabilities is extended to dynamic nonlinear optical properties in the infinite optical frequency approximation. Illustrative calculations are carried out on H2 O and N H3. The former molecule is weakly anharmonic while the latter contains a strongly anharmonic umbrella mode. The effect on vibrational (hyper)polarizabilities due to various truncations of the potential energy and property surfaces involved in the calculation are examined

Calculation of Franck-Condon factors including anharmonicity: simulation of the C2 H4+ X̃2B 3u←C2 H4X̃1 Ag band in the photoelectron spectrum of ethylene

Our new simple method for calculating accurate Franck-Condon factors including nondiagonal (i.e., mode-mode) anharmonic coupling is used to simulate the C2H4+X2B 3u←C2H4X̃1 Ag band in the photoelectron spectrum. An improved vibrational basis set truncation algorithm, which permits very efficient computations, is employed. Because the torsional mode is highly anharmonic it is separated from the other modes and treated exactly. All other modes are treated through the second-order perturbation theory. The perturbation-theory corrections are significant and lead to a good agreement with experiment, although the separability assumption for torsion causes the C2 D4 results to be not as good as those for C2 H4. A variational formulation to overcome this circumstance, and deal with large anharmonicities in general, is suggested

The aromatic fluctuation index (FLU): a new aromaticity index based on electron delocalization

In this work, the aromatic fluctuation index (FLU) that describes the fluctuation of electronic charge between adjacent atoms in a given ring is introduced as a new aromaticity measure. This new electronic criterion of aromaticity is based on the fact that aromaticity is related to the cyclic delocalized circulation of π electrons. It is defined not only considering the amount of electron sharing between contiguous atoms, which should be substantial in aromatic molecules, but also taking into account the similarity of electron sharing between adjacent atoms. For a series of rings in 15 planar polycyclic aromatic hydrocarbons, we have found that, in general, FLU is strongly correlated with other widely used indicators of local aromaticity, such as the harmonic-oscillator model of aromaticity, the nucleus independent chemical shift, and the para-delocalization index (PDI). In contrast to PDI, the FLU index can be applied to study the aromaticity of rings with any number of members and it can be used to analyze both the local and global aromatic character of rings and molecules

Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrödinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational Møller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed

Basis set and electron correlation effects on initial convergence for vibrational nonlinear optical properties of conjugated organic molecules

The level of ab initio theory which is necessary to compute reliable values for the static and dynamic (hyper)polarizabilities of three medium size π-conjugated organic nonlinear optical (NLO) molecules is investigated. With the employment of field-induced coordinates in combination with a finite field procedure, the calculations were made possible. It is stated that to obtain reasonable values for the various individual contributions to the (hyper)polarizability, it is necessary to include electron correlation. Based on the results, the convergence of the usual perturbation treatment for vibrational anharmonicity was examined

Electronic and vibrational contributions to first hyperpolarizability of donor-acceptor-substituted azobenzene

In this study we report on the electronic and vibrational (hyper)polarizabilities of donor–acceptorsubstituted azobenzene. It is observed that both electronic and vibrational contributions to the electric dipole first hyperpolarizability of investigated photoactive molecule substantially depend on the conformation. The contributions to the nuclear relaxation first hyperpolarizability are found to be quite important in the case of two considered isomers (cis and trans). Although the double-harmonic term is found to be the largest in terms of magnitude, it is shown that the total value of the nuclear relaxation contribution to vibrational first hyperpolarizability is a result of subtle interplay of higher-order contributions. As a part of the study, we also assess the performance of long-range-corrected density functional theory in determining vibrational contributions to electric dipole (hyper)polarizabilities. In most cases, the applied long-range-corrected exchange correlation potentials amend the drawbacks of their conventional counterparts

Basin-scale hydrodynamics in a Mediterranean reservoir. Implications for the phytoplankton dynamics

Procesos hidrodinámicos determinan, en un alto grado la calidad del agua en embalse, sin embargo dichos procesos han sido tradicionalmente olvidados en la gestión de embalse. En esta tesis se presentan evidencias de los principales procesos hidrodinámicos que ocurren en un embalse Mediterráneo a escala de cuenca a través de campañas experimentales y modelización numérica; y su influencia en la dinámica de poblaciones de fitoplancton. Dichos procesos son principalmente la generación de ondas internas o secas y la intrusión del río. La presencia de viento periódico genera secas forzadas, amplificando los modos cercanos al periodo del viento, de manera que modos verticales altos, considerados como raros en la naturaleza, tienden a dominar en el sistema.