Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic & molecular structure”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Numerical simulation of three-dimensional self-gravitating flow

The three-dimensional flow of a self-gravitating fluid is numerically simulated using a Fourier pseudospectral method with a logarithmic variable formulation. Two cases with zero total angular momentum are studied in detail, a 323 simulation (Run B). Other than the grid size, the primary difference between the two cases are that Run A modeled atomic hydrogen and had considerably more compressible motion initially than Run B, which modeled molecular hydrogen. The numerical results indicate that gravitational collapse can proceed in a variety of ways. In the Run A, collapse led to an elongated tube-like structure, while in the Run B, collapse led to a flatter, disklike structure.

Shebalin, John V.↗
NASA NTRS1994

An Einstein survey of the 1 keV soft X-ray background in the Galactic plane

We have analyzed 56 Einstein Observatory Imaging Proportional Counter (IPC) observations within +/- 3 deg of the Galactic plane in order to determine the low-latitude soft X-ray background flux in the 0.56-1.73 keV band. Any detected X-ray point source which fell within our regions of study was removed from the image, enabling us to present maps of the background flux as a function of Galactic latitude along 18 meridians. These maps reveal considerable structure to the background in the Galactic plane on an angular scale of approximately 1 deg. Our results are compared with those of an earlier study of the 1 keV X-ray background along l = 25 deg by Kahn & Caillault. The double-peaked structure they found is not discernible in our results, possibly because of the presence of solar backscattered flux in their data. A model which takes into account contributions to the background by extragalactic and stellar sources, the distribution of both atomic and molecular absorbing material with the Galaxy, the energy dependence of the cross section for absorption of X-rays, and the energy dependence of the detector has been constructed and fitted to these new data to derive constraints on the scale height, temperature, and volume emissivity of the unaccounted-for X-ray-emitting material. The results of this model along l = 25 deg are roughly similar to those of the model of Kahn & Caillault along the same meridian.

Stanford, John M.↗
DOE OSTI2025

Evidence for In-Plane Electrical Polarization in 3R-β’-In 2 Se 3 Thin Films Grown by Molecular Beam Epitaxy

β-In 2 Se 3 has been identified as a potential ferroelectric. This work describes the growth of β-In 2 Se 3 via molecular beam epitaxy along with a description of its properties, including the search for switchable polarization in 3R-β’-In2Se3. The thin-film morphology, crystal structure, and phase maps of β-In 2 Se 3 on Si(111) and Al 2 O 3 (0001) were evaluated by atomic force microscopy, X-ray diffraction, and Raman spectroscopy as a function of the atomic Se/In flux ratio and growth temperature. Smooth β-In 2 Se 3 thin films were successfully realized on Si(111) at a substrate temperature of 150 °C using a Se/In flux ratio of 5.7, as well as on Al 2 O 3 (0001) at 450 °C using a Se/In flux ratio of 5.5. Scanning transmission electron microscopy (STEM) confirms the 3R polytype of β-In 2 Se 3 in films on both substrates, with a minor disorder associated with the 2H polytype at the interface. Indications of in-plane Se atom displacements characteristic of the 3R-β’-In 2 Se 3 polytype were found by STEM and second harmonic generation analysis in films on Al 2 O 3 (0001), but not in films on Si(111). Furthermore, attempts at electrical polarization switching did not produce compelling evidence for ferroelectricity. Instead, electrical transport measurements demonstrated locally varying anisotropic responses with the applied electric fields along different in-plane directions, with some hysteresis associated with the trapping of charges.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electrochemical and Optical Spectroscopic Probing of Transition–Sized Au 130 (SR) 50 Nanoclusters

Ultrasmall metal nanoclusters (NCs) exhibit a quantized conduction band, hence, a distinct HOMO-LUMO gap (E g ). Such a quantized electronic structure gives rise to multiple discrete peaks in the optical absorption spectrum of the NCs. As the size grows to 130 gold atoms (Au 130 protected by ligands), electrical charging and optical behaviors seem to show certain metal-like features (hence, transition-sizes). To probe such behaviors, especially the potential ligand effect, we have devised the synthesis of Au 130 NCs protected by phenylethanethiolate and naphthalenethiolate, respectively, with the former having a nonconjugated separation between the aromatic molecular group and the metal core while the latter being in direct bonding. A careful comparison of these two Au 130 nanoclusters with the earlier reported analogues is carried out, including the structurally characterized Au 130 (pMBT) 50 and the aqueous counterpart. While all of these nanoclusters possess the same 80 free electron counts in the core, some notable differences in electrochemical and optical properties are found, which are attributed to the ligand effects. The obtained insights may stimulate further interest in the transition-sized nanoclusters and also promote their applications in optics, energy conversion, and biomedicine.

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

Materials Data on Ti2H2CSe4Cl16O5 by Materials Project

Ti2C(OCl2)3(HSeOCl)2(SeCl3)2Cl2 crystallizes in the monoclinic P2/c space group. The structure is two-dimensional and consists of two molecular chlorine molecules; four HSeOCl clusters; four SeCl3 clusters; and one Ti2C(OCl2)3 sheet oriented in the (0, 1, 0) direction. In each HSeOCl cluster, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. Se3+ is bonded in a distorted L-shaped geometry to one O2- and one Cl1- atom. The Se–O bond length is 1.83 Å. The Se–Cl bond length is 2.50 Å. O2- is bonded in a distorted water-like geometry to one H1+ and one Se3+ atom. Cl1- is bonded in a distorted single-bond geometry to one Se3+ atom. In each SeCl3 cluster, Se3+ is bonded in a single-bond geometry to one Cl1- atom. The Se–Cl bond length is 2.76 Å. There are three inequivalent Cl1- sites. In the first Cl1- site, Cl1- is bonded in a 1-coordinate geometry to one Cl1- atom. The Cl–Cl bond length is 2.30 Å. In the second Cl1- site, Cl1- is bonded in a 2-coordinate geometry to one Se3+ and one Cl1- atom. The Cl–Cl bond length is 2.32 Å. In the third Cl1- site, Cl1- is bonded in a distorted linear geometry to two Cl1- atoms. In the Ti2C(OCl2)3 sheet, Ti4+ is bonded in a 5-coordinate geometry to one O2- and four Cl1- atoms. The Ti–O bond length is 2.71 Å. There are a spread of Ti–Cl bond distances ranging from 2.15–2.47 Å. C4+ is bonded in a linear geometry to two equivalent O2- atoms. Both C–O bond lengths are 1.18 Å. There are two inequivalent O2- sites. In the first O2- site, O2- is bonded in a single-bond geometry to one Ti4+ and one C4+ atom. In the second O2- site, O2- is bonded in a bent 150 degrees geometry to two equivalent Cl1- atoms. Both O–Cl bond lengths are 1.93 Å. There are three inequivalent Cl1- sites. In the first Cl1- site, Cl1- is bonded in a single-bond geometry to one Ti4+ atom. In the second Cl1- site, Cl1- is bonded in a water-like geometry to one Ti4+ and one O2- atom. In the third Cl1- site, Cl1- is bonded in an L-shaped geometry to two equivalent Ti4+ atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on H12Cl2O13 by Materials Project

H7O3H5O2(ClO4)2 crystallizes in the monoclinic P2_1/c space group. The structure is zero-dimensional and consists of four molecular hydrogen;dihydrate molecules, eight ClO4 clusters, and four H7O3 clusters. In each ClO4 cluster, there are four inequivalent O sites. In the first O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.48 Å. In the second O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.46 Å. In the third O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.46 Å. In the fourth O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.46 Å. Cl is bonded in a tetrahedral geometry to four O atoms. In each H7O3 cluster, there are seven inequivalent H sites. In the first H site, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 1.02 Å. In the second H site, H is bonded in a linear geometry to two O atoms. There is one shorter (1.05 Å) and one longer (1.53 Å) H–O bond length. In the third H site, H is bonded in a linear geometry to two O atoms. There is one shorter (1.05 Å) and one longer (1.53 Å) H–O bond length. In the fourth H site, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 1.00 Å. In the fifth H site, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 1.00 Å. In the sixth H site, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 0.99 Å. In the seventh H site, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 0.98 Å. There are three inequivalent O sites. In the first O site, O is bonded in a 3-coordinate geometry to three H atoms. In the second O site, O is bonded in a trigonal non-coplanar geometry to three H atoms. In the third O site, O is bonded in a 3-coordinate geometry to three H atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nickel(ii) complexes based on dithiolate–polyamine binary ligand systems: crystal structures, hirshfeld surface analysis, theoretical study, and catalytic activity study on photocatalytic hydrogen generation

To ascertain the influence of binary ligand systems [1,1-dicyanoethylene-2,2-dithiolate (i-mnt -2 ) and polyamine {tetraen = tris(2-aminoethyl)amine, tren = diethylene triamine and opda = o-phenylenediamine}] on the coordination modes of the Ni(II) metal center and resulting supramolecular architectures, a series of nickel(II) thiolate complexes [Ni(tetraen)(i-mnt)](DMSO) (1), [Ni 2 (tren) 2 (i-mnt) 2 ] (2), and [Ni 2 (i-mnt) 2 (opda) 2 ] n (3) have been synthesized in high yield in one step in water and structurally characterized by single crystal X-ray crystallography and spectroscopic techniques. X-ray diffraction studies disclose the diverse i-mnt -2 coordination to the Ni +2 center in the presence of active polyamine ligands, forming a slightly distorted octahedral geometry (NiN 4 S 2 ) in 1, square planar (NiS 4 ) and distorted octahedral geometries (NiN 6 ) in the bimetallic co-crystallized aggregate of cationic [Ni(tren) 2 ] +2 and anionic [Ni(i-mnt) 2 ] -2 in 2, and a one dimensional (1D) polymeric chain along the [100] axis in 3, having consecutive square planar (NiS 4 ) and octahedral (NiN 6 ) coordination kernels. The N–H···O, N–H···S, N–H···N, N–H···S, N–H···N, and N–H···O type hydrogen bonds stabilize the supramolecular assemblies in 1, 2, and 3 respectively imparting interesting graph-set-motifs. The molecular Hirshfeld surface analyses (HS) and 2D fingerprint plots were utilized for decoding all types of non-covalent contacts in the crystal networks. Atomic HS analysis of the Ni +2 centers reveals significant Ni–N metal–ligand interactions compared to Ni–S interactions. We have also studied the unorthodox interactions observed in the solid state structures of 1–3 by QTAIM and NBO analyses. Moreover, all the complexes proved to be highly active water reduction co-catalysts (WRC) in a photo-catalytic hydrogen evolution process involving iridium photosensitizers, wherein 2 and 3 having a square planar arrangement around the nickel center(s) – were found to be the most active ones, achieving 1000 and 1119 turnover numbers (TON), respectively.

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

Vibronic Excitons and Conical Intersections in Semiconductor Quantum Dots

Surface defects and organic surface-capping ligands affect the photoluminescence properties of semiconductor quantum dots (QDs) by altering the rates of competing nonradiative relaxation processes. In this study, broadband two-dimensional electronic spectroscopy reveals that absorption of light by QDs prepares vibronic excitons, excited states derived from quantum coherent mixing of the core electronic and ligand vibrational states. Rapidly damped coherent wavepacket motions of the ligands are observed during hot-carrier cooling, with vibronic coherence transferred to the photoluminescent state. Furthermore, these findings suggest a many-electron, molecular theory for the electronic structure of QDs, which is supported by calculations of the structures of conical intersections between the exciton potential surfaces of a small ammonia-passivated model CdSe nanoparticle.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Novel synthesis approach for “stubborn” metals and metal oxides

Significance Atomically precise complex oxides containing “stubborn” elements, such as ruthenium, iridium, and platinum, hold tremendous promise as designer quantum materials for exploring novel electronic, magnetic, superconducting, and topological phases owing to their strong spin–orbit interaction. This study shows a method to synthesize such materials by eliminating the major synthesis bottleneck of low vapor pressure and difficulty in oxidation. This study serves as a “proof-of-concept” allowing us to 1) grow Pt, RuO 2 , and SrRuO 3 thin films by supplying Pt and Ru precursors at 65 to 100 °C in a low-temperature effusion cell, as opposed to the several thousand degrees Celsius needed using electron beam evaporators; 2) reveal bulk-like room-temperature resistivity; and 3) ultimately provide pathways to creating atomically precise quantum structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spectroscopy of Rydberg states in erbium using electromagnetically induced transparency

We present a study of the Rydberg spectrum in 166 Er for series connected to the 4 f 12 ( 3 H 6 )6s, J c = 13/2 and J c = 11/2 ionic core states using an all-optical detection based on electromagnetically induced transparency in an effusive atomic beam. Identifying approximately 550 individual states, we find good agreement with a multi-channel quantum defect theory (MQDT) which allows assignment of most states to ns or nd Rydberg series. We provide an improved accuracy for the lowest two ionization thresholds to E IP, J c =13/2 = 49260:750(1) cm -1 and E IP, J c =11/2 = 49701:184(1) cm -1 as well as the corresponding quantum defects for all observed series. We identify Rydberg states in five different isotopes, and states between the two lowest ionization thresholds. Our results open the way for future applications of Rydberg states for quantum simulation using erbium and exploiting its special open-shell structure.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2022

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
DOE OSTI2020

Materials Data on H7ClO3 by Materials Project

H5O2H7O3ClH2OCl crystallizes in the triclinic P1 space group. The structure is zero-dimensional and consists of one hydrochloric acid, monohydrate molecule; one molecular hydrogen;dihydrate molecule; and one H7O3Cl cluster. In the H7O3Cl cluster, there are seven inequivalent H1+ sites. In the first H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 1.00 Å. In the second H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. In the third H1+ site, H1+ is bonded in a distorted single-bond geometry to one O2- and one Cl1- atom. The H–O bond length is 1.03 Å. The H–Cl bond length is 1.88 Å. In the fourth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.98 Å. In the fifth H1+ site, H1+ is bonded in a distorted linear geometry to two O2- atoms. There is one shorter (1.02 Å) and one longer (1.61 Å) H–O bond length. In the sixth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 1.00 Å. In the seventh H1+ site, H1+ is bonded in a linear geometry to two O2- atoms. There is one shorter (1.14 Å) and one longer (1.29 Å) H–O bond length. There are three inequivalent O2- sites. In the first O2- site, O2- is bonded in a trigonal non-coplanar geometry to three H1+ atoms. In the second O2- site, O2- is bonded in a distorted water-like geometry to three H1+ atoms. In the third O2- site, O2- is bonded in a trigonal non-coplanar geometry to three H1+ atoms. Cl1- is bonded in a single-bond geometry to one H1+ atom.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unfolding of the Villin Headpiece Domain: Revealing Structural Heterogeneity with Time‐Resolved X‐Ray Solution Scattering and Markov State Modeling

Understanding protein folding pathways is crucial to deciphering the principles of protein structure and function. Here, the unfolding dynamics of the 35‐residue villin headpiece (HP35) and a norleucine‐substituted variant (2F4K) using a combination of experimental and computational techniques is investigated. Time‐resolved X‐ray solution scattering coupled with equilibrium molecular dynamics simulations and Markov state modeling reveals distinct unfolding mechanisms between the two variants: HP35 and 2F4K. Specifically, HP35 exhibits a two‐state unfolding process, whereas an intermediate state is identified for the 2F4K mutant. A Markov state model constructed from simulations is used to map atomic‐level transitions to experimental observations, providing insights into the role of sequence variations in modulating folding pathways. The findings underscore the importance of integrating experimental and computational approaches to unravel protein unfolding mechanisms between heterogenous structural ensembles.

Nijhawan, Adam K. [Department of Chemistry Northwe↗
NASA NTRS2018

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗
DOE OSTI2026

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗
DOE OSTI2025

Structure and ultrafast dynamics of tri-nuclear Ag-/Tl–Pt 2 POP 4 complexes in solution

The energetics and dynamics of ion assembly in solution has broad influence in nanomaterials and inorganic synthesis. To investigate the fundamental processes involved, we present a time-resolved x-ray solution scattering (TR-XSS) study of the trinuclear silver and thallium complexes of the diplatinum ion PtPOP [Pt 2 (H 2 P 2 O 5 )$_4^{4−}$] in aqueous solution. These complexes, their structural properties, and their electronic structure are not well understood and afford a unique opportunity to study the metal–metal bond formation that influences molecular and material assembly in solution. We present model-independent analysis of the observed dynamics as well as an analysis incorporating time-resolved structural refinements of key bond lengths with $<$100 fs time resolution. We find that upon photoexcitation, the Pt atoms contract ∼0.25 Å toward the center of both the Ag- and the Tl-PtPOP complexes, as previously observed for the PtPOP anion. For the AgPtPOP system, an ultrafast Ag-Pt bond expansion of ∼0.2 Å is observed, whereas in contrast, the TlPtPOP system exhibits a Tl-Pt bond contraction of ∼0.3 Å upon photoexcitation. For both complexes, the change in electronic state leads to coherent (“wave-packet”) oscillations along the metal–Pt coordinates. Based on these structural dynamics, we propose an electronic structure model that describes the metal–metal bonding behavior in both the ground and excited state for both complexes.

Lenzen, Philipp [Technical Univ. of Denmark, Lyngb↗
DOE OSTI2026

On the Prospect of Chemically Transferable Coarse-Grained Electronic Models for Soft Materials

Electronic coarse-graining (ECG) methods predict quantum-mechanical electronic properties directly from coarse-grained (CG) molecular configurations, enabling electronic predictions at mesoscale length scales. Here, we present a diagnostic assessment of the feasibility of chemically transferable ECG models across a broad polymer-relevant chemical space using all-atom, united-atom, and Martini-scale representations. While high-resolution ECG models achieve near-quantitative accuracy, we show that chemically transferable ECG at the Martini resolution fails because the CG force field does not sample the same configurational distribution of local molecular structure as that underlying the DFT-parameterized ECG model. We demonstrate that our proposed Element-Count-Label (ECL) representation, which augments Martini beads with explicit stoichiometric data, significantly improves chemical generalization across diverse polymer chemistries. However, we find that even with improved chemical resolution, the model cannot recover electronic property distributions that are absent from the configurational space sampled by the CG force field. These results demonstrate that chemically transferable ECG requires future Martini-like force fields to explicitly preserve quantum chemistry–compatible local molecular structure in addition to thermodynamic and structural fidelity.

Kidder, Katherine M [Department of Chemistry; Univ↗
NASA NTRS1988

Recent advances in electronic structure theory and their influence on the accuracy of ab initio potential energy surfaces

Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F+H2 yields HF+H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces.

Bauschlicher, Charles W., Jr.↗