Search NASA⌕ Search

SEARCH · Search NASA

Results for “EXTRAPOLATION”

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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Density Measurements of Molten LiF–BeF 2 and LiF–BeF 2 –LaF 3 Salt Mixtures by Neutron Radiography

The densities of eutectic (LiF) 2 –BeF 2 and mixtures of this salt (FLiBe) with LaF 3 were measured by dilatometry and by neutron attenuation from 673 K to 1,073 K. Because LaF 3 has a limited solubility in FLiBe, it was necessary to determine the amount of LaF 3 in solution before the density could be determined. The FLiBe density determination was favorably benchmarked against the literature data. A simple comparison was not available for the LaF 3 –FLiBe mixtures, so extrapolation of published data was necessary based on analysis using the Molten Salt Thermal Properties Database-Thermochemistry, or MSTDB-TC, developed by the US Department of Energy. Solubilities for LaF 3 in FLiBe ranged from 1 to 4 mol % over 673 to 1,073 K. The salt system was heated and cooled over 24 h to evaluate potential changes in composition and hysteresis during the measurement. Changes in the meniscus were observed, and these were included in the correction for density determinations. Salt surface tension may have led to supersaturation of LaF 3 in the salt because the solubility curve was nonlinear with respect to the inverse temperature, as would be expected for an ideal system. Surface tension measurements are currently underway to test this hypothesis.

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

Henry’s Solubility and Diffusion Coefficients for 29 Volatile Organic Compounds in Polydimethylsiloxane Sylgard 184 at 293 K

Two-dimensional (2D) inverse gas chromatography (IGC) enables simultaneous determination of Henry’s solubility and Fickian diffusion coefficients for volatile organic compounds (VOCs) in polymer films. This technique offers a significant advantage over traditional cylindrical column IGC by providing precise control and measurement of the film thickness (here, 0.064 ± 0.002 mm), which is the critical length scale for accurate diffusivity determination. We apply this methodology to characterize VOC transport in Sylgard 184, a widely used polydimethylsiloxane (PDMS)-based polymer containing substantial silica filler content. At room temperature (20 °C), we measured solubility and diffusion coefficients for 29 common VOCs spanning diverse chemical functionalities, including alkanes, aromatics, chlorinated solvents, ketones, esters, and alcohols. Comparison with literature data for pure PDMS reveals that VOC solubility in Sylgard 184 is generally higher; for most non-hydrogen-bonding compounds it remains within a factor of 2 of pure PDMS, whereas alcohols are enhanced by roughly 1.8 to 3.7 times, which we attribute to favorable interactions with residual silanol groups on the silanized silica filler. Diffusion coefficients range from 1.0 × 10 –6 cm 2 /s (n-undecane) to 8.9 × 10–5 cm 2 /s (acetonitrile) and align well with extrapolated literature values for PDMS. This study provides essential thermodynamic and transport data for predicting VOC permeation in Sylgard 184 while demonstrating the utility of 2D IGC as a robust technique for characterizing rubbery polymer membranes across diverse industrial applications.

organic↗
DOE OSTI2015

Hydrocarbon Liquid Production via Catalytic Hydroprocessing of Phenolic Oils Fractionated from Fast Pyrolysis of Red Oak and Corn Stover

Phenolic oils were produced from fast pyrolysis of two different biomass feedstocks, red oak and corn stover and evaluated in hydroprocessing tests for production of liquid hydrocarbon products. The phenolic oils were produced with a bio-oil fractionating process in combination with a simple water wash of the heavy ends from the fractionating process. Phenolic oils derived from the pyrolysis of red oak and corn stover were recovered with yields (wet biomass basis) of 28.7 wt% and 14.9 wt%, respectively, and 54.3% and 58.6% on a carbon basis. Both precious metal catalysts and sulfided base metal catalyst were evaluated for hydrotreating the phenolic oils, as an extrapolation from whole bio-oil hydrotreatment. They were effective in removing heteroatoms with carbon yields as high as 81% (unadjusted for the 90% carbon balance). There was nearly complete heteroatom removal with residual O of only 0.4% to 5%, while N and S were reduced to less than 0.05%. Use of the precious metal catalysts resulted in more saturated products less completely hydrotreated compared to the sulfided base metal catalyst, which was operated at higher temperature. The liquid product was 42-52% gasoline range molecules and about 43% diesel range molecules. Particulate matter in the phenolic oils complicated operation of the reactors, causing plugging in the fixed-beds especially for the corn stover phenolic oil. This difficulty contrasts with the catalyst bed fouling and plugging, which is typically seen with hydrotreatment of whole bio-oil. This problem was substantially alleviated by filtering the phenolic oils before hydrotreating. More thorough washing of the phenolic oils during their preparation from the heavy ends of bio-oil or on-line filtration of pyrolysis vapors to remove particulate matter before condensation of the bio-oil fractions is recommended.

09 BIOMASS FUELS↗
DOE OSTI2016

Resource Sharing Controls Gene Expression Bursting

Episodic gene expression, with periods of high expression separated by periods of no expression, is a pervasive biological phenomenon. This bursty pattern of expression draws from a finite reservoir of expression machinery in a highly time variant way, i.e., requiring no resources most of the time but drawing heavily on them during short intense bursts, that intimately links expression bursting and resource sharing. Yet, most recent investigations have focused on specific molecular mechanisms intrinsic to the bursty behavior of individual genes, while little is known about the interplay between resource sharing and global expression bursting behavior. In this paper, we confine Escherichia coli cell extract in both cell-sized microfluidic chambers and lipid-based vesicles to explore how resource sharing influences expression bursting. Interestingly, expression burst size, but not burst frequency, is highly sensitive to the size of the shared transcription and translation resource pools. The intriguing implication of these results is that expression bursts are more readily amplified than initiated, suggesting that burst formation occurs through positive feedback or cooperativity. When extrapolated to prokaryotic cells, these results suggest that large translational bursts may be correlated with large transcriptional bursts. This correlation is supported by recently reported transcription and translation bursting studies in E. coli. Finally, the results reported here demonstrate a strong intimate link between global expression burst patterns and resource sharing, and they suggest that bursting plays an important role in optimizing the use of limited, shared expression resources.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Growth Kinetics and Morphology of Barite Crystals Derived from Face-Specific Growth Rates

Here we investigate the growth kinetics and morphology of barite (BaSO 4 ) crystals by measuring the growth rates of the (001), (210), (010), and (100) surfaces using vertical scanning interferometry. Solutions with saturation indices 1.1, 2.1, and 3.0 without additional electrolyte, in 0.7 M NaCl, or in 1.3 mM SrCl2 are investigated. Face-specific growth rates are inhibited in the SrCl 2 solution relative to a solution without electrolyte, except for (100). Contrarily, growth of all faces is promoted in the NaCl solution. The variation of face-specific rates is solution-specific, which leads to a. change of the crystal morphology and overall growth rate of crystals. The measured face-specific growth rates are used to model the growth of single crystals. Modeled crystals have a morphology and size similar to those grown from solution. Based on the model the time dependence of surface area and growth rates is analyzed. Growth rates change with time due to surface area normalization for small crystals and large growth intervals. By extrapolating rates to crystals with large surfaces areas, time-independent growth rates are 0.783, 2.96, and 0.513 mmol∙m -2 ∙h -1 , for saturation index 2.1 solutions without additional electrolyte, NaCl, and SrCl 2 , respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Platinum nanoparticle during electrochemical hydrogen evolution: Adsorbate distribution, active reaction species, and size effect

For small Pt nanoparticles (NPs), catalytic activity is, as observed, adversely affected by size in the 1–3 nm range. We elucidate, via first-principles-based thermodynamics, the operation H* distribution and cyclic voltammetry (CV) during the hydrogen evolution reaction (HER) across the electrochemical potential, including the underpotential region (U ≤ 0) that is difficult to assess in experiment. We consider multiple adsorption sites on a 1 nm Pt NP model and show that the characteristic CV peaks from different H* species correspond well to experiment. We next quantify the activity contribution from each H* species to explain the adverse effect of size. From the resolved CV peaks at the standard hydrogen electrode potential (U = 0), we first deduce that the active species for the HER are the partially covered (100)-facet bridge sites and the (111)-facet hollow sites. Upon evaluation of the reaction barriers at operation H* distribution and microkinetic modeling of the exchange current, we find that the nearest-neighbor (100)-facet bridge site pairs have the lowest activation energy and contribute to ~75% of the NP activity. Edge bridge sites (fully covered by H*) per se are not active; however, they react with neighboring (100)-facet H* to account for ~18% of the activity, whereas (111)-facet hollow sites contribute little. As a result, extrapolating the relative contributions to larger NPs in which the ratio of facet-to-edge sites increases, we show that the adverse size effect of Pt NP HER activity kicks in for sizes below 2 nm.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2012

First-Principles Modeling of Non-Covalent Interactions in Supramolecular Systems: The Role of Many-Body Effects

Supramolecular host–guest systems play an important role for a wide range of applications in chemistry and biology. The prediction of the stability of host–guest complexes represents a great challenge to first-principles calculations due to an interplay of a wide variety of covalent and noncovalent interactions in these systems. In particular, van der Waals (vdW) dispersion interactions frequently play a prominent role in determining the structure, stability, and function of supramolecular systems. On the basis of the widely used benchmark case of the buckyball catcher complex (C60@C60H28), we assess the feasibility of computing the binding energy of supramolecular host–guest complexes from first principles. Large-scale diffusion Monte Carlo (DMC) calculations are carried out to accurately determine the binding energy for the C60@C60H28 complex (26 ± 2 kcal/mol). On the basis of the DMC reference, we assess the accuracy of widely used and efficient density-functional theory (DFT) methods with dispersion interactions. The inclusion of vdW dispersion interactions in DFT leads to a large stabilization of the C60@C60H28 complex. However, DFT methods including pairwise vdW interactions overestimate the stability of this complex by 9–17 kcal/mol compared to the DMC reference and the extrapolated experimental data. A significant part of this overestimation (9 kcal/mol) stems from the lack of dynamical dielectric screening effects in the description of the molecular polarizability in pairwise dispersion energy approaches. The remaining overstabilization arises from the isotropic treatment of atomic polarizability tensors and the lack of many-body dispersion interactions. A further assessment of a different supramolecular system – glycine anhydride interacting with an amide macrocycle – demonstrates that both the dynamical screening and the many-body dispersion energy are complex contributions that are very sensitive to the underlying molecular geometry and type of bonding. We discuss the required improvements in theoretical methods for achieving “chemical accuracy” in the first-principles modeling of supramolecular systems.

Chemistry↗
DOE OSTI2013

Excitation energies with cost-reduced variant of the active-space EOMCCSDT method: the EOMCCSDt-(3) over-bar approach

In this paper we discuss the performance of the several simplified variants of equation-of-motion coupled cluster method (EOMCC) with iterative inclusion of singles, doubles and active-space triples (EOMCCSDt). In particular, we explore simplified EOMCCSDt approaches which enable one to generate the triply excited amplitudes in on-the-fly manner. The original EOMCCSDt formulation has already demonstrated a great success in encapsulating the most important excited-state correlation effects due to triples. In analogy to the original EOMCCSDT formulation, the proposed approach can by-pass the typical bottlenecks associated with the need for storing triply-excited amplitudes. In this paper, we illustrate the performance of several approximate EOMCCSDt methods, named EOMCCSDt-3 and EOMCCSdt-3x, on typical benchmark systems including C2, N2, and the ozone molecules. These new methods yield excitation energies close to the EOMCCSDt ones. The extrapolation of excitation energies for basis sets ranging from cc-pVDZ to cc-pV6Z for N2 and C2 shows very good convergence to the experimental results for states dominated by single excitations.

Coupled cluster↗
DOE OSTI2013

Micromodel Investigation of Transport Effect on the Kinetics of Reductive Dissolution of Hematite

Reductive dissolution of hematite in porous media was investigated using a micromodel with realistic pore network structures that include distinctive advection domain, macro-pores and micro-pores created in silicon substrate. The micromodel pore surface was sputter deposited with a thin layer (230 nm) of hematite. The hematite in the micromodel was reduced by injecting pH-varying solutions containing a reduced form of flavin mononucleotide (FMNH2), a biogenic soluble electron transfer mediator produced by Shewanella species. The reduction kinetics was determined by measuring effluent Fe(II) concentration and by spectroscopically monitoring the hematite dissolution front in the micromodel. Batch experiment was also performed to estimate the hematite reduction rate under the well-mixed condition. The results showed a significant spatial variation in local redox reaction rate that was controlled by the coupled diffusion and reaction. The overall rate of the redox reaction in the micromodel required a three-domain numerical model to effectively describe with distinctive rate parameters in different pore domains. Results from this study demonstrated the important scaling effect when extrapolating geochemical or biogeochemical reaction rate from batch reactor to porous media and indicated a significant control of physical transport mechanisms on the reaction rate scaling.

Reactive transport, porous media, hematite↗
DOE OSTI2014

Effect of Subgrid Heterogeneity on Scaling Geochemical and Biogeochemical Reactions: A Case of U(VI) Desorption

The effect of subgrid heterogeneity in sediment properties on the rate of uranyl[U(VI)] desorption was investigated using a sediment collected from the US Department of Energy Hanford site. The sediment was sieved into 7 grain size fractions that each exhibited different U(VI) desorption properties. Six columns were assembled using the sediment with its grain size fractions arranged in different spatial configurations to mimic subgrid heterogeneity in reactive transport properties. The apparent rate of U(VI) desorption varied significantly in the columns. Those columns with sediment structures leading to preferential transport had much lower rates of U(VI) desorption than those with relatively homogeneous transport. Modeling analysis indicated that the U(VI) desorption model and parameters characterized from well-mixed reactors significantly over-predicted the measured U(VI) desorption in the columns with preferential transport. A dual domain model, which operationally separates reactive transport properties into two subgrid domains improved the predictions significantly. A similar effect of subgrid heterogeneity, albeit at a less degree, was observed for denitrification, which also occurred in the columns. The results imply that subgrid heterogeneity is an important consideration in extrapolating reaction rates from the laboratory to field.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2012

Miscibility gap in the U-Nd-O phase diagram: a new approach of nuclear oxides in the environment?

To some extent, rare-earth-doped UO{sub 2} is representative of an irradiated nuclear fuel. The two phases we observed previously in neodymium-doped UO{sub 2} are now interpreted as the existence of a miscibility gap in the U-Nd-O phase diagram using new results obtained with Raman spectroscopy. Extrapolating the miscibility gap in the U-Nd-O phase diagram to irradiated UO{sub 2} opens the path to a new understanding of nuclear oxides in the environment. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Interactions of Polar and Nonpolar Groups of Alcohols in Zeolite Pores

Understanding the quantitative interactions among zeolite pore walls, Bro̷nsted acid sites, and molecules with both polar and nonpolar regions is essential for scoping out the potential of zeolites as sorbents and catalysts. Purely siliceous zeolites (MFI and Beta in the present study) are hydrophobic, whereas those containing aluminum are considered hydrophilic, preferentially adsorbing organic molecules even in aqueous environments. To characterize these interactions, we use primary alcohols of increasing molecular weight, quantifying their specific interactions in the confined pore space of the alkyl (CH x ) and OH groups. Three types of interactions were identified: (i) alkyl CH x groups interacting with the zeolite pore walls (approximately 10 kJ mol −1 per carbon), (ii) alcohol OH groups interacting with the pore walls (30−35 kJ mol −1 ), and (iii) alcohol OH groups interacting with Bro̷nsted acid sites (37 kJ mol −1 ). All three interactions were well mirrored by computational simulations. The contribution of the alkyl CH x groups was inferred from the incremental increase in sorption enthalpy with increasing molecular weight; the interaction strength of the OH groups was determined by extrapolating the global adsorption enthalpy of the alcohols to a hypothetical OH group without an alkyl group. This value was identical to the adsorption enthalpy of water. The experiments demonstrated that only water has an adsorption enthalpy on zeolite pore walls lower than its condensation enthalpy (30−35 kJ mol −1 vs 45 kJ mol −1 ), limiting the concentration of water that can be adsorbed.

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

In Situ Neutron Diffraction Studies of the Ion Exchange Synthesis Mechanism of Li 2 Mg 2 P 3 O 9 N: Evidence for a Hidden Phase Transition

Motivated by predictions made using a bond valence sum difference map (BVS-DM) analysis, the novel Li-ion conductor Li 2 Mg 2 P 3 O 9 N was synthesized in this paper by ion exchange from a Na 2 Mg 2 P 3 O 9 N precursor. Impedance spectroscopy measurements indicate that Li 2 Mg 2 P 3 O 9 N has a room temperature Li-ion conductivity of about 10 –6 S/cm (comparable to LiPON), which is 6 orders of magnitude higher than the extrapolated Na-ion conductivity of Na 2 Mg 2 P 3 O 9 N at this temperature. The structure of Li 2 Mg 2 P 3 O 9 N was determined from ex situ synchrotron and time-of-flight neutron diffraction data to retain the P2 1 3 space group, though with a cubic lattice parameter of a = 9.11176(8) Å that is significantly smaller than the a = 9.2439(1) Å of Na 2 Mg 2 P 3 O 9 N. The two Li-ion sites are found to be very substantially displaced (~0.5 Å) relative to the analogous Na sites in the precursor phase. The non-molten salt ion exchange method used to prepare Li 2 Mg 2 P 3 O 9 N produces a minimal background in powder diffraction experiments, and was therefore exploited for the first time to follow a Li + /Na + ion exchange reaction using in situ powder neutron diffraction. Lattice parameter changes during ion exchange suggest that the reaction proceeds through a Na 2–x Li x Mg 2 P 3 O 9 N solid solution (stage 1) followed by a two-phase reaction (stage 2) to form Li 2 Mg 2 P 3 O 9 N. However, full Rietveld refinements of the in situ neutron diffraction data indicate that the actual transformation mechanism is more complex and instead involves two thermodynamically distinct solid solutions in which the Li exclusively occupies the Li1 site at low Li contents (stage 1a) and then migrates to the Li3 site at higher Li contents (stage 1b), a crossover driven by the different signs of the local volume change at these sites. Finally, in addition to highlighting the importance of obtaining full structural data in situ throughout the ion exchange process, these results provide insights into the general question of what constitutes a thermodynamic phase.

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

Kinetics of the Strongly Correlated CH 3 O + O 2 Reaction: The Importance of Quadruple Excitations in Atmospheric and Combustion Chemistry

Kinetics measurements on radical–radical reactions are often unavailable experimentally, and obtaining quantitative rate constants for such reactions by theoretical methods is challenging because the transition states and the reactants are often strongly correlated. Treating strongly correlated systems by coupled cluster theory limited to single, double, and triple connected excitations is often inadequate. We therefore use a new method, called GMM(P), for extrapolation to the complete configuration interaction limit to go beyond triple excitations and in particular to approximate the CCSDTQ(P)/CBS limit. Here, we present this method and use it to investigate the CH 3 O + O 2 reaction. The contribution of connected quadruple excitations to the barrier height energy is found to be -3.13 kcal/mol, and adding a quasiperturbative calculation of the effect of connected pentuple excitations brings the post-connected-triples contributions to -3.44 kcal/mol, which corresponds to Boltzmann factors that increase calculated rate constants by factors of 1.0 × 10 3 , 3.3 × 10 2 , and 18 at 250, 298, and 600 K, respectively. We present rate constants for temperatures from 250 to 2000 K, and we find that the Arrhenius activation energy increases from 0.58 to 9.68 kcal/mol over this range. We also find reasonably good accuracy for the barrier height with the MN15-L exchange–correlation functional, and we calculate rate constants by a combination of GMM(P) and MN15-L electronic structure calculations and conventional and variational transition state theory, in particular canonical variational theory with small-curvature tunneling. The present findings have broad implications for obtaining quantitative rate constants for complex reaction systems in atmospheric and combustion chemistry.

Chemistry↗
DOE OSTI2011

Trends in Ln(III) Sorption to Quartz Assessed by Molecular Dynamics Simulations and Laser Induced Flourescence Studies

Molecular dynamics simulations were performed to examine trends in trivalent lanthanide [Ln(III)] sorption to quartz surface SiOH0 and SiO- sites across the 4f period. Complementary laser induced fluorescence studies examined Eu(III) sorption to quartz at varying ionic strength such that the surface sorbed species could be extrapolated at zero ionic strength, the conditions under which the simulations are performed. This allowed for direct comparison of the data, enabling a molecular understanding of the surface sorbed species and the role of the ion surface charge density upon the interfacial reactivity. Thus, this combined theoretical and experimental approach aids in the prediction of the fate of trivalent radioactive contaminants at temporary and permanent nuclear waste storage sites. Potential of mean force molecular dynamics, as well as simulations of pre-sorbed Ln(III) species agrees with the spectroscopic study of Eu(III) sorption, indicating that strongly bound inner-sphere complexes are formed upon sorption to an SiO- site. The coordination shell of the ion contains 6-7 waters of hydration and it is predicted that surface OH groups dissociate from the quartz and bind within the inner coordination shell of Eu(III). Molecular simulations predict less-strongly bound inner2 sphere species in early lanthanides and more strongly bound species in late lanthanides, following trends in the ionic radius of the 4f ions. The participation of surface dissociated OHgroups within the inner coordination shell of the Ln(III) ion is, however, consistent across the series studied. Sorption to a fully protonated quartz surface is not predicted to be favorable by any Ln(III), except perhaps Lu.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Benchmark Theoretical Study of the π–π Binding Energy in the Benzene Dimer

We establish a new estimate for the binding energy between two benzene molecules in the parallel-displaced (PD) conformation by systematically converging (i) the intra- and intermolecular geometry at the minimum, (ii) the expansion of the orbital basis set, and (iii) the level of electron correlation. The calculations were performed at the second-order Møller–Plesset perturbation (MP2) and the coupled cluster including singles, doubles, and a perturbative estimate of triples replacement [CCSD(T)] levels of electronic structure theory. At both levels of theory, by including results corrected for basis set superposition error (BSSE), we have estimated the complete basis set (CBS) limit by employing the family of Dunning’s correlation-consistent polarized valence basis sets. The largest MP2 calculation was performed with the cc-pV6Z basis set (2772 basis functions), whereas the largest CCSD(T) calculation was with the cc-pV5Z basis set (1752 basis functions). The cluster geometries were optimized with basis sets up to quadruple-ζ quality, observing that both its intra- and intermolecular parts have practically converged with the triple-ζ quality sets. The use of converged geometries was found to play an important role for obtaining accurate estimates for the CBS limits. Our results demonstrate that the binding energies with the families of the plain (cc-pVnZ) and augmented (aug-cc-pVnZ) sets converge [within <0.01 kcal/mol for MP2 and <0.15 kcal/mol for CCSD(T)] to the same CBS limit. In addition, the average of the uncorrected and BSSE-corrected binding energies was found to converge to the same CBS limit much faster than either of the two constituents (uncorrected or BSSE-corrected binding energies). Due to the fact that the family of augmented basis sets (especially for the larger sets) causes serious linear dependency problems, the plain basis sets (for which no linear dependencies were found) are deemed as a more efficient and straightforward path for obtaining an accurate CBS limit. We considered extrapolations of the uncorrected (ΔE) and BSSE-corrected (ΔE cp ) binding energies, their average value (ΔE ave ), as well as the average of the latter over the plain and augmented sets (Δ~E ave ) with the cardinal number of the basis set n. Our best estimate of the CCSD(T)/CBS limit for the π–π binding energy in the PD benzene dimer is D e = -2.65 ± 0.02 kcal/mol. The best CCSD(T)/cc-pV5Z calculated value is -2.62 kcal/mol, just 0.03 kcal/mol away from the CBS limit. For comparison, the MP2/CBS limit estimate is -5.00 ± 0.01 kcal/mol, demonstrating a 90% overbinding with respect to CCSD(T). Finally, the spin-component-scaled (SCS) MP2 variant was found to closely reproduce the CCSD(T) results for each basis set, while scaled opposite spin (SOS) MP2 yielded results that are too low when compared to CCSD(T).

Environmental Molecular Sciences Laboratory↗
DOE OSTI2014

Effective Ion Mobility Calculations for Macromolecules by Scattering off Electron Clouds

Broad commercialization and increasing resolving power of ion mobility spectrometry/mass spectrometry (IMS/MS) platforms has engendered an explosion of IMS applications to structural characterization of gas-phase biomolecules. That has renewed interest in more accurate and rapid ion mobility calculations needed to elicit ion geometries from the measurements. An approach based on scattering off electron density isosurfaces (SEDI) that mirrors the physics of molecular collisions was proven superior to the common methods involving atomic coordinates a decade ago, but has remained impractical for large ions because of extreme computational demands. Here, we accelerate SEDI by up to ~500 times using the fragment molecular orbital (FMO) approach for surface generation and the multiplexed scattering algorithm in conjunction with the new grid extrapolation procedure for cross section evaluations. Parallelization of the code on a supercomputer has produced major further speed gains, allowing SEDI calculations for proteins (defined by over a million surface points) with the precision of <0.1% in one minute. Initial tests reveal the anticipated dependences of mobility on the ion charge state and lower cross sections in view of reduced surface roughness. Present developments are expected to lead to broad application of SEDI in IMS studies of macromolecules, enabling more accurate and reliable structural assignments.

Alexeev, Yuri↗