Molecular dynamics simulations of binary structure H hydrogen and methyl-tert-butylether clathrate hydrates

Binary structure H (sH) hydrogen and methyl-tert-butylether (MTBE) clathrate hydrates are studied with molecular dynamics simulations. Simulations on a 3 x 3 x 3 sH unit cell with up to 4.7 mass % hydrogen gas are run at pressures of 100 bars and 2 kbars at 100 and 273 K. For the small and medium cages of the sH unit cell, H2 guest molecule occupancies of 0, 1 (single occupancy), and 2 (double occupancy) are considered with the MTBE molecule occupying all of the large cages. An increase of the small and medium cage occupancies from 1 to 2 leads to a jump in the unit cell volume and configurational energy. Calculations are also set up with 13, 23, and 89 of the MTBE molecules in the large cages replaced by sets of three to six H2 molecules, and the effects on the configurational energy and volume of the simulation cell are determined. As MTBE molecules are replaced with sets of H2 guests in the large cages, the configurational energy of the unit cell increases. At the lower temperature, the energy and volume of the clathrate are not sensitive to the number of hydrogen guests in the large cages; however, at higher temperatures the repulsions among the H2 guest molecules in the large cages cause an increase in the system energy and volume.     

Rapid Estimation of Basis Set Error and Correlation Energy Based on Mulliken Charges and Mulliken Matrix with the Small 6-31g* Basis Set

Good, density functional quality (B3LYP/6-31G*) ground state total electronic energies have been approximated using single point Hartree–Fock-self consistent field (HF-SCF/6-31G*) total energies and Mulliken partial charges versus. Mulliken matrix (electrons assigned to atoms and atoms pairs from Mulliken population analysis). This is a development of our rapid estimation of basis set error and correlation energy from partial charges (REBECEP) method, published earlier (see references [21,22,30]. The development is as follows: (1) A larger set of atoms (H, C, N, O, F, Si, P and S) are considered as building blocks for closed shell, neutral, ground state molecules at their equlibrium geometry; (2) geometries near equilibrium geometry are also considered; (3) A larger set, containing 115 molecules, was used to fit REBECEP parameters; (4) most importantly, electrons belonging to chemical bonds (between atom pairs) are also considered (Mulliken matrix) in addition to the atoms (Mulliken charges), using more REBECEP parameters to fit and yielding a more flexible algorithm. With these parameters a rather accurate closed shell ground state electronic total energy can be obtained from a small basis set HF-SCF calculation in the vicinity of optimal geometry. The 3.3 kcal/mol root mean square deviation of REBECEP improves to 1.5 kcal/mol when using Mulliken matrix instead of Mulliken charges.     

Multireference coupled-cluster theory: the easy way

The multi-ionization equation-of-motion coupled-cluster (CC) method is developed for multireference (MR) problems. It is operationally single reference, depending upon a formal matrix diagonalization step to define the coefficients in the wavefunction in an unbiased way that allows for important MR character. The method is illustrated for the autoisomerization of cyclobutadiene, which has a very large multireference effect and compared to other MR-CC results. The newly implemented methods are also used to obtain the vertical double ionization (DI) potentials of several small molecules (H(2)O, CO, C(2)H(2), C(2)H(4)). Also, the performance of the new methods is analyzed by plotting the potential energy curve for twisted ethylene as a function of a dihedral angle between two methylenes. Evaluation of the total molecular energy via MR-DI-CC calculations makes it possible to avoid an unphysical cusp.     

Comparison of localization and delocalization indices obtained with Hartree-Fock and conventional correlated methods: effect of Coulomb correlation

Atomic populations and localization [lambda(A)] and delocalization [delta(A,B)] indices (LIs and DIs) are calculated for a large set of molecules at the Hartree-Fock (HF), MP2, MP4(SDQ), CISD, and QCISD levels with the 6-311++G(2d,2p) basis set. The HF method and the conventional correlation methods [MP2, MP4(SDQ), CISD, and QCISD] yield distinct sets of LIs and DIs. Yet, within the four conventional correlation methods the differences in atomic populations and LIs and DIs are small. Relative to HF, the conventional correlation methods [MP2, MP4(SDQ), CISD, QCISD] yield virtually the same LIs and DIs for molecules with large charge separations while LIs and DIs that differ significantly from the HF values--the LIs are increased and DIs decreased--are obtained for bonds with no or small charge separations. Such is the case in the archetypal homopolar molecules HC(triple bond)CH, H2C=CH2, CH3-CH3, and "protonated cyclopropane" C(3)H(7) (+), in which case the bonding may be atypical. Relative to HF, the typical effect of the conventional correlation methods is to decrease the DI between atoms.     

Molecular-dynamics simulations of binary structure II hydrogen and tetrahydrofurane clathrates

The binary structure II hydrogen and tetrahydrofurane (THF) clathrates are studied with molecular-dynamics simulations. Simulations are done at pressures of 120 and 1.013 bars for temperatures ranging from 100 to 273 K. For the small cages of the structure II unit cell, H2 guest molecule occupancies of 0, 16 (single occupancy), and 32 (double occupancy) are considered. THF occupancies of 0-8 in the large cages are studied. For cases in which THF does not occupy all large cages in a unit cell, the remaining large cages can be occupied with sets of four H2 guest molecules. The unit-cell volumes and configurational energies are compared in the different occupancy cases. Increasing the small cage occupancy leads to an increase in the unit-cell volume and thermal-expansion coefficient. Among simulations with the same small cage occupancy, those with the large cages containing 4H2 guests have the largest volumes. The THF guest molecules have a stabilizing effect on the clathrate and the configurational energy of the unit cell decreases linearly as the THF content increases. For binary THF + H2 clathrates, the substitution of the THF molecules in the large cages with sets of 4H2 molecules increases the configurational energy. For the binary clathrates, various combinations of THF and H2 occupancies have similar configurational energies.     

Variational quantum mechanical and active database approaches to the rotational-vibrational spectroscopy of ketene, H2CCO

A variational quantum mechanical protocol is presented for the computation of rovibrational energy levels of semirigid molecules using discrete variable representation of the Eckart-Watson Hamiltonian, a complete, "exact" inclusion of the potential energy surface, and selection of a vibrational subspace. Molecular symmetry is exploited via a symmetry-adapted Lanczos algorithm. Besides symmetry labels, zeroth-order rigid-rotor and harmonic-oscillator quantum numbers are employed to characterize the computed rovibrational states. Using the computational molecular spectroscopy algorithm presented, a large number of rovibrational states, up to J = 50, of the ground electronic state of the parent isotopologue of ketene, H(2) (12)C=(12)C=(16)O, were computed and characterized. Based on 12 references, altogether 3982 measured and assigned rovibrational transitions of H(2) (12)C=(12)C=(16)O have been collected, from which 3194 were validated. These transitions form two spectroscopic networks (SN). The ortho and the para SNs contain 2489 and 705 validated transitions and 1251 and 471 validated energy levels, respectively. The computed energy levels are compared with energy levels obtained, up to J = 41, via an inversion protocol based on this collection of validated measured rovibrational transitions. The accurate inverted energy levels allow new assignments to be proposed. Some regularities and irregularities in the rovibrational spectrum of ketene are elucidated.     

Calculation of Verdet constants with time-dependent density functional theory: implementation and results for small molecules

We report the implementation of a method by which to calculate Verdet constants for molecules. The method is based on gauge-including atomic orbitals (GIAOs) and density functional theory. Calculations based on this method afford magneto-optical rotations of the right magnitude for the molecules H2, N2, CO, HF, CH4, C2H2, H2O, and CS2. The results are in satisfactory agreement with experiment. We investigate the dependency of the results on the gauge origin if GIAOs are not chosen, the convergence of the results with the size of the basis set for AOs and GIAOs, and for H2O and CS2 a comparison of gas-phase and liquid phase values. For the small molecules studied here, large polarized basis sets with diffuse functions are required to obtain well converged results. The use of an asymptotically correct Kohn-Sham potential is advantageous.     

Novel method for geometry optimization of molecular clusters: application to benzene clusters

A heuristic and unbiased method for searching optimal geometries of clusters of nonspherical molecules was constructed from the algorithm recently proposed for Lennard-Jones atomic clusters. In the method, global minima are searched by using three operators, interior, surface, and orientation operators. The first operator gives a perturbation on a cluster configuration by moving molecules near the center of mass of a cluster, and the second one modifies a cluster configuration by moving molecules to the most stable positions on the surface of a cluster. The moved molecules are selected by employing a contribution of the molecules to the potential energy of a cluster. The third operator randomly changes the orientations of all molecules. The proposed method was applied to benzene clusters. It was possible to find new global minima for (C6H6)11, (C6H6)14, and (C6H6)15. Global minima for (C6H6)16 to (C6H6)30 are first reported in this article.     

Encapsulation of carbon chain molecules in single-walled carbon nanotubes

The vacuum space inside carbon nanotubes offers interesting possibilities for the inclusion, transportation, and functionalization of foreign molecules. Using first-principles density functional calculations, we show that linear carbon-based chain molecules, namely, polyynes (C(m)H(2), m = 4, 6, 10) and the dehydrogenated forms C(10)H and C(10), as well as hexane (C(6)H(14)), can be spontaneously encapsulated in open-ended single-walled carbon nanotubes (SWNTs) with edges that have dangling bonds or that are terminated with hydrogen atoms, as if they were drawn into a vacuum cleaner. The energy gains when C(10)H(2), C(10)H, C(10), C(6)H(2), C(4)H(2), and C(6)H(14) are encapsulated inside a (10,0) zigzag-shaped SWNT are 1.48, 2.04, 2.18, 1.05, 0.55, and 1.48 eV, respectively. When these molecules come inside a much wider (10,10) armchair SWNT along the tube axis, they experience neither an energy gain nor an energy barrier. They experience an energy gain when they approach the tube walls inside. Three hexane molecules can be encapsulated parallel to each other (i.e., nested) inside a (10,10) SWNT, and their energy gain is 1.98 eV. Three hexane molecules can exhibit a rotary motion. One reason for the stability of carbon chain molecules inside SWNTs is the large area of weak wave function overlap. Another reason concerns molecular dependence, that is, the quadrupole-quadrupole interaction in the case of the polyynes and electron charge transfer from the SWNT in the case of the dehydrogenated forms. The very flat potential surface inside an SWNT suggests that friction is quite low, and the space inside SWNTs serves as an ideal environment for the molecular transport of carbon chain molecules. The present theoretical results are certainly consistent with recent experimental results. Moreover, the encapsulation of C(10) makes an SWNT a (purely carbon-made) p-type acceptor. Another interesting possibility associated with the present system is the direction-controlled transport of C(10)H inside an SWNT under an external field. Because C(10)H has an electric dipole moment, it is expected to move under a gradient electric field. Finally, we derive the entropies of linear chain molecules inside and outside an open-ended SWNT to discuss the stability of including linear chain molecules inside an SWNT at finite temperatures.     

Extrapolation of electron correlation energies to finite and complete basis set targets

The electron correlation energy of two-electron atoms is known to converge asymptotically as approximately (L+1)(-3) to the complete basis set limit, where L is the maximum angular momentum quantum number included in the basis set. Numerical evidence has established a similar asymptotic convergence approximately X(-3) with the cardinal number X of correlation-consistent basis sets cc-pVXZ for coupled cluster singles and doubles (CCSD) and second order perturbation theory (MP2) calculations of molecules. The main focus of this article is to probe for deviations from asymptotic convergence behavior for practical values of X by defining a trial function X(-beta) that for an effective exponent beta=beta(eff)(X,X+1,X+N) provides the correct energy E(X+N), when extrapolating from results for two smaller basis sets, E(X) and E(X+1). This analysis is first applied to "model" expansions available from analytical theory, and then to a large body of finite basis set results (X=D,T,Q,5,6) for 105 molecules containing H, C, N, O, and F, complemented by a smaller set of 14 molecules for which accurate complete basis set limits are available from MP2-R12 and CCSD-R12 calculations. beta(eff) is generally found to vary monotonically with the target of extrapolation, X+N, making results for large but finite basis sets a useful addition to the limited number of cases where complete basis set limits are available. Significant differences in effective convergence behavior are observed between MP2 and CCSD (valence) correlation energies, between hydrogen-rich and hydrogen-free molecules, and, for He, between partial-wave expansions and correlation-consistent basis sets. Deviations from asymptotic convergence behavior tend to get smaller as X increases, but not always monotonically, and are still quite noticeable even for X=5. Finally, correlation contributions to atomization energies (rather than total energies) exhibit a much larger variation of effective convergence behavior, and extrapolations from small basis sets are found to be particularly erratic for molecules containing several electronegative atoms. Observed effects are discussed in the light of results known from analytical theory. A carefully calibrated protocol for extrapolations to the complete basis set limit is presented, based on a single "optimal" exponent beta(opt)(X,X+1,infinity) for the entire set of molecules, and compared to similar approaches reported in the literature.     

Theoretical study of isomerization and dissociation transition states of C3H5O radical isomers: ab initio characterization of the critical points and statistical transition-state theory modeling of the dynamics

I use coupled-cluster theory and a modest basis set, aug-cc-pVDZ, to calculate structures and harmonic vibrational frequencies of local minima and transition states on the C(3)H(5)O potential energy surface. Accurate energies are computed using explicitly correlated coupled-cluster methods and a large basis set, cc-pVQZ-F12, to approach the one-particle basis set limit. My computations characterize eight additional stable radical structures on the global potential energy surface for this system. Additionally, this study encompasses many more isomerization and dissociation pathways, both between previously known intermediates and ones first characterized here. Analysis of the transition states and statistical transition-state theory results shows that energetically small barriers connect many of the alkenol and epoxide intermediates to the straight-chain alkoxy isomers, leading to significant branching to these alkoxy radical intermediates. Facile isomerization to these alkoxy intermediates is significant because the barrier heights leading to H + acrolein and HCO + C(2)H(4) product channels are energetically accessible even at low vibrational energies. The low dissociation barrier heights and loose transition states of these pathways result in unimolecular dissociation as opposed to isomerization to a different C(3)H(5)O intermediate.     

Gas phase hydration and deprotonation of the cyclic C3H3+ cation. Solvation by acetonitrile, and comparison with the benzene radical cation

The binding energies of the first 5 H2O molecules to c-C3H3+ were determined by equilibrium measurements. The measured binding energies of the hydrated clusters of 9-12 kcal/mol are typical of carbon-based CH+...X hydrogen bonds. The ion solvation with the more polar CH3CN molecules results in stronger bonds consistent with the increased ion-dipole interaction. Ab initio calculations show that the lowest energy isomer of the c-C3H3+(H2O)4 cluster consists of a cyclic water tetramer interacting with the c-C3H3+ ion, which suggests the presence of orientational restraint of the water molecules consistent with the observed large entropy loss. The c-C3H3+ ion is deprotonated by 3 or more H2O molecules, driven energetically by the association of the solvent molecules to form strongly hydrogen bonded (H2O)nH+ clusters. The kinetics of the associative proton transfer (APT) reaction C3H3+ + nH2O --> (H2O)nH+ + C3H2* exhibits an unusually steep negative temperature coefficient of k = cT(-63+/-4) (or activation energy of -37 +/- 1 kcal mol(-1)). The behavior of the C3H3+/water system is exactly analogous to the benzene+*/water system, suggesting that the mechanism, kinetics and large negative temperature coefficients may be general to multibody APT reactions. These reactions can become fast at low temperatures, allowing ionized polycyclic aromatics to initiate ice formation in cold astrochemical environments.     

Accurate thermochemistry from corrected Hartree–Fock results: rapid estimation of nearly experimental quality total energy using the small 6-31G(d) basis set

Gaussian-3 ground-state total electronic energies have been approximated using single point 6-31G(d) basis set Harteee–Fock self-consistent-field (HF-SCF) total energies and partial charges based on our earlier rapid estimation of correlation energy from partial charges method. Sixty-five closed-shell neutral molecules (composed of H, C, N, O, and F atoms) of the G2/97 thermochemistry database were selected for the present study. The main feature in this work is that the␣basis set error has been treated by the least squares fit of rapid estimation of basis set error and correlation energy from partial charges (REBECEP) atomic parameters. With these parameters a rather accurate closed-shell ground-state electronic total energy can be obtained from a small basis set HF-SCF calculation in the vicinity of stationary points. The average absolute deviation of the best REBECEP enthalpies of formation from the experimental enthalpies of formation is 1.39 kcal/mol for the test set of 65 enthalpies of neutral molecules. Received: 11 December 2000 / Accepted: 6 February 2001/Published online: 11 October 2001     

H2O nucleation around Au+

First principles electronic structure calculations have been carried out to investigate the ground state geometry, electronic structure, and the binding energy of [Au(H2O)n]+ clusters containing up to 10 H2O molecules. It is shown that the first coordination shell of Au+ contains two H2O molecules forming a H2O-Au+-H2O structure with C2 symmetry. Subsequent H2O molecules bind to the previous H2O molecules forming stable and fairly rigid rings, each composed of 4 H2O molecules, and leading to a dumbbell structure at [Au(H2O)8]+. The 9th and the 10th H2O molecules occupy locations above the Au+ cation mainly bonded to one H2O from each ring, leading to structures where the side rings are partially distorted and forming structures that resemble droplet formation around the Au+ cation. The investigations highlight quantum effects in nucleation at small sizes and provide a microscopic understanding of the observed incremental binding energy deduced from collision induced dissociation that indicates that [Au(H2O)n]+ clusters with 7-10 H2O molecules have comparable binding energy. The charge on the Au+ is shown to migrate to the outside H2O molecules, suggesting an interesting screening phenomenon.     

One and two hydrogen molecules in the large cage of the structure II clathrate hydrate: quantum translation-rotation dynamics close to the cage wall

We have performed a rigorous theoretical study of the quantum translation-rotation (T-R) dynamics of one and two H2 and D2 molecules confined inside the large hexakaidecahedral (5(12)6(4)) cage of the sII clathrate hydrate. For a single encapsulated H2 and D2 molecule, accurate quantum five-dimensional calculations of the T-R energy levels and wave functions are performed that include explicitly, as fully coupled, all three translational and the two rotational degrees of freedom of the hydrogen molecule, while the cage is taken to be rigid. In addition, the ground-state properties, energetics, and spatial distribution of one and two p-H2 and o-D2 molecules in the large cage are calculated rigorously using the diffusion Monte Carlo method. These calculations reveal that the low-energy T-R dynamics of hydrogen molecules in the large cage are qualitatively different from that inside the small cage, studied by us recently. This is caused by the following: (i) The large cage has a cavity whose diameter is about twice that of the small cage for the hydrogen molecule. (ii) In the small cage, the potential energy surface (PES) for H2 is essentially flat in the central region, while in the large cage the PES has a prominent maximum at the cage center, whose height exceeds the T-R zero-point energy of H2/D2. As a result, the guest molecule is excluded from the central part of the large cage, its wave function localized around the off-center global minimum. Peculiar quantum dynamics of the hydrogen molecule squeezed between the central maximum and the cage wall manifests in the excited T-R states whose energies and wave functions differ greatly from those for the small cage. Moreover, they are sensitive to the variations in the hydrogen-bonding topology, which modulate the corrugation of the cage wall.     

Electronic state of small and large cavities for methane hydrate

It is well known that methane hydrate is aggregates of small and large hydrogen bonded water cavities (composed of 12 pentagonal faces of 20 water molecules, and 12 pentagonal and two hexagonal faces of 24 water molecules, respectively) where one methane molecule is encaged. We calculated the methane molecule in vacuum, the small and large cavities by ab initio MO method to clarify the electronic state. The proton of methane in the cavities is shown to form the weak hydrogen bond (O...H[bond]C) between methane and four water molecules, and the H-bond lengths and energies in the small and large cavities were estimated as (0.293 nm, 6.8 kJ/mol) and (0.309 nm, 5.2 kJ/mol), respectively. The calculated values of symmetric C[bond]H stretching frequencies and (13)C-NMR chemical shieldings of the methane in the two cluster cavities show good agreement with the experimental ones observed by Sum et al. and Ripmeester and coworker, respectively.     

Weakly bound complexes of N2O: an ab initio theoretical analysis toward the design of N2O receptors

Ab initio calculations at MP2/6-311++G(2d,2p) and MP2/6-311++G(3df,3pd) computational levels have been used to analyze the interactions between nitrous oxide and a series of small and large molecules that act simultaneously as hydrogen bond donors and electron donors. The basis set superposition error (BSSE) and zero point energy (ZPE) corrected binding energies of small N2O complexes (H2O, NH3, HOOH, HOO*, HONH2, HCO2H, H2CO, HCONH2, H2CNH, HC(NH)NH2, SH2, H2CS, HCSOH, HCSNH2) vary between -0.93 and -2.90 kcal/mol at MP2/6-311++G(3df,3pd) level, and for eight large complexes of N2O they vary between -2.98 and -3.37 kcal/mol at the MP2/6-311++G(2d,2p) level. The most strongly bound among small N2O complexes (HCSNH2-N2O) contains a NH..N bond, along with S-->N interactions, and the most unstable (H2S-N2O) contains just S-->N interactions. The electron density properties have been analyzed within the atoms in molecules (AIM) methodology. Results of the present study open a window into the nature of the interactions between N2O with other molecular moieties and open the possibility to design N2O abiotic receptors.     

The performance of the rapid estimation of basis set error and correlation energy from partial charges method on new molecules of the G3/99 test set

Experimental enthalpies of formation have been approximated using single-point Hartree–Fock (HF)–self-consistent-field (SCF) total energies plus the rapid estimation of basis set error and correlation energy from partial charges (REBECEP) energy corrections. The energy corrections are calculated from the HF–SCF partial atomic charges and optimized atomic energy parameters. The performance of the method was tested on 51 closed-shell neutral molecules (50 molecules from the G3/99 thermochemistry database plus urea, composed of H, C, N, O, and F atoms). The predictive force of the method is demonstrated, because these larger molecules were not used for the optimization of the atomic parameters. We used the earlier RECEP-3 [HF/6-311+G(2d,p)] and REBECEP [HF/6-31G(d)] atomic parameter sets obtained from the G2/97 thermochemistry database (containing small molecules) together with natural population analysis and Mulliken partial charges. The best results were obtained using the natural population analysis charges, although the Mulliken charges also provide useful results. The root-mean-square deviations from the experimental enthalpies of formation for the selected 51 molecules are 1.15, 3.96, and 2.92 kcal/mol for Gaussian-3, B3LYP/6-11+G(3df,2p), and REBECEP (natural population analysis) enthalpies of formation, respectively (the corresponding average absolute deviations are 0.94, 7.09, and 2.27 kcal/mol, respectively). The REBECEP method performs considerably better for the 51 test molecules with a moderate 6-31G(d) basis set than the B3LYP method with a large 6-311+G(3df,2p) basis set. Received: 10 March 2001 / Accepted: 5 July 2001 / Published online: 11 October 2001     

Time-dependent density-functional theory for molecular photoionization with noniterative algorithm and multicenter B-spline basis set: CS2 and C6H6 case studies

In this work a new direct (noniterative) algorithm to solve the time-dependent density-functional theory equations for molecular photoionization has been proposed and implemented, using a multicentric basis set expansion of B-spline functions and complete exploiting of the molecular point-group symmetry. The method has been applied to study the photoionization dynamics of CS2 and C6H6: the results confirmed the expectation of large screening effects in CS2. For C6H6 the screening effects have been found to play a minor role than in CS2, however, also in this case the quality of the final results is definitely improved. The method has proven suitable to study with confidence molecules of medium size, and there is still room for further improvement working on more elaborate treatment of the exchange-correlation functional.