Biomedical subjects

Nathan E Schultz

Publications and source records attributed to Nathan E Schultz.

7 recordsLinked to original sources

PubMed2006

Phase behavior of elemental aluminum using monte carlo simulations.

Monte Carlo simulations are presented for two models of aluminum: an embedded-atom model and an explicit many-body model. Vapor/liquid coexistence curves are determined using Gibbs ensemble Monte Carlo simulations. The normal boiling points predicted by both models are somewhat higher (by about 10%) than the experimental value. Isothermal constant-stress simulations are used to simulate solid Al from 300 K to the triple point. The solid structures are at least metastable in the face-centered cubic configuration, and the specific heat is determined to be lower than the experimental value. The melting point for the embedded-atom model determined via thermodynamic integration along a pseudo-supercritical path is approximately 20% higher than the experimental value.

Journal Article↗

PubMed2006

PdnCO (n = 1,2): accurate Ab initio bond energies, geometries, and dipole moments and the applicability of density functional theory for fuel cell modeling.

Electrode poisoning by CO is a major concern in fuel cells. As interest in applying computational methods to electrochemistry is increasing, it is important to understand the levels of theory required for reliable treatments of metal-CO interactions. In this paper we justify the use of relativistic effective core potentials for the treatment of PdCO and hence, by inference, for metal-CO interactions where the predominant bonding mechanism is charge transfer. We also sort out key issues involving basis sets and we recommend that bond energies of 17.2, 43.3, and 69.4 kcal/mol be used as the benchmark bond energy for dissociation of Pd2 into Pd atoms, PdCO into Pd and CO, and Pd2CO into Pd2 and CO, respectively. We calculated the dipole moments of PdCO and Pd2CO, and we recommend benchmark values of 2.49 and 2.81 D, respectively. Furthermore, we tested 27 density functionals for this system and found that only hybrid density functionals can qualitatively and quantitatively predict the nature of the sigma-donation/pi-back-donation mechanism that is associated with the Pd-CO and Pd2-CO bonds. The most accurate density functionals for the systems tested in this paper are O3LYP, OLYP, PW6B95, and PBEh.

Journal Article↗

PubMed2006

Critical properties of aluminum.

Gibbs ensemble Monte Carlo calculations are performed using a validated embedded-atom potential to obtain the vapor-liquid coexistence curve for elemental aluminum in good agreement with available experimental data up to the boiling point. These calculations are then extended to make a reliable prediction of the critical temperature, pressure, and density of Al, which have previously been known only with very large uncertainties. This demonstrates the ability of modern simulations to predict fundamental physical properties that are extremely difficult to measure directly.

Journal Article↗

PubMed2005

Density functionals for inorganometallic and organometallic chemistry.

We present a database of 21 bond dissociation energies for breaking metal-ligand bonds. The molecules in the metal-ligand bond energy database are AgH, CoH, CoO+, CoOH+, CrCH3+, CuOH2+, FeH, Fe(CO)5, FeO, FeS, LiCl, LiO, MgO, MnCH3NiCH2+, Ni(CO)4, RhC, VCO+, VO, and VS. We have also created databases of metal-ligand bond lengths and atomic ionization potentials. The molecules used for bond lengths are AgH, BeO, CoH, CoO+, FeH, FeO, FeS, LiCl, LiO, MgO, RhC, VO, and VS and the ionization potentials are for the following atoms: C, Co, Cr, Cu, Ni, O, and V. The data were chosen based on their diversity and expected reliability, and they are used along with three previously developed databases (transition metal dimer bond energies and bond lengths and main-group molecular atomization energies) for assessing the accuracy of several kinds of density functionals. In particular, we report tests for 42 previously defined functionals: 2 local spin density approximation (LSDA) functionals, 14 generalized gradient approximation (GGA) methods, 13 hybrid GGA methods, 7 meta GGA methods, and 8 hybrid meta GGA methods. In addition to these functionals, we also examine the effectiveness of scaling the correlation energy by testing 13 functionals with scaled or no gradient-corrected correlation energy, and we find that functionals of this kind are more accurate for metal-metal and metal-ligand bonds than any of the functionals already in the literature. We also present a readjusted GGA and a hybrid GGA with parameters adjusted for metals. When we consider these 57 functionals for metal-ligand and metal-metal bond energies simultaneously with main-group atomization energies, atomic ionization potentials, and bond lengths we find that the most accurate functional is G96LYP, followed closely by MPWLYP1M (new in this article), XLYP, BLYP, and MOHLYP (also new in this article). Four of these five functionals have no Hartree-Fock exchange, and the other has only 5%. As a byproduct of this work we introduce a convenient diagnostic, called the B1 diagnostic, for ascertaining the multireference character in a bond.

Journal Article↗

PubMed2005

Exchange-correlation functional with broad accuracy for metallic and nonmetallic compounds, kinetics, and noncovalent interactions.

By incorporating kinetic-energy density in a balanced way in the exchange and correlational functionals and removing self-correlation effects, we have designed a density functional that is broadly applicable to organometallic, inorganometallic, and nonmetallic bonding, thermochemistry, thermochemical kinetics, and noncovalent interactions as well as satisfying the uniform electron gas limit. The average error is reduced by a factor of 1.3 compared with the best previously available functionals, but even more significantly, we find a functional that has a high accuracy for all four categories of interaction.

Journal Article↗

PubMed2005

Databases for transition element bonding: metal-metal bond energies and bond lengths and their use to test hybrid, hybrid meta, and meta density functionals and generalized gradient approximations.

We propose a data set of bond lengths for 8 selected transition metal dimers (Ag(2), Cr(2), Cu(2), CuAg, Mo(2), Ni(2), V(2), and Zr(2)) and another data set containing their atomization energies and the atomization energy of ZrV, and we use these for testing density functional theory. The molecules chosen for the test sets were selected on the basis of the expected reliability of the data and their ability to constitute a diverse and representative set of transition element bond types while the data sets are kept small enough to allow for efficient testing of a large number of computational methods against a very reliable subset of experimental data. In this paper we test 42 different functionals: 2 local spin density approximation (LSDA) functionals, 12 generalized gradient approximation (GGA) methods, 13 hybrid GGAs, 7 meta GGA methods, and 8 hybrid meta GGAs. We find that GGA density functionals are more accurate for the atomization energies of pure transition metal systems than are their meta, hybrid, or hybrid meta analogues. We find that the errors for atomization energies and bond lengths are not as large if we limit ourselves to dimers with small amounts of multireference character. We also demonstrate the effects of increasing the fraction of Hartree-Fock exchange in multireference systems by computing the potential energy curve for Cr(2) and Mo(2) with several functionals. We also find that BLYP is the most accurate functional for bond energies and is reasonably accurate for bond lengths. The methods that work well for transition metal bonds are found to be quite different from those that work well for organic and other main group chemistry.

Journal Article↗

PubMed2005

Analytic potential energy functions for simulating aluminum nanoparticles.

Potential energy functions (PEFs) parametrized to bulk data are shown to perform poorly for small aluminum nanoparticles and clusters. In contrast, PEFs parametrized to a limited set of cluster and bulk data, but no nanoparticle data, perform well for nanoparticles. This validates a practical scheme for developing PEFs for nanoscale systems. Building on these findings, we optimized five PEFs by minimizing the error in the fit over a broad data set. Two of these PEFs have errors of less than or equal to 0.08 eV/atom for each of three categories of system sizes, i.e., for small clusters, for nanoparticles, and for bulk potential energies.

Journal Article↗