Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel computational chemistry”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Extending atomistic simulation timescale in solid/liquid systems: Crystal growth from solution by a parallel-replica dynamics and continuum hybrid method

Deposition of solid material from solution is ubiquitous in nature. However, due to the inherent complexity of such systems, this process is comparatively much less understood than deposition from a gas or vacuum. Further, the accurate atomistic modeling of such systems is computationally expensive, therefore leaving many intriguing long-timescale phenomena out of reach. We present an atomistic/continuum hybrid method for extending the simulation timescales of dynamics at solid/liquid interfaces. We demonstrate the method by simulating the deposition of Ag on Ag (001) from solution with a significant speedup over standard MD. The results reveal specific features of diffusive deposition dynamics, such as a dramatic increase in the roughness of the film.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Ultrafast infrared transient absorption spectroscopy of gas-phase Ni(CO) 4 photodissociation at 261 nm

Here, we employ ultrafast mid-infrared transient absorption spectroscopy to probe the rapid loss of carbonyl ligands from gas-phase nickel tetracarbonyl following ultraviolet photoexcitation at 261 nm. Here, nickel tetracarbonyl undergoes prompt dissociation to produce nickel tricarbonyl in a singlet excited state; this electronically excited tricarbonyl loses another CO group over tens of picoseconds. Our results also suggest the presence of a parallel, concerted dissociation mechanism to produce nickel dicarbonyl in a triplet excited state, which likely dissociates to nickel monocarbonyl. Mechanisms for the formation of these photoproducts in multiple electronic excited states are theoretically predicted with one-dimensional cuts through the potential energy surfaces and computation of spin–orbit coupling constants using equation of motion coupled cluster methods (EOM-CC) and coupled cluster theory with single and double excitations (CCSD). Bond dissociation energies are calculated with CCSD, and anharmonic frequencies of ground and excited state species are computed using density functional theory (DFT) and time-dependent density functional theory (TD-DFT).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

A hydrogen bonded Co(II) coordination complex and a triply interpenetrating 3-D manganese(II) coordination polymer from diaza crown ether with dibenzoate sidearms

In this paper, the diaza crown ether dicarboxylate ligand, 4,4'-((1,7-dioxa-4,10-diazacyclododecane-4,10-diyl)bis(methylene)dibenzoate), L, forms a monomeric coordination complex of Co(II) ions, CoL(H 2 O) 2 ·2H 2 O, 1, and a coordination polymer (MnL·H 2 O) n , 2, with Mn(II) ions under hydrothermal conditions. The monomeric coordination complex (structure 1) is polar with mirror symmetry and crystallizes in the non-centrosymmetric space group Cm. The mirror plane bisects the complex at the six-coordinate Co(II) ion, the two oxygen atoms of the crown moiety and the two oxygen atoms from coordinated water molecules. The coordinated water molecules take part in strong, linear hydrogen bonds with carboxylate oxygens provided by neighboring Co(II)–crown complexes resulting in a three-dimensional (1 : 1) n polar network in which the topology of the underlying 8-coordinated net is sqc3. The structure of 2 crystallizes in the orthorhombic Fdd2 space group as an infinite, polar triply interpenetrating three-dimensional network. The eight-coordinate Mn(II) ion coordinates two oxygen and two nitrogen atoms of the crown moiety, as well as single carboxylate O atom from each of two neighboring ligands. In both structures, the ligand assumes a flexed-wing bird shape, with the two benzoate sidearms of the crown moiety locked in a syn orientation. The metal ion elevated above the plane of the diaza-crown oxygen and nitrogen atoms. The diaza-crown moiety with its two benzoate sidearms has the peculiar property of forming an oriented crystal structure where the metal-crown vectors are oriented parallel to each other in the crystal. However, the structures are achiral, but the rigid crystal lattice prevents re-orientation of the structure through inversion. The coordination polymer (MnL·H 2 O) n , 2, is a 3-center uninodal net with ths (ThSi 2 ) topology. Computational results using the density functional theory (DFT) calculations explain why the flexed-wing bird shape is the preferred and most stable ligand conformation for metal ion binding. Finally, both structures demonstrate magnetic properties that are characteristics of the respective non-interacting isolated paramagnetic transition metal ions that are present.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Periodic subsystem density-functional theory

By partitioning the electron density into subsystem contributions, the Frozen Density Embedding (FDE) formulation of subsystem Density Functional Theory (DFT) has recently emerged as a powerful tool for reducing the computational scaling of Kohn–Sham DFT. To date, however, FDE has been employed to molecular systems only. Periodic systems, such as metals, semiconductors, and other crystalline solids have been outside the applicability of FDE, mostly because of the lack of a periodic FDE implementation. To fill this gap, in this work we aim at extending FDE to treat subsystems of molecular and periodic character. This goal is achieved by a dual approach. On one side, the development of a theoretical framework for periodic subsystem DFT. On the other, the realization of the method into a parallel computer code. We find that periodic FDE is capable of reproducing total electron densities and (to a lesser extent) also interaction energies of molecular systems weakly interacting with metallic surfaces. In the pilot calculations considered, we find that FDE fails in those cases where there is appreciable density overlap between the subsystems. Conversely, we find FDE to be in semiquantitative agreement with Kohn–Sham DFT when the inter-subsystem density overlap is low. We also conclude that to make FDE a suitable method for describing molecular adsorption at surfaces, kinetic energy density functionals that go beyond the GGA level must be employed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Cluster perturbation theory X: A parallel implementation of Lagrangian perturbation series for the coupled cluster singles and doubles ground-state energy through fifth order

We describe an efficient implementation of cluster perturbation and Møller-Plesset Lagrangian energy series through fifth order that target the coupled cluster singles and doubles energy utilizing the resolution of the identity approximation. We illustrate the computational performance of the implementation by performing ground state energy calculations on systems with up to 1200 basis functions using a single node and by comparison to conventional CCSD calculations. We further show that our hybrid MPI/OMP parallel implementation that also utilizes graphical processing units can be used to obtain fifth order energies on systems with almost 1200 basis functions with a 90 minute "time to solution" running on Frontier at Oak Ridge National Laboratory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

Liquid-crystal-mediated self-assembly at nanodroplet interfaces

Technological applications of liquid crystals have generally relied on control of molecular orientation at a surface or an interface. Such control has been achieved through topography, chemistry and the adsorption of monolayers or surfactants. The role of the substrate or interface has been to impart order over visible length scales and to confine the liquid crystal in a device. In this work, we report results from a computational study of a liquid-crystal-based system in which the opposite is true: the liquid crystal is used to impart order on the interfacial arrangement of a surfactant. Recent experiments on macroscopic interfaces have hinted that an interfacial coupling between bulk liquid crystal and surfactant can lead to a two-dimensional phase separation of the surfactant at the interface, but have not had the resolution to measure the structure of the resulting phases. To enhance that coupling, we consider the limit of nanodroplets, the interfaces of which are decorated with surfactant molecules that promote local perpendicular orientation of mesogens within the droplet. In the absence of surfactant, mesogens at the interface are all parallel to that interface. As the droplet is cooled, the mesogens undergo a transition from a disordered (isotropic) to an ordered (nematic or smectic) liquid-crystal phase. As this happens, mesogens within the droplet cause a transition of the surfactant at the interface, which forms new ordered nanophases with morphologies dependent on surfactant concentration. Such nanophases are reminiscent of those encountered in block copolymers, and include circular, striped and worm-like patterns.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Lateral Standing of the Pentacene Derivative 5,6,7-Trithiapentacene-13-one on Gold: A Combined STM, DFT, and NC-AFM Study

In this paper, high resolution scanning tunneling microscopy investigation of the pentacene derivative, 5,6,7-trithiapentacene-13-one (TTPO), on close-packed gold surfaces reveals nanoscale structures ranging from molecular chains at low coverage to highly ordered self-assembled molecular monolayers. Robust, thermally stable, and highly inert to photo-oxidation, TTPO is a promising organic semiconductor with potential applications in high-temperature photovoltaic devices. When adsorbed on a gold electrode, experiment and first-principle computation reveal a novel 3-D angular assembly of the TTPO molecules, with the long axis of the molecule parallel to the gold surface, distinctive from previously observed pentacene and pentacene derivative assemblies. Structures assembled are angularly dependent on TTPO molecular interactions while commensurate with the underlying gold substrate, allowing for potential tailoring of π-molecular orbital overlap through tilt-angle control. Finally, rigorous noncontact atomic force microscopy analysis was used to experimentally determine the molecular tilt of the molecule at 5 K, consistent with the calculated and observed structures at room temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

High-Pressure Electrides: A Quantum Chemical Perspective

It has long been assumed that all matter will adopt simple close-packed lattices and become metallic under pressure, in accordance with the Thomas–Fermi–Dirac (TFD) model. However, this model struggles to explain pressure-driven complex structural transitions that have been observed in elements, including sodium, challenging our conventional understanding of compressed matter. Moreover, in stark contrast to the TFD model, first-principles calculations suggest that various elements and compounds become electrides under pressure. Electrides, characterized by concentrations of charge density at interstitial regions, can be thought of as ionic compounds where electrons behave as the anions. Though ambient-pressure molecular electrides have been extensively studied via experiments and computations, high-pressure electrides (HPEs) are not well-understood. The identification and characterization of HPEs have been, to date, based purely on theory, including topological analysis of the electron density and the electron localization function. Here, we review these theoretical analysis tools and suggest guidelines that can be used to classify systems as electrides. Moreover, we describe models used to rationalize the electronic structure of HPEs, drawing parallels with ambient-pressure molecular systems, and encourage the development of experimental techniques that provide evidence for the theoretically calculated charge localization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

A hybrid framework of first principles molecular orbital calculations and a three-dimensional integral equation theory for molecular liquids: Multi-center molecular Ornstein–Zernike self-consistent field approach

In this study, we reported the development of a new quantum mechanics/molecular mechanics (QM/MM)-type framework to describe chemical processes in solution by combining standard molecular-orbital calculations with a three-dimensional formalism of integral equation theory for molecular liquids (multi-center molecular Ornstein–Zernike (MC-MOZ) method). The theoretical procedure is very similar to the 3D-reference interaction site model self-consistent field (RISM-SCF) approach. Since the MC-MOZ method is highly parallelized for computation, the present approach has the potential to be one of the most efficient procedures to treat chemical processes in solution. Benchmark tests to check the validity of this approach were performed for two solute (solute water and formaldehyde) systems and a simple S{sub N}2 reaction (Cl{sup −} + CH{sub 3}Cl → ClCH{sub 3} + Cl{sup −}) in aqueous solution. The results for solute molecular properties and solvation structures obtained by the present approach were in reasonable agreement with those obtained by other hybrid frameworks and experiments. In particular, the results of the proposed approach are in excellent agreements with those of 3D-RISM-SCF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Modeling electrokinetics in ionic liquids: General

Using direct numerical simulations we provide a thorough study on the electrokinetics of ionic liquids. In particular, the modfied Poisson-Nernst-Planck (MPNP) equations are solved to capture the crowding and overscreening effects that are the characteristics of an ionic liquid. For modeling electrokinetic flows in an ionic liquid, the MPNP equations are coupled with the Navier-Stokes equations to study the coupling of ion transport, hydrodynamics, and electrostatic forces. Specifically, we consider the ion transport between two parallel plates, charging dynamics in a 2D straight-walled pore, electro-osmotic ow in a nano-channel, electroconvective instability on a plane ion-selective surface, and electroconvective ow on a curved ion-selective surface. We discuss how the crowding and overscreening effects and their interplay affect the electrokinetic behaviors of ionic liquids in these application problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. Here, we present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis

We report that replica-exchange molecular dynamics (REMD) has been proven to efficiently improve the convergence of free energy perturbation (FEP) calculations involving considerable reorganization of their surrounding. We previously introduced the FEP/(λ,H)-REMD algorithm for ligand binding, in which replicas along the alchemical thermodynamic coupling axis λ were expanded as a series of Hamiltonian boosted replicas along a second axis to form a two-dimensional (2D) replica-exchange exchange map [Jiang, W.; Roux, B., J. Chem. Theory Comput. 2010, 6 (9), 2559-2565]. Aiming to achieve a similar performance at a lower computational cost, we propose here a modified version of this algorithm in which only the end-states along the alchemical axis are augmented by boosted replicas. The reduced FEP/(λ,H)-REMD method with one-dimensional (1D) unbiased alchemical thermodynamic coupling axis λ is implemented on the basis of generic multiple copy algorithm (MCA) module of the biomolecular simulation program NAMD. The flexible MCA framework of NAMD enables a user to design customized replica-exchange patterns through Tcl scripting in the context of a highly parallelized simulation program without touching the source code. Two Hamiltonian tempering boosting scheme were examined with the new algorithm: a first one based on potential energy rescaling of a pre-identified “solute”, and a second one via the introduction of flattening torsional free energy barriers. As two illustrative examples with reliable experiment data, the absolute binding free energies of pxylene and n-butylbenzene to the nonpolar cavity of the L99A mutant of T4 lysozyme were calculated. Lastly, the tests demonstrate that the new protocol efficiently enhances the sampling of torsional motions for backbone and side chains around the binding pocket and accelerates the convergence of the free energy computations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Origin of the magnetic couplings for the weak ferromagnet Li + [TCNE] •- (TCNE = Tetracyanoethylene)

Based upon the 16-K crystal structure, the nearest-neighbor spin couplings (J; H = –2JS a •S b ) for the weak ferromagnet (=canted antiferromagnet) Li + [TCNE] •- (TCNE = tetracyanoethylene) (T c = 21 K), which possesses two interpenetrating diamondoid sublattices with two layers of parallel [TCNE] •- s canted with respect to each other by ~60°, are computed at the B3LYP/aug-cc-pV-TZ-level. Li[TCNE] has computed interlattice, through-space interactions that exceed the intralattice, through-bonding to Li + interactions, that are antiferromagnetic. The interlattice interactions dominate the spin coupling leading to the observed weak ferromagnetic behavior. Finally, the computed J's based upon the 50-K structure have the strongest ferromagnetic intralayer interlattice interactions increase by 2.9 ± 0.6% while the antiferromagnetic interlayer interlattice interaction increases by 23% with respect to the 16K results and are in accord with the observed initial increase and subsequent decrease in the temperature dependence of the canting angle that is a consequence of the these changing opposed interactions arising from the change in structure as a function of temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

MADNESS: A Multiresolution, Adaptive Numerical Environment for Scientific Simulation

MADNESS (multiresolution adaptive numerical environment for scientific simulation) is a high-level software environment for solving integral and differential equations in many dimensions that uses adaptive and fast harmonic analysis methods with guaranteed precision based on multiresolution analysis and separated representations. Underpinning the numerical capabilities is a powerful petascale parallel programming environment that aims to increase both programmer productivity and code scalability. This paper describes the features and capabilities of MADNESS and briefly discusses some current applications in chemistry and several areas of physics.

00A72↗
DOE OSTI2013

Massively parallel implementations of coupled-cluster methods for electron spin resonance spectra. I. Isotropic hyperfine coupling tensors in large radicals

Coupled cluster (CC) methods provide highly accurate predictions of molecular properties, but their high computational cost has precluded their routine application to large systems. Fortunately, recent computational developments in the ACES III program by the Bartlett group [the OED/ERD atomic integral package, the super instruction processor, and the super instruction architecture language] permit overcoming that limitation by providing a framework for massively parallel CC implementations. In that scheme, we are further extending those parallel CC efforts to systematically predict the three main electron spin resonance (ESR) tensors (A-, g-, and D-tensors) to be reported in a series of papers. In this paper inaugurating that series, we report our new ACES III parallel capabilities that calculate isotropic hyperfine coupling constants in 38 neutral, cationic, and anionic radicals that include the {sup 11}B, {sup 17}O, {sup 9}Be, {sup 19}F, {sup 1}H, {sup 13}C, {sup 35}Cl, {sup 33}S,{sup 14}N, {sup 31}P, and {sup 67}Zn nuclei. Present parallel calculations are conducted at the Hartree-Fock (HF), second-order many-body perturbation theory [MBPT(2)], CC singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] levels using Roos augmented double- and triple-zeta atomic natural orbitals basis sets. HF results consistently overestimate isotropic hyperfine coupling constants. However, inclusion of electron correlation effects in the simplest way via MBPT(2) provides significant improvements in the predictions, but not without occasional failures. In contrast, CCSD results are consistently in very good agreement with experimental results. Inclusion of perturbative triples to CCSD via CCSD(T) leads to small improvements in the predictions, which might not compensate for the extra computational effort at a non-iterative N{sup 7}-scaling in CCSD(T). The importance of these accurate computations of isotropic hyperfine coupling constants to elucidate experimental ESR spectra, to interpret spin-density distributions, and to characterize and identify radical species is illustrated with our results from large organic radicals. Those include species relevant for organic chemistry, petroleum industry, and biochemistry, such as the cyclo-hexyl, 1-adamatyl, and Zn-porphycene anion radicals, inter alia.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗