General derivative relations for anharmonic force fields

The brace notation, introduced by Allen and Csaszar (1993, J. chem. Phys., 98, 2983), provides a simple and compact way to deal with derivatives of arbitrary non-tensorial quantities. One of its main advantages is that it builds the permutational symmetry of the derivatives directly into the formalism. The brace notation is applied to formulate the general nth-order Cartesian derivatives of internal coordinates, and to provide closed forms for general, nth-order transformation equations of anharmonic force fields, expressed as Taylor series, from internal to Cartesian or normal coordinate spaces.

The two-dimensional anharmonic oscillator I: vibration-rotation constants of fulminic acid, HCNO

The lowest-wavenumber vibration of HCNO and DCNO, ν5, is known to involve a largeamplitude low-frequency anharmonic bending of the CH bond against the CNO frame. In this paper the anomalous vibrational dependence of the observed rotational constants B(v5, l5), and of the observed l-doubling interactions, is interpreted according to a simple effective vibration-rotation Hamiltonian in which the appropriate vibrational operators are averaged in an anharmonic potential surface over the normal coordinates (Q5x, Q5y). All of the data on both isotopes are interpreted according to a single potential surface having a minimum energy at a slightly bent configuration of the HCN angle ( 170°) with a maximum at the linear configuration about 2 cm−1 higher. The other coefficients in the Hamiltonian are also interpreted in terms of the structure and the harmonic and anharmonic force fields; the substitution structure at the “hypothetical linear configuration” determined in this way gives a CH bond length of 1.060 Å, in contrast to the value 1.027 Å determined from the ground-state rotational constants. We also discuss the difficulties in rationalizing our effective Hamiltonian in terms of more fundamental theory, as well as the success and limitations of its use in practice.

Vibrational intensities VIII: CH3- and CD3- chloride, bromide and iodide

Absolute intensity measurements have been made on the fundamental vibrations of methyl chloride, bromide, and iodide, and their fully deuterated derivatives, by integrating the optical density over the absorption bands. The bands were fully pressure broadened by using up to 80 atmos of foreign gas. Band separations were made graphically. The results are analyzed in terms of the dipole moment derivatives with respect to symmetry coordinates in the molecule, (∂p/∂Si). The data on the different isotopic species are shown to yield consistent results, and this requirement of consistency has also been used as an aid in the analysis. In the E‐class vibrations the signs of the dipole moment derivatives have been determined unambiguously by assuming the permanent dipole to be directed CH3+☒X—.

A potential energy surface for the ground state of formaldehyde, H2CO(1A1)

A model potential energy function for the ground state of H2CO has been derived which covers the whole space of the six internal coordinates. This potential reproduces the experimental energy, geometry and quadratic force field of formaldehyde, and dissociates correctly to all possible atom, diatom and triatom fragments. Thus there are good reasons for believing it to be close to the true potential energy surface except in regions where both hydrogen atoms are close to the oxygen. It leads to the prediction that there should be a metastable singlet hydroxycarbene HCOH which has a planar trans structure and an energy of 2•31 eV above that of equilibrium formaldehyde. The reaction path for dissociation into H2 + CO is predicted to pass through a low symmetry transition state with an activation energy of 4•8 eV. Both of these predictions are in good agreement with recently published ab initio calculations.

ASYM20: a program for force constant and normal coordinate calculations, with a critical review of the theory involved

The theory of harmonic force constant refinement calculations is reviewed, and a general-purpose program for force constant and normal coordinate calculations is described. The program, called ASYM20. is available through Quantum Chemistry Program Exchange. It will work on molecules of any symmetry containing up to 20 atoms and will produce results on a series of isotopomers as desired. The vibrational secular equations are solved in either nonredundant valence internal coordinates or symmetry coordinates. As well as calculating the (harmonic) vibrational wavenumbers and normal coordinates, the program will calculate centrifugal distortion constants, Coriolis zeta constants, harmonic contributions to the α′s. root-mean-square amplitudes of vibration, and other quantities related to gas electron-diffraction studies and thermodynamic properties. The program will work in either a predict mode, in which it calculates results from an input force field, or in a refine mode, in which it refines an input force field by least squares to fit observed data on the quantities mentioned above. Predicate values of the force constants may be included in the data set for a least-squares refinement. The program is written in FORTRAN for use on a PC or a mainframe computer. Operation is mainly controlled by steering indices in the input data file, but some interactive control is also implemented.

Variational calculations of rovibrational states: a precise high-energy potential surface for HCN [and discussion]

We report the results of variational calculations of the rovibrational energy levels of HCN for J = 0, 1 and 2, where we reproduce all the ca. 100 observed vibrational states for all observed isotopic species, with energies up to 18000 cm$^{-1}$, to about $\pm $1 cm$^{-1}$, and the corresponding rotational constants to about $\pm $0.001 cm$^{-1}$. We use a hamiltonian expressed in internal coordinates r$_{1}$, r$_{2}$ and $\theta $, using the exact expression for the kinetic energy operator T obtained by direct transformation from the cartesian representation. The potential energy V is expressed as a polynomial expansion in the Morse coordinates y$_{i}$ for the bond stretches and the interbond angle $\theta $. The basis functions are built as products of appropriately scaled Morse functions in the bond-stretches and Legendre or associated Legendre polynomials of cos $\theta $ in the angle bend, and we evaluate matrix elements by Gauss quadrature. The hamiltonian matripx is factorized using the full rovibrational symmetry, and the basis is contracted to an optimized form; the dimensions of the final hamiltonian matrix vary from 240 $\times $ 240 to 1000 $\times $ 1000.We believe that our calculation is converged to better than 1 cm$^{-1}$ at 18 000 cm$^{-1}$. Our potential surface is expressed in terms of 31 parameters, about half of which have been refined by least squares to optimize the fit to the experimental data. The advantages and disadvantages and the future potential of calculations of this type are discussed.

Friends or relatives? Phylogenetics and species delimitation in the controversial European orchid Genus ophrys

Background and Aims Highly variable, yet possibly convergent, morphology and lack of sequence variation have severely hindered production of a robust phylogenetic framework for the genus Ophrys. The aim of this study is to produce this framework as a basis for more rigorous species delimitation and conservation recommendations. Methods Nuclear and plastid DNA sequencing and amplified fragment length polymorphism (AFLP) were performed on 85 accessions of Ophrys, spanning the full range of species aggregates currently recognized. Data were analysed using a combination of parsimony and Bayesian tree-building techniques and by principal coordinates analysis. Key Results Complementary phylogenetic analyses and ordinations using nuclear, plastid and AFLP datasets identify ten genetically distinct groups (six robust) within the genus that may in turn be grouped into three sections (treated as subgenera by some authors). Additionally, genetic evidence is provided for a close relationship between the O. tenthredinifera, O. bombyliflora and O. speculum groups. The combination of these analytical techniques provides new insights into Ophrys systematics, notably recognition of the novel O. umbilicata group. Conclusions Heterogeneous copies of the nuclear ITS region show that some putative Ophrys species arose through hybridization rather than divergent speciation. The supposedly highly specific pseudocopulatory pollination syndrome of Ophrys is demonstrably 'leaky', suggesting that the genus has been substantially over-divided at the species level.

Rovibrational energy levels of hydrogen peroxide, studied by MULTIMODE with a reaction path Hamiltonian

Recently. Carter and Handy [J. Chem. Phys. 113 (2000) 987] have introduced the theory of the reaction path Hamiltonian (RPH) [J. Chem. Phys. 72 (1980) 99] into the variational scheme MULTIMODE, for the calculation of the J = 0 vibrational levels of polyatomic molecules, which have a single large-amplitude motion. In this theory the reaction path coordinate s becomes the fourth dimension of the moment-of-inertia tensor, and must be treated separately from the remaining 3N - 7 normal coordinates in the MULTIMODE program. In the modified program, complete integration is performed over s, and the M-mode MULTIMODE coupling approximation for the evaluation of the matrix elements applies only to the 3N - 7 normal coordinates. In this paper the new algorithm is extended to the calculation of rotational-vibration energy levels (i.e. J > 0) with the RPH, following from our analogous implementation for rigid molecules [Theoret. Chem. Acc. 100 (1998) 191]. The full theory is given, and all extra terms have been included to give the exact kinetic energy operator. In order to validate the new code, we report studies on hydrogen peroxide (H2O2), where the reaction path is equivalent to torsional motion. H2O2 has previously been studied variationally using a valence coordinate Hamiltonian; complete agreement for calculated rovibrational levels is obtained between the previous results and those from the new code, using the identical potential surface. MULTIMODE is now able to calculate rovibrational levels for polyatomic molecules which have one large-amplitude motion. (C) 2003 Elsevier B.V. All rights reserved.

Bis[bis(methoxycarbimido)aminato]copper(II) 1-methylpyrrolidin-2-one disolvate

The title compound, [Cu(C4H8N3O2)(2)]center dot 2C(5)H(9)NO, consists of a neutral copper complex, in which the Cu II centre coordinates to two bis(methoxycarbimido) aminate ligands, solvated by two molecules of 1-methylpyrrolidin-2-one. The complex is planar and centrosymmetric, with the Cu II centre occupying a crystallographic inversion centre and adopting approximately square-planar geometry. N-H center dot center dot center dot O hydrogen-bonding interactions exist between the amine NH groups of the ligands and the O atoms of the 1-methylpyrrolidin-2-one molecules. The associated units pack to form sheets.

Theoretical rovibrational line intensities in the electronic ground state of ozone

First-principles calculations of absolute line intensities and rovibrational energies of ozone (O-16(3)) are reported using potential energy and electric dipole moment functions calculated by the internally contracted MRCI approach. The rovibrational energies and eigenfunctions (up to about 8500 cm(-1) and J = 64) were obtained variationally with an exact Hamiltonian in internal valence coordinates. More than 4.8 x 10(6) electric dipole transition matrix elements were calculated for the absolute rovibrational line intensities. They are compared with the values of the HITRAN database. The purely rotational absolute line intensities in the (000) state and the rovibrational intensities for the (001)-(000) band agree to within about 0.3 to 1% for the (0 10)-(000) band to within about 3 to 4%. Excellent agreement with experiment is also achieved for low-lying overtone and combination bands. Inconsistencies are found for the (100)-(000) band overlapping with the antisymmetric stretching fundamental and also for the (002)-(000) antisymmetric stretching overtone. The generated dipole moment function can be used for predicting the absorption intensities in any of the heavier isotopomers, hot bands or the rates of spontaneous emission.

Structural diversity in supramolecular architectures of Mn(II) through non-covalent interactions

Two mononuclear complexes of manganese(II), [Mn(OCN)(2)(phen)(2)] 1 and [Mn(NCO)(2)(bpy)(2)] 2 [1,10-phenanthroline (phen); 2,2'-bipyridine (bpy)], have been synthesized and characterized by single crystal X-ray analysis, infra-red spectroscopy and magnetic studies. The coordination structure of complex 2 is already reported. The cyanate anions are pendent in both the complexes. In 1, cyanate anion links manganese(II) through O-atom, whereas in 2 it coordinates through N-atom. The mononuclear fragments of 1 are built up to a supramolecular lamellar 3D architecture by pi-pi interactions only. On the other hand, mononuclear fragments of 2 are assembled to a 2D supramolecular brick-wall architecture by C-H-... pi interactions.

Large vibrational variational calculations using 'multimode' and an iterative diagonalization technique

We have favoured the variational (secular equation) method for the determination of the (ro-) vibrational energy levels of polyatomic molecules. We use predominantly the Watson Hamiltonian in normal coordinates and an associated given potential in the variational code 'Multimode'. The dominant cost is the construction and diagonalization of matrices of ever-increasing size. Here we address this problem, using pertubation theory to select dominant expansion terms within the Davidson-Liu iterative diagonalization method. Our chosen example is the twelve-mode molecule methanol, for which we have an ab initio representation of the potential which includes the internal rotational motion of the OH group relative to CH3. Our new algorithm allows us to obtain converged energy levels for matrices of dimensions in excess of 100 000.

Vibrations in the B-4 rhombic structure

A double minimum six-dimensional Potential energy surface (PES) is determined in symmetry coordinates for the most stable rhombic (D-2h) B-4 isomer in its (1)A(g) electronic ground state by fitting to energies calculated ab initio. The PES exhibits a barrier to the D-4h square structure of 255 cm(-1). The vibrational levels (J=0) are calculated variationally using an approach which involves the Watson kinetic energy operator expressed in normal coordinates. The pattern of about 65 vibrational levels up to 1600 cm-1 for all stable isotopomers is analyzed. Analogous to the inversion in ammonia-like molecules, the rhombus rearrangements lead to splittings of the vibrational levels. In B-4 it is the B-1g (D-4h mode which distorts the square molecule to its planar rhombic form. The anharmonic fundamental vibrational transitions of B-11(4) are calculated to be (splittings in parentheses): G(O) = 2352(22) cm(-1), v(1)(A(1g)) - 1136(24) cm(-1,) v(2)(B-1g)=209(144) cm(-1) v(3)(B-2g)=1198(19)cm(-1), v(4)(B-2u) = 271(24) cm(-1), and v(5) (E-u) = 1030( 166) cm(-1) (D-4h notation). Their variations in all stable isotoporners were investigated. Due to the presence of strong anharmonic resonances between the B-1g in-plane distortion and the B-2u, out-of-plane bending modes. the hiaher overtones and combination levels are difficult to assign unequivocally. (C) 2005 American Institute of Physics.

Influence of counter anions on the structures and magnetic properties of trinuclear Cu(II) complexes containing a mu(3)-OH core and Schiff base ligands

Four trinuclear Cu(II) complexes, [(CuL1)(3)(mu(3)-OH)](NO3)(2) (1), [(CuL2)(3)(mu(3)-OH)](I)(2)center dot H2O (2), [(CuL3)(3)(mu(3)-OH)](I)(2) (3) and [(CuL1)(3)(mu(3)-OH)][(CuI3)-I-1] (4), where HL1 (8-amino-4-methyl-5-azaoct-3-en-2-one), HL2 [7-amino-4-methyl-5-azaoct-3-en-2-one] and HL3 [7-amino-4-methyl-5-azahept-3-en-2- one] are the three tridentate Schiff bases, have been synthesized and structurally characterized by X-ray crystallography. All four complexes contain a partial cubane core, [(CuL)(3)(mu(3)-OH)](2+) in which the three [CuL] subunits are interconnected through two types of oxygen bridges afforded by the oxygen atoms of the ligands and the central OH- group. The copper(II) ions are in a distorted square-pyramidal environment. The equatorial plane consists of the bridging oxygen of the central OH- group together with three atoms (N, N, O) from the Schiff base. The oxygen atom of the Schiff base also coordinates to the axial position of Cu(II) of another subunit to form the cyclic trimer. Magnetic susceptibilities have been determined for these complexes over the temperature range of 2-300 K. The isotropic Hamiltonian, H = -J(12)S(1)S(2) - J(13)S(1)S(3) - J(23)S(2)S(3) has been used to interpret the magnetic data. The best fit parameters obtained are: J = - 54.98 cm(-1) g = 2.24 for 1; J = - 56.66 cm(-1), g = 2.19 for 2; J = -44.39 cm(-1), g = 2.16 for 3; J = - 89.92 cm(-1), g = 2.25 for 4. The EPR data at low temperature indicate that the phenomenon of spin frustration occurs for complexes 1-3. (c) 2007 Elsevier B.V. All rights reserved.

Lanthanide complexes of 2-hydroxynicotinic acid: synthesis, luminescence properties and the crystal structures of [Ln(HnicO)(2)(mu-HnicO)(H2O)] center dot nH(2)O (Ln = Tb, Eu)

New lanthanide complexes of 2-hydroxynicotinic acid (H(2)nicO) [Ln(HnicO)(2)(mu-HnicO)(H2O)] (.) nH(2)O (Ln = Eu, Gd, Tb, Er, Tm) were prepared. The crystal structures of the [Tb(HnicO)(2)(g-HnicO)(H2O)] (.) 1.75H(2)O(1) and [Eu(HniCO)(2)(mu-HnicO)(H2O)] (.) 1.25H(2)O (2) complexes were determined by X-ray diffraction. The 2-hydroxynicotinate ligand coordinates through O,O-chelation to the lanthanide(III) ions as shown by X-ray diffraction and the infrared, Raman and NMR spectroscopy results. Photoluminescence measurements were performed for the Eu(III) and Tb(III) complexes. Lifetimes of 0.592 +/- 0.007 and 0.113 +/- 0.002 ms were determined for the Eu3+ and Tb3+ emitting states D-5(0) and D-5(4), respectively. A value around 30% was found for the D-5(0) quantum efficiency. The energy transfer mechanisms between the lanthanide ions and the ligands are discussed and compared with those observed in similar complexes involving the 3-hydroxypicolinate ligand based on the luminescence of the respective Gd3+-based complexes. (C) 2003 Published by Elsevier Ltd.