Accurate ab initio potential energy surfaces for the (3)A '' and (3)A ' electronic states of the O(P-3) plus HBr system

In this work, we report the construction of potential energy surfaces for the (3)A '' and (3)A' states of the system O(P-3) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O(P-3) + HBr -> OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A '' electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A' surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A '' and 4.16 kcal/mol for the (3)A' state. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4705428]

The infrared fundamental intensities of some cyanopolyynes

Some cyanopolyynes, HCnN (n = 1, 3, ... , 17), are investigated by means of calculations at the MP2/cc-pVTZ and CCSD/cc-pVDZ levels. Although the MP2/cc-pVTZ results for geometries and molecular dipole moments are encouraging, the CCSD/cc-pVDZ level was superior for the study of infrared fundamental intensities. The main bands are also analyzed with a charge-charge flux-dipole flux (CCFDF) partition model based on quantities given by the Quantum Theory of Atoms in Molecules (QTAIM). The intensity of vibrations corresponding to the stretching of CH bonds (3471-3473 cm(-1)) increases in line with the number of carbon atoms (from 61 to 146 km mol(-1) between HCN and HC13N). This increase is due to the charge flux contribution while the other contributions remain roughly unaltered except for HCN. Moreover, the hydrogen atom loses an almost constant amount of electronic charge during the CH bond enlargement and a small fraction of this charge spreads to atoms farther and farther away from hydrogen as the molecule size increases. The band associated with the doubly degenerate CH bending vibrations (643-732 cm(-1)) presents approximately the same intensity in all the studied cyanopolyynes (from 67 to 76 km mol(-1)). The CCFDF/QTAIM contributions are also nearly the same for these bending modes in HC5N and larger systems. The intensity of the mode mostly identified as CN stretching (around 2378-2399 cm(-1) except for HCN) increases from HCN up to HC7N (from 0.3 to 83 km mol(-1)) and nearly stabilizes around 80-90 km mol(-1) for larger systems. The CCFDF/QTAIM contributions for this mode also change significantly up to HC7N and remain almost constant in larger systems. We also observed the appearing of a very relevant band between 2283 and 2342 cm(-1). This mode is mainly associated with the symmetric stretching of CC triple bonds near the molecule center and exhibits large charge fluxes while the other contributions are almost negligible in the largest cyanopolyynes. The two vibrational bands associated with the smallest frequencies are also studied and extrapolation equations are suggested to predict their positions in larger cyanopolyynes. (C) 2012 Elsevier B.V. All rights reserved.

The nuclear electric quadrupole moment of hafnium from the molecular method

The molecular method is used to obtain nuclear electric quadrupole moment (NQM) values for hafnium through electric field gradients (EFGs) at this nucleus in HfO and HfS. Dirac-Coulomb calculations with the Coupled Cluster approach, DC-CCSD (T) and DC-CCSD-T, were carried out to achieve the most accurate estimates of these EFGs. Higher order corrections are also added. Hence, the most reliable values for 177Hf and 179Hf determined here are 3319(33) and 3750(37) mbarn, respectively, in nice accordance with the best currently accepted NQMs for this element. (C) 2012 Elsevier B.V. All rights reserved.

Coupled Boltzmann calculation of mixed axion/neutralino cold dark matter production in the early universe

We calculate the relic abundance of mixed axion/neutralino cold dark matter which arises in R-parity conserving supersymmetric (SUSY) models wherein the strong CP problem is solved by the Peccei-Quinn (PQ) mechanism with a concommitant axion/saxion/axino supermultiplet. By numerically solving the coupled Boltzmann equations, we include the combined effects of 1. thermal axino production with cascade decays to a neutralino LSP, 2. thermal saxion production and production via coherent oscillations along with cascade decays and entropy injection, 3. thermal neutralino production and re-annihilation after both axino and saxion decays, 4. gravitino production and decay and 5. axion production both thermally and via oscillations. For SUSY models with too high a standard neutralino thermal abundance, we find the combined effect of SUSY PQ particles is not enough to lower the neutralino abundance down to its measured value, while at the same time respecting bounds on late-decaying neutral particles from BBN. However, models with a standard neutralino underabundance can now be allowed with either neutralino or axion domination of dark matter, and furthermore, these models can allow the PQ breaking scale f(a) to be pushed up into the 10(14) - 10(15) GeV range, which is where it is typically expected to be in string theory models.

Encapsulation of small magnetic clusters in fullerene cages: A density functional theory investigation within van der Waals corrections

The encapsulation of magnetic transition-metal (TM) clusters inside carbon cages (fullerenes, nanotubes) has been of great interest due to the wide range of applications, which spread from medical sensors in magnetic resonance imaging to photonic crystals. Several theoretical studies have been reported; however, our atomistic understanding of the physical properties of encapsulated magnetic TM 3d clusters is far from satisfactory. In this work, we will report general trends, derived from density functional theory within the generalized gradient approximation proposed by Perdew, Burke, and Ernzerhof (PBE), for the encapsulation properties of the TMm@C-n (TM = Fe, Co, Ni; m = 2-6, n = 60,70,80,90) systems. Furthermore, to understand the role of the van der Waals corrections to the physical properties, we employed the empirical Grimme's correction (PBE + D2). We found that both PBE and PBE + D2 functionals yield almost the same geometric parameters, magnetic and electronic properties, however, PBE + D2 strongly enhances the encapsulation energy. We found that the center of mass of the TMm clusters is displaced towards the inside C-n surfaces, except for large TMm clusters (m = 5 and 6). For few cases, e. g., Co-4 and Fe-4, the encapsulation changes the putative lowest-energy structure compared to the isolated TMm clusters. We identified few physical parameters that play an important role in the sign and magnitude of the encapsulation energy, namely, cluster size, fullerene equatorial diameter, shape, curvature of the inside C-n surface, number of TM atoms that bind directly to the inside C-n surface, and the van der Waals correction. The total magnetic moment of encapsulated TMm clusters decreases compared with the isolated TMm clusters, which is expected due to the hybridization of the d-p states, and strongly depends on the size and shape of the fullerene cages.