The calculation of frequency-dependent hyperpolarizabilities including electron correlation effects

The theory for the calculation of the frequency-dependent hyperpolarizabilities β(?2ω; 0, ω), β(?ω; 0, ω), and β(0; ω, ?ω) is discussed. New relations between these tensors are derived for those wave functions that obey the time-dependent Hellmann–Feynman theorem (e.g., the self-consistent field [SCF] or the exact wave function). Using second-order Møller–Plesset perturbation theory (MP2), expressions are obtained for the hyperpolarizabilities in terms of derivatives of appropriately defined linear polarizability tensors with respect to a static electric field. Results are presented for ammonia and formaldehyde for the optical Kerr effect and for secondharmonic generation. These results indicate that it is desirable to determine the frequency-dependent contribution to the hyperopolarizability at the MP2 rather than the SCF level of theory, in cases where the static hyperpolarizability has a large contribution from electron correlation and/or where the frequency-dependent contribution may be more significant, such as for secondharmonic generation.     

A method for selection and use of numerical integration points for molecules with cylindrical symmetry

Boys and Handy [1] have discussed the solution of the bivariational equations with restricted numerical integration. One of the weaknesses of the method was that in the numerical summations over points, some points arose with r _ij= 0 and non-zero weights. This makes the method quite impractical for the Schrodinger Hamiltonian (because of the singularity at r _ij= 0), and it cannot be advantageous for the transcorrelated Hamiltonian C^–1HC because there will be some discontinuous higher derivatives at r _ij=0. Here it is shown how the symmetry of cylindrically symmetric molecules can be used to eliminate such points, without losing any of the advantages of the overall method, such as the convergence of the eigensolutions. It is also shown how the primary numerical integration points (z _i, r_i) may be chosen in any calculation such that each is associated with an equal amount of one-electron density. The choice of the angular coordinates are governed by the removal of the r _ij=0 points and maintaining the natural orthogonality between orbitals of different symmetry types. The method has been programmed and found to be practical, although no new molecular calculations have yet been performed. It is to be hoped that these points will give a basis for new transcorrelated calculations on diatomic molecules.This paper was presented during the session on numerical integration methods for molecules of the 1970 Quantum Theory Conference in Nottingham. It has been revised in the light of the interesting discussion which followed.     

Exact solution (within a double-zeta basis set) of the schrodinger electronic equation for water

Using recently developed theoretical techniques it has been possible to achieve an exact variational solution of the Schrödinger equation within a modest basis set or one-eleltron functions. The full configuration interaction for this system included a total of 256473 ¹A₁ spin-and space-adapted configurations. Comparison with many-body perturbation theory proves to be quite interesting.     

Density-functional generalized-gradient and hybrid calculations of electromagnetic properties using Slater basis sets

In this paper we extend our density-functional theory calculations, with generalized gradient approximation and hybrid functionals, using Slater-type orbitals (STOs), to the determination of second-order molecular properties. The key to the entire methodology involves the fitting of all STO basis function products to an auxiliary STO basis, through the minimization of electron-repulsion integrals. The selected properties are (i) dipole polarizabilities, (ii) nuclear magnetic shielding constants, and (iii) nuclear spin-spin coupling constants. In all cases the one-electron integrals involving STOs were evaluated by quadrature. The implementation for (ii) involved some complexity because we used gauge-including atomic orbitals. The presence of two-electron integrals on the right-hand side of the coupled equations meant that the fitting procedure had to be implemented. For (iii) in the hybrid case, fitting procedures were again required for the exchange contributions. For each property we studied a number of small molecules. We first obtained an estimate of the basis set limit using Gaussian-type orbitals (GTOs). We then showed how it is possible to reproduce these values using a STO basis set. For (ii) a regular TZ2P quality STO basis was adequate; for (i) the addition of one set of diffuse functions (determined by Slater's rules) gave the required accuracy; for (iii) it was necessary to add a set of 1s functions, including one very tight function, to give the desired result. In summary, we show that it is possible to predict second-order molecular properties using STO basis sets with an accuracy comparable with large GTO basis sets. We did not encounter any major difficulties with either the selection of the bases or the implementation of the procedures. Although the energy code (especially in the hybrid case) may not be competitive with a regular GTO code, for properties we find that STOs are more attractive.     

Ab initio studies of internal rotation barriers and vibrational frequencies of (C2H2)2, (CO2)2, and C2H2-CO2

We have performed large-scaleab initio calculations using second order Møller-Plesset perturbation theory (MP2) on the three van der Waals dimers formed from acetylene and carbon dioxide. Intermolecular geometrical parameters are reliably computed at this level of theory. Calculations of vibrational frequencies of the van der Waals modes, currently unobtainable by experimental means, give important information about the intermolecular potential and predict significant large-amplitude motion. Zero point energy contributions are shown to be vital in assessing the relative stability of conformations which are close in energy. Our studies suggest that the barrier to interconversion tunnelling in (CO₂)₂ is significantly smaller than previously inferred and is approximately the same as in (C₂H₂)₂. The reason for the rigidity of (CO₂)₂ is the difference in monomer centre-of-mass separation between ground state and transition state. We also show that, in addition to the previously observedC _2v form, the collinear form of C₂H₂-CO₂ is a local minimum on its potential energy surface.     

An acid-free Pictet–Spengler reaction using deep eutectic solvents (DES)

Deep eutectic solvents (such as the combination of either urea or glycerol with choline chloride) are effective solvents/organocatalysts for Pictet–Spengler condensations to form carbolines. The reaction conditions are quite mild and do not require additional Bronsted or Lewis acid catalyst. Given the inexpensive, non-toxic, and recyclable nature of the DES, these reaction conditions are simple and highly environmentally friendly.     

A mild synthesis of bis(indolyl)methanes using a deep eutectic solvent

Deep eutectic solvents (such as the combination of urea and choline chloride) are effective solvents/organocatalysts for the condensation of indole and aryl or alkenyl aldehydes to form bis(indolyl)methanes. The reaction conditions are quite mild and do not require additional Bronsted or Lewis acid catalyst, though they fail with ketones or aliphatic aldehydes. Given the inexpensive, non-toxic, and recyclable nature of the DES, these reaction conditions are simple and highly environmentally friendly.     

One-pot Double Suzuki Couplings of Dichloropyrimidines

An effective one-pot, regioselective double Suzuki coupling of 2,4-dichloropyrimidine has been developed, which enables the quick and efficient synthesis of diarylated pyrimidines. The choice of solvent proved critical to the success of this reaction sequence, with alcoholic solvent mixtures affording much greater reactivity and correspondingly lower temperatures than the use of polar aprotic solvents.