Cyclic and ladder hydrogen bonded cyanamide oligomers: a density functional theory and many-body analysis approach

This work reports hydrogen bonding interaction in cyclic and ladder oligomers using density functional theory method. Many-body analysis technique has been used to study the nature of interactions between different molecules and their contribution to the binding energy of a respective hydrogen bonded oligomers. Hydrogen bonds in cyclic trimer to pentamer are stronger than those in corresponding ladder structures. Cyanamide monomer shows the lowest energy at B3LYP/aug-cc-pvdz level among different methods used here with the same basis set. The geometrical parameters for cyanamide monomer obtained at B3LYP/aug-cc-pvdz level are in excellent agreement with the experimental determinations. Cyclic structures are more stable than the ladder. In cyclic oligomers not only total two-body energies, but higher body energies also contribute significantly to the binding energy of a respective complex whereas in ladder, only total two-body energies contribute significantly and higher-body energies are almost negligible for cyanamide trimer to pentamer.     

New ab initio potential energy surface and the vibration-rotation-tunneling levels of (H2O)2 and (D2O)2

We report a new full-dimensional potential energy surface (PES) for the water dimer, based on fitting energies at roughly 30,000 configurations obtained with the coupled-cluster single and double, and perturbative treatment of triple excitations method using an augmented, correlation consistent, polarized triple zeta basis set. A global dipole moment surface based on Moller-Plesset perturbation theory results at these configurations is also reported. The PES is used in rigorous quantum calculations of intermolecular vibrational frequencies, tunneling splittings, and rotational constants for (H2O)2 and (D2O)2, using the rigid monomer approximation. Agreement with experiment is excellent and is at the highest level reported to date. The validity of this approximation is examined by comparing tunneling barriers within that model with those from fully relaxed calculations.     

Generalizing Rosenbluth's Algorithm to Include Along‐the‐Chain Intramolecular Energies

We incorporate the Boltzmann factors for inter‐monomer bending energy into the monomer growth direction choice in Rosenbluth's algorithm to model chains of arbitrary nearest‐neighbor rigidity. This allows for the consideration of compact (bent state lower in energy), free (straight and bent state equal in energy), or extended chains (bent state higher). We validate against, and compare to, various other results, showing very good agreement with known results for short chains and demonstrate the ability to model chains up to 500 segments long, far beyond the length at which the normal Rosenbluth method becomes unstable for reasonable nonzero bending energies. This approach is easily generalizable both to other energies determinable during chain growth, for example, polymers composed of more than one type of monomer with differing monomer interaction energies, as well as to other chain production algorithms. © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 , 57, 1684–1691     

On the determination of monomer dissociation energies of small water clusters from photoionization experiments

Recently, monomer dissociation energies of neutral water clusters were estimated via a thermodynamic cycle that utilized the measured appearance energies of vacuum ultraviolet (VUV) photoionized water clusters and the previously reported dissociation energies of protonated water clusters. The thermodynamic cycle incorrectly assumed that the energy difference between the (H2O)n+ --> (H2O)n-1H+ + OH* asymptotes (the relaxation energy) was zero. We show that these relaxation energies are large and cannot be neglected in the analysis. Thus, the neutral water cluster monomer dissociation energies cannot be directly determined from the measured ionization potentials because they are themselves involved in the appropriate thermodynamic cycle.     

Potential energy surfaces for small alcohol dimers I: methanol and ethanol

Potential energy landscapes for homogeneous dimers of methanol and ethanol were calculated using counterpoise (CP) corrected energies at the MP26-311+G(2df,2pd) level. The landscapes were sampled at approximately 15 dimer separation distances for different relative monomer geometries, or routes, given in terms of a relative monomer yaw, pitch, and roll and the spherical angles between the monomer centers (taken as the C atom attached to the O). The 19 different routes studied for methanol and the 22 routes examined for ethanol include 607 CP corrected energies. Both landscapes can be adequately represented by site-site, pairwise-additive models, suitable for use in molecular simulations. A modified Morse potential is used for the individual pair interactions either with or without point charges to represent the monomer charge distribution. A slightly better representation of the methanol landscape is obtained using point charges, while the potential energy landscape of ethanol is slightly better without point charges. This latter representation may be computationally advantageous for molecular simulations because it avoids difficulties associated with long-range effects of point-charge-type models.     

Elastic/inelastic and charge transfer collisions of H++H2 at collision energies of 4.67, 6, 7.3, and 10 eV

Quantum mechanical studies of vibrational and rotational state-resolved differential cross sections, integral cross sections, and transition probabilities for both the elastic/inelastic and charge transfer processes have been carried out at collision energies of 4.67, 6, 7.3, and 10 eV using the vibrational close-coupling rotational infinite-order sudden approach. The dynamics has been performed employing our newly obtained quasidiabatic potential energy surfaces which were generated using ab initio procedures and Dunning's correlation-consistent-polarized quadrupole zeta basis set. The present theoretical results for elastic/inelastic processes provide an overall excellent agreement with the available experimental data and they are also found to be almost similar to that obtained in earlier theoretical results using the ground electronic potential energy surface, lending credence to the accuracy and reliability of the quasidiabatic potential energy surfaces. The results for the complementary charge transfer processes are also presented at these energies.     

Saturated hydraulic conductivity of a volcanic ash soil affected by repulsive potential energy in a multivalent anionic system

Acid rain is supposed to influence soil structures, because soils have pH-dependent charges. The adverse effects of acid rain on soils must be assessed. Although repulsive potential energy among soil clay particles generates swelling and dispersion, thereby changing the soil’s hydraulic conductivity, the relationship between hydraulic conductivity and repulsive potential energy has not been evaluated. Moreover, research into repulsive potential energy in multivalent counterionic systems has been rare. In this paper, repulsive potential energies for a volcanic ash soil (allophanic Andisol), which is characterized by a number of pH-dependent charges, were evaluated in a multivalent counterionic system. Changes in saturated hydraulic conductivity (K) of volcanic ash soil during dilute acid leaching and their relationship with the repulsive potential energies were examined. When nitric acid at pH 3 or 4 was leached, K decreased rapidly. On the other hand, the decrease in K attenuated as the proportion of sulfate in the dilute acid increased. Electrophoretic mobilities were measured and the zeta potentials were estimated. From the zeta potentials and the calculation of repulsive potential energies between the clay particles in the NO₃–SO₄ system, we concluded that the decrease in K for an acid solution with a high proportion of nitrate was due to swelling and dispersion of the soil induced by electrostatic repulsive potential energy. Because sulfate formed complexes on the clay surface, the repulsive potential energy decreased as the proportion of sulfate in the dilute acid increased. Then, the flocculation of the soil was maintained, thereby inhibiting the decrease in K.     

6-甲基-4-羟基嘧啶单体及二聚体质子转移过程的理论研究

应用密度泛函理论的B3LYP/6-311+G(d)方法研究了6-甲基-4-羟基嘧啶单体及二聚体质子转移的异构化反应.对反应势能面的研究发现,该化含物可能存在9种单体异构体,对其最稳定的单体构型进行分析.各单体间异构化反应的过渡态共有9种,反应的活化能最小为22.06 kJ/mol,最大为356.55 kJ/mol,最可能的反应路径在室温下即可进行. 研究了2种二聚体及其异构化反应的过渡态,发现二聚体均比其对应的单体稳定,而且质子转移所需要的活化能仅为20.13 kJ/mol,比单体低很多. 氢键在这种变化中起了主要作用,由单体和二聚体的总能量计算了氢键的键能.     

Cation spectroscopy and binding energy determination for 1,4-benzodioxan-Ar1 and -Ar2 complexes

Cation vibronic spectra are measured for 1,4-benzodioxan (BZD) and van der Waals complexes of BZD with one and two Ar atoms using zero electron kinetic energy and mass analyzed threshold ionization spectroscopy. The spectra of the monomer cation were used to measure the frequencies of the two key low-frequency modes which had previously been extensively studied in the neutral S0 and S1 states. The aliphatic ring twisting mode, nu25, has an energy of 146 cm(-1) in the cation, intermediate between the values found in the S0 and S1 states. The bending, butterfly-like mode nu48 has an energy of 125 cm(-1), which is of higher frequency than either of the neutral states. The S1 spectra of the BZD-Ar1 and BZD-Ar2 complexes are recorded and observed to have modest red shifts from the monomer. The cation spectra of the complexes are also measured using mass analyzed threshold ionization spectroscopy including scans at higher energy which are used to determine the Ar binding energies. The energies for the loss of one Ar atom were determined to be 630 +/- 10 and 650 +/- 10 cm(-1) for BZD-Ar and BZD-Ar2, respectively. The similar cation spectra and similar binding energies indicate that each Ar atom in BZD-Ar2 has a similar binding geometry. Quantum chemical calculations were performed which had fair agreement with the measured binding energies and give some insight into the specific binding geometry.     

Development of a Flexible‐Monomer Two‐Body Carbon Dioxide Potential and Its Application to Clusters up to (CO2)13

A flexible‐monomer two‐body potential energy function was developed that approaches the high level CCSD(T)/CBS potential energy surface (PES) of carbon dioxide (CO₂) systems. This function was generated by fitting the electronic energies of unique CO₂ monomers and dimers to permutationally invariant polynomials. More than 200,000 CO₂ configurations were used to train the potential function. Comparisons of the PESs of six orientations of flexible CO₂ dimers were evaluated to demonstrate the accuracy of the potential. Furthermore, the potential function was used to determine the minimum energy structures of CO₂ clusters containing as many as 13 molecules. For isomers of (CO₂)₃, the potential demonstrated energetic agreement with the M06‐2X functional and structural agreement of the B2PLYP‐D functional at substantially reduced computational costs. A separate function, fit to MP2/aug‐cc‐pVDZ reference energies, was developed to directly compare the two‐body potential to the ab initio MP2 level of theory. © 2017 Wiley Periodicals, Inc.     

Many-Body Rayleigh-Schrödinger Perturbation calculations of the correlation energy of open shell molecules in the restricted Roothaan-Hartree-Fock formalism. Application to heats of reaction and energies of activation

The title technique was applied to a series of elementary chemical reactions. Second and third order contributions to the correlation energy were computed for the basis sets of the double zeta and double zeta plus polarization quality. Calculated heats of reaction and energies of activation were compared with the experimental data and the results of the bestab initio calculations reported in the literature.     

Potential energy surfaces for small alcohol dimers. II. Propanol, isopropanol, t-butanol, and sec-butanol

Potential energy landscapes for homogeneous dimers of propanol, isopropanol, tert-butanol, and sec-butanol were obtained using 735 counterpoise-corrected energies at the MP2/6-311+G(2df,2pd) level. The landscapes were sampled at 15 dimer separation distances for different relative monomer geometries, or routes, given in terms of the yaw, pitch, and roll of one monomer relative to the other and the spherical angles between the two monomer centers (taken as the C atom attached to the O). The resultant individual energy surfaces and their complex topographies were also regressed using a site-site pair potential model using a modified Morse potential that provides a mathematically simple representation of the landscapes suitable for use in molecular simulations. Generalized Morse parameters were also obtained for this model from a composite regression of these energy landscapes and those previously reported for methanol and ethanol. The quality of fit for all these energy landscapes suggests that these site parameters have transferability for possible use on other alcohols.     

Calculation of the O-H stretching vibrational overtone spectrum of the water dimer

The O-H stretching vibrational overtone spectrum of the water dimer has been calculated with the dimer modeled as two individually vibrating monomer units. Vibrational term values and absorption intensities have been obtained variationally with a computed dipole moment surface and an internal coordinate Hamiltonian, which consists of exact kinetic energy operators within the Born-Oppenheimer approximation of the monomer units. Three-dimensional ab initio potential energy and dipole moment surfaces have been calculated using the internal coordinates of the monomer units using the coupled cluster method including single, double, and perturbative triple excitations [CCSD(T)] with the augmented correlation consistent valence triple zeta basis set (aug-cc-pVTZ). The augmented correlation consistent valence quadruple zeta basis set (aug-cc-pVQZ), counterpoise correction, basis set extrapolation to the complete basis set limit, relativistic corrections, and core and valence electron correlations effects have been included in one-dimensional potential energy surface cuts. The aim is both to investigate the level of ab initio and vibrational calculations necessary to produce accurate results when compared with experiment and to aid the detection of the water dimer under atmospheric conditions.     

Ab initio study of chiral recognition in the propylene imine.hydrogen peroxide complex

High level ab initio calculations were employed to study the chiral recognition effect in several chiral molecular pairs that consist of the propylene imine and hydrogen peroxide molecules. The potential energy surfaces for the complexes formed between S-cis-1,2-propylene imine and the two enantiomeric forms of hydrogen peroxide were constructed, using the calculated interaction energies at different separations and orientations. The energy calculations were done using the MOLPRO suite of programs with CCSD(T)/cc-pVDZ. The energies were counterpoise corrected at every point to eliminate the basis set superposition error. Complete geometry optimizations were further carried out for the molecular complexes consisting of the cis- or trans-propylene imine isomers and the two enantiomeric forms of hydrogen peroxide. The geometry optimizations were done using the Gaussian 98 and 03 suites of programs, with MP2/aug-cc-pVDZ being the highest level used. Altogether, eight stable complexes were identified, and the corresponding dissociation energies were calculated with MP2/aug-cc-pVTZ. The largest chirodiastaltic energy is found at 0.74 kcal mol(-1) for the (syn)trans-propylene imine.hydrogen peroxide complexes, where hydrogen peroxide acts as a hydrogen donor and is on the opposite side of the ring from the methyl group. The rotational constants, dipole moments, and harmonic frequencies of the complexes are presented to assist future spectroscopic investigations.     

A π-stacked phenylacetylene dimer

The structure of the phenylacetylene-dimer has been elucidated using IR-UV double resonance spectroscopy in combination with high level ab initio calculations at the CCSD(T)/CBS level. The IR spectra in the acetylenic and the aromatic C-H stretching regions indicate that the two phenylacetylene moieties are in identical environments and very similar to the phenylacetylene monomer. Calculated stabilization energies and the free energies at the CCSD(T)/CBS level favor the formation of an anti-parallel π-stacked structure. The DFT-SAPT energy decomposition analysis points out that the anti-parallel π-stacked structure maximizes electrostatic as well as the dispersion components of energy. The observed IR spectra are consistent with the anti-parallel π-stacked structure.     

Gaussian multipoles in practice: Electrostatic energies for intermolecular potentials

A method is presented for calculating the total electrostatic interaction energies between molecules from ab initio monomer wave functions. This approach differs from existing methods, such as Stone's distributed multipole analysis (DMA), in including the short-range penetration energy as well as the long-range multipolar energy. The monomer charge densities are expressed as distributed series of atom-centered functions which we call Gaussian multipoles; these are analogous to the distributed point multipoles used in DMA. Our procedure has been encoded in the GMUL program. Calculations have been performed on the formamide/formaldehyde complex, a model system for N? H …? O hydrogen bonding in biological molecules, and also on guanidinium/benzene, modeling amino/aromatic interactions in proteins. We find that the penetration energy can be significant, especially in its contribution to the variation of the electrostatic energy with interaction geometry. A hybrid method, which uses Gaussian multipoles for short-range atom pair interactions and point multipoles for long-range ones, allows the electrostatic energies, including penetration, to be calculated at a much reduced cost. We also note that the penetration energy may provide the best route to an atom–atom anisotropic model for the exchange-repulsion energy in intermolecular potentials. © 1994 by John Wiley & Sons, Inc.     

Photodissociation of small group-11 metal cluster ions: fragmentation pathways and photoabsorption cross sections

Noble metal cluster ions Cu(n)(+), Ag(n)(+) and Au(n)(+) (n = 3-21) have been stored in a Penning trap and photodissociated by low intensity laser pulses of 10 ns at photon energies of 3.49 eV and 4.66 eV. The fragmentation pathways, neutral monomer and dimer evaporation, have been monitored as a function of cluster size, excitation energy and element. It is found that the behavior of the branching ratio between monomer and dimer evaporation as a function of excitation energy depends on the metal under investigation. In particular, the slope of the energy dependence is positive for silver but negative for gold and copper cluster ions. Furthermore, photoabsorption cross sections are determined from observed total fragment yields in single-photon dissociation.     

Intermolecular potential calculations for polynuclear aromatic hydrocarbon clusters

Calculations of intermolecular potentials are presented for homo-molecular and hetero-molecular clusters of 24 peri-condensed PAH spanning monomer masses ranging from 78 to 1830 Da. Binding energies of homo-molecular dimers rise rapidly with molecular size and asymptotically approach the experimentally established exfoliation energy for graphite of 5.0 kJ mol(-1) (carbon atom)(-1). Binding energies of hetero-molecular dimers correlate well with the reduced mass of the pair. From calculations of homo-molecular stacks, binding energies were observed to increase with each added molecule and rise asymptotically, approaching a limit which scales linearly with monomer molecular mass. These results are reviewed in the context of molecular growth in flames and in the context of astrophysical observations.     

Quantum and classical studies of vibrational motion of CH5+ on a global potential energy surface obtained from a novel ab initio direct dynamics approach

We report a full dimensional, ab initio based potential energy surface for CH(5) (+). The ab initio electronic energies and gradients are obtained in direct-dynamics calculations using second-order M?ller-Plesset perturbation theory with the correlation consistent polarized valence triple zeta basis. The potential energy and the dipole moment surfaces are fit using novel procedures that ensure the full permutational symmetry of the system. The fitted potential energy surface is tested by comparing it against additional electronic energy calculations and by comparing normal mode frequencies at the three lowest-lying stationary points obtained from the fit against ab initio ones. Well-converged diffusion Monte Carlo zero-point energies, rotational constants, and projections along the CH and HH bond lengths and the tunneling coordinates are presented and compared with the corresponding harmonic oscillator and standard classical molecular dynamics ones. The delocalization of the wave function is analyzed through comparison of the CH(5) (+) distributions with those obtained when all of the hydrogen atoms are replaced by (2)H and (3)H. The classical dipole correlation function is examined as a function of the total energy. This provides a further probe of the delocalization of CH(5) (+).     

Quantum dynamics of vibrationally activated OH-CO reactant complexes

The MP2 complete basis set (CBS) limit for the binding energy of the two low-lying water octamer isomers of D2d and S4 symmetry is estimated at -72.7+/-0.4 kcal/mol using the family of augmented correlation-consistent orbital basis sets of double through quintuple zeta quality. The largest MP2 calculation with the augmented quintuple zeta (aug-cc-pV5Z) basis set produced binding energies of -73.70 (D2d) and -73.67 kcal/mol (S4). The effects of higher correlation, computed at the CCSD(T) level of theory, are estimated at <0.1 kcal/mol. The newly established MP2/CBS limit for the water octamer is reproduced quite accurately by the newly developed all atom polarizable, flexible interaction potential (TTM2-F). The TTM2-F binding energies of -73.21 (D2d) and -73.24 kcal/mol (S4) for the two isomers are just 0.5 kcal/mol (or 0.7%) larger than the MP2/CBS limit.