C-H···O H-bonded complexes: how does basis set superposition error change their potential-energy surfaces?

Geometries, vibrational frequencies, and interaction energies of the CNH⋯O3 and HCCH⋯O3 complexes are calculated in a counterpoise-corrected (CP-corrected) potential-energy surface (PES) that corrects for the basis set superposition error (BSSE). Ab initio calculations are performed at the Hartree-Fock (HF) and second-order Møller-Plesset (MP2) levels, using the 6-31G(d,p) and D95++(d,p) basis sets. Interaction energies are presented including corrections for zero-point vibrational energy (ZPVE) and thermal correction to enthalpy at 298 K. The CP-corrected and conventional PES are compared; the unconnected PES obtained using the larger basis set including diffuse functions exhibits a double well shape, whereas use of the 6-31G(d,p) basis set leads to a flat single-well profile. The CP-corrected PES has always a multiple-well shape. In particular, it is shown that the CP-corrected PES using the smaller basis set is qualitatively analogous to that obtained with the larger basis sets, so the CP method becomes useful to correctly describe large systems, where the use of small basis sets may be necessary

Double proton transfer and one-electron oxidation behaviors in double H-bonded glycinamide-formamidine complex and comparison with biological base pair

The behaviors of double proton transfer (DPT) occurring in a representative glycinamide-formamidine complex have been investigated employing the B3LYP/6-311++G** level of theory. Computational results suggest that the participation of a formamidine molecule favors the proceeding of the proton transfer (PT) for glycinamide compared with that without mediator-assisted case. The DPT process proceeds with a concerted mechanism rather than a stepwise one since no zwitterionic complexes have been located during the DPT process. The barrier heights are 14.4 and 3.9 kcal/mol for the forward and reverse directions, respectively. However, both of them have been reduced by 3.1 and 2.9 kcal/mol to 11.3 and 1.0 kcal/mol with further inclusion of zero-point vibrational energy (ZPVE) corrections, where the lower reverse barrier height implies that the reverse reaction should proceed easily at any temperature of biological importance. Additionally, the one-electron oxidation process for the double H-bonded glycinamide-formamidine complex has also been investigated. The oxidated product is characterized by a distonic radical cation due to the fact that one-electron oxidation takes place on glycinamide fragment and a proton has been transferred from glycinamide to formamidine fragment spontaneously. As a result, the vertical and adiabatic ionization potentials for the neutral double H-bonded complex have been determined to be about 8.46 and 7.73 eV, respectively, where both of them have been reduced by about 0.79 and 0.87 eV relative to those of isolated glycinamide due to the formation of the intermolecular H-bond with formamidine. Finally, the differences between model system and adenine-thymine base pair have been discussed briefly.

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

The Role of Anharmonicity in Hydrogen-Bonded Systems: The Case of Water Clusters

The nature of vibrational anharmonicity has been examined for the case of small water clusters using second-order vibrational perturbation theory (VPT2) applied on second-order Møller–Plesset perturbation theory (MP2) potential energy surfaces. Using a training set of 16 water clusters (H2O)n=2–6,8,9 with a total of 723 vibrational modes, we determined scaling factors that map the harmonic frequencies onto anharmonic ones. The intermolecular modes were found to be substantially more anharmonic than intramolecular bending and stretching modes. Due to the varying levels of anharmonicity of the intermolecular and intramolecular modes, different frequency scaling factors for each region were necessary to achieve the highest accuracy. Furthermore, new scaling factors for zero-point vibrational energies (ZPVE) and vibrational corrections to the enthalpy (ΔHvib) and the entropy (Svib) have been determined. All the scaling factors reported in this study are different from previous works in that they are intended for hydrogen-bonded systems, while others were built using experimental frequencies of covalently bonded systems. An application of our scaling factors to the vibrational frequencies of water dimer and thermodynamic functions of 11 larger water clusters highlights the importance of anharmonic effects in hydrogen-bonded systems.

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.

Vibrational energy levels of water

Several quartic force fields and a full sextic anharmonic force field for H,O have been determined from high-quality ab initio calculations, the highest at the aug-cc-pVQZ CCSD(T) level of theory. These force fields have been used to determine vibrational excited state band origins up to 15 000 cm - ’ above the zero-point level, using both a perturbation-resonancea pproach and a variational approach. An optimisedq uartic force field hasb eeno btained by least squares refinement of our best ab initio results to fit the observed overtone levels of 5 symmetrically substituted isotopomers of water (Hi60, Hi70, HisO, D,O, and T,O) with an rms error of less than 10 cm-‘, using the perturbation-resonancem odel for the vibrational calculation. Predicatel east squaresr efinement was usedt o provide a loose constraint of the refined force field to the ab initio results. The results obtained prove the viability of the perturbation-resonancem odel for usei n larger molecular systemsa nd also highlight someo f its weaknesse

When the Casimir energy is not a sum of zero-point energies

We compute the leading radiative correction to the Casimir force between two parallel plates in the lambdaPhi(4) theory. Dirichlet and periodic boundary conditions are considered. A heuristic approach, in which the Casimir energy is computed as the sum of one-loop corrected zero-point energies, is shown to yield incorrect results, but we show how to amend it. The technique is then used in the case of periodic boundary conditions to construct a perturbative expansion which is free of infrared singularities in the massless limit. In this case we also compute the next-to-leading order radiative correction, which turns out to be proportional to lambda(3/2).

Factors controlling the competition among rotational and vibrational energy transfer channels in glyoxal

The state-to-state transfer of rotational and vibrational energy has been studied for S1 glyoxal (CHOCHO) in collisions with D2, N2, CO and C2H4 using crossed molecular beams. A laser is used to pump glyoxal seeded in He to its S1 zero point level with zero angular momentum about its top axis (K′ = 0). The inelastic scattering to each of at least 26 S1 glyoxal rotational and rovibrational levels is monitored by dispersed S1–S0 fluorescence. Various collision partners are chosen to investigate the relative influences of reduced mass and the collision pair interaction potential on the competition among the energy transfer channels. When the data are combined with that obtained previously from other collision partners whose masses range from 2 to 84 amu, it is seen that the channel competition is controlled primarily by the kinematics of the collisional interaction. Variations in the intermolecular potential play strictly a secondary role.

Non-Fermi liquid behavior in a mean-field study of the t-J model: Role of zero-point spin fluctuations

We present analytic results to show that the Schwinger-boson hole-fermion mean-field state exhibits non-Fermi liquid behavior due to spin-charge separation. The physical electron Green's function consists of three additive components. (a) A Fermi-liquid component associated with the bose condensate. (b) A non-Fermi liquid component which has a logarithmic peak and a long tail that gives rise to a linear density of states that is symmetric about the Fermi level and a momentum distribution function with a logarithmic discontinuity at the Fermi surface. (c) A second non-Fermi liquid component associated with the thermal bosons which leads to a constant density of states. It is shown that zero-point fluctuations associated with the spin-degrees of freedom are responsible for the logarithmic instabilities and the restoration of particle-hole symmetry close to the Fermi surface.

Zero-point entropy in 'spin ice'

Common water ice (ice I-h) is an unusual solid-the oxygen atoms form a periodic structure but the hydrogen atoms are highly disordered due to there being two inequivalent O-H bond lengths'. Pauling showed that the presence of these two bond lengths leads to a macroscopic degeneracy of possible ground states(2,3), such that the system has finite entropy as the temperature tends towards zero. The dynamics associated with this degeneracy are experimentally inaccessible, however, as ice melts and the hydrogen dynamics cannot be studied independently of oxygen motion(4). An analogous system(5) in which this degeneracy can be studied is a magnet with the pyrochlore structure-termed 'spin ice'-where spin orientation plays a similar role to that of the hydrogen position in ice I-h. Here we present specific-heat data for one such system, Dy2Ti2O7, from which we infer a total spin entropy of 0.67Rln2. This is similar to the value, 0.71Rln2, determined for ice I-h, SO confirming the validity of the correspondence. We also find, through application of a magnetic field, behaviour not accessible in water ice-restoration of much of the ground-state entropy and new transitions involving transverse spin degrees of freedom.

zero point of charge of cupric hydroxide and surface area determination by dye adsorption

zero point of charge of freshly precipitated cu(oh)2 has been determined to lie at pH 7.7 by means of microclectrophoresis technique. Day aged hydroxide shows an acid zpc shift to pH 7.3. these experimental values approximate the equivalence points of cu+ and oh_ ,which can be estimated from the solubility diagram constructed fo gu(oh)2 and cuo.

Comparison of salt titration and potentiometric titration methods for the determination of zero point of charge

the salt ritration metod was evaluated as a method to determine zpc in comparison with the potentiometric titration method for 26 soil with variable charge clays,i.e.,Oxisols and Ultisols from Thailand and Andisols from Japan. In addition to the determination of ST-pH0 as the zero point of charge, a calculation procedure was adopted here in order to acquire more information from the titration curve . fuithermore, for the purpose of cross-checking of zpc determined by the pt method, the st procedure was successively applied to the samples analyzed by the pt method.