The carbon dioxide molecule

On the performance of composite schemes in determining equilibrium molecular structures

Determination of equilibrium molecular structures is an essential ingredient in predicting spectroscopic parameters that help in identifying molecular carriers of microwave transitions. Here, the performance of two different ab initio composite approaches for obtaining equilibrium structures, "energy scheme" and "geometry scheme," is explored and compared to semi-experimental equilibrium structures. This study is performed for a set of 11 molecules which includes diatomics, linear triatomics, and a few non-linear molecules. The ab initio calculations were performed using three tiers of composite chemical recipes. The current results show that as the overall rigor of calculation is increased, the semi-experimental and the ab initio numbers agree to within 0.0003 Å for all molecules in the test set. The composite approach based on correcting the potential energy surface (energy scheme) and the one based on correcting the geometry directly (geometry scheme) show excellent agreement with each other. This work represents a step toward development of efficient and highly accurate procedures for computing ab initio equilibrium structures.

Rotational excitation of CO2 induced by He: New potential energy surface and scattering calculations

The CO2 molecule is of great interest for astrophysical studies since it can be found in a large variety of astrophysical media where it interacts with the dominant neutral species, such as He, H2, or H2O. The CO2–He collisional system was intensively studied over the last two decades. However, collisional data appear to be very sensitive to the potential energy surface (PES) quality. Thus, we provide, in this study, a new PES of the CO2–He van der Waals complex calculated with the coupled-cluster method and a complete basis set extrapolation in order to provide rotational rate coefficients that are as accurate as possible. The PES accuracy was tested through the calculations of bound state transition frequencies and pressure broadening coefficients that were compared to experimental data. An excellent agreement was globally found. Then, revised collisional data were provided for the 10–300 K temperature range. Rate coefficients were compared to previously computed ones and are found to be up to 50% greater than previously provided ones. These differences can induce non-negligible consequences for the modeling of CO2 abundance in astrophysical media.

Analytic energy gradients for the orbital-optimized third-order Møller-Plesset perturbation theory

Analytic energy gradients for the orbital-optimized third-order Møller-Plesset perturbation theory (OMP3) [U. Bozkaya, J. Chem. Phys. 135, 224103 (2011)] are presented. The OMP3 method is applied to problematic chemical systems with challenging electronic structures. The performance of the OMP3 method is compared with those of canonical second-order Møller-Plesset perturbation theory (MP2), third-order Møller-Plesset perturbation theory (MP3), coupled-cluster singles and doubles (CCSD), and coupled-cluster singles and doubles with perturbative triples [CCSD(T)] for investigating equilibrium geometries, vibrational frequencies, and open-shell reaction energies. For bond lengths, the performance of OMP3 is in between those of MP3 and CCSD. For harmonic vibrational frequencies, the OMP3 method significantly eliminates the singularities arising from the abnormal response contributions observed for MP3 in case of symmetry-breaking problems, and provides noticeably improved vibrational frequencies for open-shell molecules. For open-shell reaction energies, OMP3 exhibits a better performance than MP3 and CCSD as in case of barrier heights and radical stabilization energies. As discussed in previous studies, the OMP3 method is several times faster than CCSD in energy computations. Further, in analytic gradient computations for the CCSD method one needs to solve λ-amplitude equations, however for OMP3 one does not since λ(ab)(ij(1))=t(ij)(ab(1)) and λ(ab)(ij(2))=t(ij)(ab(2)). Additionally, one needs to solve orbital Z-vector equations for CCSD, but for OMP3 orbital response contributions are zero owing to the stationary property of OMP3. Overall, for analytic gradient computations the OMP3 method is several times less expensive than CCSD (roughly ~4-6 times). Considering the balance of computational cost and accuracy we conclude that the OMP3 method emerges as a very useful tool for the study of electronically challenging chemical systems.

Orbital-free embedding applied to the calculation of induced dipole moments in CO2⋯X (X=He, Ne, Ar, Kr, Xe, Hg) van der Waals complexes

The orbital-free frozen-density embedding scheme within density-functional theory [T. A. Wesolowski and A. Warshel, J. Phys. Chem. 97, 8050 (1993)] is applied to the calculation of induced dipole moments of the van der Waals complexes CO2...X (X = He, Ne, Ar, Kr, Xe, Hg). The accuracy of the embedding scheme is investigated by comparing to the results of supermolecule Kohn-Sham density-functional theory calculations. The influence of the basis set and the consequences of using orbital-dependent approximations to the exchange-correlation potential in embedding calculations are examined. It is found that in supermolecular Kohn-Sham density-functional calculations, different common approximations to the exchange-correlation potential are not able to describe the induced dipole moments correctly and the reasons for this failure are analyzed. It is shown that the orbital-free embedding scheme is a useful tool for applying different approximations to the exchange-correlation potential in different subsystems and that a physically guided choice of approximations for the different subsystems improves the calculated dipole moments significantly.

Least-Squares Mass-Dependence Molecular Structures

The major part of the difference between the zero-point and equilibrium moment of inertia of a molecule is a homogeneous function of degree 12 in the atomic masses. The use of explicit mathematical models of is studied here as a means of determining near-equilibrium molecular structures from the zero-point moments of inertia of a range of isotopomers. It is found that models with more than two parameters per axis generally give strongly correlated fits for the number of isotopomers typically available. The two-parameter model, where N is the number of atoms, mi are the atomic masses, and M is the molecular mass, usually gives well-conditioned fits with standard deviations in the parts-per-million range, and structures called that are often close to equilibrium structures. Laurie-type corrections for hydrogen atoms and parameters to describe isotopic rotations of the epsilon tensor can also be included. Copyright 1999 Academic Press.