Calculation of converged rovibrational energies and partition function for methane using vibrational-rotational configuration interaction

The rovibration partition function of CH4 was calculated in the temperature range of 100-1000 K using well-converged energy levels that were calculated by vibrational-rotational configuration interaction using the Watson Hamiltonian for total angular momenta J=0-50 and the MULTIMODE computer program. The configuration state functions are products of ground-state occupied and virtual modals obtained using the vibrational self-consistent field method. The Gilbert and Jordan potential energy surface was used for the calculations. The resulting partition function was used to test the harmonic oscillator approximation and the separable-rotation approximation. The harmonic oscillator, rigid-rotator approximation is in error by a factor of 2.3 at 300 K, but we also propose a separable-rotation approximation that is accurate within 2% from 100 to 1000 K. (C) 2004 American Institute of Physics.

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.

Harmonic force fields from scaled SCF calculations: Program ASYM40

We report an extended version of our normal coordinate program ASYM40, which may be used to transform Cartesian force constants from ab initio calculations to a force field in nonredundant internal (symmetry) coordinates. When experimental data are available, scale factors for the theoretical force field may then be optimized by least-squares refinement. The alternative of refining an empirical force field to fit a wide variety of data, as with the previous version ASYM20, has been retained. We compare the results of least-squares refinement of the full harmonic force field with least-squares refinement of only the scale factors for an SCF calculated force field and conclude that the latter approach may be useful for large molecules where more sophisticated calculations are impractical. The refinement of scale factors for a theoretical force field is also useful when there are only limited spectroscopic data. The program will accept ab initio calculated force fields from any program that presents Cartesian force constants as output. The program is available through Quantum Chemistry Program Exchange.

The two-dimensional anharmonic oscillator III: the CCC bending mode of C3O2

Newly observed data on the rotational constants of carbon suboxide in excited vibrational states of the low-wavenumber bending vibration ν7 have been successfully interpreted in terms of the two-dimensional anharmonic oscillator wavefunctions associated with this vibration. By combining these results with published infrared and Raman spectra the vibrational assignment has been extended and a refined bending potential for ν7 has been derived: this has a minimum at a bending angle of about 24° at the central C atom, with an energy maximum at the linear configuration some 23 cm−1 above the minimum. From similar data on the combination and hot bands of ν7 with ν4 (1587 cm−1) and ν2 (786 cm−1) the effective ν7 bending potential has also been determined in the one-quantum excited states of ν4 and ν2. The effective ν7 potential shows significant changes from the ground vibrational state; the central hump in the ν7 potential surface is increased to about 50 cm−1 in the v4 = 1 state, and decreased to about 1 cm−1 in the v2 = 1 state. In the light of these results vibrational assignments are suggested for most of the observed bands in the infrared and Raman spectra of C3O2.

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.

Wigner Transport models of the electron-phonon kinetics in quantum wires

Two quantum-kinetic models of ultrafast electron transport in quantum wires are derived from the generalized electron-phonon Wigner equation. The various assumptions and approximations allowing one to find closed equations for the reduced electron Wigner function are discussed with an emphasis on their physical relevance. The models correspond to the Levinson and Barker-Ferry equations, now generalized to account for a space-dependent evolution. They are applied to study the quantum effects in the dynamics of an initial packet of highly nonequilibrium carriers, locally generated in the wire. The properties of the two model equations are compared and analyzed.

Assessment of the RE(OH)3 Ising magnetic materials as possible candidates for the study of transverse-field-induced quantum phase transitions

The LiHoxY1−xF4 Ising magnetic material subject to a magnetic field perpendicular to the Ho3+ Ising direction has shown over the past 20 years to be a host of very interesting thermodynamic and magnetic phenomena. Unfortunately, the availability of other magnetic materials other than LiHoxY1−xF4 that may be described by a transverse-field Ising model remains very much limited. It is in this context that we use here a mean-field theory to investigate the suitability of the Ho(OH)3, Dy(OH)3, and Tb(OH)3 insulating hexagonal dipolar Ising-type ferromagnets for the study of the quantum phase transition induced by a magnetic field, Bx, applied perpendicular to the Ising spin direction. Experimentally, the zero-field critical (Curie) temperatures are known to be Tc≈2.54, 3.48, and 3.72 K, for Ho(OH)3, Dy(OH)3, and Tb(OH)3, respectively. From our calculations we estimate the critical transverse field, Bxc, to destroy ferromagnetic order at zero temperature to be Bxc=4.35, 5.03, and 54.81 T for Ho(OH)3, Dy(OH)3, and Tb(OH)3, respectively. We find that Ho(OH)3, similarly to LiHoF4, can be quantitatively described by an effective S=1/2 transverse-field Ising model. This is not the case for Dy(OH)3 due to the strong admixing between the ground doublet and first excited doublet induced by the dipolar interactions. Furthermore, we find that the paramagnetic (PM) to ferromagnetic (FM) transition in Dy(OH)3 becomes first order for strong Bx and low temperatures. Hence, the PM to FM zero-temperature transition in Dy(OH)3 may be first order and not quantum critical. We investigate the effect of competing antiferromagnetic nearest-neighbor exchange and applied magnetic field, Bz, along the Ising spin direction ẑ on the first-order transition in Dy(OH)3. We conclude from these preliminary calculations that Ho(OH)3 and Dy(OH)3 and their Y3+ diamagnetically diluted variants, HoxY1−x(OH)3 and DyxY1−x(OH)3, are potentially interesting systems to study transverse-field-induced quantum fluctuations effects in hard axis (Ising-type) magnetic materials.

Induced random fields in the LiHoxY1-xF4 quantum ising magnet in a transverse magnetic field

The LiHoxY1-xF4 magnetic material in a transverse magnetic field Bxx̂ perpendicular to the Ising spin direction has long been used to study tunable quantum phase transitions in a random disordered system. We show that the Bx-induced magnetization along the x̂ direction, combined with the local random dilution-induced destruction of crystalline symmetries, generates, via the predominant dipolar interactions between Ho3+ ions, random fields along the Ising ẑ direction. This identifies LiHoxY1-xF4 in Bx as a new random field Ising system. The random fields explain the rapid decrease of the critical temperature in the diluted ferromagnetic regime and the smearing of the nonlinear susceptibility at the spin-glass transition with increasing Bx and render the Bx-induced quantum criticality in LiHoxY1-xF4 likely inaccessible.

Switching on a two-dimensional time-harmonic scalar wave in the presence of a diffracting edge

This paper concerns the switching on of two-dimensional time-harmonic scalar waves. We first review the switch-on problem for a point source in free space, then proceed to analyse the analogous problem for the diffraction of a plane wave by a half-line (the ‘Sommerfeld problem’), determining in both cases the conditions under which the field is well-approximated by the solution of the corresponding frequency domain problem. In both cases the rate of convergence to the frequency domain solution is found to be dependent on the strength of the singularity on the leading wavefront. In the case of plane wave diffraction at grazing incidence the frequency domain solution is immediately attained along the shadow boundary after the arrival of the leading wavefront. The case of non-grazing incidence is also considered.

Approximation by harmonic polynomials in star-shaped domains and exponential convergence of Trefftz hp-dGFEM

We study the approximation of harmonic functions by means of harmonic polynomials in two-dimensional, bounded, star-shaped domains. Assuming that the functions possess analytic extensions to a delta-neighbourhood of the domain, we prove exponential convergence of the approximation error with respect to the degree of the approximating harmonic polynomial. All the constants appearing in the bounds are explicit and depend only on the shape-regularity of the domain and on delta. We apply the obtained estimates to show exponential convergence with rate O(exp(−b square root N)), N being the number of degrees of freedom and b>0, of a hp-dGFEM discretisation of the Laplace equation based on piecewise harmonic polynomials. This result is an improvement over the classical rate O(exp(−b cubic root N )), and is due to the use of harmonic polynomial spaces, as opposed to complete polynomial spaces.