Engineering topics

Boschen, Jeffery S.

Publications and source records attributed to Boschen, Jeffery S..

9 records

DOE OSTI2021

From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape

Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units. This review explores the challenges and some of the solutions in transforming software from the terascale to the petascale and now to the upcoming exascale computers. While discussing the field in general, NWChem and its redesign, NWChemEx, will be highlighted as one of the early codesign projects to take advantage of massively parallel computers and emerging software standards to enable large scientific challenges to be tackled.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Developing a Computational Chemistry Framework for the Exascale Era

Within computational chemistry the NWChem package has arguably been the de facto standard for running high-accuracy numerical simulations on the most powerful supercomputers. In order to better address the challenges presented by emerging exascale architectures, the decision has been made to rewrite NWChem. Design of the resulting package, NWChemEx, has been driven by exascale computing; however, significant additional design considerations have arisen from the team's involvement with the Molecular Sciences Software Institute (MolSSI). MolSSI is a National Science Foundation initiative focused on establishing coding and data standards for the computational chemistry community. As a result, NWChemEx is built upon a general computational chemistry framework called the Simulation Development Environment (SDE) that is designed with a focus on extensibility and interoperability. The present manuscript describes the modular approach of the SDE and how it has been used to implement the self-consistent field algorithm within NWChemEx.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Escape of anions from geminate recombination in THF due to charge delocalization

Geminate recombination of 24 radical anions (M˙ – ) with solvated protons (RH 2 + ) was studied in tetrahydrofuran (THF) with pulse radiolysis. The recombination has two steps: (1) diffusion of M˙ – and RH 2 + together to form intimate (contact and solvent separated) ion pairs, driven by Coulomb attraction; (2) annihilation of anions due to proton transfer (PT) from RH 2 + to M˙ – . The non-exponential time-dependence of the geminate diffusion was determined. For all molecules protonated on O or N atoms the subsequent PT step is too fast (<0.2 ns) to measure, except for the anion of TCNE which did not undergo proton transfer. PT to C atoms was as slow as 70 ns and was always slow enough to be observable. A possible effect of charge delocalization on the PT rates could not be clearly separated from other factors. For 21 of the 24 molecules studied here, a free ion yield (71.6 ± 6.2 nmol J –1 ) comprising ~29% of the total, was formed. This yield of “Type I” free ions is independent of the PT rate because it arises entirely by escape from the initial distribution of ion pair distances without forming intimate ion pairs. Furthermore, three anions of oligo(9,9-dihexyl)fluorenes, F n ˙ – (n = 2–4) were able to escape from intimate ion-pairs to form additional yields of “Type II” free ions with escape rate constants near 3 × 10 6 s –1 . These experiments find no evidence for an inverted region for proton transfer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Stabilization of X–Au–X complexes on the Au(111) surface: A theoretical investigation and comparison of X = S, Cl, CH 3 S, and SiH 3 S

Alnico alloys have long been used as strong permanent magnets because of their ferromagnetism and high coercivity. Understanding their structural details allows for better prediction of the resulting magnetic properties. However, quantitative three-dimensional characterization of the phase separation in these alloys is still challenged by the spatial quantification of nanoscale phases. Herein, we apply a dual tomography approach, where correlative scanning transmission electron microscopy (STEM) energy-dispersive X-ray spectroscopic (EDS) tomography and atom probe tomography (APT) are used to investigate the initial phase separation process of an alnico 8 alloy upon non-magnetic annealing. STEM-EDS tomography provides information on the morphology and volume fractions of Fe–Co-rich and Νi–Al-rich phases after spinodal decomposition in addition to quantitative information of the composition of a nanoscale volume. Subsequent analysis of a portion of the same specimen by APT offers quantitative chemical information of each phase at the sub-nanometer scale. Furthermore, APT reveals small, 2–4 nm Fe-rich α 1 phases that are nucleated in the Ni-rich α 2 matrix. From this information, we show that phase separation of the alnico 8 alloy consists of both spinodal decomposition and nucleation and growth processes. The complementary benefits and challenges associated with correlative STEM-EDS and APT are discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Correlation energy extrapolation by many-body expansion

Accounting for electron correlation is required for high accuracy calculations of molecular energies. The full configuration interaction (CI) approach can fully capture the electron correlation within a given basis, but it does so at a computational expense that is impractical for all but the smallest chemical systems. In this work, a new methodology is presented to approximate configuration interaction calculations at a reduced computational expense and memory requirement, namely, the correlation energy extrapolation by many-body expansion (CEEMBE). This method combines a MBE approximation of the CI energy with an extrapolated correction obtained from CI calculations using subsets of the virtual orbitals. The extrapolation approach is inspired by, and analogous to, the method of correlation energy extrapolation by intrinsic scaling. Benchmark calculations of the new method are performed on diatomic fluorine and ozone. Finally, the method consistently achieves agreement with CI calculations to within a few mhartree and often achieves agreement to within ~1 millihartree or less, while requiring significantly less computational resources.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Comparison of S-adsorption on (111) and (100) facets of Cu nanoclusters

In order to gain insight into the nature of chemical bonding of sulfur atoms on coinage metal surfaces, we compare the adsorption energy and structural parameters for sulfur at four-fold hollow (4fh) sites on (100) facets and at three-fold hollow (3fh) sites on (111) facets of Cu nanoclusters. Consistent results are obtained from localized atomic orbital and plane-wave based density functional theory using the same functionals. PBE and its hybrid counterpart (PBE0 or HSE06) also give similar results. 4fh sites are preferred over 3fh sites with stronger bonding by ~0.6 eV for nanocluster sizes above ~280 atoms. However, for smaller sizes there are strong variations in the binding strength and the extent of the binding site preference. In addition, we show that suitable averaging over clusters of different sizes, or smearing the occupancy of orbitals, provide useful strategies to aid assessment of the behavior in extended surface systems. From site-projected density of states analysis using the smearing technique, we show that S adsorbed on a 4fh site has similar bonding interactions with the substrate as that on a 3fh site, but with much weaker antibonding interactions.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Size Dependence of S-bonding on (111) Facets of Cu Nanoclusters

We demonstrate a strong damped oscillatory size dependence of the adsorption energy for sulfur on the (111) facets of tetrahedral Cu nanoclusters up to sizes of ~300 atoms. This behavior reflects quantum size effects. Consistent results are obtained from density functional theory analyses utilizing either atomic orbital or plane-wave bases and using the same Perdew–Burke–Ernzerhof functional. Behavior is interpreted via molecular orbitals (MO), density of states (DOS), and crystal orbital Hamilton population (COHP) analyses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Mild partial deoxygenation of esters catalyzed by an oxazolinylborate-coordinated rhodium silylene

An electrophilic, coordinatively unsaturated rhodium complex supported by borate-linked oxazoline, oxazoline-coordinated silylene, and N-heterocyclic carbene donors [{κ³-N,Si,C-PhB(Ox Me² )(Ox Me² SiHPh)Im Mes }Rh(H)CO][HB(C₆F₅)₃] (2, Ox Me² = 4,4-dimethyl-2-oxazoline; Im Mes = 1-mesitylimidazole) is synthesized from the neutral rhodium silyl {PhB(Ox Me² )₂Im Mes }RhH(SiH 2 Ph)CO (1) and B(C 6 F 5 ) 3 . The unusual oxazoline-coordinated silylene structure in 2 is proposed to form by rearrangement of an unobserved isomeric cationic rhodium silylene species [{PhB(Ox Me² )₂Im Mes }RhH(SiHPh)CO][HB(C₆F₅)₃] generated by H abstraction. Complex 2 catalyzes reductions of organic carbonyl compounds with silanes to give hydrosilylation products or deoxygenation products. The pathway to these reactions is primarily influenced by the degree of substitution of the organosilane. Reactions with primary silanes give deoxygenation of esters to ethers, amides to amines, and ketones and aldehydes to hydrocarbons, whereas tertiary silanes react to give 1,2-hydrosilylation of the carbonyl functionality. In contrast, the strong Lewis acid B(C₆F₅)₃ catalyzes the complete deoxygenation of carbonyl compounds to hydrocarbons with PhSiH₃ as the reducing agent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Accurate ab initio potential energy curves and spectroscopic properties of the four lowest singlet states of C 2

The diatomic carbon molecule has a complex electronic structure with a large number of low-lying electronic excited states. In this work, the potential energy curves (PECs) of the four lowest lying singlet states (X-1 Sigma(+)(g), A(1)Pi(u), B-1 Delta(g), and B'(1)Sigma(+)(g)) were obtained by high-level ab initio calculations. Valence electron correlation was accounted for by the correlation energy extrapolation by intrinsic scaling (CEEIS) method. Additional corrections to the PECs included core-valence correlation and relativistic effects. Spin-orbit corrections were found to be insignificant. The impact of using dynamically weighted reference wave functions in conjunction with CEEIS was examined and found to give indistinguishable results from the even weighted method. The PECs showed multiple curve crossings due to the B-1 Delta(g) state as well as an avoided crossing between the two (1)Sigma(+)(g) states. Vibrational energy levels were computed for each of the four electronic states, as well as rotational constants and spectroscopic parameters. Comparison between the theoretical and experimental results showed excellent agreement overall. Equilibrium bond distances are reproduced to within 0.05 %. The dissociation energies of the states agree with experiment to within similar to 0.5 kcal/mol, achieving "chemical accuracy." Vibrational energy levels show average deviations of similar to 20 cm(-1) or less. The B-1 Delta(g) state shows the best agreement with a mean absolute deviation of 2.41 cm(-1). Calculated rotational constants exhibit very good agreement with experiment, as do the spectroscopic constants.

36 MATERIALS SCIENCE↗