Search PubMed⌕ Search

Biomedical subjects

Emily A Carter

Publications and source records attributed to Emily A Carter.

13 recordsLinked to original sources
PubMed2006

Insight into selected reactions in low-temperature dimethyl ether combustion from Born-Oppenheimer molecular dynamics.

Dimethyl ether is under consideration as an alternative diesel fuel. Its combustion chemistry is as yet ill-characterized. Here we use Born-Oppenheimer molecular dynamics (BOMD) based on DFT-B3LYP forces to investigate the short-time dynamics of selected features of the low-temperature dimethyl ether (DME) oxidation potential energy surface. Along the chain propagation pathway, we run BOMD simulations from the transition state involving the decomposition of (*)CH(2)OCH(2)OOH to two CH(2)=O and an (*)OH radical. We predict that formaldehyde C-O stretch overtones are excited, consistent with laser photolysis experiments. We also predict that O-H overtones are excited for the (*)OH formed from (*)CH(2)OCH(2)OOH dissociation. We also investigate short-time dynamics involved in chain branching. First, we examine the isomerization transition state of (*)OOCH(2)OCH(2)OOH --> HOOCH(2)OCHOOH. The latter species is predicted to be a short-lived metastable radical that decomposes within 500 fs to hydroperoxymethyl formate (HPMF; HOOCH(2)OC(=O)H) and the first (*)OH of chain branching. The dissociation of HOOCH(2)OCHOOH exhibits non-RRKM behavior in its lifetime profile, which may be due to conformational constraints or slow intramolecular vibrational energy transfer (IVR) from the nascent H-O bond to the opposite end of the radical, where O-O scission occurs to form HPMF and (*)OH. In a few trajectories, we see HOOCH(2)OCHOOH recross back to (*)OOCH(2)OCH(2)OOH because the isomerization is endothermic, with only an 8 kcal/mol barrier to recrossing. Therefore, some inhibition of chain-branching may be due to recrossing. Second, trajectories run from the transition state leading to the direct decomposition of HPMF (an important source of the second (*)OH radical in chain branching) to HCO, (*)OH, and HC(=O)OH show that these products can recombine to form many other possible products. These products include CH(2)OO + HC(=O)OH, H(2)O + CO + HC(=O)OH, HC(=O)OH + HC(=O)OH, and HC(=O)C(=O)H + H(2)O, which (save CH(2)OO + HC(=O)OH) are all more thermodynamically stable than the original HCO + (*)OH + HC(=O)OH products. Moreover, the multitude of extra products suggest that standard statistical rate theories cannot completely describe the reaction kinetics of significantly oxygenated compounds such as HPMF. These secondary products consume the second (*)OH required for explosive combustion, suggesting an inhibition of DME fuel combustion is likely.

Journal Article↗
PubMed2006

A nanoscale mechanism of fatigue in ionic solids.

We employ periodic density functional theory to explore the effect of cyclic tensile loading on the behavior of alumina in the threshold region of crack formation. We find evidence for nanoscale fatigue when the alumina lattice is subjected to uniaxial tensile loading and unloading and tensile stresses normal to the applied load. It is possible that such atomic-scale fatigue impacts the durability of ceramics, since the highly ionic bonding requires near-ideal lattice structures in order to maintain cohesive strength.

Journal Article↗
PubMed2006

Local electronic structure around a single Kondo impurity.

The local electronic structure due to the adsorption of Co on Cu111 is studied using an embedded cluster model, in which the crystal background is taken into account via an effective density functional theory (DFT)-based potential. This approach goes beyond the usual single-impurity Anderson model, where the ground state consists of a singly occupied impurity whose moment is compensated by the background conduction electrons. Ab initio correlated wave function calculations for the embedded cluster provide an alternative picture for this ground state and indicate that the quenching of the Co magnetic moment is due to the formation of metal-metal bonds with the Cu substrate. Low-lying excitations are also studied within the embedding model, and our results are discussed in the context of scanning tunneling microscopy experiments (Manoharan, H. C.; Lutz, C. P.; Eigler, D. M. Nature 2000, 403, 512), which find a sharp tunneling resonance localized in the vicinity of the Co adatom.

Adsorption↗
PubMed2005

Atomic origin of hysteresis during cyclic loading of Si due to bond rearrangements at the crack surfaces.

The atomistic origin of fatigue failure in micron-sized silicon devices is not fully understood. Two series of density-functional theory calculations on cubic diamond Si explore the effect of surface bond formation on crack healing in systems which exhibit strong surface reconstruction. Both series introduce a separation between Si(100) layers (i.e., the crack) and allow the ions to relax to their minimum-energy configuration. The initial surface ionic positions are either bulk terminated or 2 x 1 reconstructed. A plot of the energy versus the introduced separation reveals that once the surfaces reconstruct, the crack is no longer able to return to the equilibrium configuration. Rather, the healed crack interface contains defects which places the flawed energy minimum at a finite strain of 3% and an increased energy of 1.13 Jm2 relative to the equilibrium configuration. The irreversible plastic deformation supports the mechanism proposed by Kahn et al. [Science 298 1215 (2002)] that invokes mechanically induced subcritical cracking to explain the delayed onset of failure.

Journal Article↗
PubMed2005

Ultrasoft spin-dependent pseudopotentials.

The use of the spin-dependent pseudopotentials has been shown to markedly enhance the transferability of the commonly used spin-neutral pseudopotential method for the study of the structural and magnetic properties of transition-metal-containing materials. Unfortunately, because the method was based on the rather expensive norm-conserving pseudopotential formalism, the method was limited to the study of fairly small systems. Here we present an extension of the spin-dependent pseudopotential method for the far more computationally advantageous ultrasoft formalism and show that it is very easy to add such a feature to any preexisting computer code. We benchmark our new method by comparing to previously published results and then apply it to the study of several relevant test cases: bulk Ni, Fe, and Co, as well as a Pd atomic wire.

Journal Article↗
PubMed2005

First principles local pseudopotential for silver: towards orbital-free density-functional theory for transition metals.

Orbital-free density-functional theory (OF-DFT) with modern kinetic-energy density functionals (KEDFs) is a linear scaling technique that accurately describes nearly-free-electron-like (main group) metals. In an attempt towards extending OF-DFT to transition metals, here we consider whether OF-DFT can be used effectively to study Ag, a metal with a localized d shell. OF-DFT has two approximations: use of a KEDF and local pseudopotentials (LPSs). This paper reports construction of a reasonably accurate LPS for Ag by means of inversion of the Kohn-Sham (KS) DFT equations in a bulk crystal environment. The accuracy of this LPS is determined within KS-DFT (where the exact noninteracting kinetic energy is employed) by comparing its predictions of bulk properties to those obtained from a conventional (orbital-based) nonlocal pseudopotential (NLPS). We find that the static bulk properties of fcc and hcp Ag predicted within KS-DFT using this LPS compare fairly well to those predicted by an NLPS. With the transferability of the LPS established, we then use this LPS in OF-DFT, where several approximate KEDFs were tested. We find that a combination of the Thomas-Fermi (T(TF)) and von Weizsacker (T(vW)) functionals (T(vW)+0.4T(TF)) produces better densities than those from the linear-response-based Wang-Teter KEDF. However, the equations of state obtained from both KEDFs in OF-DFT contain unacceptably large errors. The lack of accurate KEDFs remains the final barrier to extending OF-DFT to treat transition metals.

Chemistry, Physical↗
PubMed2005

Finding transition states for crystalline solid-solid phase transformations.

We present a method to identify transition states and minimum energy paths for martensitic solid-solid phase transformations, thereby allowing quantification of the activation energies of such transformations. Our approach is a generalization of a previous method for identifying transition states for chemical reactions, namely the climbing image-nudged elastic band algorithm, where here the global deformation of the crystalline lattice (volume and shape fluctuations) becomes the reaction coordinate instead of atomic motion. We also introduce an analogue to the Born-Oppenheimer approximation that allows a decoupling of nuclear motion and lattice deformation, where the nuclear positions along the path are determined variationally according to current deformation state. We then apply this technique to characterize the energetics of elemental lithium phase transformations as a function of applied pressure, where we see a validation of the Born-Oppenheimer-like approximation, small energy barriers (expected for martensitic transformations), and a pronounced pressure dependence of various properties characterizing the phase transitions.

Algorithms↗
PubMed2005

Improving the orbital-free density functional theory description of covalent materials.

The essential challenge in orbital-free density functional theory (OF-DFT) is to construct accurate kinetic energy density functionals (KEDFs) with general applicability (i.e., transferability). During the last decade, several linear-response (LR)-based KEDFs have been proposed. Among them, the Wang-Govind-Carter (WGC) KEDF, containing a density-dependent response kernel, is one of the most accurate that still affords a linear scaling algorithm. For nearly-free-electron-like metals such as Al and its alloys, OF-DFT employing the WGC KEDF produces bulk properties in good agreement with orbital-based Kohn-Sham (KS) DFT predictions. However, when OF-DFT, using the WGC KEDF combined with a recently proposed bulk-derived local pseudopotential (BLPS), was applied to semiconducting and metallic phases of Si, problems arose with convergence of the self-consistent density and energy, leading to poor results. Here we provide evidence that the convergence problem is very likely caused by the use of a truncated Taylor series expansion of the WGC response kernel. Moreover, we show that a defect in the ansatz for the first-order reduced density matrix underlying the LR KEDFs limits the accuracy of these KEDFs. By optimizing the two free parameters involved in the WGC KEDF, the two-body Fermi wave vector mixing parameter gamma and the reference density rho* used in the Taylor expansion, OF-DFT calculations with the BLPS can achieve semiquantitative results for nine phases of bulk silicon. These new parameters are recommended whenever the WGC KEDF is used to study nonmetallic systems.

Journal Article↗
PubMed2004

Importance of shear in the bcc-to-hcp transformation in iron.

Iron shows a pressure-induced martensitic phase transformation from the ground state ferromagnetic bcc phase to a nonmagnetic hcp phase at approximately 13 GPa. The exact transformation pressure (TP) and pathway are not known. Here we present a multiscale model containing a quantum-mechanics-based multiwell energy function accounting for the bcc and hcp phases of Fe and a construction of kinematically compatible and equilibrated mixed phases. This model suggests that shear stresses have a significant influence on the bcc<-->hcp transformation. In particular, the presence of modest shear accounts for the scatter in measured TPs. The formation of mixed phases also provides an explanation for the observed hysteresis in TP.

Journal Article↗
PubMed2004

First principles resonance widths for Li near an Al(001) surface: predictions of scattered ion neutralization probabilities.

By combining a first-principles periodic density functional theory calculation of adsorbate resonance widths with a many-body dynamical theory of charge transfer, we assess charge-transfer rates for ions scattering off metal surfaces. This goes beyond previous approaches, which have been limited to modeling the surfaces with either static potentials or finite clusters. Here we consider Li(+) scattering from an Al(001) surface. We show how the Li 2s orbital hybridizes with metal valence bands, near the surface, increasing the width of the 2s energy level. This in turn affects the charge-transfer rates between the ion and the metal surface. Our predictions for Li(+)-Al(001) scattering yield the correct angular dependence of the fraction of neutral Li atoms formed when compared to experiment.

Aluminum↗
PubMed2004

Oligoacenes: theoretical prediction of open-shell singlet diradical ground states.

A series of oligoacenes from benzene to decacene were studied computationally with DFT and CASSCF methods. In contrast to the common view that acenes are closed-shell systems or may have a triplet ground state, these results offer the first theoretical predictions for the singlet ground state and diradical character for oligoacenes. The nature of the ground states of these molecules arises from the disjoint nature of the NBMOs that are singly occupied in the diradical.

Journal Article↗
PubMed2004

Size extensive modification of local multireference configuration interaction.

We recently developed a reduced scaling multireference configuration interaction (MRCI) method based on local correlation in the internal (occupied) and external (virtual) orbital spaces. This technique can be used, e.g., to predict bond dissociation energies in large molecules with reasonable accuracy. However, the inherent lack of size extensivity of truncated CI is a disadvantage that in principle worsens as the system size grows. Here we implement an a priori size-extensive modification of local MRCI known as the averaged coupled pair functional (ACPF) method. We demonstrate that local MR-ACPF recovers more correlation energy than local MRCI, in keeping with trends observed previously for nonlocal ACPF. We test the size extensivity of local ACPF on noninteracting He atoms and a series of hydrocarbons. Basis set and core correlation effects are explored, as well as bond breaking in a variety of organic molecules. The local MR-ACPF method proves to be a useful tool for investigating large molecules and represents a further improvement in predictive accuracy over local MRCI.

Journal Article↗