SEARCH · Engineering Papers

Results for “Large deviations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Can Neutron Star Tidal Effects Obscure Deviations from General Relativity?

Abstract One of the main goals of gravitational-wave astrophysics is to study gravity in the strong-field regime and constrain deviations from general relativity (GR). Any such deviation affects not only binary dynamics and gravitational-wave emission but also the structure and tidal properties of compact objects. In the case of neutron stars, masses, radii, and tidal deformabilities can all differ significantly between different theories of gravity. Currently, the measurement uncertainties in neutron star radii and tidal deformabilities are quite large. However, much less is known about how the large uncertainty in the nuclear equation of state (EOS) might affect tests of GR using binary neutron star mergers. Conversely, using the wrong theory of gravity might lead to incorrect constraints on the nuclear EOS. Here, we study this problem within scalar–tensor (ST) theory. We apply the recently derived ℓ = 2 tidal Love numbers in this theory to parameter estimation of GW170817. Correspondingly, we test if physics beyond GR could bias measurements of the nuclear EOS and neutron star radii. We find that parameter inference for both the GR and ST cases returns consistent component masses and tidal deformabilities. The radius and the EOS posteriors, however, differ between the two theories, but neither is excluded by current observational limits. This indicates that measurements of the nuclear EOS may be biased and that deviations from GR could go undetected when analyzing current binary neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

The Use of the 'Preclosure' Technique for Antegrade Aspiration Thrombectomy with Large Catheters in Acute Limb Ischemia

This study was designed to assess retrospectively short- and mid-term outcomes of the use of a suture-mediated closure device to close the antegrade access in patients undergoing percutaneous aspiration thrombectomy with large catheters for acute leg ischemia. Between November 2005 and February 2010, a suture-mediated active closure system (ProGlide{sup Registered-Sign} 6F, Abbott) was placed before arterial sheath (mean 9 F, range 6-12 F) introduction in 101 patients (74 men, 73 %, mean age 70.1 {+-} 12.6 years standard deviation). Data regarding mortality, complications, and factors contributing to vascular complications at the access site was collected for 6 month after the intervention to detect device-related problems. As a coincidence, 77 patients had follow-up visits for a duplex ultrasound. There were a total of 19 vascular complications (19 %) at the puncture site, all of which were of hemorrhagic nature and none of which consisted of vessel occlusion. Two major outcome complications (2 %) occurred. A retroperitoneal hematoma and a serious inguinal bleeding required additive treatment and did not result in permanent sequelae. Nine cases involved death of which eight were not attributable to the closure and one remained unclear. Successful closure was achieved in 95 patients (94 %); additional manual compression was sufficient in the majority of the remaining patients. Numerous factors contributing to vascular complications were encountered. With acceptable short- and mid-term outcomes, the 'preclose' technique can be a reliable option for the closure of a large antegrade femoral access even for patients at a high risk of vascular complications, such as those undergoing aspiration thrombectomy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2025

Room‐Temperature Mid‐Infrared Detection Using Metasurface‐Absorber‐Integrated Phononic Crystal Oscillator

Mid-infrared (MIR) detectors find extensive applications in chemical sensing, spectroscopy, communications, biomedical diagnosis, and space exploration. Alternative to semiconductor MIR photodiodes and bolometers, mechanical-resonator-based MIR detectors show advantages in higher sensitivity and lower noise at room temperature, especially toward longer wavelength infrared. Here, uncooled room-temperature MIR detectors based on lithium niobate surface acoustic wave phononic crystal (PnC) resonators integrated with wavelength-and-polarization-selective metasurface absorber arrays are demonstrated. The detection is based on the resonant frequency shift induced by the local temperature change due to MIR absorptions. The PnC resonator is configured in an oscillating mode, enabling active readout and low-frequency noise. The 1-GHz oscillator-based MIR detector shows a relative frequency deviation of 5.24 × 10 −10 Hz −1/2 at an integration time of 50 µs, leading to an incident noise equivalent power of 197 pW Hz −1/2 when input 6-µm MIR light is modulated at 1.8 kHz, and a large dynamic range of 10 7 in incident MIR power. The device architecture is compatible with the scalable manufacturing process and can be readily extended to a broader spectral range by tailoring the absorbing wavelengths of metasurface absorbers.

Xi, Zichen [Virginia Polytechnic Inst. and State U↗

DOE OSTI2016

DNS of a turbulent lifted DME jet flame

A three-dimensional direct numerical simulation (DNS) of a turbulent lifted dimethyl ether (DME) slot jet flame was performed at elevated pressure to study interactions between chemical reactions with low-temperature heat release (LTHR), negative temperature coefficient (NTC) reactions and shear generated turbulence in a jet in a heated coflow. By conditioning on mixture fraction, local reaction zones and local heat release rate, the turbulent flame is revealed to exhibit a “pentabrachial” structure that was observed for a laminar DME lifted flame [Krisman et al., (2015)]. The propagation characteristics of the stabilization and triple points are also investigated. Potential stabilization points, spatial locations characterized by preferred temperature and mixture fraction conditions, exhibit autoignition characteristics with large reaction rate and negligible molecular diffusion. The actual stabilization point which coincides with the most upstream samples from the pool of potential stabilization points fovr each spanwise location shows passive flame structure with large diffusion. The propagation speed along the stoichiometric surface near the triple point is compared with the asymptotic value obtained from theory [Ruetsch et al., (1995)]. At stoichiometric conditions, the asymptotic and averaged DNS values of flame displacement speed deviate by a factor of 1.7. However, accounting for the effect of low-temperature species on the local flame speed increase, these two values become comparable. In conclusion, this suggests that the two-stage ignition influences the triple point propagation speed through enhancement of the laminar flame speed in a configuration where abundant low-temperature products from the first stage, low-temperature ignition are transported to the lifted flame by the high-velocity jet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Three-dimensional, Time-Resolved, Intrafraction Motion Monitoring Throughout Stereotactic Liver Radiation Therapy on a Conventional Linear Accelerator

Purpose: To investigate the time-resolved 3-dimensional (3D) internal motion throughout stereotactic body radiation therapy (SBRT) of tumors in the liver using standard x-ray imagers of a conventional linear accelerator. Methods and Materials: Ten patients with implanted gold markers received 11 treatment courses of 3-fraction SBRT in a stereotactic body-frame on a conventional linear accelerator. Two pretreatment and 1 posttreatment cone-beam computed tomography (CBCT) scans were acquired during each fraction. The CBCT projection images were used to estimate the internal 3D marker motion during CBCT acquisition with 11-Hz resolution by a monoscopic probability-based method. Throughout the treatment delivery by conformal or volumetric modulated arc fields, simultaneous MV portal imaging (8 Hz) and orthogonal kV imaging (5 Hz) were applied to determine the 3D marker motion using either MV/kV triangulation or the monoscopic method when marker segmentation was unachievable in either MV or kV images. The accuracy of monoscopic motion estimation was quantified by also applying monoscopic estimation as a test for all treatments during which MV/kV triangulation was possible. Results: Root-mean-square deviations between monoscopic estimations and triangulations were less than 1.0 mm. The mean 3D intrafraction and intrafield motion ranges during liver SBRT were 17.6 mm (range, 5.6-39.5 mm) and 11.3 mm (2.1-35.5mm), respectively. The risk of large intrafraction baseline shifts correlated with intrafield respiratory motion range. The mean 3D intrafractional marker displacement relative to the first CBCT was 3.4 mm (range, 0.7-14.5 mm). The 3D displacements exceeded 8.8 mm 10% of the time. Conclusions: Highly detailed time-resolved internal 3D motion was determined throughout liver SBRT using standard imaging equipment. Considerable intrafraction motion was observed. The demonstrated methods provide a widely available approach for motion monitoring that, combined with motion-adaptive treatment techniques, has the potential to improve the accuracy of radiation therapy for moving targets.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2025

Lithium's low-temperature phase transitions: Insights into quantum lattice dynamics and superconductivity

The large lattice dynamics of lithium, driven by its low atomic mass, results in energetically similar structures and significant isotope effects under pressure, posing challenges to current theoretical models. Above 20 GPa and at low temperatures, lithium's electronic properties deviate from simple metallic behavior, with superconductivity emerging in a complex, pressure-dependent manner, alongside an unusual isotope effect. The structural phases of 7 Li reported under these conditions are inconsistent across studies, and the structures of 6 Li remain unexamined. These gaps limit our understanding of the effects of pressure on lithium's electronic properties and the role of quantum lattice effects on its structural behavior under pressure. Here, we integrate experimental and theoretical approaches to investigate the low-temperature structural phase boundaries in lithium isotopes. We map the structural phase diagram of 7 Li from 5 to 55 GPa and 15–75 K, identifying the sequence fcc → ℎ⁡R⁢1 → cI16. A pronounced isotope effect is observed, with 6 Li shifting the fcc → ℎ⁡R⁢1 phase boundary to lower pressures at 15 K. Density functional theory calculations further clarify how these structural changes affect superconducting properties, particularly emphasizing the role of the fcc → ℎ⁡R⁢1 transition in lithium's superconductivity. Furthermore, our findings offer insights into the unique behavior of lithium isotopes under pressure.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Parity-time symmetry breaking in magnetic systems

The understanding of out-of-equilibrium physics, especially dynamic instabilities and dynamic phase transitions, is one of the major challenges of contemporary science, spanning the broadest wealth of research areas that range from quantum optics to living organisms. By focusing on nonequilibrium dynamics of an open dissipative spin system, we introduce a non-Hermitian Hamiltonian approach, in which non-Hermiticity reflects dissipation and deviation from equilibrium. The imaginary part of the proposed spin Hamiltonian describes the effects of Gilbert damping and applied Slonczewski spin-transfer torque. In the classical limit, our approach reproduces Landau-Lifshitz-Gilbert-Slonczewski dynamics of a large macrospin. Here, we reveal the spin-transfer torque-driven parity-time symmetry-breaking phase transition corresponding to a transition from precessional to exponentially damped spin dynamics. Micromagnetic simulations for nanoscale ferromagnetic disks demonstrate the predicted effect. These findings can pave the way to a general quantitative description of out-of-equilibrium phase transitions driven by spontaneous parity-time symmetry breaking.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

SU-E-T-603: Analysis of Optical Tracked Head Inter-Fraction Movements Within Masks to Access Intracranial Immobilization Techniques in Proton Therapy

Purpose: We present a quantitative methodology utilizing an optical tracking system for monitoring head inter-fraction movements within brain masks to assess the effectiveness of two intracranial immobilization techniques. Methods and Materials: A 3-point-tracking method was developed to measure the mask location for a treatment field at each fraction. Measured displacement of mask location to its location at first fraction is equivalent to the head movement within the mask. Head movements for each of treatment fields were measured over about 10 fractions at each patient for seven patients; five treated in supine and two treated in prone. The Q-fix Base-of-Skull head frame was used in supine while the CIVCO uni-frame baseplate was used in prone. Displacements of recoded couch position of each field post imaging at each fraction were extracted for those seven patients. Standard deviation (S.D.) of head movements and couch displacements was scored for statistical analysis. Results: The accuracy of 3PtTrack method was within 1.0 mm by phantom measurements. Patterns of head movement and couch displacement were similar for patients treated in either supine or prone. In superior-inferior direction, mean value of scored standard deviations over seven patients were 1.6 mm and 3.4 mm for the head movement and the couch displacement, respectively. The result indicated that the head movement combined with a loose fixation between the mask-to-head frame results large couch displacements for each patient, and also large variation between patients. However, the head movement is the main cause for the couch displacement with similar magnitude of around 1.0 mm in anterior-posterior and lateral directions. Conclusions: Optical-tracking methodology independently quantifying head movements could improve immobilization devices by correctly acting on causes for head motions within mask. A confidence in the quality of intracranial immobilization techniques could be more efficient by eliminating the need for frequent imaging.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2018

Search for new phenomena in high-mass final states with a photon and a jet from $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search is performed for new phenomena in events having a photon with high transverse momentum and a jet collected in 36.7 fb -1 of proton–proton collisions at a centre-of-mass energy of s√ = 13 TeV recorded with the ATLAS detector at the Large Hadron Collider. The invariant mass distribution of the leading photon and jet is examined to look for the resonant production of new particles or the presence of new high-mass states beyond the Standard Model. No significant deviation from the background-only hypothesis is observed and cross-section limits for generic Gaussian-shaped resonances are extracted. Excited quarks hypothesized in quark compositeness models and high-mass states predicted in quantum black hole models with extra dimensions are also examined in the analysis. The observed data exclude, at 95% confidence level, the mass range below 5.3 TeV for excited quarks and 7.1 TeV (4.4 TeV) for quantum black holes in the Arkani-Hamed–Dimopoulos–Dvali (Randall–Sundrum) model with six (one) extra dimensions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Search for new phenomena in high-mass final states with a photon and a jet from p p collisions at s = 13 TeV with the ATLAS detector

A search is performed for new phenomena in events with a high-transverse momentum photon and jets collected in 36.7/fb of pp collisions at a centre-of-mass energy of 13 TeV recorded with the ATLAS detector at the Large Hadron Collider. The invariant mass distribution of the leading photon and jet is examined to look for a resonant production of new particles or the presence of new high-mass states beyond the Standard Model. No significant deviation is observed in the mass distribution from the background-only hypothesis and cross section limits for generic Gaussian-shaped resonances are extracted. Excited quarks hypothesized in the quark compositeness model and high-mass states predicted in the quantum black hole models with extra dimensions are also examined in the analysis. The observed data exclude, at 95% confidence level, mass range below 5.3 TeV for the excited quarks and 7.1 TeV (4.4 TeV) for the quantum black holes with six (one) extra dimensions in the Arkani-Hamed-Dimopoulos-Dvali (Randall-Sundrum) model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Evaluating the Benefits of Bayesian Hierarchical Methods for Analyzing Heterogeneous Environmental Datasets: A Case Study of Marine Organic Carbon Fluxes

Large compilations of heterogeneous environmental observations are increasingly available as public databases, allowing researchers to test hypotheses across datasets. Statistical complexities arise when analyzing compiled data due to unbalanced spatial sampling, variable environmental context, mixed measurement techniques, and other reasons. Hierarchical Bayesian modeling is increasingly used in environmental science to describe these complexities, however few studies explicitly compare the utility of hierarchical Bayesian models to simpler and more commonly applied methods. Here we demonstrate the utility of the hierarchical Bayesian approach with application to a large compiled environmental dataset consisting of 5,741 marine vertical organic carbon flux observations from 407 sampling locations spanning eight biomes across the global ocean. We fit a global scale Bayesian hierarchical model that describes the vertical profile of organic carbon flux with depth. Profile parameters within a particular biome are assumed to share a common deviation from the global mean profile. Individual station-level parameters are then modeled as deviations from the common biome-level profile. The hierarchical approach is shown to have several benefits over simpler and more common data aggregation methods. First, the hierarchical approach avoids statistical complexities introduced due to unbalanced sampling and allows for flexible incorporation of spatial heterogeneitites in model parameters. Second, the hierarchical approach uses the whole dataset simultaneously to fit the model parameters which shares information across datasets and reduces the uncertainty up to 95% in individual profiles. Third, the Bayesian approach incorporates prior scientific information about model parameters; for example, the non-negativity of chemical concentrations or mass-balance, which we apply here. We explicitly quantify each of these properties in turn. We emphasize the generality of the hierarchical Bayesian approach for diverse environmental applications and its increasing feasibility for large datasets due to recent developments in Markov Chain Monte Carlo algorithms and easy-to-use high-level software implementations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Potential-pH diagrams considering complex oxide solution phases for understanding aqueous corrosion of multi-principal element alloys

The potential-pH diagram, a graphical representation of the thermodynamically predominant reaction products in aqueous corrosion, is originally proposed for the corrosion of pure metals. The original approach only leads to stoichiometric oxides and hydroxides as the oxidation products. However, numerous experiments show that non-stoichiometric oxide scales are prevalent in the aqueous corrosion of alloys. In the present study, a room temperature potential-pH diagram considering oxide solid solutions, as a generalization of the traditional potential-pH diagram with stoichiometric oxides, is constructed for an FCC single-phase multi-principal element alloy (MPEA) based on the CALculation of PHAse Diagram method. The predominant reaction products, the ions in aqueous solution, and the cation distribution in oxides are predicted. The oxide solid solution is stabilized by the mixing free energy (or mixing entropy) and the stabilizing effect becomes more significant as the temperature increases. Consequently, solid solution oxides are stable in large regions of the potential-pH diagram and the mixing free energy mostly affects the equilibrium composition of the stable oxides, while the shape of stable regions for oxides is mostly determined by the structure of the stable oxides. Agreements are found for Ni 2+ , Fe 2+ , and Mn 2+ between the atomic emission spectroelectrochemistry measurements and thermodynamic calculations, while deviations exist for Cr 3+ and Co 2+ possibly due to surface complexation with species such as Cl - and the oxide dissolution. By incorporating the solution models of oxides, the current work presents a general and more accurate way to analyze the reaction products during aqueous corrosion of MPEAs.

Materials Science↗

DOE OSTI2016

Deformation induced dynamic recrystallization and precipitation strengthening in an Mg−Zn−Mn alloy processed by high strain rate rolling

The microstructure of a high strain-rate rolled Mg−Zn−Mn alloy was investigated by transmission electron microscopy to understand the relationship between the microstructure and mechanical properties. The results indicate that: (1) a bimodal microstructure consisting of the fine dynamic recrystallized grains and the largely deformed grains was formed; (2) a large number of dynamic precipitates including plate-like MgZn{sub 2} phase, spherical MgZn{sub 2} phase and spherical Mn particles distribute uniformly in the grains; (3) the major facets of many plate-like MgZn{sub 2} precipitates deviated several to tens of degrees (3°–30°) from the matrix basal plane. It has been shown that the high strength of the alloy is attributed to the formation of the bimodal microstructure, dynamic precipitation, and the interaction between the dislocations and the dynamic precipitates. - Highlights: •A bimodal microstructure was formed in a high strain-rate rolled Mg−Zn−Mn alloy. •Plate-like MgZn{sub 2}, spherical MgZn{sub 2} and spherical Mn phases were observed. •The major facet of the plate-like MgZn{sub 2} deviated from the matrix basal plane.

36 MATERIALS SCIENCE↗

DOE OSTI2024

DRiFT current mode, trigger settings and flexible detector specifications applied to scintillator arrays

MCNP radiation transport output is post-processed by DRiFT, a Detector Response Function Toolkit to simulate detailed nuclear instrumentation response. DRiFT can be used to assess the performance and potential limitations of scintillator, gas, and semiconductor detectors under a variety of simulated conditions not easily achievable in a laboratory setting. This work describes new updates in DRiFT for scintillator simulations which focus on the capability to simulate scintillators in current mode, an expansion of trigger options, and the ability to customize individual detector properties in a simulation. These improvements are designed to facilitate the ability to model large arrays of scintillator detectors with higher fidelity than was previously possible and are demonstrated in three examples. The first shows the difference between operating DRiFT in current and pulse mode. In the second example, which is intended to demonstrate deviations in individual detector performance, each detector has properties (PMT gain, optical transport, scintillation yield, etc.) that vary between detectors and are specified in DRiFT. A final example examines how DRiFT could be used to optimize digitizer settings in high rate measurements with split signals using the new common trigger option.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2020

Reconstruction of Electron and Ion Distribution Functions in a Magnetotail Reconnection Diffusion Region

In the diffusion region of magnetotail reconnection, particle distributions are highly structured, exhibiting triangular shapes and multiple striations that deviate dramatically from the Maxwellian distribution. Fully kinetic simulations have been demonstrated to be capable of producing the essential structures of the observed distribution functions, yet are computationally not feasible for 3D global simulations. The fluid models used for large-scale simulations, on the other hand, do not have the kinetic physics necessary for describing reconnection accurately. Our study aims to bridge fully kinetic and fluid simulations by quantifying the information required to capture the non-Maxwellian features in the distributions underlying the closures used in the fluid code. Here, we compare the results of fully kinetic simulations with observed electron velocity distributions in a magnetotail reconnection diffusion region and use the maximum entropy model to reconstruct electron and ion distributions using various numbers of moments obtained from the simulation. Our results indicate that using only local moments, the maximum entropy model can reproduce many of the features of the distributions: (1) the electron outflow distribution with a tilted triangular structure is reproduced with 21 or more moments in agreement with Ng et al. and (2) counterstreaming distributions can be captured with the 35-moment model when the separation in velocity space between the populations is large.

58 GEOSCIENCES↗

DOE OSTI2024

Uncovering the phonon spectra and lattice dynamics of plastically deformable InSe van der Waals crystals

Stacking two-dimensional (2D) van der Waals (vdW) materials in a layered bulk structure provides an appealing platform for the emergence of exotic physical properties. As a vdW crystal with exceptional plasticity, InSe offers the opportunity to explore various effects arising from the coupling of its peculiar mechanical behaviors and other physical properties. Here, we employ neutron scattering techniques to investigate the correlations of plastic interlayer slip, lattice anharmonicity, and thermal transport in InSe crystals. Not only are the interlayer slip direction and magnitude well captured by shifts in the Bragg reflections, but we also observe a deviation from the expected Debye behaviour in the heat capacity and lattice thermal conductivity. Combining the experimental data with first-principles calculations, we tentatively attribute the observed evidence of strong phonon-phonon interactions to a combination of a large acoustic-optical frequency resonance and a nesting effect. These findings correlate the macroscopic plastic slip and the microscopic lattice dynamics, providing insights into the mechano-thermo coupling and modulation in 2D vdW materials.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Isotherm Model for Moisture-Controlled CO 2 Sorption

Moisture-controlled sorption of CO 2 , the basis for moisture-swing CO 2 capture from air, is a novel phenomenon observed in strong-base anion exchange materials. Prior research has shown that Langmuir isotherms provide an approximate fit to moisture-controlled CO 2 sorption isotherm data. However, this fit still lacks a governing equation derived from an analytic model. As such, in this paper we derive an analytic form for an isotherm equation from a bottom-up approach, starting with a fundamental theory for an alkali liquid. In the range of interest relevant to CO 2 capture from air, an isotherm equation for an alkali liquid reduces to a simple analytic form with a single parameter, $K_{\text{eq}}$. In the limit $K_{\text{eq}} \gg 1$, a 2nd order approximation simplifies to a Langmuir isotherm that, however, deviates from experimental data. The isotherm theory for an alkali liquid has been generalized to a strong-base anion exchange material. In a strong-base anion exchange material, water concentration inside a sorbent, [H 2 O], is not large enough to be regarded as constant, which allows us to extend $K_{\text{eq}}$ to $K_{\text{eq(AEM)eff}} = K_{\text{eq(AEM)}} \times$ [H 2 O]$^{–n}$ according to the law of mass action. The final isotherm formula has been validated by experimental data from the literature. For a moisture-controlled CO 2 sorbent, $K_{\text{eq(AEM)eff}}$ varies significantly with moisture content of the sorbent. Depending on moisture level, the observed $K_{\text{eq(AEM)eff}}$ in a specific sorbent ranges from a few times to a few thousand times the value of $K_{\text{eq}}$ of a 2 mol L –1 alkali liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Improving intermolecular interactions in DFTB3 using extended polarization from chemical-potential equalization

Semi-empirical quantum mechanical methods traditionally expand the electron density in a minimal, valence-only electron basis set. The minimal-basis approximation causes molecular polarization to be underestimated, and hence intermolecular interaction energies are also underestimated, especially for intermolecular interactions involving charged species. In this work, the third-order self-consistent charge density functional tight-binding method (DFTB3) is augmented with an auxiliary response density using the chemical-potential equalization (CPE) method and an empirical dispersion correction (D3). The parameters in the CPE and D3 models are fitted to high-level CCSD(T) reference interaction energies for a broad range of chemical species, as well as dipole moments calculated at the DFT level; the impact of including polarizabilities of molecules in the parameterization is also considered. Parameters for the elements H, C, N, O, and S are presented. The Root Mean Square Deviation (RMSD) interaction energy is improved from 6.07 kcal/mol to 1.49 kcal/mol for interactions with one charged species, whereas the RMSD is improved from 5.60 kcal/mol to 1.73 for a set of 9 salt bridges, compared to uncorrected DFTB3. For large water clusters and complexes that are dominated by dispersion interactions, the already satisfactory performance of the DFTB3-D3 model is retained; polarizabilities of neutral molecules are also notably improved. Overall, the CPE extension of DFTB3-D3 provides a more balanced description of different types of non-covalent interactions than Neglect of Diatomic Differential Overlap type of semi-empirical methods (e.g., PM6-D3H4) and PBE-D3 with modest basis sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗