Search NASA⌕ Search

SEARCH · Search NASA

Results for “polarizability”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Meta-Generalized Gradient Approximation for the Cavity-Dependent Exchange-Correlation Interaction in Strongly Coupled Light–Matter Systems

Strong light–matter coupling in optical cavities enables the manipulation of chemical and physical properties without altering molecular composition. Theoretical modeling of such phenomena requires exchange-correlation (XC) functionals that account for both electron–electron and electron–photon (ep) interactions within quantum electrodynamical density functional theory (QEDFT). In this work, we develop a meta-generalized gradient approximation (meta-GGA) specifically targeting the cavity-dependent XC interaction in strongly coupled light–matter systems. This novel approximation is built upon a new semilocal polarizability approximation, which draws from the jellium-with-a-gap model, and can be extended to a “global hybrid” variant that goes beyond the isotropic model from previous approximations. The polarizability model yields significantly improved dispersion coefficients and benchmark calculations with the cavity-dependent XC functional demonstrate improved agreement with QED Hartree–Fock (QED-HF) reference energies. Application to the regioselectivity of brominated nitrobenzene intermediates reveals the functional’s capacity to capture cavity-induced energetic shifts. In conclusion, our results advance the Jacob’s ladder of functionals for QEDFT and provide a practical tool for modeling polaritonic chemistry.

Approximation↗
DOE OSTI2020

Molecular simulation of liquid–vapor coexistence for NaCl: Full-charge vs scaled-charge interaction models

Scaled-charge models have been recently introduced for molecular simulations of electrolyte solutions and molten salts to attempt to implicitly represent polarizability. Although these models have been found to accurately predict electrolyte solution dynamic properties, they have not been tested for coexistence properties, such as the vapor pressure of the melt. In this work, we evaluate the vapor pressure of a scaled-charge sodium chloride (NaCl) force field and compare the results against experiments and a non-polarizable full-charge force field. The scaled-charge force field predicts a higher vapor pressure than found in experiments, due to its overprediction of the liquid-phase chemical potential. Reanalyzing the trajectories generated from the scaled-charge model with full charges improves the estimation of the liquid-phase chemical potential but not the vapor pressure.

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

Time propagation of electronic wavefunctions using nonorthogonal determinant expansions

The use of truncated configuration interaction in real-time time-dependent simulations of electron dynamics provides a balance of computational cost and accuracy, while avoiding some of the failures associated with real-time time-dependent density functional theory. However, low-order truncated configuration interaction also has limitations, such as overestimation of polarizability in configuration interaction singles, even when perturbative doubles are included. Increasing the size of the determinant expansion may not be computationally feasible, and so, in this work, we investigate the use of nonorthogonality in the determinant expansion to establish the extent to which higher-order substitutions can be recovered, providing an improved description of electron dynamics. Model systems are investigated to quantify the extent to which different methods accurately reproduce the (hyper)polarizability, including the high-harmonic generation spectrum of H 2 , water, and butadiene.

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

Non-resonant Raman optical activity from phase-space electronic structure theory

In order to model experimental non-resonant Raman optical activity, chemists must compute a host of second-order response tensors (e.g., the electric-dipole–magnetic-dipole polarizability) and their nuclear derivatives along a set of vibrational modes. While these response functions are almost always computed within a Born–Oppenheimer (BO) framework, here we provide a natural interpretation of the electric-dipole–magnetic-dipole polarizability within phase space electronic structure theory, a beyond-BO model whereby the electronic structure depends on nuclear momentum (P) in addition to nuclear position (R). By coupling to nuclear momentum, phase space electronic structure theory is able to capture the asymmetric response of the electronic properties to an external field, in so far as for a vibrating (non-stationary) molecule, $\frac{∂μ}{∂B}$≠$\frac{∂m}{∂F}$, where μ and m are the electrical linear and magnetic dipoles, and F and B are electric and magnetic fields. As an example, for a prototypical methyloxirane molecule, we show that phase space electronic structure theory is able to deliver a reasonably good match with experimental results in a manner that is formally invariant to gauge origin G 0 —provided that one uses a complete basis or, alternatively, gauge invariant atomic orbitals.

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

Post-remediation geophysical assessment: Investigating long-term electrical geophysical signatures resulting from bioremediation at a chlorinated solvent contaminated site

There is a growing need to assess long-term impacts of active remediation strategies on treated aquifers. A variety of biogeochemical alterations can result from interactions of the amendment with the aquifer, conceivably leading to a geophysical footprint of the long-term alteration of an aquifer. This concept of post-remediation geophysical assessment was investigated in a shallow, chlorinated solvent-contaminated aquifer six to eight years after amendment delivery. Surface resistivity imaging and cross-borehole resistivity and induced polarization (IP) imaging were performed on a transect that spanned treated and untreated zones of the aquifer. Established relationships between IP parameters and surface electrical conductivity were used to predict vertical profiles of electrolytic conductivity and surface conductivity from the inverted cross-borehole images. Aqueous geochemistry data, along with natural gamma and magnetic susceptibility logs, were used to constrain the interpretation. The electrical conductivity structure determined from surface and borehole imaging was foremost controlled by the electrolytic conductivity of the interconnected pore space, being linearly related to fluid specific conductance. The electrolytic conductivity (and thus the conductivity images alone) did not discriminate between treated and untreated zones of the aquifer. In contrast, inverted phase angles and surface conductivities did discriminate between treated and untreated zones of the aquifer, with the treated zone being up to an order of magnitude more polarizable in places. Supporting aqueous chemistry and borehole logging datasets indicate that this geophysical footprint of the long-term impact of the remediation on the aquifer is most likely associated with the formation of polarizable, dispersed iron sulfide minerals. The study opens the door to the possibility of employing time-lapse electrical geophysical measurements to assist with long-term environmental stewardship of legacy sites undergoing active remediation.

bioremediation↗
DOE OSTI2023

Chemistry and Dynamics of Supercritical Carbon Dioxide and Methane in the Slit Pores of Layered Silicates

In the mid 2010s, high-pressure diffraction and spectroscopic tools opened a window into the molecular-scale behavior of fluids at the conditions of many CO 2 sequestration and shale/tight gas reservoirs - conditions where CO 2 and CH 4 are present as variably wet supercritical fluids. Integrating high-pressure spectroscopy and diffraction with molecular modeling has revealed much about the ways that supercritical CO 2 and CH 4 behave in reservoir components, particularly in the slit-shaped micro- and meso-pores of layered silicates (phyllosilicates) abundant in caprocks and shales. This Account summarizes how supercritical CO 2 and CH 4 behave in the slit pores of swelling phyllosilicates as functions of the H 2 O activity, framework structural features, and charge-balancing cation properties at 90 bar and 323 K – conditions similar to a reservoir at ~1 km depth. Slit pores containing cations with large radii, low hydration energy, and large polarizability readily interact with CO 2 , allowing CO 2 and H 2 O to adsorb and co-exist in these interlayer pores over a wide range of fluid humidities. In contrast, cations with small radii, high hydration energy, and low polarizability weakly interact with CO 2 leading to reduced CO 2 uptake and a tendency to exclude CO 2 from interlayers when H 2 O is abundant. The reorientation dynamics of confined CO 2 depends on the interlayer pore height, which is strongly influenced by the cation properties, framework properties, and fluid humidity. The silicate structural framework also influences CO 2 uptake and behavior - for example, smectites with increasing F-for-OH substitution in the framework take up greater quantities of CO 2 . Reactions that trap CO 2 in carbonate phases have been observed in thin H 2 O films near smectite surfaces, including a dissolution-reprecipitation mechanism when the edge surface area is high and an ion exchange-precipitation mechanism when the interlayer cation can form a highly insoluble carbonate. In contrast, supercritical CH 4 does not readily associate with cations, does not react with smectites, and is only incorporated into interlayer slit mesopores when (i) the pore has a $z$-dimension large enough to accommodate CH 4 , (ii) the smectite has low charge and (iii) the H 2 O activity is low. In conclusion, the adsorption and displacement of CH 4 by CO 2 and vice-versa have been studied on the molecular scale in one shale, but opportunities remain to examine behavioral details in this more complicated, slit-pore inclusive system.

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

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2024

Binding of Phosphate Species to Ca 2+ and Mg 2+ in Aqueous Solution

Phosphate derivatives and their interaction with metal cations are involved in many important biological phenomena, so an accurate characterization of the phosphate-metal interaction is necessary to properly understand the role of phosphate-metal contacts in mediating biological function. Herein, we improved the standard 12-6 Lennard-Jones (LJ) potential via the usage of the 12-6-4 LJ model, which incorporates ion-induced dipole interactions. Via parameter scanning, we fine-tuned the 12-6-4 LJ polarizability values to obtain accurate absolute binding free energies for the phosphate anions H 2 PO 4 – , HPO 4 2– , PO 4 3– coordinating with Ca 2+ and Mg 2+ . First, we modified the phosphate 12-6-4 LJ parameters to reproduce the solvation free energies of the series of phosphate anions using the thermodynamic integration (TI) method. Then, using the potential mean force (PMF) method, the polarizability of the metal-phosphate interaction was obtained. We show that the free energy profiles of phosphate ions coordinated to Ca 2+ and Mg 2+ generally show similar trends at longer metal–phosphate distances, while the absolute binding energy values increased with deprotonation. The resulting parameters demonstrate the flexibility of the 12-6-4 LJ-type nonbonded model and its usefulness in accurately describing cation–anion interactions.

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

Probing the Dipole-Bound State in the 9-Phenanthrolate Anion by Photodetachment Spectroscopy, Resonant Two-Photon Photoelectron Imaging, and Resonant Photoelectron Spectroscopy

Valence-bound anions with a dipolar core can support dipole-bound states (DBSs) below the electron detachment threshold. The highly diffuse DBS observed is usually of σ symmetry with an s-like orbital. Recently, a π-type DBS was observed experimentally in the 9-anthrolate anion (9AT – ) and it was shown to be stabilized due to the large anisotropic polarizability of the 9AT core. To confirm the general existence of π-DBS and its structural dependence, here we report an investigation of the 9-phenanthrolate anion (9PT – ), which has a different structure and lower symmetry than 9AT – . Photodetachment spectroscopy revealed a DBS 257 cm -1 below the detachment threshold of 9PT – at 19,627 cm -1 (2.4334 eV). Resonant two-photon photoelectron imaging indeed showed a π symmetry for the DBS. Similar to that observed in 9AT – , the π-DBS in 9PT – is also stabilized by the anisotropic polarizability of the 9PT core and accessed via non-adiabatic population transfer from the initially populated σ-DBS. Photodetachment spectroscopy unveiled nine above-threshold vibrational resonances of the DBS, resulting in nine highly non-Franck-Condon resonant photoelectron spectra by tuning the detachment laser to the vibrational resonances. In conclusion, the combination of photodetachment spectroscopy and resonant photoelectron spectroscopy allowed frequencies for nine vibrational modes of the 9-phenathroxy radical to be measured, including the six lowest frequency bending modes.

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

Efficient Measurement of Length Distribution of 1D Nanoparticles in Solution via Optical Polarimetry

The efficient measurement of the length distribution of nanotubes, nanowires, and other one-dimensional (1D) nanoparticles in solution is important to enable their incorporation into materials and devices and to optimize their processing for properties of interest, such as thermal/electrical conductivity or mechanical strength, in suspensions and composites. We report an electric-field (E-field)-assisted optical-polarimetry technique to measure the length distribution of ensembles of high-aspect-ratio particles in dilute suspension. The degree of alignment of polarizable 1D particles suspended in a fluid under Brownian motion explicitly depends on the E-field strength and the particle length. We show that it is possible to extract the length distribution of 1D nanoparticles suspended in an insulating fluid by applying a range of E-fields and using optical polarimetry to measure the corresponding alignment order parameter. Notably, the method is relatively insensitive to the diameter of the 1D particles, which can be poorly known or vary within a sample. The technique is validated with silver nanowires and carbon nanotubes of known lengths, as well as polymer-depletion-length-separated single-wall carbon nanotube samples with length distributions independently measured with analytical ultracentrifugation. Finally, we demonstrate the ability of the optical-polarimetry technique to quantify changes in the length distribution of ultranarrow, sub-nanometer-diameter single-wall carbon nanotubes under different types and durations of ultrasonication. Within its range of applicability (polarizable 1D nanoparticles in the 0.5 to 15 μm length range, constrained by the voltage stability of the media and the suspended particles), the E-field-assisted optical-polarimetry method is a particularly efficient and accurate method to measure the length distribution of nanowires and nanotubes in suspension.

1D nanoparticles↗
DOE OSTI2025

Algae Asphalt to Enhance Pavement Sustainability and Performance at Subzero Temperatures

This paper evaluates the potential of algae-derived biobinders as sustainable alternatives for pavement construction. It specifically examines the physicochemical and rheological properties of biomodified binders and their potential to offset carbon emissions when used as partial replacements for conventional petroleum-based asphalt binders. Biosequestration of CO 2 using microalgal cell factories is a promising way of recycling CO 2 into biomass via photosynthesis. Our study demonstrates that incorporating algae-derived binders into asphalt can significantly reduce carbon emissions. Each 1% increase in algae-based biobinder leads to an approximate 4.5% decrease in net carbon emissions. This indicates that a blend containing about 22% biobinder has the potential to achieve carbon neutrality. Blends with higher proportions may even result in net-negative emissions, highlighting a promising strategy for environmentally responsible road construction. In terms of performance, the study shows that certain algae-derived biobinders significantly enhance the cracking resistance of asphalt, particularly under subzero temperatures, by improving its stress-relief capacity. A key contribution of this work is the introduction of polarizability as a novel molecular-level parameter for assessing the compatibility of algae-derived bio-oils with asphalt. By capturing the electronic responsiveness of bio-oil molecules, polarizability serves as a predictive indicator of their interaction potential with asphalt components, providing a new dimension for evaluating the binder performance at the molecular scale. Among the tested materials, the biobinder derived from Haematococcus pluvialis demonstrated particularly strong improvements in resistance to permanent deformation under repeated loading conditions analogous to traffic-induced stress, as well as enhanced resistance to moisture-induced damage. In conclusion, these findings advance the chemistry-driven design of biomass-based binders and highlight a promising pathway toward the development of low-carbon, high-performance, and sustainable infrastructure materials.

Algae↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

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

Molecular transport enhancement in pure metallic carbon nanotube porins

Nanofluidic channels impose extreme confinement on water and ions, giving rise to unusual transport phenomena strongly dependent on the interactions at the channel–wall interface. Yet how the electronic properties of the nanofluidic channels influence transport efficiency remains largely unexplored. Here we measure transport through the inner pores of sub-1 nm metallic and semiconducting carbon nanotube porins. We find that water and proton transport are enhanced in metallic nanotubes over semiconducting nanotubes, whereas ion transport is largely insensitive to the nanotube bandgap value. Molecular simulations using polarizable force fields highlight the contributions of the anisotropic polarizability tensor of the carbon nanotubes to the ion–nanotube interactions and the water friction coefficient. We also describe the origin of the proton transport enhancement in metallic nanotubes using deep neural network molecular dynamics simulations. Finally, these results emphasize the complex role of the electronic properties of nanofluidic channels in modulating transport under extreme nanoscale confinement.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Mutually guided light and particle beam propagation

The polarizability of atoms and molecules gives rise to optical forces that trap particles and a refractive index that guides light beams, potentially leading to a self-guided laser and particle beam propagation. In this paper, the mutual interactions between an expanding particle beam and a diffracting light beam are investigated using an axisymmetric particle-light coupled simulation. The nonlinear coupling between particles and photons is dependent on the particle beam radius, particle density, particle velocity and temperature, polarizability, light beam waist, light frequency (with respect to the resonance frequency), and light intensity. The computational results show that the maximum propagation distance is achieved when the waveguiding effect is optimized to single-mode operation. The application of the coupled beam propagation as a space propulsion system is discussed.

42 ENGINEERING↗
DOE OSTI2020

Understanding the thermally activated charge transport in NaPb m SbQ m+2 (Q = S, Se, Te) thermoelectrics: weak dielectric screening leads to grain boundary dominated charge carrier scattering

Many thermoelectric materials feature irregular electrical conductivity with thermally activated transport below ~600 K and metallic behavior at high temperatures, despite possessing degenerate carrier concentrations. The suppression of the electrical conductivity ultimately degrades the thermoelectric performance on the cold side and limits the device energy conversion efficiency. As such, establishing the origin of the low temperature scattering and developing strategies to mitigate its effect are paramount issues. To date, the correct microscopic description of the low temperature carrier scattering remains an open issue, and there is little work addressing why some thermoelectric materials are more susceptible to the deleterious behavior. Here, we use the promising thermoelectric alloys of PbQ and NaSbQ 2 (Q = S, Se, Te) as model systems to address these concerns. We directly show the thermally activated transport stems from the scattering of charge carriers by the grain boundaries (GBs), and that the expected metallic electrical conductivity is recovered by preparing large grained samples with reduced densities of GBs. We furthermore study the electrical properties NaPb m SbSe m+2 as a function of NaSbSe 2 fraction, as well as those of the chalcogenide analogues, PbTe–NaSbTe 2 and PbS–NaSbS 2 , and demonstrate that the strength of GB scattering can be understood by utilizing simple chemical principles. By considering the polarizability of the host atoms, we directly relate the magnitude of GB scattering to the relative degree of charge carrier screening in each material, and demonstrate that GB scattering is strongest in the ionic NaSbQ 2 -rich compounds and weakest in more polarizable PbQ-rich phases. We finally show how these chemical arguments elegantly explain the strong GB scattering in numerous other thermoelectric materials. By uniting the deleterious charge transport properties exhibited by many different compounds into a common picture, we discuss how our work gives design principles for proper microstructure engineering in emerging thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The critical role of configurational flexibility in facilitating reversible reactive metal deposition from borohydride solutions

Development of calcium metal batteries has been historically frustrated by a lack of electrolytes capable of supporting reversible calcium electrodeposition. In this paper, we report the study of an electrolyte consisting of Ca(BH 4 ) 2 in tetrahydrofuran (THF) to gain important insight into the role of the liquid solvation environment in facilitating the reversible electrodeposition of this highly reactive, divalent metal. Through interrogation of the Ca 2+ solvation environment and comparison with Mg 2+ analogs, we show that an ability to reversibly electrodeposit metal at reasonable rates is strongly regulated by dication charge density and polarizability. Our results indicate that the greater polarizability of Ca 2+ over Mg 2+ confers greater configurational flexibility, enabling ionic cluster formation via neutral multimer intermediates. Increased concentration of the proposed electroactive species, CaBH 4 + , enables rapid and stable delivery of Ca 2+ to the electrode interface. Finally, this work helps set the stage for future progress in the development of electrolytes for calcium and other divalent metal batteries.

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

Recent advances in coupled cluster computations of open-shell atomic nuclei

In this contribution, we report on recent progress in coupled-cluster simulations of open-shell atomic nuclei using interactions consistently derived from chiral effective field theory. In particular, we compare different coupled-cluster approaches by computing binding energies and electric dipole polarizabilities in medium-mass calcium isotopes.

Marino, Francesco [Johannes Gutenberg Univ., Mainz↗
DOE OSTI2023

Accurate prediction of global-density-dependent range-separation parameters based on machine learning

In this work, we develop an accurate and efficient XGBoost machine learning model for predicting the global-density-dependent range-separation parameter, ωGDD, for long-range corrected functional (LRC)-ωPBE. This ωGDDML model has been built using a wide range of systems (11 466 complexes, ten different elements, and up to 139 heavy atoms) with fingerprints for the local atomic environment and histograms of distances for the long-range atomic correlation for mapping the quantum mechanical range-separation values. The promising performance on the testing set with 7046 complexes shows a mean absolute error of 0.001 117 a0−1 and only five systems (0.07%) with an absolute error larger than 0.01 a0−1, which indicates the good transferability of our ωGDDML model. In addition, the only required input to obtain ωGDDML is the Cartesian coordinates without electronic structure calculations, thereby enabling rapid predictions. LRC-ωPBE(ωGDDML) is used to predict polarizabilities for a series of oligomers, where polarizabilities are sensitive to the asymptotic density decay and are crucial in a variety of applications, including the calculations of dispersion corrections and refractive index, and surpasses the performance of all other popular density functionals except for the non-tuned LRC-ωPBE. Finally, LRC-ωPBE (ωGDDML) combined with (extended) symmetry-adapted perturbation theory is used in calculating noncovalent interactions to further show that the traditional ab initio system-specific tuning procedure can be bypassed. The present study not only provides an accurate and efficient way to determine the range-separation parameter for LRC-ωPBE but also shows the synergistic benefits of fusing the power of physically inspired density functional LRC-ωPBE and the data-driven ωGDDML model.

Chemistry↗