QSID Tool: a new three-dimensional QSAR environmental tool

QSID Tool (Quantitative structure–activity relationship tool for Innovative Discovery) was developed to provide an easy-to-use, robust and high quality environmental tool for 3D QSAR. Predictive models developed with QSID Tool can accelerate the discovery of lead compounds by enabling researchers to formulate and test hypotheses for optimizing efficacy and increasing drug safety and bioavailability early in the process of drug discovery. QSID Tool was evaluated by comparison with SYBYL^® using two different datasets derived from the inhibitors of Trypsin (Böhm et al., J Med Chem 42:458, 1999) and p38-MAPK (Liverton et al., J Med Chem 42:2180, 1999; Romeiro et al., J Comput Aided Mol Des 19:385, 2005; Romeiro et al., J Mol Model 12:855, 2006). The results suggest that QSID Tool is a useful model for the prediction of new analogue activities.     

Diffusion in linear porous media with periodic entropy barriers: A tube formed by contacting spheres

The problem of transport in quasi-one-dimensional periodic structures has been studied recently by several groups [D. Reguera et al., Phys. Rev. Lett.96, 130603 (2006); P. S. Burada et al., Phys. Rev. E75, 051111 (2007); B. Q. Ai and L. G. Liu, ibid.74, 051114 (2006); B. Q. Ai et al., ibid.75, 061126 (2007); B. Q. Ai and L. G. Liu, J. Chem. Phys.126, 204706 (2007); 128, 024706 (2008); E. Yariv and K. D. Dorfman, Phys. Fluids19, 037101 (2007); N. Laachi et al., Europhys. Lett.80, 50009 (2007); A. M. Berezhkovskii et al., J. Chem. Phys.118, 7146 (2003); 119, 6991 (2003)]. Using the concept of "entropy barrier" [R. Zwanzig, J. Phys. Chem.96, 3926 (1992)] one can classify such structures based on the height of the entropy barrier. Structures with high barriers are formed by chambers, which are weakly connected with each other because they are connected by small apertures. To escape from such a chamber a diffusing particle has to climb a high entropy barrier to find an exit that takes a lot of time [I. V. Grigoriev et al., J. Chem. Phys.116, 9574 (2002)]. As a consequence, the particle intrachamber lifetime tau(esc) is much larger than its intrachamber equilibration time, tau(rel), tau(esc)>tau(rel). When the aperture is not small enough, the intrachamber escape and relaxation times are of the same order and the hierarchy fails. This is the case of low entropy barriers. Transport in this case is analyzed in the works of Schmid and co-workers, Liu and co-workers, and Dorfman and co-workers, while the work of Berezhkovskii et al. is devoted to diffusion in the case of high entropy barriers.     

Variable selection and specification of robust QSAR models from multicollinear data: arylpiperazinyl derivatives with affinity and selectivity for α<Subscript>2</Subscript>-adrenoceptors

Two QSAR models have been identified that predict the affinity and selectivity of arylpiperazinyl derivatives for alpha1 and alpha2 adrenoceptors (ARs). The models have been specified and validated using 108 compounds whose structures and inhibition constants (Ki) are available in the literature [Barbaro et al., J. Med. Chem., 44 (2001) 2118; Betti et al., J. Med. Chem., 45 (2002) 3603; Barbaro et al., Bioorg. Med. Chem., 10 (2002) 361; Betti et al., J. Med. Chem., 46 (2003) 3555]. One hundred and forty-seven predictors have been calculated using the Cerius 2 software available from Accelrys. This set of variables exhibited redundancy and severe multicollinearity, which had to be identified and removed as appropriate in order to obtain robust regression models free of inflated errors for the beta estimates - so-called bouncing betas. Those predictors that contained information relevant to the alpha2 response were identified on the basis of their pairwise linear correlations with affinity (-log Ki) for alpha2 adrenoceptors; the remaining variables were discarded. Subsequent variable selection made use of Factor Analysis (FA) and Unsupervised Variable Selection (UzFS). The data was divided into test and training sets using cluster analysis. These two sets were characterised by similar and consistent distributions of compounds in a high dimensional, but relevant predictor space. Multiple regression was then used to determine a subset of predictors from which to determine QSAR models for affinity to alpha2-ARs. Two multivariate procedures, Continuum Regression (the Portsmouth formulation) and Canonical Correlation Analysis (CCA), have been used to specify models for affinity and selectivity, respectively. Reasonable predictions were obtained using these in silico screening tools.     

GA strategy for variable selection in QSAR studies: application of GA-based region selection to a 3D-QSAR study of acetylcholinesterase inhibitors

Comparative molecular field analysis (CoMFA) with partial least squares (PLS) is one of the most frequently used tools in three-dimensional quantitative structure-activity relationships (3D-QSAR) studies. Although many successful CoMFA applications have proved the value of this approach, there are some problems in its proper application. Especially, the inability of PLS to handle the low signal-to-noise ratio (sample-to-variable ratio) has attracted much attention from QSAR researchers as an exciting research target, and several variable selection methods have been proposed. More recently, we have developed a novel variable selection method for CoMFA modeling (GARGS: genetic algorithm-based region selection), and its utility has been demonstrated in the previous paper (Kimura, T., et al. J. Chem. Inf. Comput. Sci. 1998, 38, 276-282). The purpose of this study is to evaluate whether GARGS can pinpoint known molecular interactions in 3D space. We have used a published set of acetylcholinesterase (AChE) inhibitors as a test example. By applying GARGS to a data set of AChE inhibitors, several improved models with high internal prediction and low number of field variables were obtained. External validation was performed to select a final model among them. The coefficient contour maps of the final GARGS model were compared with the properties of the active site in AChE and the consistency between them was evaluated.     

Exploration of a binding mode of benzothiazol-2-yl acetonitrile pyrimidine core based derivatives as potent c-Jun N-terminal kinase-3 inhibitors and 3D-QSAR analyses

C-Jun N-terminal kinase (JNK) is a therapeutic target for inhibitors which may provide clinical benefit in the pathogenesis of rheumatoid arthritis (RA) as well as in various apoptosis-related disorders. The benzothiazol-2-yl acetonitrile derivatives, recently reported by Pascale et al. (J. Med. Chem. 2005, 48, 4596-4607), are the first generation JNK inhibitors of this class. To understand inhibitory mechanisms and elucidate pharmacophoric properties of these derivatives molecular docking and 3D-QSAR studies were performed on a set of 44 compounds. Ligand Fit module of Cerius2 (4.9) was employed to locate the binding orientations of all the compounds within the JNK-3 ATP binding site. A good correlation (r2=0.810) between the calculated binding free energies (-PMF score) and the experimental inhibitory activities suggests that the identified binding conformations of these potential inhibitors are reliable. Based on the binding conformations, robust and highly predictive 3D-QSAR models were developed with conventional r2 0.886 and 0.802, full cross-validation r2 0.980 and 0.788, and predictive r2 0.965 and 0.968 for MFA and MSA, respectively. The interaction mode was demonstrated taking into consideration inhibitor conformation, hydrogen bonding, and electrostatic interaction. The 3D-QSAR model built in this study will provide clear guidelines for a novel inhibitor design based on the benzothiazole derivatives against JNK-3 for the treatment of inflammatory disorders.     

Theoretical study of pattern formation during the catalytic oxidation of CO on Pt{100} at low pressures

Theoretical studies have thus far been unable to model pattern formation during the reaction in this system on physically feasible length and time scales. In this paper, we derive a computational reaction-diffusion model for this system in which most of the input parameters have been determined experimentally. We model the surface on a mesoscopic scale intermediate between the microscopic size of CO islands and the macroscopic length scale of pattern formation. In agreement with experimental investigations [M. Eiswirth et al., Z. Phys. Chem., Neue Folge 144, 59 (1985)], the results from our model divide the CO and O(2) partial pressure parameter space into three regions defined by the level of CO coverage or the presence of sustained oscillations. We see CO fronts moving into oxygen-covered regions, with the 1 x 1 to hex phase change occurring at the leading edge. There are also traveling waves consisting of successive oxygen and CO fronts that move into areas of relatively high CO coverage, and in this case, the phase change is more gradual and of lower amplitude. The propagation speed of these reaction waves is similar to those observed experimentally for CO and oxygen fronts [H. H. Rotermund et al., J. Chem. Phys. 91, 4942 (1989); H. H. Rotermund et al., Nature (London) 343, 355 (1990); J. Lauterbach and H. H. Rotermund, Surf. Sci. 311, 231 (1994)]. In the two-dimensional version of our model, the traveling waves take the form of target patterns emitted from surface inhomogeneities.     

CoMFA and docking studies on glycogen phosphorylase a inhibitors as antidiabetic agents

Glycogen phosphorylase (GP(a)) is a specific target for the design of inhibitors and may prevent glycogenolysis under high glucose conditions in type II diabetes. The carboxamides first reported by Hoover D. J. et al. (J. Med. Chem. 1998, 41, 2934-2938) are one of the major classes of GP(a) inhibitors other than glucose derivatives. The recent, X-ray crystallographic analyses (Oikonomakos et al. Biochim. Biophys. Acta 2003, 1647, 325-332) have revealed a distinct mechanism of action for these inhibitors, which bind at a new allosteric site away from the inhibitory and catalytic sites. To elucidate the essential structural and physicochemical requirements responsible for binding to the GP(a) enzyme and to develop predictive models, CoMFA and docking studies have been carried out on a series of indole-2-carboxamide derivates. The CoMFA model developed using pharmacophoric alignments and hydrogen-bonding fields demonstrated high predictive ability against the training (r2 = 0.98, q2 = 0.68) and the test set (r2pred = 0.85). Further the superimposition of PLS coefficient contour maps from CoMFA with the GP(a) active site (PDB: 1lwo) has shown a high level of compatibility.     

H2 storage in isostructural UiO-67 and UiO-66 MOFs

The recently discovered UiO-66/67/68 class of isostructural metallorganic frameworks (MOFs) [J. H. Cavka et al. J. Am. Chem. Soc., 2008, 130, 13850] has attracted great interest because of its remarkable stability at high temperatures, high pressures and in the presence of different solvents, acids and bases [L. Valenzano et al. Chem. Mater., 2011, 23, 1700]. UiO-66 is obtained by connecting Zr(6)O(4)(OH)(4) inorganic cornerstones with 1,4-benzene-dicarboxylate (BDC) as linker resulting in a cubic MOF, which has already been successfully reproduced in several laboratories. Here we report the first complete structural, vibrational and electronic characterization of the isostructural UiO-67 material, obtained using the longer 4,4'-biphenyl-dicarboxylate (BPDC) linker, by combining laboratory XRPD, Zr K-edge EXAFS, TGA, FTIR, and UV-Vis studies. Comparison between experimental and periodic calculations performed at the B3LYP level of theory allows a full understanding of the structural, vibrational and electronic properties of the material. Both materials have been tested for molecular hydrogen storage at high pressures and at liquid nitrogen temperature. In this regard, the use of a longer ligand has a double benefit: (i) it reduces the density of the material and (ii) it increases the Langmuir surface area from 1281 to 2483 m(2) g(-1) and the micropore volume from 0.43 to 0.85 cm(3) g(-1). As a consequence, the H(2) uptake at 38 bar and 77 K increases from 2.4 mass% for UiO-66 up to 4.6 mass% for the new UiO-67 material. This value is among the highest values reported so far but is lower than those reported for MIL-101, IRMOF-20 and MOF-177 under similar pressure and temperature conditions (6.1, 6.2 and 7.0 mass%, respectively) [A. G. Wong-Foy et al. J. Am. Chem. Soc., 2006, 128, 3494; M. Dinca and J. R. Long. Angew. Chem., Int. Ed., 2008, 47, 6766]. Nevertheless the remarkable chemical and thermal stability of UiO-67 and the absence of Cr in its structure would make this material competitive.     

Designing novel nicotinic agonists by searching a database of molecular shapes

Summary We introduce an approach by which novel ligands can be designed for a receptor if a pharmacophore geometry has been established and the receptor-bound conformations of other ligands are known. We use the shape-matching method of Kuntz et al. [J. Mol. Biol., 161 (1982) 269–288] to search a database of molecular shapes for those molecules which can fit inside the combined volume of the known ligands and which have interatomic distances compatible with the pharmacophore geometry. Some of these molecules are then modified by interactive modeling techniques to better match the chemical properties of the known ligands. Our shape database (about 5000 candidate molecules) is derived from a subset of the Cambridge Crystallographic Database [Allen et al., Acta Crystallogr., Sect. B,35 (1979) 2331–2339]. We show, as an example, how several novel designs for nicotinic agonists can be derived by this approach, given a pharmacophore model derived from known agonists [Sheridan et al., J. Med. Chem., 29 (1986) 889–906]. This report complements our previous report [DesJarlais et al., J. Med. Chem., in press], which introduced a similar method for designing ligands when the structure of the receptor is known.     

Characterization of the structure and reactivity of monocopper-oxygen complexes supported by beta-diketiminate and anilido-imine ligands

Copper-oxygen complexes supported by beta-diketiminate and anilido-imine ligands have recently been reported (Aboelella et al., J Am Chem Soc 2004, 126, 16896; Reynolds et al., Inorg Chem 2005, 44, 6989) as potential biomimetic models for dopamine beta-monooxygenase (DbetaM) and peptidylglycine alpha-hydroxylating monooxygenase (PHM). However, in contrast to the enzymatic systems, these complexes fail to exhibit C--H hydroxylation activity (Reynolds et al., Chem Commun 2005, 2014). Quantum chemical characterization of the 1:1 Cu-O(2) model adducts and related species (Cu(III)-hydroperoxide, Cu(III)-oxo, and Cu(III)-hydroxide) indicates that the 1:1 Cu-O(2) adducts are unreactive toward substrates because of the weakness of the O--H bond that would be formed upon hydrogen-atom abstraction. This in turn is ascribed to the 1:1 adducts having both low reduction potentials and basicities. Cu(III)-oxo species on the other hand, determined to be intermediate between Cu(III)-oxo and Cu(II)-oxyl in character, are shown to be far more reactive toward substrates. Based on these results, design strategies for new DbetaM and PHM biomimetic ligands are proposed: new ligands should be made less electron rich so as to favor end-on dioxygen coordination in the 1:1 Cu-O(2) adducts. Comparison of the relative reactivities of the various copper-oxygen complexes as hydroxylating agents provides support for a Cu(II)-superoxide species as the intermediate responsible for substrate hydroxylation in DbetaM and PHM, and suggests that a Cu(III)-oxo intermediate would be competent in this process as well.     

Molecular beam optical Stark study of rhodium mononitride

The optical Stark effect in the Q(1) and R(0) lines of the [15.1]1-X (1)Sigma+ (1,0) band of rhodium mononitride (RhN) were recorded and analyzed to determine the permanent electric dipole moments mu for the X (1)Sigma+(upsilon=0) and [15.1]1(upsilon=1) states to be 2.43(5) and 1.75(1) D, respectively. The determined dipole moments are compared to predicted values obtained from density functional theory [Stevens et al., Chem. Phys. Lett. 421, 281 (2006)] and an all-electron ab initio calculation [Shim et al., J. Mol. Struct. THEOCHEM 393, 127 (1997)]. A simple single configuration molecular orbital correlation diagram is used to rationalize the relative values of mu for the 4d mononitrides and RhO. An electronic configuration for the [15.1]1 state is proposed based on the interpretation of the (103)Rh and (14)N magnetic hyperfine interactions.     

Three-dimensional potential energy surface of the Ar-OH(2Pi i) complex

Pure rotational transitions in the ground state for Ar-OH and Ar-OD [Y. Ohshima et al., J. Chem. Phys. 95, 7001 (1991) and Y. Endo et al., Faraday Discuss. 97, 341 (1994)], those in the excited states of the OH vibration, nu(s)=1 and 2, observed by Fourier-transform microwave spectroscopy in the present study, rotation-vibration transitions observed by infrared-ultraviolet double-resonance spectroscopy [K. M. Beck et al., Chem. Phys. Lett. 162, 203 (1989) and R. T. Bonn et al., J. Chem. Phys. 112, 4942 (2000)], and the P-level structure observed by stimulated emission pumping spectroscopy [M. T. Berry et al., Chem. Phys. Lett. 178, 301 (1991)] have been simultaneously analyzed to determine the potential energy surface of Ar-OH in the ground state. A Schrodinger equation, considering all the freedom of motions for an atom-diatom system in the Jacobi coordinate, R, theta, and r, was numerically solved to obtain energies of the rovibrational energy levels using the discrete variable representation method. A three-dimensional potential energy surface is determined by a least-squares fitting. In the analysis the potential parameters, obtained by ab initio calculations at the RCCSD(T) level of theory with a set of basis functions of aug-cc-pVTZ and midbond functions, are used as initial values. The determined intermolecular potential energy surface and its dependence on the OH monomer bond length are compared with those of an isovalent radical complex, Ar-SH.     

From homogeneous dispersion to micelles-a molecular dynamics simulation on the compromise of the hydrophilic and hydrophobic effects of sodium dodecyl sulfate in aqueous solution

The structural and functional diversity of surfactant systems has attracted simulation works in atomistic, coarse grain, and mesoscopic models (Bandyopadhyay, S.; et al. Langmuir 2000, 16, 942; Senapati, S.; et al. J. Phys. Chem. B 2003, 107, 12906; Maiti, P. K.; et al. Langmuir 2002, 18, 1908; Srinivas, G.; et al. J. Phys. Chem. B 2004, 108, 8153; Groot, R. D.; et al. J. Chem. Phys. 1999, 110, 9739; Rekvig, L.; et al. Langmuir 2003, 19, 8195). However, atomistic models have suffered from their tremendous computational cost and are, so far, not able to simulate the structural behaviors in sufficient spatio-temporal scales (Shelley, J. C.; Shelley, M. Y. Curr. Opin. Colloid Interface Sci. 2000, 5, 101). The other two approaches are not microscopic enough to describe the configurations of the surfactants that determine their behaviors (Shelley and Shelley). In this study, we propose to simplify atomistic models based on the observation that the compromise of the hydrophilic and hydrophobic effects (Li, J.; Kwauk, M. Chem. Eng. Sci. 2003, 58, 521-535) and molecular structures of surfactants are the dominant factors shaping their structures in the systems. With this simplification, we are able to simulate with moderate computing cost the whole process of micelle formation from an initially uniform dispersion of sodium dodecyl sulfate (SDS) in aqueous solution. The resulting micelle structures are different from those predicted by atomistic simulations that started with a predefined micelle configuration at the same surfactant concentrations. However, if we use their initial micelle configuration, micelle structures the same as theirs are obtained. Analyses show that our results are more realistic and that the results of the atomistic simulations suffer from artificial initial conditions. Therefore, our model may serve as a reasonable simplification of atomistic models in terms of the general structure of micelles.     

Comparison of different methods for calculating the paramagnetic relaxation enhancement of nuclear spins as a function of the magnetic field

The enhancement of the spin-lattice relaxation rate for nuclear spins in a ligand bound to a paramagnetic metal ion [known as the paramagnetic relaxation enhancement (PRE)] arises primarily through the dipole-dipole (DD) interaction between the nuclear spins and the electron spins. In solution, the DD interaction is modulated mostly by reorientation of the nuclear spin-electron spin axis and by electron spin relaxation. Calculations of the PRE are in general complicated, mainly because the electron spin interacts so strongly with the other degrees of freedom that its relaxation cannot be described by second-order perturbation theory or the Redfield theory. Three approaches to resolve this problem exist in the literature: The so-called slow-motion theory, originating from Swedish groups [Benetis et al., Mol. Phys. 48, 329 (1983); Kowalewski et al., Adv. Inorg. Chem. 57, (2005); Larsson et al., J. Chem. Phys. 101, 1116 (1994); T. Nilsson et al., J. Magn. Reson. 154, 269 (2002)] and two different methods based on simulations of the dynamics of electron spin in time domain, developed in Grenoble [Fries and Belorizky, J. Chem. Phys. 126, 204503 (2007); Rast et al., ibid. 115, 7554 (2001)] and Ann Arbor [Abernathy and Sharp, J. Chem. Phys. 106, 9032 (1997); Schaefle and Sharp, ibid. 121, 5387 (2004); Schaefle and Sharp, J. Magn. Reson. 176, 160 (2005)], respectively. In this paper, we report a numerical comparison of the three methods for a large variety of parameter sets, meant to correspond to large and small complexes of gadolinium(III) and of nickel(II). It is found that the agreement between the Swedish and the Grenoble approaches is very good for practically all parameter sets, while the predictions of the Ann Arbor model are similar in a number of the calculations but deviate significantly in others, reflecting in part differences in the treatment of electron spin relaxation. The origins of the discrepancies are discussed briefly.     

Hybrid one-electron/many-electron methods for ionized states of molecular clusters

To describe singly-ionized states of molecular clusters we devised an effective Hamiltonian approach that combines (1) accurate monomer ionization potentials from many-electron wave functions with (2) polarization shifts and (3) effective monomer couplings obtained from a simple one-electron approach (the superposition-of-fragment-states (SFS) method [Valeev et al., J. Am. Chem. Soc., 2006, 128, 9882]). The accuracy of the intermolecular coupling parameters evaluated with SFS Hartree-Fock (HF) and Density-Functional-Theory (DFT) variants was evaluated for several weakly-bound dimers and compared against the state-of-the-art equation-of-motion ionization-potential coupled-cluster singles and doubles (EOM-IP-CCSD) data of Krylov and co-workers. The SFS-HF method produces coupling integrals accurate to a few percent, whereas SFS-DFT predictions are substantially worse. A hybrid approach combining SFS-HF couplings and shifts with EOM-IP-CCSD ionization potentials of monomers (denoted as SFS-EOM-IP-CCSD) was applied to ionized states of two conformers of a benzene dimer and ten representative DNA base pairs. The 16 considered SFS-EOM-IP-CCSD ionization potentials of the benzene dimer differed from the reference EOM-IP-CCSD IPs of Krylov and co-workers [Pieniazek et al., J. Chem. Phys. 2007, 127, 044317; Bravaya et al., Phys. Chem. Chem. Phys. 2010, 12, 2261] by less than 0.1 eV on average, and at most by 0.2 eV. For the DNA base pairs the mean absolute (median) deviation of the SFS-EOM-IP-CCSD IPs was 0.27 (0.23) eV; several deviations for non-Koopmans states were as large as 0.9 eV. The SFS-EOM-IP-CCSD method can be readily applied to large molecular clusters with computational effort scaling cubically with the size of the cluster.     

Quantum Monte Carlo calculations of the dissociation energy of the water dimer

We report diffusion quantum Monte Carlo (DMC) calculations of the equilibrium dissociation energy D(e) of the water dimer. The dissociation energy measured experimentally, D(0), can be estimated from D(e) by adding a correction for vibrational effects. Using the measured dissociation energy and the modern value of the vibrational energy Mas et al., [J. Chem. Phys. 113, 6687 (2000)] leads to D(e)=5.00+/-0.7 kcal mol(-1), although the result Curtiss et al., [J. Chem. Phys. 71, 2703 (1979)] D(e)=5.44+/-0.7 kcal mol(-1), which uses an earlier estimate of the vibrational energy, has been widely quoted. High-level coupled cluster calculations Klopper et al., [Phys. Chem. Chem. Phys. 2, 2227 (2000)] have yielded D(e)=5.02+/-0.05 kcal mol(-1). In an attempt to shed new light on this old problem, we have performed all-electron DMC calculations on the water monomer and dimer using Slater-Jastrow wave functions with both Hartree-Fock approximation (HF) and B3LYP density functional theory single-particle orbitals. We obtain equilibrium dissociation energies for the dimer of 5.02+/-0.18 kcal mol(-1) (HF orbitals) and 5.21+/-0.18 kcal mol(-1) (B3LYP orbitals), in good agreement with the coupled cluster results.     

Estimated MP2 and CCSD(T) interaction energies of n-alkane dimers at the basis set limit: comparison of the methods of Helgaker et al. and Feller

The MP2 (the second-order M?ller-Plesset calculation) and CCSD(T) (coupled cluster calculation with single and double substitutions with noniterative triple excitations) interaction energies of all-trans n-alkane dimers were calculated using Dunning's [J. Chem. Phys. 90, 1007 (1989)] correlation consistent basis sets. The estimated MP2 interaction energies of methane, ethane, and propane dimers at the basis set limit [EMP2(limit)] by the method of Helgaker et al. [J. Chem. Phys. 106, 9639 (1997)] from the MP2/aug-cc-pVXZ (X=D and T) level interaction energies are very close to those estimated from the MP2/aug-cc-pVXZ (X=T and Q) level interaction energies. The estimated EMP2(limit) values of n-butane to n-heptane dimers from the MP2/cc-pVXZ (X=D and T) level interaction energies are very close to those from the MP2/aug-cc-pVXZ (X=D and T) ones. The EMP2(limit) values estimated by Feller's [J. Chem. Phys. 96, 6104 (1992)] method from the MP2/cc-pVXZ (X=D, T, and Q) level interaction energies are close to those estimated by the method of Helgaker et al. from the MP2/cc-pVXZ (X=T and Q) ones. The estimated EMP2(limit) values by the method of Helgaker et al. using the aug-cc-pVXZ (X=D and T) are close to these values. The estimated EMP2(limit) of the methane, ethane, propane, n-butane, n-pentane, n-hexane, n-heptane, n-octane, n-nonane, and n-decane dimers by the method of Helgaker et al. are -0.48, -1.35, -2.08, -2.97, -3.92, -4.91, -5.96, -6.68, -7.75, and -8.75 kcal/mol, respectively. Effects of electron correlation beyond MP2 are not large. The estimated CCSD(T) interaction energies of the methane, ethane, propane, and n-butane dimers at the basis set limit by the method of Helgaker et al. (-0.41, -1.22, -1.87, and -2.74 kcal/mol, respectively) from the CCSD(T)/cc-pVXZ (X=D and T) level interaction energies are close to the EMP2(limit) obtained using the same basis sets. The estimated EMP2(limit) values of the ten dimers were fitted to the form m0+m1X (X is 1 for methane, 2 for ethane, etc.). The obtained m0 and m1 (0.595 and -0.926 kcal/mol) show that the interactions between long n-alkane chains are significant. Analysis of basis set effects shows that cc-pVXZ (X=T, Q, or 5), aug-cc-pVXZ (X=D, T, Q, or 5) basis set, or 6-311G** basis set augmented with diffuse polarization function is necessary for quantitative evaluation of the interaction energies between n-alkane chains.     

Reactivity enhancement of ultracold O(3P)+H2 collisions by van der Waals interactions

The role of van der Waals forces in O((3)P)+H(2)(upsilon=1,j=0) collisions is investigated theoretically at low and ultralow temperatures. Quantum scattering calculations have been performed for zero total angular momentum using the lowest London-Eyring-Polanyi-Sato double-polynomial (3)A(") potential-energy surface reported by [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)] and its recent BMS1 and BMS2 extensions developed by [Brandao et al., J. Chem. Phys. 121, 8861 (2004)] which provide a more accurate treatment of the van der Waals interaction. Our calculations show that van der Waals forces strongly influence chemical reactivity at ultracold translational energies. The presence of a zero-energy resonance for the BMS1 surface is found to enhance reactivity in the ultracold regime and shift the Wigner threshold to lower temperatures.     

Molecular dynamics simulations of surface-initiated melting of nitromethane

The melting of nitromethane initiated at solid-vacuum interfaces has been investigated using molecular dynamics nvt simulations with a realistic force field [D. C. Sorescu et al., J. Phys. Chem. B 104, 8406 (2000)]. The calculated melting point (251+/-5 K) is in good agreement with experiment (244.73 K) and values obtained previously (approximately 255.5 and 266.5+/-8 K) using other simulation methods [P. M. Agrawal et al., J. Chem. Phys. 119, 9617 (2003)]. Analyses of the molecular orientations and diffusion during the simulations as functions of the distance from the exposed surfaces show that the melting is a direct crystal-to-liquid transition, in which the molecules first gain rotational freedom, then mobility. There is a slight dependence of the melting temperature on the exposed crystallographic face.