Biomedical subjects

Mark S Gordon

Publications and source records attributed to Mark S Gordon.

At least 19 recordsLinked to original sources

PubMed2006

The structure of the Si9H12 cluster: A coupled cluster and multi-reference perturbation theory study.

Full geometry optimizations using both singles and doubles coupled cluster theory with perturbative triple excitations, CCSD(T), and second order multi-reference perturbation theory, MRMP2, have been employed to predict the structure of Si9H12, a cluster commonly used in calculations to represent the Si(100) surface. Both levels of theory predict the structure of this cluster to be symmetric (not buckled), and no evidence for a buckled (asymmetric) structure is found at either level of theory.

Journal Article↗

PubMed2006

Dissociation potential curves of low-lying states in transition metal hydrides. 3. Hydrides of groups 6 and 7.

The dissociation curves of low-lying spin-mixed states in monohydrides of groups 6 and 7 were calculated by using an effective core potential (ECP) approach. This approach is based on the multiconfiguration self-consistent field (MCSCF) method, followed by first-order configuration interaction (FOCI) calculations, in which the method employs an ECP basis set proposed by Stevens and co-workers (SBKJC) augmented by a set of polarization functions. Spin-orbit coupling (SOC) effects are estimated within the one-electron approximation by using effective nuclear charges, since SOC splittings obtained with the full Breit-Pauli Hamitonian are underestimated when ECP basis sets are used. The ground states of group 6 hydrides have Omega = (1)/(2)(X(6)Sigma(+)(1/2)), where Omega is the z component of the total angular momentum quantum number. Although the ground states of group 7 hydrides have Omega = 0(+), their main adiabatic components are different; the ground state in MnH originates from the lowest (7)Sigma(+), while in TcH and ReH the main component of the ground state is the lowest (5)Sigma(+). The present paper reports a comprehensive set of theoretical results including the dissociation energies, equilibrium distances, electronic transition energies, harmonic frequencies, anharmonicities, and rotational constants for several low-lying spin-mixed states in these hydrides. Transition dipole moments were also computed among the spin-mixed states and large peak positions of electronic transitions are suggested theoretically for these hydrides. The periodic trends of physical properties of metal hydrides are discussed, based on the results reported in this and other recent studies.

Journal Article↗

PubMed2006

A systematic multireference perturbation-theory study of the low-lying states of SiC3.

The three known lowest-energy isomers of SiC(3), two cyclic singlets (2s and 3s) and a linear triplet (1t), have been reinvestigated using multireference second-order perturbation theory (MRPT2). The dependence of the relative energies of the isomers upon the quality of the basis sets and the sizes of the reference active spaces is explored. When using a complete-active-space self-consistent-field reference wave function with 12 electrons in 11 orbitals [CASSCF (12, 11)] together with basis sets that increase in size up to the correlation-consistent polarized core-valence quadruple zeta basis set (cc-pCVQZ), the MRPT2 method consistently predicts the linear triplet to be the most stable isomer. A new parallel direct determinant MRPT2 code has been used to systematically explore reference spaces that vary in size from CASSCF (8,8) to full optimized reaction space [FORS or CASSCF (16,16)] with the cc-pCVQZ basis. It is found that the relative energies of the isomers change substantially as the active space is increased. At the best level of theory, MRPT2 with a full valence FORS reference, the 2s isomer is predicted to be more stable than 3s and 1t by 4.7 and 2.2 kcal/mol, respectively.

Carbon Compounds, Inorganic↗

PubMed2006

Modeling styrene-styrene interactions.

This study is the first step in the systematic investigation of substituted (carboxyl) polystyrene nanoparticles. Understanding the fundamental interactions between the p-carboxyl styrene monomers, where an ethyl group is used instead of a vinyl group (referenced, for convenience, as "p-carboxyl styrene"), provides the basic information needed to construct potentials for nanoparticles composed of these monomers. In this work, low-energy isomers of p-carboxyl styrene dimer were studied. The dimer structures and their relative and binding energies were determined using both Møller-Plesset second-order perturbation theory (MP2) and the general effective fragment potential (EFP2) method. Sections of the intermolecular potential energy surface (PES) of the p-carboxylated styrene dimer in its global minimum orientation were also determined. As expected, double hydrogen bonding between the two carboxylic groups provides the strongest interaction in this system, followed by isomers with a single H-bond and strong benzene ring-benzene ring (pi-pi) type interactions. Generally, the EFP2 method reproduces the MP2 geometries and relative energies with good accuracy, so it appears to be an efficient alternative to the correlated ab initio methods, which are too computationally demanding to be routinely used in the study of the more-complex polymeric systems of interest.

Journal Article↗

PubMed2006

Scalable implementation of analytic gradients for second-order Z-averaged perturbation theory using the distributed data interface.

The analytic gradient expression for second-order Z-averaged perturbation theory is revised and its parallel implementation is described in detail. The distributed data interface is used to access molecular-orbital integral arrays stored in distributed memory. The algorithm is designed to maximize the use of local data and reduce communication costs. The iterative solution and the preconditioner used to induce the convergence of the coupled-perturbed Hartree-Fock equations are presented. Several illustrative timing examples are discussed.

Journal Article↗

PubMed2006

Parallel coupled perturbed CASSCF equations and analytic CASSCF second derivatives.

A parallel algorithm for solving the coupled-perturbed MCSCF (CPMCSCF) equations and analytic nuclear second derivatives of CASSCF wave functions is presented. A parallel scheme for evaluating derivative integrals and their subsequent use in constructing other derivative quantities is described. The task of solving the CPMCSCF equations is approached using a parallelization scheme that partitions the electronic hessian matrix over all processors as opposed to simple partitioning of the 3 N solution vectors among the processors. The scalability of the current algorithm, up to 128 processors, is demonstrated. Using three test cases, results indicate that the parallelization of derivative integral evaluation through a simple scheme is highly effective regardless of the size of the basis set employed in the CASSCF energy calculation. Parallelization of the construction of the MCSCF electronic hessian during solution of the CPMCSCF equations varies quantitatively depending on the nature of the hessian itself, but is highly scalable in all cases.

Journal Article↗

PubMed2005

Influence of multi-atom bridging ligands on the electronic structure and magnetic properties of homodinuclear titanium molecules.

The electronic structure and magnetic properties of homodinuclear titanium(III) molecules with bridging ligands from groups 14, 15, and 16 are examined. Single- and multireference methods with triple-zeta plus polarization basis sets are employed. Dynamic electron correlation effects are included via second-order multireference perturbation theory. Isotropic interaction parameters are calculated, and two of the complexes studied are predicted to be ferromagnetic based on multireference second-order perturbation (MRMP2) theory, using the TZVP(fg) basis set. Zero-field splitting parameters are determined using spin-orbit coupling obtained from complete active space (CAS) self-consistent field (SCF) and multiconfigurational quasi-degenerate perturbation theory (MCQDPT) wave functions. Three Breit-Pauli-based spin coupling methods were employed: full Breit-Pauli (HSO2), the partial two-electron method (P2E), and the semiempirical one-electron method (HSO1).

Computer Simulation↗

PubMed2005

Time-resolved gas-phase kinetic and quantum chemical studies of the reaction of silylene with oxygen.

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with O(2). The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF(6) bath gas, at five temperatures in the range 297-600 K. The second order rate constants at 10 Torr were fitted to the Arrhenius equation: [see text] The decrease in rate constant values with increasing temperature, although systematic is very small. The rate constants showed slight increases in value with pressure at each temperature, but this was scarcely beyond experimental uncertainty. From estimates of Lennard-Jones collision rates, this reaction is occurring at ca. 1 in 20 collisions, almost independent of pressure and temperature. Ab initio calculations at the G3 level backed further by multi-configurational (MC) SCF calculations, augmented by second order perturbation theory (MRMP2), support a mechanism in which the initial adduct, H(2)SiOO, formed in the triplet state (T), undergoes intersystem crossing to the more stable singlet state (S) prior to further low energy isomerisation processes leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are H2O+SiO. The decomposition of the intermediate cyclo-siladioxirane, via O-O bond fission, plays an important role in the overall process. The bottleneck for the overall process appears to be the T-->S process in H2SiOO. This process has a small spin-orbit coupling matrix element, consistent with an estimate of its rate constant of 1x10(9) s-1 obtained with the aid of RRKM theory. This interpretation preserves the idea that, as in its reactions in general, SiH2 initially reacts at the encounter rate with O2. The low values for the secondary reaction barriers on the potential energy surface account for the lack of an observed pressure dependence. Some comparisons are drawn with the reactions of CH2+O2 and SiCl2+O2.

Algorithms↗

PubMed2005

Cycloaddition isomerizations of adsorbed 1,3-cyclohexadiene on Si(100)-2x1 surface: first neighbor interactions.

The initial and subsequent surface reaction mechanisms of 1,3-cyclohexadiene on the Si(100)-2x1 surface were theoretically explored, focusing on the possible first-neighbor interactions. Five different initial reaction channels leading to nine different surface products were identified, confirming previous experimental reports of inter-dimer structures. Among the nine identified products, five of these surface products are new species that have not previously been reported. Potential energy surface studies reveal that the net reaction barriers within a given channel are very small, indicating that the final product distributions within that channel are determined by thermodynamics. On the other hand, thermal isomerizations between different channels are not expected to occur easily. Therefore, the surface product distributions among the five different channels are more likely to be determined by kinetics. As a result, understanding the relationships among the available reaction channels both kinetically and thermodynamically is essential for properly interpreting the experimental results. The current study shows that the subsequent surface chemical reactions of unsaturated initial surface products are strongly coupled with the first-neighbor interactions.

Journal Article↗

PubMed2005

Cycloaddition of benzene on Si(100) and its surface conversions.

A comprehensive ab initio study of the adsorption of benzene on the silicon(100) surface is presented. Five potential candidates ([2+2] adduct, [4+2] adduct, two tetra-sigma-bonded structures, and one radical-like structure) for the reaction product are examined to determine the lowest energy adsorption configuration. A [4+2] butterfly structure is determined to be the global minimum (-29.0 kcal/mol), although one of the two tetra-sigma-bonded structures (-26.7 kcal/mol) is similar in energy to it. Multireference perturbation theory suggests that the [4+2] addition mechanism of benzene on Si(100) is very similar to the usual Diels-Alder reaction (i.e., small or zero activation barrier), even though benzene adsorption entails the loss of benzene aromaticity during the reaction. On the other hand, the [2+2] cycloaddition mechanism is shown to require a relatively high activation barrier (17.8 kcal/mol), in which the initial step is to form a (relatively strongly bound) van der Waals complex (-8.9 kcal/mol). However, the net activation barrier relative to reactants is only 8.9 kcal/mol. Careful examination of the interconversion reactions among the reaction products indicates that the two tetra-sigma-bonded structures (that are energetically comparable to the [4+2] product) can be derived from the [2+2] adduct with activation barriers of 15.5 and 21.4 kcal/mol. However, unlike the previous theoretical predictions, it is found that the conversion of the [4+2] product to the tetra-sigma-bonded structures entails huge barriers (>37.0 kcal/mol) and is unlikely to occur. This suggests that the [4+2] product is not only thermodynamically the most stable configuration (lowest energy product) but also kinetically very stable (large barriers with respect to the isomerization to other products).

Journal Article↗

PubMed2005

Where does the planar-to-nonplanar turnover occur in small gold clusters?

Several levels of theory, including both Gaussian-based and plane wave density functional theory (DFT), second-order perturbation theory (MP2), and coupled cluster methods (CCSD(T)), are employed to study Au6 and Au8 clusters. All methods predict that the lowest energy isomer of Au6 is planar. For Au8, both DFT methods predict that the two lowest isomers are planar. In contrast, both MP2 and CCSD(T) predict the lowest Au8 isomers to be nonplanar.

Journal Article↗

PubMed2005

Multireference second-order perturbation theory: how size consistent is "almost size consistent".

A systematic study of the deviation from size consistency of the multireference second-order Moller-Plesset perturbation theory (MRMP2) method is presented. The size-consistency error is shown to depend on the number of monomers in a supermolecule calculation, size of basis set, number of correlated valence electrons, and size of active space. HF, F(2), and N(2) are used as test cases, with stretched bonds, to include simple, well-defined multireference character. This is essential in ensuring that MRMP2 is being tested as a multireference method. It is concluded that the MRMP2 and other multireference perturbation theory methods can exhibit significant size-consistency errors, and that the size of the error depends on the manner in which the perturbation theory is implemented.

Journal Article↗

PubMed2004

Fast fragments: the development of a parallel effective fragment potential method.

The Effective Fragment Potential (EFP) method for solvation decreases the cost of a fully quantum mechanical calculation by dividing a chemical system into an ab initio region that contains the solute plus some number of solvent molecules, if desired, and an "effective fragment" region that contains the remaining solvent molecules. Interactions introduced with this fragment region (for example, Coulomb and polarization interactions) are added as one-electron terms to the total system Hamiltonian. As larger systems and dynamics are just starting to be studied with the EFP method, more needs to be done to decrease the calculation time of the method. This article considers parallelization of both the EFP fragment-fragment and mixed quantum mechanics (QM)-EFP interaction energy and gradient computation within the GAMESS suite of programs. The iteratively self-consistent polarization term is treated with a new algorithm that makes use of nonblocking communication to obtain better scalability. Results show that reasonable speedup is achieved with a variety of sizes of water clusters and number of processors.

Journal Article↗

PubMed2004

Binding of propene on small gold clusters and on Au(111): simple rules for binding sites and relative binding energies.

We use density functional theory (DFT) to investigate the bonding of propene to small gas-phase gold clusters and to a Au(111) surface. The desorption energy trends and the geometry of the binding sites are consistent with the following set of rules. (1) The bond of propene to gold is formed by donation of electron density from the highest occupied molecular orbital (HOMO) of propene to one of the low-lying empty orbitals [denoted by LUMO1, LUMO2, em leader (LUMO-lowest unoccupied molecular orbital)] of the gold cluster. (2) Propene binds to a site on the Au cluster where one of the low-lying LUMOs protrudes in the vacuum. Different isomers (same cluster, but different binding sites for propene) correspond to sites where different low-lying LUMOs protrude in space. (3) The desorption energy of the lowest energy isomer correlates with the energy of the lowest empty orbital of the cluster; the lower the energy of that LUMO, the higher the desorption energy. (4) If the lowest-lying LUMO protrudes into space at two nonequivalent sites at the edge of a cluster, propene binds more strongly to the site with the lowest coordination. These rules are consistent with the calculated bond energies and geometries for [Au(n)(C(3)H(6))](q), for n=1-5 and n=8 and q=-1, 0, +1. Based on them we have made a number of predictions that have been confirmed by DFT calculations. The bond of propene to gold is strengthened as the net charge of the cluster varies from -1, to zero, to +1. Compared to a gas-phase cluster, a cluster on a support binds propene more strongly if the support takes electron density from the cluster (e.g., a Au cluster on a gold surface) and more weakly if the support donates electron density to the cluster (e.g., a Au cluster on an oxygen vacancy on an oxide surface).

Alkenes↗

PubMed2004

The effective fragment potential: small clusters and radial distribution functions.

The effective fragment potential (EFP) method for treating solvent effects provides relative energies and structures that are in excellent agreement with the analogous fully quantum [i.e., Hartree-Fock (HF), density functional theory (DFT), and second order perturbation theory (MP2)] results for small water clusters. The ability of the method to predict bulk water properties with a comparable accuracy is assessed by performing EFP molecular dynamics simulations. The resulting radial distribution functions (RDF) suggest that as the underlying quantum method is improved from HF to DFT to MP2, the agreement with the experimental RDF also improves. The MP2-based EFP method yields a RDF that is in excellent agreement with experiment.

Journal Article↗

PubMed2004

Ab initio vibrational state calculations with a quartic force field: applications to H2CO, C2H4, CH3OH, CH3CCH, and C6H6.

For polyatomic molecules, n-mode coupling representations of the quartic force field (nMR-QFF) are presented, which include terms up to n normal coordinate couplings in a fourth-order polynomial potential energy function. The computational scheme to evaluate third-and fourth-order derivatives by finite differentiations of the energy is fully described. The code to generate the nMR-QFF has been implemented into GAMESS program package and interfaced with the vibrational self-consistent field (VSCF) and correlation corrected VSCF (cc-VSCF) methods. As a demonstration, fundamental frequencies have been calculated by the cc-VSCF method based on 2MR-QFF for formaldehyde, ethylene, methanol, propyne, and benzene. The applications show that 2MR-QFF is a highly accurate potential energy function, with errors of 1.0-1.9% relative to the experimental value in fundamental frequencies. This approach will help quantitative evaluations of vibrational energies of a general molecule with a reasonable computational cost.

Journal Article↗

PubMed2004

A new hierarchical parallelization scheme: generalized distributed data interface (GDDI), and an application to the fragment molecular orbital method (FMO).

A two-level hierarchical scheme, generalized distributed data interface (GDDI), implemented into GAMESS is presented. Parallelization is accomplished first at the upper level by assigning computational tasks to groups. Then each group does parallelization at the lower level, by dividing its task into smaller work loads. The types of computations that can be used with this scheme are limited to those for which nearly independent tasks and subtasks can be assigned. Typical examples implemented, tested, and analyzed in this work are numeric derivatives and the fragment molecular orbital method (FMO) that is used to compute large molecules quantum mechanically by dividing them into fragments. Numeric derivatives can be used for algorithms based on them, such as geometry optimizations, saddle-point searches, frequency analyses, etc. This new hierarchical scheme is found to be a flexible tool easily utilizing network topology and delivering excellent performance even on slow networks. In one of the typical tests, on 16 nodes the scalability of GDDI is 1.7 times better than that of the standard parallelization scheme DDI and on 128 nodes GDDI is 93 times faster than DDI (on a multihub Fast Ethernet network). FMO delivered scalability of 80-90% on 128 nodes, depending on the molecular system (water clusters and a protein). A numerical gradient calculation for a water cluster achieved a scalability of 70% on 128 nodes. It is expected that GDDI will become a preferred tool on massively parallel computers for appropriate computational tasks.

Journal Article↗

PubMed2004

A study of the reactions of molecular hydrogen with small gold clusters.

This work presents a study of reactions between neutral and negatively charged Au(n) clusters (n=2,3) and molecular hydrogen. The binding energies of the first and second hydrogen molecule to the gold clusters were determined using density functional theory (DFT), second order perturbation theory (MP2) and coupled cluster (CCSD(T)) methods. It is found that molecular hydrogen easily binds to neutral Au(2) and Au(3) clusters with binding energies of 0.55 eV and 0.71 eV, respectively. The barriers to H(2) dissociation on these clusters with respect to Au(n)H(2) complexes are 1.10 eV and 0.59 eV for n=2 and 3. Although negatively charged Au(n) (-) clusters do not bind molecular hydrogen, H(2) dissociation can occur with energy barriers of 0.93 eV for Au(2) (-) and 1.39 eV for Au(3) (-). The energies of the Au(2)H(2) (-) and Au(3)H(2) (-) complexes with dissociated hydrogen molecules are lower than the energies of Au(2) (-)+H(2) and Au(3) (-)+H(2) by 0.49 eV and 0.96 eV, respectively. There is satisfactory agreement between the DFT and CCSD(T) results for binding energies, but the agreement is not as good for barrier heights.

Journal Article↗